---
title: Dicke Lattice Model Overview
url: https://www.emergentmind.com/topics/dicke-lattice-model
type: topic
---

# Dicke Lattice Model Overview

The Dicke lattice model denotes a class of spatial extensions of Dicke superradiance in which collective light–matter coupling is embedded in a lattice or otherwise spatially structured setting. In its standard form, each site of a cavity or resonator array hosts a bosonic mode coupled to a collective spin, while photon hopping couples neighboring sites; related usages include bosonic optical lattices whose density-wave degree of freedom couples to a single cavity mode, as well as momentum-space constructions built from timed Dicke states [2607.00557] [1511.00850] [1403.7097]. Across these variants, the defining feature is that superradiant ordering no longer occurs in a spatially featureless environment: it competes with hopping, detuning, on-site interactions, boundary conditions, and dissipation, producing normal, patterned superradiant, superfluid, and Mott-insulating regimes.

## 1. Definition and range of usages

In the cavity-array formulation, the Dicke lattice model is a spatial extension of the Dicke model to an array of coupled cavity–spin units. Each lattice site \(i\) contains a single bosonic cavity mode \(c_i\) of frequency \(\omega_c\) and an ensemble of \(N_a\) identical two-level atoms or spins described by a collective spin \(S_i\) of frequency \(\omega_a\), with neighboring cavities coupled by photon hopping \(\xi\) [2607.00557]. The competition between local Dicke coupling and inter-site hopping organizes superradiant order into lattice-symmetry sectors distinguished by momentum and sign pattern.

A second usage arises in ultracold atoms in optical cavities. There, a Bose–Einstein condensate loaded into an optical lattice inside a high-finesse cavity realizes what was termed a Dicke–Hubbard setting: a lattice of bosonic sites whose collective density-wave operator couples to a single radiation mode, generating an effective infinite-range interaction across the lattice [1511.00850]. In that setting, the Dicke degree of freedom is not a local two-level atom per site, but a coarse-grained checkerboard density-wave mode embedded in a Bose–Hubbard system.

Further related constructions broaden the term’s scope. The long-range Dicke–Ising model consists of a single global bosonic mode coupled uniformly to a spin lattice with algebraically decaying antiferromagnetic Ising interactions; it realizes a “Dicke lattice” in the sense of a spatially structured spin sector, but not an array of hopping photon modes [2503.02734]. The “superradiance lattice” goes further and realizes a tight-binding lattice in momentum space from timed Dicke states coupled by electromagnetically induced transparency (EIT) fields [1403.7097]. Taken together, these works show that “Dicke lattice model” is a family resemblance term rather than a single Hamiltonian.

## 2. Hamiltonian structure and symmetry

The standard closed Dicke lattice Hamiltonian on a periodic ring is
\[
H=\sum_{i=1}^{N} H_i^{\rm Dicke}-\xi\sum_{i=1}^{N}(c_i^\dagger c_{i+1}+c_{i+1}^\dagger c_i),
\]
with
\[
H_i^{\rm Dicke}=\omega_c c_i^\dagger c_i+\omega_a S_i^z+\frac{2g}{\sqrt{N_a}}(c_i+c_i^\dagger)S_i^x.
\]
Its dissipative extension is governed by the Lindblad equation
\[
\frac{d\rho}{dt}=-i[H,\rho]+\kappa\sum_i D[c_i]\rho,
\qquad
D[O]\rho=2O\rho O^\dagger-O^\dagger O\rho-\rho O^\dagger O,
\]
with cavity loss \(\kappa\) taken as the dominant dissipative channel in the analyses under discussion [2607.00557].

Photon hopping reduces the independent local \(Z_2\) symmetries of uncoupled sites to a single global \(Z_2\) and imposes discrete translation symmetry. Superradiant phases are therefore classified by irreducible representations of the translation group, equivalently by lattice momenta \(k\), and on finite rings by the sitewise sign structure of \(\mathrm{Re}\,\alpha_i\), where \(\alpha_i=\langle c_i\rangle\) [2607.00557]. In Dicke–Hubbard lattices with counter-rotating terms, the relevant conserved symmetry is likewise a global \(Z_2\) parity,
\[
P=\exp\!\left[i\pi\sum_n\left(a_n^\dagger a_n+J_z^n+\frac{N}{2}\right)\right],
\]
which is spontaneously broken when the photonic order parameter \(\psi=\langle a\rangle\) becomes nonzero [1604.06920].

In the cavity-coupled Bose–Hubbard realization, the matter sector is
\[
H_{BH}=-t\sum_{\langle i,j\rangle}(b_i^\dagger b_j+\text{h.c.})+\frac{U}{2}\sum_i n_i(n_i-1)-\mu\sum_i n_i,
\]
while the cavity field and light–matter coupling are
\[
H_{\rm cav}=-\delta_{\rm eff}\,a^\dagger a,
\qquad
H_{\rm int}=g(a+a^\dagger)\Theta,
\]
with
\[
\Theta=\sum_i f_i n_i,\qquad f_i=\cos(ky_i)\cos(kz_i).
\]
After adiabatic elimination of the cavity in the stationary weakly dissipative regime,
\[
\alpha\equiv\langle a\rangle=\frac{g}{\delta_{\rm eff}+i\kappa}\langle\Theta\rangle,
\]
and the atoms experience the effective interaction
\[
H_{\rm eff}=H_{BH}-\frac{g^2\delta_{\rm eff}}{\delta_{\rm eff}^2+\kappa^2}\Theta^2.
\]
For red detuning \(\delta_{\rm eff}<0\), the induced infinite-range term lowers the energy when \(|\langle\Theta\rangle|\) grows, favoring checkerboard density-wave order and superradiance [1511.00850].

## 3. Ordered phases, order parameters, and critical behavior

For dissipative cavity arrays, the normal phase is defined by \(\alpha_i=0\) and \(S_i^z\simeq -N_a/2\), while a superradiant phase has finite \(\alpha_i\) and \(S_i^x\) [2607.00557]. On a four-site ring, the principal superradiant configuration classes are the uniform pattern \([+ + + +]\), the staggered pattern \([+ - + -]\), the stripe-like pattern \([+ + - -]\), and a low-symmetry three-minus pattern \([+ - - -]\), with degeneracies \(m=2,2,4,8\), respectively. The first three were labeled HSRP, ISRP1, and ISRP2; the fourth, obtained numerically, was labeled ISRP3 [2607.00557].

The linear instability of the normal phase in the dissipative model occurs at
\[
g_c^{NP}=\min_k \frac{1}{2}\sqrt{\omega_a\omega_k\left(1+\frac{\kappa^2}{\omega_k^2}\right)},
\]
where \(\omega_k\) is the photon normal-mode frequency. In the closed model, minimizing the mean-field energy yields
\[
g_c^\pm=\frac{1}{2}\sqrt{\omega_a(\omega_c\mp 2\xi)}.
\]
For \(\xi>0\), the equilibrium ground state is uniform; for \(\xi<0\), it is staggered [2607.00557]. The same work further found that, in the dissipative problem, different nonequilibrium branches can belong to different universality classes: NP\(\rightarrow\)HSRP and NP\(\rightarrow\)ISRP1 have \(\gamma_{\rm open}=1\), whereas NP\(\rightarrow\)ISRP2 has \(\gamma_{\rm open}=1/2\). By contrast, the closed Dicke lattice shares the equilibrium Dicke exponent \(\gamma_{\rm closed}=1/2\) across distinct configurations [2607.00557].

In the optical-lattice realization, three phases were observed in the ground-state phase diagram. The homogeneous superfluid has \(\alpha=0\), \(\rho_{DW}=0\), and \(\psi\neq 0\); the self-organized superradiant superfluid has \(\alpha\neq 0\), \(\rho_{DW}\neq 0\), and \(\psi\neq 0\); and the self-organized superradiant Mott insulator has \(\alpha\neq 0\), \(\rho_{DW}\neq 0\), and \(\psi=0\) [1511.00850]. Here \(\psi=\langle b_i\rangle\) measures matter-wave coherence, \(\alpha=\langle a\rangle\) the coherent intracavity field, and \(\rho_{DW}\propto N_a^{-1}\langle\Theta\rangle\) the checkerboard density-wave order. The cavity-induced term \(-\lambda\Theta^2\) with \(\lambda=g^2\delta_{\rm eff}/(\delta_{\rm eff}^2+\kappa^2)<0\) reshapes the Bose–Hubbard phase structure by suppressing number fluctuations on the “wrong” sublattice and stabilizing a Mott insulator inside a self-organized density wave [1511.00850].

## 4. Dissipation, boundaries, and multistability

A defining nonequilibrium feature of the dissipative Dicke lattice is multistability. Photon hopping organizes the candidate superradiant steady states by lattice symmetry, and dissipation can stabilize several symmetry-distinguished configurations simultaneously [2607.00557]. For the dissipative four-site ring, the complete phase diagram contains a stable normal region, regions with only one stable superradiant phase, twofold coexistence regions, threefold coexistence regions, and a fourfold coexistence region in which HSRP, ISRP1, ISRP2, and ISRP3 are all dynamically stable. A representative parameter set exhibiting fourfold coexistence is \(\omega_a=\omega_c\equiv\omega\), \(\kappa=0.4\omega\), \(\xi=0.1\omega\), and \(g=0.8\omega\) [2607.00557]. Which branch is reached depends on initial conditions because several attracting fixed points can coexist.

Open boundary conditions qualitatively alter this picture. For the dissipative Dicke lattice with finite \(N\), open boundaries are implemented by setting the end-to-end hopping \(\lambda=0\), whereas periodic boundaries correspond to \(\lambda=\xi\). Under open boundary conditions, the photonic normal-mode dispersion changes from
\[
\omega_{P,k}=\omega_c-2\xi\cos\!\left[\frac{2\pi(k-1)}{N}\right]
\]
to
\[
\omega_{O,k}=\omega_c-2\xi\cos\!\left[\frac{\pi k}{N+1}\right].
\]
The normal-to-superradiant threshold remains
\[
g_c^{P(O),NP}=\min_k \frac{1}{2}\sqrt{\omega_a\omega_{P(O),k}\left(1+\frac{\kappa^2}{\omega_{P(O),k}^2}\right)},
\]
but the phase structure becomes strongly boundary-sensitive [2508.10296].

For \(N=3\), open boundaries generate a “zoo of superradiant phases” not present in the corresponding infinite system: O1 with \(\langle c_1\rangle=\langle c_3\rangle\neq\langle c_2\rangle\neq0\) and equal signs on all sites, O2 with \(\langle c_1\rangle=-\langle c_3\rangle\) and \(\langle c_2\rangle=0\), O3 with all three amplitudes unequal and nonzero, and O4 with equal edge amplitudes opposite in sign to the center [2508.10296]. The same study showed that a homogeneous superradiant steady state is absent under open boundary conditions for any finite \(N\), because the edge and bulk steady-state constraints cannot be satisfied by one uniform complex \(\langle c_j\rangle\). This produces monostable, bistable, and tristable regions that are absent or greatly simplified under periodic boundaries [2508.10296].

A related dissipative effect is finite-momentum instability in hybrid microwave-cavity arrays. There the mode-resolved threshold is
\[
G_k=\sqrt{\frac{\Delta_k\Delta_s}{4}\left(1+\frac{\kappa^2}{\Delta_k^2}\right)},
\qquad
\Delta_k=\Delta_c-2t\cos k,
\]
so \(G_{\rm crit}=\min_k G_k\). For \(\Delta_c-2t>\kappa\), the instability is at \(k=0\), giving a homogeneous superradiant transition. For \(0<\Delta_c-2t<\kappa<\Delta_c+2t\), the transition occurs at finite momentum, with
\[
G_{\rm crit}=\sqrt{\frac{\kappa\Delta_s}{2}},
\qquad
k_c=\arccos\!\left(\frac{\Delta_c-\kappa}{2t}\right).
\]
This finite-\(k\) selection is a distinct non-equilibrium signature of the driven open Dicke lattice model [1405.3289].

## 5. Implementations and direct observations

A direct realization of Dicke–Hubbard physics was achieved with a cigar-shaped \(^{87}\mathrm{Rb}\) Bose–Einstein condensate of \(N_a\approx 5\times 10^4\) atoms in \(|F=2,m_F=2\rangle\), overlapped with a high-finesse cavity mode along \(z\) and an optical lattice inside the cavity [1511.00850]. The cavity had waist \(\approx 32\,\mu{\rm m}\), finesse \(344{,}000\), Purcell factor \(44\), and \(\kappa=2\pi\times 4.45\,{\rm kHz}\). Two standing waves at \(\lambda=803\,{\rm nm}\) were used: a transverse pump along \(y\) and an external lattice along \(x\) at fixed depth \(14\,E_{\rm rec}\), with
\[
E_{\rm rec}=\frac{\hbar^2k^2}{2m}=2\pi\hbar\times 3.56\,{\rm kHz}.
\]
For \(N_a\approx 5\times 10^4\), the dispersive shift was \(\delta_-\approx -2\pi\times 9\,{\rm kHz}\approx -2\kappa\), placing the system in a strong cooperative-coupling regime [1511.00850].

Superradiance was detected through the intracavity photon number \(N_p=\langle a^\dagger a\rangle\) obtained from leaked photons, while matter-wave coherence and density-wave order were monitored by \(25\,{\rm ms}\) time-of-flight absorption images. Bragg peaks at \((\pm2,0)\hbar k\) arose from the pump lattice, and peaks at \((\pm1,\pm1)\hbar k\) from the cavity–pump interference lattice [1511.00850]. The HSF\(\rightarrow\)SSF transition appeared as the onset of nonzero \(N_p\) in the negative-\(\delta_{\rm eff}\) half-plane. The SSF\(\rightarrow\)SMI transition was identified by a sudden increase in the zero-momentum width \(\mathcal W\), together with continued growth of \(N_p\) and kinks in \(N_p\) versus \(\delta_{\rm eff}\) that were interpreted as reduced number fluctuations [1511.00850].

A complementary implementation was proposed for arrays of superconducting microwave cavities coupled to ensembles of nitrogen-vacancy centers in diamond. By cavity-assisted Raman transitions between two NV spin states, the scheme realizes a generalized Dicke interaction containing both Jaynes–Cummings and anti–Jaynes–Cummings terms, with tunable effective detunings and couplings [1405.3289]. The effective lattice Hamiltonian is
\[
H_{\rm DLM}
=
\Delta_c\sum_{\ell=1}^{N_L} a_\ell^\dagger a_\ell
-t\sum_{\ell=1}^{N_L-1}(a_\ell^\dagger a_{\ell+1}+a_\ell a_{\ell+1}^\dagger)
+\sum_{\ell=1}^{N_L}\Delta_s J_\ell^z
+\sum_{\ell=1}^{N_L}\frac{G}{\sqrt{\mathcal N}}(J_\ell^++J_\ell^-)(a_\ell+a_\ell^\dagger),
\]
with photon loss at rate \(\kappa\) [1405.3289]. The paper quoted realistic values \(g_0\sim 10\,{\rm Hz}\), bare collective coupling \(G_0\approx 7.5\)–\(10\,{\rm MHz}\), Raman coupling \(G\approx 1.5\,{\rm MHz}\) for \(\delta_B=100\,{\rm MHz}\) and \(\Omega\) up to \(20\,{\rm MHz}\), cavity decay \(\kappa\sim 0.1\)–\(0.5\,{\rm MHz}\), inhomogeneous broadening \(\gamma_s\sim 20\,{\rm MHz}\), and arrays of \(N_L\sim 50\)–\(100\) cavities as feasible [1405.3289].

Finite cavity or resonator arrays with engineered boundaries are also experimentally relevant in circuit QED and quantum optics. Proposed observables include site-resolved cavity fields \(\alpha_i=\langle c_i\rangle\), photon numbers \(n_i\), spin polarizations \(m_i^{x,z}\), emitted-light phase patterns, and multistability under parameter sweeps [2508.10296].

## 6. Related models, conceptual extensions, and open directions

The Dicke lattice idea has been extended in several directions that change the balance between locality, interaction range, and the role of counter-rotating terms. In the Dicke–Hubbard lattice with full counter-rotating interaction,
\[
H=\sum_n\left[\omega_c a_n^\dagger a_n+\omega_q J_z^n+\frac{2g}{\sqrt N}(a_n+a_n^\dagger)J_x^n\right]
-\kappa\sum_{\langle n,m\rangle}(a_n^\dagger a_m+a_m^\dagger a_n),
\]
mean-field plus extended coherent-state calculations found a localization–delocalization transition marked by a finite photonic order parameter \(\psi=\langle a\rangle\), spontaneous parity breaking, complete suppression of Mott lobes, and monotonic enhancement of \(\psi\) with increasing \(g\) [1604.06920]. In the large-\(N\) limit, the critical line is
\[
\lambda_c=\frac{1}{2}\sqrt{(\omega_c-z\kappa)\omega_q},
\]
and the model contrasts sharply with the rotating-wave Dicke–Hubbard or Tavis–Cummings lattice, where lobe physics survives [1604.06920].

A different extension couples a single bosonic mode to a long-range Ising lattice. In the long-range Dicke–Ising model on square and triangular lattices, the \(g=0\) limit exhibits devil’s staircase magnetization plateaux, while finite light–matter coupling produces both uniform and magnetically ordered superradiant phases with finite photon density [2503.02734]. Examples include a three-sublattice \(1/2\)–\(1/4\)–\(1/4\) superradiant phase on the square lattice and a superradiant Wigner crystal with a four-site unit cell, the \(3/4\)–\(1/4\) SR phase, on the triangular lattice [2503.02734]. There, the transition from normal to superradiant phases is second order with Dicke universality when magnetic order is preserved and first order when it changes [2503.02734].

The superradiance lattice is conceptually more distant but clarifies how Dicke collectivity can itself define a lattice. Timed Dicke states of three-level atoms form momentum-space sites, while a standing-wave EIT coupling induces nearest-neighbor hopping between them. In one dimension, detuning between the two standing-wave components produces an effective uniform force in momentum space, enabling Bloch oscillations, Wannier–Stark ladders, Bloch band collapse, and dynamic localization [1403.7097]. In two dimensions, three coupling beams generate a honeycomb superradiance lattice with graphene-like Dirac physics [1403.7097]. This construction differs from the conventional real-space Dicke lattice of coupled cavities, but it preserves the central idea that collective light–matter states can be arranged into lattice-like structures with tunable connectivity [1403.7097].

Several open directions follow directly from the surveyed literature. Multimode cavities could realize richer Dicke lattices with competing patterns and frustrated long-range interactions; higher dimensions and different geometries may host additional superradiant insulating and supersolid phases; and theoretical treatments incorporating full cavity dynamics, such as Keldysh or Langevin approaches, are expected to refine the description of hysteresis and bistability beyond static cavity elimination [1511.00850]. In finite arrays, boundary engineering already emerges as a control parameter of the same status as hopping and dissipation [2508.10296]. These developments indicate that the Dicke lattice model is best understood not as one fixed Hamiltonian, but as a broad framework for studying how collective radiation reorganizes many-body order once spatial structure is made dynamical.

Source: https://www.emergentmind.com/topics/dicke-lattice-model