Monte-Carlo simulations of superradiant lasing

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.