We simulate the superradiant dynamics of ensembles of atoms in the presence of collective and individual atomic decay processes. We apply the Monte-Carlo wave-function method and identify quantum jumps in a reduced Dicke state basis, which reflects the permutation symmetry of the system. While the number of density matrix elements in the Dicke representation increases polynomially with atom number, the quantum jump dynamics populates only a single Dicke state at the time and thus efficient simulations can be carried out for tens of thousands of atoms. The superradiant pulses from initially excited atoms agree quantitatively with recent experimental results of strontium atoms but rapid atom loss in these experiments does not permit steady-state superradiance. By introducing an incident flux of new atoms to maintain a large average atom number, our theoretical calculations predict lasing with a millihertz linewidth despite rapid atom number fluctuations.
We simulate the superradiant dynamics of ensembles of atoms in the presence of collective and individual atomic decay processes. We apply the Monte-Carlo wave-function method and identify quantum jumps in a reduced Dicke state basis, which reflects the permutation symmetry of the system. While the number of density matrix elements in the Dicke representation increases polynomially with atom number, the quantum jump dynamics populates only a single Dicke state at the time and thus efficient simulations can be carried out for tens of thousands of atoms. The superradiant pulses from initially excited atoms agree quantitatively with recent experimental results of strontium atoms but rapid atom loss in these experiments does not permit steady-state superradiance. By introducing an incident flux of new atoms to maintain a large average atom number, our theoretical calculations predict lasing with a millihertz linewidth despite rapid atom number fluctuations.
Collective emission behavior is usually described by the decay dynamics of the completely symmetric Dicke states. To study a more realistic scenario, we investigate alternative initial states inducing a more complex time evolution. Superposition states of the fully inverted Dicke state and the Dicke ground state with unequal mutual weights are studied as examples as well as superradiance stemming from atoms in clusters separated by more than one wavelength. The Monte Carlo wave function method serves as a framework to study the dynamics of quantum states, which is determined by quantum jumps, on the one hand, and continuous evolution dynamics, on the other hand. We compare this method with the classical picture of a system of rate equations written for the diagonal components of the density matrix.
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.
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 needing approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a new, recently found exact solution. Applying two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The longterm 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 even qualitatively nonphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, at a longer time scale, the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when applied for long-time evolution or reliably finding steady states.
Superradiance is an extensive cooperative spontaneous emission phenomenon, exhibited by some atomic collective-spin states. However, distinct initial states differ in their decay dynamics. Dicke states ∣j,m⟩ with distinct numbers of excitations n=m+j, driven by vacuum fluctuations, have their peak emission intensity shifted in time. Rotating Dicke states relative to the decay axis introduces an interesting parity structure that affects the pulse profile and photon correlations. Squeezed-bath prepared states undergo a squeezing-controlled crossover to the rotated Dicke states, making the emission character dependent on the amount of squeezing transferred from light to the atomic state. For semiclassical states with a macroscopic dipole moment, like the atomic coherent state, the emission intensity depends on their polarization. We present detailed results on the superradiant dynamics of a representative selection of states expanded in Dicke states to highlight the initial state as an independent dynamical control parameter. For large-N, we are able to predict fairly accurately the pulse profile in each case using the mean-field approximation, an approach based on the Fokker--Planck equation. We also present comparative results on the intensity correlation function, quantify the coherent and incoherent contributions of the emission, and contrast between small and large ensembles.
We report a time-resolved study of collective emission in dense ensembles of two-level atoms. We compare, on the same sample, the buildup of superradiance and subradiance from the ensemble when driven by a strong laser. This allows us to measure the dynamics of the population of superradiant and subradiant states as a function of time. In particular, we demonstrate the buildup in time of subradiant states through superradiant dynamics. This illustrates the dynamics of the many-body density matrix of superradiant ensembles of two-level atoms when departing from the ideal conditions of Dicke superradiance, in which symmetry forbids the population of subradiant states.
Steady-state superradiance and superradiant lasing attract significant attentions in the field of optical lattice clocks, but have not been achieved so far due to the technical challenges and atom loss problem. In this article, we propose that their counter-part may be observed in the microwave domain with solid-state spins-microwave resonator systems at room temperature with realistic technical restrictions. To validate our proposal, we investigate systematically the system dynamics and steady-state by solving quantum master equations for the multi-level and multi-process dynamic of trillions of spins. To this end, we employ a mean-field approach, and convert the mean-field dynamics of the spin ensemble into the one in a more intuitive Dicke state picture. Our calculations show that for systems with nitrogen vacancy center spins and pentacene molecular spins the superradiant Rabi oscillations occur firstly due to transitions among different Dicke states, and the subsequent continuous-wave superradiant masing can achieve a linewidth well below millihertz. Our work may guide further exploration of transient and steady-state superradiant masing with the mentioned and other solid-state spins systems, such as silicon vacancy centers in silicon carbide and boron vacancy centers in hexagonal boron nitride, where the coherent radiation with ultra-narrow linewidth may find applications in deep-space communications, radio astronomy and high-precision metrology.