Quantum Rabi Hamiltonian

• Quantum Rabi Hamiltonian is a fundamental model that describes the interaction between a two-level system and a quantized bosonic field without using the rotating wave approximation.
• It captures complex spectral phenomena such as parity symmetry, Juddian baselines, and deep-strong coupling effects through methods like the Bargmann–Fock representation and canonical transformations.
• The model underpins quantum simulation and circuit QED research, offering practical insights for numerical computation and experimental setups in light–matter interaction studies.

The Quantum Rabi Hamiltonian is a fundamental model of light–matter interaction widely used in quantum optics, cavity and circuit quantum electrodynamics (QED), and quantum information science. It describes the coupling between a two-level system (qubit) and a single mode of a quantized bosonic field without applying the rotating wave approximation, thus capturing the full dynamics including both energy-conserving and counter-rotating terms. Its mathematical structure encapsulates complex spectral, algebraic, and dynamical phenomena arising in both perturbative and deep-strong coupling regimes.

1. Formal Definition and Basic Structure

The canonical Quantum Rabi Hamiltonian (QRM), in units where $\hbar=1$, is given by

$$H = \omega\,a^\dagger a + \frac{\Omega}{2}\,\sigma_z + g\,\sigma_x (a + a^\dagger)$$

where:

• $a$, $a^\dagger$ are bosonic annihilation and creation operators for the field mode, $[a, a^\dagger]=1$;
• $\sigma_x$, $\sigma_z$ are Pauli matrices acting on the two-level system (qubit) of splitting $\Omega$;
• $\omega$ is the field mode frequency;
• $g$ is the light–matter (dipole) coupling strength.

The QRM retains both the rotating (Jaynes–Cummings) and counter-rotating (Bloch–Siegert) terms, thus remaining valid for all coupling strengths, up to the deep–strong coupling regime where $g / \omega \gtrsim 1$ (Feranchuk et al., 2016).

2. Spectral Theory and Exact Solutions

Analytical Structure and Integrability

The spectrum of the Quantum Rabi model is purely discrete and exhibits parity symmetry: the Hamiltonian commutes with the operator $\Pi = (-1)^{a^\dagger a}\sigma_z$, splitting the Hilbert space into even and odd parity chains (Reyes-Bustos et al., 2019, Vandaele et al., 2016). The Hamiltonian admits an explicit exact solution via several complementary formulations:

• Bargmann–Fock Representation: Mapping $a^\dagger \rightarrow z$, $a \rightarrow d/dz$ in the space of analytic functions yields two coupled first-order ODEs for the spinor-like wavefunctions $f_1(z), f_2(z)$ (Vandaele et al., 2016).
• Canonical (Birkhoff) Transformation: Via a Birkhoff gauge transformation, the problem reduces to a canonical system whose quantization condition is the single-valuedness of confluent hypergeometric solutions, leading to integer quantum numbers associated with Juddian baselines (Vandaele et al., 2016).
• Progressive Diagonalization Scheme (PDS): Decomposes the infinite-dimensional problem into a sequence of finite-degree polynomial equations (the so-called “leapfrog” or “rational” secular equations), efficiently returning low-lying eigenvalues and eigenstates (Pan et al., 2010).
• Transcendental Equations and Parity Sectors: The roots of transcendental functions (built from series solutions or continued fractions) in each parity sector provide the complete spectrum (Vandaele et al., 2016, Pan et al., 2010).

A salient aspect is the existence of Juddian baselines, $E_k^{(0)} = k - g^2$ (for $k \in \mathbb{N}$), where level crossings occur between parity chains, corresponding to exact solutions at isolated parameter values (Vandaele et al., 2016).

Strong-Coupling Spectral Limit and Zeta Functions

In the ultra-strong coupling limit $g \rightarrow \infty$, the spectral zeta function of the QRM converges to the Hurwitz zeta function, and the eigenvalue sequence $E_n(g) + g^2 \rightarrow n$ (i.e., the spectrum approaches that of the bare oscillator, modulo a global shift) (Hiroshima et al., 2024). The construction of the heat kernel and explicit spectral zeta allows the analysis of spectral statistics and ground-state path measures (Reyes-Bustos et al., 2019, Hiroshima et al., 2024).

3. Algebraic Symmetry and Duality

The Quantum Rabi model possesses a rich algebraic structure:

• Parity–Duality SU(2) Algebra: The set $\{ \Pi, D_1, D_2 \}$, with $D_j = \sigma_j P$ and $P = e^{i\pi a^\dagger a}$, forms a closed SU(2) algebra. This allows the QRM to be mapped, via duality transformations, to effective bosonic, fermionic, or “coupling-only” Hamiltonians (Omolo, 2021).
• Parity Sectors: Each eigenstate is labeled by its parity eigenvalue, and the spectrum is composed of paired levels of opposite parity, except at Judd points where crossings occur (Vandaele et al., 2016, Omolo, 2021).
• Symmetry Breaking and Generalizations: Extensions to multi-qubit systems, $N$-level generalizations, and inclusion of both one- and two-photon terms may break discrete symmetries or lead to new chain decompositions and approximate solutions (Albert, 2011, Xie et al., 2018, Peng et al., 2015).

4. Physical Regimes and Effective Theories

Multi-mode Embedding and Validity of Single-Mode Approximation

The QRM emerges from a more general multi-mode TLS–boson system via a canonical transformation that separates a collective (“bright”) mode from “dark” fluctuations. The single-mode model is physically accurate when the collective coupling $g$ dominates over residual fluctuations: $$\frac{\max_{k \in \Delta} |\omega_k - \omega_c|}{\sqrt{N}} \ll \bar{g}$$ where the collective frequency $\omega_c$ and coupling $\bar{g}$ are appropriately averaged (Feranchuk et al., 2016).

Rotating Wave Approximation (RWA) and Jaynes–Cummings Limit

For $g/\omega \ll 1$ and near resonance, the rotating wave approximation (RWA) reduces the model to the Jaynes–Cummings Hamiltonian by dropping counter-rotating terms, yielding analytic eigenstates and the well-known vacuum Rabi splitting. The full QRM describes correction effects such as the Bloch–Siegert shift and spectral nonlinearity as $g$ increases (Feranchuk et al., 2016, Gartner et al., 2021).

Semiclassical Limit

A rigorous semiclassical limit is obtained by taking $g \to 0$, $|\alpha| \to \infty$ with $A = g|\alpha|$ fixed (displaced Fock-state basis), so the QRM reduces unambiguously to the semiclassical driven two-level model,

$$H_{sc}(t) = \frac{1}{2}\Omega\,\sigma_z + 2A\,\sigma_x \cos(\omega_0 t)$$

This procedure is essential for resolving the quantum-to-classical crossover and controlling quantum back-action (Irish et al., 2022).

5. Extensions: Diamagnetic Term and Generalized Models

QRM with $A^2$ (Diamagnetic) Term

The complete QED Hamiltonian for light–matter coupling includes the diamagnetic $A^2$ term: $$H_{\mathrm{QRM+A^2}} = \omega\,a^\dagger a + \frac{\omega_0}{2}\,\sigma_z + g\sigma_x(a + a^\dagger) + D(a + a^\dagger)^2$$ with $D \simeq g^2/\omega$. By an exact squeezing (Bogoliubov) transformation, this maps onto a QRM with renormalized frequency $\omega' = \omega\sqrt{1+4k}$ and reduced coupling $g' = g/(1+4k)^{1/4}$ (2002.03702, Boutakka et al., 2024). Physically, the $A^2$ term shifts the ground-state energy upward, suppresses spectral crossings, and in the ultra-strong regime, limits ground-state entanglement as quantified by von Neumann entropy (Boutakka et al., 2024).

Generalized Quantum Rabi Models

The addition of two-photon ($g_2$) terms and multi-level atomic structure introduces further complexity: $$H = \frac{\Delta}{2}\sigma_z + \omega a^\dagger a + g_1 \sigma_x(a + a^\dagger) + g_2 \sigma_x(a^2 + a^{\dagger 2})$$ This breaks parity symmetry for $g_1 g_2 \neq 0$, resulting in cubic secular equations for eigenvalues, multimodal vacuum Rabi splittings, and new phenomena such as two dominant Rabi frequencies and altered sideband structure in the emission spectrum (Xie et al., 2018).

Multi-qubit and $N$-State Generalizations

• Two-Qubit QRM: Hidden algebraic structure yields quasi-exact finite-photon solutions and exact “dark states” at specific detunings; the full spectrum is given via convergent power series constructed with Bogoliubov operators (Peng et al., 2015).
• $N$-State Rabi Models: Group-theoretical analysis reveals $Z_N$ symmetry and chain decomposition, with novel parity inversion phenomena at strong coupling for $N\geq4$ (Albert, 2011).

6. Heat Kernel and Spectral Functions

The QRM's heat kernel admits an explicit analytic expansion as a 2×2 matrix-valued function, obtained via the Trotter–Kato product formula and orbit combinatorics of the symmetric group (Reyes-Bustos et al., 2019). This permits nonperturbative calculation of the partition function, spectral zeta, and dynamical propagator, providing access to thermodynamic and quantum statistical properties even at arbitrary coupling. Explicit representations for observables under the ground-state path measure are available in terms of expectations with respect to stochastic processes (Ornstein–Uhlenbeck and Poisson spin-jump process) (Hiroshima et al., 2024).

7. Computational and Experimental Relevance

The Quantum Rabi Hamiltonian serves as a benchmark for quantum simulation, circuit QED, and quantum devices operating beyond the RWA. Its exact solution and effective Hamiltonians enable precise computation of operator averages (photon number, excitations) and guide the experimental exploration of the ultra- and deep-strong coupling regimens, including the prescription for including the $A^2$ term to avoid spurious superradiant phase transitions in closed cavities (2002.03702, Gartner et al., 2021). Efficient numerical schemes based on polynomial secular equations or Bogoliubov operator expansions allow for high-precision determination of low-lying spectra with moderate computational resources (Pan et al., 2010).

The rich interplay of symmetry, integrability, and quantum-classical crossover continues to stimulate theoretical physics and practical quantum technology research.