Quantum Rabi

Author: TheoryLand

The quantum Rabi model is the minimal fully quantized light-matter Hamiltonian, while the Jaynes-Cummings model is its rotating-wave approximation [1, 2]. This demo propagates the initial state e,0|e,0\rangle in a truncated bosonic Fock basis and compares the full excited-state trace against the Jaynes-Cummings reference.

The sliders vary the boson cutoff nmaxn_{\max} and coupling gg. At weak coupling and modest times the two curves nearly coincide; as gg grows, the counter-rotating terms drive the expected departure from the simple cos2(gt)\cos^2(gt) reference.

Method

We use the resonant Hamiltonian

H=ωaa+ω02σz+g(a+a)(σ++σ),H=\omega a^\dagger a+\frac{\omega_0}{2}\sigma_z+g(a+a^\dagger)(\sigma_++\sigma_-),

with ω=ω0=1\omega=\omega_0=1 and bosonic truncation

n=0,,nmax.n=0,\dots,n_{\max}.

After removing the trivial free phases, the amplitudes satisfy nearest-neighbor ladder equations with explicit e±i(ω+ω0)te^{\pm i(\omega+\omega_0)t} counter-rotating factors. The demo integrates that truncated system with RK4, so increasing nmaxn_{\max} is the direct numerical check on how much bosonic support is needed at the selected coupling.

The comparison curve keeps only the rotating-wave terms,

HJC=ωaa+ω02σz+g(aσ++aσ),H_{\mathrm{JC}}=\omega a^\dagger a+\frac{\omega_0}{2}\sigma_z+g(a\sigma_++a^\dagger\sigma_-),

which for e,0|e,0\rangle gives

PeJC(t)=cos2(gt).P_e^{\mathrm{JC}}(t)=\cos^2(gt).

So the plotted difference is a direct visual measure of how much the neglected counter-rotating sector matters at the chosen (g,nmax)(g,n_{\max}).

References

  1. [1] E. T. Jaynes and F. W. Cummings. “Comparison of quantum and semiclassical radiation theories with application to the beam maser,” Proceedings of the IEEE 51 (1), 89--109 (1963). DOI: 10.1109/PROC.1963.1664.
  2. [2] D. Braak. “Integrability of the Rabi model,” Physical Review Letters 107 (10), 100401 (2011). DOI: 10.1103/PhysRevLett.107.100401.