Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2603.05824

Gaussian dynamics in the double Siegel disk

Giacomo Pantaleoni, Nicolas C. Menicucci

quant-phSubmitted Mar 6, 2026• Announced Mar 9, 2026• 0 citations
v1 · Mar 6, 2026v2 · Mar 9, 2026
arXivHTMLPDF

Abstract

We show that multimode deterministic (CPTP) Gaussian channels admit a symmetric-space description: by passing from the \(n\)-mode Siegel disk \(Δ_{n}\) to the \(2n\)-mode \emph{double Siegel disk} \(Δ_{2n}\), general Gaussian dynamics becomes a linear-fractional (Möbius) action on a single matrix representative. Concretely, \(Δ_{2n}\) naturally parametrizes Gaussian kernels in the Fock--Bargmann representation, and we identify an explicit physical subset \(\sspace\subsetΔ_{2n}\) corresponding to valid mixed Gaussian states. We then construct, from the standard \((X,Y)\) parametrization of a deterministic Gaussian channel, a normalized oscillator-semigroup element \(\bar E\) whose fractional action implements the channel update \(\amat\mapstoφ_{\bar E}(\amat)\) on \(\sspace\); Gaussian unitaries arise as the symplectic (isometric) subcase. The resulting calculus bridges covariance-matrix channel theory with the adjacency-matrix/symmetric-space picture, retains a simple composition law (matrix multiplication of the acting blocks), and suggests a direct route to graphical update rules beyond the pure-state setting.

Comments
27 pages, 9 figures

Related papers

A common parametrization for finite mode Gaussian states, their symmetries, and associated contractions with some applications

Tiju Cherian John, Kalyanapuram Rangachari Parthasarathy

Let Γ(H) be the boson Fock space over a finite dimensional Hilbert space H. It is shown that every Gaussian symmetry admits a Klauder–Bargmann integral representation in terms of coherent states. Furthermore, Gaussian states, Gaussian symmetries, and second quantization contractions belong to a weakly closed self-adjoint semigroup E2(H) of bounded operators in Γ(H). This yields a common parametrization for these operators. It is shown that the new parametrization for Gaussian states is a fruitful alternative to the customary parametrization by position–momentum mean vectors and covariance matrices. This leads to a rich harvest of corollaries: (i) every Gaussian state ρ admits a factorization ρ=Z1†Z1, where Z1 is an element of E2(H) and has the form Z1=cΓ(P)exp∑r=1nλrar+∑r,s=1nαrsaras on the dense linear manifold generated by all exponential vectors, where c is a positive scalar, Γ(P) is the second quantization of a positive contractive operator P in H, ar, 1 ≤ r ≤ n, are the annihilation operators corresponding to the n different modes in Γ(H), λr∈C, and [αrs] is a symmetric matrix in Mn(C); (ii) an explicit particle basis expansion of an arbitrary mean zero pure Gaussian state vector along with a density matrix formula for a general Gaussian state in terms of its E2(H)-parameters; (iii) a class of examples of pure n-mode Gaussian states that are completely entangled; (iv) tomography of an unknown Gaussian state in Γ(Cn) by the estimation of its E2(Cn) parameters using O(n2) measurements with a finite number of outcomes.

arXiv:1911.06555•Journal of Mathematical Physics•Feb 1, 2021•12 citations

Standard forms and entanglement engineering of multimode Gaussian states under local operations

Alessio Serafini, Gerardo Adesso

Abstract. We investigate the action of local unitary operations on multimode (pure or mixed) Gaussian states and single out the minimal number of locally invariant parametres which completely characterise the covariance matrix of such states. For pure Gaussian states, central resources for continuous-variable quantum information, we investigate separately the parametre reduction due to the additional constraint of global purity, and the one following by the local-unitary freedom. Counting arguments and insights from the phase-space Schmidt decomposition and in general from the framework of symplectic analysis, accompany our description of the standard form of pure n-mode Gaussian states. In particular we clarify why only in pure states with n ≤ 3 modes all the direct correlations between position and momentum operators can be set to zero by local unitary operations. For any n, the emerging minimal set of parametres contains complete information about all forms of entanglement in the corresponding states. An efficient state engineering scheme (able to encode direct correlations between position and momentum operators as well) is proposed to produce entangled multimode Gaussian resources, its number of optical elements matching the minimal number of locally invariant degrees of freedom of general pure n-mode Gaussian states.

•Jan 1, 2007•16 citations

Implementation of Two-Mode Gaussian States Whose Covariance Matrix Has the Standard Form

Gianfranco Cariolaro, Roberto Corvaja

This paper deals with the covariance matrix (CM) of two-mode Gaussian states, which, together with the mean vector, fully describes these states. In the two-mode states, the (ordinary) CM is a real symmetric matrix of order 4; therefore, it depends on 10 real variables. However, there is a very efficient representation of the CM called the standard form (SF) that reduces the degrees of freedom to four real variables, while preserving all the relevant information on the state. The SF can be easily evaluated using a set of symplectic invariants. The paper starts from the SF, introducing an architecture that implements with primitive components the given two-mode Gaussian state having the CM with the SF. The architecture consists of a beam splitter, followed by the parallel set of two single–mode real squeezers, followed by another beam splitter. The advantage of this architecture is that it gives a precise non-redundant physical meaning of the generation of the Gaussian state. Essentially, all the relevant information is contained in this simple architecture.

Symmetry•Jul 20, 2022•4 citations

Geometry and dynamics of Gaussian wave packets and their Wigner transforms

Tomoki Ohsawa, Cesare Tronci

We find a relationship between the dynamics of the Gaussian wave packet and the dynamics of the corresponding Gaussian Wigner function from the Hamiltonian/symplectic point of view. The main result states that the momentum map corresponding to the natural action of the symplectic group on the Siegel upper half space yields the covariance matrix of the corresponding Gaussian Wigner function. This fact, combined with Kostant’s coadjoint orbit covering theorem, establishes a symplectic/Poisson-geometric connection between the two dynamics. The Hamiltonian formulation naturally gives rise to corrections to the potential terms in the dynamics of both the wave packet and the Wigner function, thereby resulting in slightly different sets of equations from the conventional classical ones. We numerically investigate the effect of the correction term and demonstrate that it improves the accuracy of the dynamics as an approximation to the dynamics of expectation values of observables.

arXiv:1611.04081•Journal of Mathematical Physics•Sep 1, 2017•12 citations

Gaussian pure states in quantum mechanics and the symplectic group

R. Simon, E. C. G. Sudarshan, Narasimhaiengar Mukunda

Gaussian pure states of systems with n degrees of freedom and their evolution under quadratic Hamiltonians are studied. The Wigner-Moyal technique together with the symplectic group Sp(2n,openR) is shown to give a convenient framework for handling these problems. By mapping these states to the set of n\ifmmode\times\else\texttimes\fi{}n complex symmetric matrices with a positive-definite real part, it is shown that their evolution under quadratic Hamiltonians is compactly described by a matrix generalization of the M\"obius transformation; the connection between this result and the ``abcd law'' of Kogelnik in the context of laser beams is brought out. An equivalent Poisson-bracket description over a special orbit in the Lie algebra of Sp(2n,openR) is derived. Transformation properties of a special class of partially coherent anisotropic Gaussian Schell-model optical fields under the action of Sp(4, openR) first-order systems are worked out as an example, and a generalization of the ``abcd law'' to the partially coherent case is derived. The relevance of these results to the problem of squeezing in multimode systems is noted.

Physical Review A•Apr 1, 1988•140 citations

From particle counting to Gaussian tomography

Kalyanapuram Rangachari Parthasarathy, Ritabrata Sengupta

The momentum and position observables in an [Formula: see text]-mode boson Fock space [Formula: see text] have the whole real line [Formula: see text] as their spectrum. But the total number operator [Formula: see text] has a discrete spectrum [Formula: see text]. An [Formula: see text]-mode Gaussian state in [Formula: see text] is completely determined by the mean values of momentum and position observables and their covariance matrix which together constitute a family of [Formula: see text] real parameters. Starting with [Formula: see text] and its unitary conjugates by the Weyl displacement operators and operators from a representation of the symplectic group [Formula: see text] in [Formula: see text], we construct [Formula: see text] observables with spectrum [Formula: see text] but whose expectation values in a Gaussian state determine all its mean and covariance parameters. Thus measurements of discrete-valued observables enable the tomography of the underlying Gaussian state and it can be done by using five one-mode and four two-mode Gaussian symplectic gates in single and pair mode wires of [Formula: see text]. Thus the tomography protocol admits a simple description in a language similar to circuits in quantum computation theory. Such a Gaussian tomography applied to outputs of a Gaussian channel with coherent input states permit a tomography of the channel parameters. However, in our procedure the number of counting measurements exceeds the number of channel parameters slightly. Presently, it is not clear whether a more efficient method exists for reducing this tomographic complexity. As a byproduct of our approach an elementary derivation of the probability generating function of [Formula: see text] in a Gaussian state is given. In many cases the distribution turns out to be infinitely divisible and its underlying Lévy measure can be obtained. However, we are unable to derive the exact distribution in all cases. Whether this property of infinite divisibility holds in general is left as an open problem.

arXiv:1504.07054•Infinite Dimensional Analysis Quantum Probability and Related Topics•Dec 1, 2015•25 citations
We show that multimode deterministic (CPTP) Gaussian channels admit a symmetric-space description: by passing from the (n)-mode Siegel disk (Δ_{n}) to the (2n)-mode \emph{double Siegel disk} (Δ_{2n}), general Gaussian dynamics becomes a linear-fractional (Möbius) action on a single matrix representative. Concretely, (Δ_{2n}) naturally parametrizes Gaussian kernels in the Fock--Bargmann representation, and we identify an explicit physical subset (\sspace\subsetΔ_{2n}) corresponding to valid mixed Gaussian states. We then construct, from the standard ((X,Y)) parametrization of a deterministic Gaussian channel, a normalized oscillator-semigroup element (\bar E) whose fractional action implements the channel update (\amat\mapstoφ_{\bar E}(\amat)) on (\sspace); Gaussian unitaries arise as the symplectic (isometric) subcase. The resulting calculus bridges covariance-matrix channel theory with the adjacency-matrix/symmetric-space picture, retains a simple composition law (matrix multiplication of the acting blocks), and suggests a direct route to graphical update rules beyond the pure-state setting.
Let Γ(H) be the boson Fock space over a finite dimensional Hilbert space H. It is shown that every Gaussian symmetry admits a Klauder–Bargmann integral representation in terms of coherent states. Furthermore, Gaussian states, Gaussian symmetries, and second quantization contractions belong to a weakly closed self-adjoint semigroup E2(H) of bounded operators in Γ(H). This yields a common parametrization for these operators. It is shown that the new parametrization for Gaussian states is a fruitful alternative to the customary parametrization by position–momentum mean vectors and covariance matrices. This leads to a rich harvest of corollaries: (i) every Gaussian state ρ admits a factorization ρ=Z1†Z1, where Z1 is an element of E2(H) and has the form Z1=cΓ(P)exp∑r=1nλrar+∑r,s=1nαrsaras on the dense linear manifold generated by all exponential vectors, where c is a positive scalar, Γ(P) is the second quantization of a positive contractive operator P in H, ar, 1 ≤ r ≤ n, are the annihilation operators corresponding to the n different modes in Γ(H), λr∈C, and [αrs] is a symmetric matrix in Mn(C); (ii) an explicit particle basis expansion of an arbitrary mean zero pure Gaussian state vector along with a density matrix formula for a general Gaussian state in terms of its E2(H)-parameters; (iii) a class of examples of pure n-mode Gaussian states that are completely entangled; (iv) tomography of an unknown Gaussian state in Γ(Cn) by the estimation of its E2(Cn) parameters using O(n2) measurements with a finite number of outcomes.
Abstract. We investigate the action of local unitary operations on multimode (pure or mixed) Gaussian states and single out the minimal number of locally invariant parametres which completely characterise the covariance matrix of such states. For pure Gaussian states, central resources for continuous-variable quantum information, we investigate separately the parametre reduction due to the additional constraint of global purity, and the one following by the local-unitary freedom. Counting arguments and insights from the phase-space Schmidt decomposition and in general from the framework of symplectic analysis, accompany our description of the standard form of pure n-mode Gaussian states. In particular we clarify why only in pure states with n ≤ 3 modes all the direct correlations between position and momentum operators can be set to zero by local unitary operations. For any n, the emerging minimal set of parametres contains complete information about all forms of entanglement in the corresponding states. An efficient state engineering scheme (able to encode direct correlations between position and momentum operators as well) is proposed to produce entangled multimode Gaussian resources, its number of optical elements matching the minimal number of locally invariant degrees of freedom of general pure n-mode Gaussian states.
This paper deals with the covariance matrix (CM) of two-mode Gaussian states, which, together with the mean vector, fully describes these states. In the two-mode states, the (ordinary) CM is a real symmetric matrix of order 4; therefore, it depends on 10 real variables. However, there is a very efficient representation of the CM called the standard form (SF) that reduces the degrees of freedom to four real variables, while preserving all the relevant information on the state. The SF can be easily evaluated using a set of symplectic invariants. The paper starts from the SF, introducing an architecture that implements with primitive components the given two-mode Gaussian state having the CM with the SF. The architecture consists of a beam splitter, followed by the parallel set of two single–mode real squeezers, followed by another beam splitter. The advantage of this architecture is that it gives a precise non-redundant physical meaning of the generation of the Gaussian state. Essentially, all the relevant information is contained in this simple architecture.
We find a relationship between the dynamics of the Gaussian wave packet and the dynamics of the corresponding Gaussian Wigner function from the Hamiltonian/symplectic point of view. The main result states that the momentum map corresponding to the natural action of the symplectic group on the Siegel upper half space yields the covariance matrix of the corresponding Gaussian Wigner function. This fact, combined with Kostant’s coadjoint orbit covering theorem, establishes a symplectic/Poisson-geometric connection between the two dynamics. The Hamiltonian formulation naturally gives rise to corrections to the potential terms in the dynamics of both the wave packet and the Wigner function, thereby resulting in slightly different sets of equations from the conventional classical ones. We numerically investigate the effect of the correction term and demonstrate that it improves the accuracy of the dynamics as an approximation to the dynamics of expectation values of observables.
Gaussian pure states of systems with n degrees of freedom and their evolution under quadratic Hamiltonians are studied. The Wigner-Moyal technique together with the symplectic group Sp(2n,openR) is shown to give a convenient framework for handling these problems. By mapping these states to the set of n\ifmmode\times\else\texttimes\fi{}n complex symmetric matrices with a positive-definite real part, it is shown that their evolution under quadratic Hamiltonians is compactly described by a matrix generalization of the M\"obius transformation; the connection between this result and the ``abcd law'' of Kogelnik in the context of laser beams is brought out. An equivalent Poisson-bracket description over a special orbit in the Lie algebra of Sp(2n,openR) is derived. Transformation properties of a special class of partially coherent anisotropic Gaussian Schell-model optical fields under the action of Sp(4, openR) first-order systems are worked out as an example, and a generalization of the ``abcd law'' to the partially coherent case is derived. The relevance of these results to the problem of squeezing in multimode systems is noted.
The momentum and position observables in an [Formula: see text]-mode boson Fock space [Formula: see text] have the whole real line [Formula: see text] as their spectrum. But the total number operator [Formula: see text] has a discrete spectrum [Formula: see text]. An [Formula: see text]-mode Gaussian state in [Formula: see text] is completely determined by the mean values of momentum and position observables and their covariance matrix which together constitute a family of [Formula: see text] real parameters. Starting with [Formula: see text] and its unitary conjugates by the Weyl displacement operators and operators from a representation of the symplectic group [Formula: see text] in [Formula: see text], we construct [Formula: see text] observables with spectrum [Formula: see text] but whose expectation values in a Gaussian state determine all its mean and covariance parameters. Thus measurements of discrete-valued observables enable the tomography of the underlying Gaussian state and it can be done by using five one-mode and four two-mode Gaussian symplectic gates in single and pair mode wires of [Formula: see text]. Thus the tomography protocol admits a simple description in a language similar to circuits in quantum computation theory. Such a Gaussian tomography applied to outputs of a Gaussian channel with coherent input states permit a tomography of the channel parameters. However, in our procedure the number of counting measurements exceeds the number of channel parameters slightly. Presently, it is not clear whether a more efficient method exists for reducing this tomographic complexity. As a byproduct of our approach an elementary derivation of the probability generating function of [Formula: see text] in a Gaussian state is given. In many cases the distribution turns out to be infinitely divisible and its underlying Lévy measure can be obtained. However, we are unable to derive the exact distribution in all cases. Whether this property of infinite divisibility holds in general is left as an open problem.