PHYSICAL REVIEW A 95, 013847 (2017)
Deriving Lindblad master equations with Keldysh diagrams: Correlated gain and loss in higher order perturbation theory
Clemens Muller¨ * and Thomas M. Stace† ARC Centre of Excellence for Engineered Quantum Systems, School of Mathematics and Physics, The University of Queensland, Saint Lucia, Queensland 4072, Australia (Received 12 August 2016; revised manuscript received 14 November 2016; published 27 January 2017)
Motivated by correlated decay processes producing gain, loss, and lasing in driven semiconductor quantum dots [Phys. Rev. Lett. 113, 036801 (2014); Science 347, 285 (2015); Phys.Rev.Lett.114, 196802 (2015)], we develop a theoretical technique by using Keldysh diagrammatic perturbation theory to derive a Lindblad master equation that goes beyond the usual second-order perturbation theory. We demonstrate the method on the driven dissipative Rabi model, including terms up to fourth order in the interaction between the qubit and both the resonator and environment. This results in a large class of Lindblad dissipators and associated rates which go beyond the terms that have previously been proposed to describe similar systems. All of the additional terms contribute to the system behavior at the same order of perturbation theory. We then apply these results to analyze the phonon-assisted steady-state gain of a microwave field driving a double quantum dot in a resonator. We show that resonator gain and loss are substantially affected by dephasing-assisted dissipative processes in the quantum-dot system. These additional processes, which go beyond recently proposed polaronic theories, are in good quantitative agreement with experimental observations.
DOI: 10.1103/PhysRevA.95.013847
I. INTRODUCTION One problem with relying on a canonical transformation as the basis for a perturbative expansion is that it is tailored to a Microwave-driven double quantum dots (DQD) have specific frame, which emphasizes some processes over others. demonstrated a rich variety of quantum phenomena, including population inversion [1–5], gain [6–8], masing [9–14] and It is therefore not guaranteed to find all dissipative processes Sysiphus thermalization [15]. These processes are well un- that occur at a given order in perturbation theory. derstood in quantum optical systems; however, mesoscopic In this paper we present a systematic approach to derive a electrostatically defined quantum dots exhibits additional Lindblad-type master equation at higher perturbative orders, complexity not typically seen in their optical counterparts, which contains all relevant correlated dissipative processes at arising from coupling to the phonon environment. a given order. We use the Keldysh real-time diagrammatic A notable experimental example of this, which motivates technique to calculate the perturbative series and explicitly our work, is an electronically open DQD system coupled to evaluate all fourth-order diagrams that generate correlated a driven resonator, as illustrated in Figs. 1(a) and 1(b) [6]. dissipation. In doing so, we make explicit the link between Substantial gain in the resonator field was observed when the Keldysh self-energy and the Lindblad superoperators. the DQD is blue-detuned with respect to the resonator and We then apply this technique to the experimental situation capacitively biased to induce substantial population inversion. described above, showing that the additional terms we find The observed gain is attributed to correlated emission of a quantitatively explain anomalous gain and loss profiles in resonator photon and a phonon into the semiconductor medium driven DQD-resonator systems. in which the DQD system is defined. This process ensures This paper is structured as follows: conservation of energy, since the phonon carries the energy Section II introduces the model system that motivates this − difference,h ¯(ωq ωr ), between the energy of the qubit and work; namely, a single artificial two-level system (qubit) cou- the energy of the resonator phonon. pled to a resonator as well as a dissipative environment [6,11]. In some experimental regimes of Ref. [6], the observed Section III derives the Keldysh self-energy superoperator gain is well described by a theory based on a canonical in terms of irreducible Keldysh diagrams and decomposes it transformation to a polaron frame [16]. In this frame, conven- into a sum of Lindblad superoperators to form the Keldysh– tional quantum-optics techniques and approximations (Born– Lindblad master equation. The complete set of fourth-order Markov, secular, etc.) are used to derive dissipative Lindblad Keldysh–Lindblad dissipators, presented in Eq. (20), is the superoperators that are quadratic in both the qubit-phonon central theory result of this paper. bath coupling strength β and in the qubit-resonator coupling j Section IV introduces resonator driving and solves the strength g. However, the same theory fails to describe the steady-state Keldysh–Lindblad master equation for a driven, substantial loss (subunity gain) in other experimental regimes, which strongly suggests that there are additional dissipative DQD-resonator system in a mean-field approximation. processes that are not captured in the polaron frame. Section V shows that the additional Lindblad dissipators that arise in the Keldysh analysis of the Dyson series quantitatively account for the of gain and loss in a driven *[email protected] DQD-resonator system across all experimental regimes of †[email protected] Ref. [6].
2469-9926/2017/95(1)/013847(24) 013847-1 ©2017 American Physical Society CLEMENS MULLER¨ AND THOMAS M. STACE PHYSICAL REVIEW A 95, 013847 (2017)
This paper therefore has the subsidiary objective of making explicit the connection between the Keldysh and Lindblad approaches to open quantum systems. As such, we have provided extensive appendixes giving the technical details of our calculations. As we show, Keldysh perturbation theory is a powerful approach for deriving consistent Lindblad dissipators at arbitrary order. It is especially well suited for the treatment of open quantum systems in which correlated dissipative processes, which go beyond the usual second order in perturbation theory, are significant. As such, in addition to the specific problem of correlated gain and loss in a driven DQD resonator, we foresee immediate applications of this technique in numerous contexts, such as cotunneling in open quantum dots [32], charge transport through chains of FIG. 1. (a) Schematic showing a double quantum dot (DQD) Josephson junctions [33], or the calculation of correlated decay coupled to a cavity resonator mode, and to source and drain leads. rates in multilevel circuit-QED systems [34]. (b) An electron (black dot) tunnels from the source lead to the DQD state |L , which couples to the DQD state |R with matrix element II. DRIVEN DISSIPATIVE RABI MODEL q , and then tunnels to the drain, leaving the DQD in the empty state |∅ . Also shown are interdot bias q and coupling rates L and R to The Rabi Hamiltonian describe a single two-level system metallic leads, with chemical potentials μL and μR. (c) Dissipative or qubit coupled to a harmonic mode. It is one of the processes due to phonon emission (dashed arrows) responsible for main workhorses of modern quantum physics and describes rates in Eq. (22). Each processes is correlated with resonator photon such diverse situations as natural atoms interacting with creation (downward wiggly arrows) or annihilation (upward wiggly visible light [18], superconducting artificial atoms coupled to (ωq ∓ωr ) | | arrows); γ↓± correspond to DQD relaxation from e to g , while microwave resonators [35], or semiconductor double quantum (ωr ) | γϕ− corresponds to DQD dephasing, leaving the populations of e dots coupled to superconducting resonators [36]. Here we and |g unchanged. additionally consider the situation where the two-level system is coupled to an environment, inducing dissipation in the dynamics. Here we adopt the language of a semiconducting We finish in Sec. VI with a discussion of our results; open double quantum dot, but the technique presented here is specifically, the correlated decay rates we obtain, which applicable to many other situations that are described by the describe the dynamics of qubit and resonator. Rabi model. The Keldysh and Lindblad formalisms are each well used The Hamiltonian of a DQD coupled to a superconducting in their respective research communities; however, the link resonator, expressed in the DQD position basis {|R ,|L },is = + + between them is usually not made explicit. Equations of H HS HB HI , with the Lindblad type are the generators of completely posi- = † − (p) + (p) tive, trace-preserving dynamics [17] and a wide variety of HS ωr a a q σz /2 q σx /2, (1) techniques have been developed to analyze the dynamics † HB = ωj b bj , (2) of open quantum systems under the actions of Lindblad modes j j master equations, including the input-output [18] and cascaded = (p) + † + (p) scattering matrix, Lindbladian and Hamiltonian (SLH) for- HI gσz (a a )/2 σz X/2, (3) malisms [19], quantum trajectories and measurement [20–24], (p) =| |−| | and phase-space methods [25]. There are widely used but where σz L L R R is the dipole operator of the (p) =| |+| | somewhat heuristic techniques in quantum optics to derive DQD(alsoreferredtoasaqubit),σx R L L R Lindblad master equations at second order in the bath induces transitions between its two charge states, and a is coupling [25]. Going beyond this approximation requires a an annihilation operator for microwave photons in the cavity well-formulated perturbative approach. mode at frequency ωr . Here the position states indicate the The Keldysh diagrammatic technique [26] is widely presence of a single charge either on the right or on the left used in condensed-matter physics as a tool to calculate dot. The DQD has an asymmetry energy between its position Green’s functions of nonequilibrium systems [27]. It has states of q , and they are tunnel coupled with strength q . been previously applied with great success for calculating The coupling between DQD and resonator is described by the the properties of mesoscopic quantum systems [13,28–31]. dipolar interaction, with coupling strength g. Here, as in the = However, in these previous cases, the resulting master rest of the paper, we use the conventionh ¯ 1. equations are not typically expressed explicitly in the The DQD couples to a bath of environmental modes bj with Lindblad form. Furthermore, these derivations are either eigenfrequencies ωj via the system-bath coupling term limited to second order or include only a small class of † X = β (b + b ). (4) diagrams for which closed-form expressions are available. j j j In contrast, here we develop a systematic approach to ac- j count for all Keldysh diagrams at a particular order, with which For descriptive purposes we will assume the environment to to derive a master equation that is explicitly in Lindblad form. be a bath of phonons because this is generally the dominant
013847-2 DERIVING LINDBLAD MASTER EQUATIONS WITH KELDYSH . . . PHYSICAL REVIEW A 95, 013847 (2017) source of dissipation for charge-based quantum dots [2]. In In particular, L2 includes dispersive and dissipative terms at general, the DQD environment could include other electronic second order in HI, L4 includes dissipative terms at fourth fluctuations, such as charged intrinsic two-level defects [37,38] order in HI, Lres accounts for the resonator decay, and Lleads or charge traps on surfaces [39]. describes a transport bias across the DQD, as in Ref. [6]. In addition to the coupling of the DQD to a bath, the The superoperators Lm are composed of a linear combi- microwave resonator will also be coupled to a dissipative en- nation of coherent evolution processes, expressed as com- vironment leading to loss and spurious excitation of resonator mutators with some effective Hamiltonian, and dissipative photons. Also, in the experiments motivating this work, the evolution, expressed in the form of Lindblad dissipators quantum dot was subject to an external bias, leading to charge † † † transport through the DQD as illustrated in Fig. 1(b).Both D[oˆ]ρ = oρˆ oˆ − (oˆ oρˆ + ρoˆ oˆ)/2. (9) these processes will be included in the master equation and we describe them in Secs. III B and III C. For the purposes of this paper, we will refer to any term of the We write the interaction Hamiltonian HI as a sum of form in Eq. (9) as a Lindblad dissipator, although we note that, the intrasystem interaction HI,S (i.e., between the qubit and at higher order, we find such terms with negative coefficients. resonator), and the interaction between the system and phonon Keldysh perturbation theory provides a systematic way to environment HI,B, evaluate self-energy terms in the Dyson equation, in order to derive the dynamical master equation for the system density = + HI HI,S HI,B matrix. To determine the system evolution, we evaluate the 1 Keldysh self-energy superoperator , which is written as a = g(cos θσ + sin θσ)(a + a†) 2 z x perturbation series in the interaction Hamiltonian. Terms in this series are expressed in the graphical language of Keldysh 1 † + (cos θσ + sin θσ) β (b + b ), (5) diagrams in order to keep track of all relevant contributions 2 z x j j j j and avoid double counting. To begin the analysis, the Keldysh master equation for the (p) = + where we have decomposed σz cos θσz sin θσx into the density matrix in the Schrodinger¨ picture, ρ(S)(t), is given by =| |+| | =| |−| | eigenoperators, σx e g g e and σz e e g g of the bare system-bath Hamiltonian t (S) (S) (S) ρ˙ (t) =−i[H0,ρ (t)] + dt1ρ (t1) (t1,t), (10) H = H + H (6) t0 0 S B = † − + † ωr a a ωq σz/2 ωj bj bj , (7) where the superoperator is the self-adjoint self-energy, j which acts on the state from the right. A derivation of Eq. (10) can be found in Appendix A1. In the interaction picture defined where the bare qubit level splitting is ω = (2 + 2 )1/2 and q q q by H , this becomes the qubit mixing angle tan θ = / , with the convention 0 q q 0 θ<π. The qubit eigenstates are given by t ρ˙(I)(t) = dt ρ(I)(t ) (t ,t). (11) θ θ 1 1 1 |g = cos |L − sin |R , t0 2 2 θ θ Hereafter, we drop the interaction-picture label (I). In what |e = sin |L + cos |R . 2 2 follows, we seek to express the self-energy superoperator on the right-hand side (RHS) of Eq. (11) in the form of dispersive In the following we perform perturbation theory in the terms (i.e., commutators with the Hamiltonian) and dissipative interaction Hamiltonian HI. We are thus assuming that terms (i.e., Lindblad superoperators), as in Eq. (8). the qubit resonator coupling is weak compared with their detuning, g<|ωq − ωr |. Additionally, we make the usual weak-coupling approximation between the system and the A. Keldysh self-energy phonon environment. The standard situation in experiments As described in Appendix A1, the self-energy superoper- is ω>g βj , where ω represents all relevant system and ator (t1,t) is given by the sum of irreducible superoperator- bath frequencies. valued Keldysh diagrams, each of which evolves the state from an early time t1 to the current time t. We evaluate (t1,t)to III. LINDBLAD SUPEROPERATORS FROM KELDYSH a given order of perturbation theory by truncating the sum at SELF-ENERGY that order. The objective of this paper is to derive a Lindblad master Keldysh superoperators are represented diagrammatically equation as two parallel lines representing the time evolution of the forward (·|) and reverse (|· ) components of the density ρ˙ = Lmρ matrix. Vertices placed on the lines correspond to the action m of the interaction Hamiltonian, so at mth order of perturbation theory there are m interaction vertices. These act at specific = L2ρ + L4ρ + Lresρ + Lleadsρ, (8) interaction times tj over which we integrate. Tracing over for the evolution of the (slowly varying) system density matrix bath degrees of freedom corresponds to contracting vertices ρ, where the index m labels orders of perturbation theory. together, indicated by dashed lines.
013847-3 CLEMENS MULLER¨ AND THOMAS M. STACE PHYSICAL REVIEW A 95, 013847 (2017)
To evaluate the time integral appearing in the Keldysh (a) master equation (11), we note that each Keldysh diagram takes the form of a convolution kernel with oscillatory time- dependence, so that the Laplace transform in the time-domain becomes a product of terms in Laplace space, i.e., ρ(t) → ρs . In Laplace space, Eq. (11) becomes (b) − = ¯ sρs ρ(0) ρs+iωj s , (12)
ωj where ωj labels a set of system frequencies. One approach to the solution of these equations was presented in Ref. [4] starting from the ansatz FIG. 2. (a) Diagrams representing all second-order terms in the Keldysh self-energy. The red dashed line indicates that bath operators ρ¯0 ρ¯j in the vertices connected by it are contracted together to form a ρs = + . (13) s s + iωj bath spectral function. (b) Two examples of irreducible fourth-order ωj Keldysh diagrams. All 32 irreducible diagrams at this order can be In brief, the unknown residues in Eq. (13) are found by generated from the two shown by swapping subsets of interaction vertices between the lower and upper lines. comparing the poles on both sides of Eq. (12) and imposing a self-consistency condition thereupon. This self-consistent solution ultimately truncates the summation in Eq. (13), which these diagrams in Appendix A4. At second order we write, corresponds to a secular or rotating-wave approximation. L = L − Further details are given in Appendix A3. 2ρ¯ 2,dissρ¯ i[H2,ρ¯], (15) In the experiments that motivate this work, the steady-state where we specify L2,diss and H2 below. response of the system determines the phenomenology, so we Second-order dissipative terms. If the vertices in Fig. 2(a) focus here on the slow dynamics of the system. represent the interaction with the underlying phonon bath Retaining only the pole at s = 0 corresponds to the conven- HI,B, the diagrams yield the usual dissipative terms in tional secular approximation made in deriving the quantum the master equation for a two-level system coupled to an optical master equation [4], so we use the quasistatic ansatz environment [18,40]. They evaluate to (cf. Appendix A4a) ρ = ρ/s,¯ (14) s L2,dissρ¯ = γ↓,2D[σ−]¯ρ + γ↑,2D[σ+]¯ρ + γϕ,2D[σz]¯ρ, (16) so thatρ ¯ captures the quasistatic evolution associated with where the second-order rates are the weak coupling to the bath. 2 γ↓ = θCω / , We now describe the results of evaluating Lm at each order. ,2 sin ( q ) 2 2 γ↑,2 = sin θC(−ωq )/2, (17) 1. First order 2 γϕ,2 = cos θC(0)/2, There are no dissipative terms surviving the perturbative treatment at first order in perturbation theory, since for any and in which the spectral function of the phonon environment † thermal state of the environment one finds b =b = is j th j th 0[18,25]. This is generally true for any odd moments of bath ∞ iω(t1−t2) operators. C(ω) = dte Xˆ (t1)Xˆ (t2) −∞ In our case, the perturbative interaction Hamiltonian HI 1 also contains coherent terms between the qubit and resonator, = J (ω)[n (ω) + θ(ω)], i.e., terms that couple two parts of the system to each other. 2 th In this case the above reasoning no longer applies and some − † ˆ = iωj t + iωj t terms might survive even in first order of perturbation theory. where X(t) j βj (bj e bj e ) is the bath-coupling However, these terms carry a rapidly rotating time dependence. operator in the interaction picture, J (ω) is the spectral density, † − = = βωj − 1 Making a rotating-wave approximation [consistent with our and bj bj nth(ωj ) (e 1) is the Bose distribution = approximate ansatz, Eq. (14)] eliminates these terms also. It with β 1/kB T . θ(ω) is the step function. See Appendix A2 follows that L1ρ¯ = 0. for more details. Second-order coherent terms. If the vertices in Fig. 2(a) represent dot-resonator interactions from H , most diagrams 2. Second order I,S carry fast oscillatory terms, so that they vanish in a rotating- In second-order perturbation theory there are four classes wave approximation. The residual, stationary terms contribute of diagrams to consider, depicted in Fig. 2(a). Evaluating these to the imaginary part of the self-energy, leading to a renormal- diagrams generates a dissipative contribution arising from ization of the Hamiltonian. We find (cf. Appendix A4b) system-phonon interactions, as well as a coherent contribution ω from the dispersive, coherent terms arising from qubit- H = χ q σ + a†aσ , 2 + ( z 2 z) (18) resonator coupling. We provide worked examples evaluating ωq ωr
013847-4 DERIVING LINDBLAD MASTER EQUATIONS WITH KELDYSH . . . PHYSICAL REVIEW A 95, 013847 (2017)
2 2 where the dispersive shift is χ = g sin θ/(4ωq − 4ωr ), An example diagram at fourth order is evaluated in similar to the standard treatment of the Jaynes–Cummings Appendix A5afor pedagogical purposes. After evaluating all model in the dispersive regime [35]. The difference between the relevant fourth-order Keldysh diagrams, we are left with our result (18) and the usual Jaynes–Cummings dispersive a large number of superoperators that contribute to L4,corrρ¯. shift lies in the inclusion of counter-rotating terms in the Collecting these into Lindblad dissipators requires careful system-bath coupling (5), leading to the additional ratio of analysis, which we describe in Appendix A6. system frequencies in Eq. (18). Our result does not make The full set of correlated photon-phonon decay terms the rotating-wave approximation in this coupling and thus that appear at fourth order of Keldysh–Lindblad perturbation includes both the usual dispersive as well as the Bloch–Siegert theory are shift [41,42] in one analytic expression. L4,corrρ¯ † † 3. Third order = ↓+D[σ−a ]¯ρ + ↓−D[σ−a]¯ρ + ↑+D[σ+a ]¯ρ † Third-order diagrams can again contribute due to the + ↑−D[σ+a]¯ρ + ϕ+D[σza ]¯ρ + ϕ−D[σza]¯ρ † system-system coupling terms in our perturbative interaction + −D[a]¯ρ + +D[a ]¯ρ + ↑↓(D[σ+]¯ρ + D[σ−]¯ρ) ∼ 3 ∼ 2 Hamiltonian. Their contributions can be at order g or gβ , † † + n(D[a a]¯ρ − D[σz + a a]¯ρ) + ϕ,4D[σz]¯ρ i.e., they are at even order in HI,B and odd order in HI,S. These † † terms contribute Hamiltonian corrections, i.e., small energy + ±,ϕ±(D[a + σza]¯ρ + D[a + σza ]¯ρ) shifts and changes of the bare system parameters, so we ignore † † + (D[σ a a]¯ρ − D[σ + σ a a]¯ρ) them here. Instead, we assume that these contributions are ϕn z z z † † implicitly included in renormalized experimental parameters. + ↓n(D[σ−a a]¯ρ − D[σ− + σ−a a]¯ρ) Thus we take L3ρ¯ = 0. † † + ↑n(D[σ+a a]¯ρ − D[σ+ + σ+a a]¯ρ). (20)
4. Fourth order Each of the rates j are sums of terms depending on the bath spectral function, C(ω), evaluated at frequencies ω ∈{0, ± In fourth-order perturbation theory, there is a total of 32 ± ± − ± + } irreducible diagrams. Two examples are depicted in Fig. 2(b) ωq , ωr , (ωq ωr ), (ωq ωr ) , cf. Table II. Addi- tionally, there are contributions to the j that are proportional and the full set is shown in Appendix A5. These 32 diagrams = to the derivative of the bath spectral function, C (ω), evaluated all have four vertices, each of which represents terms in HI =± H + H . When evaluating diagrams, we first expand each at frequency ω ωq . Those originate from diagrams where I,S I,B the time evolution involves degenerate frequencies, leading to one into all possible combinations of HI,S and HI,B over the vertices. Within this expansion, terms with an odd number of double poles in Laplace space, as has been noticed before [43]. Their physical interpretation is uncertain. One possibility is instances of HI,B vanish, since odd moments of bath operators vanish. Terms with an even number of instances of both H that they are first-order terms in a Taylor expansion of the I,B bath spectral function and as such indicative of an implicit and HI,S are further classified into three distinct categories of diagrams: renormalization of system frequencies. We have verified that (1) Correlated photon-phonon decay processes: diagrams the derivative terms do not contribute significantly in the in which two vertices come from the intrasystem interaction experimental parameter regime we consider below and leave Hamiltonian H and the remaining two are from the system- their detailed understanding to future work. I,S Analyzing all 21 dissipators above in detail is not the main bath interaction HI,B. (2) Coherent dispersive photon-photon processes: dia- concern of this paper. However, the first six dissipators in grams in which all four vertices originate from the coherent Eq. (20) are easy to interpret and have dominant contributions γj to the corresponding rates j , which can be expressed intrasystem interaction HI,S. (3) Two phonon decay processes: diagrams in which all simply. These are four vertices come from the system-bath interaction term HI,B. (0) (ωq −ωr ) † (ωq +ωr ) L ρ¯ = γ D[σ−a ]¯ρ + γ D[σ−a]¯ρ We thus write 4,corr ↓+ ↓− (−ωq −ωr ) † (−ωq +ωr ) + D + + D + L4ρ¯ = L4,corrρ¯ − i[H4,ρ¯] + L4,2phononρ,¯ (19) γ↑+ [σ a ]¯ρ γ↑− [σ a]¯ρ − + (ωr )D + ( ωr )D † where the order of the terms corresponds to the numbered list γϕ− [σza]¯ρ γϕ+ [σza ]¯ρ, (21) above. where ↑ and ↓ denote qubit flipping processes, ϕ indicates In this paper we are only concerned with calculating the qubit dephasing, and + and − photon excitation and loss, self-energy diagrams leading to correlated photon-phonon respectively. The rates are decay processes, L4,corr. We therefore evaluate neither H4 nor 2 2 L4,2phonon here, but leave this to future work. In passing, we − 1 ω sin θ γ (ωq ωr ) = g2 cos2 θ q C(ω − ω ), note that H4 will renormalize the bare Hamiltonian terms, and ↓+ 2 2 q r 2 ω (ωq − ωr ) L will renormalize the rates in Eq. (17), and potentially r 4,2phonon 2 2 + 1 ω sin θ include dispersive shifts as well. γ (ωq ωr ) = g2 cos2 θ q C(ω + ω ), (22) ↓− 2 + 2 q r Fourth-order correlated photon-phonon decay terms. 2 ωr (ωq ωr ) Evaluating correlated photon-phonon decay self-energy di- ω2 sin2 θ (ωr ) 1 2 2 q agrams at fourth order results in a total of 21 individual γ − = g sin θ C(ω ). ϕ 2 (ω − ω )2(ω + ω )2 r dissipators, each associated with rates j . q r q r
013847-5 CLEMENS MULLER¨ AND THOMAS M. STACE PHYSICAL REVIEW A 95, 013847 (2017)
The remaining three rates in Eq. (21) are related by thermal We therefore apply a displacement transformation to the occupation factors, resonator modes [45,46], a → a˜ + α. This transformation (−ω +ω ) (ω −ω ) − − forms the basis of a semiclassical, mean-field approximation q r = q r β(ωq ωr ) γ↑− γ↓+ e , for dissipative steady-state equations for qubit and resonator, (−ω −ω ) (ω +ω ) − + from which we establish a self-consistency condition that γ q r = γ q r e β(ωq ωr ), ↑+ ↓− determines α. − − ( ωr ) = (ωr ) βωr The canonical unitary for the displacement transformation γϕ+ γϕ− e . is In Sec. V, we analyze in detail the contributions from Eqs. (22) = { † − ∗ } to gain and loss in the parameter regime of recent experiments D(α) exp αa α a , (26) in semiconductor double quantum dots coupled to a microwave so that D(α)†aD(α) = a˜ + α. The Hamiltonian in the dis- resonator. There we show that the six dissipators in Eq. (22) placed frame is transformed as quantitatively explain recent experimental results [6], whereas † † including the full set of rates from Eq. (20) does not change H˜ = D(α) HD(α) − iD(α) D˙ (α) the picture qualitatively apart from a change in numerical † † ∗ = D(α) HD(α) − i(a˜ α˙ − a˜α˙ ), (27) parameters. where α may be time dependent. B. Resonator decay When applying the displacement operation, we have to take care to include dissipative processes containing resonator The resonator couples weakly to the external electromag- operators in the transformation. For example, when applying netic environment [44], and we include this as an additional the transformation to the dissipator describing direct photon decay process in the usual way [25] through the Lindblad decay ∼D[a]ρ we find dissipators D † =D − ∗ − † † [D aD]¯ρ [a˜]¯ρ i[i(α a˜ αa˜ )/2,ρ¯], (28) Lresρ¯ = κ−,r D[a]¯ρ + κ+,r D[a ]¯ρ, (23) where the second term describes a coupling between the with the resonator relaxation and excitation rates κ∓,r and the residual mode a˜ and the coherent field amplitude α.This = − bare resonator linewidth κ κ−,r κ+,r . term, and similar ones from other dissipators, contributes to a self-consistent equation for the coherent field amplitude α. C. External quantum dot bias Correlated dissipators transform similarly. For example, in D † In the experiments that motivate this work, the open the displaced frame, [σ−a ]ρ transforms as † † † 2 DQD coupled to external leads, inducing a charge-discharge D[D σ−a D]¯ρ = D[σ−a˜ ]¯ρ + |α| D[σ−]¯ρ transport cycle, illustrated in Fig. 1(b). To model this process, † ∗ we extend the DQD basis to include the empty state |Ø ,in + ασ−a˜ ρσ¯ + + α σ−ρ¯aσ˜ + which the DQD is uncharged. As electrons tunnel between the 1 † ∗ − (ασ+σ−a˜ ρ¯ + α aσ˜ +σ−ρ¯ leads and the DQD, it passes transiently through the empty 2 † ∗ state. This process is described by the incoherent Lindblad + αρσ¯ +σ−a˜ + α ρ¯aσ˜ +σ−). (29) superoperator [2,17,18,21,22] A. Separable and semiclassical approximations Lleadsρ = LD[|L Ø|]ρ + RD[|Ø R|]ρ. (24) In principle, we could solve the steady-state (or dynamical For simplicity, we assume = = in the following. R L behavior) of the master equation (8), accounting fully for Depending on the sign of and the strength of , the DQD q correlations between the resonator and qubit. However, to population may become inverted in the steady state. make analytical progress we now make a simplifying mean- field approximation that enables us to calculate the qubit and IV. DRIVEN RESONATOR resonator steady states. The additional dissipators tend to drive the system to its That is, we assume the resonator to be in the semiclassical | ground state. To see nontrivial behavior arising from the master regime, so that it is well approximated by a coherent state α . equation, the system needs to be driven out of equilibrium. We will self-consistently determine the coherent amplitude α Depending on the specific system, there are a variety of ways by finding a displaced resonator frame which is effectively to achieve this. Here we consider resonator driving, in which undriven. an external microwave field is imposed on the resonator. To In this approximation, the resonator and qubit are separable, this end we introduce an additional resonator driving term in ρ¯ = ρ¯r ⊗ ρ¯q . (30) the system Hamiltonian, H = HS + HD, where † − Under this approximation, Eq. (8) decomposes into two = iωd t + iωd t HD d (a e ae )/2, (25) coupled mean-field equations: one for the resonator state, ρ¯r =|α α|, and one for the qubit state,ρ ¯q , in which each with the amplitude of the drive d and frequency ωd . In writing Eq. (25) we have assumed a rotating-wave approximation, subsystem experiences the mean field of the other. We find ≈ relevant to near-resonant driving ωd ωr . Under driving, we anticipate that the resonator will tend ρ¯˙ = Tr D†L Dρ¯ , (31) | r q m to a steady state which is close to a coherent state α . m
013847-6 DERIVING LINDBLAD MASTER EQUATIONS WITH KELDYSH . . . PHYSICAL REVIEW A 95, 013847 (2017)