---
title: Subradiant Correlations in Quantum Emitter Arrays
url: https://www.emergentmind.com/topics/subradiant-correlations
type: topic
---

# Subradiant Correlations in Quantum Emitter Arrays

Subradiant correlations are collective correlation structures in open many-emitter systems whose radiative decay is suppressed by destructive interference in a common optical environment. The term is used in several technically distinct but related senses. In Liouvillian formulations, it denotes correlation modes with eigenvalues \(\lambda\) satisfying \(|\Re\lambda|\ll\gamma\), so that temporal correlations decay much more slowly than the single-emitter rate. In bad-cavity lasers, it is quantified by negative collective pair coherence and a negative subradiance factor. In ordered atomic arrays, it is realized as delocalized spin-wave or phase-imprinted excitations whose collective linewidth is smaller than the single-atom linewidth. Across waveguide QED, free-space lattices, Dicke-type models, metamaterials, and disordered open systems, the common element is a many-body excitation or eigenoperator that couples weakly to the radiation continuum [2509.09993; 2103.07402; 1606.06403].

## 1. Definitions and conceptual scope

The phrase “subradiant correlations” does not refer to a single universal observable. Its precise meaning depends on the level of description, the excitation sector, and whether one studies states, eigenoperators, or steady-state correlation functions.

| Context | Definition | Hallmark |
|---|---|---|
| Driven waveguide QED | Liouvillian modes with small \(|\Re\lambda|\) | Slow decay of temporal correlations |
| Bad-cavity laser | Negative \(S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle\) | Suppressed cavity output |
| Phase-imprinted arrays | Collective single-excitation spin waves | Small collective linewidth |
| Random driven Dicke model | Long-lived traceless Liouvillian eigenoperators | Small Liouvillian gap, possibly \(\Im\lambda\neq 0\) |

In the Liouvillian language, one writes
\[
\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},
\]
so that \(-\Re\lambda_\alpha\) is the decay rate and \(\Im\lambda_\alpha\) is the oscillation frequency. In this sense, subradiant correlations are not density matrices of physical states but traceless eigenoperators governing slow components of expectation values and two-time correlators [2509.09993; 2507.19467].

A different but compatible definition is used in the bad-cavity laser, where subradiance is encoded directly in collective observables. The subradiance factor
\[
S_f=\frac{1}{N}\left[\langle \hat{J}^+\hat{J}^- \rangle - \left(\frac{N}{2} + \langle \hat{J}^z \rangle\right)\right]
\]
can be rewritten as \(S_f=(N-1)\langle \hat{\sigma}_1^+ \hat{\sigma}_2^- \rangle\). Negative \(S_f\) means that collective emission is suppressed relative to independent atoms, and the macroscopic singlet yields the maximally subradiant value \(S_f=-1/2\) [2103.07402].

In single-excitation array problems, subradiant correlations are often represented as nonlocal phase patterns. In one- and three-dimensional lattices, “De Moivre” states distribute one excitation over all atoms with site-dependent phases; some of these states overlap predominantly with eigenmodes whose collective decay rates are far below the single-atom rate, thereby realizing long-lived spin-wave correlations [1603.00996; 1606.06403].

## 2. Microscopic origin: common channels, interference, and collective operators

The microscopic source of subradiant correlations is coherent coupling of many emitters to shared radiative channels. In a waveguide-QED array of \(N\) two-level atoms at positions \(z_j\), the guided environment enters through collective jump operators
\[
c_L = \sum_{j=1}^N e^{i\phi_j}\sigma_j,\qquad
c_R = \sum_{j=1}^N e^{-i\phi_j}\sigma_j,
\]
with \(\phi_j=k z_j\), and the master equation
\[
\partial_t \rho = -i[H_{\rm I}+H_{\rm II},\rho] + \gamma\mathcal{D}[c_L]\rho+\gamma\mathcal{D}[c_R]\rho.
\]
Here \(H_{\rm I}\) is the coherent drive and \(H_{\rm II}\) is the waveguide-mediated exchange interaction. The off-diagonal structure of \(c_{L,R}\) means that emission is a collective process rather than a sum of independent jumps [2509.09993].

Closely related structures appear in cavity-mediated models. In the bad-cavity laser, adiabatic elimination of the cavity yields a purely dissipative master equation with a collective jump \(\hat{J}^-=\sum_j \hat{\sigma}_j^-\), supplemented by individual decay and repumping. There, subradiance is the suppression of \(\langle \hat{J}^+\hat{J}^- \rangle\) relative to the independent-emitter contribution, and the collective Dicke-manifold structure provides a natural basis for describing dark and nearly dark correlations [2103.07402].

In disordered open systems, the same logic reappears in effective non-Hermitian Hamiltonians. In the open 3D Anderson-Dicke model, coherent disorder is combined with a rank-1 opening,
\[
H_{\rm eff}=H_0 - i\frac{\gamma}{2}Q,\qquad Q_{ij}=1,
\]
which produces one bright collective state and \(N-1\) dark or weakly radiative states. The opening simultaneously induces long-range hopping and a collective decay channel, so subradiant correlations arise from amplitudes arranged to minimize overlap with the bright vector [1304.5451].

An analogous eigenmode picture holds in metamaterial arrays and free-space atomic arrays, where the basic object is a complex interaction matrix built from dyadic Green’s functions. Its eigenvectors are collective current or dipole patterns, and its eigenvalues furnish collective linewidths and shifts. Subradiant correlations are then the phase-coherent eigenvectors whose radiative linewidth is strongly reduced [1611.01509; 1609.08350].

## 3. Liouvillian spectra, oscillatory branches, and driven regimes

A central development in recent work is the shift from subradiant states to subradiant Liouvillian modes. In a strongly driven waveguide-QED array, the Liouvillian spectrum organizes into bands labelled by an integer \(m\), with
\[
\Im\lambda \approx 2m|\Omega|,
\]
and the drive defines
\[
H_{\rm I}=-2\Omega J_x.
\]
For the oscillatory branches \(m\neq 0\), the exact bound
\[
|\Re\lambda| \ge \frac{m\gamma}{2}
\]
rules out oscillating subradiant correlations in the strong-drive regime \(|\Omega|\gg N\gamma\). The result is independent of \(N\): oscillatory modes can never become arbitrarily long-lived with increasing array size, whereas the nonoscillating branch \(m=0\) can still host subradiant correlations with \(|\Re\lambda|\ll\gamma\) [2509.09993].

This strong-drive prohibition is model-specific rather than universal. In a random driven Dicke model with collective decay, strong coherent drive \(\Omega\gg\delta\omega\) dynamically suppresses the effect of inhomogeneous broadening and restores a large manifold of long-lived Liouvillian modes. In the drive-dominated basis, these modes take the form
\[
\rho(j)_d^{\nu,\nu'}=\frac{1}{2j+1}\sum_{m_x=-j}^{j}\ket{j,m_x,\nu}\bra{j,m_x,\nu'},
\]
and they remain long-lived even when conventional dark states are destroyed by frequency disorder. When nearest-neighbor dipole-dipole interactions are added, some of these long-lived modes acquire nonzero \(\Im\lambda\), producing slowly decaying oscillatory correlations in finite-size systems [2507.19467].

A related expansion of the subradiant domain occurs in driven anti-Bragg waveguide QED. For periodic arrays with \(d=\lambda/4\) or \(3\lambda/4\), there are no such states at low driving powers, but strong coherent driving generates strongly subradiant eigenstates of the master equation and directly manifests them in long-living quantum correlations between qubit excitations [2202.10138].

Taken together, these results distinguish three regimes. First, conventional subradiance can arise already in the linear or weak-drive spectrum. Second, strong driving can create new long-lived Liouvillian modes not present in the linear regime. Third, strong driving can also forbid a particular class of oscillatory subradiant correlations, as in the waveguide array with left- and right-propagating collective jumps. A common misconception is therefore that “strong driving either always destroys or always stabilizes subradiance”; the literature instead shows that the answer depends on Liouvillian structure, symmetry, and the definition of subradiance being used [2509.09993; 2507.19467; 2202.10138].

## 4. Mathematical structures and analytical methods

The mathematical analysis of subradiant correlations spans spectral theory, combinatorics, group representation theory, and Green-function methods. In the strong-drive waveguide array, the key step is to project the Liouvillian into fixed-\(m\) sectors and represent the dissipative block as
\[
\mathcal{Q}=-\mathcal{F}-\mathcal{A},
\]
where \(\mathcal{F}\) is diagonal and strictly positive,
\[
\mathcal{F}_\rho=\frac{1}{4}|a_w-a_{w'}|^2+\frac{m}{2},
\qquad
a_w=\sum_{j=1}^N (-1)^{w_j}e^{i\phi_j},
\]
and \(\mathcal{A}=\mathcal{B}^T\mathcal{B}\) is a positive semidefinite Laplacian built from a weighted incidence operator on a ranked poset of word pairs \((w,w')\). The positivity of \(\mathcal{A}\), combined with the anti-Hermitian character of \(-i[H_{\rm II},\cdot]\), yields the spectral bound on \(\Re\lambda\) via the Bendixson inequality. The proof makes no use of the detailed form of \(H_{\rm II}\), only of the Liouvillian block structure and the poset/Laplacian decomposition [2509.09993].

In the random driven Dicke model, exact diagonalization of the \(4^N\times4^N\) Liouvillian is supplemented by group representation theory. Depending on the Hamiltonian, the relevant symmetry is \(S_N\), \(D_N\), or \(C_s\), and the decomposition into irreducible representations determines how many dark or long-lived modes survive and which of them become oscillatory when additional coherent interactions split degenerate sectors. This approach yields explicit counting rules for both dark manifolds and oscillation frequencies [2507.19467].

In array scattering problems outside the full Liouvillian setting, the principal mathematical object is often an effective non-Hermitian Hamiltonian or Green-tensor kernel. In waveguide-coupled qubits, 3D phase-imprinted atomic lattices, and planar 2D arrays, its complex eigenvalues encode collective shifts and radiative widths, while its eigenvectors furnish the spatial phase patterns responsible for destructive interference. This framework underlies the identification of spin-wave subradiance, fermionized multi-excitation subradiance, and collective storage modes [1803.02115; 1606.06403; 1609.08350].

## 5. Spatial structures, many-body organization, and entanglement

Subradiant correlations are not restricted to a single spatial form. In one- and three-dimensional ordered arrays, they can be engineered by phase imprinting. The De Moivre states
\[
|\phi_m\rangle=\frac{1}{\sqrt{N}}\sum_{\mu=1}^{N}
e^{i\mathbf{k}\cdot \mathbf{r}_\mu}
e^{i\frac{2\pi m}{N}(\mu-1)}
|e_\mu\rangle |g\rangle^{\otimes(N-1)}
\]
form a complete orthonormal basis of the single-excitation manifold. Some of these states overlap predominantly with eigenmodes whose decay rates are much smaller than the free-space decay rate, and their fluorescence can exhibit decayed Rabi-like oscillations with beating frequency set by differences of cooperative Lamb shifts. For one hundred atoms, lifetimes up to hundred milliseconds were predicted in the one-dimensional optical-lattice proposal, and in a \(3\times 3\times 10\) three-dimensional array the lifetime of a subradiant De Moivre state reaches \(\sim 2\) milliseconds [1603.00996; 1606.06403].

In low-excitation one-dimensional waveguide QED, the spatial structure becomes strongly constrained by the hard-core nature of spin excitations. The most subradiant multi-excitation eigenstates are well approximated by fermionic or antisymmetrized combinations of single-excitation eigenstates, so that excitations avoid one another and avoid the boundaries. In this regime the most subradiant single-excitation modes obey the universal scaling
\[
\Gamma_\xi/\Gamma_{\rm 1D}\propto \xi^2/N^3,
\]
and multi-excitation subradiant eigenstates inherit strong real-space anti-bunching and long-lived temporal photon correlations [1803.02115].

In chiral waveguide QED with nonreciprocal couplings, subradiant correlations can bind multiple excitations into shape-preserving dimers and trimers. The long-time dynamics at \(\xi=\pi\) displays persistent connected density-density correlations
\[
\langle G^{(2)}(r)\rangle
\]
and modified third-order correlations
\[
\langle G^{(3)}\rangle,
\]
with ballistic but shape-preserving propagation. The diffusion speed depends on the initial coherence between the excited atoms and is robust to relative phase fluctuations [2102.03757].

Steady-state subradiance can also organize the entire many-body Hilbert space. In the bad-cavity laser, there is a dissipative phase transition at \(w_c=\gamma\) between two distinct subradiant phases. Both have negative \(S_f\) and approach the maximally subradiant value \(S_f\to -1/2\), but they differ qualitatively: one is concentrated at low \(M\) and extensive \(J\), the other near the singlet corner with \(J/N\to 0\). Near the critical region, the generalized spin-squeezing parameter
\[
\xi^2=\frac{(\Delta \hat{J}^x)^2+(\Delta \hat{J}^y)^2+(\Delta \hat{J}^z)^2}{N/2}
\]
satisfies \(\xi^2\ll 1\), with \(\xi^2_{\min}\propto N^{-\alpha}\) and \(\alpha\approx 0.34\) from exact diagonalization, indicating macroscopic entanglement and a vanishing fraction of unentangled atoms in the large-\(N\) limit [2103.07402].

Disorder does not simply eliminate subradiant correlations; it can also reshape them. In the open 3D Anderson-Dicke model, the subradiant hybrid regime combines an Anderson-localized core with an extended plateau of height \(\sim 1/N\). The participation ratio remains size-independent even though the state has a weak global background, showing that coherent opening and disorder can cooperate to produce hybrid subradiant states rather than purely localized or purely extended ones [1304.5451].

## 6. Spectral signatures, applications, and limitations

A recurring signature of subradiant correlations is the appearance of unusually narrow spectral features. In a planar 2D atomic array, high-fidelity preparation of a collective subradiant mode normal to the plane produces sharp transmission resonances and can be described by an effective two-mode model coupling a broad in-plane mode to a narrow perpendicular mode. In periodic one-dimensional arrays, extremely subradiant states produce very narrow transmission and reflection features and interaction-induced transparency in a narrow spectral range [1609.08350; 1906.07423].

The same narrowness can be turned into a metrological resource. In waveguides and subwavelength free-space arrays, subradiant collective states generate sharp transmission features that enhance sensitivity to global and spatially varying perturbations. The precision estimate
\[
\Delta\omega \sim \frac{\Gamma_{\rm sub}}{\sqrt{pN\Gamma_0\tau}}
\]
shows explicitly that the effective linewidth \(\Gamma_{\rm sub}\) of the subradiant feature replaces the bare single-emitter linewidth in the metrological scaling. Proposed applications include atomic clock operation, imaging of emitter positions, and detection of global or spatially varying detunings such as electromagnetic fields or gravitational gradients [2512.09050].

Subradiant correlations also support storage protocols. In phase-imprinted one- and three-dimensional arrays, a single photon can be mapped into a long-lived subradiant spin wave and later reconverted into a bright mode for readout. In 2D arrays, light storage is realized by transferring population from a bright collective mode into a highly subradiant mode through Zeeman-induced mixing. In 1D waveguide arrays, measurement protocols based on on-site readout and photon correlations provide direct access to real-space and temporal signatures of multi-excitation subradiant states [1603.00996; 1606.06403; 1609.08350; 1803.02115].

EIT-like constructions provide another route to observation. A superradiant state can act as the excited level and a subradiant state as the metastable level of an effective \(\Lambda\) scheme, so that the transparency point reveals the collective energy splitting and can be used for subwavelength metrology. In that setting the relevant frequency is \(\nu=2\Delta_c\), where \(\Delta_c\) is the collective Lamb shift [1701.08175].

The limitations are as instructive as the applications. Positional disorder in ordered atomic arrays broadens the narrowest linewidths and destroys the most extreme subradiant scaling. In the strong-drive waveguide array of uniformly driven emitters, oscillatory subradiant correlations are forbidden by the \( |\Re\lambda|\ge m\gamma/2 \) bound. In disordered Liouvillian systems, frequency disorder can destroy conventional dark states, although strong drive can dynamically restore long-lived Liouvillian correlations. A plausible implication is that “subradiant correlations” should be treated as a family of interference-protected many-body structures rather than as a single universal phase: their stability depends on geometry, symmetry, disorder, and the specific spectral object—state, eigenmode, or eigenoperator—under consideration [2509.09993; 2507.19467; 1304.5451].

Source: https://www.emergentmind.com/topics/subradiant-correlations