PHYSICAL REVIEW A 102, 032207 (2020)
Perturbation theory for Lindblad superoperators for interacting open quantum systems
V. Yu. Shishkov,1,2,3 E. S. Andrianov,1,2,3 A. A. Pukhov,1,2,3 A. P. Vinogradov ,1,2,3 and A. A. Lisyansky 4,5 1Dukhov Research Institute of Automatics, 22 Sushchevskaya, Moscow 127055, Moscow Region, Russia 2Moscow Institute of Physics and Technology, 9 Institutskiy Pereulok, Dolgoprudny 141700, Moscow Region, Russia 3Institute for Theoretical and Applied Electromagnetics, 13 Izhorskaya, Moscow 125412, Moscow Region, Russia 4Department of Physics, Queens College of the City University of New York, Flushing, New York 11367, USA 5The Graduate Center, City University of New York, New York, New York 10016, USA
(Received 4 December 2019; revised 2 May 2020; accepted 10 August 2020; published 8 September 2020)
We study the dynamics of interacting open quantum subsystems. Interactions merge these subsystems into a unified system. It is often believed that the Lindblad superoperator of the unified system is the sum of the Lindblad superoperators of the separate subsystems. Such an approach, however, results in a violation of the second law of thermodynamics. To avoid a cumbersome direct derivation of the correct superoperator from first principles, we develop a perturbation theory based on using Lindblad superoperators of separate subsystems, which are assumed to be known. Using interacting two-level systems as an example, we show that, starting with a certain order of the perturbation theory, the second law of thermodynamics holds. We demonstrate that the theory developed can be applied to the problems of quantum transport.
DOI: 10.1103/PhysRevA.102.032207
I. INTRODUCTION unified system grows exponentially with the number of sub- systems. Therefore, for a unified system, the direct approach Recently, the study of open quantum systems has experi- is impractical. enced rapid growth connected to studies of interacting arrays It is desirable to express the Lindblad superoperator of of qubits [1], ensembles of quantum dots and molecules [2], the unified system through the Lindblad superoperators of and systems of atoms coupled with resonators [3–5]. Such subsystems. The key issue is that the eigenstates of the uni- systems are open in the sense that they interact with the fied system are no longer a direct product of the eigenstates environment. of isolated subsystems. Consequently, the transition opera- The exact dynamics of an open quantum system can tors between eigenstates of the unified system that appear only be obtained for a small number of problems [6–10]. in the Lindblad superoperator are no longer equal to direct Therefore, various approximate methods are used [11–13]. products of the transition operators of subsystems, and the One of the most common is the method of the Lindblad- Lindblad superoperator of the unified system is not equal Gorini-Kossakowski-Sudarshan (LGKS) equation [14,15]. to the sum of the Lindblad superoperators for isolated sub- This method implies that both the system and the reservoir systems [17–20]. However, in many cases it is believed that have separated Hamiltonians and that the coupling between even if the whole Hamiltonian includes the interactions be- them can be expressed as the interaction Hamiltonian. The tween subsystems, all the relaxation processes in the unified derivation of the LGKS equation involves the use of the system are the same as for isolated subsystems [2,21–32]. von Neumann equation for the system-reservoir density ma- Such an approximation is called the local approach. Math- trix and the subsequent elimination of the reservoir degrees ematically, it implies that the Lindblad superoperator of the of freedom. The states of the latter are assumed to be un- unified system is the sum of Lindblad superoperators of the changeable during the evolution of the system. After such isolated subsystems. One might suppose that the local ap- elimination, the evolution operator of the system’s density proach will only lead to minor quantitative differences from matrix includes the Hamiltonian and Lindblad superoperators the exact solution. Unfortunately, this is not the case. In [16]. These Lindblad superoperators are constructed from particular, the local approach may result in the violation of the product of the transition operators between the system the second law of thermodynamics [20,33]. Recently, consid- eigenstates. The transitions occur under the influence of the erable attention has been devoted to a detailed comparison reservoirs. The transition rates are calculated according to the between the dynamics predicted by the local approach and Fermi golden rule [17]. the exact Lindblad superoperators for a number of systems Although the general theory of deducing the LGKS equa- [34–41], including interacting two-level systems [42], har- tion is well developed [16], a direct way of constructing monic oscillators [43], and optomechanical systems [44]. In Lindblad superoperators requires knowledge of the exact [45] it was shown that the second law of thermodynamics system eigenstates and eigenenergies [18,19]. For a unified can be restored if the environment is divided into an ensem- system containing interacting subsystems, the complexity of ble of identically prepared auxiliary systems, which interact the calculation of the eigenstates and eigenenergies of the
2469-9926/2020/102(3)/032207(17) 032207-1 ©2020 American Physical Society V. YU. SHISHKOV et al. PHYSICAL REVIEW A 102, 032207 (2020) sequentially with an individual subsystem for an infinitesi- frequency of the atom, ωk is the frequency of a photon with mally short time. However, for the permanent interaction of the wave vector k and the polarization η, and kη is the the system and reservoirs the fulfillment of the second law is coupling constant between the atom and the photon with not guaranteed. the wave number k and polarization η. The operatorσ ˆ is the In this paper we develop a perturbative method for building transition operator from the excited state to the ground state the Lindblad superoperator of a system comprised of inter- of the two-level atom; the operatora ˆkη is the annihilation acting subsystems. The perturbation theory developed for the operator of the photon with the wave vector k and polarization η ˆ σ † + σ λ Lindblad superoperator is based on the construction of the . In this example the operator Sα is ( ˆ ˆ ), α is k0η0 , and † transition operators that appear in this Lindblad superoperator. ˆα η + η / η η R is k (ˆakη aˆk ) k0 0 , where k0 and 0 are the wave We demonstrate that the transition operators can be found vector and the polarization of the free-space mode for which = from the interaction of the Hamiltonian of the reservoir and k0η0 0, respectively. the system by employing the Heisenberg equation of motion. The von Neumann equation for the density matrixρ ˆ (t )of At each step of the perturbation theory the system of the the unified system and the reservoir follows from the Hamil- Heisenberg equations for the transition operator is linear and tonian (1), closed, can be solved exactly. The developed method ensures the fulfillment of the principles of thermodynamics. This ap- ∂ρˆ (t ) i = [ˆρ(t ), Hˆ + Hˆ + Hˆ ], (3) proach permits us to express the Lindblad operator of the ∂t h¯ S SR R unified system through the linear combination of products of the transition operators of the separate subsystems. Using in- where we introduced the Hamiltonian of the unified system teracting two-level systems as an example, we find the lowest ˆ = ˆ (is) + ε ˆ HS HS W . The exclusion of the reservoir degrees of order of the perturbation series that ensures the fulfillment of freedom according to the standard procedure [16,20,46,47] the second law of thermodynamics. demands the transition of the equation into the interaction representation, with the new density matrix defined as II. PERTURBATION THEORY APPROACH FOR OBTAINING THE LINDBLAD SUPEROPERATORS Hˆ + Hˆ Hˆ + Hˆ ρ˜ˆ (t ) = exp i S R t ρˆ (t )exp −i S R t . (4) We consider a system described by the Hamiltonian h¯ h¯ Hˆ = Hˆ (is) + εWˆ + Hˆ + Hˆ , S SR R (1) The von Neumann equation for the new density matrix ρ˜ˆ (t ) ˆ (is) = ˆ (is) can be obtained from Eq. (3), where HS α HSα is the Hamiltonian of the system com- prised of isolated subsystems, which are enumerated by the α ˆ (is) ∂ρˆ subscript and described by the Hamiltonians HSα ; each ˜ (t ) i = ρˆ t , H˜ˆ α t , αth subsystem interacts with its own reservoir described by ∂ ˜ ( ) SR ( ) (5) t h¯ α the operator HˆRα, so the Hamiltonian for the reservoir of the whole system is Hˆ = Hˆ . The Hamiltonians Hˆ (is) and R α Rα Sα ˆ ˆ where the operator H˜ α (t ) is the Hamiltonian Hˆ α, given by HRα depend only on the subsystem and reservoir degrees of SR SR freedom, respectively. The interaction of the αth subsystem Eq. (2), in the interaction representation ˆ and its reservoir is described by the operator HSRα;soforthe whole system Hamiltonian for the interaction with reservoirs ˆ + ˆ ˆ + ˆ ˆ = HS HR − HS HR ˆ = ˆ H˜ SRα (t ) exp i t HˆSRα exp i t is HSR α HSRα. The interaction between the subsystems is h¯ h¯ governed by the Hamiltonian εWˆ = ε α α Wˆα α , where ε is 1 2 1 2 = λ ˆ / ˆ − ˆ / ˆ / ˆ a small dimensionless parameter. For example, in the case of h¯ α exp(iHSt h¯)Sα exp( iHSt h¯)exp(iHRt h¯)Rα ε interacting two-level systems, is the dimensionless coupling × exp(−iHˆRt/h¯). (6) constant between them. For simplicity, we assume that ε is the same for all subsystems. We assume that the Hamiltonian The standard procedure implies the averaging over the ˆ (is) α HSα of the th isolated subsystem has known nondegenerate reservoir degrees of freedom in Eq. (5)[18,19]. This aver- | (is) ω(is) eigenstates kαi and eigenfrequencies αi enumerated by aging procedure is similar to the well-known rotating-wave the subscript i. approximation in quantum optics [48,49]. The averaging over Usually, it is assumed that the Hamiltonian of the interac- the reservoir degrees of freedom in Eq. (5) leads to the master tion between a system and a reservoir has the form [16,17] equation for the system density matrix of the unified system ρ˜ˆS (t ) = TrR[ρ˜ˆ (t )] [18,19]:
Hˆ α = h¯λαSˆαRˆα, (2) SR ρˆ d ˜ S (t ) = ρ . where λα is the interaction constant and Sˆα and Rˆα are the Lˆex[˜ˆS (t )] (7) dt operators of the subsystem and its reservoir, respectively. α = For example, for a single two-level atom, 1, interacting The influence of the reservoir on the unified system in with the electromagnetic field of free space, the Hamiltoni- Eq. (7) reduces to the Lindblad superoperator Lˆ [ρ˜ˆ (t )]. To (is) † † ex S ans are Hˆ α = h¯ω0σˆ σˆ , HˆRα = ,η h¯ωkaˆ aˆkη, and Hˆ α = S k kη SR derive this Lindblad superoperator Lˆex[ρ˜ˆS (t )], according to † + σ † + σ ω k,η h¯ kη(ˆakη aˆkη )( ˆ ˆ ), where 0 is the transition the standard algorithm [18,19], we should find the explicit
032207-2 PERTURBATION THEORY FOR LINDBLAD … PHYSICAL REVIEW A 102, 032207 (2020)
expression for the operators exp(iHˆSt/h¯)Sˆα exp(−iHˆSt/h¯)in time, the near-field interaction between dipole moments of the form two TLSs, which is the strongest electromagnetic interaction between subsystems at the nanoscale, is ∼ (ˆσ †σˆ + σˆ †σˆ ), ˆ / ˆ − ˆ / 1 2 2 1 exp(iHSt h¯)Sα exp( iHSt h¯) i.e., can be presented in the form (13). In Appendix A we discuss the extension of the theory developed for the general = k |Sˆα|k |k k | exp(i ω t ), (8) i j i j ij case. ij To develop the perturbation approach for determination where ωij = ω j − ωi and ω j are eigenfrequencies corre- of the Lindblad superoperators we represent the Hamiltonian sponding to the eigenstates |k of the Hamiltonian Hˆ . Once ˜ˆ j S HSRα (t )[Eq.(6)] as the series (8) is found, one can apply the standard procedure ˜ˆ = λ ˆ ˆ . and obtain the Lindblad superoperator HSRα (t ) h¯ α Sαij(t )Rα (t ) (14) ij ˆ ρˆ = λ2 ω ˆ , ρˆ ˆ† Lex[˜S (t )] αGα ( αij) S ω ˜ S(t )S ω αij αij ˆ = α ij In Eq. (14) we introduce the operators Rα (t ) (is) exp(iHˆ t/h¯)Rˆα exp(−iHˆ t/h¯) and Sˆ α (t ) = exp(iHˆ t/h¯)Sˆ † R R ij S αij + Sˆ ω ρ˜ˆ (t ), Sˆ , (9) (is) α S ωα − ˆ / ˆ ij ij exp( iHSt h¯), with Sαij defined by Eq. (12), and | (is) we use the fact that the eigenvectors kαi of the where we introduced the operators ˆ (is) Hamiltonian HS of isolated subsystems form the complete basis. The perturbation theory for determination of the Sˆ ω = k |Sˆα|k |k k | (10) αij i j i j Lindblad superoperators developed here is based on searching for the approximate expansion for the operators ω ˆ = ˆ / ˆ(is) − ˆ / and the function Gα ( αij) is defined as Sαij(t ) exp(iHSt h¯)Sαij exp( iHSt h¯) in the form +∞ ˆ = εn ˆ (n) . ω = ˆ† ˆ / ˆ − ˆ / Sαij(t ) Sαij(t ) (15) Gα ( αij) TrR[Rα exp(iHRt h¯)Rα exp( iHRt h¯)] 0 n
× exp(i ωα t )dt. (11) To construct the Lindblad equation it is neces- ij (n) sary to represent the operator Sˆ α (t ) in the form Hereafter we call the Lindblad superoperator (9) the exact ij Sˆ (n) exp(−i ω(n)t ), where the operators Sˆ (n) Lindblad superoperator. Generally, finding the exact Lindblad ω(n) ω(n) ω(n) ω(n) superoperator (9) analytically or even numerically is a very do not depend on time and are some yet unknown difficult problem because the exact Lindblad superoperator frequencies. Once such a representation is found, by applying the standard procedure we can arrive at the series expansion Lˆex[ρ˜ˆS (t )] includes the product of the transition operators Sˆ ω Sˆ ω of the Lindblad superoperator αij. Finding the transition superoperators αij [Eq. (10)] is a very cumbersome problem because it demands the so- +∞ ˆ ρ = εn (n) ρ . lution of the eigenproblem for the unified Hamiltonian HS. Lˆ[˜ˆS (t )] Lˆ [˜ˆ S (t )] (16) At the same time, the complexity of this eigenproblem grows n=0 exponentially with the number of subsystems. To find the coefficients Sˆ (n) and frequencies ω(n),we As an alternative, we develop a perturbation approach for ω(n) ˆ the analytical determination of the Lindblad superoperators. solve the Heisenberg equation for Sαij(t ): For future derivations, we introduce the operators dSˆ α (t ) i ij = ˆ , ˆ . ˆ(is) = (is) ˆ (is) (is) (is) . [HS Sαij(t )] (17) Sαij kαi Sα kα j kαi kα j (12) dt h¯ We assume that there is only pair interaction between sub- The initial condition for Eq. (17) follows from the defini- ˆ systems (which is true for the vast majority of problems in tion of Sαij(t ): quantum optics [49]) such that the interaction Hamiltonian Wˆ (is) Sˆ αij(0) = Sˆα . (18) between subsystems has the form ij We solve Eq. (17) using the perturbation method with the (is) (is) Wˆ = h¯ wα α Sˆ Sˆ . (13) ε ˆ 1 2i1 j1i2 j2 α1i1 j1 α2i2 j2 small parameter . To do this we expand the operator Sαij(t ) α =α 1 2 i1 j1i2 j2 in a power series of ε [Eq. (15)]. We substitute this series into Eq. (17) and collect all the terms with the same order of ε.As For the simpler and clearer presentation we assume that ε the interaction between the subsystems (13) is defined by a result, the terms proportional to the zeroth order of give the same operators that introduce each system-environment the equation for the zeroth order of the perturbation theory interaction in Eq. (2). This assumption is true for the most (see Appendix A), optical systems. For example, the interactions of each two- dSˆ (0) (t ) αij =− ω(is) ˆ (0) , level system (TLS) with dissipative and dephasing reservoirs i αijSαij(t ) (19) † † † † dt are ∼ (ˆσ , bˆk + bˆ σˆ1,2 ) and ∼ σˆ , σˆ1,2(bˆk + bˆ ), re- k 1 2 k k 1 2 k (is) (is) (is) † ω = ω − ω σ σ where αij α j αi . We obtain the initial condition spectively ( 1,2 and 1,2 are lowering and raising operators ˆ ˆ† for Eq. (19) from Eq. (18): between TLS states and bk and bk are annihilation and creation ˆ (0) = ˆ(is). operators of kth reservoir mode, respectively). At the same Sαij(0) Sαij (20)
032207-3 V. YU. SHISHKOV et al. PHYSICAL REVIEW A 102, 032207 (2020)
The solution to Eq. (19) with the initial condition (20) is the zeroth-order Lindblad superoperator. Using the Lindblad superoperator (22), we can write the LGKS equation for the ˆ (0) = ˆ(is) − ω(is) . ρˆ Sαij(t ) Sαij exp i αijt (21) density matrix of the system ˜S (t ), dρ˜ˆ (t ) The direct application of the standard procedure for S = Lˆ (0)[ρ˜ˆ (t )]. (23) ˆ ≈ ˆ (0) dt S Sαij(t ) Sαij(t ) leads to the Lindblad superoperator It can be seen from Eqs. (22) and (9) that the Lindblad (0) ˆ (0) ρˆ = λ2 ω(is) ˆ(0) , ρˆ ˆ(0)† superoperator Lˆ [ρ˜ˆS (t )] coincides with the one obtained with L [˜S (t )] αGα αij Sαij ˜ S (t )Sαij α ij the standard procedure for the isolated subsystems in the absence of the interaction between the subsystems (ε = 0). + ˆ(0) ρˆ , ˆ(0)† , Sαij ˜ S (t ) Sαij (22) The first order of the perturbation theory for Eq. (17) ˆ = may be obtained by the substitution of the series Sαij(t ) ˆ(0) = ˆ(is) λ2 ω(is) εn ˆ (n) where Sαij Sαij and αGα ( αij) is defined by (11). We n Sαij(t ) into Eq. (17) and selecting the terms that are call the Lindblad superoperator (22) the Lindblad superoper- proportional to ε. As a result, we obtain the equation for ˆ (1) ator of the zeroth order of the perturbation theory, or in short Sαij(t ) (see Appendix A), ⎡ ⎤ dSˆ (1) (t ) αij (is) (1) ⎣ (0) (0) (0) ⎦ =−i ωα Sˆ α (t ) + i wα α i j i j Sˆ α (t )Sˆ α (t ), Sˆ α (t ) . (24) dt ij ij 1 2 1 1 2 2 1i1 j1 2i2 j2 ij α1 =α2 i1 j1i2 j2
The initial condition for Eq. (24) follows from the initial condition (18) for the operator Sˆ αij(t ) and takes the form ˆ (1) = . Sαij(0) 0 (25) The solution of Eq. (24) with the initial condition (25) is (see Appendix A) (0) (0) (0) wα α i j i j Sˆα Sˆα , Sˆα Sˆ (1) (t ) = 1 2 1 1 2 2 1i1 j1 2i2 j2 ij 1 − exp − i ω(is) + ω(is) t exp − i ω(is)t . (26) αij (is) (is) α1i1 j1 α2i2 j2 αij ωα + ωα α1 =α2 i1 j1i2 j2 1i1 j1 2i2 j2 Combining the solutions (21) and (26), we obtain the solution to Eq. (17) up to the first order of the perturbation theory (0) (1) Sˆ α (t ) ≈ Sˆ α (t ) + εSˆ α (t ), ij ij ij (0) (0) (0) wα α i j i j Sˆα Sˆα , Sˆα ˆ ≈ ˆ(0) + ε ˆ(1) − ω(is) − ε 1 2 1 1 2 2 1i1 j1 2i2 j2 ij Sαij(t ) Sαij Sαij exp i αijt (is) (is) ωα + ωα α1 =α2 i1 j1i2 j2 1i1 j1 2i2 j2 × exp − i ω(is) + ω(is) + ω(is) t , (27) α1i1 j1 α2i2 j2 αij
where superoperator. Using the Lindblad superoperator (29), we can (0) (0) (0) write the LGKS equation for the density matrix of the system wα α i j i j Sˆα Sˆα , Sˆα Sˆ(1) = 1 2 1 1 2 2 1i1 j1 2i2 j2 ij . (28) as follows: αij ω(is) + ω(is) α =α α α 1 2 i1 j1i2 j2 1i1 j1 2i2 j2 ρˆ d ˜ S (t ) (0) (1) = Lˆ [ρ˜ˆS (t )] + εLˆ [ρ˜ˆS (t )]. (30) The direct application of the standard procedure for the dt ˆ ≈ ˆ (0) + εˆ (1) operator Sαij(t ) Sαij(t ) Sαij(t ) leads to the Lindblad It is important to note that the Lindblad superoperator (0) (1) (0) (0) ρ + superoperator Lˆ [ρ˜ˆS(t )] + εLˆ [ρ˜ˆS(t )], with Lˆ [ρ˜ˆS (t )] de- of the first order of the perturbation theory Lˆ [˜ˆS (t )] (1) ε (1) ρ fined by Eq. (22) and Lˆ [ρ˜ˆS (t )] equal to Lˆ [˜ˆS (t )] is nontrivial in the sense that, unlike the Lindblad superoperator (22), it cannot be represented as a sum of the ˆ (1) ρˆ = λ2 ω(is) ˆ(1) , ρˆ ˆ(0)† L [˜ S (t )] αGα αij Sαij ˜ S (t )Sαij Lindblad superoperators of isolated subsystems. Nevertheless, α (0) (1) ij the Lindblad superoperator Lˆ [ρ˜ˆS (t )] + εLˆ [ρ˜ˆS (t )] is con- + ˆ(1) ρˆ , ˆ(0)† + λ2 ω(is) structed from the product of the transition operators of the Sαij ˜ S (t ) Sαij αGα αij (0) subsystems Sˆα that can be seen from (28) and (29). Together, α ij ij these two facts indicate that the first-order Lindblad superop- × ˆ(0) , ρˆ ˆ(1)† + ˆ(0) ρˆ , ˆ(1)† . Sαij ˜ S (t )Sαij Sαij ˜ S (t ) Sαij erators of the perturbation theory takes into account the cross (29) relaxation of the subsystems. For practical applications it is more convenient to write the (0) We call the Lindblad superoperator Lˆ [ρ˜ˆS (t )] + LGKS equation for the density matrixρ ˆS (t ) rather than for the (1) εLˆ [ρ˜ˆS (t )] the Lindblad superoperator of the first order density matrix in the interaction representation ρ˜ˆS (t ). These of the perturbation theory, or in short the first-order Lindblad two density matrices for the unified system are connected by
032207-4 PERTURBATION THEORY FOR LINDBLAD … PHYSICAL REVIEW A 102, 032207 (2020) the equation that follows from Eq. (4): The relaxation dynamics, described by the LGKS equation (23), may violate the second law of thermodynamics and ˆ ˆ HS HS incorrectly predicts the dynamics of the system [20,33]. ρ˜ˆS (t ) = exp i t ρˆS (t )exp −i t . (31) h¯ h¯ In the next section we show that although the Lindblad superoperator of the lower orders of the perturbation theory This transition from the density matrix ρ˜ˆ (t ) to the density S may violate the second law of thermodynamics, the Lindblad matrixρ ˆ (t ) in the LGKS equation (7) can be easily done S superoperator of the higher order of the perturbation theory because restores the fulfillment of the second law of thermodynamics. HˆS HˆS exp −i t Lˆex[ρ˜ˆS (t )] exp i t = Lˆex[ˆρS (t )], (32) h¯ h¯ III. FULFILLMENT OF THE SECOND LAW OF THERMODYNAMICS which follows from the definitions of the exact Lindblad superoperator Lˆex[ρ˜ˆS (t )] [see Eq. (9)] and the transition op- Below, for simplicity, we consider all the reservoirs to Sˆ ω have zero temperature Tα = 0, which is a good approxima- erators αij [see Eq. (10)]. As a result, the LGKS equation with the exact Lindblad superoperator for the density matrix tion for quantum optics. In this limit, the second law of the ρˆS (t ) takes the form thermodynamics is equivalent to the demand that the energy flow between the system and reservoir is directed from the dρˆS (t ) i = [ˆρS (t ), HˆS] + Lˆex[ˆρS (t )]. (33) system to the reservoir. The energy flow in our case (when dt h¯ ˆ the Hamiltonian of the unified system HS does not depend on The same transition can be done in the LGKS equations time) is (23) and (30) with the zeroth-order Lindblad superoper- dρˆ (t ) ˆ (0) ρˆ = ˆ S . ator L [˜S (t )] and the first-order Lindblad superoperator J Tr HS (39) (0) (1) dt Lˆ [ρ˜ˆS (t )] + εLˆ [ρ˜ˆS (t )]. To do this transition we use the approximate equalities The positive sign of the energy flow J indicates that the energy flows from the reservoir to the system. The negative Hˆ Hˆ − S ˆ (0) ρ S ≈ ˆ (0) ρ , sign of the energy flow J indicates that the energy flows exp i t L [˜ˆ S(t )] exp i t L [ˆS(t )] (34) h¯ h¯ from the system to the reservoir. Therefore, the second law ˆ ˆ HS (0) (1) HS of thermodynamics at zero temperature takes the form exp −i t {Lˆ [ρ˜ˆS (t )] + εLˆ [ρ˜ˆS (t )]} exp i t h¯ h¯ J 0. (40) ≈ ˆ (0) ρ + ε ˆ (1) ρ . L [ˆS (t )] L [ˆS (t )] (35) It is worth emphasizing that if there is no interaction be- tween subsystems then, as was pointed out after Eq. (37), the The difference between the left-hand side and the right- LGKS equation (33) with the exact Lindblad superoperator hand side of the approximate equality (34) is proportional to takes the form ε, which follows directly from Eq. (22) and initial condition dρˆS (t ) i = ρˆ (t ), Hˆ (is) + Lˆ (0)[ˆρ (t )]. (41) (20) (see Appendix B). The difference between the left-hand dt h¯ S S S side and the right-hand side of the approximate equality (35)is This means that in the case of ε = 0 the Hamiltonian of proportional to ε2, which follows directly from Eqs. (29) and (is) the unified system is HˆS = Hˆ and the exact Lindblad su- (27) (see Appendix B). Therefore, we can omit these terms S peroperator Lˆex[ˆρS (t )] is equal to the zeroth-order Lindblad on the right-hand sides of Eqs. (34) and (35). As a result, we ˆ (0) ρ ρ superoperator L [ˆS (t )]. In this case, at zero temperature obtain the LGKS equation for the density matrix ˆS (t ) with it is easy to prove that the second law of thermodynamics the zeroth-order Lindblad superoperator holds. Indeed, the substitution of Eq. (41) and the equality (is) dρˆ (t ) i HˆS = Hˆ into Eq. (39) and the use of the commutation S = ρ , ˆ + ˆ (0) ρ S [ˆS (t ) HS] L [ˆS (t )] (36) ˆ (is), ˆ(0) =− ω(is) ˆ(0) dt h¯ relation [HαS Sαij] h¯ αijSαij leads to and the LGKS equation for the density matrixρ ˆ (t ) with the (is) 2 (is) S J =−2¯h ω λ Gα ω first-order Lindblad superoperator 0 αij α αij α ωαij0 ρ d ˆS (t ) i (0) (1) (0)† (0) = [ˆρ (t ), Hˆ ] + Lˆ [ˆρ (t )] + εLˆ [ˆρ (t )]. (37) × Tr Sˆα Sˆα ρˆS (t ) 0. (42) dt h¯ S S S S ij ij In Eq. (42) the subscript 0 in J means that the energy We note that the LGKS equation (36) is equal to one used 0 flow (39) is calculated when there is no interaction between in the local approach [17–20]. Indeed, the Hamiltonian of (is) (0) (is) subsystems (ε = 0). The definitions of ω and Sˆ are the unified system, Hˆ = Hˆ + εWˆ , contains the interaction αij αij S S given after Eq. (22). The sum in Eq. (42) includes only non- between the subsystems εWˆ . At the same time the Lindblad (is) negative frequencies ω because at zero temperature and superoperator Lˆ (0)[ˆρ (t )] coincides with the one obtained with αij S ω(is) ω(is) the exact standard procedure, as it would be in the absence of negative αij, Gα ( αij) is equal to zero (this means that the interaction between the subsystems, only transitions with the energy flowing out of the system are possible) [17]. Thus, in this case the second law of thermo- (0) Lˆ [ˆρS (t )] = Lˆsubsystem[ˆρS (t )]. (38) dynamics holds (for a rigorous proof of the second law for subsystem arbitrary temperature see [18,19,46]).
032207-5 V. YU. SHISHKOV et al. PHYSICAL REVIEW A 102, 032207 (2020)
When the subsystems interact with each other, it turns out ˆ deph = σ †σ ˆ Hamiltonian HSR ˆ ˆ Rdeph, is dephasing. A reservoir that in the lower order of perturbation theory the second law of ˆ dis = which interacts with a TLS through the Hamiltonian HSR thermodynamics may be violated, whereas the higher order of deph (ˆσ † + σˆ )Rˆ is dissipative. Indeed, the Hamiltonian Hˆ = the perturbation theory restores its fulfillment. Indeed, let us dis SR σ †σ ˆ ˆ (is) = consider the zeroth order of perturbation theory. To obtain the ˆ ˆ Rdeph commutes with the TLS Hamiltonian HS ω σ †σ ˆ dis = σ † + σ ˆ energy flow between the unified system and the reservoirs, one h¯ TLS ˆ ˆ , while the Hamiltonian HSR (ˆ ˆ )Rdis does should substitute Eq. (36) into Eq. (39). As a result we obtain not. (local approach) When the interaction of an isolated subsystem with a dissi- pative reservoir permits transitions between system states and (0) (is) (0) (is) J = Tr Hˆ + εWˆ Lˆ [ˆρS (t )] = J0 + ε J0, (43) ω = < S at least some αij 0, then according to Eq. (8), J0 0, (0) and the second law of thermodynamics holds until ε| J0| where the (0) in J means that the energy flow (39)is | | calculated with the zeroth-order Lindblad superoperator. The J0 . This puts a limit on the perturbation parameter. In this energy flow J on the right-hand side of Eq. (43) is defined case, even using local Lindblad superoperators does not cause 0 a violation of the second law of thermodynamics. by Eq. (42). The energy flow J0 on the right-hand side of Eq. (43)is The situation changes when a subsystem interacts only with a dephasing reservoir [20,33]. In this case, all the fre- = { ˆ ˆ (0) ρ }. ω(is) = J0 Tr W L [ˆS (t )] (44) quencies αij 0 and consequently J0 is equal to zero [see Eq. (42)]. Then the second law of thermodynamics in As mentioned in the preceding section, the zeroth-order zeroth order (local approach) holds only if ε = 0orε J < Lindblad superoperator Lˆ (0)[ˆρ (t )] does not take into account 0 S 0. In a general case, ε J may be positive and the second the cross-relaxation processes between the subsystems. This 0 law of thermodynamics would be violated. The addition of means that the application of the zeroth-order Lindblad super- a dissipative reservoir may not help because now we re- operator Lˆ (0)[ˆρ (t )] leads to the additional “wrong” term ε J S 0 ε| deph + dis| | dis| in the energy flow (43) between the system and reservoirs. We quire that J0 J0 J0 . Thus, if the subsystems interact with each other, ε = 0, the local Lindblad superop- showed after Eq. (42) that the term J0, is nonpositive, but the erator, which does not take into account the cross-relaxation wrong term ε J0 may have a positive or a negative sign. As processes, may lead to a violation of the second law of ther- a result, when ε J0 > 0 and |J0| <ε| J0| the second law of thermodynamics is violated, because the energy starts to flow modynamics [20,33] even if the system interacts with both from the reservoir to the system, and J (0) > 0. dephasing and dissipative reservoirs. For example, the local Let us discuss how to determine the lowest order of the Lindblad superoperator may lead to the pumping of the system perturbation series that ensures the fulfillment of the second by the dephasing reservoir with zero temperature [20,33]. law of thermodynamics. To be more concrete, we consider Application of the first-order Lindblad superoperator ˆ (0) ρ + ε ˆ (1) ρ two types of reservoirs: a dephasing reservoirs that does not L [ˆS (t )] L [ˆS (t )] may fix this problem. Indeed, the substitution of Eq. (37) into the energy flow (39)givesthe change the energies of the isolated subsystems but affects (1) the phases of the nondiagonal elements of the each of the energy flow J calculated with the zeroth-order Lindblad subsystem density matrices and a dissipative reservoirs that superoperator drains the energy from the isolated subsystems. For example, (1) = ˆ (is) + ε ˆ { ˆ (0) ρ + ε ˆ (1) ρ } the reservoir Rˆ interacting with a two-level system described J Tr HS W L [ˆS (t )] L [ˆS (t )] ˆ (is) = ω σ †σ by the Hamiltonian HS h¯ TLS ˆ ˆ (where the operator = J + ε2 J , (45) σˆ =|g e| is the transition operator between the ground |e 1 1 | and the excited g states of the TLS), through the coupling where we introduce the energy flows