Generalized Tavis-Cummings Model

• The generalized Tavis-Cummings model is a family of light-matter Hamiltonians that extends the standard two-level emitter framework with collective coupling and symmetry-based structures.
• It incorporates additional interactions such as direct spin-spin coupling, anharmonic multilevel emitters, Kerr-nonlinearity, and multi-mode effects to enrich the excitation manifold and quantum dynamics.
• These extensions enable practical applications in quantum metrology, entanglement control, nonlinear spectroscopy, and digital quantum simulations for exploring many-body quantum phenomena.

The generalized Tavis-Cummings model is a family of light-matter Hamiltonians that extend the standard Tavis-Cummings description of $N$ two-level systems collectively coupled to a single bosonic mode. In the standard resonant, homogeneous form, the model is

$$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$

or, in collective notation, $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$. The generalized models surveyed in recent work preserve the basic cavity-mediated collective coupling while adding direct spin-spin interactions, multilevel anharmonicity, intensity dependence, Kerr terms, multiple bosonic modes, explicit driving, dissipation, disorder, or reduced-subspace constructions. As a result, the phrase “generalized Tavis-Cummings model” refers not to a single canonical modification but to a broad research program organized around symmetry, integrability, excitation-manifold structure, and collective quantum dynamics (0906.4005, Dukalski et al., 2013).

1. Canonical structure and symmetry backbone

The standard Tavis-Cummings model is the multi-emitter generalization of the Jaynes-Cummings model: a set of identical two-level atoms or qubits interacts with a single cavity mode under the rotating-wave approximation. In homogeneous settings, the collective spin operators $J_\pm=\sum_i \sigma_i^\pm$ and $J_z$ expose the permutation symmetry of the emitter sector and reduce the many-qubit problem to collective pseudospin sectors. A central consequence is that the full $2^N$-dimensional qubit Hilbert space can be block-diagonalized into independent subsystems labeled by a pseudospin $j$, with each block behaving as a higher-pseudospin Jaynes-Cummings model; the largest block has size $N+1$ (Dukalski et al., 2013).

This symmetry backbone also governs the excitation-manifold structure. Invariant subspace constructions based on Dicke states and fixed total excitation number allow the dynamics to remain within compact sectors whose dimension scales with $\min\{N_b,n_0\}+1$ rather than exponentially, a fact exploited both analytically and algorithmically in later generalized models (Yang et al., 2023). The homogeneous model additionally exhibits collective dynamical phenomena such as collapse and revival, attractor states, and the migration of qubit entanglement into the field. For multi-qubit states within a “basin of attraction,” the reduced qubit state reaches a universal attractor spin-coherent state at special fractions of the revival time, while states outside that basin can display sudden birth and sudden death of entanglement (0906.4005).

A recurrent feature of generalizations is therefore not the abandonment of the original TC structure, but the controlled modification of its conserved quantities, collective operators, and symmetry sectors. Some extensions retain excitation conservation and simply enrich the matter sector; others break $U(1)$ symmetry, introduce dissipation, or replace exact solvability by effective descriptions.

2. Interacting-spin generalizations

A major line of generalization replaces noninteracting emitters by an interacting spin medium. In “Quantum metrology enhanced by the $$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$0 spin interaction in a generalized Tavis-Cummings model” (Su et al., 2023), the emitters form a one-dimensional spin chain with nearest-neighbor $$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$1-type interaction,

$$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$2

with

$$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$3

Under large detuning, a time-averaged method yields

$$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$4

which directly links metrological sensitivity to collective spin fluctuations. In this model, a finite $$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$5 is necessary and sufficient for Heisenberg-scaling quantum Fisher information. The noninteracting case $$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$6 yields standard-quantum-limit scaling, whereas anisotropic $$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$7 or Ising interactions $$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$8 produce enhanced fluctuations in the ordered phase and permit Heisenberg scaling. The isotropic $$\hat H=\hbar\omega \hat a^\dagger \hat a+\frac{\hbar\omega}{2}\sum_{i=1}^{N_q}\hat\sigma_i^z+\hbar\lambda\sum_{i=1}^{N_q}\left(\hat a\hat\sigma_i^+ + \hat a^\dagger \hat\sigma_i^- \right),$$9 limit $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$0 does not achieve this enhancement, and increasing $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$1 increases the QFI monotonically (Su et al., 2023).

A closely related development maps the generalized $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$2-interacting TC model to a central spin model. In “Dynamical multipartite entanglement in a generalized Tavis-Cummings model with XY spin interaction” (Su et al., 9 May 2025), strong detuning and adiabatic elimination lead to an effective Hamiltonian containing both a $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$3 term and a $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$4 correction. After a Holstein-Primakoff transformation, the cavity mode is represented by a collective central spin, producing

$\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$5

Within this reduced description, strong coupling and low temperature are necessary conditions for genuine multipartite entanglement in the $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$6 model, and the magnetic field modulates both the period and amplitude of the QFI dynamics. The work makes explicit that generalized TC models can serve as bridges between cavity QED and central-spin many-body physics (Su et al., 9 May 2025).

Another interacting-spin extension adds all-to-all qubit interactions and a parametric cavity drive. In “Parametrically Driven Superradiance of an Interacting Tavis-Cummings Model” (Geng et al., 31 Aug 2025), the Hamiltonian includes a squeezing term $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$7 and a uniform interaction term proportional to $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$8. The parametric drive breaks the $\hat{\mathcal V}=g(J_+\hat a+J_-\hat a^\dagger)$9 symmetry of the standard TC model down to $J_\pm=\sum_i \sigma_i^\pm$0, allowing a superradiant transition with squeezed cavity fields. Repulsive interactions shift the phase boundary upward and suppress superradiance; attractive interactions generally enhance it, although re-entrant normal regions (“fingers”) appear at discrete attractive couplings. The low-energy behavior is captured by effective models involving only a few Dicke states, underscoring how generalized TC physics often remains collective even when the microscopic Hamiltonian becomes more elaborate (Geng et al., 31 Aug 2025).

3. Multilevel emitters, anharmonicity, and intensity dependence

A second major generalization abandons the two-level approximation for the matter subsystem. In “Generalization of the Tavis-Cummings model for multi-level anharmonic systems” (Campos-Gonzalez-Angulo et al., 2021), each emitter has $J_\pm=\sum_i \sigma_i^\pm$1 discrete, possibly anharmonic levels, and only adjacent transitions couple to the cavity: $J_\pm=\sum_i \sigma_i^\pm$2 Because the emitters are identical, the spectrum decomposes into symmetry-classified submanifolds associated with irreducible representations of the permutation group. Operationally, the number of emitters becomes only a parameter: the block sizes to be diagonalized are independent of $J_\pm=\sum_i \sigma_i^\pm$3. This formalism exposes multiple classes of dark or subradiant states with no single-molecule counterpart, together with delocalized modes containing more than one excitation per molecule, which have no equivalent in the standard two-level TC model (Campos-Gonzalez-Angulo et al., 2021).

The second excitation manifold was developed in greater detail in “Generalization of the Tavis-Cummings model for multi-level anharmonic systems: insights on the second excitation manifold” (Campos-Gonzalez-Angulo et al., 2022). There, a brute-force $J_\pm=\sum_i \sigma_i^\pm$4 diagonalization is reduced to matrices of at most $J_\pm=\sum_i \sigma_i^\pm$5, and explicit results are given for harmonic and anharmonic emitters. Two forms of anharmonicity are distinguished: mechanical anharmonicity $J_\pm=\sum_i \sigma_i^\pm$6, which shifts level spacings, and electrical anharmonicity $J_\pm=\sum_i \sigma_i^\pm$7, which modifies successive transition dipoles. The analysis identifies resonant conditions between bipolaritons and anharmonic transitions under which two-photon absorption is enhanced. At the same time, the induced polaritonic energy shifts scale away with increasing $J_\pm=\sum_i \sigma_i^\pm$8, so calculations based on a single or a few emitters qualitatively fail to represent the nonlinear optical response of the collective strong-coupling regime (Campos-Gonzalez-Angulo et al., 2022).

A different route to generalized matter structure appears in “Entangling operations in nonlinear two-atom Tavis-Cummings models” (Gómez-Rosas et al., 2021), where the interaction is intensity dependent: $J_\pm=\sum_i \sigma_i^\pm$9 The standard Tavis-Cummings model is recovered by $J_z$0, while $J_z$1 yields the Buck-Sukumar case. For large coherent initial fields, the time-dependent state vector admits an approximate analytical form in terms of atomic Bell states and field coherent states. This permits simple formulas for concurrence and purity at quarter and half revival times, and it supports measurement-induced Bell-state projectors, a complete Bell measurement protocol, and an effective three-qubit gate that can generate GHZ and W states. The same generalized structure arises naturally in an ion-trap implementation where the bosonic mode is the quantized center-of-mass motion rather than a cavity photon field (Gómez-Rosas et al., 2021).

4. Multimode, Kerr, and explicitly time-dependent extensions

Generalization also occurs on the bosonic side. A three-mode trilinear variant,

$J_z$2

was treated algebraically in “Algebraic approach to the Tavis-Cummings model with three modes of oscillation” (Choreño et al., 2017). Bogoliubov transformations and normal-mode operators solve special cases, while $J_z$3 and $J_z$4 tilting transformations diagonalize the most general form. The resulting eigenfunctions combine ordinary oscillator states with Perelomov number coherent states, and the spectra reveal both exact solvability and, in the $J_z$5 realization, a regime of complex energies that signals dynamical instability (Choreño et al., 2017).

The time-dependent three-mode extension was used in “Berry phase of the Tavis-Cummings model with three modes of oscillation” (Choreño et al., 2019), where

$J_z$6

By rewriting the Hamiltonian in $J_z$7 and $J_z$8 forms and applying Lewis-Riesenfeld invariant theory with displacement operators, explicit Berry phases were obtained for the trilinear model and, more generally, for time-dependent Hamiltonians linear in the generators of those groups (Choreño et al., 2019).

Kerr-type nonlinear generalizations remain within a single-mode cavity but alter the field spectrum. In “Dynamic correlation funtions of the generalized Tavis-Cummings model” (Bogoliubov et al., 2017), the Hamiltonian

$J_z$9

is treated by the quantum inverse method. Bethe states, Bethe equations, Slavnov determinants, and determinant expressions for transition elements produce exact dynamic correlation functions. Relative to the $2^N$0 limit, finite Kerr nonlinearity sharpens and localizes the collapse-and-revival structure of observables such as the atomic inversion, while also making the dynamics more irregular as the nonlinear parameter increases (Bogoliubov et al., 2017).

A distinct explicitly time-dependent generalization is the Landau-Zener or driven Tavis-Cummings model,

$2^N$1

where the bosonic mode frequency is swept linearly in time. “Landau-Zener extension of the Tavis-Cummings model: structure of the solution” (Sun et al., 2016) gives exact state-to-state transition probabilities and coarse-grained emitted-boson distributions expressed through $2^N$2-Pochhammer symbols and $2^N$3-binomial coefficients. This places a generalized TC model inside the class of integrable multistate Landau-Zener problems and connects its probabilities to $2^N$4-deformed binomial statistics (Sun et al., 2016).

5. Open-system, disordered, and reduced-complexity formulations

Open-system generalizations replace unitary TC evolution by Lindblad dynamics. In the single-excitation manifold, “Singly-excited resonant open quantum system Tavis-Cummings model with quantum circuit mapping” (Marinkovic et al., 2022) gives an analytical solution for arbitrary $2^N$5 with linear complexity. The associated Q-MARINA mapping uses $2^N$6 qubits, one for each atom plus one for a lossy cavity/environment qubit, and encodes the dynamics in $2^N$7 entangling gates. Because the model is restricted to a single excitation, the circuit depth and qubit count scale linearly with $2^N$8, whereas the classical master-equation description scales exponentially with the full Hilbert space size (Marinkovic et al., 2022).

A more general digital-simulation treatment appears in “Digital Quantum Simulations of the Non-Resonant Open Tavis-Cummings Model” (Sims et al., 30 Jan 2025). There, the model is non-resonant, inhomogeneous, dissipative, optionally driven, and allows up to three cavity excitations. The Hamiltonian includes distinct emitter frequencies $2^N$9, couplings $j$0, and a coherent pump, while the open dynamics is governed by a Lindblad master equation with cavity and emitter loss. Two quantum algorithms are implemented: a Split $j$1-Matrix method with gate complexity $j$2, and a wave matrix Lindbladization method with gate complexity

$j$3

These constructions are benchmarked against classical solvers and are used to access non-resonant, inhomogeneous regimes that the paper describes as classically intractable (Sims et al., 30 Jan 2025).

Disorder supplies another axis of generalization. In “Integrability, multifractality, and two-photon dynamics in disordered Tavis-Cummings models” (Wierzchucka et al., 2023), the Hamiltonian contains both local energy shifts $j$4 and inhomogeneous couplings $j$5. The paper establishes exact solvability, in the two-excitation sector, not only for homogeneous coupling with disordered energies, but also for homogeneous energies with disordered couplings. Three classes of eigenstates are identified—triplet polaritons, doublet polaritons, and dark states—with distinct multifractal and semilocalized properties. This suggests that integrability in generalized TC models can survive certain kinds of disorder, although only in specific excitation sectors or symmetry limits (Wierzchucka et al., 2023).

Open molecular versions show how dissipation couples bright and dark sectors. In “The role of dephasing for dark state coupling in a molecular Tavis-Cummings model” (Davidsson et al., 2023), a CO molecule with vibrational and electronic structure is embedded alongside an ensemble of two-level systems and evolved under a Lindblad equation with pure dephasing and photon decay. In the ideal symmetric Hamiltonian, dark states remain optically decoupled, but dephasing breaks the symmetry and transfers population into the dark-state subspace. The paper reports that both the dephasing rate and the number of particles have a significant influence on the dark-state population, showing that in generalized molecular TC settings the environment is not merely a source of loss but also a mechanism for symmetry breaking and state redistribution (Davidsson et al., 2023).

6. Collective phenomena, applications, and conceptual scope

Generalized Tavis-Cummings models are used to study metrology, multipartite entanglement, nonlinear spectroscopy, quantum information processing, quantum batteries, and scalable simulation. In metrological settings, spin anisotropy in the $j$6-interacting model is the essential resource behind Heisenberg-scaling QFI for weak-field estimation, and the transition from Heisenberg to standard-quantum-limit scaling tracks the underlying many-body phase boundary (Su et al., 2023). In entanglement dynamics, the nonlinear two-atom model yields Bell-state-selective field measurements and three-qubit gate constructions, while the $j$7-interacting many-body model identifies strong coupling and low temperature as necessary for genuine multipartite entanglement (Gómez-Rosas et al., 2021, Su et al., 9 May 2025).

Quantum thermodynamic applications supply a different interpretation of collectivity. In “Optimal energy storage in the Tavis-Cummings quantum battery” (Yang et al., 2023), the TC Hamiltonian is recast as a charger-battery system in which $j$8 atoms interact collectively with an initial $j$9-photon field. An invariant subspace method reduces the problem to a tridiagonal matrix of dimension $N+1$0. In the limits $N+1$1 or $N+1$2, an emergent $N+1$3 symmetry appears, yielding optimal energy storage at a finite charging time. The paper also establishes a negative relationship between battery-charger entanglement and stored energy, and identifies $N+1$4 as an asymptotically optimal regime (Yang et al., 2023).

Nonlinear optical and spectroscopic consequences are equally prominent. Multilevel anharmonic generalizations produce bipolaritonic resonances and two-photon absorption channels absent from the two-level TC model, while disordered two-excitation models display boson bunching and a two-photon blockade tied directly to the localization properties of the eigenstates (Campos-Gonzalez-Angulo et al., 2022, Wierzchucka et al., 2023). The multi-level group-theoretic framework was developed specifically with nonlinear spectroscopy and polaritonic chemistry in view, where dark submanifolds, multi-excitation delocalization, and weakly $N+1$5-dependent anharmonic shifts all matter for interpretation (Campos-Gonzalez-Angulo et al., 2021).

A consistent theme across these works is that performance enhancements are conditional rather than generic. Superradiance requires symmetry breaking beyond the standard $N+1$6 TC conservation law in the parametrically driven case; Heisenberg scaling requires finite spin fluctuations generated by anisotropy; genuine multipartite entanglement requires strong coupling and low temperature; and exact solvability, when available, is usually tied to special excitation manifolds, collective sectors, or algebraic structures. The generalized Tavis-Cummings model is therefore best understood as a structured landscape of extensions of the TC paradigm, unified by collective light-matter coupling but differentiated by which symmetries are preserved, broken, or replaced (Geng et al., 31 Aug 2025, Sun et al., 2016, Wierzchucka et al., 2023).