Benchmarking the Cumulant Expansion Method Using Dicke Superradiance

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.