Entanglement and non-Markovianity of a multi-level atom decaying in a cavity
Fan Zi-Long , Ren Yu-Kun , Zeng Hao-Sheng †,
Key Laboratory of Low-Dimensional Quantum Structures and Quantum Control of Ministry of Education,Department of Physics, Hunan Normal University, Changsha 410081, China

 

† Corresponding author. E-mail: hszeng@hunnu.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11275064 and 11075050), the Specialized Research Fund for the Doctoral Program of Higher Education, China (Grant No. 20124306110003), and the Construct Program of the National Key Discipline, China.

Abstract
Abstract

We present a paradigmatic method for exactly studying non-Markovian dynamics of a multi-level V-type atom interacting with a zero-temperature bosonic bath. Special attention is paid to the entanglement evolution and the dynamical non-Markovianity of a three-level V-type atom. We find that the entanglement negativity decays faster and non-Markovianity is smaller in the resonance regions than those in the non-resonance regions. More importantly, the quantum interference between the dynamical non-Markovianities induced by different transition channels is manifested, and the frequency domains for constructive and destructive interferences are found.

1. Introduction

Quantum non-Markovian dynamics, due to its wide existence in the quantum optical system, [ 1 ] quantum dot, [ 2 ] superconductor system, [ 3 ] quantum chemistry, [ 4 ] and biological system [ 5 ] and some possible applications in quantum metrology [ 6 ] and quantum communication, [ 7 10 ] has received much attention in recent years. Several proposals [ 11 17 ] for the measure of non-Markovianity have been presented and various dynamical properties [ 18 39 ] of non-Markovian processes have been investigated. Experimentally, the simulation of non-Markovianity [ 40 , 41 ] under controlled environments has been realized.

The nature of a quantum non-Markovian process is the flow of lost information from the environment back to the open system, which leads to the oscillations of some pivotal physical quantities such as quantum trace distance, [ 11 ] quantum correlation, [ 12 , 14 ] and quantum Fisher information. [ 13 ] Although the properties of these oscillations are not in complete agreement, each of them can serve as the signature of non-Markovian dynamics. We will study the dynamical non-Markovianity in terms of the entanglement negativity. [ 42 , 43 ]

The study of non-Markovian dynamics of open quantum systems is typically very involved. Most of the previous works about non-Markovian dynamics mainly concentrated on the two-level systems. The high-dimensional open quantum systems, especially the high-dimensional dissipative quantum systems, due to their complexity, have been seldom investigated. In this paper, we present a paradigmatic method for exactly studying the non-Markovian dynamics of a multi-level V-type atom decaying in a zero-temperature bosonic bath. As a special example, we will emphatically study the case of a three-level V-type atom. The evolution of entanglement, the non-Markovianity of the quantum process, and the interference effect between different transition channels will be the main contents investigated.

The paper is organized as follows. In Section 2, we introduce the microscopic model for a multi-level V-type atom interacting with a zero-temperature bosonic bath and present the schematic method for treating the problem. In Section 3 and Section 4, we study respectively the evolution of entanglement (measured by negativity) and the non-Markovianity of the dynamics for the special case of a three-level V-type atom. Finally, the conclusion is given in Section 5.

2. Microscopic model

Consider an atom A embedded in a zero-temperature bosonic bath modeled by an infinite chain of quantum harmonic oscillators (see Fig.  1(a) ). The atom has V-type ( n + 1) levels, labeled as |0〉, |1〉, …, | n 〉, with transition frequencies ω i , i = 1,2,…, n (Fig.  1(b) ). The Hamiltonian of this system in the Schrödinger picture is given by

where ω k , a k , and are respectively the frequency, annihilation, and creation operators for the k -th harmonic oscillator of the bath. For simplicity but without loss of generality, we assume the coupling strengths g ik to be real.

(a) Schematic diagram of the model. The atom A (main system) is embedded in a bosonic bath and entangled initially with another identical ancillary atom B. The total initial state is | Ψ (0)〉. (b) Energy levels of the multi-level V-type atom.

For the following study, we now assume another identical ancillary atom B which is free from the environment but initially entangled with atom A. The initial state of the two atoms plus the environment is assumed to be

where | 0 〉 = |00···0〉 denotes the vacuum state of the bath. Due to the feature of Hamiltonian (1), the evolved state at any time t may be written as

where | 1 k 〉 = |00···1 k ···0〉 denotes the state that a photon is in mode k of the bath. As the state of atom B is unchangeable, we only need to consider the conservation of the total excitation number of atom A plus its bath. The first term in the above equation has total excitation number of zero; the second term has total excitation number of one, in which the state of atom A has altered from level j A to i A through the interaction with the bath; the last term also has total excitation number of one but the excitation in atom A has been transferred to the bath.

The state | Ψ ( t )〉 must satisfy the Schrodinger equation governed by Hamiltonian (1), which leads to the following equations about coefficients α i j ( t ) and β jk ( t ):

Integrating Eq. ( 5 ) formally and plugging it into Eq. ( 4 ), we have

where the correlation function is defined as

In the above derivation, we have used the continuum limit ∑ k g ik g i′k → ∫ d ωj ii′ ( ω ), with j ii′ ( ω ) being the transition-dependent spectral density. In this paper, we assume the spectral density to be a Lorentzian distribution

with central frequency ω 0 and width λ . The parameter Ω i describes the transition-dependent strength. For this spectral distribution, the correlation function may be further written as

Plugging this correlation function into Eq. ( 6 ) and performing the Laplace transformation, we have

where i , j = 1, …, n . This is a set of linear equations with respect to α i j ( p ), which can be solved through programmatic procedures. Furthermore, we find that the inverse Laplace transformation on α i j ( p ) is in principle feasible, because they can be written in the form of rational functions (see Eq. ( 12 ) in the next section). Combining with Eq. ( 4 ), we can thus obtain all the coefficients α i j ( t ) and β jk ( t ) needed for the study of system dynamics.

3. Entanglement negativity

In this section, we take the three-level V-type atom as a paraphrastic example to demonstrate the procedure for exactly solving the dynamics. We study the entanglement evolution in terms of negativity. The time evolution of entanglement in a one-sided noisy environment is an important model, which may be used to simulate the procedure of a real entanglement distribution. Suppose Alice wants to establish an entanglement resource with distant Bob. She first produces a local entangled state of two atoms A and B, and then sends atom B to Bob through a quantum channel and keeps atom A in her hand. If we suppose only the sent atom suffered from noise, then the procedure corresponds to the time evolution of entanglement in a one-sided environment. The non-local entanglement of high-dimensional quantum systems is very relevant because of the security in quantum cryptography. [ 44 , 45 ]

Tracing over the bath degrees of freedom, from Eq. ( 3 ), we obtain the reduced density matrix for atoms A and B

where the basis has been set as {|0 A j B 〉; |1 A j B 〉; |2 A j B 〉; | j = 0,1,2}. All the matrix elements may be determined through Eqs. ( 5 ) and ( 10 ). Firstly, equation ( 10 ) with the initial conditions α 11 (0) = α 22 (0) = 1 and α 12 (0) = α 21 (0) = 0 leads to

where

and ν 1 , ν 2 , ν 3 are the roots of the following equation with respect to ν :

The inverse Laplace transforms of Eq. ( 12 ) produce the following compact expressions for α i j ( t ):

with

for α 11 ( t );

for α 22 ( t ); and

for α 12 ( t ) and α 21 ( t ).

Secondly, let and , then equation ( 5 ) becomes

Integrating this equation with the initial condition β jk (0) = 0, we have

Multiplying this equation by its complex conjugation and summing over k , we obtain

Performing the continuum limit ∑ k g ik g i′k → ∫ d ω J ii′ ( ω ) and employing Eqs. ( 7 ) and ( 8 ), we finally arrive at

This expression applies in principle to the system of arbitrary high-dimensional V-type atoms, and for the case of three-level V-type atoms, we should let i , i′ , j , j′ = 1,2. Now, all the elements in density matrix (11) can be evaluated directly through Eqs. ( 16 ) and ( 20 ).

From Eq. ( 11 ), we can then study the evolution of entanglement for atom A with its ancillary atom B. We will employ the entanglement negativity [ 42 , 43 ] defined by

where and are respectively the negative and all eigenvalues of the partial transpose of ρ AB with respect to subsystem A. Through a tedious but straightforward calculation, we find and with , and the remaining four roots are defined by the equation

For simplifying the expressions of , we have used the relation α 12 = α 21 .

In Fig.  2 , we show the time evolution of entanglement negativity for different spectral widths. For all the numerical simulations involved in this paper, we always set , which has the frequency dimension. The dimensionless time is defined by τ = Ωt . It is shown for both non-degenerate ( ω 1 ω 2 , Fig.  2(a) ) and degenerate ( ω 1 = ω 2 , Fig.  2(b) ) transitions that the less the λ is, the severer the curve oscillates and the slower the negativity decreases. The former suggests that narrower spectral width can increase the non-Markovianity of open quantum process, and the latter suggests that stronger non-Markovianity can effectively prevent the decay of entanglement. These properties are similar to those observed in open two-level quantum systems. [ 26 ]

Time evolution of negativity for (a) non-degenerate case with ω 1 = 5 Ω , ω 2 = 7 Ω , and ω 0 = 6 Ω , and (b) degenerate case with ω 1 = ω 2 = ω 0 = 6 Ω .

Figure  3 shows the time evolution of entanglement negativity in the case of degenerate transition and for different detuning between the transition and the bath frequencies, where we choose ω 1 = ω 2 = 6 Ω , λ = 0.4 Ω , and the dimensionless detuning is defined by Δ = ( ω 0 ω 1 )/ Ω . We see that in the resonant regime ( ω 0 ω 1 ), the entanglement negativity decreases very quickly (Fig.  3(a) ), because the resonance is beneficial for the dissipation of the open quantum system. In order to see clearly the oscillation, we present a plane graph in Fig.  3(b) . It is shown that larger detuning between the transition and the bath frequencies can not only slow the decay of the entanglement negativity, but also improve the non-Markovianity of the quantum process. Similar results also have been obtained in the open two-level quantum systems. [ 11 , 26 ]

(a) Time evolution of negativity for degenerate transition and for different detuning from the bath frequency, where ω 1 = ω 2 = 6 Ω , λ = 0.4 Ω , and the dimensionless detuning is defined by Δ = ( ω 0 ω 1 )/ Ω . (b) A plane view.

In Fig.  4 , we show the time evolution of entanglement negativity in the case of non-degenerate transition and for different detuning between the transition and the bath frequencies, where the parameters are set as ω 1 = 5 Ω , ω 2 = 7 Ω , and λ = 0.4 Ω . It behaves similarly, i.e., the negativity in the two resonant regimes decays faster and the non-Markovianity in the non-resonant regimes becomes larger.

Time evolution of negativity for non-degenerate transition and for different detuning from the bath frequency, where ω 1 = 5 Ω , ω 2 = 7 Ω , and λ = 0.4 Ω .

4. Interference of non-Markovianity

We now discuss the interference between non-Markovianities induced by different transition channels. We will use the non-Markovian measure [ 12 ] based on the entanglement with an ancillary particle. Suppose a system of interest is initially prepared in a maximally entangled state with an ancillary particle, where only the system is influenced by a noise environment and the ancillary particle is noise-free. Then the quantity for measuring the non-Markovianity of the quantum process is defined as [ 12 ]

where Ė ( t ) denotes the time derivative of some kind of entanglement between the system and the ancillary particle, and the integration is extended over all time intervals in which Ė ( t ) is positive. For two-level systems, one usually measures the entanglement E ( t ) by concurrence, while for higher-dimensional systems, the entanglement negativity may be a suitable choice.

The main and ancillary systems in our problem are respectively atoms A and B, which are initially in the maximally entangled state given by Eq. ( 2 ). For the system of three-level V-type atom, the state at any time t is given by Eq. ( 11 ) and its entanglement negativity Eq. ( 21 ) has already been worked out in the above section. Thus the non-Markovianity of Eq. ( 23 ) can be calculated numerically. In Fig.  5 , we plot the non-Markovianity as a function of the bath frequency for different transition frequencies, where λ = 0.4 Ω . We find that the non-Markovianity (the bar charts) in the resonance regimes ( ω i ω 0 ) is smaller and becomes larger when away from the resonance regimes, which agrees with the analysis of the above section.

In order to see the interference effect of non-Markovianity, we also plot the non-Markovianities that only the first transition channel ω 1 or only the second channel ω 2 is opened, i.e., corresponding to a two-level system (The red and green lines in Figs.  5(a) and 5(b) respectively correspond to the non-Markovianity induced by channels ω 1 and ω 2 . In the case of resonance in Fig.  5(c) , the two lines overlap completely and thus only the red line is shown). The black lines are the sum of the red and green lines. It clearly shows that when the bath frequency is between the two transitions, the non-Markovianity causes constructive interference, and outside this frequency domain, it causes destructive interference. When the two transition frequencies are close to each other, the frequency domain of constructive interference becomes narrower and narrower, and finally disappears for degenerate transitions.

Non-Markovianity as a function of the bath frequency with λ = 0.4 Ω : (a) ω 1 = 4 Ω , ω 2 = 8 Ω ; (b) ω 1 = 5 Ω , ω 2 = 7 Ω ; and (c) ω 1 = ω 2 = 6 Ω . The red (green) lines correspond to the non-Markovianity that only the first (second) transition exists and the black lines are the their sum. Note that the red and green lines in panel (c) overlap completely and the green line is omitted. The bar charts correspond to the non-Markovianity of the real three-level atom.

5. Conclusion

We have presented a paradigmatic method for exactly studying non-Markovian dynamics of a multi-level V-type atom interacting with a zero-temperature bosonic bath. For the special case of a three-level V-type atom, the exact analytical expressions for entanglement negativity have been obtained and the dynamical non-Markovianity has been investigated. We have found that the entanglement negativity decays faster and non-Markovianity is smaller in the resonant or quasi-resonant regions compared with that in the non-resonant case. In addition, narrower spectral width is beneficial for improving the dynamical non-Markovianity and preventing the decay of entanglement. These properties are similar to those in the open two-level quantum systems. A distinctive phenomenon observed in the case of three-level V-type atoms is the quantum interference between the non-Markovianities induced by different transition channels. We have demonstrated the constructive and destructive interferences and the corresponding frequency domains for the two kinds of interferences.

Quantum correlation and interference are very important phenomena in quantum physics. Our method for exactly solving dynamics applies in principle to the case of more higher-dimensional open quantum systems, particularly suitable for the numerical simulation of high-dimensional open systems. We expect that non-Markovian dynamics of high-dimensional open quantum systems will show richer features.

Reference
1 Breuer H P Petruccione F 2007 The Theory of Open Quantum Systems Oxford Oxford University Press
2 Kubota Y Nobusada K 2009 J. Phys. Soc. Jpn. 78 114603
3 Ji Y H Hu J J 2010 Chin. Phys. B 19 060304
4 Shao J 2004 J. Chem. Phys. 120 5053
5 Chin A W Datta A Caruso F Huelga S F Plenio M B 2010 New J. Phys. 12 065002
6 Chin A W Huelga S F Plenio M B 2012 Phys. Rev. Lett. 109 233601
7 Bylicka B Chruściński D Maniscalco S 2014 Sci. Rep. 4 5720
8 Laine E M Breuer H P Piilo J 2014 Sci. Rep. 4 4620
9 Vasile R Olivares S Paris M G A Maniscalco S 2011 Phys. Rev. A 83 042321
10 Tang N Fan Z L Zeng H S 2015 Quant. Inform. Comput. 15 0568
11 Breuer H P Laine E M Piilo J 2009 Phys. Rev. Lett. 103 210401
12 Rivas Á Huelga S F Plenio M B 2010 Phys. Rev. Lett. 105 050403
13 Lu X M Wang X G Sun C P 2010 Phys. Rev. A 82 042103
14 Luo S Fu S Song H 2012 Phys. Rev. A 86 044101
15 Hou S C Liang S L Yi X X 2015 Phys. Rev. A 91 012109
16 Zeng H S Tang N Zheng Y P Wang G Y 2011 Phys. Rev. A 84 032118
17 Xu Z Y Yang W L Feng M 2010 Phys. Rev. A 81 044105
18 Chruściński D Wudarski F A 2015 Phys. Rev. A 91 012104
19 Gu W J Li G X 2012 Phys. Rev. A 85 014101
20 Yin X Ma J Wang X Nori F 2012 Phys. Rev. A 86 012308
21 Mäkelä H 2015 Phys. Rev. A 91 012108
22 Deffner S Lutz E 2013 Phys. Rev. Lett. 111 010402
23 Shen H Z Qin M Yi X X 2013 Phys. Rev. A 88 033835
24 Zou C L Chen X D Xiong X Sun F W Zou X B Han Z F Guo G C 2013 Phys. Rev. A 88 063806
25 He Z Zou J Li L Shao B 2011 Phys. Rev. A 83 012108
26 Zeng H S Zheng Y P Tang N Wang G Y 2013 Quantum Inf. Process 12 1637
27 Xue S B Wu R B Zhang W M Zhang J Li C W Tarn T J 2012 Phys. Rev. A 86 052304
28 Tang N Cheng W Zeng H S 2014 Eur. Phys. J. D 68 278
29 Zhang Y J Yang X Q Han W Xia Y J 2013 Chin. Phys. B 22 090307
30 Wang X Y Ding B F Zhao H P 2013 Chin. Phys. B 22 020309
31 Wang X Y Ding B F Zhao H P 2013 Chin. Phys. B 22 040308
32 Xue P Zhang Y S 2013 Chin. Phys. B 22 070302
33 Xie D Wang A M 2014 Chin. Phys. B 23 040302
34 Liu X Wu W 2014 Chin. Phys. B 23 070303
35 Fan Z L Tian J Zeng H S 2014 Chin. Phys. B 23 060303
36 Wang G Y Tang N Liu Y Zeng H S 2015 Chin. Phys. B 24 050302
37 Liu Y Cheng W Gao Z Y Zeng H S 2015 Opt. Express 23 23021
38 Tian L J Ti M M Zhai X D 2015 Chin. Phys. B 24 100305
39 Zou H M Fang M F 2015 Chin. Phys. B 24 080304
40 Liu B H Li L Huang Y F Li C F Guo G C Laine E M Breuer H P Piilo J 2011 Nat. Phys. 7 931
41 Tang J S Li C F Li Y L Zou X B Gou G C Breuer H P Laine E M Piilo J 2012 Europhys. Lett. 97 10002
42 Peres A 1996 Phys. Rev. Lett. 77 1413
43 Horodečki P 1997 Phys. Lett. A 232 333
44 Bruß D Macchiavello C 2002 Phys. Rev. Lett. 88 127901
45 Cerf N J Bourennane M Karlsson A Gisin N 2002 Phys. Rev. Lett. 88 127902