Gaussian quantum states of bosonic systems are an important class of states. In particular, they play a key role in quantum optics as all processes generated by Hamiltonians up to second order in the field operators (i.e. linear optics and quadrature squeezing) preserve Gaussianity. A powerful approach to calculations and analysis of Gaussian states is using phase-space variables and symplectic transformations. The purpose of this note is to serve as a concise reference for performing phase-space calculations on Gaussian states. In particular, we list symplectic transformations for commonly used optical operations (displacements, beam splitters, squeezing), and formulae for tracing out modes, treating homodyne measurements, and computing fidelities.
Gaussian quantum states of bosonic systems are an important class of states. In particular, they play a key role in quantum optics as all processes generated by Hamiltonians up to second order in the field operators (i.e. linear optics and quadrature squeezing) preserve Gaussianity. A powerful approach to calculations and analysis of Gaussian states is using phase-space variables and symplectic transformations. The purpose of this note is to serve as a concise reference for performing phase-space calculations on Gaussian states. In particular, we list symplectic transformations for commonly used optical operations (displacements, beam splitters, squeezing), and formulae for tracing out modes, treating homodyne measurements, and computing fidelities.
Gaussian states, operations, and measurements are central building blocks for continuous-variable quantum information processing which paves the way for abundant applications, especially including network-based quantum computation and communication. To make the most use of the Gaussian processes, it is required to understand and utilize suitable mathematical tools such as the symplectic space, symplectic algebra, and Wigner representation. Applying these mathematical tools to practical quantum scenarios, we developed various schemes for quantum transduction, interference-based bosonic mode permutation and bosonic sensing. We demonstrated that generic coupler characterized by Gaussian unitary process can be transformed into a high-fidelity transducer, assuming the access to infinite squeezing and adaptive feedforward with homodyne measurements. To address the practical limitation of finite squeezing, we explored the interference-based protocols. These protocols let us freely permute bosonic modes only assuming the access to single-mode Gaussian operations and multiple uses of a given multi-mode Gaussian process. Thus, such a scheme not only enables universal decoupling for bosonic systems, which is useful for suppressing undesired coupling between the system and the environment, but also faithful bidirectional single-mode quantum transduction. Moreover, noticing that the Gaussian processes are appropriate theoretical models for optical sensors, we studied the quantum noise theory for optical parameter sensing and its potential in providing great measurement precision enhancement. We also extended the Gaussian theories to discrete variable systems, with several examples such as quantum (gate) teleportation. All the analyses originated from the fundamental quantum commutation relations, and therefore are widely applicable.
Covariance matrices lie at the heart of the powerful and well-established symplectic framework for describing continuous-variable Gaussian quantum states. However, since this framework only relies on first- and second-order moments, it is not sufficient for the analysis of non-Gaussian states because their higher-order moments are essential to capture some of their key properties. Here, we define higher-order covariance matrices -- more precisely, covariance matrices built from higher-order quadrature monomials -- which provide a simple way to evaluate the effect of Gaussian transformations on non-Gaussian states and can be used, for example, to address nonlinear squeezing or non-Gaussian nullifiers. Higher-order covariance matrices can be estimated from homodyne measurement data using only a limited number of quadrature angles, which involves matrices of moderate dimension compared with a full simulation in the Fock basis. The dimension of a higher-order covariance matrix does not depend on the span of the quantum states in Fock basis and, furthermore, scales only polynomially with the number of modes.
Quantum continuous variable systems are quantum systems obeying the canonical commutation relations. As we have seen in the previous chapter, such systems always require the adoption of infinite-dimensional Hilbert spaces, even when a finite number of degrees of freedom is considered. Besides making their dynamics potentially intractable, this fact also implies that the study of properties related to the information and correlation content of continuous variable systems turns out to be particularly difficult. A widely known way to make continuous variable systems more tractable is to consider the restriction to Gaussian states. Already well trodden in older research on quantum optics, the study of Gaussian systems has more recently become a central theme of continuous variable quantum information. If open quantum systems, coupled to an external environment, are considered, the Gaussian regime applies when the system couples linearly to its environment and when system and environment are governed by Hamiltonians which are at most quadratic in the canonical operators (or, equivalently, when the overall Hamiltonian of system and environment is at most quadratic). These conditions may seem rather restrictive but, as we shall point out in the course of the book, are actually ordinarily met by a vast number of existing experimental set-ups in the areas of quantum optics, trapped ions, optomechanics, atomic ensembles and certain superconducting degrees of freedom (all of these platforms are reviewed in Chapter 9 ). What we will refer to as Gaussian measurements – POVMs that preserve the Gaussian character of initial states – are also customarily carried out with comparatively high efficiency, while the description of quantum noise can be incorporated seamlessly into the Gaussian picture. What is perhaps even more important to note is that the restriction to the Gaussian realm, while not exhaustive of the full wealth of dynamics and situations allowed in the whole Hilbert space, still lets one include most of the processes relevant to quantum technologies. Squeezing, whereby certain canonical quadratures have uncertainties below the vacuum state noise, and can hence be used in precision measurements, quantum entanglement, and cooling, where the entropy is drained out of a system in order to reach a pure quantum state and thus initialise information protocols or low noise experiments, may all be described within the Gaussian regime. Within the quantum information discourse, the Gaussian restriction has been often criticised on the grounds that no real genuinely quantum effects can ever be observed without leaving it at some point, since the statistics obtained from Gaussian states through Gaussian measurements can be reproduced by classical probability distributions. While technically correct, the reader should be reminded that such critiques often disregard the fact that the object described by a Gaussian state does entail a genuinely quantum description in a Hilbert space. Hence, for instance, the preparation of a pure Gaussian quantum state is the preparation of a pure vector of the Hilbert space. Likewise, a highly squeezed Gaussian state prepared for a metrological protocol is a state with very low noise in a certain physical observable, regardless of whether it can be mimicked by a classical distribution or not. The testing and exploitation of quantum non-locality (stronger than classical correlations that imply the presence of quantum entanglement) with Gaussian states does instead require a departure from Gaussian measurements. However, as we will see, certain non-Gaussian measurements are customarily implemented with current technology, and it is hence still relevant to study in detail the creation of such a resource for Gaussian states. This chapter will introduce Gaussian states from a somewhat unusual angle: as ground and thermal states of at most second-order Hamiltonians. Our treatment will highlight the privileged role played by quadratic canonical ( i.e. , ‘symplectic’) transformations in this context, which we will take care to describe in full generality, and will hinge on the reduction of Gaussian states into their normal modes, which allows for their systematic diagonalisation. We will then proceed to explore an equivalent characterisation in terms of statistical moments of the canonical operators, and to introduce the implications of non-commutativity on such statistics in the form of the Robertson–Schrödinger uncertainty relation.
Gaussian unitaries, generated by quadratic Hamiltonians, are fundamental in quantum optics and continuous-variable computing. Their structures correspond to symplectic (bosons) and orthogonal (fermions) groups, but physical realizations give rise to their respective double covers, introducing phase and sign ambiguities. The homogeneous (quadratic-only) case has been resolved through a parameterization constructed in a recent work [arXiv:2409.11628]. We extend the previous framework to inhomogeneous Gaussian unitaries parameterized by (M,z,Ψ). The Baker-Campbel-Hausdorff formula allows us then to factor any Gaussian unitary into a squeezing and a displacement transformation, from which we derive the group multiplication law.
Gaussian operations, as one usually refers to CP-maps that send Gaussian states to Gaussian states, are apt to model a wide range of situations of practical interest and may often be controlled and applied on demand in the laboratory, in a variety of set-ups (that will be the subject matter of Chapter 9 ). Furthermore, they can be classified in a hierarchy of compact, powerful and far-reaching theoretical descriptions, which will be introduced by degrees and discussed in the current chapter. Note that here we are dealing with discrete quantum operations, rather than dynamical equations, which will be dealt with in the next chapter. We shall start from Gaussian unitary operations, corresponding to the symplectic group already introduced in Chapter 3 , and will also clarify how other basic manipulations such as tensoring and partial tracing carry over to the Gaussian picture. Then, we shall broaden our view to include all deterministic (trace-preserving) maps resulting from Gaussian unitary dilations, which are often known as the set of ‘bosonic Gaussian channels’ and encompass open dynamics subject to noise and ‘decoherence’ (the loss of quantum coherence due to the interaction with an environment). Next, we will move on to consider POVMs corresponding to the well-known homodyne and heterodyne detection schemes which, when acting on a subsystem, preserve the Gaussian character of the initial state both when the system is filtered by recording the measurement outcomes, and when acting unconditionally (that is, when the outcome of the measurement is not recorded). Our description will be further enlarged to the wider class of ‘general-dyne’ detections, and we will also introduce methods to account for imperfect detections. The ability to accommodate noise in the analytical description, often daunting with discrete variable and, more generally, outside the Gaussian regime, is probably one of the greatest strengths of the formalism under discussion. Finally, we shall introduce the Gaussian version of the Choi–Jamiolkowski isomorphism, and attain a unified description of all the Gaussian maps previously defined. Throughout our exposition, we will refer to quantum optical systems to exemplify the practical relevance of the operations encountered.
Gaussian unitaries are specified by a second order polynomial in the bosonic operators, that is, by a quadratic polynomial and a linear term. From the Hamiltonian other equivalent representations of the Gaussian unitaries are obtained, such as Bogoliubov and real symplectic transformations. The paper develops an alternative representation, called complex symplectic transformation, which is more compact and is comprehensive of both Bogoliubov and real symplectic transformations. Moreover, it has other advantages. One of the main results of the theory, not available in the literature, is that the final displacement is not simply given by the linear part of the Hamiltonian, but depends also on the quadratic part. In particular, it is shown that by combining squeezing and rotation, it is possible to achieve a final displacement with an arbitrary amount.