Dicke model


The Dicke model is a fundamental model of quantum optics, which describes the interaction between light and matter. In the Dicke model, the light component is described as a single quantum mode, while the matter is described as a set of two-level systems. When the coupling between the light and matter crosses a critical value, the Dicke model shows a mean-field phase transition to a superradiant phase. This transition belongs to the Ising universality class and was realized in cavity quantum electrodynamics experiments. Although the superradiant transition bears some analogy with the lasing instability, these two transitions belong to different universality classes.



The Dicke model is a quantum mechanical model that describes the coupling between a single-mode cavity and   two-level systems, or equivalently   spin-1/2 degrees of freedom. The model was first introduced in 1973 by K. Hepp and E. H. Lieb.[1] Their study was inspired by the pioneering work of R. H. Dicke on the superradiant emission of light in free space [2] and named after him.

Like any other model in quantum mechanics, the Dicke model includes a set of quantum states (the Hilbert space) and a total-energy operator (the Hamiltonian). The Hilbert space of the Dicke model is given by (the tensor product of) the states of the cavity and of the two-level systems. The Hilbert space of the cavity can be spanned by Fock states with   photons, denoted by  . These states can be constructed from the vacuum state   using the canonical ladder operators,   and  , which add and subtract a photon from the cavity, respectively. The states of each two-level system are referred to as up and down and are defined through the spin operators  , satisfying the spin algebra  . Here   is the reduced Planck constant and   indicates a specific two-level system.[3]

The Hamiltonian of the Dicke model is



Here, the first term describes the energy of the cavity and equals to the product of the energy of a single cavity photon   (where   is the cavity frequency), times the number of photons in the cavity,  . The second term describes the energy of the two-level systems, where   is the energy difference between the states of each two-level system. The last term describes the coupling between the two-level systems and the cavity and is assumed to be proportional to a constant,  , times the inverse of the square root of the number of two-level systems. This assumption allows one to obtain a phase transition in the limit of   (see below). The coupling can be written as the sum of two terms: a co-rotating term that conserves the number of excitations and is proportional to   and a counter-rotating term proportional to  , where   are the spin ladder operators.

The Hamiltonian in Eq. 1 assumes that all the spins are identical (i.e. have the same energy difference and are equally coupled to the cavity). Under this assumption, one can define the macroscopic spin operators  , with  , which satisfy the spin algebra,  . Using these operators, one can rewrite the Hamiltonian in Eq. 1 as


This notation simplifies the numerical study of the model because it involves a single spin-S with  , whose Hilbert space has size  , rather than   spin-1/2, whose Hilbert space has size  .

The Dicke model has one global symmetry,



Because   squares to unity (i.e. if applied twice, it brings each state back to its original state), it has two eigenvalues,   and  . This symmetry is associated with a conserved quantity: the parity of the total number of excitations,  , where



This parity conservation can be seen from the fact that each term in the Hamiltonian preserves the excitation number, except for the counter-rotating terms, which can only change the excitation number by  . A state of the Dicke model is said to be normal when this symmetry is preserved, and superradiant when this symmetry is spontaneously broken.


The Dicke model is closely related to other models of quantum optics. Specifically, the Dicke model with a single two-level system,  , is called the Rabi model. In the absence of counter-rotating terms, the model is called Jaynes-Cummings for   and Tavis-Cummings for  . These two models conserve the number of excitations   and are characterized by a   symmetry. The spontaneous breaking of this symmetry gives rise to a lasing state (see below).

The relation between the Dicke model and other models is summarized in the table below [4]

Model's name Counter-rotating terms? symmetry Number of two-level systems
Jaynes-Cummings no    
Tavis-Cummings no    
Rabi model yes    
Dicke yes    

Superradiant phase transition

Schematic plot of the order parameter of the Dicke transition, which is zero in the normal phase and finite in the superradiant phase. The inset shows the free energy in the normal and superradiant phases, see Eq. 5.

Early studies of the Dicke model considered its equilibrium properties.[1] These works considered the limit of   (also known as the thermodynamic limit) and assumed a thermal partition function,  , where   is the Boltzmann constant and   is the temperature. It was found that, when the coupling   crosses a critical value  , the Dicke model undergoes a second-order phase transition, known as the superradiant phase transition. In their original derivation, Hepp and Lieb[1] neglected the effects of counter-rotating terms and, thus, actually considered the Tavis-Cummings model (see above). Further studies of the full Dicke model found that the phase transition still occurs in the presence of counter-rotating terms, albeit at a different critical coupling.[5]

The superradiant transition spontaneously breaks the parity symmetry,  , defined in Eq. 3. The order parameter of this phase transition is  . In the thermodynamic limit, this quantity tends to zero if the system is normal, or to one of two possible values, if the system is superradiant. These two values correspond to physical states of the cavity field with opposite phases (see Eq. 3 and, correspondingly, to states of the spin with opposite   components). Close to the superradiant phase transition, the order parameter depends on   as  . This dependence corresponds to the mean-field critical exponent  .

Mean-field description of the transition


The simplest way to describe the superradiant transition is to use a mean-field approximation, in which the cavity field operators are substituted by their expectation values. Under this approximation, which is exact in the thermodynamic limit, the Dicke Hamiltonian of Eq. 1 becomes a sum of independent terms, each acting on a different two-level system, which can be diagonalized independently. At thermal equilibrium (see above), one finds that the free energy per two-level system is[6]



The critical coupling of the transition can be found by the condition  , leading to



For  ,   has one minimum, while for  , it has two minima. In the limit of   one obtains an expression for the critical coupling of the zero-temperature superradiant phase transition,  .

Semiclassical limit and chaos


Semiclassical limit


A phase space for the Dicke model in the symmetric atomic subspace with   may be constructed by considering the tensor product of the Glauber coherent states


where   is the displacement operator and   is the photon vacuum Fock state, and the SU(2) coherent states


where   is the rotation operator in the Bloch sphere,   and   is the state with all atoms in their ground state. This yields a four-dimensional phase space with canonical coordinates   and  .

A classical Hamiltonian is obtained by taking the expectation value of the Dicke Hamiltonian given by Eq. 2 under these states,[7][8]

Percentage of classical trajectories with positive Lyapunov exponent as a function of the energy per particle   and the coupling parameter   (divided by the critical coupling  ). The parameters are  .

In the limit of  , the quantum dynamics given by the quantum Hamiltonian of Eq. 2 and the classical dynamics given by Eq. 9 coincide. For a finite system size, there is a classical and quantum correspondence that breaks down at the Ehrenfest time, which is inversely proportional to  .

Quantum chaos


The Dicke model provides an ideal system to study the quantum-classical correspondence and quantum chaos.[9]

The classical system given by Eq. 9 is chaotic or regular depending on the values of the parameters  ,  , and   and the energy  .[8][10] Note that there may be chaos in both the normal and superradadiant regimes.

It was recently found that the exponential growth rate of the out-of-time-order correlator coincides with the classical Lyapunov exponents[11][12] in the chaotic regime and at unstable points of the regular regime. In addition, the evolution of the survival probability (i.e. the fidelity of a state with itself at a later time) of initial coherent states highly delocalized in the energy eigenbasis is well-described by random matrix theory,[13][14] while initial coherent states strongly affected by the presence of quantum scars display behaviors that break ergodicity.[15][16]

Open Dicke model


The Dicke model of Eq. 1 assumes that the cavity mode and the two-level systems are perfectly isolated from the external environment. In actual experiments, this assumption is not valid: the coupling to free modes of light can cause the loss of cavity photons and the decay of the two-level systems (i.e. dissipation channels). It is worth mentioning, that these experiments use driving fields (e.g. laser fields) to implement the coupling between the cavity mode and the two-level systems. The various dissipation channels can be described by adding a coupling to additional environmental degrees of freedom. By averaging over the dynamics of these external degrees of freedom one obtains equations of motion describing an open quantum system. According to the common Born-Markov approximation, one can describe the dynamics of the system with the quantum master equation in Lindblad form [17]



Here,   is the density matrix of the system,   is the Lindblad operator of the decay channel  , and   the associated decay rate. When the Hamiltonian   is given by Eq. 1, the model is referred to as the open Dicke model.

Some common decay processes that are relevant to experiments are given in the following table:

- Cavity decay Atomic decay Atomic dephasing Collective decay
Decay rate        

In the theoretical description of the model, one often considers the steady state where  . In the limit of  , the steady state of the open Dicke model shows a continuous phase transition, often referred to as the nonequilibrium superradiant transition. The critical exponents of this transition are the same as the equilibrium superradiant transition at finite temperature (and differ from the superradiant transition at zero temperature).

Superradiant transition and Dicke superradiance

Schematic representation of the difference between Dicke superradiance and the superradiant transition of the open Dicke model.

The superradiant transition of the open Dicke model is related to, but differs from, Dicke superradiance.

Dicke superradiance is a collective phenomenon in which many two-level systems emit photons coherently in free space.[2][18] It occurs if the two-level systems are initially prepared in their excited state and placed at a distance much smaller than the relevant photon's wavelength. Under these conditions, the spontaneous decay of the two-level systems becomes much faster: the two-level systems emit a short pulse of light with large amplitude. Under ideal conditions, the pulse duration is inversely proportional to the number of two-level systems,  , and the maximal intensity of the emitted light scales as  . This is in contrast to the spontaneous emission of   independent two-level systems, whose decay time does not depend on   and where the pulse intensity scales as  .

As explained above, the open Dicke model rather models two-level systems coupled to a quantized cavity and driven by an external pump. In the normal phase, the intensity of the cavity field does not scale with the number of atoms  , while in the superradiant phase, the intensity of the cavity field is proportional to  .

The scaling laws of Dicke superradiance and of the superradiant transition of the Dicke model are summarized in the following table:

Dicke superradiance[2] Superradiant transition of the Dicke model[1]
Environment Free space Cavity
Duration Transient Steady state
Intensity of the field (normal)    
Intensity of the field (superradiant)    

Experimental realizations

Schematic representation of two schemes to experimentally realize the Dicke model: on the left, the equilibrium approach based on the dipole coupling between the two levels and, on the right, the nonequilibrium approach based on two-photon processes, namely stimulated Raman scattering. Only the latter scheme is used to realize the Dicke model.

The simplest realization of the Dicke model involves the dipole coupling between two-level atoms in a cavity. In this system, the observation of the superradiant transition is hindered by two possible problems: (1) The bare coupling between atoms and cavities is usually weak and insufficient to reach the critical value  , see Eq. 6.[19] (2) An accurate modelling of the physical system requires to consider   terms that according to a no-go theorem, may prevent the transition. Both limitations can be circumvented by applying external pumps on the atoms and creating an effective Dicke model in an appropriately rotating frame.[20][21]

In 2010, the superradiant transition of the open Dicke model was observed experimentally using neutral Rubidium atoms trapped in an optical cavity.[22] In these experiments, the coupling between the atoms and the cavity is not achieved by a direct dipole coupling between the two systems. Instead, the atoms are illuminated by an external pump, which drives a stimulated Raman transition. This two-photon process causes the two-level system to change its state from down to up, or vice versa, and emit or absorb a photon into the cavity. Experiments showed that the number of photons in the cavity shows a steep increase when the pump intensity crosses a critical threshold. This threshold was associated with the critical coupling of the Dicke model.

In the experiments, two different sets of physical states were used as the down and up states. In some experiments,[23][22][24] the two states correspond to atoms with different velocities, or momenta: the down state had zero momentum and belonged to a Bose-Einstein condensate, while the up state had a momentum equal to sum of the momentum of a cavity photon and the momentum of a pump photon.[25] [26] In contrast, later experiments[27][28] used two different hyperfine levels of the Rubidium atoms in a magnetic field. The latter realization allowed the researchers to study a generalized Dicke model (see below). In both experiments, the system is time-dependent and the (generalized) Dicke Hamiltonian is realized in a frame that rotates at the pump's frequency.

Generalized model and lasing


The Dicke model can be generalized by considering the effects of additional terms in the Hamiltonian of Eq. 1.[6] For example, a recent experiment[28] realized an open Dicke model with independently tunable rotating and counter-rotating terms. In addition to the superradiant transition, this generalized Dicke model can undergo a lasing instability, which was termed inverted lasing or counter-lasing.[6] This transition is induced by the counter-rotating terms of the Dicke model and is most prominent when these terms are larger than the rotating ones.

The nonequilibrium superradiant transition and the lasing instability have several similarities and differences. Both transitions are of a mean-field type and can be understood in terms of the dynamics of a single degree of freedom. The superradiant transition corresponds to a supercritical pitchfork bifurcation, while the lasing instability corresponds to a Hopf instability. The key difference between these two types of bifurcations is that the former gives rise to two stable solutions, while the latter leads to periodic solutions (limit cycles). Accordingly, in the superradiant phase the cavity field is static (in the frame of the pump field), while it oscillates periodically in the lasing phase.[6]

See also




  This article was adapted from the following source under a CC BY 4.0 license (2020) (reviewer reports): Mor M Roses; Emanuele Dalla Torre (4 September 2020). "Dicke model". PLOS One. 15 (9): e0235197. doi:10.1371/JOURNAL.PONE.0235197. ISSN 1932-6203. PMID 32886669. Wikidata Q98950147.

  1. ^ a b c d Hepp, Klaus; Lieb, Elliott H (1973). "On the superradiant phase transition for molecules in a quantized radiation field: the dicke maser model". Annals of Physics. 76 (2): 360–404. Bibcode:1973AnPhy..76..360H. doi:10.1016/0003-4916(73)90039-0. ISSN 0003-4916.
  2. ^ a b c Dicke, R. H. (1954). "Coherence in Spontaneous Radiation Processes". Physical Review. 93 (1): 99–110. Bibcode:1954PhRv...93...99D. doi:10.1103/PhysRev.93.99. ISSN 0031-899X.
  3. ^ Note that the spin operators are often represented by Pauli matrices  , through the relation  . In some References, the Hamiltonian of the Dicke model is represented in terms of Pauli matrices, rather than spin operators.
  4. ^ Larson, Jonas; Irish, Elinor K (2017). "Some remarks on 'superradiant' phase transitions in light-matter systems". Journal of Physics A: Mathematical and Theoretical. 50 (17): 174002. arXiv:1612.00336. Bibcode:2017JPhA...50q4002L. doi:10.1088/1751-8121/aa65dc. ISSN 1751-8113. S2CID 119474228.
  5. ^ See Garraway, B. M. (2011). "The Dicke model in quantum optics: Dicke model revisited". Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences. 369 (1939): 1137–1155. Bibcode:2011RSPTA.369.1137G. doi:10.1098/rsta.2010.0333. ISSN 1364-503X. PMID 21320910. and references therein.
  6. ^ a b c d See Kirton, Peter; Roses, Mor M.; Keeling, Jonathan; Dalla Torre, Emanuele G. (2018). "Introduction to the Dicke Model: From Equilibrium to Nonequilibrium, and Vice Versa". Advanced Quantum Technologies. 2 (1–2): 1800043. arXiv:1805.09828. doi:10.1002/qute.201800043. hdl:10023/18678. ISSN 2511-9044. S2CID 51695881. and references therein.
  7. ^ de Aguiar, M.A.M; Furuya, K; Lewenkopf, C.H; Nemes, M.C (1992). "Chaos in a spin-boson system: Classical analysis". Annals of Physics. 216 (2): 291–312. Bibcode:1992AnPhy.216..291D. doi:10.1016/0003-4916(92)90178-O.
  8. ^ a b Bastarrachea-Magnani, Miguel Angel; López-del-Carpio, Baldemar; Lerma-Hernández, Sergio; Hirsch, Jorge G (2015). "Chaos in the Dicke model: quantum and semiclassical analysis". Physica Scripta. 90 (6): 068015. Bibcode:2015PhyS...90f8015B. doi:10.1088/0031-8949/90/6/068015. ISSN 0031-8949. S2CID 123306936.
  9. ^ Emary, Clive; Brandes, Tobias (2003). "Quantum Chaos Triggered by Precursors of a Quantum Phase Transition: The Dicke Model". Physical Review Letters. 90 (4): 044101. arXiv:cond-mat/0207290. Bibcode:2003PhRvL..90d4101E. doi:10.1103/PhysRevLett.90.044101. PMID 12570425. S2CID 18617320.
  10. ^ Chávez-Carlos, Jorge; Bastarrachea-Magnani, Miguel Angel; Lerma-Hernández, Sergio; Hirsch, Jorge G. (2016). "Classical chaos in atom-field systems". Physical Review E. 94 (2): 022209. arXiv:1604.00725. Bibcode:2016PhRvE..94b2209C. doi:10.1103/PhysRevE.94.022209. PMID 27627300. S2CID 18123430.
  11. ^ Chávez-Carlos, Jorge; López-del-Carpio, Baldemar; Bastarrachea-Magnani, Miguel Angel; Stránský, Pavel; Lerma-Hernández, Sergio; Santos, Lea F.; Hirsch, Jorge G. (2019). "Quantum and Classical Lyapunov Exponents in Atom-Field Interaction Systems". Physical Review Letters. 122 (2): 024101. arXiv:1807.10292. Bibcode:2019PhRvL.122b4101C. doi:10.1103/PhysRevLett.122.024101. PMID 30720302.
  12. ^ Pilatowsky-Cameo, Saúl; Chávez-Carlos, Jorge; Bastarrachea-Magnani, Miguel A.; Stránský, Pavel; Lerma-Hernández, Sergio; Santos, Lea F.; Hirsch, Jorge G. (2020). "Positive quantum Lyapunov exponents in experimental systems with a regular classical limit". Physical Review E. 101 (1): 010202. arXiv:1909.02578. Bibcode:2020PhRvE.101a0202P. doi:10.1103/PhysRevE.101.010202. PMID 32069677. S2CID 210023416.
  13. ^ Lerma-Hernández, S.; Villaseñor, D.; Bastarrachea-Magnani, M. A.; Torres-Herrera, E. J.; Santos, L. F.; Hirsch, J. G. (2019). "Dynamical signatures of quantum chaos and relaxation time scales in a spin-boson system". Physical Review E. 100 (1): 012218. arXiv:1905.03253. Bibcode:2019PhRvE.100a2218L. doi:10.1103/PhysRevE.100.012218. PMID 31499773.
  14. ^ Villaseñor, David; Pilatowsky-Cameo, Saúl; Bastarrachea-Magnani, Miguel Angel; Lerma-Hernández, Sergio; Santos, Lea F.; Hirsch, Jorge G. (2020). "Quantum vs classical dynamics in a spin-boson system: manifestations of spectral correlations and scarring". New Journal of Physics. 22 (6): 063036. arXiv:2002.02465. Bibcode:2020NJPh...22f3036V. doi:10.1088/1367-2630/ab8ef8. ISSN 1367-2630.
  15. ^ Aguiar, M. A. M. de; Furuya, K; Lewenkopf, C. H; Nemes, M. C (1991). "Particle-Spin Coupling in a Chaotic System: Localization-Delocalization in the Husimi Distributions". Europhysics Letters (EPL). 15 (2): 125–131. Bibcode:1991EL.....15..125D. doi:10.1209/0295-5075/15/2/003. ISSN 0295-5075. S2CID 250898244.
  16. ^ Pilatowsky-Cameo, Saúl; Villaseñor, David; Bastarrachea-Magnani, Miguel A; Lerma, Sergio; Santos, Lea F; Hirsch, Jorge G (2021). "Quantum scarring in a spin-boson system: fundamental families of periodic orbits". New Journal of Physics. 23 (3): 033045. arXiv:2009.08523. Bibcode:2021NJPh...23c3045P. doi:10.1088/1367-2630/abd2e6. ISSN 1367-2630.
  17. ^ Scully, Marlan O.; Zubairy, M. Suhail (1997). Quantum Optics. Cambridge University Press. doi:10.1017/CBO9780511813993. ISBN 9780521435956.
  18. ^ Gross, M.; Haroche, S. (1982). "Superradiance: An essay on the theory of collective spontaneous emission". Physics Reports. 93 (5): 301–396. Bibcode:1982PhR....93..301G. doi:10.1016/0370-1573(82)90102-8. ISSN 0370-1573.
  19. ^ Frisk Kockum, Anton; Miranowicz, Adam; De Liberato, Simone; Savasta, Salvatore; Nori, Franco (2019). "Ultrastrong coupling between light and matter". Nature Reviews Physics. 1 (1): 19–40. arXiv:1807.11636. Bibcode:2019NatRP...1...19F. doi:10.1038/s42254-018-0006-2. ISSN 2522-5820. S2CID 51963566.
  20. ^ Dimer, F.; Estienne, B.; Parkins, A. S.; Carmichael, H. J. (2007). "Proposed realization of the Dicke-model quantum phase transition in an optical cavity QED system". Physical Review A. 75 (1): 013804. arXiv:quant-ph/0607115. Bibcode:2007PhRvA..75a3804D. doi:10.1103/PhysRevA.75.013804. ISSN 1050-2947. S2CID 5513429.
  21. ^ Nagy, D.; Kónya, G.; Szirmai, G.; Domokos, P. (2010). "Dicke-Model Phase Transition in the Quantum Motion of a Bose-Einstein Condensate in an Optical Cavity". Physical Review Letters. 104 (13): 130401. arXiv:0912.3260. Bibcode:2010PhRvL.104m0401N. doi:10.1103/PhysRevLett.104.130401. ISSN 0031-9007. PMID 20481867. S2CID 36968566.
  22. ^ a b Baumann, Kristian; Guerlin, Christine; Brennecke, Ferdinand; Esslinger, Tilman (2010). "Dicke quantum phase transition with a superfluid gas in an optical cavity". Nature. 464 (7293): 1301–1306. arXiv:0912.3261. Bibcode:2010Natur.464.1301B. doi:10.1038/nature09009. ISSN 0028-0836. PMID 20428162. S2CID 205220396.
  23. ^ Black, Adam T.; Chan, Hilton W.; Vuletić, Vladan (2003). "Observation of Collective Friction Forces due to Spatial Self-Organization of Atoms: From Rayleigh to Bragg Scattering". Physical Review Letters. 91 (20): 203001. Bibcode:2003PhRvL..91t3001B. doi:10.1103/PhysRevLett.91.203001. ISSN 0031-9007. PMID 14683358.
  24. ^ Klinder, Jens; Keßler, Hans; Wolke, Matthias; Mathey, Ludwig; Hemmerich, Andreas (2015). "Dynamical phase transition in the open Dicke model". Proceedings of the National Academy of Sciences. 112 (11): 3290–3295. arXiv:1409.1945. Bibcode:2015PNAS..112.3290K. doi:10.1073/pnas.1417132112. ISSN 0027-8424. PMC 4371957. PMID 25733892.
  25. ^ Ritsch, Helmut; Domokos, Peter; Brennecke, Ferdinand; Esslinger, Tilman (2013). "Cold atoms in cavity-generated dynamical optical potentials". Reviews of Modern Physics. 85 (2): 553–601. arXiv:1210.0013. Bibcode:2013RvMP...85..553R. doi:10.1103/RevModPhys.85.553. ISSN 0034-6861. S2CID 118314846.
  26. ^ Mivehvar, Farokh; Piazza, Francesco; Donner, Tobias; Ritsch, Helmut (2 January 2021). "Cavity QED with quantum gases: new paradigms in many-body physics". Advances in Physics. 70 (1): 1–153. arXiv:2102.04473. Bibcode:2021AdPhy..70....1M. doi:10.1080/00018732.2021.1969727. ISSN 0001-8732. S2CID 231855430.
  27. ^ Zhiqiang, Zhang; Lee, Chern Hui; Kumar, Ravi; Arnold, K. J.; Masson, Stuart J.; Parkins, A. S.; Barrett, M. D. (2017). "Nonequilibrium phase transition in a spin-1 Dicke model". Optica. 4 (4): 424. arXiv:1612.06534. Bibcode:2017Optic...4..424Z. doi:10.1364/OPTICA.4.000424. ISSN 2334-2536. S2CID 119232097.
  28. ^ a b Zhang, Zhiqiang; Lee, Chern Hui; Kumar, Ravi; Arnold, K. J.; Masson, Stuart J.; Grimsmo, A. L.; Parkins, A. S.; Barrett, M. D. (2018). "Dicke-model simulation via cavity-assisted Raman transitions". Physical Review A. 97 (4): 043858. arXiv:1801.07888. Bibcode:2018PhRvA..97d3858Z. doi:10.1103/PhysRevA.97.043858. ISSN 2469-9926. S2CID 119232888.