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

# Open Dicke Lattice Model

Searching arXiv for recent papers on the open Dicke lattice model and closely related open Dicke criticality.
The **open Dicke lattice model** denotes a class of driven or dissipative many-body light–matter systems in which cavity or resonator modes are arranged on a lattice, each site is locally coupled to a collective spin or atomic ensemble by a Dicke-type interaction, and openness enters through environmental coupling such as cavity loss. In the most direct lattice realizations, photon hopping endows the photonic sector with a dispersion or band structure, while the on-site light–matter coupling retains the full Dicke form with rotating and counter-rotating contributions [1405.3289][2508.10296][2607.00557]. The topic sits at the intersection of nonequilibrium quantum optics, cavity QED, circuit QED, and dissipative many-body physics, and it is distinguished from the single-mode open Dicke model by the presence of spatial structure, boundary conditions, and momentum- or configuration-selective instabilities [1405.3289][2508.10296].

## 1. Definition and model classes

A standard open Dicke lattice consists of an array of coupled resonators, each locally coupled to a collective spin degree of freedom. One explicit formulation is the Dicke lattice Hamiltonian
\[
H_{\rm DLM}=\Delta_c \sum_{\ell=1}^{N_L} a^\dagger_{\ell} a_{\ell} -t\sum_{\ell=1}^{N_L-1}(a^\dagger_\ell a_{\ell+1}+a_\ell a^\dag_{\ell+1})
+ \sum_{\ell=1}^{N_L} \Delta_s  J^z_\ell + \sum_{\ell=1}^{N_L} \frac{G}{\sqrt{\mathcal N}}(J^+_{\ell}+J^-_\ell)(a_\ell+a_\ell^\dag)
\]
with local cavity operators \(a_\ell\), collective spin operators \(J_\ell^{\pm,z}\), hopping \(t\), cavity detuning \(\Delta_c\), spin detuning \(\Delta_s\), and collective coupling \(G\) [1405.3289]. A closely related finite-size formulation uses
\[
H= \sum_{i=1}^{N} H_i^{\rm Dicke} - \sum_{i=1}^{N-1} \xi \left( c_{i}^{\dagger} c_{i+1} + c_{i} c_{i+1}^{\dagger} \right) - \lambda \left( c_{1}^{\dagger} c_{N} + c_{1} c_{N}^{\dagger} \right),
\]
with local Dicke blocks
\[
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
\]
and boundary parameter \(\lambda=\xi\) for periodic boundary conditions and \(\lambda=0\) for open boundary conditions [2508.10296].

A periodic one-dimensional ring formulation is
\[
H= \sum_{i=1}^{N} H_i^{\rm Dicke} - \xi \sum_{i=1}^{N} \left( c_i^\dagger c_{i+1} + c_{i+1}^\dagger c_i \right),
\]
with \(c_{N+1}=c_1\) and
\[
H_i^{\rm Dicke} = \omega_c c_i^\dagger c_i + \omega_a S_i^z + \frac{2g}{\sqrt{N_a}} \left( c_i+c_i^\dagger \right) S_i^x .
\]
This form underlies a configuration-based analysis of nonequilibrium superradiant phases in finite rings [2607.00557].

These models are “open” because the dynamics is not unitary. In the most common treatments, photon loss is included through a Lindblad master equation,
\[
\frac{d\rho}{dt} = -i[H,\rho] + \kappa \sum_i \left(2 c_i \rho c_i^\dagger - c_i^\dagger c_i \rho - \rho c_i^\dagger c_i\right),
\]
or its site-indexed variants [2508.10296][2607.00557][1405.3289]. By contrast, several influential works on single-mode open Dicke systems treat non-Markovian baths or additional matter dissipation; these are not lattice models, but they supply mechanisms that can be imported into lattice settings [1611.09378][1503.04672].

## 2. Canonical ingredients: hopping, symmetry, and photonic structure

The defining lattice ingredient is photon hopping. In the superconducting-cavity/NV-center implementation, hopping produces the photon band
\[
\Delta_k = \Delta_c - 2t\cos k,
\]
and after Holstein–Primakoff reduction the quadratic normal-phase Hamiltonian becomes
\[
H_{\rm DLM}= \sum_k\Delta_s b_k^\dagger b_k+\sum_k\Delta_k a^\dagger_k a_k +G\sum_k\big( b^\dagger_k a_k+b_k a^\dagger_k + b_k^\dagger a^\dagger_{-k}+b_k a_{-k} \big)
\]
[1405.3289]. In finite lattices, the photonic normal-mode frequencies depend on boundary conditions:
\[
\omega_{P,k} = \omega_c - 2\xi \cos\left[ \frac{2\pi(k-1)}{N} \right], 
\]
\[
\omega_{O,k} = \omega_c - 2\xi \cos\left[ \frac{\pi k}{N+1} \right], \qquad k=1,2,\dots,N,
\]
with positivity condition \(\omega_{P(O),k}>0\) [2508.10296]. Periodic lattices therefore support plane-wave modes, whereas open chains support standing-wave modes [2508.10296].

The Dicke coupling preserves a discrete parity structure. In the single-site model discussed in the non-Markovian setting, the model has a \(Z_2\) parity symmetry \((a,S^x)\to(-a,-S^x)\) [1611.09378]. In the lattice setting, this local Dicke symmetry coexists with lattice translation or boundary symmetries, and the ordered phases are best classified by the spatial pattern of the local cavity order parameters. In the four-site ring, the real part of the local cavity field,
\[
R_j = \frac{\mathrm{Re}\,\langle c_j\rangle}{\sqrt{N_a}},
\]
serves as the local order parameter, and the possible superradiant configurations are organized into sign-pattern classes such as \([+\,+\,+\,+]\), \([+\,-\,+\,-]\), \([+\,+\,-\,-]\), and \([+\,-\,-\,-]\) [2607.00557]. This classification reflects the reduction of independent local parities at \(\xi=0\) to a global \(\mathbb Z_2\) parity once hopping is switched on [2607.00557].

## 3. Dissipation, steady states, and superradiant thresholds

Open Dicke lattice criticality is formulated as a stability problem for nonequilibrium steady states rather than as energy minimization. In the homogeneous single-cavity open Dicke model with cavity loss, the dissipative threshold is shifted from the closed-system value to
\[
G_{\rm crit} = \sqrt{ \frac{\Delta_c\Delta_s}{4}\left(1+\frac{\kappa^2}{\Delta_c^2}\right)}
\]
[1405.3289]. In the lattice, each momentum mode \(k\) has its own threshold,
\[
G_k=\sqrt{\frac{\Delta_k\Delta_s}{4}\left(1+\frac{\kappa^2}{\Delta_k^2}\right)},
\]
and the actual instability is determined by
\[
G_{\rm crit}=\min_k G_k
\]
[1405.3289]. In finite lattices with site-resolved mean-field equations, the normal-phase threshold is
\[
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) }
\]
[2508.10296], while for the four-site periodic ring the normal phase loses stability at
\[
g_c^{\rm NP} = \min_k \frac{1}{2} \sqrt{ \omega_a \omega_k \left( 1+\frac{\kappa^2}{\omega_k^2} \right) } .
\]
The isolated-site threshold remains
\[
g_c=\frac{1}{2}\sqrt{\omega_a\omega_c\left(1+\frac{\kappa^2}{\omega_c^2}\right)}
\]
when \(\xi=0\) [2607.00557].

The normal phase is characterized by vanishing cavity fields and fully polarized spins. In the three-site boundary-sensitive model,
\[
\langle c_j \rangle = \langle S_j^- \rangle = 0,\qquad \langle S_j^z\rangle=-\frac{N_a}{2}
\]
[2508.10296]. In the four-site ring,
\[
R_j=I_j=R_j^s=I_j^s=0,\qquad S_j^z=-\frac12
\]
[2607.00557]. The superradiant phase is identified by nonzero coherent cavity amplitudes on at least one site, \(\langle c_j\rangle\neq 0\) [2508.10296].

A central nonequilibrium point is that the instability is controlled by the complex spectrum of damped modes. In the cavity-array implementation, the transition occurs when one eigenvalue of the linearized dynamical matrix acquires positive real part [1405.3289]. This differs from closed equilibrium reasoning based only on softening of a Hermitian excitation. A closely related lesson emerges in the experimental single-mode open Dicke model: the phase boundary is not determined solely by the vanishing of the excitation frequency, because one must also track the damping or growth rate \(\gamma\) [2502.12155]. This suggests that lattice critical manifolds should generally be located by the full complex spectrum of momentum-resolved modes rather than by real-frequency softening alone.

## 4. Spatial ordering, finite-momentum phases, multistability, and boundaries

One of the defining features of the open Dicke lattice is the possibility of ordered phases beyond homogeneous superradiance. In the cavity-array model, if \(\Delta_c-2t>\kappa\), the first unstable mode is at \(k=0\), giving a homogeneous superradiant phase with \(\langle a_\ell\rangle=\alpha\neq 0\) [1405.3289]. In the parameter window \(0<\Delta_c-2t<\kappa<\Delta_c+2t\), however, the instability occurs at a finite wavevector \(k_c\) satisfying
\[
\Delta_{k_c}=\kappa, \qquad k_c=\arccos\!\left(\frac{\Delta_c-\kappa}{2t}\right),
\]
with threshold
\[
G_{\rm crit}=\sqrt{\frac{\kappa\Delta_s}{2}} .
\]
The corresponding ordered state is spatially modulated,
\[
\langle a_\ell\rangle \simeq \alpha \cos(\phi_0+k_c\ell),
\]
and spontaneously breaks translation symmetry [1405.3289]. For larger losses \(\kappa>\Delta_c+2t\), the preferred instability is at \(k_c=\pi\), corresponding to an antiferromagnetic or staggered pattern [1405.3289].

Finite systems display an additional layer of structure. For the dissipative four-site ring, the complete phase diagram contains a normal phase and four superradiant phase types: a homogeneous superradiant phase and three inhomogeneous superradiant phases [2607.00557]. The homogeneous phase corresponds to
\[
R_1^{\rm ss}=R_2^{\rm ss}=R_3^{\rm ss}=R_4^{\rm ss}\equiv \mathcal R_1,
\]
the alternating phase to
\[
R_1^{\rm ss}=-R_2^{\rm ss}=R_3^{\rm ss}=-R_4^{\rm ss}\equiv \mathcal R_2,
\]
and the domain-wall-like phase to
\[
R_1^{\rm ss}=R_2^{\rm ss}=-R_3^{\rm ss}=-R_4^{\rm ss}\equiv \mathcal R_3
\]
[2607.00557]. The paper further identifies coexistence regions with up to four stable superradiant phases, so multistability is interpreted as simultaneous dynamical stabilization of several symmetry-allowed configuration classes rather than as unstructured nonlinear complexity [2607.00557].

Open boundary conditions alter the phase structure even more dramatically. For a three-site dissipative Dicke lattice, open boundaries support a “zoo of superradiant phases with broken translational symmetry,” including monostable, bistable, and tristable regions [2508.10296]. Representative open-boundary patterns are
\[
\langle c_1 \rangle = \langle c_3 \rangle \neq \langle c_2 \rangle \neq 0,
\]
\[
\langle c_1 \rangle = -\langle c_3 \rangle,\qquad \langle c_2 \rangle = 0,
\]
and
\[
\langle c_1 \rangle \neq \langle c_2 \rangle \neq \langle c_3 \rangle \neq 0
\]
[2508.10296]. A central analytic result is that, for any finite \(N\), a spatially homogeneous steady-state solution is impossible under open boundary conditions because edge and bulk sites satisfy incompatible balance equations:
\[
-\kappa\,\mathrm{Re}(\langle c_j \rangle) + (\omega_c - \xi)\,\mathrm{Im}(\langle c_j \rangle) = 0 \quad \text{for } j=1,N,
\]
whereas
\[
-\kappa\,\mathrm{Re}(\langle c_j \rangle) + (\omega_c - 2\xi)\,\mathrm{Im}(\langle c_j \rangle) = 0 \quad \text{for } j=2,\ldots,N-1.
\]
This incompatibility forbids a fully homogeneous superradiant steady state under open boundary conditions for finite \(N\) [2508.10296].

These results establish that the open Dicke lattice is not simply the periodic Dicke lattice with edges attached. Boundaries, standing-wave mode structure, and dissipative selection reshape the available stationary phases [2508.10296].

## 5. Universality, fluctuations, and the role of dissipation

Dissipation in Dicke systems does not merely shift thresholds; it can also alter fluctuation criticality and universality. In the open four-site ring, discrete truncated Wigner calculations show that different superradiant configurations can belong to different nonequilibrium universality classes. The photon fluctuation scaling
\[
\Delta_j^\beta \propto |g-g_c^\beta|^{-\gamma_{\rm open}}
\]
yields \(\gamma_{\rm open}=-1\) for the normal-to-homogeneous and normal-to-ISRP1 transitions, but \(\gamma_{\rm open}=-1/2\) for the normal-to-ISRP2 transition [2607.00557]. By contrast, in the corresponding closed lattice the lowest excitation gap scales with exponent \(\gamma_{\rm closed}=1/2\), and the equilibrium universality class is shared by different superradiant configurations [2607.00557].

Single-mode open Dicke studies clarify the microscopic reasons dissipation can matter so strongly. In the Markovian open Dicke model, the soft mode is complex, and criticality is governed by the vanishing damping gap rather than only by the real-frequency softening [1206.5131][1107.4323]. In the driven open Dicke model with a sub-Ohmic matter bath, the critical photon-number exponent is \(1\) without spin-bath coupling and becomes less than \(1\) when the spin-like mode couples to a sub-Ohmic reservoir, decreasing monotonically as the bath exponent \(s\) is reduced below \(1\) [1503.04672]. In a distinct non-Markovian construction, coupling the cavity displacement to a zero-temperature power-law bath yields a retarded hybridization function
\[
\Delta^R(\omega)=\Lambda(\omega)-iJ(\omega),
\]
a bath-renormalized photon frequency
\[
\tilde{\omega}_0(\kappa)=\omega_0-\frac{2\kappa}{\pi s},
\]
and a large-\(N\) critical coupling
\[
\lambda^{\infty}_c(\kappa)=\sqrt{\tilde{\omega}_0(\kappa)\omega_q}/2,
\]
so that dissipation can promote rather than suppress superradiance [1611.09378]. These are single-mode results, not lattice results, but they isolate mechanisms—frequency-dependent self-energies, bath-controlled infrared behavior, and dissipative universality changes—that are directly relevant when generalized to momentum-dependent lattice propagators [1611.09378][1503.04672].

A further single-mode result of broad relevance is that local Markovian spin decay can enhance collective quantum correlations in the open Dicke model. In the large-\(N\) fluctuation treatment with cavity loss \(\kappa\) and local spin decay \(\gamma\), the critical coupling is
\[
\lambda_c=\sqrt{\frac{\left[\omega_m^2+\left(\frac{\kappa}{2}\right)^2\right]\left[\omega_z^2+\left(\frac{\gamma}{2}\right)^2\right]}{4 \omega_z \omega_m} ,
\]
and local dissipation can lead to enhancement of logarithmic negativity near the critical line and to a superradiant phase with nonzero spin-boson entanglement [2110.13191]. This suggests that local dissipation in lattice settings may also modify fluctuation structure in ways that are not reducible to simple decoherence.

## 6. Implementations, methods, and scope of the field

A concrete implementation of the Dicke lattice was proposed in arrays of superconducting microwave cavities coupled to NV-center ensembles in diamond. The microscopic model uses cavity-assisted Raman transitions to engineer an effective Dicke interaction, with cavity loss providing openness [1405.3289]. The starting Hamiltonian contains a cavity mode, NV-center ground-state triplet levels \(|m_s=0\rangle, |m_s=\pm1\rangle\), and two microwave drives; after adiabatic elimination of the intermediate state, the effective two-level model becomes
\[
H_{\rm eff}= \Delta_c a^\dag a +  \sum_i  \left[  \frac{\Delta_s^i}{2}+ \lambda_i  a^\dag a \right] (\sigma_z^i+1)
 + \sum_i  \left(g_1^i  a  +  g_2^i a^\dag \right)\sigma_-^i + {\rm H.c.}
\]
[1405.3289]. In the homogeneous limit \(g_1^i=g_2^i=g_i\), this reduces to an effective Dicke form and can be extended to a cavity array with hopping [1405.3289]. The same work emphasizes robustness of the superradiant transition to substantial inhomogeneous broadening, encapsulated in the generalized instability condition
\[
\lim_{\epsilon\rightarrow 0}  \sum_\mu  \frac{4 G_\mu^2 \Delta_c \Delta_\mu }{(\Delta_c^2+\kappa^2)(\Delta_\mu^2+\epsilon^2)}=1
\]
and, for a Lorentzian spin distribution, the broadened critical coupling
\[
G_{\rm crit} = \sqrt{\frac{\Delta_c\bar \Delta_s}{4}\left(1+\frac{\kappa^2}{\Delta_c^2}\right)\left(1+\frac{\gamma_s^2}{4\bar \Delta_s^2}\right) }
\]
[1405.3289].

Methodologically, the field combines mean-field theory, Holstein–Primakoff bosonization, linear stability analysis, Routh–Hurwitz criteria, Keldysh functional integrals, discrete truncated Wigner approximation, and, in single-mode problems, Liouvillian spectral analysis [2508.10296][2607.00557][1611.09378][1503.04672][2307.05675]. The Liouvillian-spectral work on the single-mode open Dicke model is not a lattice study, but it is notable for showing that complex spectral windows can display 2D Poisson statistics in regular regimes and GinUE statistics in chaotic ones, together with an eigenstate-based convergence criterion for bosonic Liouvillians [2307.05675]. This suggests a possible route for future open Dicke lattice studies of nonequilibrium chaos, though the spatial case is substantially more difficult because of hopping, larger Hilbert spaces, and momentum structure.

A persistent source of confusion is the use of “open Dicke” for models that are not lattices. Several important papers treat a single cavity mode coupled collectively to matter and then add photon loss, matter baths, or local spin dissipation [1611.09378][1503.04672][2110.13191]. These are not Dicke lattice models, because they lack photon hopping, site-resolved cavity fields, and lattice momentum structure. Their relevance lies instead in identifying local or mode-resolved mechanisms—bath-induced cavity softening, dissipative critical exponents, local-dissipation-induced sector selection, or open-system correlation structure—that can serve as building blocks for genuine open Dicke lattices [1611.09378][1503.04672][2505.12696][2110.13191].

In current usage, then, the open Dicke lattice model is best understood not as a single universal Hamiltonian but as a family of dissipative lattice light–matter theories defined by three ingredients: local Dicke coupling, spatially extended photonic dynamics, and explicit openness. Across implementations and formulations, the field’s most characteristic results are the emergence of finite-momentum and boundary-induced superradiant phases, the organization of steady states into symmetry-related configuration classes, and the fact that dissipation can create qualitatively new nonequilibrium phase structure rather than merely perturb the equilibrium Dicke transition [1405.3289][2508.10296][2607.00557].

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