Thermal state structure in the Tavis--Cummings model and rapid simulations in mesoscopic quantum ensembles

TL;DR

This paper investigates the thermal behavior of the Tavis-Cummings model in mesoscopic quantum ensembles, revealing a lower temperature transition where the Dicke subspace approximation becomes invalid, and introduces efficient methods for analyzing such systems.

Contribution

It demonstrates the limitations of the Dicke subspace approximation at certain temperatures and develops scalable perturbative methods for thermal property calculations in the Tavis-Cummings model.

Findings

Dicke subspace approximation fails above a certain temperature.

Efficient $O(\sqrt{N})$ methods for thermal property estimation.

Transition temperature lower than previously thought.

Abstract

Hybrid quantum systems consisting of a collection of N spin-1/2 particles uniformly interacting with an electromagnetic field, such as one confined in a cavity, are important for the development of quantum information processors and will be useful for metrology, as well as tests of collective behavior. Such systems are often modeled by the Tavis-Cummings model and having an accurate understanding of the thermal behaviors of this system is needed to understand the behavior of them in realistic environments. We quantitatively show in this work that the Dicke subspace approximation is at times invoked too readily, in specific we show that there is a temperature above which the degeneracies in the system become dominant and the Dicke subspace is minimally populated. This transition occurs at a lower temperature than priorly considered. When in such a temperature regime, the key constants of the motion are the total excitation count between the spin system and cavity and the collective angular momentum of the spin system. These enable perturbative expansions for thermal properties in terms of the energy shifts of dressed states, called Lamb shifts herein. These enable efficient numeric methods for obtaining certain parameters that scale as $O(\sqrt{N})$, and is thus highly efficient. These provide methods for approximating, and bounding, properties of these systems as well as characterizing the dominant population regions, including under perturbative noise. In the regime of stronger spin-spin coupling the perturbations outweigh the expansion series terms and inefficient methods likely are needed to be employed, removing the computational efficiency of simulating such systems. The results in this work can also be used for related systems such as coupled-cavity arrays, cavity mediated coupling of collective spin ensembles, and collective spin systems.