Keldysh Field Theory for Driven Open Quantum Systems
L. M. Sieberer,1 M. Buchhold,2 and S. Diehl2, 3, 4 1Department of Condensed Matter Physics, Weizmann Institute of Science, Rehovot 7610001, Israel 2Institute of Theoretical Physics, TU Dresden, D-01062 Dresden, Germany 3Kavli Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106, USA 4Institut f¨urTheoretische Physik, Universit¨atzu K¨oln,D-50937 Cologne, Germany Recent experimental developments in diverse areas - ranging from cold atomic gases over light-driven semi- conductors to microcavity arrays - move systems into the focus, which are located on the interface of quantum optics, many-body physics and statistical mechanics. They share in common that coherent and driven-dissipative quantum dynamics occur on an equal footing, creating genuine non-equilibrium scenarios without immediate counterpart in condensed matter. This concerns both their non-thermal flux equilibrium states, as well as their many-body time evolution. It is a challenge to theory to identify novel instances of universal emergent macro- scopic phenomena, which are tied unambiguously and in an observable way to the microscopic drive conditions. In this review, we discuss some recent results in this direction. Moreover, we provide a systematic introduction to the open system Keldysh functional integral approach, which is the proper technical tool to accomplish a merger of quantum optics and many-body physics, and leverages the power of modern quantum field theory to driven open quantum systems.
CONTENTS II Applications 32
I. Introduction2 III. Non-equilibrium stationary states: spin models 32 A. New phenomena2 A. Ising spins in a single-mode cavity 32 B. Theoretical concepts and techniques3 1. Large-spin Holstein-Primakoff representation 33 C. Experimental platforms5 2. Effective Ising spin action for the single-mode cavity 35 1. Cold atoms in an optical cavity, and microcavity arrays: driven-dissipative spin-boson models5 B. Ising spins in a random multi-mode cavity 36 1. Keldysh action and saddle point equations 37 2. Exciton-polariton systems: driven open interacting bosons7 2. Non-equilibrium glass transition 38 3. Thermalization of photons and atoms in the glass 3. Cold atoms in optical lattices: heating phase 39 dynamics9 4. Quenched and Markovian bath coupling in the D. Outline and scope of this review 10 Keldysh formalism 40
IV. Non-equilibrium stationary states: bosonic models 40 I Theoretical background 11 A. Density-phase representation of the Keldysh action 42 II. Keldysh functional integral for driven open systems 11 B. Critical dynamics in 3D 44 A. From the quantum master equation to the Keldysh 1. Effective action for driven-dissipative functional integral 11 condensation 45 B. Examples 16 2. Non-equilibrium FRG flow equations 46 1. Single-mode cavity, and some exact properties 16 3. Scaling solutions and critical behavior 48 2. Driven-dissipative condensate 20 C. Absence of algebraic order in 2D 49 C. Semiclassical limit of the Keldysh action 21 D. KPZ scaling in 1D 51 D. Symmetries of the Keldysh action 23 1. Thermodynamic equilibrium as a symmetry of V. Universal heating dynamics in 1D 52 the Keldysh action 24 A. Heating an interacting Luttinger liquid 53 2. Semiclassical limit of the thermal symmetry 26 B. Kinetic equation 54 3. The Noether theorem in the Keldysh formalism 27 C. Self-consistent Born approximation 56 4. Closed systems: energy, momentum and particle 1. Kinetics for small momenta 58 number conservation 28 2. Scaling of the self-energy 58 5. Extended continuity equation in open systems 29 D. Heating and universality 58 6. Spontaneous symmetry breaking and the 1. Phonon densities 58 Goldstone theorem 30 2. Non-equilibrium scaling 59 E. Open system functional renormalization group 31 3. Observability 59 2
VI. Conclusions and outlook 60 techniques to be used, and calls for the development of new theo- retical tools. The physical framework sparks broader theoretical VII. Acknowledgements 60 questions on the existence of new phases of bosonic [18, 19] and fermionic [20–22] matter, the nature of phase transitions in such A. Functional differentiation 61 driven systems [10, 23–26], and the observable consequences of quantum mechanics at the largest scales [27]. Beyond stationary B. Gaussian functional integration 61 states, a fundamental challenge is set by the time evolution of interacting quantum systems, which is currently explored theo- References 62 retically [28–33] and experimentally in cold atomic [34–38] and photonic systems [39]. A key goal is to identify universal dy- namical regimes that hold beyond specific realizations or precise I. INTRODUCTION initial conditions. Combining idealized closed system evolution with the intrinsic open system character of any real world ex- periment takes this setting to the next stage, and exhibits emer- Understanding the quantum many-particle problem is one of gent dynamics markedly different from closed systems both for the grand challenges of modern physics. In thermodynamic short [40, 41] and long evolution times [42–47]. equilibrium, the combined effort of experimental and theoretical The interplay of coherent and driven-dissipative dynamics can research has made tremendous progress over the last decades, be a natural consequence of the driving necessary to maintain revealing the key concepts of emergent phenomena and univer- a certain many-body state. Going one step further, it is possi- sality. This refers to the observation that the relevant degrees ble to exploit and further develop the toolbox of quantum optics of freedom governing the macrophysics may be vastly differ- for the driven-dissipative manipulation of many-body systems. ent from those of the microscopic physics, but on the other Recently, it has been recognized that the concept of dissipative hand are constrained by basic symmetries on the short distance state preparation in quantum optics [48, 49] can be developed scale, restoring predictive power. Regarding the role and power into a many-body context, both theoretically [50–53] and exper- of these concepts in out-of-equilibrium situations, there is a imentally [54–56]. Suitably tailored dissipation then does not large body of work in the context of classical near equilibrium necessarily act as an adversary to subtle quantum mechanical and non-equilibrium many-body physics and statistical mechan- correlations such as phase coherence or entanglement [57–62]. ics [1–5]. However, analogous scenarios and theoretical tools In contrast, it can even create these correlations, and dissipa- for non-equilibrium quantum systems are much less developed. tion then represents the dominant resource of many-body dy- This review addresses recent theoretical progress in an impor- namics. In particular, even topologically ordered states in spin tant and uprising class of dynamical non-equilibrium phases of systems [63] or of fermionic matter [64, 65] can be induced dis- quantum matter, which emerge in driven open quantum systems sipatively ([52] and [66–68], respectively). These developments where a Hamiltonian is not the only resource of dynamics. This open up a new arena for many-body physics, where the quan- concerns both non-equilibrium stationary states, but also the dy- tum mechanical microscopic origin is of key importance despite namics of such ensembles. Strong motivation for the exploration a dominantly dissipative dynamics. of non-thermal stationary states comes from a recent surge of Summarizing, this is a fledging topical area, where first results experiments in diverse areas: in cold atomic gases [6–8], hy- underpin the promise of these systems to exhibit genuinely new brid light-matter systems of Bose-Einstein condensates placed physics. In the remainder of the introduction, we will discuss in optical cavities are created [9, 10], or driven Rydberg ensem- in some more detail the three major challenges which emerge in bles are prepared [11]; light-driven semiconductor heterostruc- these systems. The first one concerns the identification of novel tures realize Bose condensation of exciton-polaritons [12, 13]; macroscopic many-body phenomena, which witness the micro- coupled microcavity arrays are pushed to a regime demonstrat- scopic driven open nature of such quantum systems. Second, we ing both strong coupling of light and matter [14] and scalabil- anticipate the theoretical machinery which allows us to perform ity [15]; large ensembles of trapped ions implement varieties of the transition from micro- to macrophysics in a non-equilibrium driven spin models [16, 17]. All those systems are genuinely context in practice. Third, we describe some representative ex- made of quantum ingredients, and share in common that coher- perimental platforms, which motivate the theoretical efforts, and ent and driven-dissipative dynamics occur on equal footing. This in which the predictions can be further explored. creates close ties to typical setups in quantum optics. But on the other hand, they exhibit a continuum of spatial degrees of free- dom, characteristic for many-body physics. Systems located at this new interface are not guaranteed to thermalize, due to the A. New phenomena absence of energy conservation and the resulting breaking of detailed balance by the external drive at the microscale. They As pointed out above, one of the key goals of the research re- rather converge to flux-equilibrium states of matter, creating viewed here is the identification of new macroscopic many-body scenarios without counterpart in condensed, equilibrium matter. phenomena, which can be uniquely traced back to the micro- This rules out conventional theoretical equilibrium concepts and scopic driven open nature of such systems, and do not have an 3 immediate counterpart in equilibrium systems. cal behavior is characterized by a fine structure in a new and in- The driven nature common to the systems considered here can dependent critical exponent, which measures decoherence, and always be associated to the fact that the underlying microscopic whose value distinguishes equilibrium and non-equilibrium dy- Hamiltonian is time dependent, with a time dependence relating namics, see Sec.IVB. to external driving fields such as lasers. When such an ensemble Instead of emergent thermal behavior indicating the fadeout of (i) in addition has a natural partition into a “system” and a “bath” non-equilibrium conditions upon coarse graining, also the oppo- — a continuum of modes well approximated by harmonic oscil- site behavior is possible. For example, low dimensional (d ≤ 2) lators with short memory —, (ii) the system-bath coupling is bosonic systems at low noise level, such as exciton-polaritons weak compared to a typical energy scale of the “system” Hamil- well above threshold, are not attracted to the equilibrium fixed tonian and (iii) linear in the bath operators, then an effective point as their three dimensional counterparts, but rather flow to (still microscopic) description in terms of combined Hamilto- the non-equilibrium fixed point of the Kardar-Parisi-Zhang [69] nian and driven-dissipative Markovian quantum dynamics of the universality class, see Sec.IVC. This can be interpreted as a “system” ensues. The “system” dynamics obtains by tracing out universal and indefinite increase of the non-equilibrium strength, the bath variables. which is triggered even if the violation of detailed balance at the The resulting effective microdynamics is “non-equilibrium”, microscopic level is very small. in a sense sharpened Sec. II D 1. This not only concerns the Universal non-equilibrium phenomena can also occur in the time evolution, but also holds for the stationary flux equilibrium time evolution of driven open systems. For example, intriguing states. More precisely, the above situation implies an explicit scaling laws describing algebraic decoherence [43], anomalous breaking of detailed balance, since the “system” energy is not diffusion [70], or glass-like behavior [42, 71–73] have been iden- conserved. In an operational sense, this principle states that there tified in the long time asymptotics of driven spin systems close is a partition invariance for the temperature (or, more generally, to the stationary states. Conversely, the short time behavior of the noise level) present in the system: any arbitrary bipartition driven open lattice bosons shows universal scaling laws directly of the system can be chosen, one part can be traced out, and the witnessing the non-equilibrium drive, see Sec.V. This scaling resulting subsystem will be at the same temperature (noise level) can be related to a strongly pronounced non-equilibrium shape as the total system. This situation is the condition for a globally of the time-dependent distribution function in the early stages of well-defined temperature characteristic for systems in thermal evolution, and be traced back to conservation laws of the driven- equilibrium. More formally, thermal equilibrium can be detected dissipative generator of dynamics. by means of so-called fluctuation-dissipation relations (FDRs). The above discussion mainly focuses on the difference be- These connect the two fundamental observables in physical sys- tween equilibrium and non-equilibrium systems on the macro- tems — correlation and response functions (see Sec.IIB). In scopic level of observation. Another direction, still much less the case of thermal equilibrium, the connection is dictated by developed, concerns the distinction between classical and quan- the particle statistics alone. It is then given by the Bose- and tum effects. Again, although the quantum mechanical descrip- Fermi-distributions, respectively. Deviations from this univer- tion is necessary at a microscopic level, the persistence of quan- sal form, which has only two free parameters (temperature and tum effects at the macroscale is not guaranteed. This is mainly chemical potential, relating to the typical conserved quantities due to the Markovian noise level inherent to such quantum sys- energy and particle number), provide a necessary requirement tems. Nevertheless, systems with suitably engineered driven- for non-equilibrium conditions. dissipative dynamics show typical quantum mechanical phe- In non-equilibrium stationary states, no such general form ex- nomena such as phase coherence [50, 51, 57, 59, 60, 74], en- ists. We will encounter a concrete and simple example in the tanglement [54–56, 58], or topological order [52, 62, 66–68]. context of the driven Dicke model (a cavity photon coupled to Especially fermionic systems, which do not possess a classical a collective spin) in Sec. III A, where the form of the FDR de- limit, are promising in this direction. pends on the observable we are choosing (e.g., the position or momentum correlations and responses). On the other hand, thermal FDRs can emerge at long wave- B. Theoretical concepts and techniques length, even though the microscopic dynamics manifestly breaks detailed balance [26]. In particular, in three dimensions and The development of theoretical tools needed to perform the close to a critical point, a degeneracy of critical exponents in- transition from micro- to macro-physics in driven open quantum dicates a universal asymptotic thermalization, in the sense of an systems is still in progress, as a topic of current research. A rea- emergent thermal FDR. A similar phenomenology is observed son for the preliminary status of the theory lies in the fact that in a disordered multimode extension of the Dicke model, see two previously rather independent disciplines — quantum op- Sec. III B. This underlines the strongly attractive nature of the tics and many-body physics — need to be unified on a technical thermal equilibrium fixed point at low frequencies. Still in these level. systems, non-equilibrium conditions leave their traces in the dy- Quantum optical systems are well described microscopically namical response of the system, in terms of information that in terms of Markovian quantum master equations, which treat does not enter the FDR at leading order. For example, the criti- coherent Hamiltonian and driven-dissipative dynamics on equal 4 footing. To solve such equations both for their dynamics dation in terms of the renormalization group already mentioned and their stationary states, powerful techniques have been de- above, yet another tool developed and most clearly formulated vised. This comprises efficient exact numerical techniques for in a functional language. Finally, the functional integral based small enough systems, such as the quantum trajectories ap- on a single scalar quantity — the system’s action which encodes proach [45, 75, 76]. But it also includes analytical approaches all the dynamics on the microscopic scale — is a convenient such as perturbation theory for quantum master equations, e.g., framework when it comes to the classification of symmetries and in the frame of the Nakajima-Zwanzig projection operator tech- associated conservation laws, and their use in devising approxi- nique [77, 78], or mappings to P, W, or Q representations [79], mation schemes respecting them. casting the problem from a second quantized formulation into To put it short, while the driven open many-body systems partial differential equations. are well described by microscopic master equations, the tradi- A characteristic feature of traditional quantum optical sys- tional techniques of quantum optics cannot be used efficiently tems is the finite spacing of the few energy levels which play – at least not in the case where the generic complications of a role. When considering systems with a spatial continuum of many-body systems start to play a role. Conversely, their driven degrees of freedom instead, the energy levels become contin- open character makes it impossible to approach these problems uous. This does not mean that the microscopic modelling in in the framework of equilibrium many-body physics. This sit- terms of a quantum master equation is inappropriate: for the uation calls for a merger of the disciplines of quantum optics driven-dissipative terms in such an equation, the assumption of and many-body physics on a technical level. On the numerical spatially independent dissipative processes (such as atomic loss side, progress has been made in one spatial dimension recently or spontaneous emission) is still valid as long as the emitted by combining the method of quantum trajectories with power- wavelength of radiation is well below the spatial resolution at the ful density matrix renormalization group algorithms [83–86], scale where the microscopic model is defined. Indeed, in this sit- see [45] for an excellent review on the topic. For more analyti- uation, destructive interference of radiation justifies the descrip- cal approaches, the Keldysh functional or path integral [87–92] tion of driven dissipation in terms of incoherent processes. How- is ideally suited (but see [93] for a recent systematic perturbative ever, under these circumstances the smallness of a microscopic approach to lattice Lindblad equations with extensions to sophis- expansion parameter no longer guarantees the smallness of the ticated resummation schemes [94], and [95–97] for advanced associated perturbative correction. Here the reason is that in per- variational techniques). Conceptually, the latter captures the turbation theory, one is summing over intermediate states with most general situation in many-body physics — the dynamics propagation amplitudes down to the longest wavelengths. This of a density matrix under an arbitrary temporally local generator can lead to infrared divergences in naive perturbation theory. of dynamics. We refer to [98] for an introduction to the Keldysh This circumstance found its physical interpretation and techni- functional integral. In our context, it can be derived by a di- cal remedy in equilibrium in terms of the renormalization group rect functional integral quantization of the Markovian quantum [80, 81]. We emphasize that it is precisely this situation of long master equation. This procedure results in a simple translation wavelength dominance which underlies much of the universality, table from the master equation to the key player in the associated i.e., insensitivity to microscopic details, which is encountered Keldysh functional integral, the Markovian action. At first sight, when moving from the microscale to macroscopic observables the complexity of the non-equilibrium Keldysh functional inte- in many-particle systems. gral is increased by the characteristic “doubling” of degrees of The modern framework to understand many-particle prob- freedom compared to thermal equilibrium. However, it should lems in thermodynamic equilibrium is in terms of the functional be noted that it is precisely this feature which relates the Keldysh integral formulation of quantum field theory. The spectrum of functional integral much closer to real time (Heisenberg-von its application ranges over a remarkable range of energy scales, Neumann) evolution familiar from quantum mechanics. We will from ultracold atomic gases over condensed matter systems with demonstrate this in Sec.IIA, making it a rather intuitive object strong correlations to quantum chromodynamics and the quan- to work with. Furthermore, when properly harnessing symme- tum theory of gravity. It provides us with a well-developed tool- tries, the complexity of calculations can often be made compara- box of techniques, such as diagrammatic perturbation theory ble to thermodynamic equilibrium. Most importantly however, including sophisticated resummation schemes. But it also en- the powerful toolbox of quantum field theory is opened in this compasses non-perturbative approaches, which often capitalize way. It is thus possible to leverage the full power of sophisti- on the flexibility of the functional integral when it comes to pick- cated techniques from equilibrium field theory to driven open ing the relevant degrees of freedom for a given problem. This is many-body systems. the challenge of emergent phenomena, whose solution typically Relation to classical dynamical field theories — The quan- is strongly scale dependent [82]. Familiar examples include tum mechanical Keldysh formulation reduces to the so-called an efficient description of emergent Cooper pair or molecular Martin-Siggia-Rose (MSR) functional integral in the (semi-) degrees of freedom in interacting fermion systems, or vortices classical limit, in turn equivalent to a stochastic Langevin equa- which conveniently parameterize the long-wavelength physics tion formulation [98]. This statement will be made precise and of interacting bosons in two dimensions. The description of the discussed in Sec.IIC. A large amount of work has been dedi- change of physics with scale was given its mathematical foun- cated to this limit in the past. On the one hand, this concerns 5 dynamical aspects of equilibrium statistical mechanics, and we bosonic degrees of freedom. However, it is also possible to ad- refer to the classic work by Hohenberg and Halperin [1] for an dress spin systems in terms of functional integrals [106], and overview. This work shows that, while the static universal crit- simple models systems have been analyzed in this way, see [107] ical behavior is determined by the symmetries and the dimen- and Sec. III A. More sophisticated approaches to spin systems sionality of the problem, the dynamical critical behavior is sen- were elaborated in the context of multimode optical cavities sitive to additional dynamical conservation laws. This leads to in [108], and systematically for various symmetries for lattice a fine structure, defining dynamical universality classes spanned systems in [109]. Fermi statistics is also conveniently imple- by models A–J. These models also provide a convenient frame- mented in the functional integral formulation. This is relevant, work to describe the statistics of work, summarized in Jarzyn- e.g., for driven open Fermi gases in optical cavities [21, 22, 110] ski’s work theorem [99, 100] and Crooks’ relation [101]. On or lattices [111], or dissipatively stabilized topological fermion the other hand, non-equilibrium situations are captured as well. matter [66–68]. Here, we highlight in particular genuine non-equilibrium univer- sality classes, which are not smoothly connected to the equilib- rium models. Among them is the problem of reaction-diffusion C. Experimental platforms models [102] including directed percolation [4], which is rele- vant to certain chemical processes (for an implementation of this The progress in controlling, manipulating, detecting and scal- universality class with driven Rydberg gases, see [47]). Another ing up driven open quantum systems to many-body scenarios has key example is surface growth, described by the Kardar-Parisi- been impressive over the last decade. Here we sketch the basic Zhang equation [69], giving rise to a non-equilibrium universal- physics of three representative platforms, and indicate the rele- ity class which is at the heart of driven phenomena such as the vant microscopic theoretical models in the frame of the Marko- growth of bacterial colonies or the spreading of fire. vian quantum master equation. In later sections, we will trans- In the same class of approaches ranges the so-called Doi-Peliti late this physics into the language of the Keldysh functional inte- functional integral [103, 104], which is a functional representa- gral. For each of these platforms, excellent reviews exist, which tion of classical master equations, and may be viewed as an MSR we refer to at the end of Sec.ID together with further litera- theory with a specific, highly non-linear appearance of the field ture on open systems. The purpose of this section is to give variables. It reduces to the conventional MSR form in a leading an overview only, and to put the respective platforms into their order Taylor expansion of the field non-linearities. A compre- overarching context as driven open quantum systems with many hensive overview of models, methods, and physical phenomena degrees of freedom. in the (semi-)classical limit is provided in [5]. We also note that the usual mean field theory, where corre- lation functions are factorized into products of field amplitudes 1. Cold atoms in an optical cavity, and microcavity arrays: and which is often used as an approximation to the quantum driven-dissipative spin-boson models master equation in the literature, corresponds to a further for- mal simplification of the semi-classical limit. Here the effects of Cavity quantum electrodynamics (cavity QED), with its fo- noise are neglected completely. Conversely, the semi-classical cus on strong light-matter interactions, is a growing field of re- limit represents a systematic extension of mean field theory, search, which has experienced several groundbreaking advances which includes the Markovian noise fluctuations. This level of in the past few years. Historically, these systems were devel- approximation is referred to as optimal path approximation in oped as few or single atom experiments, detecting the radiation the literature on MSR functional integrals [5, 98]. properties of atoms, which are strongly coupled to a quantized In many cases, even though the microscopic description is light field. The focus has recently been shifted towards loading in terms of a quantum master equation, at long distances the more and more atoms inside a cavity. Thereby, not only single Keldysh field theory reduces to a semi-classical MSR field the- particle dynamics in strong radiation fields can be probed, but ory. The reason is the finite Markovian noise level that such also collective, macroscopic phenomena, which are driven by systems exhibit generically, as explained in Sec.IIC. The prefix light-matter interactions. For an excellent review on this topic, “semi” refers to the fact that phase coherence may still persist covering important experimental and theoretical developments, in such circumstances — the situation is comparable to a Bose- see [10]. In the experiments, cold atoms are loaded inside an Einstein condensate at finite temperature. Recently however, sit- optical or microwave cavity, for which the coherent interaction uations have been identified where the drastic simplifications of between the atomic internal states and a single cavity mode dom- the semi-classical limit do not apply. In particular, this occurs inates over dissipative processes [112]. The atoms absorb and in systems with dark states — pure quantum states which are emit cavity photons, thereby changing their internal states. Due dynamical fixed points of driven-dissipative evolution [50, 105]. to this process, the spatial modulation of the intra-cavity light In these cases, classical dynamical field theories are inappropri- field induces a coupling of the cavity photons to the atomic in- ate, which calls for the development of quantum dynamical field ternal state as well as their motional degree of freedom [113]. In theories [27]. These developments are just in their beginnings. this way, cavity photons mediate an effective atom-atom interac- In this review, we concentrate on systems composed of tion, which leads to a back-action of a single atom on the motion 6 of other atoms inside the cavity. This cavity mediated interaction of the pump laser was chosen such that the atoms effectively re- is long-ranged in space and represents the source of collective ef- main in the internal ground state, but acquire a characteristic re- fects in cold atomic clouds within cavity QED experiments. One coil momentum when scattering with a cavity photon. This scat- hallmark of collective dynamics in cavity QED has been the ob- tering creates a collective, motionally excited state, which re- servation of self-organization of a Bose-Einstein condensate in places the role of an individual, internally excited atom. The ex- an optical cavity. This is accompanied by a a Dicke phase tran- perimental realization of a superradiance transition in the Dicke sition via breaking of a discrete Z2-symmetry of the underlying model is usually inhibited, since the required coupling strength model [9, 114–116]. by far exceeds the available value of the atomic dipole coupling. Although the dominant dynamics in these systems is coherent, However, for the BEC in the cavity, the energy scales of the ex- there are dissipative effects which cannot be discarded. These cited modes are much lower than the optical scale of the atomic are the loss of cavity photons due to imperfections in the cav- modes. In this way, the superradiance transition indeed became ity mirrors, or spontaneous emission processes of atoms, which experimentally accessible. This was inspired by a theoretical emit photons into transverse modes. These effects modify the proposal using two balanced Raman channels between different dynamics of the system on the longest time scales. Therefore, internal atomic states inside an optical cavity, which reduced the they become relevant for macroscopic phenomena, such as phase effective level splitting of the internal states to much lower en- transitions and collective dynamics. For instance, dissipative ef- ergy scales [124]. fects have been shown theoretically [117, 118] and experimen- In addition to the unitary dynamics represented by the Dicke tally [119] to modify the critical exponent of the Dicke transi- model, the cavity is subject to permanent photon loss due to im- tion compared to its zero temperature value. This illustrates that perfections in the cavity mirrors. For high finesse cavities, the for the analysis of collective phenomena in cavity QED experi- coupling of the intra cavity photons to the surrounding vacuum ments, the dissipative nature of the system has to be taken into radiation field is very weak, and the latter can be eliminated in a account properly [107, 120]. Born-Markov approximation [77]. This results in a Markovian Typically, for a cavity field which has a very narrow spectrum, quantum master equation for the system’s density matrix the atomic internal degrees of freedom can be reduced to two − L internal states, whose transitions are near resonant to the pho- ∂tρ = i[HD, ρ] + ρ, (2) ton frequency. The operators acting on these two internal states where ρ is the density matrix for the intra cavity system and L can be represented by Pauli matrices, making them equivalent adds dissipative dynamics to the coherent evolution of the Dicke to a spin-1/2 degree of freedom. A very important model in the model. For a vacuum radiation field, it is given by framework of cavity QED is the Dicke model [121–123], which † 1 † describes N atoms (i.e., two-level systems) coupled to a single Lρ = κ LρL − 2 {L L, ρ} . (3) quantized photon mode. This is expressed by the Hamiltonian (here and in the following we set ~ = 1) The Lindblad operator L acting on the density matrix describes pure photon loss (L = a) with an effective loss rate κ; the lat- N N † ωz X z g X x † ter depends on system specific parameters, but is typically the HD = ω0a a + σ + √ σi a + a . (1) 2 i N smallest scale in the master equation (2)[10]. i=1 i=1 In generic cavity experiments, there are also atomic sponta- z Here, ω0 is the photon frequency, σ = |1ih1| − |0ih0| repre- neous emission processes. The atoms scatter a laser or cavity sents the splitting of the two atomic levels with energy differ- photon out of the cavity, and this represents a source of deco- x ence ωz. σ = |0ih1| + |1ih0| describes the coherent excitation herence. This process has been considered in Ref. [107], and and de-excitation of the atomic state proportional to the atom- leads to an effective decay rate of the atomic excited state and photon coupling strength g. The Dicke model features a dis- therefore to a dephasing of the atoms, described by additional crete Z Ising symmetry: it is invariant under the transforma- z 2 Lindblad operators Li = σi [125]. However, these losses are † x † x tion (a , a, σi ) 7→ (−a , −a, −σi ). In the thermodynamic limit, at least three orders of magnitude smaller than the cavity de- for N → ∞, it features a phase transition, which spontaneously cay rate and typically not considered [9, 120]. The basic model breaks the Ising symmetry. Crossing the transition, the system for cavity QED with cold atoms is therefore represented by the enters a superradiant phase, characterized by finite expectation Dicke model with dissipation, formulated in terms of the master x values hai , 0, hσi i , 0. This describes condensation of the equation (2). Its dynamics is discussed in Sec. III, including the cavity photons, i.e., the formation of a macroscopically occu- Dicke superradiance transition. pied, coherent intra-cavity field, and a “ferromagnetic” ordering The Dicke model Hamiltonian takes the form HD = Hc + Hs + of the atoms in the x-direction. Hcs, where Hc,s,cs represent the bosonic cavity, spin, and spin- Although the Dicke model is a standard model√ for cavity QED boson sectors, respectively. There are many directions to go be- experiments in the ultra-strong coupling limit Ng > ωzω0, it yond the Dicke model with single collective spin, still keeping has been realized only very recently in cold atom experiments, the basic feature of coupling spin to boson (cavity photon) de- where an entire BEC was placed inside an optical cavity [9]. It grees of freedom. One direction – relevant to future cold atom has been shown that this setup maps to a Dicke model, with a experiments – is to consider multimode cavities instead of a sin- “collective” spin degree of freedom [115]. Here, the detuning gle one. In particular, in conjunction with quenched disorder, 7 intriguing analogies to the physics of quantum glasses can be the excited spin states, the latter degrees of freedom can be established in this way [108, 126]. Here, the global coupling of integrated out. In this limit, their physical effect is to gen- all spins to a single mode gi ≡ g is replaced by random couplings erate Kerr-type bosonic non-linearities, giving rise to driven gi,`, where the index ` now refers to a collection of cavity modes. open variants of the celebrated Bose-Hubbard model. These The many-body physics of such an open system is discussed in models can even be brought into the correlation dominated Sec. III. regime [24, 127, 128, 133]. Oftentimes, these approximations on The basic building block of the Dicke model is the spin-boson the spin sector actually apply, and it is both useful and interesting † x term of the form Hsc ∼ (a + a )σ , i.e. a Rabi type non- to study these effective low frequency bosonic theories instead linearity that preserves the Z2 symmetry of Hc,s. In the con- of the full many-body spin-boson problems, see Refs. [134–138] text of circuit quantum electrodynamics, a natural many-body for recent work in this direction. generalization of Hamiltonians with a spin-boson interaction is to consider entire arrays of microcavities (instead of consider- ing many modes within a single cavity). These cavities can be 2. Exciton-polariton systems: driven open interacting bosons coupled to each other by single photon tunnelling processes be- tween adjacent cavities, giving rise to Hubbard-type hopping † Exciton-polaritons are an extremely versatile experimental terms ∼ Jai a j + h. c., where i, j now label the spatial index platform, which is documented by the richness of physical phe- of the cavities. In cirquit QED, strong non-linearities can be nomena that have been studied in these systems both in the- generated, e.g., by coupling to adjacent qubits made of Cooper ory and experiment. For a comprehensive account of the sub- pair boxes [127, 128]. This gives rise many-body variants of ject, we refer to a number of excellent review articles [13, 139, the Rabi model [129], whose phase diagrams have been studied 140]. A Keldysh functional integral approach is discussed in recently [130]. Furthermore, for the implementation of lattice Refs. [141, 142], which provides both a microscopic deriva- Dicke models with large collective spins, the use of hybrid quan- tion of an exciton-polariton model and a mean field analysis tum systems consisting of superconducting cavity arrays cou- including Gaussian fluctuations. At this point, we content our- pled to solid-state spin ensembles have been proposed [131]. selves with a short introduction, with the aim of showing that In many physical situations (away from the ultra-strong cou- in a suitable parameter regime, exciton-polaritons very natu- pling limit), a form of the spin-cavity interaction alternative rally provide a test-bed to study Bose condensation phenom- to the Rabi term is more appropriate, with non-linear building ena out of thermal equilibrium. Similar physics can also arise + † − block Hcs ∼ aσ + a σ . In fact, this form results naturally in a variety of other systems, including condensates of pho- from the weak coupling rotating wave approximation of a driven tons [143], magnons [144], and potentially excitons [145]. Re- spin-cavity problem. For a single cavity mode and spin, the re- markably, even cold atoms could be brought to condense in a sulting model is the Jaynes-Cummings model, which in contrast non-equilibrium regime, where continuous loading of atoms bal- to the Dicke model possesses a continuous U(1) phase rotation ances three-body losses [146], or in atom laser setups [147–149]. − − symmetry under a → eiθa, σ → eiθσ . A basic experimental setup for exciton-polaritons consists of Clearly, when such systems are driven coherently via a Hamil- a planar semiconductor microcavity embedding a quantum well † tonian Hd = Ω(a + a ) or suitable multimode generalizations (see Fig.1 (a)). This setting allows for a strong coupling of cav- thereof, both the U(1) and even the Z2 symmetries of the above ity light and matter in the quantum well, as originally proposed models are broken explicitly. Coherent drive is usually the sim- in [150]. The free dynamics of the elementary excitations of this plest way to compensate for unavoidable losses due to cavity system — i.e., of cavity photons and Wannier-Mott excitons — leakage, although incoherent pumping schemes are conceivable is described by the quadratic Hamiltonian [13] using multiple qubits [132]. An advantage of such schemes is that the symmetries of the underlying dynamics are preserved H0 = HC + HX + HC−X, (4) or less severely corrupted in this way (cf. the discussion in Sec.IID). In other platforms, such as exciton-polariton systems where the parts of the Hamiltonian involving only photons and (cf. the subsequent section), incoherent pumping is more natural excitons, respectively, take the same form, which is given by from the outset. (here the index α labels cavity photons, α = C, and excitons, All the systems discussed here represent genuine instances of α = X, respectively)1 driven open many-body systems. Besides the coherent drive, Z dq X they undergo dissipative processes, which have to be taken into H ω q a† q a q . α = 2 α( ) α,σ( ) α,σ( ) (5) account for a proper understanding of their time evolution and (2π) σ stationary states. A generic feature of these processes is their locality. For the effective spin degrees of freedom, the typical − processes are qubit decay (local Lindblad operators Li = σi ) and z dephasing (Li = σi ). For the bosonic component, local single- 1 In Ref. [141, 142], a different model for excitons is used: they are assumed photon loss is dominant, Li = ai. to be localized by disorder, and interactions are included by imposing a hard- Under various circumstances, such as a low population of core constraint. 8
† Field operators aα,σ(q) and aα,σ(q) create and destroy, respec- (a) tively, a photon or exciton (note that both are bosonic excita- photon tions) with in-plane momentum q and polarization σ (there are exciton two polarization states of the exciton which are coupled to the e cavity mode [13]). For simplicity, we neglect polarization ef- fects leading to an effective spin-orbit coupling [13]. Due to the h confinement in the transverse (z) direction, i.e., along the cav- ity axis, the motion of photons in this direction is quantized as q = πn/l , where n is a positive integer, and l is the length z,n z z Braggmirror QW Braggmirror of the cavity. In writing the Hamiltonian (5), we are assuming that only the lowest transverse mode is populated, which leads to (b) upper ! a quadratic dispersion as a function of the in-plane momentum polariton q q = |q| = q2 + q2: x y photon q q2 excitation ω (q) = c q2 + q2 = ω0 + + O(q4). (6) relaxation C z,1 C 2m C exciton 0 Here, c is the speed of sound, ω = cqz, , and the effective mass C 1 lower of the photon is given by mC = qz,1/c. Typically, the value of the photon mass is orders of magnitude smaller than the mass of the polariton exciton, so that the dispersion of the latter appears to be flat on emission the scale of Fig.1 (b). q Upon absorption of a photon by the semiconductor, an exciton is generated. This process (and the reverse process of the emis- Figure 1. (a) Schematic of two Bragg mirrors forming a microcavity, sion of a photon upon radiative decay of an exciton) is described in which a quantum well (QW) is embedded. In the regime of strong by light-matter interaction, the cavity photon and the exciton hybridize and form new eigenmodes, which are called exciton-polaritons. (b) Energy Z dq X H a† q a q . . , dispersion of the upper and lower polariton branches as a function of in- C−X = ΩR 2 X,σ( ) C,σ( ) + H c (7) (2π) σ plane momentum q. In the experimental scheme illustrated in this figure (cf. Ref. [12]), the incident laser is tuned to highly excited states of the where ΩR is the rate of the coherent interconversion of photons quantum well. These undergo relaxation via emission of phonons and into excitons and vice versa. The quadratic Hamiltonian (4) can scattering from polaritons, and accumulate at the bottom of the lower be diagonalized by introducing new modes — the lower and polariton branch. In the course of the relaxation process, coherence is upper exciton-polaritons, ψLP,σ(q) and ψUP,σ(q) respectively, quickly lost, and the effective pumping of lower polaritons is incoher- which are linear combinations of photon and exciton modes. The ent. dispersion of lower and upper polaritons is depicted in Fig.1 (b). In the regime of strong light-matter coupling, which is reached when ΩR is larger than both the rate at which pho- peculiarity of exciton-polaritons is that they are excitations with tons are lost from the cavity due to mirror imperfections and the relatively short lifetime. In turn, this necessitates continuous re- non-radiative decay rate of excitons, it is appropriate to think of plenishment of energy in the form of laser driving in order to exciton-polaritons as the elementary excitations of the system. maintain a steady population. In Fig.1 (b), we consider the In experiments, it is often sufficient to consider only lower case in which the excitation laser is tuned to energies well above polaritons in a specific spin state, and to approximate the disper- the lower polariton band. The thus created high-energy excita- sion as parabolic [13]. Interactions between exciton-polaritons tions are deprived of their excess energy via phonon-polariton originate from various physical mechanisms, with a dominant and stimulated polariton-polariton scattering. Eventually, they contribution stemming from the screened Coulomb interactions accumulate at the bottom of the lower polariton band. As a con- between electrons and holes forming the excitons. Again, in the sequence of multiple scattering processes, the coherence of the low-energy scattering regime, this leads to an effective contact incident laser field is quickly lost, and the effective pumping of interaction between lower polaritons. As a result, the low-energy lower polaritons is incoherent. description of lower polaritons takes the form (in the following A phenomenological model for the dynamics of the lower po- we drop the subscript indices in ψLP,σ)[13] lariton field, which accounts for both the coherent dynamics gen- Z " 2 ! # erated by the Hamiltonian (8) and the driven-dissipative one de- † 0 ∇ † 2 2 HLP = dx ψ (x) ω − ψ(x) + ucψ (x) ψ(x) . (8) scribed above, was introduced in Ref. [151]. It involves a dissi- LP 2m LP pative Gross-Pitaevskii equation for the lower polariton field that While this Hamiltonian is quite generic and arises also, e.g., in is coupled to a rate equation for the reservoir of high-energy ex- cold bosonic atoms in the absence of an external potential, the citations. However, for the study of universal long-wavelength 9 behavior in Sec.IV, this degree of microscopic modeling is ac- reached, for which the dissipative coupling to the environment tually not required: indeed, any (possibly simpler) model that becomes visible in experimental observables. Such dissipation possesses the relevant symmetries (see the discussion at the be- may even hinder the system from relaxation towards a well de- ginning of Sec.IV) will yield the same universal physics. Such fined steady state. a model can be obtained by describing incoherent pumping and A relevant example of a dissipative coupling is decoherence losses of lower polaritons by means of a Markovian master equa- of an atomic cloud induced by spontaneous emission of atoms tion:2 in the lattice [158, 159]. In this way, the many-body system is heated up, and therefore driven away from the low entropy ∂tρ = −i[HLP, ρ] + Lρ, (9) state in which it was prepared initially. A detailed discussion of the microscopic physics and its long time dynamics can be where Lρ encodes incoherent single-particle pumping and found in [42–44], see also the review article [45]. For bosonic losses, as well as two-body losses: atoms in optical lattices, the coherent dynamics is described by Z the Bose-Hubbard Hamiltonian [160, 161] † 2 Lρ = dx γpD[ψ(x) ]ρ + γlD[ψ(x)]ρ + 2udD[ψ(x) ]ρ , X U X H = −J b†b + n (n − 1), (12) (10) BH l m 2 l l hl,mi l where which models bosonic atoms in terms of the creation and annihi- † † 1 † D[L]ρ = LρL − {L L, ρ} (11) lation operators [bl , bm] = δlm in the lowest band of a lattice with 2 site indices l, m. The atoms hop between neighboring lattice sites with an amplitude J, and experience an on-site repulsion U. The reflects the Lindblad form, and γp, γl, 2ud are the rates of single- particle pumping, single-particle loss, and two-body loss, re- lattice potential V(x), which leads to the second quantized form spectively. The inclusion of the non-linear loss term ensures of the Bose-Hubbard model, is created by the superposition of saturation of the pumping. An analogous mechanism is im- counter-propagating laser beams in each spatial dimension. The plemented in the above-mentioned Gross-Pitaevskii description. coherent laser field couples two internal atomic states via stimu- More precisely, in the spirit of universality, the above quantum lated absorption and emission, which leads to the single particle master equation (9) and the above mentioned phenomenologi- Hamiltonian pˆ 2 ∆ Ω(xˆ) cal model reduce to precisely the same low frequency model for H = − σz − σx , (13) bosonic degrees of freedom upon taking the semiclassical limit atom 2m 2 2 in the Keldysh path integral associated to Eq. (9) (see Sec.IIC where pˆ is the atomic momentum operator, ∆ is the detuning of for its implementation), and integrating out the upper polariton the laser from the atomic transition frequency, and Ω(xˆ) is the reservoir in the phenomenological model. laser field at the atomic position. For large detuning, the excited state of the atom can be traced out, which leads to the lattice Hamiltonian 3. Cold atoms in optical lattices: heating dynamics pˆ 2 |Ω(xˆ)|2 Heff = + . (14) atom 2m 4∆ In recent years, experiments with cold atoms in optical lat- tices have shown remarkable progress in the simulation of many- This describes a lattice potential for the ground state atoms body model systems both in and out of equilibrium. A par- generated by the spatially modulated AC Stark shift. Adding ticular strength of cold atom experiments is the unprecedented an atomic interaction potential and expanding both the single tuneability of model parameters, such as the local interaction particle Hamiltonian (14) and the interaction in terms of Wan- strength and the lattice hopping amplitude. This becomes possi- nier states, the leading order Hamiltonian is the Bose-Hubbard ble by, e.g., manipulation of the lattice laser and external mag- model (12)[161]. netic fields. It comes along with a very weak coupling of the In the semi-classical treatment of the atom-laser interaction, system to the environment, such that the dynamics can often be spontaneous emission events are neglected, as their probabil- seen as isolated on relevant time scales for typical measurements ity is typically very small (see below for a more precise state- of static, equilibrium correlations. However, more and more ex- ment). They can be taken into account on the basis of optical periments start to investigate the realm of non-equilibrium phe- Bloch equations [158], which leads, after elimination of the ex- nomena with cold atoms, e.g., by letting systems prepared in a cited atomic state, to an additional, driven-dissipative term in non-equilibrium initial condition relax in time towards a steady the atomic dynamics. It describes the decoherence of the atomic state [34, 36, 155–157]. With these experiments, time scales are state due to spontaneous emission, i.e., position dependent ran- dom light scattering. The leading order contribution to the dy- namics in the basis of lowest band Wannier states is captured by the master equation 2 X While this approach captures the universal behavior, we note that non- ∂ ρ = −i[H , ρ] − γ [n , [n , ρ]], (15) Markovian effects can be of key importance for other properties [152–154]. t BH l l l 10
Sec. II B 2: there we discuss the mean field theory of conden- sation in a bosonic many-body system with particle losses and pumping. The semi-classical limit of this model and its valid- ity are the content of Sec.IIC. This is followed by a discussion of symmetries and conservation laws in the Keldysh formalism in Sec.IID. In particular, we point out a symmetry that allows one to distinguish equilibrium from non-equilibrium conditions. Finally, an advanced field theoretical tool — the open system Figure 2. Illustration of the coherent (a) and the incoherent (b) con- functional renormalization group — is introduced in Sec.IIE. tribution to the dynamics stemming from the coupling of the atoms to a laser with amplitude |Ω| and large detuning ∆. (a) Via the AC Stark PartII harnesses this formalism to generate an understanding effect, stimulated absorption and emission lead to a coherent periodic of the physics in different experimental platforms. We begin in |Ω|2 Sec. III with simple but paradigmatic spin models with discrete potential with amplitude 4∆ . (b) Stimulated absorption and subsequent spontaneous emission lead to effective decoherence of the atomic state Ising Z2 symmetry in driven flux equilibrium states. In partic- 2 with rate Γ |Ω| , where Γ is the microscopic spontaneous emission rate. ular, we discuss the physics of the driven open Dicke model in 4∆ Sec. III A. This is followed by an extended variant of the lat- ter in the presence of disorder and a multimode cavity, which hosts an interesting spin and photon glass phase in Sec. III B. where n = b†b is the local atomic density, and ρ is the many- l l l Sec.IV is devoted to the non-equilibrium stationary states of |Ω|2 body density matrix. For a red detuned laser the rate γ = Γ 4∆2 is bosons with a characteristic U(1) phase rotation symmetry: the proportional to the microscopic spontaneous emission rate Γ and driven-dissipative condensates introduced in Sec. II B 2. After the laser amplitude |Ω|. Note the suppression of the scale γ by some additional technical developments relating to U(1) sym- a factor Γ/∆ 1 for large detuning, compared to the strength metry in Sec.IVA, we discuss critical behavior at the Bose eff of Hatom. The coherent and incoherent contribution of the atom- condensation transition in three dimensions. In particular, we laser coupling to the dynamics is illustrated in Fig.2. show the decrease of a parameter quantifying non-equilibrium The dissipative term in the master equation (15) leads to an strength in this case (Sec.IVB). The opposite behavior is ob- energy increase hHBHi(t) ∼ t linear in time, and thus to heat- served in two (Sec.IVC) and one (Sec.IVD) dimensions. Fi- ing. Furthermore, it introduces decoherence in the number state nally, leaving the realm of stationary states, in Sec.V we dis- basis, i.e., it projects the local density matrix on its diagonal in cuss an application of the Keldysh formalism to the time evolu- Fock space and leads to a decrease of the coherences in time [42– tion of open bosonic systems in one dimension, which undergo 45]. Starting from a low entropy state at t = 0, the heating leads number conserving heating processes. We set up the model in to a crossover from coherence dominated dynamics at short and Sec.VA, derive the kinetic equation for the distribution func- intermediate times [40] to a decoherence dominated dynamics at tion in Sec.VB, and discuss the relevant approximations and long times [42–44]. In one dimension, both regimes have been physical results in Secs.VC andVD, respectively. Conclusions analyzed extensively both numerically (with a focus on deco- are drawn in Sec.VI. Finally, brief introductions to functional herence dominated dynamics) as well as analytically, and dis- differentiation and Gaussian functional integration are given in play several aspects of non-equilibrium universality, see Sec.V. AppendicesA andB. Therefore, heating in interacting lattice systems represents a cru- Reflecting the bipartition of this review, the scope of it is cial example for universality in out-of-equilibrium dynamics, twofold. On the one hand, it develops the Keldysh functional in- which can be probed by cold atom experiments. tegral approach to driven open quantum systems “from scratch”, in a systematic and coherent way. It starts from the Marko- vian quantum master equation representation of driven dissipa- D. Outline and scope of this review tive quantum dynamics [77, 78, 162], and introduces an equiv- alent Keldysh functional integral representation. Direct contact The remainder of this review is split into two parts. is made to the language and typical observables of quantum op- PartI develops the theoretical framework for the e fficient de- tics. It does not require prior knowledge of quantum field the- scription of driven open many-body quantum systems. In Sec.II, ory, and we hope that it will find the interest of — and be use- we begin with a direct derivation of the open system Keldysh ful for — researchers working on quantum optical systems with functional integral from the many-body quantum master equa- many degrees of freedom. On the other hand, this review doc- tion (Sec.IIA). We then discuss in Sec. II B 1 in detail a simple uments some recent theoretical progresses made in this concep- example: the damped and driven optical cavity. This allows the tual framework in a more pedagogical way than the original lit- reader to familiarize with the functional formalism. In particu- erature. We believe that this not only exposes some interesting lar, the key players in terms of observables — correlation and re- physics, but also demonstrates the power and flexibility of the sponse functions — are described. We also point out a number of Keldysh approach to open quantum systems. exact structural properties of Keldysh field theories, which hold This work is complemetary to excellent reviews putting more beyond the specific example. Another example is introduced in emphasis on the specific experimental platforms partially men- 11 tioned above: The physics of driven Bose-Einstein condensates the norm tr(ρ) of the system density matrix, this term must be ac- in optical cavities is reviewed in [10]. A general overview of companied by fluctuations. The corresponding term, where the driven ultracold atomic systems, specifically in optical lattices, Lindblad operators act from both sides onto the density matrix, is provided in [45], and systems with engineered dissipation are is referred to as recycling or quantum jump term. Dissipation described in [53]. Detailed accounts for exciton-polariton sys- occurs at rates γα which are non-negative, so that the density tems are given in [13, 139, 140]; specifically, we refer to the matrix evolution is completely positive, i.e., the eigenvalues of review [141] working in the Keldysh formalism. The physics of ρ remain positive under the combined dynamics generated by H microcavity arrays is discussed in [14, 24, 129, 133], and trapped and D [168]. If the index α is the site index in an optical lattice ions are treated in [17, 163]. We also refer to recent reviews or in a microcavity array, or even a continuous position label (in on additional upcoming platforms of driven open quantum sys- which case the sum is replaced by an integral), in a translation tems, such as Rydberg atoms [164] and opto-nanomechanical invariant situation there is just a single scale γα = γ for all α settings [165]. For a recent exposition of the physics of quan- associated to the dissipator. tum master equations and to efficient numerical techniques for The quantum master equation (16) provides an accurate de- their solution, see [45]. scription of a system-bath setting with a strong separation of scales. This is generically the case in quantum optical systems, which are strongly driven by external classical fields. More pre- Part I cisely, there must be a large energy scale in the bath (as com- pared to the system-bath coupling), which justifies to integrate Theoretical background out the bath in second-order time-dependent perturbation the- ory. If in addition the bath has a broad bandwidth, the combined II. KELDYSH FUNCTIONAL INTEGRAL FOR DRIVEN Born-Markov and rotating-wave approximations are appropri- OPEN SYSTEMS ate, resulting in Eq. (16). A concrete example for such a set- ting is provided by a laser-driven atom undergoing spontaneous emission. Generic condensed matter systems do not display such In this part, we will be mainly concerned with a Keldysh field a scale separation, and a description in terms of a master equa- theoretical reformulation of the stationary state of Markovian tion is not justified. However, the systems discussed in the intro- many-body quantum master equations. As we also demonstrate, duction, belong to the class of systems which permit a descrip- this opens up the powerful toolbox of modern quantum field the- tion by Eq. (16). We refer to them as driven open many-body ory for the understanding of such systems. quantum systems. The quantum master equation, examples of which we have Due to the external drive these systems are out of thermody- already encountered in Sec.IC, describes the time evolution of namic equilibrium. This statement will be made more precise a reduced system density matrix ρ and reads [77, 78] in Sec. II D 1 in terms of the absence of a dynamical symme- ! try which characterizes any system evolving in thermodynamic X † 1 † ∂tρ = Lρ = −i[H, ρ] + γα LαρL − {L Lα, ρ} , (16) equilibrium, and which is manifestly violated in dynamics de- α 2 α α scribed by Eq. (16). Its absence reflects the lack of energy con- servation and the epxlicit breaking of detailed balance. As stated where the operator L acts on the density matrix ρ “from both in the introduction, the main goal of this review is point out the sides” and is often referred to as Liouville superoperator or Li- macroscopic, observable consequences of this microscopic vio- ouvillian (sometimes this term is reserved for the second con- lation of equilibrium conditions. tribution on the RHS of Eq. (16) alone). There are two contri- butions to the Liouvillian: first, the commutator term, which is familiar from the von Neumann equation, describes the coher- ent dynamics generated by a system Hamiltonian H; the second A. From the quantum master equation to the Keldysh functional part, which we will refer to as the dissipator D 3, describes the integral dissipative dynamics resulting from the interaction of the sys- tem with an environment, or “bath.” It is defined in terms of a In this section, starting from a many-body quantum master set of so-called Lindblad operators (or quantum jump operators) equation Eq. (16) in the operator language of second quantiza- Lα, which model the coupling to that bath. The dissipator has a tion, we derive an equivalent Keldysh functional integral. We characteristic Lindblad form [166, 167]: it contains an anticom- focus on stationary states, and we discuss how to extract dynam- mutator term which describes dissipation; in order to conserve ics from this framework in Sec.V. Our derivation of the Keldysh functional integral applies to a theory of bosonic degrees of free- dom. If spin systems are to be considered, it is useful to first per- form the typical approximations mapping them to bosonic fields, 3 We use the term “dissipation” here for all kinds of environmental influences and then proceed along the construction below (see Sec. III and on the system which can be captured in Lindblad form, including effects of Refs. [107, 109]). Clearly, this amounts to an approximate treat- decay and of dephasing/decoherence. ment of the spin degrees of freedom; for an exact (equilibrium) 12 functional integral representation for spin dynamics, taking into a) account the full non-linear structure of their commutation rela- (t0) tions, we refer to Ref. [169]. For fermionic problems, the con- t t t0 | i struction is analogous to the bosonic case presented here. How- U ever, a few signs have to be adjusted to account for the fermionic b) anticommutation relations [98, 169]. The basic idea of the Keldysh functional integral can be de- ⇢(t0) t t veloped in simple terms by considering the Schrodinger¨ vs. the U U † t von Neumann equation, c) + contour
i∂t|ψ(t)i = H|ψ(t)i ⇒ |ψ(t)i = U(t, t0)|ψ(t0)i, ⇢(tf ) ⇢(t0)
† tf =+ - contour t = ∂tρ(t) = −i[H, ρ(t)] ⇒ ρ(t) = U(t, t0)ρ(t0)U (t, t0), 1 0 1 (17) Figure 3. Idea of the Keldysh functional integral. a) According to the −iH(t−t0) where U(t, t0) = e is the unitary time evolution opera- Schrodinger¨ equation, the time evolution of a pure state vector is de- −iH(t−t0) tor. In the first case, a real time path integral can be constructed scribed the unitary operator U(t, t0) = e . In the Feynman func- along the lines of Feynman’s original path integral formulation tional integral construction, the time evolution is chopped into infinites- of quantum mechanics [170]. To this end, a Trotter decomposi- imal steps of length δt, and completeness relations in terms of coherent tion of the evolution operator states are inserted in between consecutive time steps. This insertion is signalled by the red arrows. b) In contrast, if the state is mixed, a den- −iH(t−t0) N e = lim (1 − iδtH) , (18) sity matrix must be evolved, and thus two time branches are needed. As N→∞ explained in the text, the dynamics need not necessarily be restricted to
t−t0 unitary evolution. The most general time-local (Markovian) dynamics with δt = N , is performed. Subsequently, in between the is generated by a Liouville operator in Lindblad form. c) For the anal- factors of the Trotter decomposition, completeness relations in ysis of the stationary state, we are interested in the real time analog of terms of coherent states are inserted in order to make the (nor- a partition function Z = tr ρ(t), starting from t0 = −∞ and running until mal ordered) Hamilton operator a functional of classical field t f = +∞. The trace operation connects the two time branches, giving variables. This is illustrated in Fig.3 a), and we will perform rise to the closed Keldysh contour. these steps explicitly and in more detail below in the context of open many-body systems. Crucially, we only need one set of field variables representing coherent Hamiltonian dynamics, So far, we have concentrated on closed systems which evolve which corresponds to the forward evolution of the Schrodinger¨ according to purely Hamiltonian dynamics. However, we can state vector. It is also clear that — noting the formal analogy of allow for a more general generator of dynamics and still proceed the operators e−iH(t−t0) and e−βH — this construction can be lever- along the two-branch strategy. The most general (time local) aged over to the case of thermal equilibrium, where the “Trotter- evolution of a density matrix is given by the quantum master ization” is done in imaginary instead of real time. equation (16). Its formal solution reads In contrast to these special cases, the von Neumann equation (t−t0)L N ρ(t) = e ρ(t0) ≡ lim (1 + δtL) ρ(t0). (19) for general mixed state density matrices cannot be rewritten in N→∞ terms of a state vector evolution, even in the case of purely co- herent Hamiltonian dynamics.4 Instead, it is necessary to study The last equality gives a meaning to the formal solution in terms the evolution of a state matrix, which transforms according to of the Trotter decomposition: at each infinitesimal time step, the the integral form of the von Neumann equation in the second exponential can be expanded to first order, such that the action line of Eq. (17). Therefore, we have to apply the Trotter formula of the Liouvillian superoperator is just given by the RHS of the and coherent state insertions on both sides of the density matrix. quantum master equation (16); At finite times, the evolved state This leads to the doubling of degrees of freedom, characteristic is given by the concatenation of the infinitesimal Trotter steps. of the Keldysh functional integral. Moreover, time evolution can If we restrict ourselves to the stationary state of the system 5 now be interpreted as occurring along the closed Keldysh con- , but want to evaluate correlations at arbitrary time differences, tour, consisting of a forward and a backward branch (cf. Fig.3 we should extend the time branch from an initial time in the dis- b)). Indeed, this is an intuitive and natural feature of evolving tant past t0 → −∞ to the distant future, t f → +∞. In analogy to matrices instead of vectors. thermodynamics, we are then interested in the so-called Keldysh partition function Z = tr ρ(t) = 1. The trace operation contracts the indices of the time evolution operator as depicted in Fig.3
4 Of course, the evolution of an M × M matrix, where M is the dimension of the Hilbert space, can be formally recast into the evolution of a vector of length M × M. While such a strategy is often pursued in numerical approaches [45], 5 We assume that it exists. We thus exclude scenarios with dynamical limit it does in general not allow for a physical interpretation. cycles, for simplicity. 13 c), giving rise to the closed time path or Keldysh contour. Con- first on a single time step, as in the usual derivation of the co- servation of probability in the quantum mechanical system is herent state functional integral [171]. That is, we decompose the reflected in the time independent normalization of the partition time evolution from t0 to t f into a sequence of small steps of du- function. In order to extract physical information, again in anal- ration δt = (t f − t0)/N, and denote the density matrix after the ogy to statistical mechanics, below we introduce sources in the n-th step, i.e., at the time tn = t0 + δtn, by ρn = ρ(tn). We then partition function. This allows us to compute the correlation and have response functions of the system by taking suitable variational derivatives with respect to the sources. δtL 1 2 ρn+1 = e ρn = ( + δtL) ρn + O(δt ). (20) After this qualitative discussion, let us now proceed with the explicit construction of the Keldysh functional integral for open As anticipated above, we proceed to represent the density matrix systems, starting from the master equation (16). As mentioned in the basis of coherent states. For instance, ρn at the time tn can above, the Keldysh functional integral is an unraveling of Li- be written as ouvillian dynamics in the basis of coherent states, and we first collect a few important properties of those. Coherent states are Z ∗ ∗ dψ+,ndψ dψ−,ndψ− ∗ ∗ +,n ,n −ψ ψ+,n−ψ− ψ−,n defined as (for simplicity, in the present discussion we restrict ρn = e +,n ,n ourselves to a single bosonic mode b) |ψi = exp(ψb†) |Ωi, where π π |Ωi represents the vacuum in Fock space. A key property of co- × hψ+,n|ρn|ψ−,ni |ψ+,ni hψ−,n| . (21) herent states is that they are eigenstates of the annihilation oper- ator, i.e., b |ψi = ψ |ψi, with the complex eigenvalue ψ. Clearly, As a next step, we would like to express the matrix element † ∗ 6 this implies the conjugate relation hψ| b = hψ| ψ . The overlap hψ+,n+1|ρn+1|ψ−,n+1i, which appears in the coherent state repre- ψ∗φ of two (non-normalized) coherent states is given by hψ|φi = e , sentation of ρn+1, in terms of the corresponding matrix element 1 R dψdψ∗ −ψ∗ψ and the completeness relation reads = π e |ψi hψ|. at the previous time step tn. Inserting Eq. (21) in Eq. (20), we The starting point of the derivation is Eq. (19), and we focus find that this requires us to evaluate the “supermatrixelement”