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 ${\displaystyle N}$ two-level systems, or equivalently ${\displaystyle N}$ 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 ${\displaystyle n}$ photons, denoted by ${\displaystyle |n\rangle }$. These states can be constructed from the vacuum state ${\displaystyle |n=0\rangle }$ using the canonical ladder operators, ${\displaystyle a^{\dagger }}$ and ${\displaystyle a}$, 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 ${\displaystyle {\vec {\sigma }}_{j}=(\sigma _{j}^{x},~\sigma _{j}^{y},~\sigma _{j}^{z})}$, satisfying the spin algebra ${\displaystyle [\sigma _{j}^{x},~\sigma _{k}^{y}]=i\hbar \sigma _{j}^{z}\delta _{j,k}}$. Here ${\displaystyle \hbar }$ is the reduced Planck constant and ${\displaystyle j=(0,1,2,...,N)}$ 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 ${\displaystyle \hbar \omega _{c}}$ (where ${\displaystyle \omega _{c}}$ is the cavity frequency), times the number of photons in the cavity, ${\displaystyle n_{c}=a^{\dagger }a}$. The second term describes the energy of the two-level systems, where ${\displaystyle \hbar \omega _{z}}$ 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, ${\displaystyle \lambda }$, 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 ${\displaystyle N\to \infty }$ (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 ${\displaystyle a\sigma ^{+}+a^{\dagger }\sigma ^{-}}$ and a counter-rotating term proportional to ${\displaystyle a\sigma ^{-}+a^{\dagger }\sigma ^{+}}$, where ${\displaystyle \sigma ^{\pm }=\sigma ^{x}\pm i\sigma ^{y}}$ 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 ${\displaystyle S^{\alpha }=\sum _{j=1}^{N}\sigma _{j}^{\alpha }}$, with ${\displaystyle \alpha =x,y,z}$, which satisfy the spin algebra, ${\displaystyle [S^{x},S^{y}]=i\hbar S^{z}}$. Using these operators, one can rewrite the Hamiltonian in Eq. 1 as

$${\displaystyle H=\hbar \omega _{c}a^{\dagger }a+\omega _{z}S^{z}+{\frac {2\lambda }{\sqrt {N}}}(a+a^{\dagger })S^{x}.}$$

(2)

This notation simplifies the numerical study of the model because it involves a single spin-S with ${\displaystyle S\leq N/2}$, whose Hilbert space has size ${\displaystyle 2S+1}$, rather than ${\displaystyle N}$ spin-1/2, whose Hilbert space has size ${\displaystyle 2^{N}}$.

The Dicke model has one global symmetry,

Because ${\displaystyle {\mathcal {P}}}$ squares to unity (i.e. if applied twice, it brings each state back to its original state), it has two eigenvalues, ${\displaystyle 1}$ and ${\displaystyle -1}$. This symmetry is associated with a conserved quantity: the parity of the total number of excitations, ${\displaystyle P=(-1)^{N_{ex}}}$, 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 ${\displaystyle \pm 2}$. 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, ${\displaystyle N=1}$, is called the Rabi model. In the absence of counter-rotating terms, the model is called Jaynes-Cummings for ${\displaystyle N=1}$ and Tavis-Cummings for ${\displaystyle N>1}$. These two models conserve the number of excitations ${\displaystyle N_{ex}}$ and are characterized by a ${\displaystyle U(1)}$ 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	${\displaystyle U(1)}$	${\displaystyle N=1}$
Tavis-Cummings	no	${\displaystyle U(1)}$	${\displaystyle N>1}$
Rabi model	yes	${\displaystyle {\mathcal {P}}}$	${\displaystyle N=1}$
Dicke	yes	${\displaystyle {\mathcal {P}}}$	${\displaystyle N>1}$

Early studies of the Dicke model considered its equilibrium properties.^[1] These works considered the limit of ${\displaystyle N\to \infty }$ (also known as the thermodynamic limit) and assumed a thermal partition function, ${\displaystyle Z=\exp(-H/k_{B}T)}$, where ${\displaystyle k_{B}}$ is the Boltzmann constant and ${\displaystyle T}$ is the temperature. It was found that, when the coupling ${\displaystyle \lambda }$ crosses a critical value ${\displaystyle \lambda _{c}}$, 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, ${\displaystyle {\mathcal {P}}}$, defined in Eq. 3. The order parameter of this phase transition is ${\displaystyle \langle {a}\rangle /{\sqrt {N}}}$. 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 ${\displaystyle x}$ components). Close to the superradiant phase transition, the order parameter depends on ${\displaystyle \lambda }$ as ${\displaystyle \langle {a}\rangle /{\sqrt {N}}\sim (\lambda _{c}-\lambda )^{-1/2}}$. This dependence corresponds to the mean-field critical exponent ${\displaystyle \beta =1/2}$.

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 ${\displaystyle dF/d\alpha (\alpha =0)=0}$, leading to

For ${\displaystyle \lambda <\lambda _{c}}$, ${\displaystyle F}$ has one minimum, while for ${\displaystyle \lambda >\lambda _{c}}$, it has two minima. In the limit of ${\displaystyle T\to 0}$ one obtains an expression for the critical coupling of the zero-temperature superradiant phase transition, ${\displaystyle \lambda _{c}={\sqrt {\omega _{c}\omega _{z}}}/2}$.

A phase space for the Dicke model in the symmetric atomic subspace with ${\displaystyle S=N/2}$ may be constructed by considering the tensor product of the Glauber coherent states

$${\displaystyle \left\vert q,p\right\rangle =D(q,p)\left\vert 0\right\rangle ,}$$

(7)

where ${\textstyle D(q,p)=e^{-{\frac {S}{4}}\left(q^{2}+p^{2}\right)}\exp \left({\sqrt {\frac {S}{2}}}\left(q+ip\right)a^{\dagger }\right)}$ is the displacement operator and ${\textstyle \left\vert 0\right\rangle }$ is the photon vacuum Fock state, and the SU(2) coherent states

$${\displaystyle \left\vert Q,P\right\rangle =R(Q,P)\left\vert S,-S\right\rangle ,}$$

(8)

where ${\textstyle R(Q,P)=\left(1-{\frac {Q^{2}+P^{2}}{4}}\right)^{S}\exp \left({\sqrt {\frac {1}{4-Q^{2}-P^{2}}}}\left(Q+iP\right){\frac {S^{+}}{\hbar }}\right)}$ is the rotation operator in the Bloch sphere, ${\textstyle Q^{2}+P^{2}\leq 4,}$ and ${\displaystyle \left\vert {S,-S}\right\rangle }$ is the state with all atoms in their ground state. This yields a four-dimensional phase space with canonical coordinates ${\textstyle (q,p)}$ and ${\displaystyle (Q,P)}$.

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

$${\displaystyle H_{\text{cl}}(q,p,Q,P)={\frac {\hbar S\omega _{c}}{2}}\left(q^{2}+p^{2}\right)+{\frac {\hbar S\omega _{z}}{2}}\left(Q^{2}+P^{2}\right)+2\hbar S\lambda \,Qq{\sqrt {1-{\frac {Q^{2}+P^{2}}{4}}}}-\hbar S\omega _{z}.}$$

(9)

In the limit of ${\displaystyle N\to \infty }$, 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 ${\displaystyle N}$.

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 ${\textstyle \lambda }$, ${\textstyle \omega _{c}}$, and ${\textstyle \omega _{z}}$ and the energy ${\textstyle E}$.^[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]

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, ${\displaystyle \rho }$ is the density matrix of the system, ${\displaystyle L_{\alpha }}$ is the Lindblad operator of the decay channel ${\displaystyle \alpha }$, and ${\displaystyle \gamma _{\alpha }}$ the associated decay rate. When the Hamiltonian ${\displaystyle H}$ 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
Lindbladian	${\displaystyle L=a}$	${\displaystyle L=\sigma _{i}^{-}}$	${\displaystyle L=\sigma _{i}^{z}}$	${\displaystyle L=\sum _{i}\sigma _{i}^{-}}$
Decay rate	${\displaystyle \kappa }$	${\displaystyle \gamma _{\downarrow }}$	${\displaystyle \gamma _{\phi }}$	${\displaystyle \Gamma _{\downarrow }}$

In the theoretical description of the model, one often considers the steady state where ${\displaystyle d\rho /dt=0}$. In the limit of ${\displaystyle N\to \infty }$, 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).

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, ${\displaystyle N}$, and the maximal intensity of the emitted light scales as ${\displaystyle N^{2}}$. This is in contrast to the spontaneous emission of ${\displaystyle N}$ independent two-level systems, whose decay time does not depend on ${\displaystyle N}$ and where the pulse intensity scales as ${\displaystyle N}$.

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 ${\displaystyle N}$, while in the superradiant phase, the intensity of the cavity field is proportional to ${\displaystyle \langle a^{\dagger }a\rangle \sim N}$.

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)	${\displaystyle N}$	${\displaystyle 1}$
Intensity of the field (superradiant)	${\displaystyle N^{2}}$	${\displaystyle N}$

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 ${\displaystyle \lambda _{c}}$, see Eq. 6.^[19] (2) An accurate modelling of the physical system requires to consider ${\displaystyle A^{2}}$ 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.

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]

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.