---
title: Extended Jaynes-Cummings Model (eJCM) Overview
url: https://www.emergentmind.com/topics/extended-jaynes-cummings-model-ejcm
type: topic
---

# Extended Jaynes-Cummings Model (eJCM) Overview

Searching arXiv for the cited papers and related uses of “extended Jaynes-Cummings model.”
arXiv search query: "extended Jaynes-Cummings model"
The extended Jaynes–Cummings model (eJCM) denotes, across the cited literature, a family of non-equivalent extensions of the Jaynes–Cummings Hamiltonian used to describe light–matter coupling beyond the textbook single-mode rotating-wave setting. In this usage, “extended” may refer to inclusion of counter-rotating terms, nonlinear and multi-boson exchange, squeezed or displaced bosonic modes, coherent driving, collective atomic modes, traveling-wave pulses, continuum-coupled autoionizing states, or lattice generalizations of Jaynes–Cummings–Hubbard type. A unifying feature is that the two-level subsystem remains coupled to bosonic or bosonized degrees of freedom, while the spectral, dynamical, or many-body structure departs from the standard JCM in a controlled way [1011.3280], [2010.13867], [1808.06727], [2405.15876], [2010.01940].

## 1. Exact non-RWA formulation

A canonical meaning of eJCM is the Jaynes–Cummings model without the rotating-wave approximation. In laboratory variables the Hamiltonian is
\[
H_0=\frac{\Delta}{2}\sigma_z+\omega a^\dagger a+g(a+a^\dagger)\sigma_x,
\]
and after a \(90^\circ\) rotation about \(y\), with energies rescaled by \(\omega\), it becomes
\[
H=-\frac{\Delta}{2}\sigma_x+a^\dagger a+g(a+a^\dagger)\sigma_z.
\]
This Hamiltonian commutes with the parity operator
\[
\Pi=\sigma_x e^{i\pi a^\dagger a},
\]
so the Hilbert space splits into even and odd parity sectors. The parity decomposition is central because counter-rotating terms mix all photon numbers, whereas parity remains conserved [1011.3280].

Chen et al. solve this non-RWA problem by an extended bosonic coherent-state ansatz built from two coherent states of opposite displacements \(\pm\alpha\), each dressed by a polynomial in \(a^\dagger\) of degree \(M\). After projection onto \((a^\dagger)^m e^{\pm \alpha a^\dagger}\lvert 0\rangle\), the stationary Schrödinger equation reduces to a three-term recurrence for coefficients \(c_m\), with \(c_0=1\) and \(c_1=0\). The truncation condition \(c_{M+1}=0\) yields a single-variable polynomial equation of degree \(M\) in \(\alpha\), and each real root \(\alpha_k\) gives an eigenenergy
\[
E_k^\pm=\alpha_k g \mp \frac{\Delta}{2}.
\]
The associated recurrence determines all coefficients \(c_m(\alpha_k)\), and hence the full eigenfunction. In this construction, the infinite-dimensional bosonic problem is mapped exactly onto a one-variable secular equation whose real roots determine the complete spectrum for arbitrary coupling \(g\) and detuning \(\Delta\) [1011.3280].

The method is both analytic and numerical. For given \((g,\Delta)\), one increases \(M\) until the real roots stabilize; a relative tolerance \(10^{-8}\) is explicitly cited as a convergence criterion. Typical \(M\approx 40\ldots 60\) suffices to converge the first \(\sim 20\) levels even at ultra-strong coupling \(g/\omega\sim 1\), and in the ultra-strong regime the required truncation can fall to \(M\lesssim 20\), whereas direct Fock-space diagonalization would require thousands of photon states for \(g\gtrsim 1\). The construction was motivated in part by circuit-QED systems reaching \(g/\omega\gtrsim 0.1\) and even \(0.5\ldots 1\), where Bloch–Siegert shifts, squeezing, and ground-state entanglement become nonperturbative [1011.3280].

## 2. Systematic approximations beyond the rotating-wave approximation

A complementary exact formulation uses an extended Swain ansatz truncated directly in the Fock basis,
\[
\lvert\psi;\pm\rangle=
\begin{pmatrix}
\sum_{n=0}^M c_n\lvert n\rangle\\[2mm]
\pm\sum_{n=0}^M (-1)^n c_n\lvert n\rangle
\end{pmatrix},
\]
again exploiting the conserved parity \(\Pi=\sigma_x e^{i\pi a^\dagger a}\). Insertion into the Schrödinger equation yields the recursion
\[
g\sqrt{m+1}\,c_{m+1}+\bigl[m\mp (\Delta/2)(-1)^m-E\bigr]c_m+g\sqrt{m}\,c_{m-1}=0,
\]
with \(c_{-1}=0\), \(c_0=1\), and boundary condition \(c_{M+1}=0\). As \(M\to\infty\), the roots of the resulting degree-\(M\) polynomial in \(E\) converge to the exact spectrum [1203.2410].

This formulation supports a hierarchy of analytic approximations. The zero-order approximation gives
\[
E_m^{(0)}=m\mp \frac{\Delta}{2}(-1)^m,
\]
with ground-state energy \(E_{GS}^{(0)}=-\Delta/2\), and contains no counter-rotating-term effects. The first-order \(2\times 2\) block approximation reproduces exactly the standard RWA spectrum,
\[
E_{n,\pm}^{(1)}=n+\frac12\pm \frac12\sqrt{\delta^2+4g^2(n+1)},
\qquad \delta\equiv \Delta-1.
\]
The second-order \(3\times 3\) block approximation produces the first counter-rotating corrections: at resonance, the \(O(g^2)\) shifts and three-state admixtures are identified explicitly as the first signatures of counter-rotating terms. The second-order energies remain accurate with errors \(\lesssim 10^{-3}\) up to \(g\approx 0.2\,\omega\), and therefore cover the \(g/\omega\sim 0.1\) regime cited for modern circuit QED [1203.2410].

The same second-order treatment modifies vacuum Rabi splitting. Under the RWA, the probe operator excites the ground state into an equal-weight superposition of \(\lvert 1,\pm\rangle\), producing two peaks of equal height separated by \(2g\). With counter-rotating terms, the dressed ground state acquires \(\lvert g,2\rangle\) and \(\lvert e,1\rangle\) admixtures, the peak frequencies shift to
\[
\Delta\omega_{1,2}\simeq \pm g \mp \frac14 g^3+O(g^4),
\]
and the total splitting becomes
\[
\Delta\omega_{\mathrm{Rabi}}=2g-\frac14 g^3+\cdots.
\]
The peak heights are unequal, with \(h_1-h_2=O(g)\); at \(g=0.1\,\omega\), the splitting shift is approximately \(0.0025\,\omega\), whereas the height ratio \(h_1/h_2\approx 1.11\) is predicted to be experimentally visible [1203.2410].

## 3. Generalized rotating-wave extensions

Beyond the non-RWA single-mode problem, the label eJCM is also applied to several generalized Hamiltonians that retain the rotating-wave approximation but enlarge the interaction structure. One broad formulation includes Stark shifts, Kerr-type nonlinearities, multi-boson exchange, and intensity-dependent coupling:
\[
H_{eJCM}
=\omega \hat n+\frac{\omega_0}{2}\sigma_z+\sigma_z F(\hat n)+G(\hat n)
+g\bigl[\sigma_+ a^k f(\hat n)+\sigma_- f(\hat n)a^{\dagger k}\bigr].
\]
Here \(F(\hat n)\) generates qubit level shifts, \(G(\hat n)\) nonlinear boson dispersion, \(k\in\mathbb N\) is the number of exchanged bosons, and \(f(\hat n)\) makes the coupling intensity dependent. The model possesses a graded Lie algebra built from nilpotent supercharges \(\mathcal Q,\mathcal Q^\dagger\), a SUSY Hamiltonian \(\mathcal H=\{\mathcal Q^\dagger,\mathcal Q\}\), and excitation operators \(\mathcal B,\mathcal N\). A unitary transformation \(D(\beta)\) diagonalizes the effective Hamiltonian sector by sector, yielding a generalized Rabi frequency \(\Omega(\mathcal N)\), closed-form eigenstates, and a closed-form propagator. Standard JC, Kerr, Buck–Sukumar, Stark-shifted two-photon, and algebraically deformed models appear as special cases [2010.13867].

A different generalized eJCM mixes Jaynes–Cummings and anti-Jaynes–Cummings terms linearly,
\[
H_{\rm eJCM}
=\hbar\bigl[\sigma_-(g^*a+f a^\dagger)+\sigma_+(g a^\dagger+f^* a)\bigr]+mc^2\sigma_z.
\]
Using the \(su(1,1)\) generators \(K_0,K_\pm\), squaring the Hamiltonian leads to a Klein–Gordon-type operator that is diagonalized by a squeezing, or tilting, transformation \(D(\xi)=\exp(\xi K_+-\xi^*K_-)\). The exact spectrum is
\[
E_{n,\pm}=\pm\sqrt{\hbar^2(|g|^2-|f|^2)n+m^2c^4},
\]
and the eigenfunctions are squeezed number states. The same model is connected in the source to the relativistic parametric amplifier and to a single trapped ion driven simultaneously on red and blue motional sidebands [1411.1974].

Related modified models are obtained by replacing the cavity operators by squeezed or squeezed-coherent operators. One construction introduces
\[
b=S(\zeta)D(\alpha)\,a\,D^\dagger(\alpha)S^\dagger(\zeta)
=\cosh r\,a+e^{i\theta}\sinh r\,a^\dagger-\alpha,
\]
so that the transformed Hamiltonian contains two-photon field terms \(b^2,b^{\dagger 2}\), counter-rotating atom–field terms, and an atomic-exchange term proportional to \(\alpha\). Under the RWA and with \(\alpha=0\), the interaction reduces to an effective coupling \(g_{\rm eff}=g\cosh r\), with dressed-state blocks identical in form to the ordinary JCM but sampled over altered squeezed-coherent photon statistics. A closely related squeezed-photon model writes the interaction in terms of \(B=S(r)aS^\dagger(r)=\cosh r\,a+\sinh r\,a^\dagger\), thereby interpolating between the JCM at \(r=0\) and the Rabi Hamiltonian in the \(r\to\infty\) limit [2212.12021], [2405.15876].

| Variant | Defining modification | Source |
|---|---|---|
| Non-RWA single-mode eJCM | Counter-rotating term and parity-resolved exact polynomial | [1011.3280], [1203.2410] |
| Nonlinear generalized eJCM | \(F(\hat n)\), \(G(\hat n)\), \(f(\hat n)\), \(k\)-boson exchange | [2010.13867] |
| JC/anti-JC mixture | Simultaneous \(a\) and \(a^\dagger\) couplings with \(g,f\) | [1411.1974] |
| Squeezed or displaced exchange | \(b=S D\,a\,D^\dagger S^\dagger\) or \(B=S(r)aS^\dagger(r)\) | [2212.12021], [2405.15876] |

## 4. Relativistic correspondences and analog mappings

Several eJCMs admit exact reformulations in relativistic language. The driven Jaynes–Cummings–Rabi Hamiltonian,
\[
H=\omega a^\dagger a+\frac{\omega_0}{2}\sigma_z+g(\sigma_++\sigma_-)(a+a^\dagger)+\epsilon(a+a^\dagger),
\]
is mapped to a massless charged Dirac particle in constant homogeneous electric and magnetic fields. In the Dirac picture, the field invariant \(F^2=B^2-E^2\) separates two regimes. For \(B>E\), corresponding to subcritical drive, the quasienergies are discrete,
\[
\varepsilon_{n,\pm}=\pm \hbar\lambda\,[1-(2\epsilon/\lambda)^2]^{3/4}\sqrt{n+1},
\]
with \(\lambda=g\), and the spectrum collapses as \(\epsilon\to g/2\). For \(E>B\), corresponding to supercritical drive, the spectrum becomes continuous. The critical amplitude is therefore \(\epsilon_c=g/2\). In the generalized case with counter-rotating coupling ratio \(\eta=g'/g\), a squeezing transform leads to \(\epsilon_c(\eta)=\lambda'(1+\eta)/2\simeq g(1+\eta)/2\). Below threshold, the eigenstates are Hermite–Gaussian spinors; above threshold, Weber-function spinors appear [1808.06727].

The generalized linear \(g,f\) Hamiltonian discussed above is also presented as a relativistic parametric amplifier. In its nonrelativistic limit, the model reproduces the degenerate parametric-amplifier Hamiltonian
\[
H_{\rm NR}
=\hbar\omega\Bigl(a^\dagger a+\frac12\Bigr)
-\frac{\hbar\chi}{2}\bigl(a^2 e^{-i\phi}+a^{\dagger 2}e^{i\phi}\bigr).
\]
The same source gives an exact trapped-ion realization under Lamb–Dicke and rotating-wave approximations, where simultaneous red- and blue-sideband driving reproduces the full eJCM coupling structure [1411.1974].

A further mapping relates the \(2+1\) Dirac–Moshinsky oscillator coupled to an external isospin field to a two-atom single-mode eJCM,
\[
H_{\rm eJCM}
=n(\sigma_+ a+\sigma_- a^\dagger)+x(\tau_+ a+\tau_- a^\dagger)
+\Delta_\sigma\sigma_z+\Delta_\tau\tau_z.
\]
The conserved quantity \(I=a^\dagger a+\tfrac12(\sigma_z+\tau_z)\) partitions the Hilbert space into invariant manifolds. Closed-form time-dependent states, reduced density matrices, concurrence or von Neumann entropy, and population inversion are given for both initial number states and initial coherent states. In the abstracted conclusion of that study, the coherent state is reported to give a better description than the number state for entanglement and population inversion [1812.00261].

## 5. Collective, traveling-wave, continuum, and lattice realizations

The eJCM can emerge as an effective low-energy theory for collective matter degrees of freedom. In a system consisting of a single control atom, a cavity, and a sample of \(N\) identical two-level atoms, adiabatic elimination in the dispersive regime \(\Delta_c\gg g\sqrt{N(\bar n+1)}\), \(\Delta_L\gg \Omega\) yields a collective bosonic mode
\[
b=\frac1{\sqrt N}\sum_{j=1}^N S_j^-,
\qquad [b,b^\dagger]\simeq 1
\quad\text{when}\quad N\gg \langle n_b\rangle.
\]
The effective Hamiltonian takes the Jaynes–Cummings form
\[
H_{\rm eJCM}
=\hbar\omega_b\,b^\dagger b+\frac{\hbar}{2}\Delta_{\rm eff}\sigma^z
+\hbar g_{\rm eff}(\sigma^+ b+\sigma^- b^\dagger),
\]
with \(g_{\rm eff}=\sqrt N\,g^2/\Delta_c\), \(\omega_b=N\,g^2/\Delta_c\), and \(\Delta_{\rm eff}=\Omega^2/\Delta_L+g^2/\Delta_c\). The resulting dynamics includes vacuum-Rabi oscillations, collapse and revival for an initial coherent state of the collective mode, stepwise Fock-state generation, Schrödinger-cat preparation, and Wigner-function reconstruction through selective displacements and atomic readout [1205.0979].

For traveling quantum pulses, a single cavity mode is insufficient. The cascaded-systems construction introduces “virtual” cavity modes \(\hat a_u\) and \(\hat a_v\) whose leak rates are chosen to emit or absorb prescribed temporal profiles \(u(t)\) and \(v(t)\). In the emission-only model,
\[
H_{\rm JCM-I}(t)=\frac{i}{2}\bigl[\sqrt\gamma\,g_u(t)\,\hat a_u^\dagger\hat\sigma^-
-\sqrt\gamma\,g_u^*(t)\,\hat a_u\hat\sigma^+\bigr],
\]
with jump operator
\[
L_{\rm JCM-I}(t)=g_u^*(t)\hat a_u+\sqrt\gamma\,\hat\sigma^-.
\]
The input–output version adds a second virtual cavity and a cross term \(\propto g_u(t)g_v^*(t)\hat a_u^\dagger\hat a_v\). In the interaction picture for \(v(t)=u(t)\), an auxiliary orthogonal mode survives explicitly, carrying non-Markovian memory. The source states directly that this differs from a “naïve” time-dependent JCM with coupling \(\chi(t)\propto u(t)\), because the correct couplings depend on temporal integrals of the pulse shape and the Lindblad operators collect radiation into all modes not explicitly retained [2403.04383].

Autoionizing polaritons furnish a continuum-coupled eJCM. The model contains a ground state \(\lvert g\rangle\), bright and dark autoionizing states \(\lvert a\rangle,\lvert b\rangle\), continua of opposite parity, and a quantized infrared mode \(\eta\). The interaction
\[
H_{\rm int}=g(\lvert a,n\rangle\langle b,n-1\rvert+\lvert b,n-1\rangle\langle a,n\rvert)
\]
couples the two autoionizing states, while Fano configuration interaction supplies complex widths. Diagonalization of the effective non-Hermitian \(2\times 2\) polariton Hamiltonian yields complex eigenenergies \(\hat E_{1,2}=E_0\pm \Delta/2-(i/2)\Gamma_{1,2}\). The widths \(\Gamma_\pm\) depend on mixing coefficients, photoionization rates, autoionization rates, and interference terms. At resonance, destructive interference can minimize \(\Gamma_-\), producing stabilization relative to the field-free width [2210.01712].

The lattice extension most directly called the extended Jaynes–Cummings–Hubbard model adds photon hopping and nearest-neighbor repulsion between atomic excitations:
\[
H=\sum_i\bigl[U(a_i^\dagger\sigma_i+\sigma_i^\dagger a_i)-\mu(n_i^a+n_i^\sigma)\bigr]
-t\sum_{\langle i,j\rangle}(a_i^\dagger a_j+\mathrm{h.c.})
+V\sum_{\langle i,j\rangle}n_i^\sigma n_j^\sigma.
\]
Large-scale worm quantum Monte Carlo on one-dimensional chains, square lattices, and triangular lattices shows that for \(V/U=0.4\), bipartite geometries exhibit a first-order solid–superfluid transition and no stable supersolid in the thermodynamic limit, whereas the triangular lattice supports a stable supersolid. At fixed \(\mu/U=-0.77\), the supersolid satisfies \(S(Q)>0\) and \(\rho_s>0\) simultaneously for \(0.018\lesssim t/U\lesssim 0.0245\) as \(L\to\infty\). The algorithm is reported to operate in continuous imaginary time, to measure stiffness from winding-number fluctuations, to reach \(\beta\sim 1500\), and to remain free of a sign problem in this bosonic setting [2010.01940].

## 6. Critical phenomena and computational complexity

The squeezed-photon extension provides an explicit interpolation between Jaynes–Cummings and Rabi limits. With
\[
H_{\rm eJCM}
=\hbar\omega a^\dagger a+\frac{\hbar\Omega}{2}\sigma_z
+\hbar g\bigl[B\sigma_++B^\dagger\sigma_-\bigr],
\qquad
B=\cosh r\,a+\sinh r\,a^\dagger,
\]
the Hamiltonian in the original basis contains both rotating and counter-rotating couplings weighted by \(\cosh r\) and \(\sinh r\). In the limit \(r\to 0\), the ordinary JCM is recovered exactly; in the limit \(r\to\infty\), one obtains the quantum Rabi Hamiltonian with effective coupling \(g_{\rm eff}=g e^r/2\). After a unitary transformation and Bogoliubov diagonalization, the source identifies a soft-mode condition \(A(A+4C)=0\), giving the critical coupling
\[
g_c^2(r)=\omega\bigl[\Omega+\omega e^{-2r}\bigr]e^{-2r}.
\]
At \(r\to\infty\) this reduces to \(g_{c,{\rm Rabi}}^2=\omega\Omega\), while at \(r\to 0\) it yields \(g_{c,{\rm JCM}}^2=\omega(\Omega+\omega)\). The associated phase structure is described as a second-order superradiant transition, with \(\langle a+a^\dagger\rangle\) serving as order parameter [2405.15876].

The driven Jaynes–Cummings–Rabi correspondence to the Dirac equation identifies a different critical phenomenon: quasienergy collapse at \(\epsilon_c=g/2\). Below critical drive, the spectrum is discrete and Landau-like; above critical drive, it is continuous. Taken together with the squeezed-photon model, these results show that the phrase “phase transition” in eJCM contexts can denote distinct mechanisms: a soft-mode superradiant transition in one construction and a discrete-to-continuous quasienergy collapse in another. This suggests that the term eJCM is structurally broad rather than tied to a single universal critical theory [1808.06727], [2405.15876].

The same breadth appears in quantum-simulation complexity analysis for multimode eJCMs under the rotating-wave approximation. For \(N_F\) truncated bosonic modes,
\[
H_{eJCM}=H_{\rm Photon}+H_{\rm Atom}+H_{\rm Photon\mbox{-}Atom},
\]
with
\[
H_{\rm Photon}=\sum_{m=1}^{N_F}\omega_m a_m^\dagger a_m,\quad
H_{\rm Atom}=\frac12\omega\sigma^z,\quad
H_{\rm Photon\mbox{-}Atom}=\sum_{m=1}^{N_F}\gamma_m(a_m\otimes \sigma^+ + a_m^\dagger\otimes \sigma^-).
\]
After Pauli decomposition, first- and second-order product formulas yield explicit operator-norm error bounds. The first-order scheme requires
\[
N_T\gtrsim \frac{T^2}{\varepsilon}[\cdots],
\]
while the second-order scheme satisfies \(\|e^{-iHT}-S_2(T)\|\le T^3 C/N_T^2\) and therefore \(N_T\gtrsim (C T^3/\varepsilon)^{1/2}\). The paper also treats pure-state binary encoding, mixed-state vectorization through the Choi–Jamiolkowski isomorphism, and fault-tolerant compilation on a rotated surface-code architecture. For the parameter set \(N_F=5\), \(k=4\) (\(n=15\)), \(\varepsilon=10^{-2}\), the reported estimate is approximately \(10^{12}\) \(T\) gates, about \(10^6\) physical qubits, and execution in \(\sim 10^8\) surface-code cycles [2509.13546].

In aggregate, the cited literature presents the eJCM less as a single model than as a technical umbrella for Jaynes–Cummings-type Hamiltonians enlarged by non-RWA physics, generalized couplings, additional modes, lattice structure, or effective reductions from more complex systems. The common payoff is exact or controlled access to regimes in which parity, squeezing, coherent displacement, collective bosonization, continuum interference, many-body order, or relativistic analogies materially change the spectrum and dynamics.

Source: https://www.emergentmind.com/topics/extended-jaynes-cummings-model-ejcm