---
title: Closed Dicke Lattice Dynamics
url: https://www.emergentmind.com/topics/closed-dicke-lattice
type: topic
---

# Closed Dicke Lattice Dynamics

Searching arXiv for recent and foundational papers on "closed Dicke lattice" and related usages.
“Closed Dicke lattice” is not a single standardized object but a family of closely related constructions built from Dicke or Tavis–Cummings light–matter models under Hamiltonian, non-dissipative dynamics. In one usage, “closed” means a Hermitian Dicke Hamiltonian and “lattice” means a Peres lattice, namely a plot of eigenstate expectation values versus energy used to diagnose integrability, chaos, and excited-state singularities. In another, it denotes ordered eigenstate lattices in the integrable Tavis–Cummings limit, where monodromy produces topological defects in energy–momentum or Peres plots. In a third, it refers to genuine spatial lattices or rings of coupled Dicke units, whose equilibrium superradiant phases are selected by photon hopping and lattice symmetry. Further uses occur in momentum-space “superradiance lattices” built from timed Dicke states and in ordered atomic arrays where superradiant avalanches effectively close competing decay channels [1312.2555] [1702.07224] [2607.00557] [1403.7097] [2304.00093].

## 1. Terminological scope and defining features

Across these usages, “closed” always denotes the absence of dissipation, baths, or Lindblad terms. In the single-mode Dicke studies, it means unitary dynamics of a Hermitian Hamiltonian. In coupled-cavity Dicke lattices it is implemented by setting the photon loss rate to zero, so that phase structure is determined by equilibrium energy minimization and Gaussian fluctuations rather than by steady-state balance. In the Tavis–Cummings monodromy setting, it similarly refers to an isolated Hamiltonian system whose regular lattices of eigenstate expectation values are meaningful precisely because one is working with Hamiltonian eigenstates rather than Liouvillian modes [1312.2555] [1702.07224] [2607.00557].

The word “lattice” is used in at least three technically distinct senses. First, it can mean a Peres lattice, i.e. a visual grid-like arrangement of points obtained by plotting $I_O(n)=\langle \psi_n|O|\psi_n\rangle$ against the eigenenergy $E_n$. Second, it can mean a joint spectral or expectation-value lattice in the integrable Tavis–Cummings model, typically organized by energy and the conserved excitation number $M$. Third, it can mean a genuine spatial lattice of coupled cavity–spin units, such as a one-dimensional ring with periodic boundary conditions, a trimer with flux, or a Dicke–Ising chain or square lattice. The superradiance-lattice literature introduces yet another usage: a tight-binding lattice in momentum space whose sites are timed Dicke states generated by EIT standing waves, while the ordered-array literature uses “closed” in the different sense that cooperative decay effectively closes weaker internal transitions [1403.7097] [2304.00093].

This semantic multiplicity is not merely lexical. It separates three different research programs: spectral diagnostics of chaos and ESQPTs, topological defects in integrable spectra, and equilibrium many-body phase structure of spatially extended Dicke systems. The same phrase therefore names different objects depending on whether the “lattice” is a visualization, a quantum-number grid, a momentum-space construction, or a real-space array.

## 2. Peres lattices in the closed Dicke model

In the single-mode Dicke setting, the closed model is studied in the symmetric atomic subspace $j=N/2$ with
$$
H=\omega a^\dagger a+\omega_0 J_z+\frac{2\gamma}{\sqrt{N}}(a+a^\dagger)J_x,
$$
with $[a,a^\dagger]=1$ and $[J_i,J_j]=i\epsilon_{ijk}J_k$. The Hamiltonian commutes with the parity operator $\Pi=\exp(i\pi\Lambda)$, where $\Lambda=a^\dagger a+J_z+j$, so parity is used to block-diagonalize the spectrum. Because $J_x$ flips parity, $J_x$ itself is not used as a Peres operator; $J_x^2$ is used instead. Integrable limits occur for $\gamma\to0$ and $\omega_0\to0$, and the rotating-wave Tavis–Cummings approximation restores conservation of $\Lambda$. In the thermodynamic limit the normal-to-superradiant critical coupling is
$$
\gamma_c=\frac{1}{2}\sqrt{\omega\omega_0},
$$
which becomes $\gamma_c=0.5$ at resonance $\omega=\omega_0=1$ [1312.2555].

A Peres lattice plots, for each eigenstate $|\psi_n\rangle$, the expectation value
$$
I_O(n)=\langle\psi_n|O|\psi_n\rangle
$$
against $E_n$, equivalently the long-time average of $O$ for that eigenstate. The observables used are $O\in\{J_z,J_x^2,n\}$ with $n=a^\dagger a$, and lattices are constructed separately in each parity sector. Ordered point patterns indicate regular dynamics, while scrambled point clouds indicate quantum chaos; local distortions inside otherwise ordered regions reveal mixed regular–chaotic structure [1312.2555].

At resonance and for $N=40$ atoms ($j=20$), the Peres lattices are highly regular near $\gamma=0.01\gamma_c$, consistent with integrability. The eigenenergies follow $E(n,j,m)\approx n+m$, and degeneracies appear as vertical point stacks at identical energies but different expectation values. For nonzero couplings below the superradiant threshold, regular lattices persist at lower energies while irregular regions appear at higher energies, notably for $E/j\gtrsim0$. The key point is that chaos appears well before the equilibrium superradiant transition. For $\gamma>\gamma_c$, global irregularity increases, but an ordered low-energy island develops below $E/j\approx-1$ in all three lattices and grows with increasing coupling; at $\gamma=2\gamma_c$ and $3\gamma_c$ the plots show a sharp visual separation between this ordered region and chaotic higher-energy sectors [1312.2555].

The same Peres geometry exposes two ESQPTs. In the $J_z$ lattice, structural slope changes occur at $E/j=+1$ and $E/j=-1$. The static ESQPT at $E/j=+1$ is associated with saturation of the atomic excitation manifold and is independent of $\gamma$. The dynamic ESQPT at $E/j=-1$ depends on $\gamma$ and is tied to the fact that the configuration $(n=0,m=-j)$, which is the ground state in the normal phase, ceases to be the ground state in the superradiant phase; strong coupling populates many-photon, many-excitation states with $E/j<-1$. The emergence and growth of the regular island below $E/j=-1$ are visually linked to that dynamic ESQPT [1312.2555].

The numerical diagonalization uses an extended displaced bosonic basis defined by
$$
A=a+\frac{2\gamma}{\sqrt{N\omega}}J_x,
$$
with basis states $|N;j,m\rangle$ satisfying $A^\dagger A|N;j,m\rangle=N|N;j,m\rangle$, and parity-adapted combinations
$$
|N;j,m;\pm\rangle=\frac{|N;j,m\rangle+(-1)^N|N;j,-m\rangle}{\sqrt{2(1+\delta_{m,0})}}.
$$
Typical truncations are $N_{\max}=250$–$300$, and the reported truncation error bound is $\Delta P\lesssim10^{-30}$ for the displayed spectral windows. The method yields large sets of converged excited states and makes the fine structure of regular, chaotic, and ESQPT-organized sectors directly visible [1312.2555].

## 3. Energy–momentum lattices, monodromy, and the integrable Tavis–Cummings sector

A different closed Dicke-lattice construction appears in the integrable limit of the extended Dicke model with counter-rotating strength $\delta$. Setting $\delta=0$ gives the Tavis–Cummings Hamiltonian, for which the $U(1)$ excitation number
$$
M=b^\dagger b+J_3+j
$$
is conserved. In this regime, eigenstates can be organized simultaneously by energy and $M$, and both energy–momentum maps and Peres lattices form ordered arrays with topological defects related to monodromy rather than to chaos. The corresponding classical large-$j$ limit has two degrees of freedom, and after a canonical transformation the dynamics at fixed $\mathcal M=M/N^*$ reduces to an effective one-degree-of-freedom Hamiltonian [1702.07224].

At $\mathcal M=1$ and above the coupling threshold
$$
\lambda'_c=\frac{1}{2}|\Delta\omega|,
$$
the point $(x',p')=(0,0)$ becomes an unstable stationary point at energy $\mathcal E'_c=\omega_0/2$. The critical classical fiber $(\mathcal M,\mathcal E)=(1,\mathcal E'_c)$ is a singular pinched torus of focus–focus type. Near the north pole of the Bloch sphere the motion spirals, takes infinite time to pass the focus, and realizes the closed-system analogue of dynamic superradiance: pulse-like exchange of excitations between field and atoms, slow near the poles and fastest near the equator. In the tuned case $\Delta\omega=0$, the photon emission or absorption rate scales linearly with $N^*$ in this closed setting [1702.07224].

Quantum mechanically, the joint spectrum $(E,M)$ forms a lattice with a point defect at $(M,E)=(N^*,N^*\mathcal E'_c)$. Encircling this defect transforms the lattice cell by
$$
\begin{pmatrix}\vec M'\\ \vec k'\end{pmatrix}
=
\begin{pmatrix}
1 & 1\\
0 & 1
\end{pmatrix}
\begin{pmatrix}\vec M\\ \vec k\end{pmatrix},
$$
which is the transpose of the classical monodromy matrix for one focus–focus singularity. Peres lattices built from $\langle b^\dagger b\rangle_k$ or $\langle J_3\rangle_k$ versus $E_k$ show corresponding dislocations, although the local topology of the defect depends on the chosen observable. In the integrable case, chains at fixed principal quantum number across varying $M$ make the kinked structure near the defect especially clear [1702.07224].

The same energy marks an ESQPT. Within the critical $M=N^*$ subspace, the unstable stationary point produces an ESQPT at
$$
E_c=N^*\mathcal E'_c=\omega_0 j.
$$
The semiclassical density of states in that subspace has a logarithmic divergence at $E=E_c$, while in the full two-degree-of-freedom spectrum the first derivative of the level density has a downward jump. The alignment of the ESQPT energy with the monodromy defect increases the local spectral density and amplifies the visibility of the defect in both energy–momentum and Peres lattices [1702.07224].

Once counter-rotating terms are restored, $U(1)$ symmetry is broken and only parity remains. Classically, the pinched-torus neighborhood becomes chaotic already for very small $\delta>0$, and quantum lattices rapidly lose their point-defect structure. The monodromy point lies on a line of minimal level spacing, so generic perturbations induce strong level repulsion, tear the lattice along a vertical break rooted at the monodromy energy, and disorder the Peres plots [1702.07224].

## 4. Closed Dicke lattices as coupled cavity–spin arrays

In the spatial sense, a closed Dicke lattice is a ring or array of local Dicke units coupled by photon hopping. A representative Hamiltonian 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,
$$
periodic boundary conditions $c_{N+1}=c_1$, and large-$N_a$ equilibrium analysis based on mean field plus Holstein–Primakoff theory. In the superradiant regime each site develops a real coherent field $\alpha_i=\langle c_i\rangle/\sqrt{N_a}$, and the sign pattern of the $\alpha_i$ defines the “configuration” class. For $N=4$ the symmetry-inequivalent classes are $[++++]$, $[+---]$, $[++--]$, and $[+-+-]$, with degeneracies $2$, $8$, $4$, and $2$, respectively [2607.00557].

The closed-lattice ground state is determined by minimizing
$$
\frac{E}{N_a}=\sum_{j=1}^N\left[\omega_c\alpha_j^2-\frac{1}{2}\sqrt{\omega_a^2+16g^2\alpha_j^2}-2\xi\alpha_j\alpha_{j+1}\right].
$$
The on-site term favors nonzero $\alpha_j$ above threshold, while the hopping term enforces sign alignment or alternation. For $\xi>0$, the energy is minimized by the uniform ferromagnetic-like configuration $[++++]$; for $\xi<0$, it is minimized by the staggered antiferromagnetic-like configuration $[+-+-]$. The Hessian eigenvalues around the normal phase are
$$
\lambda_k=2\omega_c-\frac{8g^2}{\omega_a}-4\xi\cos(2\pi k/N),
$$
which gives the mode-resolved threshold
$$
g_c(k)=\frac{1}{2}\sqrt{\omega_a\left[\omega_c-2\xi\cos(2\pi k/N)\right]}.
$$
For $N=4$, this reduces to $g_c^\pm=\frac{1}{2}\sqrt{\omega_a(\omega_c\mp2\xi)}$. The normal-to-superradiant transition is continuous, whereas scanning $\xi$ through zero at fixed $g>g_c$ produces a first-order transition between the ferromagnetic and antiferromagnetic superradiant configurations. The lowest excitation gap and the ground-state photon fluctuations scale with exponent $1/2$, matching the equilibrium mean-field Dicke universality class [2607.00557].

A closely related closed Dicke lattice model arises in hybrid arrays of superconducting microwave cavities coupled to ensembles of NV centers. The effective Hamiltonian is
$$
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^\dagger_{\ell+1})
+\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 cavity dispersion $\Delta_k=\Delta_c-2t\cos k$. In the closed limit $\kappa=0$, the momentum-resolved instability threshold is
$$
G_{\rm crit}(k)=\frac{1}{2}\sqrt{\Delta_s\Delta_k},
$$
and the global threshold is the minimum over $k$. For a 1D nearest-neighbor chain with $t>0$, the minimum lies at $k=0$, so the equilibrium superradiant phase is uniform. The open-system finite-$k$ instabilities emphasized in the same work are explicitly described as non-equilibrium effects absent in the closed Hamiltonian limit [1405.3289].

The trimer provides the minimal odd ring. Its closed Hamiltonian,
$$
H=\sum_{n=1}^3 H_n^{\rm Dicke}
+J\sum_{n=1}^3\left(e^{i\varphi}a_n^\dagger a_{n+1}+e^{-i\varphi}a_{n+1}^\dagger a_n\right),
$$
contains unbalanced rotating and counter-rotating couplings through $\eta_\pm=(1\pm\eta)/2$. At $\eta=0$, the model restores the $U(1)$ symmetry of the Tavis–Cummings limit and the superradiant phase contains a Goldstone mode. The classical ground-state functional supports a translationally invariant non-frustrated superradiant phase and a translation-breaking frustrated superradiant phase. The normal-phase Hessian has six exact eigenvalues, yielding non-frustrated and frustrated critical branches $g_c^{\rm nf}$ and $g_c^{\rm f}$; the actual threshold is their minimum. At time-reversal-symmetric $\varphi=0$, two soft modes emerge at the normal-phase boundary with exponent $1/2$, while on the frustrated side the low-energy mode scales as $| \delta g|^{3/2}$. For $0<\varphi<\varphi_{\rm tr}$, the normal-phase soft mode instead scales as $|\delta g|^{1}$ and exhibits anomalous finite critical fluctuations because the conjugate quadratic stiffnesses vanish with the same exponent [2303.11758].

## 5. Closed equilibrium Dicke lattices with additional matter interactions

The equilibrium Dicke–Ising model extends the closed Dicke lattice by adding short-range Ising couplings among the spins while retaining a single quantized cavity mode. The Hamiltonian is
$$
H=\omega a^\dagger a+\epsilon\sum_{i=1}^N\sigma_i^z+\frac{g}{\sqrt N}(a+a^\dagger)\sum_{i=1}^N\sigma_i^x+J\sum_{\langle ij\rangle}\sigma_i^z\sigma_j^z.
$$
It is studied on 1D chains and 2D square lattices with periodic boundary conditions at zero temperature using sign-problem-free QMC. For $\epsilon=0$ the model has $Z_2\times Z_2$ symmetry: Ising spin-flip in the $z$ sector and Dicke parity in the light–matter sector. The natural order parameters are the superradiant amplitude $\alpha=\langle a+a^\dagger\rangle/\sqrt N$, the photon density $n_{\rm ph}=\langle a^\dagger a\rangle/N$, the uniform magnetization $m$, and the staggered magnetization $m_s$ [2409.15082].

For ferromagnetic $J<0$, any $\epsilon\neq0$ explicitly breaks the spin $Z_2$, so only one phase boundary remains, separating a normal phase from a superradiant phase. At small longitudinal field $\epsilon\ll |J|$ the transition is first order, whereas at large $\epsilon\gg|J|$ it is continuous and follows Dicke criticality. The line separating these regimes ends at a multicritical point at $\epsilon\approx|J|$. In the continuous regime, QMC agrees with the infinite-dimensional variational mean-field condition
$$
\epsilon=2\left(\frac{g^2}{\omega}+Jd\right),
$$
while the first-order segment at small $\epsilon$ is a genuine beyond-mean-field effect [2409.15082].

For antiferromagnetic $J>0$, both $Z_2$ symmetries remain intact even at finite $\epsilon$, giving four phases: PN, AN, PS, and AS. The AS phase carries simultaneous superradiant and AFM order and is identified as the light–matter analogue of a lattice supersolid. At $\epsilon=0$, the AN–PS transition is first order and persists to finite $\epsilon$ until it intersects the mean-field AN–AS line, where it splits and opens an intermediate AS region. On the square lattice, the AS–PS boundary is first order for $g>0.30(2)$ and continuous with 3D Ising universality for $g<0.30(2)$; on the chain it remains first order throughout the explored parameter range [2409.15082].

All continuous superradiant transitions in this closed model are in the Dicke universality class,
$$
\alpha=0,\quad \beta=\frac12,\quad \gamma=1,\quad \delta=3,\quad \nu=\infty,\quad z\nu=\frac12,
$$
with the modified finite-size scaling expected above the upper critical dimension. The effective exponent entering drift and rounding is $\nu'=3/2$, leading to
$$
g_c(N)-g_c(\infty)\sim N^{-2/3},\qquad
\Delta_c(N)\sim N^{-1/3},\qquad
\alpha_c(N)\sim N^{-1/6}.
$$
The QMC formulation integrates out the photon exactly and yields a retarded all-to-all transverse interaction among spins, which is sampled by a wormhole generalization of the directed-loop algorithm. This provides unbiased equilibrium benchmarks for closed Dicke lattices with competing local and global interactions [2409.15082].

## 6. Alternative but related constructions

The phrase also appears in settings that are not real-space Hamiltonian Dicke lattices. In the superradiance-lattice construction, $N$ fixed $\Lambda$-type atoms support timed Dicke states
$$
|e_{\mathbf k}\rangle=\frac{1}{\sqrt N}\sum_{j=1}^N e^{i\mathbf k\cdot \mathbf r_j}|g_1\cdots e_j\cdots g_N\rangle,
\qquad
|s_{\mathbf k}\rangle=\frac{1}{\sqrt N}\sum_{j=1}^N e^{i\mathbf k\cdot \mathbf r_j}|g_1\cdots s_j\cdots g_N\rangle,
$$
which act as momentum-space lattice sites. A standing-wave EIT coupling generates nearest-neighbor hopping in the timed-Dicke basis,
$$
H=\sum_m(\Delta_e|e_m\rangle\langle e_m|+\Delta_s|s_m\rangle\langle s_m|)
+\sum_m\left(\frac{\Omega_+}{2}|e_m\rangle\langle s_{m+1}|+\frac{\Omega_-}{2}|e_m\rangle\langle s_{m-1}|+\text{h.c.}\right).
$$
Detuning asymmetry produces a uniform force in momentum space with Bloch frequency $\omega_B=2\delta$, leading to Bloch oscillations and Wannier–Stark ladders, while periodic modulation renormalizes the hopping through Bessel functions and yields band collapse and dynamic localization. The system is operationally “closed” because EIT confines dynamics to a dark-state manifold, suppressing spontaneous emission from $|e\rangle$ and leaving ground-state dephasing as the principal residual decoherence [1403.7097].

A different usage appears in ordered two-dimensional arrays of multilevel alkaline-earth(-like) atoms. There, the array is not a Hamiltonian Dicke lattice but a free-space cooperative emitter governed by a channel-resolved Lindblad master equation. “Closed” refers to transition closing: a superradiant avalanche can funnel emission overwhelmingly into one dominant branch, effectively suppressing weaker fine-structure or Zeeman channels. For a channel $\alpha$, the effective branching ratio is
$$
B_\alpha^{\rm eff}=\frac{\Gamma_\alpha^{\rm col}}{\sum_\beta \Gamma_\beta^{\rm col}},
$$
with $B_{\rm bright}^{\rm eff}\to1$ under strong cooperative enhancement. In square arrays at $d=0.2\lambda_0$, the integrated bright-branch share increases from approximately $60\%$ for a single Sr atom in the quoted multilevel example to approximately $70\%$ for a $12\times12$ array, while predicted burst intensities scale superlinearly, with exponents such as $N^{1.29}$–$N^{1.47}$ depending on species, spacing, and polarization. Directional superradiance can persist even when global superradiance is weak, and shows revivals near $d\approx \lambda/2$ and $d\approx\lambda$ due to geometric resonances [2304.00093].

These constructions remain adjacent to the main closed-Dicke-lattice program because they retain collective Dicke-state organization while relocating the notion of “lattice” from real space to momentum space, or the notion of “closed” from Hamiltonian isolation to effective suppression of alternative decay channels. They therefore extend, rather than unify, the semantics of the term.

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