The Tavis-Cummings model (the Dicke model treated in the rotating-wave approximation), describing many two-level systems coupled to a single bosonic mode, has been long known to show collective semiclassical oscillations when prepared in an inverted state, with all two-level systems excited and the bosonic mode is empty. This paper discusses how the quantum dynamics approaches this semiclassical result for large numbers of two-level systems, focusing on how the eigenvalues approach their semiclassical limit. The approach to the semiclassical result is found to be slow, scaling like a power of the logarithm of the system size. Considering also the effect of weak detuning between the two-level system and the bosonic field, quantum corrections are again found to decay slowly with system size, such that for a fixed detuning, the quantum effects of detuning are greater than the classical effect.
The Tavis-Cummings model (the Dicke model treated in the rotating-wave approximation), describing many two-level systems coupled to a single bosonic mode, has been long known to show collective semiclassical oscillations when prepared in an inverted state, with all two-level systems excited and the bosonic mode is empty. This paper discusses how the quantum dynamics approaches this semiclassical result for large numbers of two-level systems, focusing on how the eigenvalues approach their semiclassical limit. The approach to the semiclassical result is found to be slow, scaling like a power of the logarithm of the system size. Considering also the effect of weak detuning between the two-level system and the bosonic field, quantum corrections are again found to decay slowly with system size, such that for a fixed detuning, the quantum effects of detuning are greater than the classical effect.
The Tavis-Cummings model for more than one qubit interacting with a common oscillator mode is extended beyond the rotating wave approximation (RWA). We explore the parameter regime in which the frequencies of the qubits are much smaller than the oscillator frequency and the coupling strength is allowed to be ultrastrong. The application of the adiabatic approximation introduced by Irish et al. [Phys. Rev. B 72, 195410 (2005)] for a single-qubit system is extended to the multiqubit case. For a two-qubit system, we identify three-state manifolds of close-lying dressed energy levels and obtain results for the dynamics of intramanifold transitions that are incompatible with results from the familiar regime of the RWA. We exhibit features of two-qubit dynamics that are different from the single-qubit case, including calculations of qubit-qubit entanglement. Both number state and coherent state preparations are considered, and we derive analytical formulas that simplify the interpretation of numerical calculations. Expressions for individual collapse and revival signals of both population and entanglement are derived.
Confined electromagnetic modes strongly couple to collective excitations in ensembles of quantum emitters, producing light–matter hybrid states known as polaritons. Under such conditions, the discrete multilevel spectrum of molecular systems offers an appealing playground for exploring multiphoton processes. This work contrasts predictions from the Tavis–Cummings model in which the material is a collection of two-level systems, with the implications of considering additional energy levels with harmonic and anharmonic structures. We discuss the exact eigenspectrum, up to the second excitation manifold, of an arbitrary number N of oscillators collectively coupled to a single cavity mode in the rotating-wave approximation. Elaborating on our group-theoretic approach [New J. Phys. 23, 063081 (2021)], we simplify the brute-force diagonalization of N2 × N2 Hamiltonians to the eigendecomposition of, at most, 4 × 4 matrices for arbitrary N. We thoroughly discuss the eigenstates and the consequences of weak and strong anharmonicities. Furthermore, we find resonant conditions between bipolaritons and anharmonic transitions where two-photon absorption can be enhanced. Finally, we conclude that energy shifts in the polaritonic states induced by anharmonicities become negligible for large N. Thus, calculations with a single or few quantum emitters qualitatively fail to represent the nonlinear optical response of the collective strong coupling regime. Our work highlights the rich physics of multilevel anharmonic systems coupled to cavities absent in standard models of quantum optics. We also provide concise tabulated expressions for eigenfrequencies and transition amplitudes, which should serve as a reference for future spectroscopic studies of molecular polaritons.
An ensemble of N two-level molecules prepared in its ground state and sharing a lossless cavity with a weak field does not simply absorb: the intensity collapses and then recurs, in a train of collective photon echoes. Working from the exact solution of the Tavis-Cummings model in the few-photon regime nˉ<N, we confirm the echo time τE=4πN+Δ/g numerically at resonance from N=5 to 400, the small-N end discriminating this form from the alternative 4πN−nˉ+Δ/g in its favor. The initial state requires no preparation, being the ground state. The echo time is independent of the initial photon distribution: coherent, thermal, squeezed and oscillatory-squeezed distributions spanning variances from 2 to 34 all recur together, as does a controlled pair with identical mean and variance differing only in the shape of ρnn. The echo amplitude is not: it varies at the tens-of-percent level at fixed mean, including a factor of 1.8 with the squeezing phase at fixed squeezing strength. Detuning organizes the dynamics into two clean regimes separated by a fragmented crossover, the dispersive-branch echo time approaching one-half the resonant one, and sufficient detuning removes the dependence on the initial Dicke state. For arbitrary initial Dicke state, emission replaces absorption at m≃−N/2+nˉ, and the echo envelope acquires one component per step up the ladder: a single-component echo occurs only at the ground state. A feasibility analysis against a five-qubit superconducting device, with Lindblad simulations of cavity decay and dephasing and full-Hilbert-space disorder simulations, shows the first echoes observable at N∼5-20 on existing hardware: the echo survives the dominant loss channel with contrast ∼e−κτE/2, photons being shielded from cavity decay while resident in the emitters.
Resonant interaction of a collection of two-level atoms with a single-mode coherent cavity field is considered in the framework of the Dicke model. We focus on the role of collective atomic effects in the phenomenon of collapses and revivals of the Rabi oscillations. It is shown that the behavior of the system strongly depends on the initial atomic state. In the case of the initial half-excited Dicke state we account for a number of interesting phenomena. The correlations between the atoms result in a suppression of the revival amplitude, and the revival time is halved, compared to the uncorrelated fully excited and ground states. The phenomenon of squeezing of the radiation field in the atom-field interaction is also discussed. For the initial fully excited and ground atomic states, the field is squeezed on the short-time scale, and squeezing can be enhanced by increasing the number of atoms. Some empirical formulas are found which describe the behavior of the system in excellent agreement with numerical results. For the half-excited Dicke state, the field can be strongly squeezed on the long-time scale in the case of two atoms. This kind of squeezing is enhanced by increasing the intensity of the initial coherent field and is of the same nature as revival-time squeezing in the Jaynes-Cummings model. The appearance of this long-time squeezing can be explained using the factorization approximation for semiclassical atomic states.
The Tavis-Cummings model is a paradigmatic central-mode model in which a set of two-level quantum emitters (spins) are coupled to a collective cavity mode. Here we study the eigenstate spectrum, its localization properties, and the effect on dynamics, focusing on the two-excitation sector relevant for nonlinear photonics. These models admit two sources of disorder: in the coupling between the spins and the cavity, and in the energy shifts of the individual spins. While this model was known to be exactly solvable in the limit of a homogeneous coupling and inhomogeneous energy shifts, we establish here the solvability in the opposite limit of a homogeneous energy shift and inhomogeneous coupling, presenting the exact solution and corresponding conserved quantities. We identify three different classes of eigenstates, exhibiting different degrees of multifractality and semilocalization closely tied to the integrable points, and we study their stability to perturbations away from these solvable points. The dynamics of the cavity occupation number away from equilibrium, exhibiting boson bunching and a two-photon blockade, is explicitly related to the localization properties of the eigenstates, and it illustrates how these models support a collective spin description despite the presence of disorder. Published by the American Physical Society 2024
We study the dynamics of two qubits interacting with a single mode of a harmonic oscillator beyond the rotating wave approximation in the ideally degenerate regime. Exact analytic expressions are obtained for state properties of interest, including qubit entanglement for a certain class of initial states of the oscillator and the qubits. Qualitative differences and similarities in the evolution of the qubits in the degenerate regime when the oscillator is treated quantum mechanically and classically are discussed.