Symmetric geometric measure and dynamics of quantum discord
Abstract
A symmetric measure of quantum correlation based on the HilbertSchmidt distance is presented in this paper. For twoqubit states, we simplify considerably the optimization procedure so that numerical evaluation can be performed efficiently. Analytical expressions for the quantum correlation are attained for some special states. We further investigate the dynamics of quantum correlation of the system qubits in the presence of independent dissipative environments. Several nontrivial aspects are demonstrated. We find that the quantum correlation can increase even if the system state is suffering dissipative noise. Sudden changes occur, even twice, in the time evolution of quantum correlation. There is certain correspondence between the evolution of quantum correlation in the systems and that in the environments, and the quantum correlation in the systems will be transferred into the environments completely and asymptotically.
pacs:
03.65.Ta, 03.67.–aI introduction
Quantum systems can be correlated in ways inaccessible to classical objects. There are quantum states that cannot be prepared with the help of local operations and classical communication, or that cannot be represented as a mixture of product states. These states are called entangled states. Entanglement is certainly a kind of quantum correlation, and moreover is by far the most famous and best studied kind of quantum correlation Horodecki et al. (2009). One reason for this situation is the fact that quantum entanglement plays an important role in much of the research of quantum information science Nielsen and Chuang (2000). However, quantum entanglement is not the only aspect of nonclassicality of correlations and does not account for all nonclassical properties of quantum correlations. A typical example is that quantum nonlocality can arise without entanglement Bennett et al. (1999). More importantly, although it is well known that entanglement is essential for certain kinds of quantuminformation tasks like teleportation Bennett et al. (1993) and superdense coding Bennett and Wiesner (1992), there is no definite answer as to whether all quantum algorithms that outperform their best known classical counterparts require entanglement as a resource. Indeed, there are several instances where we see a quantum improvement in the absence or near absence of entanglement (see e.g. Braunstein et al. (1999); Meyer (2000); Vidal (2003); Biham et al. (2004)). In particular, no entanglement is present in the computational model referred to as “the power of one qubit” with the acronym DQC1 Knill and Laflamme (1998); Laflamme et al. (2002). Despite this, the mixed separable states can create an advantage for computational tasks over their classical counterparts Laflamme et al. (2002); Datta et al. (2005); *Datta.PhysRevA.75.042310.2007.
On the other hand, many works have been devoted to understanding and quantifying the quantum correlations beyond entanglement Henderson and Vedral (2001); Ollivier and Zurek (2001); Oppenheim et al. (2002); Koashi and Winter (2004); Groisman et al. (2005); Horodecki et al. (2005); Luo (2008); RodríguezRosario et al. (2008); *Shabani.PhysRevLett.102.100402.2009; Piani et al. (2008); *Piani.PhysRevLett.102.250503.2009; Modi et al. (2010); Dakić et al. (2010); Adesso and Datta (2010); Luo and Fu (2011); Cavalcanti et al. (2011); *Madhok.PhysRevA.83.032323.2011; Streltsov et al. (2011a); *Piani.PhysRevLett.106.220403.2011. As the total correlation can be split into a classical part and a quantum part Henderson and Vedral (2001), various measures of quantum correlations are proposed by considering different notions of classicality and operational means to quantify quantumness. Amongst them, quantum discord Ollivier and Zurek (2001) has attracted much attention. Quantum discord has been defined as the mismatch between two quantum analogues of classically equivalent expression of the mutual information. It can be also expressed as the difference between total correlation, measured by quantum mutual information, and the classical correlation defined in Henderson and Vedral (2001). The notion of quantum discord goes beyond entanglement: separable states can have nonzero discord. In particular, it is believed that quantum discord is the figure of merit for the DQC1 model of quantum computation Datta et al. (2008). Generally speaking, quantum discord plays an important role in quantum information processing Merali (2011).
In order to obtain the classical correlation (and thereby the quantum discord) of a bipartite quantum system, one has to perform measurement on one subsystem to extract the information about the other subsystem, i.e., locally accessible information. Hence this sort of onesided measurements implies that classical correlation and quantum discord are not symmetric under the permutation of subsystems. As a measure of correlation, whether classical or nonclassical, one would expect that it is symmetric. Some symmetric measures of correlations have been proposed in Refs. Terhal et al. (2002); Horodecki et al. (2005); Luo (2008); Modi et al. (2010). The key step in obtaining these measures is to consider the classicalclassical (CC) states. The state is called a CC state, if it can be expressed as the mixture of locally distinguishable states, namely, , where is a joint probability distribution and local states and span an orthonormal basis. The set of product states is the subset of the set of CC states, while the set of CC states in turn is the subset of the set of separable states. A CC state is in fact the embedding of a classical probability distribution in the formalism of quantum theory and as such has no quantumness. Given a state , one can find the closest CC state, , to it. It is natural to regard this minimal distance as the measure of quantum correlation. Such a measure has a transparent geometrical meaning. The distance can be measured with trace distance, HilbertSchmidt distance, Bures distance, relative entropy and so on. In Modi et al. (2010), the authors used the relative entropy as the distance measure to provide the unified view of quantum and classical correlation.
We present in this paper a symmetric geometric measure of quantum correlation based on HilbertSchmidt (HS) distance. By applying von Neumann measurement on each subsystem, any bipartite state will become a CC state, which we call a measurementinduced CC (MICC) state. Given a bipartite state , we define the quantum correlation, denoted by , as the squared HS distance between and the closest MICC state. The evaluation of requires an optimization procedure over the set of all local von Neumann measurements on two subsystems, and thus attacking the general case is a formidable task. However, for 2qubit states, we are able to simplify the optimization over twosided measurements to that over onesided measurements. In other words, the number of measurement parameters over which the optimization procedure is performed is reduced from four to two. Thus we are able to evaluate the geometric measure of quantum correlation by the efficient numerical method. We further find the exact analytical expressions for some special states.
Moveover, we will discuss the dynamics of quantum correlation, quantified by the geometric measure , when quantum states undergo a noisy channel. It has been shown that the evolution of quantum correlation, whether in terms of quantum entanglement or of quantum discord, may behave in a “sudden” way. Quantum entanglement can evolve to sudden death or birth Yu and Eberly (2004); *Eberly.Sci.316.555.2007; *Yu.Sci.323.598.2009. The decoherence regime of quantum discord may change suddenly to that of classical correlation Mazzola et al. (2010). Also, the dynamics of quantum discord has attracted much attention Maziero et al. (2009); *Maziero.PhysRevA.81.022116.2010; Werlang et al. (2009); Wang et al. (2010); *Fanchini.PhysRevA.81.052107.2010. This “sudden” behavior has not been thoroughly understood yet. It is natural to enquire as to whether and how such a geometric measure of quantum correlation will exhibit some sort of sudden changes. To this end, we consider the case that a 2qubit state is effected by the action of two independent nonunital channels, e.g., amplitude damping (AD) channel, and calculate analytically , the time evolution of the geometric measure developed in this paper. We observe the nontrivial phenomenon that the function may be neither monotone nor smooth. In other words, during the time evolution, the quantum correlation can increase under the influence of AD channel, and the rate of evolution can exhibit sudden changes at some critical times. Despite these novel aspects, the asymptotical behavior of is in accordance with the coherence decay. In addition, we see that the quantum correlation in the system qubits is completely transferred to the environments after sufficiently long time.
Ii Definitions and notations
This section is a prelude providing the definitions and notations that will be used throughout the whole paper.
ii.1 HilbertSchmidt distance
For any 2qubit state shared by two parties and , let’s define a matrix as with the elements given by for , where is identity matrix and () are usual Pauli matrices. We write as
(1) 
where and are the Bloch vectors (in row form) of reduced state and respectively, and is a matrix and usually called correlation matrix. The superscript denotes transposition.
HilbertSchmidt (HS) norm (also called Frobenius norm) is defined as for any bounded operator on complexdimensional Hilbert space. It follows that the HS distance between two operators and is given by
(2) 
In the following, we will omit the subscript .
Let and be the matrix associated with 2qubit states and respectively. The HS distance between and can be expressed in terms of the elements of matrix:
(3) 
Now suppose that we perform local von Neumann measurements on both qubits and . The measurement operators for and are given by
respectively, where and are unit vectors in threedimensional real space. After measurements, we obtain a MICC state , that is,
(4) 
Let and . Both of them are real symmetric matrices. The matrix of can be written as
(5) 
From (3), the squared distance between and is given by
where and so forth, and in the last line we have used and . Let and . We rewrite as
(6) 
If only one qubit, say , is measured, the resulting state is a classicalquantum (CQ) state. A general CQ state has the form of . In the case we are considering, the measurementinduced CQ state is expressed as
(7) 
The corresponding matrix is
(8) 
It follows that the squared distance between and is given by
(9) 
If only qubit is measured, the squared distance can be defined similarly.
ii.2 Quantum channel
We will give only a brief introduction to the quantum channel, and refer the reader to books, say, Nielsen and Chuang (2000); Bengtsson and Życzkowski (2006) for a detailed analysis.
A quantum channel is a trace preserving completely positive (TPCP) map. A linear map is completely positive if and only if it is of the form
(10) 
where are the operators on the state space of the system and are known as Kraus operator. Eq. (10) is called operatorsum representation of a quantum operation. A completely positive map is called tracepreserving if and only if
(11) 
If a quantum channel leaves the maximal mixed state invariant, it is called a bistochastic map or unital channel, that is, . Otherwise, the channel is nonunital.
The action of a quantum channel on a quantum system can be described by a unitary evolution on an extended system, the system plus the environment , followed by a partial trace over . This description is called the unitary representation of the quantum channel. This representation is not unique since many different unitary evolutions will lead to the same effect of the channel.
Now let us consider a concrete quantum channel, amplitude damping (AD) channel. The AD channel is used to describe the evolution of the system state in the presence of a dissipative environment. In this process the system interacts with a thermal bath at zero temperature. This process could be described as spontaneous emission of a twostate atom (system ) coupled with the vacuum modes of the ambient electromagnetic field (environment ) which leads the atom state to the ground state (see, e.g., Breuer and Petruccione (2002) for more indepth discussion).
Considering the behavior of a twolevel atom in a mode cavity, the interaction between them gives rise to the phenomenological map
(12)  
(13) 
Here, is the ground state and is the excited state of the atom, while is the vacuum state of the cavity and describe the cavity state with only one excitation distributed over all modes. The amplitude converges to in the limit of . The parameter is usually called coupling strength. Then the atom and environment evolve as an effective 2qubit system.
From (12) and (13), we can write the Kraus operators of the AD channel:
(14) 
The evolution of the system state is then given by .
Let us further consider the case of a 2qubit system being affected by their independent dissipative environments. We assume that at time the systemplusenvironment state is described by
(15) 
where is the vacuum state of two environments.
Iii Geometric measures of quantum correlation
As stated in the Introduction, we can define the quantum correlation in a 2qubit state as the geometric distance, measured by HS norm, between and the closest MICC state . We will present below a detailed analysis in this direction.
iii.1 Onesided measure
To begin with, we discuss a simpler case in which only onesided measurement is performed on qubit . The geometric measure of quantum correlation based on such a sort of onesided measurements has been proposed in Dakić et al. (2010) and discussed in Luo and Fu (2010). We recover these results below for completeness.
Given a 2qubit state , onesided von Neumann measurement on induces the CQ state given by (7). With the squared distance given by (9), the geometric measure of quantum correlation is defined as
(18) 
where the minimization is performed over all von Neumann measurements on qubit .
For convenience, we use Dirac braket notation to express the vectors in real threedimensional space, that is, and . We will use alternatively both notations in the present work.
It can be seen from (9) that to find the minimal value of is equivalent to calculate the maximal value of . By noting that , we have that among all unit ,
(19) 
where is the largest eigenvalue of . Thus we have
(20) 
It is exactly the result given in Dakić et al. (2010). Similarly, if qubit is measured, we have
(21) 
with the largest eigenvalue of . Generally, .
In Dakić et al. (2010), the authors have used general, rather than measurementinduced, CQ state in the derivation of (20). Subsequently, it is pointed out in Luo and Fu (2010) that such a geometric measure is essentially measurementoriented, meaning that the optimal CQ state in the most general sense is indeed the sort of measurementinduced CQ state.
iii.2 Twosided measure
Now we proceed to the case of twosided measurements, which will lead us to a symmetric measure of quantum correlation. Twosided measurements will result in MICC states. The squared HS distance between the state and the corresponding MICC state has been given by (6). We define the geometric measure of quantum correlation, based on twosided measurement, as
(22) 
where the minimization is performed over all , or equivalently, over all measurement directions . As a CC state does not incorporate any quantum correlation, the twosided quantifier can be regarded as a more strict geometric measure of quantum correlation.
By referring to (6), we need only focus on the maximal value of for all allowed and . Let’s define the following two matrices:
(23)  
(24) 
where denotes the identity matrix. It follows that
(25) 
Then the problem is to calculate the maximal value of or over all unit vectors .
By noting that the eigenvalue of is equal to the sum of and the eigenvalue of , we see that the eigenvalue of is the function of , which we denote by . Let be the maximal one over all unit vector . It follows that the maximal value of is exactly the . Similarly, let be the eigenvalue of and the maximal one. The maximal value of is . Also, we have . The directions of the optimal measurement, denoted by and , are the eigenvectors of and referring to the eigenvalues and , respectively. Thus we see that the problem of twosided optimization is indeed reduced to that of the onesided optimization.
To proceed further, we will use the following lemma, which can be proved by direct calculation.
Lemma 1.
For any two vectors and (not necessarily normalized) in , the largest eigenvalue of the matrix is given by
(26) 
with and . The corresponding normalized eigenvector is
where and are the unit vector along and respectively, and is the normalization factor.
Let and . It follows from the Lemma 1 that
(27)  
(28) 
The remaining problem is to find the maximal values or . The geometric measure is then given by
(29) 
Referring to (27), we see that involves two parameters, that is, two angles indicating the direction of . It is not difficult to attain the by efficient numerical method. For some special states, the exact results are available, which will be stated as follows.
iii.2.1 States with twosided maximally mixed marginals
In this case, . We can always transform the correlation matrix into the diagonal form by local unitary operations. Eq. (27) becomes . The maximum is given by .
In this case there is no difference between the onesided measure and the twosided measure .
iii.2.2 States with onesided maximally mixed marginal
In this case, only one reduced state, say , is maximally mixed. Then . It follows from (27) that and the is given by the largest eigenvalue of the matrix . On the other hand, if is maximally mixed, we can consider with . Similar results are easily obtained.
iii.2.3 X states with the identical local purity
We call a 2qubit state the X state if the only nonzero elements in the density matrix lie along the diagonal or skew diagonal. By local unitary transformations, the entries in the matrix can always be written as
The exact expression of for any X state is in fact available. However, the derivation is too lengthy and too tedious. So we present the results referring to a restricted class of X states, that is, the X states with , i.e., with the identical local purity. The concrete analytical results are presented in Appendix.
Iv Dynamics of quantum correlation
Now we are in the position to analyze the dynamics of quantum correlation, in terms of the geometric measure developed in this paper. We focus on a nonunital channel, AD channel. It is assumed that the identical AD channel is applied on each qubit independently and simultaneously. In order for the analysis to be precise, we take the system’s initial state as the X state with the same local purity. The time evolution of is given by (16). The environment state can be attained by considering the total state and then tracing out and . We see that at any time both and are the sort of the X state with the identical local purity, which allows for an analytical computation of the geometric quantum correlation .
Let us consider the following two examples.
Suppose that two 2qubit X states and are given by, in terms of the elements in the matrix,
(30) 
and
(31) 
respectively. The coupling strength of the AD channel is set to be . We calculate and (or simply and ) with , and plot the results in Fig. 1. We observe several nontrivial aspects.
(i) The time evolution of quantum correlation, whether for system qubits or for the environments , is neither monotone nor smooth. There exist critical times on which the evolution rate exhibits a sudden change in behavior. The sudden change may even happen twice. More importantly, the quantum correlation in the system state can increase during the time evolution, even if the state is suffering a dissipative process.
(ii) By comparing with , we see that there is a symmetry between them. Roughly speaking, the evolution of from to is similar to the “inverse” evolution of , namely from to . Precisely, we have the following relationship: , , , and , where and () are the critical time for and , respectively (see Fig.1). In other words, there is a correspondence at the initial, critical and asymptotical point. A similar situation occurs in the comparison of and .
(iii) The fact that means that the initial quantum correlation between the system qubits is transferred completely to the asymptotical quantum correlation between the ancillary qubits.
Let us make some remarks about the above observations and related topics.
For a 2qubit system interacting with two independent environments, the dynamics and the sudden change behavior of quantum discord have been analyzed in Ref. Maziero et al. (2009); Mazzola et al. (2010); Maziero et al. (2010). Subsequently, these issues are discussed in Lu et al. (2010) in terms of another quantifier of quantum correlation, onesided HS distance measure Dakić et al. (2010). The results we present here are referring to the twosided (i.e., symmetric) HS distance measure. One significant phenomenon is that quantum correlation can increase under the influence of the separate independent dissipative environments.
Recently, some authors show that the quantum discord can increase under a local amplitude damping channel with a CC input state Ciccarello and Giovannetti (2011). A more general conclusion is drawn in Streltsov et al. (2011b), that is, any local channel which is nonunital and not semiclassical can in principle create quantum correlations, independently of the considered measure, out of a CC state. Our work provides a concrete evolution process, in particular, with respect to the state , in which we see a rapid increase of the quantum correlation in system qubits. It should be noted that the input state is a separable but not a CC state. Meanwhile, we also see that the quantum correlation of the environments undergoes a similar evolution: increases in the beginning and then changes suddenly to decrease.
It should be mentioned that quantum discord can increase in the nonMarkovian environments Wang et al. (2010); Fanchini et al. (2010), whereas the AD channel corresponds to a Markovian process.
Concerning the asymptotical behavior of , we see that after the second sudden change (if existing). This behavior is qualitatively identical to the decay of the skew diagonal elements of the density matrix , e.g., . It means that the asymptotical evolution of quantum correlation is closely related the coherence decay (in the case discussed above, they are qualitatively identical). Our discussion, from the viewpoint of symmetric distance quantifier, provides an evidence for the robustness of quantum discord Werlang et al. (2009) and the claim that “almost all quantum states have nonclassical correlations” Ferraro et al. (2010).
V Conclusion
In conclusion, we have introduced in this paper a symmetric geometric measure of quantum correlation. By performing twosided von Neumann measurements on bipartite state, we obtain a MICC state. The geometric measure is defined as the HS distance between the given state and the closest MICC state. For 2qubit system, we simplify the optimization procedure considerably, that is, the twosided optimization is reduced to the onesided one. Hence the numerical computation can be performed efficiently. Moreover, the analytical results are available for some special class of 2qubit states.
Using this quantifier, we have studied the dynamics of quantum correlation under the action of AD channel. We present the nontrivial aspects which may be exhibited during the time evolution: (i) the quantum correlation can increase; (ii) the quantum correlation can change suddenly, even twice. As for the environments, we see that the quantum correlation therein increases from zero at the beginning, and then evolves asymptotically to the value of the initial quantum correlation in the system qubits, after one or two sudden changes.
The geometric measure developed in this paper provides a symmetric viewpoint to study the quantum correlation. It is also a measurementoriented measure, since the minimal HS distance is referring to the MICC, rather general CC, states. Then a question follows: Is there more rigorous quantifier, which is referring to the general CC states? This question has been solved in the case of onesided geometric measure Luo and Fu (2010). However, it remains open for the twosided measure.
Acknowledgements.
This work was supported by National Nature Science Foundation of China, the CAS, and the National Fundamental Research Program 2007CB925200.Appendix A for X state with the same local purity
Let with . We assume that all is nonzero for the sake of simplicity. Referring to (27), we rewrite as
(32) 
Define the function as
(33) 
In the following, we discuss two cases separately: and .
a.1 Case of
When , we have . So both functions and depend only on . By taking derivative of with respect to , and solving the equation for , we have the following results.
If and , we have
(34) 
If and , we have
(35) 
Some remarks are needed here.
(i) The conditions for (34) and (35) come from two considerations: one is the nonnegativity of ; the other is the requirement of .
(ii) In solving the equation , it is assumed that . In fact, if , we have and . It follows that . We will see below that cannot be the maximal value of .
(iii) Eqs. (34) and (35) require that and , respectively. In fact, if , we have and . It follows from (32) that
It is a monotone increasing function for . Then when . This result can be contained in (39). Similarly for the case of . This analysis also holds for the derivation process in the next Case.
Inserting (34) and (35) into (32) respectively, we get two candidates for :
(36)  
(37) 
We have to take the end points of into consideration, i.e., and . When , we have . Let’s define
(38) 
When , we have
(39) 
Combining the above analysis, we conclude as follows.
If and , the maximal is given by
(40) 
If and , the maximal is given by
(41) 
a.2 Case of
Let’s prove the following lemma.
Lemma 2.
If , then at least one of is zero.
With given by (32), we introduce Lagrange multiplier and take partial derivative of with respect to . Three equations follow.
(42)  
(43)  
(44) 
Assume that all are nonzero. If so, we can delete in each equation. By noting that all are nonzero (as assumed at the beginning of this Subsection), we see from (42) and (43) that the Lagrange multiplier must be zero. Then Eqs. (42) or (43) reduce to . It follows that , which contradicts the assumption. Lemma 2 is proved.
Subsequently, we will discuss one by one the measurements allowed by Lemma 2. In each case, we obtain a candidate for . The largest one is what we want.
If and .
It follows that , and that both and are the functions of only. By the approach similar to that presented in Section A.1, we have the following results.
If and , we have
(45) 
If and , we have
(46) 
If and .
In this case, . The results are very similar to that presented in the last paragraph, that is:
If , we have
(47) 
If , we have
(48) 
If .
Here it is not required that both and are nonzero. Inserting into the expression of , i.e., Eq. (32), we have
The maximal value of the above expression is given by
(49) 
If .
In this case, . It easily follows that
(50) 
It is not difficult to see that the above four cases cover all allowed measurements. We summary the results obtained in this Subsection in Table 1.
Conditions  

, and  , ^{1}^{1}1, , , .  , (45) 
,  , (46)  
, and  ,  , (47) 
,  , (48)  
, (49)  
, (50) 
Now let’s use an example to show how to calculate . Given an X state with , we write its matrix. If we see that and , then we face the following possibilities:

If and , then .

If and , then .

If and , then .

If and , then