“…Presently nanolasers are theoretically modeled either by rate equations as in [37,38], by numerical solution of the density matrix equations as in [22,[39][40][41], or by systems of equations for correlations as in the cluster expansion [25,42] or cumulant expansion [43,44] methods. Numerical analysis of superradiant emission and lasing has recently led to new and interesting results, such as mechanical effects in photon-atom interactions [45], lasing with a millihertz linewidth and rapid emitter number fluctuations [46], Wigner functions for semiconductor heterostructures [47], transition from superradiance to regular lasing by varying the coherent and incoherent driving [44], sub-and superradiance in multimode optical waveguides [48], and photon-antibunching in the fluorescence from an optical nanofiber-tip [49]. However, complementary analytical methods to model nanolasers without adiabatic elimination of polarization, that would apply to superradiant nanolasers, are not well developed.…”