We study the time evolution of quantum discord and entanglement of a two-qubit system in thermal reservoirs. We find that there are no simple ordering relations between entanglement and quantum discord for the dynamical evolution behavior; that is, quantum discord may be smaller or larger than entanglement in the evolution process. However, a strong correlation between changes of the ordering relations and the mean photon number is found. In addition, it also shows that entanglement is not the strongest form of nonclassicality.
Considerable efforts have been devoted to study quantum correlations in the last two decades, mainly due to the general belief that they are a fundamental resource for quantum information processing tasks. Among them, entanglement is the most famous and best studied kind of quantum correlation and is generally regarded as a necessary resource in quantum computation and communication [1]. The situation started to change after a so-called DQC1 (deterministic quantum computation with one qubit) model was presented which may provide the speedup in the deterministic quantum computation with one pure qubit [2]. The fact is that no entanglement is present in this model; however, other kinds of nonclassical correlations are responsible for the quantum computational efficiency of DQC1. Such correlations are characterized as quantum discord [3], which is believed to be more general and more fundamental than entanglement [4–6] and which can be viewed as the amount of entanglement consumed in the quantum-state merging [7, 8].
Besides its application in DQC1, quantum discord has also been used in studies of quantum phase transition [9, 10], Maxwell’s demon [11], and relativistic effect [12, 13]. In addition, concerning biological systems, it has been suggested that such correlations could play an important role in photosynthesis [14]. In particular, quantum discord has been predicted to show peculiar dynamics under decoherence [15]. Considering two noninteracting qubits, it was shown that, under certain conditions, the decay rate of the correlations may suffer a sudden change [16, 17]. Furthermore, by analyzing various dissipative channels, some authors have shown that, in all cases where entanglement suddenly disappears, quantum discord vanishes only asymptotically [18–20]. In this sense, quantum discord is more robust against decoherence than entanglement so that quantum algorithms based only on quantum discord correlations may be more robust than those based on entanglement. Despite intense research over the last decade, the relation between entanglement and quantum discord is not clear, and there are still a number of challenging outstanding problems that drive much effort in the field.
In this paper, we consider a two-qubit system interacting with two independent thermal reservoirs. We will explore the relation between the dynamics of entanglement and that of quantum discord of two atoms in their local reservoirs. The common belief is that measures of quantum discord should be more than those of entanglement. However, our results will show that quantum discord may be smaller or larger than entanglement in the evolution process, and there exists an interesting variation in the ordering relations between entanglement and quantum discord. Moreover, it is also shown that entanglement is not the strongest form of nonclassicality, as pointed out in [21].
For a two-qubit quantum system, the quantum discord is defined as the difference between the quantum mutual information and the classical correlation,
(1)D(ρab)=I(ρab)-C′(ρab).
Here, the quantum mutual information of two subsystems is given by
(2)I(ρab)=S(ρa)+S(ρb)-S(ρab),
and S(ρj) is the entropy of the quantum system, which is given by the von Neumann entropy,
(3)S(ρj)=-Trj(ρjlog2ρj)=-∑iλjilog2λji,
where the subscript j indicates either the subsystem a(b) or the total system and {λji} are the nonzero eigenvalues of the quantum state ρj. C′(ρab) is the classical correlation between the two subsystems [4]. As discussed in [15], the classical correlation is given by
(4)C′(ρab)=S(ρa)-min{Πk}[S(ρab∣{Πk})],
where the minimum is taken over the set of projective measurements {∏k} and S(ρab∣{∏k})=∑kpkS(ρk) is the quantum conditional entropy of a given the complete measurements on subsystem b, with ρk=Trb(∏kρab∏k)/pk and pk=Trab(ρab∏k).
Then, quantum discord is given by
(5)D(ρab)=S(ρb)-S(ρab)+min{Πk}S(ρab∣{Πk}).
The evaluation of quantum discord is a hard task due to the minimization over all possible measurements. Although this cannot be done in full generality analytically even for two qubits, the problem is tractable using analytical tools for certain subsets of states (such as X states) [22, 23].
The density matrix of a two-qubit X state written in the standard basis is given by
(6)ρ(t)=(ρ1100ρ140ρ22ρ230ρ32ρ330ρ4100ρ44).
The eigenvalues of the matrix ρ(t) in (6) can be obtained as follows:
(7)λ0,1=12[(ρ11+ρ44)±(ρ11-ρ44)2+4|ρ14|2],λ2,3=12[(ρ22+ρ33)±(ρ22-ρ33)2+4|ρ23|2],
and the von Neumann entropies of the reduced density matrix of the ρ(t) can be easily obtained as
(8)S(ρa)=-[(ρ11+ρ22)log2(ρ11+ρ22)+(ρ33+ρ44)log2(ρ33+ρ44)],S(ρb)=-[(ρ11+ρ33)log2(ρ11+ρ33)+(ρ22+ρ44)log2(ρ22+ρ44)].
The quantum conditional entropy in (4) is given by
(9)S(ρab∣{Πk})=p0S(ρ0)+p1S(ρ1),
where
(10)p0=[(ρ11+ρ33)k+(ρ22+ρ44)l],p1=[(ρ11+ρ33)l+(ρ22+ρ44)k],
and the entropies are given by
(11)S(ρ0)=-1-θ2log21-θ2-1+θ2log21+θ2,S(ρ1)=-1-θ′2log21-θ′2-1+θ′2log21+θ′2,
with
(12)θ=[(ρ11-ρ33)k+(ρ22-ρ44)l]2+Θ[(ρ11+ρ33)k+(ρ22+ρ44)l]2,θ′=[(ρ11-ρ33)l+(ρ22-ρ44)k]2+Θ[(ρ11+ρ33)l+(ρ22+ρ44)k]2,
where Θ=4kl[|ρ14|2+|ρ23|2+2Re(ρ14ρ23)]-16mRe(ρ14ρ23)+16nIm(ρ14ρ23) and Re(z) and Im(z) are the real and imaginary parts of the complex number z.
The minimum value of S(ρab∣{Πk}) in (4) depends on parameters k, l, m, and n, where k∈[0,1], m∈[0,1/4], n∈[-1/8,1/8], and k+l=1 [22].
To investigate entanglement dynamics, we use the concurrence as the measure [24]. For a pair of qubits, the concurrence is defined as
(13)C(ρ)=max{0,λ1′-λ2′-λ3′-λ4′},
where the λi′ are the eigenvalues (λ1′ being the largest one) of non-Hermitian matrix ρρ~ and ρ~ is defined as
(14)ρ~=(σy⊗σy)ρ*(σy⊗σy),
with ρ* being the complex conjugate of matrix ρ. Here, σy stands for the off-diagonal pure imaginary Pauli matrix. The concurrence C=0 corresponds to an unentangled state and C=1 corresponds to a maximally entangled state.
In the case of the density matrix with the X type (6), one easily finds [25] that the concurrence is given by
(15)C=max{0,Λ1,Λ2},
where Λ1=2(ρ23ρ32-ρ11ρ44), Λ2=2(ρ14ρ41-ρ22ρ33).
Next, we investigate the dynamics of quantum discord and entanglement of the two atoms which are each trapped separately in one of the two spatially separated thermal reservoirs. According to the general quantum reservoir theory [26], we can derive the following equation of motion for the reduced density matrix of the atoms interacting with their local thermal reservoirs of mean thermal photon numbers m and n:
(16)dρdt=-12γ1(m+1)[σ1+σ1ρ-2σ1ρσ1++ρσ1+σ1]-12γ1m[σ1σ1+ρ-2σ1+ρσ1+ρσ1σ1+]-12γ2(n+1)[σ2+σ2ρ-2σ2ρσ2++ρσ2+σ2]-12γ2n[σ2σ2+ρ-2σ2+ρσ2+ρσ2σ2+],
where γi is the spontaneous emission rate of atom i and σi+, σi are the usual Pauli raising and lowering operators of atom i. Here, we assume that the correlation time between the atoms and the reservoirs is much shorter than the characteristic time of the dynamic evolution of the atoms so that the Markov approximation is valid.
In what follows, we will focus on the initial Bell states since such states provide important test ground for the understanding of nonclassicality beyond entanglement and its interplay with entanglement itself.
Now, we would like to consider Bell state |Φ+〉=(1/2)(|10〉+|01〉) as the initial state. For simplicity, we assume m=n and γ1=γ2=γ. The solution of master equation (16) is given by
(17)ρ11=12(2n+1)2{2n2+2ne-(2n+1)γt-2n(n+1)e-2(2n+1)γt},ρ22=ρ33=12(2n+1)2{2n(n+1)+e-(2n+1)γt+2n(n+1)e-2(2n+1)γt},ρ44=12(2n+1)2{2(n+1)2-2(n+1)e-(2n+1)γt-2n(n+1)e-2(2n+1)γt},ρ23=ρ32=12e-(2n+1)γt.
Using (17), we can give the eigenvalues of the matrix discussed:
(18)λ0,1=12[(ρ11+ρ44)±|ρ11-ρ44|],λ2,3=ρ22±ρ23.
Based on the discussion of the above, we can now calculate the classical correlation and quantum discord. The quantum discord is given by
(19)D(ρab)=S(ρb)-S(ρab)+min{Πk}{S(ρ0)|θ1,S(ρab∣{Πk})|θ2,θ3},
where
(20)S(ρ0)|θ1=-1+θ12log1+θ12-1-θ12log21-θ12,S(ρab∣{Πk})|θ2,θ3=(ρ22+ρ44)(-1+θ22log21+θ22-1-θ22log21-θ22)+(ρ11+ρ33)(-1+θ32log21+θ32-1-θ32log21-θ32),
with θ1=(ρ11-ρ44)2+4|ρ23|2, θ2=(|ρ22-ρ44|/(ρ22+ρ44)), and θ3=(|ρ11-ρ33|/(ρ11+ρ33)).
Using (15) and (19), the expressions of the entanglement and quantum discord can be numerically calculated. In Figure 1, the dynamics of entanglement and quantum discord are plotted as a function of the scaled time γt and mean photon number n. Firstly, we note that quantum discord only disappears in the asymptotic limit, while the entanglement suffers a sudden death at a certain parameterized time td in a thermal reservoir of nonzero mean photon number. Considering the fact that the entanglement sudden death (ESD) was firstly reported in [27], thereafter, many works have been devoted to the related topics [28–31]. We also see that both time of ESD and vanishing of quantum discord decrease as the mean thermal photon number becomes large, as expected.
Dynamics of the quantum discord (a) and entanglement (b) as a function of the scaled time γt and mean photon number n.
In order to explore the relation between entanglement and quantum discord, we plot the evolution of both correlations for different values of the mean photon numbers n=1, 0.1, 0.01, 0 in Figure 2. The abscissa of the vertical dotted line corresponds to the scaled time γt′ of the intersection point of the two evolution curves. From Figures 2(a)–2(c), one can see a similar phenomenon that the two correlations have a different order as functions of time γt, with concurrence being initially larger than quantum discord for 0≤γt≤γt′, but, for γt≥γt′, quantum discord becomes larger than concurrence. It means that entanglement is not the strongest form of nonclassicality. Although there are no simple ordering relations between the two correlations, it is instructive to note that the position of the scaled time γt′ of the intersection point of the two evolution curves moves to the right along the γt-axis with the decreases of the mean photon number n. Moreover, one can observe from Figure 2(d) that, for n=0, γt′→∞. It can also be found that the dynamical behavior of the entanglement tends to agree with that of the quantum discord in time in the vacuum reservoirs. Meanwhile, it also implies that the ESD vanishes as the mean thermal photon number is zero. These results show that further and deeper study of the dynamical behaviors of quantum discord and entanglement in thermal environments might help to better understand the nature of quantum correlations.
Time dependence of quantum discord (dot-dashed line) and concurrence (solid line) for mean photon number: (a) n=1, (b) n=0.1, (c) n=0.01, and (d) n=0 (color online). The abscissa of the vertical dotted line corresponds to the scaled time γt′ of the intersection point of the two evolution curves.
In summary, we have investigated the dynamics of the entanglement and the quantum discord for two atoms coupled to independent thermal reservoirs. We found that there are no simple ordering relations between entanglement and quantum discord for the dynamical evolution behavior; however, the position of the intersection point of the two correlation curves moves to the right along the γt-axis as the mean photon number becomes small. We also observed that, in thermal environments where entanglement suddenly disappears, quantum discord vanishes only in the asymptotic limit and the dynamical behavior of quantum discord tends to agree with that of entanglement in time in the vacuum reservoir. Our results also show that entanglement is not the strongest form of nonclassicality. It suggests that quantum algorithms based on quantum discord, different from those based on entanglement, are more resistant to external perturbations and, therefore, introduce new hope of implementing an efficient quantum algorithm and quantum information processing.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
The authors are grateful to Dr. Z. Qi for careful reading of the original paper.
NielsenM. A.ChuangI. L.2000Cambridge, UKCambridge University PressKnillE.LaflammeR.Power of one bit of quantum information19988125567256752-s2.0-11544307709OllivierH.ZurekW. H.Quantum discord: a measure of the quantumness of correlations20028810179012-s2.0-0037033316CéleriL. C.MazieroJ.SerraR. M.Theoretical and experimental aspects of quantum discord and related measures201197-8183718732-s2.0-8305516725910.1142/S0219749911008374TathamR.MištaL.AdessoG.KorolkovaN.Nonclassical correlations in continuous-variable non-Gaussian Werner states20128522-s2.0-8485754504810.1103/PhysRevA.85.022326022326ModiK.BrodutchA.CableH.PaterekT.VedralV.The classical-quantum boundary for correlations: discord and related measures20128441655170710.1103/RevModPhys.84.1655CavalcantiD.AolitaL.BoixoS.ModiK.PianiM.WinterA.Operational interpretations of quantum discord20118332-s2.0-7995795967310.1103/PhysRevA.83.032324032324MadhokV.DattaA.Interpreting quantum discord through quantum state merging20118332-s2.0-7996063741510.1103/PhysRevA.83.032323032323SarandyM. S.Classical correlation and quantum discord in critical systems20098022-s2.0-6894913543310.1103/PhysRevA.80.022108022108WerlangT.TrippeC.RibeiroG. A. P.RigolinG.Quantum correlations in spin chains at finite temperatures and quantum phase transitions201010592-s2.0-7795628273710.1103/PhysRevLett.105.095702095702ZurekW. H.Quantum discord and Maxwell's demons20036710123202-s2.0-0037240898DattaA.Quantum discord between relatively accelerated observers20098052-s2.0-7044909650010.1103/PhysRevA.80.052304052304WangJ. C.DengJ. F.JingJ. L.Classical correlation and quantum discord sharing of Dirac fields in noninertial frames20108152-s2.0-7795321545310.1103/PhysRevA.81.052120052120BrádlerK.WildeM. M.VinjanampathyS.UskovD. B.Identifying the quantum correlations in light-harvesting complexes20108262-s2.0-7865098634510.1103/PhysRevA.82.062310062310MazieroJ.CéleriL. C.SerraR. M.VedralV.Classical and quantum correlations under decoherence20098042-s2.0-7034986844810.1103/PhysRevA.80.044102044102MazzolaL.PiiloJ.ManiscalcoS.Sudden transition between classical and quantum decoherence2010104202-s2.0-7795253406410.1103/PhysRevLett.104.200401200401HeQ. L.XuJ. B.YaoD. X.ZhangY. Q.Sudden transition between classical and quantum decoherence in dissipative cavity QED and stationary quantum discord20118422-s2.0-8005163747710.1103/PhysRevA.84.022312022312MazieroJ.WerlangT.FanchiniF. F.CéleriL. C.SerraR. M.System-reservoir dynamics of quantum and classical correlations20108122-s2.0-7724909647810.1103/PhysRevA.81.022116022116YanX. Q.YueZ. L.Dynamics of quantum and classical correlations of a two-atom system in thermal reservoirs20135711712210.1016/j.chaos.2013.09.006TanaśR.Evolution of quantum correlations in a two-atom system2013T15301405910.1088/0031-8949/2013/T153/014059PianiM.GharibianS.AdessoG.CalsamigliaJ.HorodeckiP.WinterA.All nonclassical correlations can be activated into distillable entanglement2011106222204032204072-s2.0-7996062935910.1103/PhysRevLett.106.220403AliM.RauA. R. P.AlberG.Quantum discord for two-qubit X states20108142-s2.0-7865097917610.1103/PhysRevA.82.069902042105LuoS.Quantum discord for two-qubit systems20087742-s2.0-4174911024710.1103/PhysRevA.77.042303042303WoottersW. K.Entanglement of formation of an arbitrary state of two qubits19988010224522482-s2.0-4243684526IkramM.LiF.ZubairyM. S.Disentanglement in a two-qubit system subjected to dissipation environments20077562-s2.0-3454726803410.1103/PhysRevA.75.062336062336ScullyM. O.ZubairyM. S.1997England, CambridgeCambridge UniversityYuT.EberlyJ. H.Finite-time disentanglement via spontaneous emission2004931411404042-s2.0-1964440064910.1103/PhysRevLett.93.140404RoszakK.HorodeckiP.HorodeckiR.Sudden death of effective entanglement20108142-s2.0-7795112163410.1103/PhysRevA.81.042308042308YanX. Q.Entanglement sudden death of two atoms successive passing a cavity2009414164516502-s2.0-6734928304310.1016/j.chaos.2008.07.007BellomoB.Lo FrancoR.CompagnoG.Non-Markovian effects on the dynamics of entanglement200799162-s2.0-3534901588610.1103/PhysRevLett.99.160502160502FicekZ.TanaśR.Dark periods and revivals of entanglement in a two-qubit system20067422-s2.0-3374866990610.1103/PhysRevA.74.024304024304