This work presents the development of a computational framework for simulating superfluorescence dynamics in a wide range of materials, distinguished by the explicit and exact treatment of emitter-emitter interactions and reduced reliance on conventional approximation methods, allowing for the capture of cooperative emission phenomena with high fidelity in a variety of regimes, such as high excitation density, non-classical initial system states and strong dipole–dipole coupling. This enables the further study of material advancements in the pursuit of optimizing coherent photon emission, as well as in the converse photovoltaic superabsorption.

This model stands out from the other state-of-the-art models by approaching the issue of collective emission phenomena with an ab-initio methodology, avoiding any unnecessary approximation methods such as the conventionally used mean-field approach^1,2, material-related restrictive approximations³ or derived phenomenological equations of motion for the emitter system⁴. By treating the problem explicitly, the system is limited to a rather small size, but enables the computations of all pairwise and higher-order coupling interactions, including interactions with the radiation field, thus allowing one to derive the quantum phenomena of collective spontaneous emission directly from the ab-initio dynamics of the system.

The model employs a Markovian time-evolution algorithm that iteratively updates the density matrixes of the emitter systems, allowing emergent coherence and phase synchronization to arise naturally from first principles. This novel approach is in direct comparison with several recent models which utilize equations of motion^1,4 which already include parameters for cooperative effects based on approximations such as the mean-field approximation. By minimizing said approximations, our framework avoids common sources of systematic error and reveals the true nature of a coherence buildup which is often assumed and simplified in other models.

The model demonstrates strong robustness across a wide parameter space, including variations in emitter geometry, initial system conditions, and excitation conditions. The resulting values from the model can be directly compared with experimentally accessible observables⁵, such as emission intensity profiles and second-order correlation values.

While the model allows for the user-defined arrangement of emitters, there is currently progress being made in creating unit cell structures simulating perovskite quantum dots (QD), which will enable the study of superfluorescence dynamics and coherence buildup in perovskite QD systems. This facilitates straightforward experimental validation using existing measurement techniques in perovskite systems, and may also be used to advance our understanding of superabsorption and its role in enhancing the efficiency of perovskite QD photovoltaics⁶.

Overall, the presented approach establishes a rigorous, extensible and developing platform for studying collective light–matter interactions in perovskite systems based on exact numerical dynamics, offering both enhanced predictive capability and increased assistance to experimental investigations of superfluorescent and superabsorption phenomena, helping to enable advancements in both the quantum optics and photovoltaics fields respectively.