Unitary matrix decompositions for optimal and modular linear optics architectures

We introduce procedures for decomposing N × N unitary matrices into smaller M × M unitary matrices. Our procedures enable designing modular and optimal architectures for implementing arbitrary discrete unitary transformations on light. Such architectures rely on systematically combining the M-mode linear optical interferometers together to implement a given N-mode transformation. Thus this work enables the implementation of large linear optical transformations using smaller modules that act on the spatial or the internal degrees of freedom of light such as polarization, time or orbital angular momentum. The architectures lead to a rectangular gate structure, which is optimal in the sense that realizing arbitrary transformations on these architectures needs a minimal number of optical elements and minimal circuit depth. Moreover, the rectangular structure ensures that each of the different optical modes incurs balanced optical losses, so the architectures promise substantially enhanced process fidelities as compared to existing schemes. (paper)