Ordered atomic arrays with subwavelength lattice spacing emit light collectively. For fully inverted atomic arrays, this results in an initial burst of radiation and a fast buildup of coherences between the atoms at initial times. Based on a cumulant expansion of the equations of motion, we derive exact analytical expressions for the emission properties and numerically analyze the full many-body problem resulting in the collective decay process for unprecedented system sizes of up to a few hundred atoms. We benchmark the cumulant expansion approach and show that it correctly captures the cooperative dynamics resulting in superradiance. For fully inverted arrays, this allows us to extract the scaling of the superradiant peak with particle number. For partially excited arrays where no coherences are shared among atoms, we also determine the critical number of excitations required for the emergence of superradiance in one- and two-dimensional geometries. In addition, we study the robustness of superradiance in the case of non-unit filling and position disorder.
Ordered atomic arrays with subwavelength lattice spacing emit light collectively. For fully inverted atomic arrays, this results in an initial burst of radiation and a fast buildup of coherences between the atoms at initial times. Based on a cumulant expansion of the equations of motion, we derive exact analytical expressions for the emission properties and numerically analyze the full many-body problem resulting in the collective decay process for unprecedented system sizes of up to a few hundred atoms. We benchmark the cumulant expansion approach and show that it correctly captures the cooperative dynamics resulting in superradiance. For fully inverted arrays, this allows us to extract the scaling of the superradiant peak with particle number. For partially excited arrays where no coherences are shared among atoms, we also determine the critical number of excitations required for the emergence of superradiance in one- and two-dimensional geometries. In addition, we study the robustness of superradiance in the case of non-unit filling and position disorder.
Superradiance and subradiance are collective effects that emerge from coherent interactions between quantum emitters. Due to their many-body nature, theoretical studies of extended samples with length larger than the atomic transition wavelength are usually restricted to their early-time behavior or to the few-excitation limit. We use herein a mean-field approach to reduce the complex many-body system to an effective two-atom master equation that includes all correlations up to second order and that can be numerically propagated in time. We find that three-dimensional and two-dimensional inverted atomic arrays sustain superradiance below a critical lattice spacing and quantify the scaling of the superradiant peak for both dimensionalities. Finally, we study the late-time dynamics of the system and demonstrate that a subradiant phase appears before the system finally relaxes.
The interaction between light and matter has captivated physicists for centuries, from early studies of vision and refraction in ancient Greece to the development of quantum mechanics and quantum electrodynamics in the past century. While the response of a single quantum emitter to light is well understood, the radiative properties of an ensemble of closely spaced emitters are far more intricate. Coupling to a shared electromagnetic environment induces coherent and dissipative interactions between emitters, giving rise to a collective response that cannot be captured by treating them independently. In the regime of few excitations, the system hosts delocalized subradiant states, that is, coherent superpositions that are largely decoupled from the electromagnetic field and thus decay at suppressed rates. While this weak coupling makes subradiant states attractive for quantum technologies, it also renders them difficult to manipulate. At higher excitation densities, the intrinsic nonlinearity of emitters and the exponential growth of the Hilbert space make theoretical and numerical descriptions of the system and its dynamics increasingly challenging. This thesis explores two fundamental questions: How can subradiant and dark states be selectively accessed and harnessed for practical applications in quantum technologies? And how can interacting ensembles of quantum emitters be efficiently simulated to uncover their many-body physics? The first part of the thesis develops protocols for controlling and addressing dark states in free-space and waveguide-coupled atomic arrays, demonstrating their utility in quantum storage and the deterministic generation of entangled photonic states. Incorporating atomic motion, we further show that collective subradiant states can enhance cooling in dense atomic arrays, offering new avenues for controlling motional dynamics. In the second part, we introduce cumulant expansions of the equations of motion as a powerful tool to analytically and numerically investigate collective decay in the many-body regime. We first examine the collective decay of fully excited atomic arrays in free space, characterizing the onset and scaling of the superradiant burst across different geometries. In collaboration with experiments on ultracold erbium atoms in optical lattices, we provide the first direct observations of many-body collective effects in free-space ordered arrays, including early-time superradiant bursts, late-time subradiant tails, and the emergence of atomic correlations throughout the dynamics. Finally, we theoretically and numerically explore the transient formation of multi-excitation subradiant states, and demonstrate how the existence of multiple dissipation channels suppresses steady-state superradiance in extended arrays.
In quantum optics, superradiance is a phenomenon in which a system of N fully excited quantum emitters radiate intense flashes of light during collective decay. However, computing its peak intensity exactly for many spatially separated emitters remains challenging due to the exponential growth of the underlying Hilbert space with system size N. Based on third-order cumulant expansion methods, we present general scaling laws for the expononent of the peak emission rate as a function of the emitter number in free-space emitter arrays and arrays coupled to one-dimensional waveguide reservoirs. We find, that for 1D chains in free-space the peak emission rate scales linearly with N, while for 2D and 3D arrays with finite emitter spacing it scales superlinearly but sub-quadratically. For emitter chains coupled to waveguide reservoirs we find that the peak emission rate scales quadratically with N.
When quantum emitters couple indistinguishably to light, they can synchronize into a collective light matter system with radiative properties profoundly different from those of independent particles. To date, the resulting collective effects have largely been confined to point like or homogeneous ensembles. Here, we open access to a qualitatively new collective regime by realizing geometrically ordered, spatially extended atom arrays with subwavelength spacing. This establishes a fundamentally new platform in which collective emission is no longer confined to a single Dicke mode but instead emerges from an ordered network of photon mediated interactions. We find that 2D atom arrays undergo strong super and subradiant emission. Despite subwavelength spacing, we achieve site resolved imaging and directly observe the buildup of spatial correlations, demonstrating the transformation of cooperative decay into a strongly correlated many-body process. We observe extensive scaling of superradiance, uncover superradiant revivals, and reveal the ferromagnetic nature of superradiance and the antiferromagnetic nature of subradiance. Our results realize a novel programmable platform for exploring and utilizing dissipative many-body quantum physics, opening new possibilities for photon capture, storage, and atom photon entanglement.
Determining the peak photon emission time and rate for an ensemble of N quantum systems undergoing collective superradiant decay typically requires tracking the time evolution of the density operator. Generally, the dimension of the density operator grows exponentially (∼2N) with the number of emitters, in the absence of any symmetries such as in Dicke superradiance with full or partial permutational symmetry. We present a detailed study of the superradiant peak emission rate and time for initially fully excited quantum emitter ensembles, for one-, two- and three-dimensional arrays in free-space and emitter chains coupled to waveguide reservoirs. For few emitters (N≲14) we utilize the full quantum master equation, and for mesoscopic emitter numbers (N≲400) we use a second- and third-order cumulant expansion of the operator averages to track the time evolution of the system.
Collective superradiant decay of a tightly packed inverted quantum emitter ensemble is among the most intensely studied phenomena in quantum optics. Since the seminal work of Dicke more than half a century ago, a plethora of theoretical calculations in quantum many-body physics have followed. Widespread experimental efforts range from the microwave to the X-ray regime. Nevertheless, accurate calculations of the time dynamics of the superradiant emission pulse still remain a challenging task requiring approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a newly found exact solution. The application of two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The long-term population evolution is only correctly captured for low emitter numbers, where an individual spin-based cumulant expansion proves more reliable than the collective spin-based variant. Surprisingly, odd orders show unphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, on longer time scales the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when used for long-time evolution or reliably finding steady states.