Next Article in Journal
What Features on Routine Panoramic Radiographs Could Help Orthodontists to Estimate the Occurrence of Condylar Hyperplasia from Other Mandibular Asymmetries—Retrospective Analysis Study
Next Article in Special Issue
Implementation of Two-Mode Gaussian States Whose Covariance Matrix Has the Standard Form
Previous Article in Journal
Magnetic Bands Producing a Monoclinic Magnetic Structure in NiO, FeO, MnO, and a Tetragonal One in CoO
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gaussian States: Evaluation of the Covariance Matrix from the Implementation with Primitive Component

by
Gianfranco Cariolaro
1,
Roberto Corvaja
1,* and
Filippo Miatto
2
1
Department of Information Engineering, University of Padova, Via G. Gradenigo 6/B, 35131 Padova, Italy
2
Xanadu Quantum Technologies, 777 Bay Street, Toronto, ON M5G 2C8, Canada
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(7), 1286; https://doi.org/10.3390/sym14071286
Submission received: 25 May 2022 / Revised: 9 June 2022 / Accepted: 15 June 2022 / Published: 21 June 2022
(This article belongs to the Special Issue Theory and Applications of Special Functions II)

Abstract

:
Quantum Gaussian states play a fundamental role in quantum communications and in quantum information. This paper deals with the implementation of multimode, and particularly of two-mode Gaussian unitaries and Gaussian states with primitive components (phase shifters, single-mode real squeezers, displacements, and beam splitters). The architecture thus obtained allows one to obtain an insight into the physical meaning of each variable involved. Moreover, following the implementation architecture, it is possible to formulate an easy algebra (radical free) for the main operations and transformations of Gaussian states.

1. Introduction

In recent years, quantum systems based on continuous variables have attracted great attention for the development of quantum information transmission and processing [1,2,3]. A large amount of theoretical work has been dedicated to multimode quantum states and transformations, with the prominence of Gaussian states and Gaussian transformations which can be more easily implemented and manipulated. In particular, the multimode case is attracting a lot of interest, since it may exhibit entanglement, which is a key resource for several quantum protocols and applications, such as teleportation, computing, and cryptography.
An efficient approach dealing with Gaussian unitaries and Gaussian states is based on the decomposition of the so-called fundamental Gaussian unitaries (FGUs), which are given by displacement, rotation, and squeezing, as shown in Figure 1.
Combination of these three FGUs allows for the implementation of the whole class of multimode Gaussian unitaries and the generation of the whole class of multimode Gaussian states, arriving at explicit closed-form results [4]. Another important approach is obtained resorting to the Bloch–Messiah reduction [1,5] in which the multimode squeezer is decomposed into the parallel of local and real single-mode squeezers, as shown in Figure 2.
There are several tools to efficiently describe Gaussian unitaries/states. For Gaussian unitaries, the main tools are given by Bogoliubov and symplectic transformations; in the N-mode, the former are specified by two complex matrices of order N, while the latter are specified by a 2 N × 2 N real matrix. For Gaussian states, the main tool is given by the covariance matrix, a 2 N × 2 N symmetric matrix with real elements (the mean vector is often neglected in the analysis).
A key role is played by the specification of the variables involved, which determine the degree of freedom of Gaussian unitary/states. For the decomposition of Figure 1, the variables are: the rotation matrix ϕ , the squeeze matrix z , and the displacement vector α ; this is the usual specification for a unitary matrix, which we call for convenience the algebraic specification. For the decomposition of Figure 2, the variables are: two rotation matrices, ϕ and γ , a real diagonal squeeze matrix r D , and the displacement vector α which we call structural or physical specification. Since a Gaussian state is obtained by processing a multimode thermal state through a Gaussian unitary, in order to characterize a Gaussian state, we have to add the N single mode thermal states in the specification.
However, with the algebraic specification, the first step is the evaluation of the polar decomposition of the squeeze matrix z . The formula contains radicals of radicals, which “propagate” on all the subsequent elaborations, leading to cumbersome formulas. This is completely avoided if we start from the structural specification.
The target of the paper is the implementation of Gaussian unitaries/states developing the cascade of Figure 2 with primitive components, starting from the structural specification and using the technique of Bloch–Messiah reduction. The primitive components are rotations, squeezers, displacements in the single mode, and beam splitters [6]. The architecture thus obtained allows one to obtain an insight into the physical meaning of each component involved. Moreover, following the implementation architecture, it is possible to formulate an (radical–free) easy algebra for the main operations and transformations of Gaussian states.
This paper is organized as follows. In Section 2, we formulate the Gaussian unitaries and their decompositions in FUGs, according to the Bloch–Messiah (BM) reduction. In Section 3, we formulate Gaussian states according to Williamson’s theorem. Section 4 deals with the implementation of Gaussian unitaries with primitive components, where the squeeze matrix is decomposed according to the Takagi factorization [7]. The implementation could be considered for multimode [6], but for simplicity, it is carried out for the two-mode one. Continuing with the two-mode, from the architecture with primitive components in Section 5, we evaluate the symplectic transformations and in Section 6, the covariance matrix. As we will realize, all the results are simple and radical free (for comparison, the same evaluation is carried out in Appendix A using the algebraic approach). In the final part of the paper, we focus our attention on the covariance matrix, in which the two-mode depends on 10 real variables. Examples of application are outlined by fixing the variables to specific values. Furthermore, in Section 8, we consider the so-called standard form of the covariance, which only depends on four real variables (symplectic invariants) and contains all the relevant information on the Gaussian state, particularly the entanglement property. We outline a simple method for the evaluation of the symplectic invariants of arbitrary Gaussian states in the two mode.

2. Gaussian Unitaries

2.1. Gaussian Unitary in Terms of Fundamental Unitaries (FGUs)

We remind that a quantum transformation is Gaussian when it transforms Gaussian states into Gaussian states. When the Gaussian transformation is defined by a unitary map, it is called Gaussian unitary.
It is known [2] that Gaussian unitaries can be expressed in the form U = exp ( i H / 2 ) , where the Hamiltonian H is a second-order polynomial in the bosonic operators a i and a i * or in the field operators q i , p i . When these unitaries are applied to the annihilators a = [ a 1 , , a N ] T , we obtain a Bogoliubov transformation, and, in terms of field operators, a symplectic transformation. For the purpose of describing Gaussian unitaries, the two types of transformations are equivalent.
Any Gaussian unitary can be expressed by FGUs. The fundamental unitaries are defined in terms of the column vectors a * and a of size N, collecting the bosonic operators a m * and a m , as:
  • N-mode displacement operator
    D ( α ) : = e α T a * α * a , α = [ α 1 , , a n ] T C N
    where α is called the displacement vector.
  • N-mode rotation operator
    R ( ϕ ) : = e i a * ϕ a , ϕ N × N Hermitian matrix
    where ϕ is called the rotation matrix.
  • N-mode squeeze operator
    S ( z ) : = e 1 2 ( a * z a * a T z * a ) , z N × N symmetric matrix
where z is called the squeeze matrix. The polar decomposition applied to z gives z = r e i θ where r is positive semidefinite and θ is Hermitian symmetric. The matrix r is uniquely determined by z as r = z × z . If r is positive definite, the phase matrix θ associated with the squeeze operator is determined as e i θ = r 1 z . Then [4],
Theorem 1.
The most general Gaussian unitary is obtained by the cascade of the three fundamental Gaussian unitaries D ( α ) , S ( z ) , and R ( ϕ ) , where the order is arbitrary, that is, S ( z ) D ( α ) R ( ϕ ) , R ( ϕ ) D ( α ) S ( z ) , etc.
One can pass from one cascade to any other one using the so-called switching rules, which provide a slight change in the parameters [4]. In the following, without loss of generality, we refer to the cascade
U g = D ( α ) S ( z ) R ( ϕ )
which is illustrated in Figure 1.

2.2. Gaussian Unitaries Decomposed according to the Block–Messiah Reduction

The Block–Messiah reduction [1,5] was reconsidered in [7] in terms of the Takagi factorization [8].
Theorem 2.
The most general Gaussian unitary can be decomposed into the cascade of a rotation operator R ( ψ ) , a local squeeze operator S ( r D ) with a real diagonal matrix r D , another rotation operator R ( γ ) , and finally a displacement operator D ( α ) .
The theorem, illustrated in Figure 2, is essentially based on the Takagi factorization [8] of the squeeze matrix
z = U z r D U z T = e i β r D e i β
which leads to the decomposition of the squeeze operator as S ( z ) = R ( β ) S ( r D ) R ( β ) . Then, the general unitary given by (4), becomes
U g = D ( α ) R ( β ) S ( r D ) R ( ψ )
where R ( ψ ) = R ( β ) R ( ϕ ) (with ψ uniquely determined by e i ψ = e β e i ϕ ). From (6), one can easily obtain the original form (4), in particular the squeeze matrix
z = e i γ r D e i γ T
Starting from other cascades in the application of Theorem 1, the Bloch–Messiah reduction leads to other architectures, but all are equivalent in representing the most general Gaussian unitary.

Specification and Degree of Freedom

Theorem 1 states that the specification of any N-mode Gaussian unitary is given by the matrices α , z , ϕ . Considering their symmetry, the associated degrees of freedom are
2 N 2 + 3 N real variables
As a verification, from Theorem 2 and the BM architecture of Figure 2, for the two-mode, one finds that the number of degrees of freedom is given by 14, which is to be expected, due to the equivalence between the two architectures.

3. Gaussian States

Gaussian states can be obtained from a Gaussian unitary driven by thermal states; in particular, pure Gaussian states are obtained when the thermal states degenerate into vacuum states. To formalize this statement, it is convenient to recall that a Gaussian state is completely described by the covariance matrix and the mean vector. In particular, the covariance matrix can always be written in the form (Williamson’s theorem)
V = S V S T
where S is an N-mode symplectic matrix and
V = diag [ ν 1 , ν 1 , , ν N , ν N ]
corresponds to a tensor product of N thermal states with an average number of thermal photons
N k = ν k 1 2 , k = 1 , 2 , , N
The quantities { ν k } are called the symplectic eigenvalues of the CM V and the operation performed by the matrix S is said to be the symplectic diagonalization of V .
With reference to the previous decompositions of the Gaussian unitaries shown in Figure 1 and Figure 2, we find that, when the input is driven by N thermal states, at the output, we obtain the most general N-mode Gaussian state.
Note that a decomposition into unitaries of the Lie group S p ( 2 N , R ) has also been used in [9,10] for the dynamical symmetry group describing the vibronic transitions in polyatomic molecules.

3.1. Specification and Degree of Freedom

To the specification of the Gaussian unitary, one has to add the specification of the input thermal noise, given by N real variables, as in (9). On the other hand, we have a reduction of the degrees of freedom because a rotation operator is irrelevant for the input thermal noise, such as the rotation R ( ϕ ) of Figure 1, or the rotation R ( ψ ) of Figure 2. Then, the degree of freedom of an N-mode Gaussian state results in
2 N 2 + 3 N + N N 2 = N 2 + 4 N real variables

3.2. Gaussian States in the Two-Mode

Since we focus on two-mode Gaussian states, we review in detail their specification, which is given by the complex matrices (see Figure 1)
α = α 1 α 2 , ϕ = ϕ 11 ϕ 1 , 2 ϕ 12 * ϕ 22 , z = z 11 z 12 z 12 z 22 .
and two thermal noises N 1 , N 2 .
From (8), a two-mode Gaussian unitary has 14 real variables as degrees of freedom, while from (12), a two-mode Gaussian state has 12 real variables degrees of freedom.
We want to evaluate the quantities involved in the two-mode in terms of the parameters given by (13).

4. Implementation with Primitive Components

The implementation of Gaussian unitaries and the generation of Gaussian states with primitive components could be carried out in the general N-mode [6], but hereafter, it will be developed in great detail in the two mode.

4.1. Primitive Components

They are: (1) single-mode displacement; (2) single-mode rotation operators, briefly shifters; (3) single-mode real squeezers; and (4) beam splitters (BSs). Figure 3 shows the graphical symbols we use for these components.
A shifter is specified by a phase β [ 0 , 2 π ) , which leads to the 1 × 1 exponential matrix e i β . A single-mode squeezing operator is specified by the squeeze factor r R . A free-phase BS is defined by a two-mode rotation operator corresponding to the phase matrix with antidiagonal entries
ϕ bs = 0 i ϕ i ϕ 0 R ( ϕ bs ) : = e ϕ ( a 1 * a 2 a 2 a 1 * ) .
where ϕ describe the reflectivity ( s : = sin ϕ ) and the transmittance ( c : = cos ϕ ) of the BS. An equivalent specification is given by the rotation matrix
U bs : = e i ϕ = exp 0 ϕ ϕ 0 = cos ϕ sin ϕ sin ϕ cos ϕ = c s s c .

4.2. Implementation in the Two Mode

The target is the implementation in the two-mode of the architecture of Figure 2 with primitive components. The squeezer is already decomposed into primitive components. The implementation of the two-mode displacement operator D ( α ) is trivial, since it is given by two parallel single-mode displacement operators D ( α 1 ) and D ( α 2 ) . For the rotation operators, we have (for the proof, see the appendixes of [6]).
Proposition 1.
An arbitrary two-mode rotation operator defined by the unitary matrix U = e i ψ : = [ ρ h k e i γ h k ] can be implemented by the cascade of (1) two phase shifters with phases ψ 11 and ψ 12 , (2) a BS with reflectivity s = ρ 12 , and (3) a phase shifter with phase μ = ψ 22 ψ 12 , as illustrated as in Figure 4. (The scheme of Figure 4 is sometimes called beam splitter with phase [11,12].)
The matrix U = e i ψ can be written as
U = e i ψ = c e i ψ 11 s e i ψ 12 s e i ( ψ 11 + μ ) c e i ( ψ 12 + μ ) with μ = ψ 22 ψ 11
and also
U = 1 0 0 e i μ c s s c e i ψ 11 0 0 e i ψ 12
The interpretation of the above factorization proves the proposition.
We are now ready to draw the Gaussian unitary of Theorem 1 with primitive components (see Figure 5). Note that this architecture can be used to specify the whole class of two-mode Gaussian unitaries. It is given by six shifters, two BSs, two real squeezers, and two displacements, corresponding to 14 real variables, in agreement with (8).
The second target is the generation of two-mode Gaussian states. They are obtained when the architecture of Figure 5 is driven by two thermal noises (Figure 6). Now, the shifters γ 11 and γ 12 are removed because they are irrelevant when driven by thermal states. Furthermore, the final displacements were removed because they do not give a contribution to the covariance matrix V (whose evaluation is the main target of the paper).

4.3. Strategy for the Evaluation

The above architectures (Figure 5 and Figure 6) represent the key basis of our theory, which has the goal to understand the action of each primitive component as well as obtain simple and meaningful formulas for the main descriptions of Gaussian unitaries/states, in particular to be “radical free”. This possibility relies ultimately on the Bloch–Messiah reduction.
In the next sections, following the architectures, we will evaluate the symplectic matrix and the covariance matrix. This is called for convenience the structural approach or physical approach, in competition with the algebraic approach more often considered in the literature (in special cases, e.g., with EPR states, the algebraic approach turns out to be easy, but not in the general case).
Here, we develop only the first step of the algebraic approach to see how the complication is originated. From the algebraic specification given by (13), to proceed, one has to obtain the polar decomposition of the squeeze matrix z and the singular values (SVs) of z , which are given by
σ 1 , 2 = 1 2 z 11 2 + z 22 2 + 2 | z 12 | 2 P
where
P = 4 z 11 z 12 * + z 12 z 22 * 2 + z 11 2 z 22 2 2
To proceed, we have to distinguish between the cases of distinct and coincident SVs because the matrix functions have different formulas. This is a further complication. The details of the algebraic approach are given in Appendix A.
Example 1.
We consider squeezing with matrix
z = 5 8 r 0 3 8 3 e i ϵ r 0 3 8 3 e i ϵ r 0 1 8 e 2 i ϵ r 0
In the algebraic approach, the squeeze matrix is one of the data and the first step requires the evaluation of the polar decomposition z = r e i θ , which results in
r = 7 r 0 8 1 8 3 e i ϵ r 0 1 8 3 e i ϵ r 0 5 r 0 8 , e i θ = 1 2 1 2 3 e i ϵ 1 2 3 e i ϵ 1 2 e 2 i ϵ
Additionally, the SVs of z must be evaluated. They result in σ 1 = r 0 / 2 and σ 2 = 2 r 0 .
With the structural approach, the squeeze matrix (20) is obtained with the data: r 1 = r 0 , r 2 = r 0 / 2 , γ 11 = π / 2 , γ 12 = 0 and balanced BSs, so that the free variables are only r 0 and ϵ. However, neither the polar decomposition nor the SVs are needed to proceed.

5. Evaluation of Symplectic Matrices

The Symplectic Matrix (SM)

Here, we consider the real interlaced SM  S , where the phase-space variables are arranged in the form X : = [ q 1 , p 1 , q 2 , p 2 ] T . A symplectic transformation is of the form
X S X + d
where S is a 4 × 4 real matrix and d R 4 . The commutation relations require the condition
S Ω S T = Ω with Ω = i = 1 2 0 1 1 0
A matrix S is called symplectic if it verifies this condition. Following the architecture of Figure 5, we find
S = S rot ( γ ) S sq ( r D ) S rot ( ψ )
where we can use the general formulas
S rot ( ϕ ) = Π cos ϕ sin ϕ sin ϕ cos ϕ Π T
S sq ( r e i θ ) = Π cosh r + sinh r cos θ sinh r sin θ sinh r sin θ cosh r sinh r cos θ Π T
with Π the permutation matrix
Π = 1 0 0 0 0 0 1 0 0 1 0 0 0 0 0 1
For the central squeezer, considering that it is real and diagonal, we find that
cosh r + sinh r 0 0 cosh r sinh r = e r 0 0 e r
where e r = diag e r 1 , e r 2 . Then,
S sq = Π e r 0 0 e r Π T = e r 0 0 e r = e r 1 0 0 0 0 e r 1 0 0 0 0 e r 2 0 0 0 0 e r 2
This completes the evaluation of the global symplectic matrices.
We developed the structural approach for the description of Gaussian unitaries based on the primitive-components architecture of Figure 5. Clearly all the formulas are “radical free”.

6. Evaluation of the Covariance Matrix (CM)

This is an important topic considered by several authors, e.g., [13,14,15,16,17]. In particular, [15] uses a similar approach to characterize Gaussian states, based on Williamson’s theorem and Bloch-Messiah decomposition, describing the symplectic transformations that correspond to the fundamental operators of squeezing, rotation and beam splitters. However, the final results are only provided for the single-mode and a couple of specific examples of two-mode states. In [13], the integration within an ordered product (IWOP) is used to express the covariance matrix of N-mode in terms of the squeeze matrix. Note, however, that the hyperbolic tangent of the squeeze matrix must be evaluated. The final derivation is obtained in two specific examples for the two-mode case. Ref. [14] derives the elements of the covariance matrix of a two-mode Gaussian state, with a generalized squeezing, by using the technique of integration within an ordered product (IWOP). The objective is to characterize the entanglement and non-locality properties of the corresponding state.
Here, we explicitly evaluate the CV V in the two-mode following the structural approach, where V is evaluated from the interlaced SM by adding the information on thermal noise (see (9))
V = S V S T , V = diag [ ν 1 , ν 1 , ν 2 , ν 2 ]
For the explicit evaluation of the trigonometric matrices in the two-mode appearing in the SM S , we use Proposition 1 and (16), with μ = ψ 22 ψ 12 , starting from the exponential
e i ψ = c e i ψ 11 s e i ψ 12 s e i ( ψ 11 + μ ) c e i ( ψ 12 + μ )
Then
( cos ( ψ ) = c cos ( ψ 11 ) s cos ( ψ 12 ) s cos ( ψ 11 + μ ) c cos ( ψ 12 + μ ) sin ( ψ ) = c sin ( ψ 11 ) s sin ( ψ 12 ) s sin ( ψ 11 + μ ) c sin ( ψ 12 + μ )
Analogously
( cos ( γ ) = q cos ( γ 11 ) p cos ( γ 12 ) p cos ( γ 11 + ϵ ) q cos ( γ 12 + ϵ ) sin ( γ ) = q sin ( γ 11 ) p sin ( γ 12 ) p sin ( γ 11 + ϵ ) q sin ( γ 12 + ϵ )
where p is the reflectivity of the second BS and q = 1 p 2 .
In order to characterize the CM, it is more convenient to express the result partitioned into 2 × 2 blocks. Letting
S = S 11 S 12 S 21 S 22
we obtain:
V = A C C T B
where
( A = ν 1 S 11 S 11 T + ν 2 S 12 S 12 T B = ν 1 S 21 S 21 T + ν 2 S 22 S 22 T C = ν 1 S 11 S 21 T + ν 2 S 12 S 22 T
Considering that the CM refers to Gaussian states, as remarked above in Figure 6, the dependence is on the 10 real parameters ν 1 , ν 2 , s , μ , r 1 , r 2 , γ 11 , γ 12 , p , ϵ , notwithstanding that the S i j also depend on the phases ψ 11 , ψ 12 .
A Gaussian state then depends on 10 real + 2 complex parameters, namely ν 1 , ν 2 , s, μ , r 1 , r 2 , γ 11 , γ 12 , p, ϵ , α 1 , α 2 .

6.1. Standard Form of the Covariance Matrix

We remind that, for any two-mode Gaussian states, there exists a local symplectic operation S l that brings the covariance matrix to the standard form, Simon2000, [18]
S l T V S l = V s f = a 0 c + 0 0 a 0 c c + 0 b 0 0 c 0 b
where the correlation terms a , b , c + and c are obtained from the ordinary CM V by the four local symplectic invariants
det V = ( a b c + 2 + ) ( a b c 2 ) , det A = a 2 , det B = b 2 , det C = c + c
The importance of the standard form and of the symplectic invariants lies on the fact that they concentrate all the relevant information on the two-mode Gaussian states, particularly that concerning the entanglement [19]. This means that the essential degree of freedom is reduced to 4 real variables, instead of the 10 variables listed above.
Evaluation of “invariants” ( a , b , c ± ) from the ordinary CM V . The invariants can be obtained from the blocks of the ordinary (nonstandard) covariance matrix (see (30)).
Proposition 2.
From the blocks of the ordinary covariance matrix V , evaluate the quantities
Δ V : = det V , Δ C : = det C , R = a 2 b 2 Δ V + Δ C 2 2 4 a 2 b 2 Δ C 2
Then, the invariants result in: a = det A , b = det B , and
c + = R + a 2 b 2 Δ V + Δ C 2 2 a b , c = 2 a b Δ C R + a 2 b 2 Δ V + Δ C 2
The proof is in Appendix B.
We give a numerical example. With the data
( ν 1 = 3.2 , ν 2 = 1.6 , μ = 0.999 , γ 11 = 0.9 , γ 12 = 0.66 , ϵ = 0.7 , r 1 = 0.66 , r 2 = 0.18 , p = 0.3 , q = 0.953939 , s = c = 0.707
A = 3.523 3.538 3.538 4.829 ,   B = 1.821 0.160 0.160 4.565 ,   C = 0.682 1.835 0.695 1.860
Then
det A = 4.49884 , det B = 8.28805 , det C = 0.00655408 , det V = 26.2144
a = 2.12105 , b = 2.8789
The four solutions are
c + 1.34656 , c 0.00486726 , c + 0.00486726 , c 1.34656
c + 0.00486726 , c 1.34656 , c + 1.34656 , c 0.00486726

6.2. Physical Analysis of the Global Architectures

The architectures of Figure 5 and Figure 6 allow us to obtain a physical insight into the dependence of the global performance, as stated by the symplectic matrix and by the covariance matrix, from each primitive component. In the limit, we can remove a single component from the architecture to see the effect on the ordinary covariance matrix V . In a forthcoming paper [20], we will show that a quantum state, having as ordinary covariance just the standard covariance matrix V s f , can be obtained by removing all the shifters. This shows that the shifters have less relevance with respect to the other primitive components.

7. Examples

By choosing the 14 parameters listed above, one can evaluate the description of all Gaussian unitaries in the two-mode. Analogously, by choosing the 10 parameters of Figure 6, one can evaluate the description of all two-mode Gaussian states. Here, we outline a few cases of Gaussian states focusing the attention on the covariance matrix.
Case 1: All-zero phases
A first case is considered with all zero phases, with parameters as in Table 1.
The squeeze matrix results in
z = r 2 p 2 + q 2 r 1 p q r 2 p q r 1 p q r 2 p q r 1 r 1 p 2 + q 2 r 2
The blocks which constitute the covariance matrix V are given by
A = ν 1 c q e r 1 p s e r 2 2 + c p e r 2 + q s e r 1 2 ν 2 0 0 ν 1 c q e r 1 p s e r 2 2 + c p e r 2 + q s e r 1 2 ν 2
B = ν 2 c q e r 2 p s e r 1 2 + c p e r 1 + q s e r 2 2 ν 1 0 0 ν 2 c q e r 2 p s e r 1 2 + c p e r 1 + q s e r 2 2 ν 1
C = c 1 0 0 c 2
where
( c 1 = c e r 2 q e r 1 p s c e r 2 p + e r 1 q s ν 2 c e r 1 q e r 2 p s c e r 1 p + e r 2 q s ν 1 c 2 = c e r 1 r 2 ( p 2 q 2 ) s ( ν 1 ν 2 ) + e 2 r 2 p q ν 2 c 2 + s 2 ν 1 e 2 r 1 p q ν 1 c 2 + s 2 ν 2
Case 2: All–zero phases, r 1 = r 2 = r 0 , and the second BS balanced
We consider another case, with all zero phases, the same squeezing on each mode and balanced BS, as in Table 2.
The squeeze matrix results in
z = r 0 0 0 r 0
The blocks of the covariance matrix V are given by
A = 1 2 e 2 r 0 ν 1 ( c s ) 2 + ( c + s ) 2 ν 2 0 0 1 2 e 2 r 0 ν 1 ( c s ) 2 + ( c + s ) 2 ν 2
B = 1 2 e 2 r 0 ν 2 ( c s ) 2 + ( c + s ) 2 ν 1 0 0 1 2 e 2 r 0 ν 2 ( c s ) 2 + ( c + s ) 2 ν 1
C = 1 2 e 2 r 0 ( c s ) ( c + s ) ( ν 1 ν 2 ) 0 0 1 2 e 2 r 0 ( c s ) ( c + s ) ( ν 1 ν 2 )
Case 3: EPR state with noise
The EPR unitary is a squeeze with the following matrix
z = 0 e i θ r 0 e i θ r 0 0
To obtain this unitary, the primitive-component architecture must have the following data r 1 = r 2 = r 0 , balanced BSs p = q = s = c = 1 2 and the phases indicated in Table 3.
The blocks of the covariance matrix V result in
( A = ν 1 cosh 2 ( r 0 ) + ν 2 sinh 2 ( r 0 ) 0 0 ν 1 cosh 2 ( r 0 ) + ν 2 sinh 2 ( r 0 ) B = ν 2 cosh 2 ( r 0 ) + ν 1 sinh 2 ( r 0 ) 0 0 ν 2 cosh 2 ( r 0 ) + ν 1 sinh 2 ( r 0 ) C = cos ( θ ) cosh ( r 0 ) ( ν 1 + ν 2 ) sinh ( r 0 ) cosh ( r 0 ) ( ν 1 + ν 2 ) sin ( θ ) sinh ( r 0 ) cosh ( r 0 ) ( ν 1 + ν 2 ) sin ( θ ) sinh ( r 0 ) cos ( θ ) cosh ( r 0 ) ( ν 1 + ν 2 ) sinh ( r 0 )
In particular, for θ = 0 , the block C results in
C = cosh ( r 0 ) ( ν 1 + ν 2 ) sinh ( r 0 ) 0 0 cosh ( r 0 ) ( ν 1 + ν 2 ) sinh ( r 0 )
in agreement with the result indicated in [19].

8. Separability and Entanglement

The separability and the entanglement properties of any two-mode Gaussian state can be handled from both the ordinary CM V and the standard CM V s f . Now, a Gaussian state is separable if its standard CM is diagonal. Hence, from (32), the separability condition is
c + = 0 , c = 0
When the standard CM is evaluated from the ordinary CM using Proposition 2, Equation (35) is not useful because it leads to degeneracy. Instead of (33), namely det V = ( a b c + 2 + ) ( a b c 2 ) and det C = c + c , one finds the conditions on the symplectic invariants
det V s f = a 2 b 2 , det C = 0
On the other hand, the separability based on the ordinary CM is that V must be diagonal. Then, from (30) and (31), the conditions read
A , B diagonal , C = 0
Recalling that ν 1 1 , ν 2 1 , from (31), we find the conditions
( S 11 S 11 T , S 12 S 12 T , S 21 S 21 T , S 22 S 22 T diagonal S 11 S 21 T = 0 , S 12 S 22 T = 0

9. Conclusions

For the description of Gaussian unitaries/states, we considered the algebraic and structural approaches. We remark that the class of noisy Gaussian states generated by the algebraic specification and by the structural specification coincide, since both approaches generate the whole class of noisy Gaussian states. In other words, one can choose the first class as well as the second class as the definition of noisy Gaussian states. In terms of performance, both the methods and also alternative solutions proposed in [13,14,15] achieve the same covariance matrix, since no approximations are introduced, but at different costs, in terms of complexity and often only in some particular cases. In particular, in [15], the solution is provided only for a couple of specific examples. In [13], the covariance matrix is obtained considering only a squeezed state, not a general one, by means of integration within an ordered product (IWOP), starting from the knowledge of the correlation operator. Note also that the hyperbolic tangent of the squeeze matrix must be evaluated. Furthermore, in [14], only the squeezing is considered, again with the same complexity of IWOP, and the objective is to characterize the entanglement of the corresponding state.
The main target of the paper is to show the advantages of the structural approach in the description of the most generic Gaussian states, namely:
  • The structural approach is completely radical free compared to the algebraic approach and requires several matrix operations leading to results that contain the radicals of radicals. Note that the key to avoiding radicals is the following: the quantities which exhibit radicals in the algebraic approach become independent variables (data) in the structural approach.
  • The structural approach is completely general, while the algebraic approach exhibits several degeneracies (mainly coincident eigenvalues) concerning some very important cases (see EPR states). These cases should be treated separately with ad hoc procedures. Such a distinction is not required in the structural approach.
  • In the structural approach, all the variables have a precise physical meaning, related to the corresponding components of the architecture, i.e., squeezers, beam splitters, phase shifters, and one can choose the specific variables to achieve the desired properties of the covariance matrix, for example, the entanglement.
Finally, we note that the theory, here developed in detail for the two-mode, could be extended to higher modes, of course with the penalty of complication at the increase in the order. We developed the three–mode (not reported here) with the structural approach and without difficulty, resulting, however, in long formulas (which were radical free). With the algebraic approach, the three mode results to be very complicated, mainly for the large number of particular cases to be treated separately.

Author Contributions

Formal analysis, G.C., R.C. and F.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors are especially indebted to Giuseppe Dattoli for their numerous comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Algebraic Approach

The purpose of this appendix is mainly to show the complexity of the algebraic approach.
For the evaluation of the symplectic matrices S rot and S sq (see (25) and (26)) needed in the covariance formula, several preliminary evaluations are requested.

Appendix A.1. Polar Decomposition of the Squeeze Matrix z

This requires the evaluation of the SVs of z. With the form (13), the SVs are
σ 1 , 2 = 1 2 z 11 2 + z 22 2 + 2 | z 12 | 2 P
where
P = 4 z 11 z 12 * + z 12 z 22 * 2 + z 11 2 z 22 2 2
We have to distinguish between the cases of distinct and coincident SVs because the matrix functions have different formulas. From (18), the conditions of coincident eigenvalues are | z 11 | = | z 22 | and z 11 z 12 * + z 12 z 22 * = 0 . We suppose that the matrix z is not diagonal (because the diagonal case is trivial with matrix functions). Hence, with z 12 0 , the condition of coincident SVs becomes z 22 = z 11 * z 12 / z 12 * .
The polar decomposition r e i θ is obtained considering that z is complex symmetric, r is Hermitian PSD with eigenvalues equal to the SVs of z , and e i θ is unitary. Then, with distinct SVs, one finds
r = 1 σ 1 + σ 2 z 11 2 + z 12 2 + σ 1 σ 2 z 11 z 12 * + z 12 z 22 * z 12 z 11 * + z 22 z 12 * z 12 2 + z 22 2 + σ 1 σ 2
e i θ = a 11 a 11 2 2 a 12 2 + J a 12 2 a 22 * σ 1 σ 2 ( σ 1 + σ 2 ) ( J M ) a 12 a 11 a 22 a 12 * σ 1 σ 2 ( σ 1 + σ 2 ) ( J M ) a 12 a 11 a 22 a 12 * σ 1 σ 2 ( σ 1 + σ 2 ) a 22 2 a 12 2 a 22 2 + J a 12 2 a 11 * σ 1 σ 2 ( σ 1 + σ 2 )
where
J = σ 1 2 + σ 2 σ 1 + σ 2 2 , M = z 11 2 + z 12 2 + z 22 2
With coincident SVs and the nondiagonal z , the PD reads
r = σ 1 0 0 σ 1 , e i θ = z 11 σ 1 z 12 σ 1 z 12 σ 1 z 12 z 11 * z 12 * σ 1
where the SV is given by σ 1 = z 11 2 + z 12 2 .

Appendix A.2. Evaluation of the Matrices cosh(r) and sinh(r)

Considering that the eigenvalues of r are the singular values of z, that is, σ 1 , σ 2 , and using the expression of r given by (A3), with distinct SVs, one obtains
G ( r ) = G ( σ 2 ) a 11 2 a 12 2 + σ 1 2 + G ( σ 1 ) a 11 2 + a 12 2 σ 2 2 σ 1 2 σ 2 2 a 11 a 12 * + a 12 a 22 * ( G ( σ 1 ) G ( σ 2 ) ) σ 1 2 σ 2 2 a 12 a 11 * + a 22 a 12 * ( G ( σ 1 ) G ( σ 2 ) ) σ 1 2 σ 2 2 G ( σ 2 ) a 12 2 a 22 2 + σ 1 2 + G ( σ 1 ) a 12 2 + a 22 2 σ 2 2 σ 1 2 σ 2 2
Hence, with G ( x ) = cosh ( x ) and G ( x ) = sinh ( x ) , we can obtain the desired matrix functions. With coincident eigenvalues, one obtains
cosh ( r ) = cosh ( σ 1 ) I 2 , sinh ( r ) = sinh ( σ 1 ) I 2 .

Appendix A.3. Evaluation of the Phase Matrix ϕ

We first evaluate the eigenvalues of ϕ , which are given by
μ 1 , 2 = 1 2 4 ϕ 12 2 + ( ϕ 11 ϕ 22 ) 2 + ϕ 11 + ϕ 22 .
Assuming distinct eigenvalues, we find
e i ϕ = 1 ψ 11 ψ 21 e i ψ 21 ( ψ 11 ϕ 11 ) + e i ψ 11 ( ϕ 11 ψ 21 ) e i ψ 11 e i ψ 21 ϕ 12 e i ψ 11 e i ψ 21 ϕ 12 * e i ψ 21 ( ψ 11 ϕ 22 ) + e i ψ 11 ( ϕ 22 ψ 21 )
The eigenvalues (A9) become coincident when ϕ 11 = ϕ 22 and ϕ 12 = 0 . Then, ϕ = ϕ 11 I 2 and the rotation matrix becomes
e i ϕ = e i ϕ 11 I 2

Appendix A.4. Evaluation of Matrices cosϕ and sinϕ

We can start from (A10) and use the identities
cos ϕ = 1 2 ( e i ϕ + e i ϕ ) , sin ϕ = 1 2 i ( e i ϕ e i ϕ )
From (18), the conditions of coincident eigenvalues are | z 11 | = | z 22 | and z 11 z 12 * + z 12 z 22 * = 0 . We suppose that the matrix z is not diagonal (because the diagonal case is trivial with matrix functions). Hence, with z 12 0 , the condition of coincident SVs becomes z 22 = z 11 * z 12 / z 12 * .

Appendix B

Proof of Proposition 2.
The equations are given by
det V = ( a b c + 2 ) ( a b c 2 ) , det A = a 2 , det B = b 2 , det C = c + c
where the determinants can be calculated from the data. With det V = H , det C = K , the equations are
H = a b c 2 a b d 2 , K = c d
where c = c + and d = c are the unknowns. Letting
R = a 2 b 2 H + K 2 2 4 a 2 b 2 K 2
we find four solutions given by
c = R + a 2 b 2 H + K 2 2 a b , d = K c
c = R + a 2 b 2 H + K 2 2 a b , d = K c

References

  1. Braunstein, S.L.; van Look, P. Quantum information with continuous variables. Rev. Mod. Phys. 2005, 77, 513. [Google Scholar] [CrossRef] [Green Version]
  2. Weedbrook, C.; Pirandola, S.; Garcia-Patron, R.; Cherf, N.J.; Ralf, T.C.; Shapiro, J.H.; Lloyd, S. Gaussian Quantum Information. Rev. Mod. Phys. 2012, 84, 621. [Google Scholar] [CrossRef]
  3. Babusci, D.; Dattoli, G.; Riciardi, S.; Sabia, E. Mathematical Methods for Physicists; World Scientific Publishing: Singapore, 2020. [Google Scholar]
  4. Ma, X.; Rhodes, W. Multimode squeeze operators and squeezed states. Phys. Rev. A 1990, 41, 4625–4631. [Google Scholar] [CrossRef] [PubMed]
  5. Bloch, C.; Messiah, A. The canonical form of an antisymmetric tensor and its application to the theory of superconductivity. Nucl. Phys. 1964, 39, 95–106. [Google Scholar] [CrossRef]
  6. Cariolaro, G.; Pierobon, G. Implementation of multimode Gaussian unitaries using primitive components. Phys. Rev. A 2018, 98, 032111. [Google Scholar] [CrossRef]
  7. Cariolaro, G.; Pierobon, G. Bloch-Messiah reduction of Gaussian unitaries by Takagi factorization. Phys. Rev. A 2016, 94, 062109. [Google Scholar] [CrossRef]
  8. Horn, R.A.; Johnson, C.J. Matrix Analysis; Cambridge University Press: New York, NY, USA, 1985. [Google Scholar]
  9. Doktorov, E.V.; Malkin, I.A.; Man’ko, V.I. Dynamical symmetry of vibronic transitions in polyatomic molecules and the Franck-Condon principle. J. Mol. Spectrosc. 1975, 56, 1–20. [Google Scholar] [CrossRef]
  10. Doktorov, E.V.; Malkin, I.A.; Man’ko, V.I. Dynamical symmetry of vibronic transitions in polyatomic molecules and the Franck-Condon principle. J. Mol. Spectrosc. 1977, 64, 302–326. [Google Scholar] [CrossRef]
  11. Campos, R.A.; Saleh, B.E.A.; Teich, M.C. Quantum-mechanical lossless beam splitter: SU(2) symmetry and photon statistics. Phys. Rev. A 1989, 40, 1371–1384. [Google Scholar] [CrossRef] [PubMed]
  12. Scheel, J.; Welsch, D.G. Entanglement generation and degradation by passive optical devices. Phys. Rev. A 2001, 64, 063811. [Google Scholar] [CrossRef] [Green Version]
  13. He, R. New algorithm to calculate the covariance matrix of an arbitrary form of Gaussian state. Quantum Inf. Process. 2015, 14, 3971–3981. [Google Scholar] [CrossRef]
  14. Hu, L.-Y.; Xu, X.-X.; Guo, Q.; Fan, H.-Y. Entanglement and nonlocality of one- and two-mode combination squeezed state. Opt. Commun. 2010, 283, 5074–5080. [Google Scholar] [CrossRef] [Green Version]
  15. Adesso, G.; Ragy, S.; Lee, A.R. Continuous Variable Quantum Information: Gaussian States and Beyond. Open Syst. Inf. Dyn. 2014, 21, 1440001. [Google Scholar] [CrossRef] [Green Version]
  16. Paris, M.G.A.; Illuminati, F.; Serafini, A.; de Siena, S. Purity of Gaussian states: Measurement schemes and time evolution in noisy channels. Phys. Rev. A 2003, 68, 012314. [Google Scholar] [CrossRef] [Green Version]
  17. Xiang, S.-H.; Shao, B.; Song, K.-H. Calculation of Covariance Matrix for Multi-mode Gaussian States in Decoherence Processes. Chin. Phys. Lett. 2009, 26, 030304. [Google Scholar]
  18. Duan, L.-M.; Giedke, G.; Cirac, J.I.; Zoller, P. Inseparability Criterion for Continuous vaiable Systems. Phys. Rev. Lett. 2000, 84, 2722. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Laurat, J.; Keller, G.; Oliveira-Huguenin, J.A.; Fabre, C.; Coudreau, T.; Serafini, A.; Adesso, G.; Illuminati, F. Entanglement of two-mode Gaussian states: Characterization and experimental production and manipulation. J. Opt. B Quantum Semiclass. Opt. 2005, 7, S577–S587. [Google Scholar] [CrossRef] [Green Version]
  20. Cariolaro, G.; Corvaja, R. Implementation of two–mode Gaussian states whose covariance matrix has the standard form. Symmetry 2022. submitted. [Google Scholar]
Figure 1. Decomposition of a Gaussian unitary in terms of fundamental Gaussian unitaries.
Figure 1. Decomposition of a Gaussian unitary in terms of fundamental Gaussian unitaries.
Symmetry 14 01286 g001
Figure 2. Decomposition of a Gaussian unitary obtained by the Bloch–Messiah reduction.
Figure 2. Decomposition of a Gaussian unitary obtained by the Bloch–Messiah reduction.
Symmetry 14 01286 g002
Figure 3. Graphical representation of the primitive components.
Figure 3. Graphical representation of the primitive components.
Symmetry 14 01286 g003
Figure 4. The beam splitter with phase.
Figure 4. The beam splitter with phase.
Symmetry 14 01286 g004
Figure 5. Implementation with primitive components of the most general two-mode Gaussian unitary, according to the Bloch–Messiah reduction and the Takagi factorization of the squeeze matrix.
Figure 5. Implementation with primitive components of the most general two-mode Gaussian unitary, according to the Bloch–Messiah reduction and the Takagi factorization of the squeeze matrix.
Symmetry 14 01286 g005
Figure 6. Generation with primitive components of a general two-mode Gaussian state, starting from two thermal states N 1 and N 2 . The two initial rotations ψ 11 and ψ 12 and the final displacements α 1 , α 2 are not introduced because they do not influence the covariance matrix.
Figure 6. Generation with primitive components of a general two-mode Gaussian state, starting from two thermal states N 1 and N 2 . The two initial rotations ψ 11 and ψ 12 and the final displacements α 1 , α 2 are not introduced because they do not influence the covariance matrix.
Symmetry 14 01286 g006
Table 1. Parameters of the architecture in Case 1. In red the free variables.
Table 1. Parameters of the architecture in Case 1. In red the free variables.
thermal states ν 1 ν 2
phases μ = 0 γ 11 = 0 γ 12 = 0 ϵ = 0
squeezers r 1 r 2
beam splitterspqsc
Table 2. Parameters of the architecture in Case 2. In red the free variables.
Table 2. Parameters of the architecture in Case 2. In red the free variables.
thermal states ν 1 ν 2
phases μ = 0 γ 11 = 0 γ 12 = 0 ϵ = 0
squeezers r 1 = r 0 r 2 = r 0
beam splitters p = 1 / 2 q = 1 / 2 sc
Table 3. Parameters of the architecture in Case 3. In red the free variables.
Table 3. Parameters of the architecture in Case 3. In red the free variables.
states ν 1 ν 2
phases μ = π / 2 γ 11 = π / 4 γ 12 = π / 4 ϵ = π / 2
squeezers r 1 = r 0 r 2 = r 0
beam splitters s = 1 / 2 c = 1 / 2 p = 1 / 2 q = 1 / 2
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cariolaro, G.; Corvaja, R.; Miatto, F. Gaussian States: Evaluation of the Covariance Matrix from the Implementation with Primitive Component. Symmetry 2022, 14, 1286. https://doi.org/10.3390/sym14071286

AMA Style

Cariolaro G, Corvaja R, Miatto F. Gaussian States: Evaluation of the Covariance Matrix from the Implementation with Primitive Component. Symmetry. 2022; 14(7):1286. https://doi.org/10.3390/sym14071286

Chicago/Turabian Style

Cariolaro, Gianfranco, Roberto Corvaja, and Filippo Miatto. 2022. "Gaussian States: Evaluation of the Covariance Matrix from the Implementation with Primitive Component" Symmetry 14, no. 7: 1286. https://doi.org/10.3390/sym14071286

APA Style

Cariolaro, G., Corvaja, R., & Miatto, F. (2022). Gaussian States: Evaluation of the Covariance Matrix from the Implementation with Primitive Component. Symmetry, 14(7), 1286. https://doi.org/10.3390/sym14071286

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop