A mathematical preprint predicts that emission from a strongly coupled light–matter system moves toward the energy of its cavity as the excitation count rises. At high excitation, the calculated emission becomes increasingly cavity-like rather than staying near the lower-polariton energy used as the model’s single-excitation reference.
In the calculation, emission energies in higher-m sectors move near the cavity once the excitation count exceeds half the emitter count. After the count reaches roughly three-quarters of the emitter count, all calculated energies overtake the single-excitation lower-polariton energy; at still higher excitation, they collapse toward the cavity energy.
A symmetry-based route through the calculation
The work is a modeling study of the Tavis–Cummings model, which describes two-level quantum emitters interacting with one electromagnetic cavity mode. Under the rotating-wave approximation, total excitation number is conserved, separating the model’s states into finite-dimensional excitation manifolds.
The main computational simplification comes from permutation symmetry. It divides the Hilbert space into irreducible symmetry sectors and removes direct dependence on the number of emitters from the computational workload.
The paper derives analytic formulas for the dimensions and multiplicities of those sectors. Their reduced interaction matrices are real, symmetric and tridiagonal, with simple real eigenvalues arranged symmetrically around zero.
An illustrative calculation used 1,000 emitters and excitation counts of 150, 300, 400 and 500. As the count rose, dark-polariton multiplicities became comparable to and then exceeded dark-state multiplicities.
The framework also calculates each symmetry sector’s average photonic and two-level-system excitation content. These quantities are basis-independent and have exact formulas within the model.
What the model predicts for emission
The model’s radiative transitions are constrained by symmetry. Cavity and collective-spin operators preserve the symmetry sector unless symmetry-breaking processes act. In the large-emitter-limit calculation, the lowest-multipolariton’s photonic content equals its squared transition probability, with both evaluated as X/2 in the paper’s notation.
For slow dephasing, the paper places predicted emission peaks at energy differences between the lowest-energy states of the relevant symmetry sectors. It weights those peaks by the sectors’ multiplicities, using that distribution as a model-based proxy for the density of states.
A separate fast-internal-relaxation calculation for 1,000 emitters predicts a sublinear blue shift—an upward move in energy—in the lowest-multipolariton emission. The shift approaches the cavity energy after the excitation count exceeds the emitter count.
An idealized prediction, not an experiment
The analysis is restricted to zero detuning, uniform coupling and a single cavity mode. The authors propose broader regimes as extensions of the present treatment.
The emission calculations use simplified dephasing and relaxation scenarios, so competing processes could affect whether the predicted peaks are observable. The multiplicity-based weights are model-derived density-of-states proxies, not complete measured intensities.
The supplied document is an arXiv version-1 preprint dated 20 August 2026, and no experimental validation is reported in the supplied analysis. The authors state that code for generating the data and figures is available in a GitHub repository.
Paper data and sources
Original title: Group-theoretic treatment of strong light-matter coupling with an arbitrary number of excitations
Authors: Antti Peltola, Olli Siltanen, Kimmo Luoma, Konstantinos S. Daskalakis
Journal/Repository: arXiv
Status: Preprint, not yet peer-reviewed
First online: 2026-08-20
DOI: Not available
Original paper · Full text