Arxiv:1603.06676V2 [Quant-Ph]

Home , Quantum discord

arXiv:1603.06676v2 [quant-ph] 14 Apr 2016 et hrfr ti eyipratadncsayt cre- to necessary and important disentangle- very is as it well therefore as result ment, decoherence which systems’ environments quantum external in in- inevitably their systems with any teracting to due prone effect environmental and systems. to fragile are realistic correlations in quantum manipulated Furthermore, cre- be and to maintained, difficult very at- ated, usually much are however, attracted cor- they has tention, quantum models of and tasks computational study communication the quantum quantum decade, to past applying the relations Over quantum provides that 6]. tasks it proven [5, information and quantum been some entanglement for has advantages without quantum It exist can Ol- [4]. discord by Zurek introduced and been dis- liver correlations initially quantum quantum has hand, entanglement, general other than more the captures impor- On which very cord 3]. be [2, quantum to and quantumcomputations networks proven communication of been quantum type for has tant special it a and cru- is pro- correlations, Entanglement the information as quantum [1]. in proposed cessing applications been resources have cial discord, quantum and ∗ § ‡ † [email protected] [email protected] mff[email protected] [email protected] unu orltos uha unu entanglement quantum as such correlations, Quantum eeainadpoeto fsed-tt unu correla quantum steady-state of protection and Generation 3 e aoaoyo o-iesoa unu tutrsadQ and Structures Quantum Low-Dimensional of Laboratory Key n eateto hsc,HnnNra nvriy Changs University, Normal Hunan Physics, of Department and 4 unu orltoscnb rsre taln iei h li the in time long a asym at the preserved t in be show only We can displays correlations state. discord quantum entangled quantum an while o in that disappears, prepared with initially entanglement non of is exists dynamics system there the that compare shows we also Besides, result The state. entangled channel raini eae otemmr offiin fchannel of coefficient memory ca the memory to with related channel is damping creation amplitude to entanglemen steady-state quan subject that, of glement show dynamics results on an The memory Some with detail. amplit channels in account. of as effect into such The taken memory, sented. been with have channels channels successiv noisy depolarizing where of channels types quantum Different through correlations tum ASnmes 36.m 36.x 36.d 85.35.Be 03.65.Ud, 03.67.Hx, 73.63.Nm, numbers: PACS olg fSine ua nvriyo ehooy Zhuzhou Technology, of University Hunan Science, of College ehv rpsdashm ftegnrto n preservation and generation the of scheme a proposed have We .INTRODUCTION I. µ o-egGuo, You-neng s h oeteetnlmn raini,adteerirth earlier the and is, creation entanglement the more the is, hnsaUiest,Cagh,Hnn402,Pol’ Rep People’s 410022, Hunan Changsha, University, Changsha hnsaUiest,Cagh,Hnn402,Pol’ Rep People’s 410022, Hunan Changsha, University, Changsha 2 ua rvneKyLbrtr fApidEvrnetlPho Environmental Applied of Laboratory Key Province Hunan 1 eateto lcrncadCmuiainEngineering, Communication and Electronic of Department ,2, 1, ∗ a-aFang, Mao-fa hneswt memory with channels 3, † u-o Wang, Guo-you eosrtd[34] nadto,teehv ensev- been have there addition, been In has decoherence [43–48]. damping demonstrated amplitude the correla- from quantum tions protect to measurement measurement reversal combat- weak its and the in using role Recently, ef- significant a decoherence. memory ting play the environment and the studies, environment of previous fect the to these of de- back-action above used tackling the mentioned be by As also correlations [42]can coherence. quantum effect Zeno bipartite quantum protect [39–41] as detuning decoupling well on using dynamical as based Besides, or strategies [38]. 37], modulation 35], [36, subspaces [34, decoherence-free codes such decoherence, error-correcting from as correlations quantum protect to entan- 31]. induce [29, to subsystems able among is environment glement field common Recently, electromagnetic dissipative thermal [28]. a or even environment that bath happens heat immers- by it common created a be in also entangle-ing can the qubits two instance, between For ment other. indirect each their interac- with to whose lead interactions subsystems system, decoupled ancillary the former with simple The tions a introduce [27–33]. to systems create is systems quantum to coupled quantum methods direct coupled of or indirect use either in people entanglement Usually, [21–26] experimentally years. and recent [7–20] studied widely theoretically been both have correlations quantum tipartite of field the in correlations information. quantum quantum the preserve and ate ttesm ie hr r aymtosdeveloped methods many are there time, same the At mul- or bipartite of preservation and generation The µ ewe w nta uiswtotentan- without qubits initial two between t h togrtemmr offiin of coefficient memory the stronger The . oaiyi h bec fentanglement. of absence the in locality sso h hnesaecorrelated. are channels the of uses e atmCnrlo iityo Education, of Ministry of Control uantum a nageetdnmc suddenly dynamics entanglement hat ttclmt utemr,two-qubit Furthermore, limit. ptotic a408,Pol’ eulco China of Republic People’s 410081, ha ltclo ueia eut r pre- are results numerical or alytical i of mit unu icr hnatwo-qubit a when discord quantum f 108 epesRpbi fChina of Republic People’s 412008, u orltoshsbe discussed been has correlations tum d apn,paedmig and phase-damping, damping, ude egnrtd h entanglement The generated. be n 4, ftoqbtsed-tt quan- steady-state two-qubit of µ ‡ eaal tt eoe the becomes state separable e → n eZeng Ke and bi fChina of ublic bi fChina of ublic 1. in u oquantum to due tions tocatalysis, ,2, 1, § 2 eral interesting works [49, 50] that focused on the pro- lations, this can happen when environmental correlation tection of the long time limit of quantum correlations time is longer than the time between consecutive uses, so through the addition of qubits. that the channel acts dependently on each channel input, ⊗N In this paper, we aim to study how the correlated noise εN = ε . These kinds of channels are called memory 6 acting on consecutive uses of channel influences the quan- channels. tum correlations dynamics. In real physical quantum In what follows we consider N channel uses. Given an transmission channels, the noise over consecutive uses input state ρ, a quantum channel ε is defined as a com- provide a natural theoretical framework for the study pletely positive, trace-preserving map from input-state of any noisy quantum communication since they were density matrices to output-state density matrices, first considered by C. Macchiavello and G.M. Palma [51] † to study the problem of classical capacity. Later, the ε(ρ)= EiρEi (1) study of information transmission over a correlated chan- Xi nel with an arbitrary degree of memory in quantum information process has attracted much attention [52–57]. where Ei = Pi1...iN Ai are the Kraus operators of the Recently, a report about the effect of correlated noise on channel whichp satisfy the completeness relationship, and the entanglement of X-type state of the Dirac fields in the i Pi1...iN = 1. Here Pi1...iN can be interpreted as the non-inertial frame has been investigated by M. Ramzan probabilityP that a random sequence of operations is ap- [54]. We here, have analyzed the generation and preser- plied to the sequence of N qubits transmitted through vation of two-qubit quantum correlations through quan- the channel. For a memoryless channel, these operations tum channels where successive uses of the channels are Ai are independent, and we have Pi1 ...iN = Pi1 Pi2 ...PiN . correlated. Different types of noisy channels with mem- However, for a channel with memory, these operations ory, such as amplitude damping, phase-damping, and de- are time-correlated, Pi1...iN = Pi1 Pi2|i1 ...PiN |iN−1 , here polarizing channels have been taken into account. The PiN |iN−1 is the conditional probability for that operation. effect of channels with memory on dynamics of quantum For simplicity, we consider the Kraus operators for two correlations has been discussed in detail. The results consecutive uses of a channel with partial memory are show that, steady-state entanglement between two independent qubits initially without entanglement subject to Ei,j = Pi[(1 µ)Pj + µδi,j ]Ai Aj (2) amplitude damping channel with memory can be gener- q − ⊗ ated. We observe that there exists nonlocality quantified where 0 µ 1 is the memory coefficient of channel. by quantum discord in the absence of entanglement. Be- Following,≤ ≤ we will focus on the noisy channels (e.g. am- sides, we compare the dynamics of entanglement with plitude damping, phase-damping, and depolarizing chan- that of quantum discord when a two-qubit system is ini- nels) with time-correlated Markov noise for two consec- tially prepared in an entangled state. We show that en- utive uses. Based on the Kraus operator approach, for tanglement dynamics suddenly disappears, whlie quan- any initial state ρ, the finial state under noise is given by tum discord vanishes only in the asymptotic limit. Two- [51, 55] qubit quantum correlations can be preserved at a long time due to the channels with prefect memory. † † ε(ρ)=(1 µ) Ei,j ρE + µ Ek,kρE (3) The layout is as follows: In Sec. II, we illustrate the − i,j k,k Xi,j Xk initial states and noise channels. In Sec. III, we devote to examining generation and protection of two-qubit quan- One can understand the above expression that the same tum correlations in different types of noisy channels with operation is applied to both qubits with probability µ memory. Finally, we give the conclusion in Sec. IV. while with probability 1 µ both operations are uncorrelated. −

II. INITIAL STATES AND NOISE MODEL A. Amplitude damping channel with memory We begin with a brief description of quantum-memory channels. As we all know, there are two diﬀerent types Amplitude damping channel which is used to charac- of quantum channels including memoryless channels and terize spontaneous emission describes the energy dissipa- memory channels. The simplest models for quantum tion from a quantum system. The Kraus operators for a channels are memoryless, when environmental correla- single qubit are given by tion time is smaller than the time between consecutive uses, so that at each channel use the environment back √1 p 0 A0 = − (4) action can be negligible. Namely, the system under- 0 1 takes the same quantum channel ε, in which independent noise acts on each use. Suppose N times uses of this channel, then we have ε = ε⊗N . However, real sys- 0 0 N A = (5) tems among subsequent channel uses exhibit some corre- 1 √p 0 3

0.30 0.10 Μ = 1.0 ( a ) 0.25 0.08 Μ = 1.0 0.20 Μ = 0.7 0.06 Μ = 0.7 0.15 QE QD 0.04 Μ = 0.4 0.10 Μ = 0.4 0.05 0.02 Μ = 0.1 Μ = 0.1 Μ = 0.0 0.00 Μ = 0.0 0.00 0 5 10 15 20 0 1 2 3 4 5 Γt Γt 0.10 FIG. 1. (Color online) Quantum entanglement of a two-qubit ( b ) system initially prepared in an unentangled state under am- 0.08 plitude damping channels with memory is plotted against γt Μ = 1.0 for different µ. Other parameters: α = 0.5,r = 0.3. 0.06 Μ = 0.7 QD 0.04 Μ = 0.4 0.02 Μ = 0.1 Μ = 0.0 0.00 0 1 2 3 4 5 0.3 Γt QE 0.2 1.0 0.1 0.10 ( c ) 0.0 0.08 Μ = 1.0 0.0 0.5Μ 0.1 0.06 r 0.2 Μ = 0.7 0.3 0.0 QD 0.04 Μ = 0.4 FIG. 2. (Color online) Generation of quantum entanglement 0.02 Μ = 0.1 is plotted against r and µ under amplitude damping channels Μ = 0.0 with memory. Other parameters: α = 0.5,p = 0.95. 0.00 0 1 2 3 4 5 Γt where p 1 exp( γt) is the single qubit damping rate. ≡ − − FIG. 3. (Color online) Quantum discords of a two-qubit sys- Suppose that two qubits pass through the quantum de- tem initially prepared in an unentangled state under differ- phasing channel, the environmental correlation time is ent quantum channels with memory are plotted against γt smaller than the time between consecutive uses, and then for different µ. (a) Amplitude damping channels, (b) phase the quantum amplitude damping channel with uncorre- damping channels, (c) depolarizing damping channels. Other lated noise can similarly be defined as the following Kraus parameters: α = 0.5,r = 0.3. operators

Ei,j = Ai Aj , (i, j =0, 1) (6) 0 000 ⊗ 0 000 E =   (8) Based on previous analysis [51, 55], an amplitude- 11 0 000   damping channel with ﬁnite memory for two consecutive  √p 0 0 0  uses, the task of constructing Kraus operators Ek,k for the amplitude-damping channel is given by Consider the following initial states 1 r ρ(0) = r Ψ Ψ + − I (9) √1 p 0 0 0 | ih | 4 0− 100 E =   (7) 00 0 010 where Ψ = √1 α2 01 + α 10 corresponds to the Bell-  0 001  like states with− 0 |α i 1, and| i r indicates the purity of   ≤ ≤ 4

0.5 0.20 ( a ) ( a ) Μ = 1.0 0.4 Μ = 1.0 0.15 Μ = 0.7 0.3 Μ = 0.7 QD 0.10 Μ = 0.4 QE 0.2 Μ = 0.1 Μ = 0.4 0.05 0.1 Μ = 0.0 Μ = 0.1 0.00 0.0 Μ = 0.0 0 1 2 3 4 5 0 1 2 3 4 5 Γt Γt 0.20 ( b ) 0.20 Μ = 1.0 ( b ) Μ = 1.0 0.15 0.15 Μ = 0.7

QD 0.10 0.10 Μ = 0.7 Μ = 0.4 QE 0.05 Μ = 0.1 0.05 Μ = 0.4 Μ = 0.0 Μ = 0.0 Μ = 0.1 0.00 0 1 2 3 4 5 0.00 0 1 2 3 4 5 Γt Γt 0.20 ( c ) 0.20 Μ = 1.0 ( c ) Μ = 1.0 0.15 0.15 Μ = 0.7

Μ = 0.1 QD 0.10 Μ = 0.4 0.10 QE 0.05 Μ = 0.0 Μ = 0.1 0.05 Μ = 0.0 0.00 0 1 2 3 4 5 Μ = 0.4 Μ = 0.7 0.00 0 1 2 3 4 5 Γt

Γt FIG. 5. (Color online)Quantum discords of a two-qubit system initially prepared in an entangled state under different FIG. 4. (Color online)Quantum entanglements of a two-qubit quantum channels with memory are plotted against γt for dif- system initially prepared in an entangled state under differ- ferent µ. (a) Amplitude damping channels, (b) phase damp- ent quantum channels with memory are plotted against γt ing channels, (c) depolarizing damping channels. Other pa- for different µ. (a) Amplitude damping channels, (b) phase rameters: α = 0.5,r = 0.5. damping channels, (c) depolarizing damping channels. Other parameters: α = 0.5,r = 0.5. whose the density matrix elements are 1 ρ = (1 p)(1 r)[1 p(1 µ)], the initial states. According to Eq. (3), through calcula- 11 4 − − − − tions, it is not difficult to get the density matrix of a two qubits system under an amplitude-damping channel with 1 ρ = [1+r(3 4α2) 4pr(1 α2)(1 µ) p2(1 r)(1 µ)], time-correlated Markov noise for two consecutive uses, 22 4 − − − − − − −

1 2 2 2 ρ33 = [1 r(1 4α ) p (1 r)(1 µ) 4prα (1 µ)], ρ11(t)0 0 0 4 − − − − − − − 0 ρ22(t) ρ23(t) 0 ρAB(t)=   (10) 0 ρ32(t) ρ33(t) 0 1 2   ρ44 = [1 r + p (1 r)(1 µ)+2p pµ + pr(2 3µ)],  0 0 0 ρ44(t)  4 − − − − − 5

∗ 2 ρ23 = ρ32 = rα 1 α [1 p(1 µ)]. (11) Here we will consider the case of two consecutive uses p − − − of a channel with partial memory, the task of constructing Kraus operators Ek,k for the non-Pauli depolarizing B. Phase-damping channel with memory channel is given [57]

E = P σ σ , (k =0, 1, 2, 3) (17) Phase-damping channel describes a quantum noise k,k k k k p ⊗ with loss of quantum phase information without loss of According to Eq. (3), the density matrix elements of energy. The Kraus operators for a single qubit are de- a two qubits system under a depolarizing channel with ﬁned in terms of the Pauli operators σ0 = I and σ3. time-correlated Markov noise for two consecutive uses are Suppose that two qubits pass through the quantum 1 2 2 2 dephasing channel, the channel with uncorrelated noise ρ11 = ρ44 = [9+r(36α 9+16p (1 µ) 24p(2α µ))], can similarly be deﬁned by the following Kraus operators 36 − − − −

Ei,j = PiPj σi σj (12) 1 ⊗ ρ = [9+r(27 36α2 +16p2(1 µ)+24p(2α2 +µ 2))], p 22 36 − − − and two consecutive uses of the dephasing channel with partial memory, the Kraus operators E is given as [56] k,k 1 ρ = [9 + r(36α2 9+16p2(1 µ)+24p(µ 2α2))], 33 36 − − − Ek,k = Pkσk σk (13) p ⊗ 1 where i, j, k = (0, 3), and P0 = 1 p, P3 = p, where ρ = ρ∗ = rα 1 α2[(3 4p)2(1 µ)+9µ]. (18) 1 − 23 32 9 − − − p 2 [1 exp( γt)]. Similarly according to Eq. (3), the p density≡ − matrix− elements of a two qubits system under a phase-damping channel with time-correlated Markov III. GENERATION AND PROTECTION OF noise for two consecutive uses are TWO-QUBIT QUANTUM CORRELATIONS 1 UNDER QUANTUM CHANNELS WITH ρ = ρ = (1 r), MEMORY 11 44 4 − Before investigating the generation and protection of 1 two-qubit quantum correlations under quantum chan- ρ = [1 + r(3 4α2)], 22 4 − nels with memory, let us review the correlation mea- sures quantified by quantum discord (QD) [4] which is defined by subtracting the classical correlation C(ρAB ) 1 2 ρ33 = [1 r(1 4α )], from the total amount of correlation I(ρA:B), namely, 4 − − QD(ρAB) = I(ρA:B) C(ρAB ). The total correlation is quantified by the− quantum mutual information I(ρA:B) = S(ρA) + S(ρB) + S(ρAB), where ρA(ρB ) is ρ = ρ∗ = rα 1 α2[(1 p)2(1 µ)+ µ]. (14) 23 32 − − − the reduced matrix of ρAB by tracing out B(A). The p classical correlation C(ρAB) is defined as C(ρAB) = S(ρA) min B S(ρA B), where S(ρ) = tr(ρ log ρ) − {Qk } | − 2 C. Depolarizing channel with memory is the von Neumann entropy. Note that the minimum is taken over the set of positive operator valued measure- Depolarizing channel is another important type of ment B on subsystem B, S(ρ ) is the conditional { k } A|B quantum noise, it describes the process in which the den- entropyQ for the subsystem A. sity matrix is dynamically replaced by the state I/2. I It is worthy pointing out that the calculation of classi- denoting identity matrix of a qubit. The Kraus operators cal correlation involves a potentially complex optimiza- for a single qubit are given by tion process even numerically, usually there is no general analytical expression of discord except for the simplest Ai = Piσi, (i =0, 1, 2, 3) (15) case of two-qubit state. Such as for a bipartite quantum p X-state described by the density matrix ρAB, an expres- 1 sion of the quantum discord is given [58] where P0 = 1 p, P1 = P2 = P3 = p/3, and p 2 [1 exp( γt)]. − ≡ − − Assume that the environmental correlation time is QD(ρAB) = min(Q1,Q2) (19) smaller than the time between consecutive uses, a chan- where nel is memoryless and corresponds to its Kraus operators 4 Qi = H(ρ11 + ρ33)+Σi=1ǫi log2 ǫi + Dj, 2 2 1+√[1−2(ρ33+ρ44)] +4(|ρ14|+|ρ23|) Ei,j = PiPj σi σj , (i, j =0, 1, 2, 3) (16) D = H( ), ⊗ 1 2 p 6

D = Σiρii log ρii H(ρ + ρ ), quantum channels with memory. The results show that 2 − 2 − 11 33 H(x)= x log2 x (1 x) log2(1 x). in a case where there is initially no entanglement between Besides,− to investigate− − the two-qubit− quantum correla- two-qubit system but there exists quantum discord. This tions dynamics, we also use concurrence as the quantifier means quantum discord captures more general quantum [59]. For the density matrix of a bipartite system has a correlations than entanglements. It is proven that the form given by Eq. (10), the concurrence reads states with nonzero quantum discord but not entanglement are responsible for the efficiency of a quantum com- QE(ρ ) = 2 max(0, ρ ρ ρ , ρ ρ ρ ) AB 23 √ 11 44 14 √ 22 33 puter [5, 6]. Besides, quantum correlations quantified by | |− | |− (20) quantum discord exhibit in the asymptotic limit, and can Evidence suggests entanglement quantified by concur- be well protected in the presence of channels with mem- rence varies from QE = 0 for a separable state to QE =1 ory. Especially, as the memory coefficient of channel µ for a maximally entangled state. increases, quantum discord can be protected more effec- In the following, we adopt both quantum discord and tively. In the limit of µ 1, the long-living quantum concurrence as a measure of quantum correlations for a discord preservation can be→ observed. two uncoupled qubits system initially prepared in Eq.(9) Finally, we investigate the dynamics of quantum corre- successive uses of noisy channels with memory, such as lations for two initial entangled qubits (e.g.α = 0.5, r = amplitude damping, Phase-damping, and depolarizing 0.5) in different types of quantum channels with memory, channels. One can easily check the initial state given and compare the dynamics of entanglement with that by Eq.(9), the concurrence QE(ρ ) = 0 for the purity 0 of quantum discord. The results show that, the entan- 0 r 1/3, while the purity 1/3 < r 1 corresponds glement dynamics disappears in a finite time under the to≤ the≤ concurrence QE(ρ ) > 0. ≤ 0 influence of noise as shown Fig. 4. This behaviour is Firstly, we consider the case where the initial state named entanglement sudden death, and the results are given by Eq.(9) is prepared in an unentangled state consistent with Refs. [54]. While quantum discord ex- (e.g.α =0.5, r =0.3), when two qubits pass through the hibits only in the asymptotic limit in Fig. 5. This in- quantum channels with memory. Due to the amplitude dicates quantum discord is more robustness against the damping channels with memory, there is no entanglement noise than entanglement. Besides, the amount of ini- at earlier times QE(ρ ) = 0, and at some time the quan- 0 tial entanglement QE(0) 0.18 possesses less general tum entanglement of two-qubit system starts to build up quantum correlations than≈ the initial quantum discord (see Fig. 1). This phenomenon is called delayed sudden QD(0) 0.20. Take quantum channels with memory birth of entanglement [31]. Note that the degree of en- into consideration,≈ as we can see that, with the memory tanglement depends on the memory coefficient of channel coefficient of channel µ increasing, it can completely cir- µ. The stronger the memory coefficient of channel µ is, cumvent the entanglement sudden death. Particularly, in the more the entanglement creation is, and the earlier the limit of µ 1, the memory effect can help preserve the separable state becomes the entangled state. This entanglement and→ quantum discord at a long time. More- implies the amplitude damping channels with memory over, we find, when two qubits pass through amplitude can locally create entanglement, and it can be preserved damping channels with memory, the amount of entan- an unentangled state at a steady entangled state in the glement creation is more than the initial one, as shown presence of the quantum channels with memory. Sur- Fig.4(a). prisingly, even though Phase-damping, and depolarizing channel with memory, two-qubit system can not be induced entanglement. To get a better understanding of the effects of the am- IV. CONCLUSION plitude damping channel with memory on the entanglement creation, we plot Fig.2 to show the entanglement In conclusion, we have proposed a scheme of the gener- creation as a function of initial state purity r and mem- ation and preservation of two-qubit steady state quantum ory coefficient of channel µ for p =0.95. One can see that correlations through quantum channels where successive the entanglement creation is not only depended on the uses of the channels are correlated. Different types of memory coefficient of channel µ but also on the initial noisy channels with memory, such as amplitude damp- state chosen. However, when chosen p = 1, we find that, ing, phase-damping, and depolarizing channels have been 2 for any α and 0 r 1/3, QE(ρt)= rα√1 α µ 0, taken into account. The effect of channels with mem- which means that≤ two≤ uncoupled qubits system− initially≥ ory on dynamics of quantum correlations has been dis- prepared in any separable states given by Eq.(9) can be cussed in detail. We conclude the results as follows: induced entanglement with the help of the amplitude Firstly, steady state entanglement between two indepen- damping channels with memory. Our result indicates dent qubits without entanglement subject to amplitude that the amplitude damping channels with memory can damping channel with memory can be generated, but be locally created entanglement, which is attributed to µ Phase-damping, and depolarizing channel with memory as well as p. can not. Secondly, the amplitude damping channels with Fig.3 shows quantum discord dynamics for two initial memory can locally create entanglement, the stronger the unentangled qubits (e.g.α = 0.5, r = 0.3) in different memory coefficient of channel µ is, the more the entan- 7 glement creation is, and the earlier the separable state efficient of channel µ increases, the phenomena of entan- becomes the entangled state. Thirdly, there is initially glement sudden death can be completely circumvented. no entanglement between two-qubit system but there ex- Particularly, two-qubit quantum correlations can be pre- ists quantum discord in different quantum channels with served at a long time in the limit of µ 1 memory. Finally, we compare the dynamics of entan- → glement with that of quantum discord when a two-qubit system is prepared in an entangled state. We show that ACKNOWLEDGMENTS entanglement dynamics suddenly disappears, namely entanglement sudden death occurs, while quantum discord This work is supported by the National Natural Sci- displays only in the asymptotic limit. As the memory co- ence Foundation of China (Grant Nos. 11374096).

[1] M.A. Nielsen, I.L. Chuang, Quantum Computation and [25] J.C. Lee, Y.C. Jeong, Y.S. Kim, Y.H. Kim, Opt. Express Quantum Information (2000) 19 (2011) 16309. [2] R. Horodecki, P. Horodecki, M. Horodecki, K. Horodecki, [26] B. Bellomo, M. Antezza, Europhys. Lett. 104 (2013) Rev. Mod. Phys. 81 (2009) 865. 10006. [3] R. Lo Franco, G. Compagno, Sci. Rep. 6 (2016) 20603. [27] D. Braun, Phys. Rev. Lett. 89 (2002) 277901. [4] H. Ollivier, W.H. Zurek, Phys. Rev. Lett. 88 (2001) [28] F. Benatti, R. Floreanini, M. Piani, Phys. Rev. Lett. 91 017901. (2003) 070402. [5] A. Datta, S.T. Flammia, C.M. Caves, Phys. Rev. A 72 [29] Y. Li, J. Zhou, H. Guo, Phys. Rev. A 79 (2009) 012309. (2005) 042316. [30] N. Ba An, J. Kim, K. Kim, Phys. Rev. A 84 (2011) [6] A. Datta, A. Shaji, C.M. Caves, Phys. Rev. Lett. 100 022329. (2008) 050502. [31] Z. Ficek, R. Tana´s, Phys. Rev. A 77 (2008) 054301. [7] J. Wang, H.M. Wiseman, G.J. Milburn, Phys. Rev. A 71 [32] A. González-Tudela et al, Phys. Rev. Lett. 106 (2011) (2005) 042309. 020501. [8] F.F. Fanchini, R.d.J. Napolitano, Phys. Rev. A 76 (2007) [33] K. Almutairi, R. Tana´s, Z. Ficek, Phys. Rev. A 84 (2011) 062306. 013831. [9] J. Zhang, R.B. Wu, C.W. Li, T.J. Tarn, J. Phys. A: [34] P.W. Shor, Phys. Rev. A 52 (1995) 2493. Math. Theor. 42 (2009) 035304. [35] A.M. Steane, Phys. Rev. Lett. 77 (1996) 793. [10] T.J. G. Apollaro, A. Cuccoli, C. Di Franco, M. Pater- [36] D.A. Lidar, I.L. Chuang, K.B. Whaley, Phys. Rev. Lett. nostro, F. Plastina, P. Verrucchi, New J. Phys. 12 (2010) 81 (1998) 2594. 3615-3620. [37] P.G. Kwiat, A.J. Berglund, J.B. Altepeter, A.G. White, [11] S.Y. Lin, B. L. Hu, Phys. Rev. D 81 (2010) 045019. Science 290 (2000) 498. [12] H.S. Xu, J.b. Xu, J. Opt. Soci. America B 29 (2012) [38] X. Xiao, Y. Li, K. Zeng, C. Wu, J.Phys. B: At. Mol. Opt. 2074-2079. Phys. 42 (2009) 235502. [13] M. Siomau, Ali A. Kamli, S. A. Moiseev, B.C. Sanders, [39] L. Viola, E. Knill, S. Lloyd, Phys. Rev. Lett. 82 (1999) Phys. Rev. A 85 (2012) 050303. 2417. [14] C. Creatore, R.T. Brierley, R.T. Phillips, P.B. Little- [40] J.R. West, D.A. Lidar, B.H. Fong, M.F. Gyure, Phys. wood, P.R. Eastham, Phys. Rev. B 86 (2012) 5505-5511. Rev. Lett. 105 (2010) 230503. [15] B. Bellomo, M. Antezza, New J.Phys. 15 (2013) 113052. [41] Z.X. Man, Y.J. Xia, and R. Lo Franco, Phys. Rev. A 92 [16] J. Cheng, W.Z. Zhang, L.Zhou, W.p. Zhang, arXiv: (2015) 012315. Quant-Ph (2015)1511.01267 [42] S. Maniscalco, F. Francisca, R. Zaffino, N. Gullo, F. Plas- [17] J.L. Guo, J.L. Wei, Annals of Physics, 354 (2015) 522- tina, Phys. Rev. Lett. 100 (2008) 090503. 533. [43] A.N. Korotkov, K. Keane, Phys. Rev. A 81 (2010) [18] A. D’Arrigo, R. Lo Franco, G. Benenti, E. Paladino, G. 040103. Falci, Ann. Phys. 350 (2014) 211. [44] Q. Sun, M. Al-Amri, M.S. Zubairy, Phys. Rev. A 80 [19] M. Cianciaruso, T. R. Bromley, W. Roga, R. Lo Franco, (2009) 033838. G. Adesso, Sci. Rep. 5 (2015) 10177. [45] Y.S. Kim, J.C. Lee, O. Kwon, Y.H. Kim, Nat. Phys. 8 [20] Z.X. Man, Y.J. Xia, and R. Lo Franco, Sci. Rep. 5 (2015) (2012) 117. 13843. [46] X. Xiao, Y.L Li, Eur. Phys. J. D 67 (2013) 204. [21] J.S. Xu, K. Sun, C.F. Li, X.Y. Xu, G.C. Guo, E. An- [47] S.C. Wang, Z.W. Yu, W.J. Zou, X.B. Wang, Phys. Rev. dersson, R. Lo Franco G. Compagno, Nat. Commun. 4 A 89 (2014) 022318. (2013) 2851. [48] J. Yune, K.H. Hong, H.T. Lim, J.C. Lee, O. Kwon, [22] A. Orieux, G. Ferranti, A. D’Arrigo, R. Lo Franco, G. S.W. Han, Y.S.Kim, S. Moon, Y.H. Kim Opt. Express Benenti, E. Paladino, G. Falci, F. Sciarrino, P. Mataloni, 23 (2015) 26012-26022. Sci. Rep. 5 (2015) 8575. [49] M. Merkli, G. P. Breman, F. Borgonovi, V.I. Tsifrinovich, [23] X.Q. Shao, T.Y. Zheng, C.H. Oh, S. Zhang, Phys. Rev. Advance in Math. Phys. (2012) 375182. A 89 (2014) 012319. [50] M.M. Flores, E.A. Galapon, Annals of Physics, 354 [24] R. Blattmann, H.J. Krenner, S. Kohler, P. Hänggi, Phys. (2015) 21-30. Rev. A 89 (2014) 012327. 8

[51] C. Macchiavello, G. M. Palma, Phys. Rev. A 65 (2002) [56] A. DArrigo, G. Benenti, G. Falci, New J. Phys. 9 (2007) 050301. 310. [52] A. DArrigo, G. Benenti, G. Falci, C. Macchiavello, Phys. [57] C. Macchiavello, G. M. Palma, S. Virmani, Phys. Rev. A Rev. A 88 (2013) 042337 69 (2004) 010303. [53] A. D’Arrigo, G. Benenti, G. Falci, C. Macchiavello, Phys. [58] C.Z. Wang, C.X. Li, L.Y. Nie, J.F. Li, J.Phys. B: At. Rev. A 92 (2015) 062342 Mol. Opt. Phys. 44 (2011) 015503. [54] M. Ramzan, Quantum Inf. Process. 12 (2013) 83-95. [59] W. K. Wootters, Phys. Rev. Lett. 80 (1998) 2245. [55] Y. Yeo, A. Skeen, Phys. Rev. A 67 (2003) 064301.