Arxiv:2103.12042V3 [Quant-Ph] 25 Jul 2021 Approximation)
Total Page:16
File Type:pdf, Size:1020Kb
Unified Gorini{Kossakowski{Lindblad{Sudarshan quantum master equation beyond the secular approximation Anton Trushechkin1, 2 1Steklov Mathematical Institute of Russian Academy of Sciences, Moscow 119991, Russia 2National University of Science and Technology MISIS, Moscow 119049, Russia∗ (Dated: July 27, 2021) Derivation of a quantum master equation for a system weakly coupled to a bath which takes into account nonsecular effects, but nevertheless has the mathematically correct Gorini{Kossakowski{ Lindblad{Sudarshan form (in particular, it preserves positivity of the density operator) and also satisfies the standard thermodynamic properties is a known long-standing problem in theory of open quantum systems. The nonsecular terms are important when some energy levels of the system or their differences (Bohr frequencies) are nearly degenerate. We provide a fully rigorous derivation of such equation based on a formalization of the weak-coupling limit for the general case. I. INTRODUCTION particular, does not preserve positivity of the density op- erator, thus leading to unphysical predictions. Also, this Quantum master equations are at the heart of theory of equation does not have the mentioned thermodynamic open quantum systems [1{3]. They describe the dynam- properties. ics of the reduced density operator of a system interact- Derivation of a mathematically correct quantum mas- ing with the environment (\bath") and are widely used ter equation which takes into account nonsecular effects in quantum optics, condensed matter physics, charge and is actively studied. Several heuristic approaches turning energy transfer in molecular systems, bio-chemical pro- the Redfield equation into an equation of the GKLS form cesses [4], quantum thermodynamics [5,6], etc. The without the (full) secular approximation have been pro- Redfield and Davies quantum master equations are well- posed. One possibility is the time coarse graining [22{27]. known microscopically derived equations for a system Another method is a partial secular approximation fol- weakly coupled to a bath and are crucial for understand- lowed by the approximation of slow variation of the spec- ing many physical phenomena. tral density [28{31]. Also, in certain cases, the so-called The Davies quantum master equation [7,8] is derived local approach is considered as an alternative to the sec- in a mathematically rigorous way and has the Gorini{ ular approximation and is largely debated [32{40]. But, Kossakowski{Lindblad{Sudarshan (GKLS) form [9{11], like the Redfield equation, these equations do not satisfy which guarantees that the corresponding dynamics of all the mentioned thermodynamic properties. Moreover, the reduced density operator is well defined. The Davies it has been explicitly shown that the local master equa- equation also satisfies a number of properties important tion violates the second law of thermodynamics [41]. for thermodynamics: stationarity of the Gibbs state, the detailed balance condition, a covariance law (related to In this paper, we derive a unified master equation for the first law of thermodynamics [12{14]), non-negativity the weak-coupling regime in a mathematically rigorous of the entropy production (i.e., the second law of ther- and systematic way, which leads to the GKLS form and modynamics) [1,5,6, 15{17]. all the desired thermodynamic properties. The derivation However, this equation assumes that all distinct Bohr is based on a rigorous formalization of the weak-coupling frequencies (differences between the energy levels) of the limit for the general (nonsecular) case. The unified equa- system are well separated from each other. In other tion has a simple and intuitive structure similar to that words, all differences between distinct Bohr frequencies of the Davies equation. are much higher than the dissipation rates (the secular arXiv:2103.12042v3 [quant-ph] 25 Jul 2021 approximation). This assumption is not satisfied if some Interestingly, this equation coincides with the refined energy levels or Bohr frequencies are nearly degenerate (thermodynamically consistent) form of the local mas- (but not exactly degenerate), which is often the case in ter equation when the local approach is expected to be physical systems. valid. Thus, we rigorously justify the correct form of the The Redfield master equation [18] does not adopt the local master equation (popular due to its simplicity) and secular approximation and thus is more general. It is complete the results of Ref. [42]. widely used in various physical applications, the role of Note that the general idea was also proposed in Ref. [1] the nonsecular terms is studied, e.g., in Refs. [19{21]. and master equations for a particular system were rig- But, unfortunately, it is not of the GKLS form and, in orously derived in Ref. [43]. General properties of the dynamics of particular models of open quantum systems with nearly degenerate spectrum were also rigorously es- tablished in Refs. [44, 45]. Here we derive a master equa- ∗ [email protected] tion for the general situation. 2 II. REDFIELD EQUATION AND SECULAR are bath correlation functions. APPROXIMATION The matrix γαβ(!; !0) [with two double indices i = (α; !0) and j = (β; !)] is, in general, not positive semidef- Consider the system-bath Hamiltonian inite, so, the equation is not of the GKLS form and de- fines the dynamics which violates positivity. i(!0 !)τ/λ2 H = HS + HB + λHI ; (1) As λ 0, the exponents e − for ! = !0 rapidly oscillate! and the corresponding terms can be6 ne- where HS is the isolated system Hamiltonian, HB is the glected (the secular approximation). After this, the equa- isolated bath Hamiltonian, H = A B is the I α α ⊗ α tion becomes of the first standard GKLS form [1,2] since, interaction Hamiltonian (with Aα being system opera- P for each !, the matrix γαβ(!; !) γαβ(!) (with the in- tors and Bα being bath operators), and λ is a formal dices α and β) is well known to be≡ positive semidefinite. small dimensionless parameter. Let HS have a purely We will refer to this equation as the secular or the Davies discrete spectrum: HS = j "jPj, where "j are the dis- master equation. tinct eigenvalues and Pj are the corresponding eigenpro- P However, as pointed out above, the secular approxima- jectors. The differences "j0 "j are called the Bohr fre- tion is not always valid. If so, the formal mathematical quencies. Denote the set of− all (positive, negative, and F limit in the present form does not correspond to physics. zero) Bohr frequencies. There is no limit in a concrete physical system, but all Let the initial system-bath state be ρ σ , where ρ 0 ⊗ B 0 physical quantities have concrete values. The limit λ 0 is the initial state of the system and σB is a reference is just a mathematical expression of the fact that the! dis- iHB t iHB t state of the bath such that e− σBe = σB for all sipative dynamics caused by the coupling to the bath is t (e.g., a thermal state). The dynamics of the reduced much slower than all other time scales, but this is not density operator of the system is given by the case if two different Bohr frequencies ! and !0 are iHt iHt close to each other. Obviously, nearly degenerate energy ρ(t) = Tr e− (ρ σ )e ; B 0 ⊗ B levels is a particular case. where TrB is the partial trace over the bath. The density operator in the interaction picture isρ ~(t) = iH t iH t III. UNIFIED MASTER EQUATION e S ρ(t)e− S : The derivation of the quantum master equation is based on the idea of the separation of different time A. Derivation scales. In the rigorous derivation, this intuition is formal- ized by the Bogolyubov{van Hove limit: λ 0, t , If we claim that, for a given physical system, a differ- 2 ! ! 1 i(!0 !)τ/λ2 λ t = τ = const (see Refs. [3,7,8, 46]). The standard ence ! !0 is small and, hence, the term e − \physical" derivation [2,3] leads to the Redfield equation is not rapidly− oscillating (with respect to the time scale (we express it in a GKLS-like form [27]): of the dissipation), then, in the formal derivation, the difference !0 ! should be treated as infinitesimal and, d − 2 ρ~(τ) = i[HLS(τ); ρ~(τ)] + γαβ(!; !0) moreover, of the order of λ . This should be explicitly dt − !;!0 α,β formalized in the mathematical language of the limit. X2F X 0 τ 1 Let us express the system Hamiltonian HS as i(! !) 2 e − λ A ρ~(τ)Ay 0 Ay 0 A ; ρ~(τ) ; (2) × β! αω − 2 αω β! H = H(0) + λ2δH ; (7) where S S S i(!0 !) τ (0) λ2 where [H ; δHS] = 0 and all nearly degenerate Bohr HLS(τ) = Sαβ(!; !0)e − Aαωy 0 Aβ!; (3) S (0) !;!0 α,β frequencies in H are exactly degenerate in H . In X2F X S S other words, all distinct Bohr frequencies of H(0) are well (the subindex LS stands for the Lamb shift), S separated. Proportionality of the remaining part to λ2 mathematically expresses the fact that some oscillations Aαω = Pj0 AαPj; [HS;Aαω] = !Aαω; (4) i(!0 !)τ/λ2 − e − occur on the same time scale as the dissipa- j;j0 : " " 0 =! Xj − j tion. (0) γ (!; !0) = Γ (!) + Γ∗ (!0); αβ αβ βα Due to commutativity and since HS may be more (0) 1 degenerate than HS, its spectral decomposition is H = S (!; !0) = Γ (!) Γ∗ (!0) S αβ 2i αβ − βα (0) (0) (0) k "k Pk , where "k are the distinct eigenvalues and 1 (0) i!s each eigenprojector P is either one of P or a sum of Γαβ(!) = ds e Cαβ(s) (5) P k j Z0 several Pj. (0) (0) (we assume that the last integral converges).