Subradiant Correlations in Quantum Emitter Arrays

• Subradiant correlations are collective interference phenomena in many-emitter systems that suppress radiative decay through engineered destructive interference.
• They are characterized by long-lived eigenmodes with decay rates significantly lower than individual emitter rates, observable in settings like bad-cavity lasers and ordered atomic arrays.
• Analytical methods such as Liouvillian spectral theory and non-Hermitian Hamiltonians are employed to study these effects for applications in metrology and photon storage.

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 (Shi et al., 12 Sep 2025, Shankar et al., 2021, Jen, 2016).

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 (Shi et al., 12 Sep 2025, Leppenen et al., 25 Jul 2025).

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 $|\Re\lambda|\ll\gamma$0 means that collective emission is suppressed relative to independent atoms, and the macroscopic singlet yields the maximally subradiant value $|\Re\lambda|\ll\gamma$1 (Shankar et al., 2021).

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 (Jen et al., 2016, Jen, 2016).

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 $|\Re\lambda|\ll\gamma$2 two-level atoms at positions $|\Re\lambda|\ll\gamma$3, the guided environment enters through collective jump operators

$|\Re\lambda|\ll\gamma$4

with $|\Re\lambda|\ll\gamma$5, and the master equation

$|\Re\lambda|\ll\gamma$6

Here $|\Re\lambda|\ll\gamma$7 is the coherent drive and $|\Re\lambda|\ll\gamma$8 is the waveguide-mediated exchange interaction. The off-diagonal structure of $|\Re\lambda|\ll\gamma$9 means that emission is a collective process rather than a sum of independent jumps (Shi et al., 12 Sep 2025).

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 $|\Re\lambda|$0, supplemented by individual decay and repumping. There, subradiance is the suppression of $|\Re\lambda|$1 relative to the independent-emitter contribution, and the collective Dicke-manifold structure provides a natural basis for describing dark and nearly dark correlations (Shankar et al., 2021).

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,

$|\Re\lambda|$2

which produces one bright collective state and $|\Re\lambda|$3 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 (Biella et al., 2013).

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 (Jenkins et al., 2016, Facchinetti et al., 2016).

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 $|\Re\lambda|$4, with

$|\Re\lambda|$5

and the drive defines

$|\Re\lambda|$6

For the oscillatory branches $|\Re\lambda|$7, the exact bound

$|\Re\lambda|$8

rules out oscillating subradiant correlations in the strong-drive regime $|\Re\lambda|$9. The result is independent of $$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$0: oscillatory modes can never become arbitrarily long-lived with increasing array size, whereas the nonoscillating branch $$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$1 can still host subradiant correlations with $$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$2 (Shi et al., 12 Sep 2025).

This strong-drive prohibition is model-specific rather than universal. In a random driven Dicke model with collective decay, strong coherent drive $$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$3 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

$$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$4

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 $$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$5, producing slowly decaying oscillatory correlations in finite-size systems (Leppenen et al., 25 Jul 2025).

A related expansion of the subradiant domain occurs in driven anti-Bragg waveguide QED. For periodic arrays with $$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$6 or $$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$7, 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 (Poddubny, 2022).

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 (Shi et al., 12 Sep 2025, Leppenen et al., 25 Jul 2025, Poddubny, 2022).

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-$$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$8 sectors and represent the dissipative block as

$$S_f=(N-1)\langle \hat{\sigma}_1^+\hat{\sigma}_2^-\rangle$$9

where $\Im\lambda\neq 0$0 is diagonal and strictly positive,

$\Im\lambda\neq 0$1

and $\Im\lambda\neq 0$2 is a positive semidefinite Laplacian built from a weighted incidence operator on a ranked poset of word pairs $\Im\lambda\neq 0$3. The positivity of $\Im\lambda\neq 0$4, combined with the anti-Hermitian character of $\Im\lambda\neq 0$5, yields the spectral bound on $\Im\lambda\neq 0$6 via the Bendixson inequality. The proof makes no use of the detailed form of $\Im\lambda\neq 0$7, only of the Liouvillian block structure and the poset/Laplacian decomposition (Shi et al., 12 Sep 2025).

In the random driven Dicke model, exact diagonalization of the $\Im\lambda\neq 0$8 Liouvillian is supplemented by group representation theory. Depending on the Hamiltonian, the relevant symmetry is $\Im\lambda\neq 0$9, $$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$0, or $$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$1, 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 (Leppenen et al., 25 Jul 2025).

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, Jen, 2016, Facchinetti et al., 2016).

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

$$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$2

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 $$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$3 three-dimensional array the lifetime of a subradiant De Moivre state reaches $$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$4 milliseconds (Jen et al., 2016, Jen, 2016).

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

$$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$5

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 $$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$6 displays persistent connected density-density correlations

$$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$7

and modified third-order correlations

$$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$8

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 (Jen, 2021).

Steady-state subradiance can also organize the entire many-body Hilbert space. In the bad-cavity laser, there is a dissipative phase transition at $$\mathcal{L}\rho_\alpha=\lambda_\alpha \rho_\alpha,\qquad \rho_\alpha(t)\propto e^{\lambda_\alpha t},$$9 between two distinct subradiant phases. Both have negative $-\Re\lambda_\alpha$0 and approach the maximally subradiant value $-\Re\lambda_\alpha$1, but they differ qualitatively: one is concentrated at low $-\Re\lambda_\alpha$2 and extensive $-\Re\lambda_\alpha$3, the other near the singlet corner with $-\Re\lambda_\alpha$4. Near the critical region, the generalized spin-squeezing parameter

$-\Re\lambda_\alpha$5

satisfies $-\Re\lambda_\alpha$6, with $-\Re\lambda_\alpha$7 and $-\Re\lambda_\alpha$8 from exact diagonalization, indicating macroscopic entanglement and a vanishing fraction of unentangled atoms in the large-$-\Re\lambda_\alpha$9 limit (Shankar et al., 2021).

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 $\Im\lambda_\alpha$0. 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 (Biella et al., 2013).

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 (Facchinetti et al., 2016, Kornovan et al., 2019).

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

$\Im\lambda_\alpha$1

shows explicitly that the effective linewidth $\Im\lambda_\alpha$2 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 (Zafra-Bono et al., 9 Dec 2025).

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 (Jen et al., 2016, Jen, 2016, Facchinetti et al., 2016, 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 $\Im\lambda_\alpha$3 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 $\Im\lambda_\alpha$4, where $\Im\lambda_\alpha$5 is the collective Lamb shift (Feng et al., 2017).

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 $\Im\lambda_\alpha$6 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 [(Shi et al., 12 Sep 2025); (Leppenen et al., 25 Jul 2025); (Biella et al., 2013)].