Nonclassicality by Local Gaussian Unitary Operations for Gaussian States

A measure of nonclassicality N in terms of local Gaussian unitary operations for bipartite Gaussian states is introduced. N is a faithful quantum correlation measure for Gaussian states as product states have no such correlation and every non product Gaussian state contains it. For any bipartite Gaussian state ρAB, we always have 0≤N(ρAB)<1, where the upper bound 1 is sharp. An explicit formula of N for (1+1)-mode Gaussian states and an estimate of N for (n+m)-mode Gaussian states are presented. A criterion of entanglement is established in terms of this correlation. The quantum correlation N is also compared with entanglement, Gaussian discord and Gaussian geometric discord.


Introduction
The presence of correlations in bipartite quantum systems is one of the main features of quantum mechanics. The most important one among such correlations is entanglement [1]. However, recently much attention has been devoted to the study and the characterization of quantum correlations that go beyond the paradigm of entanglement, being necessary but not sufficient for its presence. Non-entangled quantum correlations also play important roles in various quantum communications and quantum computing tasks [2][3][4][5].
For the last two decades, various methods have been proposed to quantify quantum correlations, such as quantum discord (QD) [6,7], geometric quantum discord [8,9], measurement-induced nonlocality (MIN) [10] and measurement-induced disturbance (MID) [11] for discrete-variable systems. It is also important to develop new simple criteria for witnessing correlations beyond entanglement for continuous-variable systems. In this direction, Giorda, Paris [12] and Adesso, Datta [13] independently introduced the definition of Gaussian QD for Gaussian states and discussed its properties. Adesso and Girolami in [14] proposed the concept of Gaussian geometric discord (GD) for Gaussian states. Measurement-induced disturbance of Gaussian states was studied in [15], while MIN for Gaussian states was discussed in [16]. For other related results, see [17,18] and the references therein. Note that not every quantum correlation defined for discrete-variable systems has a Gaussian analogy for continuous-variable systems [16]. On the other hand, the values of Gaussian QD and Gaussian GD are very difficult to be computed and the known formulas are only for some (1 + 1)-mode Gaussian states. Little information is revealed by Gaussian QD and GD. The purpose of this paper is to introduce a new ρ − (U A ⊗ I)ρ(U A ⊗ I) † F , where A F = Tr(A † A) denotes the Frobenius norm and U A is any unitary operator satisfying [ρ A , U A ] = 0. This quantity demands that the reduced density matrix of the subsystem A is invariant under this unitary transformation. However, the global density matrix may be changed after such local unitary operation, and therefore d U A (ρ) may be non-zero for some U A . Then, Datta, Gharibian, et al. discussed respectively in [20,21] the properties of d U A (ρ) and revealed that max U A d U A (ρ) can be used to investigate the nonclassical effect.
Motivated by the works in [19][20][21], we can consider an analogy for continuous-varable systems. In the present paper, we introduce a quantity N in terms of local Gaussian unitary operations for (n + m)-mode quantum states in Gaussian systems. Different from the finite dimensional case, besides the local Gaussian unitary invariance property for quantum states, we also show that N (ρ AB ) = 0 if and only if ρ AB is a Gaussian product state. This reveals that the quantity N is a kind of faithful measure of the nonclassicality for Gaussian states that a state has this nonclassicality if and only if it is not a product state. In addition, we show that 0 ≤ N (ρ AB ) < 1 for each (n + m)-mode Gaussian state ρ AB and the upper bound 1 is sharp. An estimate of N for any (n + m)-mode Gaussian states is provided and an explicit formula of N for any (1 + 1)-mode Gaussian states is obtained. As an application, a criterion of entanglement for (1 + 1)-mode Gaussian states is established in terms of N by numerical approaches. Finally, we compare N with Gaussian QD and Gaussian GD to illustrate that it is a better measure of the nonclassicality.

Gaussian States and Gaussian Unitary Operations
Recall that, for arbitrary state ρ in an n-mode continuous-variable system, its characteristic function χ ρ is defined as where z = (x 1 , y 1 , · · · , x n , y n ) T ∈ R 2n with R the field of real numbers and (·) T the transposition, and W(z) = exp(iR T z) is the Weyl operator. Let R = (R 1 , R 2 , · · · , R 2n ) T = (Q 1 ,P 1 , · · · ,Q n ,P n ) T . As usual,Q i andP i stand respectively for the position and momentum operators for each i ∈ {1, 2, · · · , n}. They satisfy the Canonical Commutation Relation (CCR) in natural units (h = 1) Gaussian states: ρ is called a Gaussian state if χ ρ (z) is of the form where d = ( R 1 , R 2 , . . . , R 2n ) T = (Tr(ρR 1 ), Tr(ρR 2 ), . . . , Tr(ρR 2n )) T ∈ R 2n is called the mean or the displacement vector of ρ and Γ = ( [22][23][24]). Here, M l×k (R) stands for the set of all l-by-k real matrices and, when l = k, we write M l×k (R) as M l (R). Note that the CM Γ of a state is symmetric and must satisfy the uncertainty principle Γ + i∆ ≥ 0, where From the diagonal terms of the above inequality, one can easily derive the usual Heisenberg uncertainty relation for position and momentum V(Q i )V(P i ) ≥ 1 with V(R i ) = (∆R i ) 2 [25]. Now assume that ρ AB is any (n + m)-mode Gaussian state. Then, the CM Γ of ρ AB can be written as where A ∈ M 2n (R), B ∈ M 2m (R) and C ∈ M 2n×2m (R). Particularly, if n = m = 1 , by means of local Gaussian unitary (symplectic at the CM level) operations, Γ has a standard form: [26][27][28][29]). Gaussian unitary operations. Let us consider an n-mode continuous-variable system with R = (Q 1 ,P 1 , · · · ,Q n ,P n ) T . For a unitary operator U, the unitary operation ρ → UρU † is said to be Gaussian if its output is a Gaussian state whenever its input is a Gaussian state, and such U is called a Gaussian unitary operator. It is known that a unitary operator U is Gaussian if and only if for some vector m in R 2n and some S ∈ Sp(2n, R), the symplectic group of all 2n × 2n real matrices S that satisfy Thus, every Gaussian unitary operator U is determined by some affine symplectic map (S, m) acting on the phase space, and can be denoted by U = U S,m ( [23,24]).
The following well-known facts for Gaussian states and Gaussian unitary operations are useful for our purpose.

Lemma 1 ([23]
). For any (n + m)-mode Gaussian state ρ AB , write its CM Γ as in Equation (1). Then, the CMs of the reduced states ρ A = Tr B ρ AB and ρ B = Tr A ρ AB are matrices A and B, respectively.
where Γ, Γ A and Γ B are the CMs of ρ AB , σ A and σ B , respectively. Lemma 3 ([23,24]). Assume that ρ is any n-mode Gaussian state with CM Γ and displacement vector d, and U S,m is a Gaussian unitary operator. Then, the characteristic function of the Gaussian state σ = UρU † is of the form exp(− 1 4

Quantum Correlation Introduced by Gaussian Unitary Operations
Now, we introduce a quantum correlation N by local Gaussian unitary operations in the continuous-variable system.

Definition 1. For any (n + m)-mode quantum state ρ AB ∈ S(H
where the supremum is taken over all Gaussian unitary operators U ∈ B(H B ) satisfying Uρ B U † = ρ B , and ρ B = Tr A (ρ AB ) is the reduced state. Here, B(H B ) is the set of all bounded linear operators acting on H B .
Observe that N (ρ AB ) = 0 holds for every product state. Thus, the product state contains no such correlation.

Remark 1.
For any Gaussian state ρ AB , there exist many Gaussian unitary U so that Uρ B U † = ρ B . This ensures that the definition of the quantity N (ρ AB ) makes sense for each Gaussian state ρ AB .
To see this, we need Williamson Theorem ( [31]), which states that, for any n-mode Gaussian state and v i s are called respectively the Williamson form and the symplectic eigenvalues of Γ ρ . By the Williamson Theorem, there exists a Gaussian unitary operator . Then, S θ is a symplectic matrix, and the corresponding We first prove that N is local Gaussian unitary invariant for all quantum states.

Proposition 1 (Local Gaussian unitary invariance). If ρ AB ∈ S(H
Then, W is also a Gaussian unitary operator and satisfies The next theorem shows that N (ρ AB ) is a faithful nonclassicality measure for Gaussian states.
Proof of Theorem 1. By Definition 1, the "if" part is apparent. Let us check the "only if" part. Since the mean of any Gaussian state can be transformed to zero under some local Gaussian unitary operation, it is sufficient to consider those Gaussian states whose means are zero by Proposition 1. In the sequel, assume that ρ AB is an (n + m)-mode Gaussian state with zero mean vector and CM Γ = A C C T B as in Equation (1), so that N (ρ AB ) = 0.
By Lemma 1, the CM of ρ B is B. According to the Williamson Theorem, there exists a symplectic matrix S 0 such that . It follows from Proposition 1 that N (σ AB ) = N (ρ AB ) = 0. Obviously, σ AB has the CM of form: For any θ i ∈ [0, π 2 ] for i = 1, 2, · · · , m, let S θ be the symplectic matrix as in Remark 1. Then, , and hence they must have the same CMs, that is, Note that I − S T θ is an invertible matrix if we take θ i ∈ (0, π 2 ) for each i. Then, it follows from C = C S T θ that we must have C = 0. Thus, σ AB is a product state by Lemma 2, and, consequently, is also a product state.
We can give an analytic formula of N (ρ AB ) for (1+1)-mode Gaussian state ρ AB . Since N is locally Gaussian unitary invariant, it is enough to assume that the mean vector of ρ AB is zero and the CM is standard.
Proof of Theorem 2. By Proposition 1, we may assume that the mean vector of ρ AB is zero. Let U S,m be a Gaussian unitary operator such that U S,m ρ B U † S,m = ρ B . Then, S and m meet the conditions SB 0 S T = B 0 and Sd B + m = d B = 0. It follows that m = 0. Thus, we can denote U S,m by U S .
For the general (n + m)-mode case, it is difficult to give an analytic formula of N (ρ AB ) for all (n + m)-mode Gaussian states ρ AB . However, we are able to give an estimate of N (ρ AB ).

Theorem 3. For any (n + m)-mode Gaussian state ρ AB with CM Γ =
A C C T B as in Equation (1), Particularly, when ρ AB is pure, . Moreover, the upper bound 1 in the inequality Proof of Theorem 3. By Proposition 1, without loss of generality, we may assume that the mean of Note that, for any n-mode Gaussian states ρ, σ with CMs V ρ , V σ and means d ρ , d σ , respectively, it is shown in [32] that Hence, , by Fischer's inequality (p. 506, [33]), we have . If ρ AB is a pure state, then Notice that, by Equation (6), we have 1 < 1 since det A > 0 and det B > 0, that is, the inequality (5) is true.
To see that the upper bound 1 is sharp, consider the two-mode squeezed vacuum state 2 ) is the two-mode squeezing operator with squeezed number r ≥ 0 and |00 is the vacuum state ( [24]). The CM . By Theorem 2, it is easily calculated that N (ρ(r)) = 1 − 8 6 + exp(−4r) + exp(4r) .

Comparison with Other Quantum Correlations
Entanglement is one of the most important quantum correlations, being central in most quantum information protocols [1]. However, it is an extremely difficult task to verify whether a given quantum state is entangled or not. Recall that a quantum state ρ AB ∈ S(H A ⊗ H B ) is said to be separable if it belongs to the closed convex hull of the set of all product states ρ A ⊗ ρ B ∈ S(H A ⊗ H B ). Note that a state ρ AB is separable if and only if it admits a representation ρ AB = X ρ A (x) ⊗ ρ B (x)π(dx), where π(dx) is a Borel probability measure and ρ A(B) (x) is a Borel S(H A(B) )-valued function on some complete, separable metric space X [34]. One of the most useful separability criteria is the positive partial transpose (PPT) criterion, which can be found in [35,36]. The PPT criterion states that if a state is separable, then its partial transposition is positive. For discrete systems, the positivity of the partial transposition of a state is necessary and sufficient for its separability in the 2 ⊗ 2 and 2 ⊗ 3 cases. However, it is not true for higher dimensional systems [36]. For continuous systems, in [27,37], the authors extended the PPT criterion to (n + m) -mode continuous systems. It is remarkable that, for any (1 + n)-mode Gaussian state, it has PPT if and only if it is separable. Furthermore, for the (1 + 1)-mode case, it is shown that a (1 + 1)-mode Gaussian state ρ AB is separable if and only ifv − ≥ 1, wherev − is the smallest symplectic eigenvalue of the CM of the partial transpose ρ T B AB [24,29]. Comparing N with the entanglement, we conjecture that there exists some positive number d < 1 such that N (ρ AB ) ≤ d for any (n + m)-mode separable Gaussian state ρ AB , that is, If this is true, then ρ AB is entangled when N (ρ AB ) > d. This will give a criterion of entanglement for (n + m)-mode Gaussian states in terms of correlation N . Though we can not give a mathematical proof, we show that this is true for (1 + 1)-mode separable Gaussian states with d ≤ 1 10 by a numerical approach (Firstly, we randomly generated one million, five million, ten million, fifty million, one hundred million, five hundred million separable Gaussian states with a, b, |c|, |d| ranging from 1 to 2, respectively. We found that the maximum of N is smaller than 0.09. Secondly, we used the same method and extended the range to 5. Then, the maximum of N is smaller than 0.1. Thirdly, using the same method and extending the range to 10, 100, 1000, 10000, respectively, we found that the maximum of N is still smaller than 0.1. We repeated the above computations ten times, and the result is just the same). Proposition 2. N (ρ AB ) ≤ 0.1 for any (1 + 1)-mode separable Gaussian state ρ AB .
It is followed from Theorem 1 that the quantum correlation N exists in all entangled Gaussian states and almost all separable Gaussian states except product states. In addition, Proposition 2 can be viewed as a sufficient condition for the entanglement of two-mode Gaussian states: if N (ρ AB ) > 0.1, then ρ AB is entangled.
To have an insight into the behavior of this quantum correlation by N and to compare it with the entanglement and the discords, we consider a class of physically relevant states-squeezed thermal state (STS). This kind of Gaussian state is used by many authors to illustrate the behavior of several interesting quantum correlations [12,13]. Recall that a two-mode Gaussian state ρ AB is an STS if ρ AB = S(r)ν 1 (n 1 ) ⊗ ν 2 (n 2 )S(r) † , where ν i (n i ) = ∑ kn k i (1+n i ) k+1 |k k| is the thermal state with thermal photon numbern i (i = 1, 2) and S(r) = exp{r(â † 1â † 2 −â 1â2 )} is the two-mode squeezing operator. Particularly, whenn 1 =n 2 = 0, ρ AB is a pure two-mode squeezed vacuum state, also known as an Einstein-Podolski-Rosen (EPR) state [24]. Whenn 1 > 0 orn 2 > 0, ρ AB is a mixed Gaussian state.
We first discuss the relation between N and the entanglement by considering SSTS. Regard N (ρ AB ) as a function of µ andv − . From Figure 1a, for separable states, we see that the value N at the separable SSTS is always smaller than 0.06, which supports positively Proposition 2. From Figure 1b, for fixed purity µ, N turns out to be a decreasing function ofv − . However, for fixedv − , N tends to 0 when µ increases. For the entangled SSTS, one sees from Figure 2a,b that the value of N is from 0 to 1. This reveals that, for some entangled SSTSs, N can be smaller than 1 10 . Thus, Proposition 2 is only a necessary condition for a Gaussian state to be separable. For fixed purity µ, from Figure 1b and 2b, N (ρ AB ) increases when entanglement increases (that is,v − → 0) and lim µ→1,v − →0 N = 1. However, for fixed v − , the behavior of N on µ is more complex. Regarding N as a function of r andn, Figure 3 shows that N (ρ AB ) is an increasing function of r and a decreasing function ofn, respectively. The value of N (ρ AB ) always gains the maximum at n = 0, that is, at pure states. Figure 3b also shows that N (ρ AB ) almost depends only onn when r is large enough because the curves for r = 5, 10, 20 are almost the same. Recall that an n-mode Gaussian positive operator-valued measure (GPOVM) is a collection of positive operators Π = {Π(z)} satisfying z Π(z)dz = I, where Π(z) = W(z)ωW † (z), z ∈ R 2n with W(z) the Weyl operators and ω an n-mode Gaussian state, which is called the seed of the GPOVM Π [38,39]. Let ρ AB be a (n + m)-mode Gaussian state and Π = {Π(z)} be a GPOVM of the subsystem B. Denote by ρ A (z) = 1 p(z) Tr B (ρ AB I ⊗ Π(z)) the reduced state of the system A after the GPOVM Π performed on the system B, where p(z) = Tr(ρ AB I ⊗ Π(z)). Write the von Neumann entropy of a state ρ as S(ρ), that is, S(ρ) = −Tr(ρ log ρ). Then, the Gaussian QD of ρ AB is defined as D(ρ AB ) = S(ρ B ) − S(ρ AB ) + inf Π dzp(z)S(ρ A (z)) [12,13], where the infimum takes over all GPOVMs Π performed on the system B. It is known that a (1 + 1)-mode Gaussian state has zero Gaussian QD if and only if it is a product state; in addition, for all separable (1 + 1)-mode Gaussian states, D(ρ AB ) ≤ 1; if the standard form of the CM of a (1 + 1)-mode Gaussian state ρ AB is as in Equation (2), then where the infimum takes over all one-mode Gaussian states ω, f (x) = x+1 2 log x+1 2 − x−1 2 log x−1 2 , v − and v + are the symplectic eigenvalues of the CM of ρ AB , E ω = A 0 − C 0 (B 0 + Γ ω ) −1 C T 0 with Γ ω the CM of ω. Let α = det A 0 , β = det B 0 , γ = det C 0 , δ = det Γ 0 , then we have [13] inf otherwise.
In [14], the quantum GD D G is proposed. Consider an (n + m)-mode Gaussian state ρ AB , its Gaussian GD is defined by D G (ρ AB ) = inf Π ||ρ AB − Π(ρ AB )|| 2 2 , where the infimum takes over all GPOVM Π performed on system B, ||· || 2 stands for the Hilbert-Schmidt norm and Π(ρ AB ) = dz(I ⊗ Π(z))ρ AB (I ⊗ Π(z)). If ρ AB is a (1 + 1)-mode Gaussian state with the CM Γ as in Equation (1) and Π is an one-mode Gaussian POVM performed on mode B with seed ω B , then It is known from [14] that Now it is clear that, for (1 + 1)-mode Gaussian state ρ AB , D G (ρ AB ) = 0 if and only if ρ AB is a product state. By Theorem 1 and the results mentioned above, D, D G and N describe the same quantum correlation for (1 + 1)-mode Gaussian states. However, from the definitions, D, D G use all GPOVMs, while N only employs Gaussian unitary operations, which is simpler and may consume less physical resources. Moreover, though an analytical formula of D is given for two-mode Gaussian states, the expression is more complex and more difficult to calculate (Equations (8) and (9)). D G is not handled in general and there is no analytical formula for all (1 + 1)-mode Gaussian states (Equation (10)). As far as we know, there are no results obtained on D, D G for general (n + m)-mode case.
To have a better insight into the behavior of N and D G , we compare them in scale with the help of two-mode STS. Note that D G of any two-mode STS ρ AB is given by [14] Clearly, our formula (7) for N is simpler then formula (11) for D G . Figures 4 and 5 are plotted in terms of photo numbern and squeezing parameter r. Figure 4 shows that, for the case of SSTS and for 0 < r ≤ 2.5, we have D G (ρ AB ) < N (ρ AB ). This means that N is better than D G when they are used to detect the correlation that they describe in the SSTS with r < 2.5. Figure 5a reveals that, for the case of nonsymmetric STS and for r = 0.5, we have D G (ρ AB ) < N (ρ AB ); that is, N is better in this situation too. However, for r = 5, N and D G can not be compared with each other globally, which suggests that one may use max{N (ρ AB ), D G (ρ AB )} to detect the correlation.

Conclusions
In conclusion, we introduce a measure of quantum correlation by N for bipartite quantum states in continuous-variable systems. This measure is introduced by performing Gaussian unitary operations to a subsystem and the value of it is invariant for all quantum states under local Gaussian unitary operations. N exists in all (n + m)-mode Gaussian states except product ones. In addition, N takes values in [0, 1) and the upper bound 1 is sharp. An analytical formula of N for any (1 + 1)-mode Gaussian states is obtained. Moreover, for any (n + m)-mode Gaussian states, an estimate of N is established in terms of its covariance matrix. Numerical evidence shows that the inequality N (ρ AB ) ≤ 0.1 holds for any (1 + 1)-mode separable Gaussian states ρ AB , which can be viewed as a criterion of entanglement. It is worth noting that Gaussian QD, Gaussian GD and N measure the same quantum correlation for (1 + 1)-mode Gaussian states. However, N is easer to calculate and can be applied to any (n + m)-mode Gaussian states.