Implementation Of Multimode Gaussian Unitaries Using Primitive Components

In the field of continuous-variable systems a fundamental role is played by Gaussian unitaries, that is, operators that preserve the Gaussian character of the states. The paper deals with the implementation of an arbitrary multimode Gaussian unitary by means of primitive components, namely, single-mode operators and simple two-mode operators (beam splitters). A partial solution to this problem is provided by the Bloch-Messiah reduction, which gives an architecture consisting of a multimode displacement, a diagonal squeezer, and two rotation operators. In this architecture displacement and squeezing are already formed by primitive components, so that it remains to find the implementation of the rotation operators. To this end, a suitable factorization of unitary matrices allows one to complete the desired implementation with primitive components. The theory is illustrated with an example of application to a four-mode Gaussian unitary.