Keldysh approach for nonequilibrium phase transitions in quantum optics: Beyond the Dicke model in optical cavities
The Harvard community has made this article openly available. Please share how this access benefits you. Your story matters
Citation Torre, Emanuele, Sebastian Diehl, Mikhail Lukin, Subir Sachdev, and Philipp Strack. 2013. “Keldysh approach for nonequilibrium phase transitions in quantum optics: Beyond the Dicke model in optical cavities.” Physical Review A 87 (2) (February 21). doi:10.1103/ PhysRevA.87.023831. http://dx.doi.org/10.1103/PhysRevA.87.023831.
Published Version doi:10.1103/PhysRevA.87.023831
Citable link http://nrs.harvard.edu/urn-3:HUL.InstRepos:11688797
Terms of Use This article was downloaded from Harvard University’s DASH repository, and is made available under the terms and conditions applicable to Open Access Policy Articles, as set forth at http:// nrs.harvard.edu/urn-3:HUL.InstRepos:dash.current.terms-of- use#OAP Keldysh approach for non-equilibrium phase transitions in quantum optics: beyond the Dicke model in optical cavities
1, 2, 3 1 1 1 Emanuele G. Dalla Torre, ⇤ Sebastian Diehl, Mikhail D. Lukin, Subir Sachdev, and Philipp Strack 1Department of Physics, Harvard University, Cambridge MA 02138 2Institute for Quantum Optics and Quantum Information of the Austrian Academy of Sciences, A-6020 Innsbruck, Austria 3Institute for Theoretical Physics, University of Innsbruck, A-6020 Innsbruck, Austria (Dated: January 16, 2013) We investigate non-equilibrium phase transitions for driven atomic ensembles, interacting with a cavity mode, coupled to a Markovian dissipative bath. In the thermodynamic limit and at low-frequencies, we show that the distribution function of the photonic mode is thermal, with an e↵ective temperature set by the atom-photon interaction strength. This behavior characterizes the static and dynamic critical exponents of the associated su- perradiance transition. Motivated by these considerations, we develop a general Keldysh path integral approach, that allows us to study physically relevant nonlinearities beyond the idealized Dicke model. Using standard di- agrammatic techniques, we take into account the leading-order corrections due to the finite number of atoms N. For finite N, the photon mode behaves as a damped, classical non-linear oscillator at finite temperature. For the atoms, we propose a Dicke action that can be solved for any N and correctly captures the atoms’ depolarization due to dissipative dephasing.
I. INTRODUCTION portional to the atom-photon interaction strength. The static and dynamic critical exponents at the high-temperature phase Much interest has recently been directed towards under- transitions are analogous to those in a conventional laser (or, standing many-body dynamics in open systems away from more precisely, optical parametric oscillator) threshold. thermal equilibrium. This subject is not new, as the analogies To treat the interplay between external driving, dissipation between threshold phenomena in dynamical systems, such as and many-body interactions in a more general setting, we next the laser, and the conventional phase transitions have been rec- develop a unified approach to describe phase transitions in ognized over 40 years ago. However, recent experiments with open quantum systems based on Keldysh path integrals [15– ultracold atoms in optical cavities o↵er intriguing possibili- 21]. This approach is used to analyze the driven Dicke model ties to explore the physics of strongly interacting atom-photon in the presence of finite-size e↵ects and atomic dissipative systems far away from thermal equilibrium from a new van- processes. Both perturbations are non-linear and cannot be tage point. Many fundamental concepts of condensed matter treated by the usual quantum-optical methods. Instead, we ap- physics, ranging from classification of phase transitions to the ply non-perturbative techniques, specific to the path-integral universal behavior of correlation functions in the vicinity of approach. We find that the low-frequency dynamics is ther- quantum critical points in the presence of driving and dissipa- mal even in this case, allowing for an e↵ective equilibrium tion, need to be revisited in light of these developments. description. In this paper, we investigate non-equilibrium phase tran- We expect the Keldysh approach to be directly applicable sitions for driven atomic ensembles interacting with a cavity to other dissipative models such as the recently discussed cen- mode that is subject to dissipation, focusing specifically on tral spin model [22] or fermionic lattice models [23–26]. We the dynamical superradiance transitions and associated self- believe the Keldysh calculations are not more involved, and organization of the atoms observed in Refs. [1–3]. Due to the sometimes simpler, than those of the usual quantum optics interplay of external driving, Hamiltonian dynamics and dissi- frameworks [27, 28]. At the same time, they facilitate an pative processes, the observed Dicke superradiance transitions easy comparison to other phase transitions of condensed mat- [1, 2] exhibit several properties [4–7] which are not present in ter physics. One of the objectives of the present paper is to the closed Dicke model [8–11]. Recently, it was shown that make the Keldysh approach more accessible to the broader other more interesting quantum many-body phases, such as quantum optics community. At the same time, we hope that quantum spin and charge glasses with long-range, random in- the Keldysh perspective will be helpful for condensed mat- teractions [12–14] mediated by multiple photon modes, could ter physicists to understand driven dissipative atom-photon potentially be simulated with many-body cavity QED. systems–especially in view of the qualitatively di↵erent en- arXiv:1210.3623v2 [cond-mat.quant-gas] 15 Jan 2013 In what follows we first review the physics of the non- ergy scales and bath properties in quantum optics. equilibrium Dicke transition in optical cavities with con- The paper is organized as follows. In Sec. II we intro- ventional techniques of quantum optics. Using linearized duce the Dicke model and perform a brief analysis with lin- Heisenberg-Langevin equations, we demonstrate that at low earized Heisenberg-Langevin equations pointing out that the frequencies, close to the phase transition, the system evolves relevant low-frequency correlations are thermal. In Sec. III, into a thermal state with a high e↵ective temperature, pro- we map the operators of the master equation and the associ- ated Liouvillians for dissipative processes to the field content of an equivalent real-time, dissipative Keldysh action S [a⇤, a] with a⇤, a the photon field variables. In Sec. IV we introduce ⇤ [email protected] the atomic degrees of freedom and study the thermodynamic 2 limit N of the open Dicke model. We define a non- 2. Dissipation: In addition to the coherent dynamics gener- equilibrium!1 distribution function F(!) and compare it to the ated by the Hamiltonian Eq. (1), there is a dissipative contri- equilibrium case, finding that both diverge at low-frequencies bution consisting of cavity loss and dissipative processes for as 1/!. In Sec. V we study the e↵ects of a finite size N. Com- the atoms. In the rotating frame, this vacuum is e↵ectively bining analytic and numerical methods, we derive the critical out-of-equilbrium, and leads to the non-equilibrium Marko- scaling of the photon number as function of N, and find it to be vian master equation (cf. Appendix A), equivalent to an equilibrium system at finite temperature and distinct from the zero temperature case. In Sec. VI we propose @t⇢ = i[H,⇢] + ⇢, (2) L an e↵ective method to describe the e↵ects of single atom de- Here ⇢ is the density matrix and the Liouville operator in cay across the phase transition of the Dicke model. We again L find that the distribution function is thermal, but with renor- Lindblad form malized couplings and, in particular, a di↵erent critical cou- ⇢ = ↵ 2L↵⇢L↵† L↵† L↵,⇢ , (3) pling gc. Sec. VII concludes the paper with a summary of our L ↵ { } main results and some final remarks. X ⇣ ⌘ where the curly brackets , denotes the anti-commutator and { } L↵ is a set of Lindblad or quantum jump operators. In the II. THERMAL NATURE OF THE OPEN DICKE present work we consider two types of disspative processes: TRANSITION cavity photon loss and single atom dissipative dephasing. The former is modeled by the Liouvillian: The Dicke Hamiltonian [8, 9] describes N two-level sys- x z tems or “qubits” represented by Pauli matrix operators , cav ⇢ = (2ˆa⇢aˆ† aˆ†aˆ,⇢ ) , (4) i i L { } coupled to a quantized photon mode represented by bosonic creation and annihilation operatorsa ˆ†,ˆa: where is an e↵ective decay rate (inverse lifetime) of a cavity photon of the order a few MHz [2]. Modelling the dissipative N N dynamics of the atoms depends on the specific implementa- !z z g x H = !0aˆ†aˆ + + aˆ† + aˆ . (1) tion of the driven Dicke model usually involving local pro- 2 i i i=1 pN i=1 X X ⇣ ⌘ cesses of each two-level atom separately. In Sect. VI, we account for dissipative dephasing of the atoms in an approxi- Here ! is the photon frequency, ! the level-splitting of the 0 z mate way by resorting to a simplified e↵ective low frequency qubits, and g the qubit-photon coupling, assumed to couple model. all qubits uniformly to the photon. Eq. (1) is invariant under x x an Ising-type Z2-transformation,a ˆ aˆ and i i . In the thermodynamic limit N , and! for su ciently! strong qubit-photon coupling g, the!1 ground state of Eq. (1) sponta- A. Heisenberg-Langevin analysis neously breaks this Ising symmetry and exhibits a phase tran- sition to a “superradiant” phase with a photon condensate aˆ . We now study the above non-equilibrium Dicke model h i In the context of ultracold dilute gases in optical cavities, using conventional quantum optical techniques, namely, the Dimer et al. [4] proposed to implement the qubits using two Heisenberg-Langevin equations of motion. We will later re- hyperfine states of the atoms and showed that, close to the peat and extend these calculations using the Keldysh path in- transition, the relevant Hilbert space can be exactly mapped tegral approach, in the following Sections. The master equa- to Eq. (1). Inspired by the work of Dimer et al. [4], the tion (2–3) with Hamiltonian (1) and cavity dissipation (4) is qubit states of the Dicke model were realized using two collec- equivalent to the equations of motion: tive motional degrees of the Bose-Einstein condensate of the ˙ ig x atoms in the cavity [2, 29], see [30] for a review. In that case aˆ = i!0aˆ aˆ i + , pN F !z becomes a collective recoil frequency and the two Dicke states are components of a dynamically forming charge den- + + ig z ˙ i = i!z i (ˆa + aˆ†), sity wave. pN i This open realization of the Dicke model in optical cavities z ig + with pumped atoms is di↵erent from the closed system Dicke ˙ i = 2 ( i i )(ˆa + aˆ†). (5) pN model Eq. (1), due to the interplay of coherent drive and dis- sipation: Here the force = (t) is a stochastic Markovian operator F F 1. Coherent drive: The photon-atom coupling g de- satisfying (t) †(t0) = 2 (t t0) and (t)† (t) = 0. hF F i hF F i scribes the scattering of pump-photons and rotates, as func- This term is needed in order to preserve the commutation re- tion of time, at the pump frequency !p. To obtain the time- lation [ˆa(t), aˆ†(t)] = 1, which would otherwise exponentially independent Dicke model (1), one has to move to a rotat- decrease. See for example Ref. [32] for a detailed study of the ing frame, where the explicit time-dependence of the origi- single-atom case, N = 1. nal Hamiltonian is “gauged” away[31]. In this frame, the pa- To analyze the dynamics below the superradiance threshold rameter !0 appearing in Eq. (1) is the cavity-pump detuning in the limit of N , we assume that the atoms are fully ! = ! ! , where ! the bare cavity frequency. polarized S z = 1!1 z N/2, and neglect non-linear 0 c p c 2 i i ⇡ P 3 terms in the equations of motion. The resulting operator of the field to an equilibrium situation in the lab frame. equations can be solved exactly in the Fourier domain as has been done in detail in the work of Dimer et al. [4] and we will not repeat their calculations here. We just add B. Photon flux exponent one simple point to their comprehensive analysis, namely that the relevant low-frequency dynamics of the photons The Langevin equation (6) becomes dynamically unstable occurs in the presence of a finite e↵ective temperature. at the Dicke transition, correspondent to the point where ↵ This is most easily seen from the equation of motion for vanishes, or equivalently to the critical coupling the real-valued phase-space coordinate x(!) = (ˆa(!) + aˆ ( !))/ p2! . Defining stochastic force operators f (!) = † 0 !2 + 2 1 (!) ( i(! + !)) + ( !) ( + i(! !)) , we g = 0 ! . p2! 0 † 0 c z (11) 0 F F s 4!0 can writeh the equation of motion as i 2 Upon approaching the Dicke transition, the number of pho- 4g !z!0 ⌫x ( + i!)2 + !2 x(!) = f (!) . (6) tons diverges as g gc , where ⌫x is called the “photon 0 !2 !2 | | z ! flux exponent” [5]. For the present non-equilbrium Dicke transition, it was found [5, 6] that ⌫x = 1, in contrast to The force operator satisfies the equilibrium case of a quantum phase transition at zero ⌫ = / 2 + !2 + !2 temperature[11], where x 1 2. We now explain that this 1 0 discrepancy is due to the finite e↵ective temperature of the 2 f (!) f (!0) + f (!0) f (!) = (! + !0) . (7) h i !0 low-frequency fluctuations of the system. The photon number is related to the fluctuations of x by At low frequencies,we can neglect high-order terms in !. Eq. (6) becomes identical to the Langevin equation of a classi- 2 2 n + 1 = 2! x2 + p2 = 2! 1 + x2 , (12) cal particle in a harmonic potential with oscillation frequency h i 0 h i h i 0 !2 h i ↵, defined by 0 0 1 ⇣ ⌘ B C @B AC 2 Here we defined p = i(ˆa aˆ†)/ p2!0 and, by repeating 2 2 2 4g !0 ↵ = + !0 , (8) the above derivation of the Langevin equation, observed that !z 2 2 2 2 2 p = ( /!0) x . We can compute x using an equi- libriumh i partitionh functioni equivalent [33,h 34]i to the (low- and friction constant 2. In the same low-frequency approxi- frequency limit) of the Langevin equation (6): mation, the correlation function of the stochastic force opera- tors on the right-hand-side of Eq. (7) becomes identical to the F 1 “noise correlations” provided by an equilibrium classical bath Z = exp , with F = ↵2 x2 . (13) T e↵ 2 at a non-zero temperature x ! !2 + 2 2 Performing the Gaussian integral we obtain e↵ 0 gc Tx = = . (9) 4!0 !z 2 T e↵ 1 + = + ! x , 2 n 1 2 1 2 0 2 (14) In contrast to the temperature in an equilibrium problem, h i ! ↵ ⇠ g gc 0 0 1 | | the low-frequency e↵ective temperature here is not a global B C B C property of the system, but is in general observable depen- leading to the correct photon@ fluxA exponent ⌫x = 1. dent. In Sec. IV D, we present a systematic, generalizable The above results indicate that, from the point of view of way to extract low-frequency e↵ective temperatures based on phase transitions, it is incorrect to call the driven Dicke tran- observable-dependent fluctuation-dissipation relations. We sition a quantum phase transition even though it is “made already here quote the e↵ective temperature for the atoms (see of quantum ingredients” (two collective motional states of a Sec. VI for the details of the computation): BEC [2]). Instead, it should be regarded as a classical phase transition belonging to the dynamical universality class of 2 2 e↵ + !z the classical Ising model with no conserved quantities and T = , (10) infinite-range interactions, a mean-field version of the “Model 4!z A” of Hohenberg and Halperin [35]. The e↵ect that dissi- where !z is the recoil energy and an e↵ective single atom de- pation induces a finite e↵ective temperature is not new. In e↵ e↵ cay rate. Because T , Tx , the di↵erent parts of the driven several other condensed matter systems [36–38], the coupling system thus do not equilibrate to each other and, although the to a non-equilibrium bath typically admixes the pure many- dominant low-frequency correlations are thermal, the system body states of the closed system, transforming pure quan- is not in a global thermal state. We remark that our defini- tum phase transitions into thermal phase transitions (see also tion of e↵ective temperature does not coincide with the one [16, 17, 22, 39, 40]). What is perhaps more surprising is that commonly used in laser theory, as discussed in Sec. IV D: the neither the low-frequency e↵ective temperature for the pho- former relates the fluctuations of the field to its response in tons, nor for the atoms, is set by the cavity loss rate , but the rotating frame, while the latter compares the fluctuations instead set by the atom-photon interaction g and the e↵ective 4 single atom parameters, respectively. III. KELDYSH APPROACH FOR CAVITY VACUUM As a side remark, we note that, being complex-valued ob- jects, the photons have two normal modes. One quadrature In this section, we introduce the real-time Keldysh formal- is thermally amplified and diverges at the Dicke transition. ism and fix our notation by considering the case of a single Its orthogonal quadrature is quantum squeezed [4], remains electromagnetic mode in an open cavity (without atoms). We gapped at the transition and therefore does not influence the also indicate how to relate the operators of the master equa- thermal nature of the phase transition. These attenuated and tion (2) to the field content and choice of time contour in a amplified quadratures arise naturally as the eigenmodes of the Keldysh action (see also Refs. [15, 18, 42]). photon correlation function (see Sec. IV B). The decay of a single boson (the cavity photon) into a con- tinuum of modes (the external vacuum) is described by the master equation
@ ⇢ = i[! aˆ†aˆ,⇢] + (2ˆa⇢aˆ† aˆ†aˆ,⇢ ). (20) t 0 { } C. Dynamic critical exponent This equation results from the (fully unitary) Heisenberg equation for the coupled system-bath setting, where the sys- In addition to the photon flux exponent, we identify a sec- tem is described by the degrees of freedoma ˆ, aˆ† and the bath ond indicator of criticality, the dynamical exponent. This ex- by a continuum of harmonic modes. Eq. (20) is obtained by ponent governs the decay of the two-time correlations close eliminating (“integrating-out”) the bath in the Born-Markov to criticality. Going back to the Langevin equation (6) and and rotating wave approximation (cf. e.g. [27]). The integra- keeping the !2 terms we obtain: tion of the bath variables gives rise to an e↵ective evolution including dissipative terms. Performing the same program in V!2 + 2i! + ↵2 x(!) = f (!), (15) path integral formulation, we obtain the Markovian dissipa- 2 2 2 tive action ⇣ !0⌘ + + ! f (!) f (!0) = (! + !0) , (16) h i ! 1 0 S a = dt a+⇤ (i@t !0)a+ a⇤ (i@t !0)a (21) where we defined a dimensionless parameter V = 1 + Z 1 ⇣ 2 3 i[2a+a⇤ (a+⇤ a+ + a⇤ a )] . 4!0g /! . For simplicity we further approximate V 1 and z ⇡ obtain the correlation function [41] In the path integral formalism, the quantum mechanical⌘ opera- tors are replaced by fluctuating, time-dependent and complex- xˆ(t), xˆ(0) = iGK (t) (17) h{ }i xx valued fields (we omit the time argument for notational sim- d! (!2 + 2 + !2) plicity). The fact that the density matrix can be acted on from = i!t 0 i e 2 2 2 2 2⇡ 2!0[(2!) + (↵ ! ) ] both sides, as reflected in the Heisenberg commutator struc- Z t ture of the original evolution equation, finds its counterpart in e = m2 !2 + 2 + ↵2 cos pm2t the presence of a forward (+) and backward (-) component of ! m2 ↵2 0 8 0 the fields. The former is associated to an action on the den- h ⇣ ⌘ ⇣ ⌘ + pm2 !2 + 2 ↵2 sin pm2t , sity matrix from the left, and the latter to the right. Indeed, in 0 the first of line of Eq. (21) there is a relative minus sign be- ⇣ ⌘ ⇣ ⌘i where we defined m2 = ↵2 2. In the vicinity of the transition, tween the terms involving the two components, reflecting the for 0 <↵<, the frequency m becomes purely imaginary and Heisenberg commutator structure of Eq. (2). The terms in the the oscillatory behavior in the above expression disappears. second line instead display the characteristic Lindblad form; As already mentioned, this is a generic feature of dissipative the “jump” or “recycling” term is represented by an explicit phase transitions (see App. C). For su ciently large times and coupling of the two contours. approaching the transition ↵ 0 (where only the closest pole ! It is convenient [19–21] to introduce “center-of-mass” and to zero contributes), we then obtain “relative” field coordinates, acl = (a+ + a )/ p2, aq = (a+ !2+2 a )/ p2. These new coordinates are often referred to as “clas- 0 t/⇠t xˆ(t), xˆ(0) = 2 e . (18) h{ }i 8!0↵ sical” and “quantum” fields, because the first can acquire an The correlation time: expectation value while the second one cannot in the absence of sources. In this basis, and going to frequency space, we 2 1 write ⇠t = (19) ↵2 ⇠ g g ⌫t c A 1 | | 0[G ] (!) acl S = a⇤ , a⇤ , is governed by a dynamical exponent ⌫ = 1. a ( cl q) R 1 K (22) t ! [G ] (!) D (!) aq Z 0 1 ! The remainder of the paper is dedicated to the development B C @B 1 d! AC of a unified, and generalizable, Keldysh approach. The Dicke where we used the notation ! = 2⇡ , and acl,q(t) = i!t ! 1 model will be used as a prototypical test object and we com- ! e acl,q( ). This classical-quantumR R basis is often referred pare our results to those of other approaches, where available. to as RAK basis: the entries are the inverse Retarded (lower R 5 left) and Advanced (upper right) Green’s functions, and the the bosonic degrees of freedom inverse Keldysh component. The RAK action Eq. (22) can be easily inverted to deliver the photonic Green’s functions d! (!) = [ˆa, aˆ†] = 1 . (29) 2⇡ Aaa† h i K R A 1 1 Z G (!) G (!) 0[G ] (!) A = 1 , This “sum rule” is an exact property of the theory valid in– G (!)0 [GR] (!) DK(!) ! ! and out–of– equilibrium. In our example of one cavity mode, (23) 2 where the Keldysh Green’s functions is a matrix product (!) = . (30) Aaa† (! ! )2 + 2 0 GK(!) = GR(!)DK(!)GA(!) . (24) For the open cavity of Eq. (21) the RAK inverse Green’s func- B. Cavity correlation function tions are
1 R/A R/A K K The correlation function encodes the system’s internal cor- [G ] = ! !0 +⌃ , D =⌃ , (25) a a relations, for example the frequency-resolved photon spec- with the “self-energies” trum of the intracavity photon fields. In the steady state, the photon correlation function is related to the Keldysh Green’s ⌃A = i, ⌃R =+i, ⌃K = 2i. (26) function by a a a K It is a key property of a Markovian system that the Keldysh aa† (t) = aˆ(t), aˆ†(0) = aˆ(t)ˆa†(0) + aˆ†(0)ˆa(t) = iG (t). K C h{ }i h i component ⌃a, in Eq. (26) is frequency independent. As can (31) be seen from App. A, this is due to a separation of scales be- Here the last identity is valid only in the specific case of a tween (i) the large pump (!p) and cavity frequency (!c), both optical frequencies in the Tera Hertz range ( 1014 Hz cor- scalar Keldysh Green’s function. At equal times this relation responding to temperatures T 104K), and⇠ (ii) the charac- results in ⇠ teristic frequencies of the electromagnetic vacuum outside the K d! K 12 (0) = 2 aˆ†aˆ + 1 = iG (t = 0) = i G (!) .(32) cavity (. 10 Hz corresponding to temperatures T . 300K). Caa† h i 2⇡ Z In the literature, it is often argued that a frequency inde- For the single decaying cavity mode, characterized by pendent, nonzero inverse Keldysh component indicates an ef- Eqs. (22-26), it is easy to show that (!) = (!) and fective finite temperature state. The Markovian lossy cavity aa† aa† the frequency integral over the KeldyshC Green’sA function is is a simple counterexample: Even though the inverse Keldysh unity yielding aˆ aˆ = 0. As expected, the steady state corre- component is constant 2i, the state is pure and the e↵ective † sponds to the cavityh i vacuum. temperature zero, as we⇠ will argue below. We next introduce the key propagators that encode the sys- tems’ response and correlation functions. In equilibrium, the C. Comparison with a closed system at equilibrium two are rigidly related by the Bose (or Fermi) distribution function; out-of-equilibrium, no such a priori knowledge is available, and it is important to distinguish them. In the absence of dissipation, the cavity becomes an iso- lated harmonic oscillator. Its inverse Green’s functions are still given by (22) with the self-energies serving only as regu- larization parameters. At equilibrium,
⌃A = i✏, ⌃R =+i✏, A. Cavity spectral response function a,EQ a,EQ ! ⌃K (!) = 2i✏ coth . (33) a,EQ 2T The spectral response function encodes the system’s response to active, external perturbations such as time- Here ✏ 0 at the end of the calculation and T is the actual modulated external fields coupling to spin or charge operators, temperature.! Note that the Keldysh component is odd with for example. The spectral response function is the di↵erence K K respect to the frequency ⌃a,EQ( !) = ⌃a,EQ(!), while in the between the retarded and advanced Green’s function: Markovian system (26) it is even. Using Eq. (33), we obtain (!) = i(GR(!) GA(!)). (27) A EQ 2✏ In the scalar case considered here, we have (!) = lim = 2⇡ (! !0), (34) Aaa† ✏ 0 (! ! )2 + ✏2 ! 0 R ✏ aa† (!) = 2ImG (!), (28) EQ 2 ! ! A (!) = lim coth T = 2⇡coth T (! !0), Caa† ✏ 0 (! ! )2 + ✏2 2 2 The frequency-integrated spectral response function is nor- ! 0 malized to unity, because of the exact commutator relation of In this noninteracting case, the spectral response function is 6 fully centered at the isolated mode with frequency !0.We In case of a decaying cavity, the Green’s functions are observe that formally, the thermodynamic equilibrium limit scalars and we can easily invert (35) to obtain can be seen as a situation with an infinitesimal loss, replacing K ✏, and the replacement in the inverse Keldysh component G (!) aa† (!) ! ! F(!) = = C = 1. (37) 2i 2i✏ coth . GR(!) GA(!) (!) ! 2T Aaa† confirming that the cavity vacuum has a zero e↵ective tem- perature. Moreover, as for a pure quantum state in the equilib- ↵ D. Cavity distribution function and low-frequency e ective rium case at T = 0, the distribution function here also squares temperature to a unit matrix, F2(!) = 1.
The response and correlations allow us to define a An important di↵erence between the zero temperature equi- fluctuation-dissipation relation, by introducing the distribu- librium case and Markovian case appears in the sign of the tion function F(!): distribution function. In the former case (as for any equilib- rium distribution), F(!) is anti-symmetric with respect to the K R A G (!) = G (!)F(!) F(!)G (!) (35) frequency !. On the contrary, for a Markovian bath the distri- 1 1 DK(!) = [GR(!)] F(!) F(!)[GA(!)] , bution function is symmetric with respect to !; one signature , of a strongly-out-of-equilibrium system. where the equivalence holds due to Eq. (24). At thermal equi- librium the distribution F is universal and equals to the unit matrix times
EQ ! ! F (!) = coth 2T = 2nB( T ) + 1, EQ FT=0(!) = sign(!), IV. KELDYSH APPROACH FOR PHOTON OBSERVABLES EQ 2T F! T (!) + ... (36) ⌧ ⇡ ! 1 with the Bose distribution nB(x) = (exp x 1) . The unit ma- We now analyze the Dicke model Eq. (1) with the path trix in field space signals detailed balance between all subparts integral approach explained in the previous section. We in- of the system. clude cavity photon loss but defer the inclusion of dissipa- In the present case, the system is out-of-equilibrium due to tive processes for the atoms to Sec. VI. Assuming homoge- its driven and dissipative nature. A notion of a temperature neous qubit-cavity coupling, one can use the large-N strategy is not a priori meaningful: Neither must the driven system of Emary and Brandes [11]. We introduce collective large-N = 1 N z x = 1 N + + equilibrate to an external heat bath with temperature T (in our spin operators S z 2 i=1 i and S 2 i=1( i i ) , to write the Dicke model (1) in terms of one large spin coupled case, due to the separation of scales underlying the Markov P P approximation, this temperature would be e↵ectively T = 0 to the cavity photon mode, compared to the system scales), nor do the di↵erent subparts 2g x of the system have to equilibrate with respect to each other. In H = !0aˆ†aˆ + !zS z + S aˆ† + aˆ . (38) this work, we argue that a notion of a temperature nevertheless pN ⇣ ⌘ emerges as a universal feature of the low frequency domain We then express the spin in terms of a Holstein-Primako↵ of Markovian systems. It is introduced by computing the F [11, 43] boson operator bˆ, defined by S = N/2 + bˆ†bˆ, matrix through Eq. (35) and comparing the low-frequency be- z S + = pN nˆbˆ pN(1 + nˆ/(2N))bˆ , and S x = (S + + S )/2. havior of its eigenvalues with the equilibrium result of Eq. † † Neglecting unimportant⇡ constants we obtain the normal or- (36). In particular, if F has a thermal infrared enhancement dered Hamiltonian: 1/! for small frequencies, its dimensionful coe cient is ⇠ ↵ identified as an e ective temperature. This notion of a “low H = !0aˆ†aˆ + !zbˆ†bˆ (39) frequency e↵ective temperature” (LET) becomes particularly 1 +g aˆ + aˆ† bˆ + bˆ† bˆ†(bˆ† + bˆ)bˆ . relevant in the vicinity of a phase transition, where the spec- 2N R, A tral weight encoded in G G is concentrated near zero fre- At N , the⇣ last, non-quadratic⌘⇣ term vanishes⌘ and the quency. Below, we will use this concept to establish a con- problem!1 reduces to a linear system of two coupled bosonic de- nection of Markovian quantum systems to the classical the- grees of freedom one of which (the cavity modea ˆ) decays into ory of dynamical universality classes according to Hohenberg a Markovian bath. As outlined in section III, we can trans- and Halperin [35]. Moreover, we find that, while all governed form the Liouvillian Eq. (4), with Hamiltonian Eq. (39), into /! by the 1 divergence in the distribution function, di↵erent an equivalent Keldysh action with subparts of the system exhibit di↵erent LETs. In contrast to the global temperature present in thermodynamic equilibrium, S = S a + S b + S ab (40) the LET is not an external parameter but rather a system im- manent quantity, determined by the interplay of unitary and 0 ! !z bcl S b = (bcl⇤ , bq⇤) ! !z 0 bq dissipative dynamics. Z! ! ! 7 where S a is given by Eq. (22) and the interaction in terms of with the action the “classical” and “quantum” fields reads 1 A 1 0[G4 4] (!) S N = V†(!) 1 ⇥ V8(!) . 8 R K 2 ! [G4 4] (!) D4 4 S ab = g [(aq + aq⇤)(bcl + bcl⇤ ) + (acl + acl⇤ )(bq + bq⇤)] Z 0 ⇥ ⇥ 1 ! B C (47) Z @B AC 1 (a + a⇤ )(b + b⇤) + (a + a⇤)(b + b⇤ ) b⇤ b + b⇤b The subscript N stands for normal phase and the dagger de- cl cl q q q q cl cl cl cl q q 4N notes transposition and complex conjugation. The block† en- ( h ih i tries are 4 4 Green’s functions given by + (acl + acl⇤ )(bcl + bcl⇤ ) + (aq + aq⇤)(bq + bq⇤) bcl⇤ bq + bq⇤bcl . ⇥ ) ! ! + i 0 g g h ih (41)i 0 1 0 ! ! i g g R 0 [G4 4] (!) = 0 1 We now demonstrate that the static saddle point solutions of ⇥ g g ! !z 0 B C this action reproduce the results of other approaches [4]. Vary- B g g 0 ! !z C B C ing S with respect to the quantum components of the fields K B C D4 4 = 2iB diag(,, 0, 0). (48) C ⇥ @ A and substituting acl(t) = p2a0, bcl(t) = p2b0, aq = 0, bq = 0, we obtain the coupled equations To obtain the photon-only action, we now integrate out the Holstein-Primako↵ field b and get a Keldysh func- @S 1 K = ( ! + i)a g 1 b2 2b = 0 , (42) tional integral that goes only over the photon fields: Z = 0 0 0 0 iS [a ,a] @a⇤ 2N D a , a e photon ⇤ . with the photon-only action q ! { ⇤ } @S 3 2 1 = ! b g 1 b (a + a⇤) = 0 . (43) R A z 0 0 0 0 0[G2 2] (!) @bq⇤ 2N S photon[a⇤, a] = A†(!) 1 ⇥ A4(!) . ! 4 R K ! [G2 2] (!) D2 2(!) Z 0 ⇥ ⇥ 1 where we chose b = b⇤. These saddle-point equations ad- B C 0 0 B C (49) mit solutions with non-zero “ferromagnetic” moment b0, and @ A super-radiant photon condensate a0, The photon four-vector collects the classical and quantum field components 2 2 2 2 N g gc g gc b = a = p2N (44) acl(!) 0 ± 2 g2 0 ± ! i r s p 0 a⇤ ( !) A (!) = cl . (50) 4 0 a (! ) 1 for atom-photon couplings larger than a critical value B q C B a⇤( !) C B q C B C 2 2 B C !0 + and the block entries are 2 2@ photon Green’sA functions which g > gc = !z , (45) ⇥ s 4!0 we now analyze one-by-one. in agreement with the results known from the literature[5, 6]. We are now going to integrate-out the atomic field b and ob- tain an e↵ective action describing the photons in the normal A. Photon spectral response function phase (g < gc), where a0 = b0 = 0. In App. B, we give the corresponding expressions in the superradiant phase. The inverse retarded Green’s function of the photons is
1 R [G2 2] (!) = ⇥ ! ! + i +⌃R(!) ⌃R(!) 0 , ⌃R( !) ⇤ ! ! i + ⌃R( !) ⇤ 0 ! In the limit of N , we can safely neglect the terms h i h i (51) proportional to 1/N in!1 (40). We can re-write Eq. (40) as 8 8 ⇥ matrix multiplying the 8-component fields: where the interaction induced photon-self energy reads
2 acl(!) 2g ! ⌃R(!) = z . (52) a⇤ ( !) 2 2 cl ! !z 0 bcl(!) 1 B C B bcl⇤ ( !) C The characteristic frequencies of the system are defined by V8(!) = B C (46) 1 B aq(!) C R B C the zeros of the determinant [G2 2] (!), corresponding to the B ! C ⇥ B aq⇤( ) C poles of the response function GR (!). Due to the symmetry B C 2 2 B bq(!) C ⇥ B C B bq⇤( !) C R ⇤ R B C x G2 2( !) x = G2 2(!) , (53) B C ⇥ ⇥ @B AC h i 8
15 � � = � � � ≪ �� Im[�] Im[�] 10 5 −� −� � � −� −� � � " Ω !
Re[�] Re[�] ! 0 aa A !5 −� −� !10 !15 � � ≾ �� � � = �� !2 !1 0 1 2 Im[�] Im[�] Ω
−� −� � � −� −� � � Re[�] Re[�] 15
−� −� 10 " Ω ! ! aa
C 5 FIG. 1. (Color online) Schematic plot of the position of the poles of the retarded Green’s function Eq. (54). (a) At zero coupling g = 0, two poles can be associated with the photonic mode ! = ! + i 0 ± 0 and two with the atomic mode ! = ! . (b) In the presence of ± z !2 !1 0 1 2 a finite coupling g, the modes hybridize and the corresponding fre- quencies are shifted in opposite directions. (c) When approaching the Ω transition, two solutions become purely imaginary and correspond to damped modes. (d) At the transition point g = gc one of the poles A ! FIG. 2. (Color online) Photon spectral response function aa† ( ) approaches zero, making the system dynamically unstable. C ! ! and correlation function aa† ( ) as a function of real frequencies . Numerical parameters: ! = ! = 1,= 0.2 (leading to g 0.51), 0 z c ⇡ and g=0 (dotted), 0.25 (dashed), 0.5 (solid). the poles come in pairs, such that = , meaning that { } { ⇤} they either are pure imaginary or come in pairs with opposite B. Photon correlation function real part. The explicit solution of
1 2g2! 2 2g2! 2 The inverse Keldysh component of the action Eq. (49) is 0 = = ! + z (!+i)2 z R 0 !2 !2 !2 !2 det[G2 2(!)] z z ⇥ ! ! 2i 0 (54) DK = (55) 2 2 02i yields four poles, schematically plotted in Fig. 1. Note that, in ⇥ ! the vicinity of the phase transition, two poles become purely K ! = imaginary. This phenomenon seems to apply to generic dissi- and the Keldysh Green’s function G2 2( ) R K A ⇥ pative transitions as we further describe in App. C. Indeed it G2 2(!)D2 2G2 2(!) is a 2 2 matrix. ⇥ ⇥ ⇥ ⇥ has been previously observed for a dissipative critical central K = K T The first diagonal element of G2 2, aa† i(1 0)G2 2(1 0) spin model [22]. Overdamping of collective modes, due to a corresponds to the photon correlation⇥ C function defined⇥ in similar mechanism, has also been found in dissipative multi- Eq. (31) and is plotted in Fig. 2. As noted above in Eq. (32), mode systems in symmetry broken phases [16, 17]. its frequency integral gives the steady-state photonic occupa- tion and it can be shown to diverge at the Dicke transition The imaginary part of the first diagonal element of GR (!) 2 2 according to Eq. (12). This result will be explicitly derived in corresponds to the photon spectral response function ⇥ , aa† Sec. V C using a low-frequency e↵ective description of GK. defined in (27), and is plotted in Fig. 2 for di↵erent valuesA The photon number diverges in steady state despite the fact of the coupling g. In the absence of atom-photon coupling that the system undergoes photon loss, and no explicit photon g = 0 (dotted curve) there is a single resonance peak at fre- pumping occurs within the model. The reason is that the cou- quency ! = ! , broadened by the cavity decay rate . A finite 0 pling constant g is an e↵ective parameter, which in any con- coupling g (dashed curve) “collectively Rabi-splits” the reso- crete physical realization microscopically involves a coherent nance [44] in two distinct peaks, corresponding to two distinct laser drive process compensating for the loss. poles of the system. Upon approaching the Dicke transition (solid curve), the spectral weight is shifted towards the low- The matrix structure of the Keldysh Green’s function K frequency pole; a precursor to the superradiant cavity mode. G2 2(!) can be conveniently exploited to compute the quadra- ⇥ 9 ture fluctuations for a general phase angle ✓ ! of the order of the atomic detuning !z and vanishes for large frequencies. i✓ i i✓ i✓ K e x (!)x ( !) = e e G (!) , (56) Nevertheless, also in the driven Dicke model, the experi- ✓ ✓ 2 2 i✓ out h i 4!0 ⇥ e mentally relevant homodyne spectrum G (!) of the cavity ! 2 2 ⇣ ⌘ output field shows noise reduction below the⇥ vacuum level in where x✓ is defined by the ✓⇤ quadrature [4]. Following the standard “input-output theory” [45, 47], it is possible to show that the homodyne 1 i✓ i✓ x✓ = e a + e a† . (57) spectrum can be linked to the Keldysh response function via p2!0 ⇣ ⌘ The corresponding equal time, frequency-integrated fluctua- 2 2i 0 R tions x✓(t)x✓(t) = x✓(!)x✓( !) are plotted in Fig. 3 and out ! i G2 2(!) = G2 2(!) 12 2 12 2, (59) h i h i 0 2i ⇥ ⇥ ⇥ diverge at the Dicke transition for all angles ✓ except for ✓⇤ ⇥ ! R defined by 2 where M M† M. By applying the transformation (56) 1 out|| || ⌘ ✓⇤ = ⇡ tan (!0/) . (58) to G2 2, it is then possible to compute the fluctuations of the output⇥ quadrature x (!)x ( !) . This quantity has h out,✓ out,✓ i The angle ✓⇤ can also be obtained as the phase angle for the been studied in detail in Ref. [4]: at the Dicke transition, non-diverging eigenmode of the zero-frequency limit of the the zero-frequency component of ✓ = ✓⇤ tends to the maxi- Keldysh correlation function, thereby naturally yielding the mally attenuated value of 1/(4!0). It should be noted, how- attenuated and amplified quadratures alluded to in Sec. II. ever, that the equal-time, frequency integrated fluctuations
equal-time xout,✓⇤ (t)xout,✓⇤ (t) are zero and therefore not attenuated below Here in the case of the driven Dicke model, the h i photon fluctuations of x✓ inside the cavity are independent of vacuum level. the atom-photon coupling⇤ and in particular are not attenuated below the vacuum noise level (the g = 0 limit, without atoms in the cavity). This is di↵erent from the case of the optical C. Comparison with a closed system at equilibrium parametric oscillator [45, 46] where, at threshold, the equal- time fluctuations of the non-diverging intra-cavity quadrature are reduced to 50% of the vacuum level. This di↵erence can In the case of a closed system at equilibrium, the retarded be traced back the di↵erent frequency dependencies of the ef- Green’s function is the same as Eq. (51), up to the replace- fective driving term in both situations. In the parametric oscil- ment ✏, and letting ✏ 0. The corresponding spectral response! function (!) is! similar to the one shown for the lator, the driving of the cavity occurs via a classically-treated Aaa† photon pump laser, and the driving amplitude is typically set Markovian case, with the narrow peaks substituted by delta- to a constant coherent field amplitude. In the present case of functions at the resonant frequencies. The Keldysh compo- the driven Dicke model, the e↵ective driving of the cavity is nent of the inverse Green’s function reads, at zero tempera- mediated by the atoms via virtual absorption and emission of ture: photons from the pump laser into the cavity. The correspond- 2 2 2 K 2i✏sign(!)0 ing driving term g !z/(!z ! ) is maximal for frequencies D2 2, EQ = . (60) ⇠ ⇥ 0 2i✏sign(!) ! Again note the di↵erent symmetry under frequency reflection 1.0 with respect to (55). 0.8
0.6 # D. Photon distribution function and low-frequency e↵ective 2 Θ
x temperature
" 0.4
0.2 We now show that the above mentioned di↵erence between the Markovian case and the equilibrium case leads to the gen- 0.0 eration of a “low frequency e↵ective temperature” (LET) for 0.0 0.5 1.0 1.5 2.0 2.5 3.0 the former. For this purpose, we calculate the distribution ma- Θ trix F, defined in Eq. (35). To this end, recall that at thermal equilibrium, F = coth(!/2T)1: it exponentially approaches FIG. 3. Frequency integrated equal-time correlations (x (t))2 as unity at high frequencies ( ! T), and diverges as 2T/! at h ✓ i | | function of the angle ✓, for di↵erent values of the coupling strength low frequencies ( ! T). For our Markovian problem, using | |⌧ approaching the transition at gc 0.51: g=0 (dotted), 0.25 (dashed), Eq. (49) we find: ⇡ 1 0.5 (solid). The fluctuations in the quadrature ✓⇤ = ⇡ tan (! /) 0 ⇡ 1.768 are independent of the coupling strength. Numerical parame- 2 g2! F = + z , (61) ters: !0 = !z = 1.0, = 0.2 z ! !2 !2 x z 10
4 spectral properties in this regime (cf. Fig. 2). At higher fre- quencies, it finally tends to one, following the non-equilibrium 3 curve f Te↵/(! !z). The approach to the quantum regime f 1 is⇡ polynomial, unlike the exponential approach in the " equilibrium⇡ case. In App. D, we show that the thermal 1/! Ω
! 2 #
f divergence, leading to a finite LET, is generic for Markovian systems. 1 We note that our definition of an e↵ective temperature is not the one commonly used in the context of laser theory [28]. 0 In this context, the e↵ective temperature is used only to de- !2 !1 0 1 2 scribe the fluctuations of the photonic field, as compared to the Ω equilibrium fluctuations in the lab frame. As a consequence, the divergence of the photon number at the phase transition FIG. 4. (Color online) Positive eigenvalue of the distribution func- is always associated to a diverging e↵ective temperature. In tion F(!) for the same parameters as in Fig. 2. At low frequencies contrast, our low-frequency e↵ective temperature (LET) de- the distribution diverges as 2Te↵ /!, where the finite e↵ective tem- scribes the ratio between the fluctuations and the response of 2 perature is proportional to the photon-atom interaction, Te↵ g . the system in the rotating frame. It is finite at the transition ⇠ and, as we will see, allows us to map the Dicke transition to an existing dynamical universality class of equilibrium systems. where z and x are Pauli matrices. The F matrix is hermitian and traceless, so its two eigenvalues are real and opposite: V. FINITE-N CORRECTIONS FROM A KELDYSH AND 2g2/! 1 2 LANGEVIN PERSPECTIVE ! = + z . f ( ) 1 2 (62) ± ± s ! 1 (!/!z) ! We now move beyond the quadratic theory, by consider- Fig. 4 shows the behavior of the positive eigenvalue in two ing the e↵ects of finite N. Based on the formalism developed points of the phase diagram, nearby and far away from the in the previous section, we approach this problem by scal- transition. In both cases, at low frequencies the eigenvalue ing analysis, diagrammatic technique, mapping into a low- diverges as 1/!. Exploiting the analogy to the equilibrium frequency e↵ective Langevin equation. As will be seen be- case, we obtain the LET low, these methods show quantitative agreement with a Monte
2 Carlo solution of the original Master equation, highlighting g the utility of the present formalism. Te↵ = . (63) !z
We find that Te↵ is not proportional to the decay rate . The temperature is rather proportional to the e↵ective interaction A. Scaling analysis between the spin and the photon. g is the scale that leads to a competition of unitary and dissipative dynamics, and the Up to this point we have considered only the thermody- LET is a measure of this: For g = 0, the steady state of the namic limit N . In this limit, the resulting theory is dissipative part of the dynamics (empty cavity) is an eigenstate quadratic and can!1 be studied by Keldysh means as well as by of the Hamiltonian, while this is no longer the case for any the Heisenberg-Langevin method. Corrections due to a finite finite g. In the limiting case of g 0, the LET goes to zero. ! N introduce non-quadratic terms into the problem and require A closer inspection of (62) reveals an important di↵erence a more careful study. The present path-integral approach al- between the Markovian bath and thermal equilibrium, related lows us to develop a diagrammatic approach and to resum all to the presence of a second energy scale, !z. If !z Te↵, this leading-order corrections in an organized fashion. A similar energy scale does not a↵ect the crossover between the quan- approach has been used to study the instability of an optical tum and classical regimes, which then proceeds monotonously parametric oscillator in Ref.s [46, 48, 49], where however the similar to an equilibrium problem. If, on the other hand, emerging low-frequency thermal nature of the problem has !z Te↵, the quantum-classical crossover of the Marko- not been discussed. vian⌧ bath occurs in an unusual way, highlighting the non- Leading 1/N corrections to the Hamiltonian of the Dicke equilibrium nature of the problem. Starting from a divergence model are easily obtained by retaining the first-order terms in at zero frequency, the distribution function (62) decreases as the Holstein-Primako↵ approximation (see Eq. (39)). These Te↵/!, in analogy to an equilibrium system at finite temper- terms are expressed using the Keldysh formalism in Eq. (41) ature. Then, instead of monotonously decreasing towards the and contain products of four fields. In a diagrammatic descrip- quantum regime where f 1, it exhibits a second divergence tion (see Fig. 5), they correspond to four-point vertices. These ⇡ at ! = !z. Since the spectral weight vanishes su ciently fast vertices can be “classical” if they contain only one quantum in this regime, the correlation functions still remain finite; the field (either bq or aq), or “quantum” if they contain three of pole in F accounts for a di↵erent scaling of correlations and them. The former type can be casted into a semi-classical 11
that the appropriate transformation is: (a) (c) quantum (d) Σ (b) classical × N N N0 = (68) � 2 � � !