Dicke model

Last updated

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.

Contents

Description

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 nameCounter-rotating terms?symmetryNumber of two-level systems
Jaynes-Cummings no
Tavis-Cummingsno
Rabi modelyes
Dickeyes

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. Superradiant transition.gif
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
E
/
S
{\textstyle E/S}
and the coupling parameter
l
{\textstyle \lambda }
(divided by the critical coupling
l
c
=
o
c
o
z
/
2
{\textstyle \lambda _{c}={\sqrt {\omega _{c}\omega _{z}}}/2}
). The parameters are
o
c
=
o
z
=
1
{\displaystyle \omega _{c}=\omega _{z}=1}
. Percentage of classical chaos in the Dicke Model.svg
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 superradiant 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 decayAtomic decayAtomic dephasingCollective decay
Lindbladian
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. Superradiance disambiguation.gif
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]
EnvironmentFree spaceCavity
DurationTransientSteady 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. Experimental realization of the Dicke model.gif
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

Related Research Articles

<span class="mw-page-title-main">Loop quantum gravity</span> Theory of quantum gravity, merging quantum mechanics and general relativity

Loop quantum gravity (LQG) is a theory of quantum gravity that incorporates matter of the Standard Model into the framework established for the intrinsic quantum gravity case. It is an attempt to develop a quantum theory of gravity based directly on Albert Einstein's geometric formulation rather than the treatment of gravity as a mysterious mechanism (force). As a theory, LQG postulates that the structure of space and time is composed of finite loops woven into an extremely fine fabric or network. These networks of loops are called spin networks. The evolution of a spin network, or spin foam, has a scale on the order of a Planck length, approximately 10−35 meters, and smaller scales are meaningless. Consequently, not just matter, but space itself, prefers an atomic structure.

The Ising model, named after the physicists Ernst Ising and Wilhelm Lenz, is a mathematical model of ferromagnetism in statistical mechanics. The model consists of discrete variables that represent magnetic dipole moments of atomic "spins" that can be in one of two states. The spins are arranged in a graph, usually a lattice, allowing each spin to interact with its neighbors. Neighboring spins that agree have a lower energy than those that disagree; the system tends to the lowest energy but heat disturbs this tendency, thus creating the possibility of different structural phases. The model allows the identification of phase transitions as a simplified model of reality. The two-dimensional square-lattice Ising model is one of the simplest statistical models to show a phase transition.

<span class="mw-page-title-main">Nonlinear Schrödinger equation</span> Nonlinear form of the Schrödinger equation

In theoretical physics, the (one-dimensional) nonlinear Schrödinger equation (NLSE) is a nonlinear variation of the Schrödinger equation. It is a classical field equation whose principal applications are to the propagation of light in nonlinear optical fibers and planar waveguides and to Bose–Einstein condensates confined to highly anisotropic, cigar-shaped traps, in the mean-field regime. Additionally, the equation appears in the studies of small-amplitude gravity waves on the surface of deep inviscid (zero-viscosity) water; the Langmuir waves in hot plasmas; the propagation of plane-diffracted wave beams in the focusing regions of the ionosphere; the propagation of Davydov's alpha-helix solitons, which are responsible for energy transport along molecular chains; and many others. More generally, the NLSE appears as one of universal equations that describe the evolution of slowly varying packets of quasi-monochromatic waves in weakly nonlinear media that have dispersion. Unlike the linear Schrödinger equation, the NLSE never describes the time evolution of a quantum state. The 1D NLSE is an example of an integrable model.

In probability theory and mathematical physics, a random matrix is a matrix-valued random variable—that is, a matrix in which some or all of its entries are sampled randomly from a probability distribution. Random matrix theory (RMT) is the study of properties of random matrices, often as they become large. RMT provides techniques like mean-field theory, diagrammatic methods, the cavity method, or the replica method to compute quantities like traces, spectral densities, or scalar products between eigenvectors. Many physical phenomena, such as the spectrum of nuclei of heavy atoms, the thermal conductivity of a lattice, or the emergence of quantum chaos, can be modeled mathematically as problems concerning large, random matrices.

Superconducting quantum computing is a branch of solid state physics and quantum computing that implements superconducting electronic circuits using superconducting qubits as artificial atoms, or quantum dots. For superconducting qubits, the two logic states are the ground state and the excited state, denoted respectively. Research in superconducting quantum computing is conducted by companies such as Google, IBM, IMEC, BBN Technologies, Rigetti, and Intel. Many recently developed QPUs use superconducting architecture.

The quantum Heisenberg model, developed by Werner Heisenberg, is a statistical mechanical model used in the study of critical points and phase transitions of magnetic systems, in which the spins of the magnetic systems are treated quantum mechanically. It is related to the prototypical Ising model, where at each site of a lattice, a spin represents a microscopic magnetic dipole to which the magnetic moment is either up or down. Except the coupling between magnetic dipole moments, there is also a multipolar version of Heisenberg model called the multipolar exchange interaction.

<span class="mw-page-title-main">Jaynes–Cummings model</span> Model in quantum optics

In quantum optics, the Jaynes–Cummings model is a theoretical model that describes the system of a two-level atom interacting with a quantized mode of an optical cavity, with or without the presence of light. It was originally developed to study the interaction of atoms with the quantized electromagnetic field in order to investigate the phenomena of spontaneous emission and absorption of photons in a cavity. It is named after Edwin Thompson Jaynes and Fred Cummings in the 1960s and was confirmed experimentally in 1987.

In quantum mechanics, a weak value is a quantity related to a shift of a measuring device's pointer when usually there is pre- and postselection. It should not be confused with a weak measurement, which is often defined in conjunction. The weak value was first defined by Yakir Aharonov, David Albert, and Lev Vaidman, published in Physical Review Letters 1988, and is related to the two-state vector formalism. There is also a way to obtain weak values without postselection.

A quantum bus is a device which can be used to store or transfer information between independent qubits in a quantum computer, or combine two qubits into a superposition. It is the quantum analog of a classical bus.

<span class="mw-page-title-main">Landau–Zener formula</span> Formula for the probability that a system will change between two energy states.

The Landau–Zener formula is an analytic solution to the equations of motion governing the transition dynamics of a two-state quantum system, with a time-dependent Hamiltonian varying such that the energy separation of the two states is a linear function of time. The formula, giving the probability of a diabatic transition between the two energy states, was published separately by Lev Landau, Clarence Zener, Ernst Stueckelberg, and Ettore Majorana, in 1932.

The Ghirardi–Rimini–Weber theory (GRW) is a spontaneous collapse theory in quantum mechanics, proposed in 1986 by Giancarlo Ghirardi, Alberto Rimini, and Tullio Weber.

In computational complexity theory, QMA, which stands for Quantum Merlin Arthur, is the set of languages for which, when a string is in the language, there is a polynomial-size quantum proof that convinces a polynomial time quantum verifier of this fact with high probability. Moreover, when the string is not in the language, every polynomial-size quantum state is rejected by the verifier with high probability.

Circuit quantum electrodynamics provides a means of studying the fundamental interaction between light and matter. As in the field of cavity quantum electrodynamics, a single photon within a single mode cavity coherently couples to a quantum object (atom). In contrast to cavity QED, the photon is stored in a one-dimensional on-chip resonator and the quantum object is no natural atom but an artificial one. These artificial atoms usually are mesoscopic devices which exhibit an atom-like energy spectrum. The field of circuit QED is a prominent example for quantum information processing and a promising candidate for future quantum computation.

<span class="mw-page-title-main">Jaynes–Cummings–Hubbard model</span> Model in quantum optics

The Jaynes–Cummings–Hubbard (JCH) model is a many-body quantum system modeling the quantum phase transition of light. As the name suggests, the Jaynes–Cummings–Hubbard model is a variant on the Jaynes–Cummings model; a one-dimensional JCH model consists of a chain of N coupled single-mode cavities, each with a two-level atom. Unlike in the competing Bose–Hubbard model, Jaynes–Cummings–Hubbard dynamics depend on photonic and atomic degrees of freedom and hence require strong-coupling theory for treatment. One method for realizing an experimental model of the system uses circularly-linked superconducting qubits.

<span class="mw-page-title-main">Superradiant phase transition</span> Process in quantum optics

In quantum optics, a superradiant phase transition is a phase transition that occurs in a collection of fluorescent emitters, between a state containing few electromagnetic excitations and a superradiant state with many electromagnetic excitations trapped inside the emitters. The superradiant state is made thermodynamically favorable by having strong, coherent interactions between the emitters.

In quantum mechanics, weak measurement is a type of quantum measurement that results in an observer obtaining very little information about the system on average, but also disturbs the state very little. From Busch's theorem any quantum system is necessarily disturbed by measurement, but the amount of disturbance is described by a parameter called the measurement strength.

In quantum computing, Mølmer–Sørensen gate scheme refers to an implementation procedure for various multi-qubit quantum logic gates used mostly in trapped ion quantum computing. This procedure is based on the original proposition by Klaus Mølmer and Anders Sørensen in 1999-2000.

The quantum Fisher information is a central quantity in quantum metrology and is the quantum analogue of the classical Fisher information. It is one of the central quantities used to qualify the utility of an input state, especially in Mach–Zehnder interferometer-based phase or parameter estimation. It is shown that the quantum Fisher information can also be a sensitive probe of a quantum phase transition. The quantum Fisher information of a state with respect to the observable is defined as

A quantum battery is a type of electric battery that uses the principles of quantum mechanics to store energy. They have the potential to be more efficient and powerful than traditional batteries.

In quantum optics, fhe Tavis–Cummings model is a theoretical model to describe an ensemble of identical two-level atoms coupled symmetrically to a single-mode quantized bosonic field. The model extends the Jaynes–Cummings model to larger spin numbers that represent collections of multiple atoms. It differs from the Dicke model in its use of the rotating-wave approximation to conserve the number of excitations of the system.

References

Open Access logo PLoS transparent.svg 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. 1 2 3 4 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. 1 2 3 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. 1 2 3 4 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. 1 2 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. 1 2 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. 1 2 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.