quantum electrodynamics cavity QED Jaynes Cummings dressed states purcell effect
# Cavity Quantum Electrodynamics: Jaynes-Cummings Hamiltonian, Dressed Atom States, and Purcell Enhancement
## 1. Introduction and Historical Context
Cavity Quantum Electrodynamics (CQED) represents the theoretical and experimental domain where a single quantum emitter—typically an atom, ion, or artificial atom (e.g., superconducting qubit, quantum dot)—interacts coherently with the quantized electromagnetic field confined within an optical or microwave cavity. The foundational framework emerges from the interplay between spontaneous emission (a purely quantum phenomenon absent in classical electromagnetism) and the finite density of electromagnetic modes in a confined geometry. In a bulk medium, an excited atom radiates into a continuum of modes, releasing energy monotonically. In contrast, inside a high-quality-factor cavity with mode volume V and quality factor Q, the density of electromagnetic modes near the cavity resonance frequency ω_c is enhanced, fundamentally altering the emitter's decay channels and enabling coherent light-matter interactions.
The Jaynes-Cummings model, introduced by Jaynes and Cummings in 1963, provides the canonical theoretical description of a two-level atom interacting with a single cavity mode. This model has proven instrumental in understanding phenomena such as vacuum Rabi splitting (observable as an avoided crossing in the cavity transmission spectrum when the atom is tuned through resonance), Rabi oscillations of excited state population under resonant driving, and the Purcell effect (modification of spontaneous emission rates in cavities). Beyond pedagogical value, the Jaynes-Cummings model serves as the foundation for quantum optics experiments ranging from micromasers and optical cavities to circuit QED platforms, where superconducting qubits replace atoms and lumped LC circuits replace physical cavities.
## 2. Quantization of the Electromagnetic Field in a Cavity
Begin with the quantized electromagnetic field confined within a perfectly conducting cavity. The full Hamiltonian for a mode with frequency ω_c reads
$$H_{ ext{cavity}} = \hbar \omega_c \left( a^\dagger a + \frac{1}{2} ight),$$
where $a$ and $a^\dagger$ are the standard photon annihilation and creation operators obeying the commutation relation $[a, a^\dagger] = 1$. The number state $|n
angle$ represents a state with exactly n photons in the cavity, with energy eigenvalue $(n + 1/2)\hbar\omega_c$.
In atomic units or when referencing energies relative to the ground state, the zero-point energy $\hbar\omega_c/2$ is often omitted, yielding $H_{ ext{cavity}} = \hbar\omega_c a^\dagger a$.
The cavity mode is characterized by a spatial profile $\mathbf{E}_0(\mathbf{r})$ satisfying Maxwell's equations with appropriate boundary conditions. For simplicity, the interaction with a two-level atom situated at position $\mathbf{r}_0$ (typically near the antinode of the standing wave pattern to maximize coupling) depends on the electric field amplitude $E_0(\mathbf{r}_0)$ at the emitter location.
## 3. Two-Level Atom Hamiltonian and Dipole Interaction
A two-level atom (or generic two-level system) possesses a ground state $|g
angle$ with energy 0 and an excited state $|e
angle$ with energy $\hbar\omega_e$. The atomic Hamiltonian in the bare-state basis is
$$H_{ ext{atom}} = \frac{\hbar\omega_e}{2} \sigma_z,$$
where $\sigma_z = |e
angle\langle e| - |g
angle\langle g|$ is the Pauli z-matrix. The raising and lowering operators are $\sigma_+ = |e
angle\langle g|$ (promoting ground to excited) and $\sigma_- = |g
angle\langle e|$ (demoting excited to ground).
The interaction of the two-level system with the cavity's electric field is mediated by the electric dipole moment $\mathbf{d} = d(\sigma_+ + \sigma_-)$, where d is the matrix element of the dipole operator between the ground and excited states. Under the standard rotating-wave approximation (RWA), which neglects rapidly oscillating counter-rotating terms $a\sigma_+$ and $a^\dagger\sigma_-$ when the atom and cavity are near-resonant, the interaction Hamiltonian becomes
$$H_{ ext{int}} = \hbar g (a^\dagger \sigma_- + a \sigma_+),$$
where $g$ is the vacuum Rabi frequency (also called the single-photon Rabi frequency or coupling strength):
$$g = \frac{\omega_e d E_0(\mathbf{r}_0)}{\hbar V^{1/2}}.$$
Here $V$ is the quantization volume (the cavity volume in the idealized limit). The coupling strength g scales inversely with cavity volume: smaller cavities produce stronger coupling.
## 4. Jaynes-Cummings Hamiltonian
The complete Jaynes-Cummings Hamiltonian, omitting the zero-point energy term, reads
$$H_{ ext{JC}} = \hbar \omega_c a^\dagger a + \frac{\hbar \omega_e}{2} \sigma_z + \hbar g (a^\dagger \sigma_- + a \sigma_+).$$
This Hamiltonian conserves the excitation number $\mathcal{N} = a^\dagger a + \sigma_+ \sigma_-$, which counts the total number of excitations (photons + atomic excitations). This conservation law permits block-diagonalizing the Hamiltonian into sectors labeled by N = 0, 1, 2, ....
### 4.1 Hilbert Space Structure and Block Diagonalization
Within the N-excitation sector, the eigenstates are superpositions of the product states. For N = 1, the relevant basis states are $|1, g
angle$ (one photon, atom in ground state) and $|0, e
angle$ (zero photons, atom in excited state). The restriction of H_JC to this subspace is the 2×2 matrix
$$H_1 = \begin{pmatrix} \hbar\omega_c & \hbar g \\ \hbar g & \frac{\hbar\omega_e}{2} \end{pmatrix}.$$
The eigenvalues are obtained by solving the secular equation:
$$\det \begin{pmatrix} \hbar\omega_c - E & \hbar g \\ \hbar g & \frac{\hbar\omega_e}{2} - E \end{pmatrix} = 0.$$
Let $\Delta = \omega_e - \omega_c$ (the atom-cavity detuning). The eigenvalues of the N = 1 sector are
$$E_{+,1} = \frac{\hbar}{2}(\omega_e + \omega_c) + \frac{\hbar}{2}\sqrt{\Delta^2 + 4g^2},$$
$$E_{-,1} = \frac{\hbar}{2}(\omega_e + \omega_c) - \frac{\hbar}{2}\sqrt{\Delta^2 + 4g^2}.$$
The energy splitting between the two dressed states is
$$\Delta E_1 = E_{+,1} - E_{-,1} = \hbar\sqrt{\Delta^2 + 4g^2}.$$
On resonance ($\Delta = 0$), this splitting simplifies to $\Delta E_1 = 2\hbar g$, defining the vacuum Rabi splitting.
For the general N-excitation manifold, the eigenvalues are
$$E_{\pm, N} = \hbar\left[ \omega_c \left(N + \frac{1}{2} ight) + \frac{\omega_e}{2} ight] \pm \frac{\hbar}{2}\sqrt{\Delta^2 + 4g^2(N+1)}.$$
The key observation is that the Rabi splitting scales as $\sqrt{N+1}$, so higher-photon-number states experience larger splittings.
### 4.2 Dressed Atom States
The eigenstates in the N = 1 manifold are the dressed states (also called polariton states), which are coherent superpositions of the bare cavity and atomic basis states:
$$|\psi_{+,1}
angle = \sin heta |1,g
angle + \cos heta |0,e
angle,$$
$$|\psi_{-,1}
angle = \cos heta |1,g
angle - \sin heta |0,e
angle,$$
where the mixing angle $ heta$ is defined by
$$ an(2 heta) = \frac{2g}{\Delta}.$$
On resonance ($\Delta = 0$, $ heta = 45°$), the dressed states are equal superpositions:
$$|\psi_{+,1}
angle = \frac{1}{\sqrt{2}}(|1,g
angle + |0,e
angle),$$
$$|\psi_{-,1}
angle = \frac{1}{\sqrt{2}}(|1,g
angle - |0,e
angle).$$
These symmetric and antisymmetric combinations represent bright and dark polariton eigenmodes of the coupled system.
## 5. Vacuum Rabi Oscillations and Coherent Dynamics
Starting from the dressed state $|\psi_{+,1}
angle$ at t = 0, the system evolves under the Jaynes-Cummings Hamiltonian as
$$|\psi(t) angle = e^{-iH_{ ext{JC}}t/\hbar} |\psi_{+,1} angle.$$
Decomposing this initial state in the bare basis and allowing each component to acquire its phase factor:
$$|\psi(t) angle = e^{-i(E_{+,1} + E_{-,1})t/(2\hbar)} \left[ e^{-i(E_{+,1} - E_{-,1})t/(2\hbar)} |\psi_{+,1} angle + e^{-i(E_{-,1} - E_{+,1})t/(2\hbar)} |\psi_{-,1} angle ight].$$
The relative phase between the two dressed states oscillates at the Rabi frequency $\Omega_R = \sqrt{\Delta^2 + 4g^2}/(2)$, causing the expectation values of observables (e.g., atomic excitation $\sigma_+\sigma_-$) to oscillate periodically. On resonance, $\Omega_R = g$, and the population of the excited state oscillates as
$$P_e(t) = \sin^2(gt).$$
This is the Rabi oscillation, a hallmark coherent phenomenon where atomic population oscillates back and forth between ground and excited states (equivalently, between having zero and one photon in the cavity) without net energy dissipation.
## 6. Vacuum Rabi Splitting in Transmission Spectra
Experimentally, the vacuum Rabi splitting is observed in the cavity transmission spectrum. When a weak probe laser at frequency $\omega_p$ is incident on the cavity, the transmission exhibits a dip (or doublet structure when the atom and cavity are tuned near resonance). At resonance ($\omega_e = \omega_c$), the transmission spectrum displays two symmetric peaks separated by the vacuum Rabi splitting $2\hbar g$, corresponding to absorption into the upper and lower dressed states. This spectral splitting is a direct signature of strong light-matter coupling and is a diagnostic feature of the cavity QED regime.
## 7. Purcell Effect and Spontaneous Emission Modification
In free space, a two-level atom excited state decays via spontaneous emission at the Einstein A coefficient rate:
$$\Gamma_0 = \frac{\omega_e^3 d^2}{3\pi\epsilon_0\hbar c^3}.$$
Inside a cavity, this decay rate is modified by the local density of electromagnetic states (LDOS). Purcell's theory predicts that the decay rate becomes
$$\Gamma(\omega_e) = \Gamma_0 + \Gamma_{ ext{cavity}},$$
where the cavity contribution $\Gamma_{ ext{cavity}}$ depends on the cavity's quality factor Q, mode volume V, and the atom's position within the cavity.
The Purcell factor (or enhancement factor) quantifies the ratio of the cavity-enhanced decay to the vacuum decay:
$$F_P = 1 + \frac{3}{4\pi^2} \left(\frac{\lambda}{n} ight)^3 \frac{Q}{V},$$
where $\lambda = 2\pi c/\omega_e$ is the wavelength of the atomic transition in vacuum, n is the refractive index of the cavity material, and Q/V is the quality factor per unit volume.
For cavities with high Q and small V (e.g., photonic crystal nanocavities with Q ~ 10⁴ and V ~ (λ/n)³), the Purcell factor can exceed 100, dramatically enhancing spontaneous emission rates. This enhanced decay can be viewed as a consequence of the cavity mode's density of states at the atomic transition frequency being much larger than the free-space density.
When the Purcell factor is large (F_P >> 1), spontaneous emission is predominantly into the cavity mode (not into radiation). This regime is advantageous for quantum information processing because it enables rapid, high-fidelity single-photon generation and efficient atom-photon entanglement.
## 8. Strong and Weak Coupling Regimes
The light-matter interaction regime is characterized by the ratio of the single-photon Rabi frequency g to the dissipation rates. Define:
- κ: the cavity mode decay rate (inverse cavity photon lifetime)
- γ: the atomic spontaneous emission rate
The strong coupling regime occurs when $g \gg \kappa, \gamma$, so the coherent interaction g dominates over dissipation. In this regime, vacuum Rabi splitting is clearly resolved, and coherent oscillations (Rabi oscillations) persist for many cycles before being washed out by decay.
The weak coupling regime occurs when $g \ll \kappa, \gamma$, meaning dissipation dominates. Here, the atom decays before coherent splitting can develop, and vacuum Rabi splitting is not observed.
The crossover between these regimes is often quantified by the cooperativity parameter:
$$C = \frac{g^2}{\kappa \gamma}.$$
Strong coupling requires C > 1; weak coupling corresponds to C << 1.
## 9. Lindblad Master Equation for Open-System Dynamics
The Jaynes-Cummings model assumes an isolated system with no dissipation. Real cavities experience photon loss (characterized by the cavity decay rate κ), and atoms undergo spontaneous emission (rate γ) into modes outside the cavity. These dissipative processes are incorporated via the Lindblad master equation:
$$\frac{d ho}{dt} = -\frac{i}{\hbar}[H_{ ext{JC}}, ho] + \mathcal{L}_\kappa[ ho] + \mathcal{L}_\gamma[ ho],$$
where the Lindblad jump operators are
$$\mathcal{L}_\kappa[
ho] = \kappa \left( a
ho a^\dagger - \frac{1}{2}\{a^\dagger a,
ho\}
ight),$$
$$\mathcal{L}_\gamma[
ho] = \gamma \left( \sigma_-
ho \sigma_+ - \frac{1}{2}\{\sigma_+ \sigma_-,
ho\}
ight).$$
The κ term models cavity photon loss (via the annihilation operator a), while the γ term models atomic decay (via σ₋). The anti-commutator terms ensure that the density matrix trace is conserved.
For weak dissipation (g >> κ, γ), the dressed states maintain partial coherence, and Rabi oscillations are damped but still observable. The dissipation-limited visibility of Rabi oscillations is approximately
$$V = \frac{4g^2}{(2g)^2 + (\kappa + \gamma)^2}.$$
## 10. Experimental Realizations and Applications
### 10.1 Optical Cavities
Traditional CQED uses high-finesse optical cavities (Q ~ 10⁴ to 10⁶) containing neutral atoms or trapped ions. These systems have demonstrated vacuum Rabi splitting, Rabi oscillations, and strong coupling. A landmark example is the observation of Rabi oscillations in a Fabry-Perot cavity with ultracold rubidium atoms, where g reached values comparable to κ and γ.
### 10.2 Photonic Crystal Cavities
Photonic crystal nanocavities achieve extremely high Q and very small V by confining light in a structure with a photonic bandgap. These cavities can be integrated with quantum dots or other solid-state emitters, enabling strong coupling with compact, chip-scale devices. The Purcell factor in such cavities routinely exceeds 100.
### 10.3 Circuit QED
Superconducting qubits acting as artificial atoms interact with microwave resonators formed by lumped LC circuits or coplanar waveguides. Circuit QED has become the leading platform for quantum computing, where the controllable qubit-resonator coupling (g ~ 1-100 MHz) and engineered dissipation rates (κ, γ ~ 1-100 kHz) allow exploration of the full range from weak to strong coupling and even ultrastrong coupling (g ~ ω_c).
### 10.4 Quantum Information Processing
Strong coupling enables high-fidelity quantum gates and single-photon generation. In particular, cavity-mediated interactions between qubits provide an alternative to direct gate implementations and have been instrumental in demonstrating quantum error correction and quantum simulation in superconducting processors.
## 11. Cavity QED Beyond the Jaynes-Cummings Model
The simplest Jaynes-Cummings model assumes a single cavity mode, a two-level emitter, and weak driving. Real systems often require extensions:
1. Multi-level atoms: More realistic atoms possess hyperfine structure or multiple excited states, leading to additional decay channels and mode coupling.
2. Multiple cavity modes: Cavities support multiple resonant frequencies. Emitters can couple to several modes simultaneously, complicating the spectral response.
3. Strong driving: When the driving field amplitude becomes comparable to the Rabi frequency, nonlinear effects emerge, including multi-photon transitions and frequency shifts.
4. Ultrastrong coupling: At extremely high coupling strengths (g ~ ω_c), the rotating-wave approximation fails, and counter-rotating terms become important, leading to qualitatively different phenomena.
5. Cavity-mediated interactions: Multiple emitters in the same cavity interact via the cavity field, enabling collective effects like superradiance and subradiance.
## 12. Numerical Solver: Jaynes-Cummings Energy Spectrum and Rabi Oscillations
The following Python code constructs the Jaynes-Cummings Hamiltonian for multiple excitation manifolds, diagonalizes it to extract dressed state energies, and simulates the time-dependent dynamics under the full Hamiltonian (without dissipation) as well as under the Lindblad equation (with dissipation):
import numpy as np
import matplotlib.pyplot as plt
from scipy.integrate import odeint
from scipy.linalg import eigh
def jaynes_cummings_hamiltonian(n_max, omega_c, omega_e, g, delta_detuning=0):
"""
Construct the Jaynes-Cummings Hamiltonian matrix in the Fock basis.
Basis ordering: |n, ground>, |n, excited> for n = 0, 1, ..., n_max
Total Hilbert space dimension: 2*(n_max + 1)
"""
dim = 2 * (n_max + 1)
H = np.zeros((dim, dim), dtype=complex)
omega_e_actual = omega_c + delta_detuning
for n in range(n_max + 1):
# |n, g> state index
idx_g = 2 * n
# |n, e> state index
idx_e = 2 * n + 1
# Cavity energy
H[idx_g, idx_g] = omega_c * n
H[idx_e, idx_e] = omega_c * n + omega_e_actual / 2
# Interaction terms
if n < n_max:
# a^\dagger \sigma_- term: |n, e> -> |n+1, g>
H[idx_e, 2*(n+1)] += g * np.sqrt(n + 1)
# a \sigma_+ term: |n+1, g> -> |n, e>
H[2*(n+1), idx_e] += g * np.sqrt(n + 1)
return H
def solve_jaynes_cummings_spectrum(n_max, omega_c, omega_e, g, detuning=0):
"""
Diagonalize the Jaynes-Cummings Hamiltonian and return eigenvalues and eigenvectors.
"""
H = jaynes_cummings_hamiltonian(n_max, omega_c, omega_e, g, detuning)
eigenvalues, eigenvectors = eigh(H)
return eigenvalues, eigenvectors, H
def plot_energy_ladder(n_max, omega_c, omega_e, g, detunings):
"""
Plot the dressed state energy ladder as a function of detuning.
"""
energies_vs_detuning = []
for delta in detunings:
eigenvalues, _, _ = solve_jaynes_cummings_spectrum(n_max, omega_c, omega_e, g, delta)
energies_vs_detuning.append(eigenvalues)
energies_vs_detuning = np.array(energies_vs_detuning)
fig, ax = plt.subplots(figsize=(10, 6))
for i in range(energies_vs_detuning.shape[1]):
ax.plot(detunings / g, energies_vs_detuning[:, i] / (np.hbar * g), 'b-', linewidth=0.8)
ax.set_xlabel('Detuning Δ / g', fontsize=12)
ax.set_ylabel('Energy (in units of ℏg)', fontsize=12)
ax.set_title('Jaynes-Cummings Energy Spectrum: Dressed State Ladder', fontsize=13)
ax.grid(True, alpha=0.3)
plt.tight_layout()
return fig
def rabi_oscillation_coherent(n_max, omega_c, omega_e, g, t_values, initial_state_index=1):
"""
Compute time-dependent population oscillations (Rabi oscillations) without dissipation.
Initial state: eigenstate of H at delta=0.
"""
eigenvalues, eigenvectors, _ = solve_jaynes_cummings_spectrum(n_max, omega_c, omega_e, g, detuning=0)
# Use the first excited dressed state as initial condition
psi_0 = eigenvectors[:, initial_state_index]
# Time evolution operator: U(t) = exp(-i H t / hbar)
populations_over_time = []
for t in t_values:
phases = np.exp(-1j * eigenvalues * t)
psi_t = np.sum([phases[i] * eigenvectors[:, i] * np.conj(psi_0[i]) for i in range(len(eigenvalues))], axis=0)
# Compute excited state population
excited_pop = 0
for n in range(n_max + 1):
idx_e = 2 * n + 1
excited_pop += np.abs(psi_t[idx_e])**2
populations_over_time.append(excited_pop)
return np.array(populations_over_time)
def lindblad_master_equation(rho, t, H, kappa, gamma, n_max):
"""
Right-hand side of the Lindblad master equation for density matrix evolution.
rho: density matrix (flattened to 1D array)
H: Hamiltonian matrix
kappa: cavity decay rate
gamma: atomic spontaneous emission rate
"""
dim = 2 * (n_max + 1)
rho_mat = rho.reshape((dim, dim))
# Coherent evolution
drho_dt = -1j * (np.dot(H, rho_mat) - np.dot(rho_mat, H))
# Cavity loss: L_kappa = kappa * a rho a^\dagger - kappa/2 * {a^\dagger a, rho}
a_op = np.zeros((dim, dim), dtype=complex)
for n in range(n_max):
for m in [0, 1]: # ground and excited states
idx_from = 2*n + m
idx_to = 2*(n+1) + m
if idx_to < dim:
a_op[idx_to, idx_from] = np.sqrt(n + 1)
cavity_term = kappa * (np.dot(a_op, np.dot(rho_mat, a_op.conj().T))
- 0.5 * (np.dot(a_op.conj().T, np.dot(a_op, rho_mat))
+ np.dot(rho_mat, np.dot(a_op.conj().T, a_op))))
drho_dt += cavity_term
# Atomic decay: L_gamma = gamma * sigma_- rho sigma_+ - gamma/2 * {sigma_+ sigma_-, rho}
sigma_minus = np.zeros((dim, dim), dtype=complex)
for n in range(n_max + 1):
sigma_minus[2*n, 2*n + 1] = 1 # |n,g> <n,e|
atomic_term = gamma * (np.dot(sigma_minus, np.dot(rho_mat, sigma_minus.conj().T))
- 0.5 * (np.dot(sigma_minus.conj().T, np.dot(sigma_minus, rho_mat))
+ np.dot(rho_mat, np.dot(sigma_minus.conj().T, sigma_minus))))
drho_dt += atomic_term
return drho_dt.flatten()
def rabi_oscillation_with_dissipation(n_max, omega_c, omega_e, g, kappa, gamma, t_values):
"""
Compute Rabi oscillations including cavity decay and spontaneous emission.
"""
eigenvalues, eigenvectors, H = solve_jaynes_cummings_spectrum(n_max, omega_c, omega_e, g, detuning=0)
# Initialize in the upper dressed state |+, 1>
psi_0 = eigenvectors[:, 1]
rho_0 = np.outer(psi_0, psi_0.conj()).flatten()
# Solve Lindblad equation
rho_evolution = odeint(lindblad_master_equation, rho_0, t_values, args=(H, kappa, gamma, n_max))
# Extract excited state population
excited_populations = []
for rho_flat in rho_evolution:
rho_mat = rho_flat.reshape((2*(n_max+1), 2*(n_max+1)))
excited_pop = 0
for n in range(n_max + 1):
idx_e = 2 * n + 1
excited_pop += np.real(rho_mat[idx_e, idx_e])
excited_populations.append(excited_pop)
return np.array(excited_populations)
# Parameters
hbar = 1 # Natural units
omega_c = 1.0 # Cavity frequency (normalized)
omega_e = 1.0 # Atomic transition frequency (on resonance)
g = 0.1 # Coupling strength
kappa = 0.01 # Cavity decay rate
gamma = 0.01 # Atomic spontaneous emission rate
n_max = 5 # Maximum photon number
# Detuning scan for energy ladder
detunings = np.linspace(-0.5, 0.5, 100) * g
fig_spectrum = plot_energy_ladder(n_max, omega_c, omega_e, g, detunings)
# Time evolution for Rabi oscillations
t_max = 50 / g # Time in units of inverse Rabi frequency
t_values = np.linspace(0, t_max, 500)
# Coherent dynamics (no dissipation)
pop_coherent = rabi_oscillation_coherent(n_max, omega_c, omega_e, g, t_values, initial_state_index=1)
# Dynamics with dissipation
pop_dissipative = rabi_oscillation_with_dissipation(n_max, omega_c, omega_e, g, kappa, gamma, t_values)
# Plot results
fig, axes = plt.subplots(2, 2, figsize=(14, 10))
# Panel 1: Energy spectrum
ax = axes[0, 0]
for i in range(min(6, len(eigenvalues))):
eigenvalues, _, _ = solve_jaynes_cummings_spectrum(n_max, omega_c, omega_e, g, detuning=0)
ax.plot(detunings / g, (eigenvalues[1] - eigenvalues[0]) * np.ones_like(detunings) / (np.hbar * g), 'ro-', markersize=4, label='Vacuum Rabi splitting')
ax.set_xlabel('Detuning Δ / g', fontsize=11)
ax.set_ylabel('Splitting (ℏg)', fontsize=11)
ax.set_title('Vacuum Rabi Splitting', fontsize=12)
ax.grid(True, alpha=0.3)
# Panel 2: Coherent Rabi oscillations
ax = axes[0, 1]
ax.plot(t_values * g, pop_coherent, 'b-', linewidth=2, label='Excited state')
ax.set_xlabel('Rabi Time (gt)', fontsize=11)
ax.set_ylabel('P_e(t)', fontsize=11)
ax.set_title('Coherent Rabi Oscillations (no dissipation)', fontsize=12)
ax.grid(True, alpha=0.3)
ax.legend()
# Panel 3: Dissipative Rabi oscillations
ax = axes[1, 0]
ax.plot(t_values * g, pop_dissipative, 'r-', linewidth=2, label='With dissipation')
ax.plot(t_values * g, pop_coherent, 'b--', linewidth=1.5, alpha=0.7, label='Without dissipation')
ax.set_xlabel('Rabi Time (gt)', fontsize=11)
ax.set_ylabel('P_e(t)', fontsize=11)
ax.set_title('Damped Rabi Oscillations', fontsize=12)
ax.grid(True, alpha=0.3)
ax.legend()
# Panel 4: Cooperativity and decay rates
ax = axes[1, 1]
cooperativity = g**2 / (kappa * gamma)
ax.text(0.5, 0.7, f'Cooperativity C = g² / (κγ) = {cooperativity:.2f}',
ha='center', fontsize=12, transform=ax.transAxes, bbox=dict(boxstyle='round', facecolor='wheat'))
ax.text(0.5, 0.5, f'Coupling g / κ = {g/kappa:.2f}',
ha='center', fontsize=11, transform=ax.transAxes)
ax.text(0.5, 0.35, f'Coupling g / γ = {g/gamma:.2f}',
ha='center', fontsize=11, transform=ax.transAxes)
ax.text(0.5, 0.2, f'Cavity quality factor (estimated) Q ~ {omega_c / kappa:.0f}',
ha='center', fontsize=11, transform=ax.transAxes)
ax.axis('off')
ax.set_title('Cavity QED Regime Parameters', fontsize=12)
plt.tight_layout()
plt.show()
print("Jaynes-Cummings Hamiltonian Analysis Complete")
print(f"Vacuum Rabi splitting (resonance): 2g = {2*g:.4f}")
print(f"Cooperativity C = {cooperativity:.4f}")
print(f"Regime: {'Strong coupling' if cooperativity > 1 else 'Weak coupling'}")## 13. Coupling Strength Engineering and Purcell Factor Optimization
To achieve strong coupling (C >> 1), one must maximize g while minimizing κ and γ. The coupling strength g depends on the dipole matrix element d (which is material-dependent), the field strength E₀ (proportional to the inverse square root of mode volume), and the quantization volume. For a given emitter, the strategy is to engineer cavities with:
1. High quality factor Q: Reduces κ = ω_c / Q, extending the photon lifetime.
2. Small mode volume V: Increases the field strength E₀, thereby increasing g.
3. Optimal emitter placement: Position the emitter at the field antinode to maximize coupling.
The Purcell factor and cooperativity are related: F_P ≈ 2C for cavities with Q >> 1. Thus, engineering for high cooperativity naturally leads to enhanced spontaneous emission rates into the cavity.
## 14. Extensions: Cavity QED with Multiple Emitters
When multiple atoms occupy the same cavity, their interaction is mediated by the cavity field. This leads to superradiance (collective enhancement of emission when atoms are in-phase) and subradiance (suppression when out-of-phase). For N identical atoms in a cavity:
- Superradiant state: All atoms oscillate in-phase; collective decay rate scales as N·Γ.
- Subradiant states: Atoms oscillate with relative phase differences; decay rates scale as Γ or lower.
This collective behavior forms the basis for quantum simulation and collective quantum gates in cavity QED platforms.
## 15. Recent Advances and Future Perspectives
Modern cavity QED experiments have pushed beyond the simple two-level model to explore:
1. Nonlinear cavity QED: Using optical Kerr nonlinearities to achieve photon blockade and generate non-classical light.
2. Hybrid quantum systems: Coupling different quantum systems (e.g., superconducting qubits and mechanical oscillators) via cavities.
3. Topological cavity QED: Exploiting topological protection in photonic structures to enhance cavity properties.
4. Distributed quantum networks: Using cavities as nodes in quantum networks for quantum repeaters and distributed quantum computing.
The field continues to expand toward practical quantum technologies, including quantum computers based on cavity QED (e.g., Google's quantum supremacy experiments with superconducting qubits), quantum metrology at the Heisenberg limit, and quantum simulation of many-body physics.
## Conclusion
The Jaynes-Cummings model provides the foundational theoretical framework for understanding light-matter interaction in cavities. Vacuum Rabi splitting, Rabi oscillations, and Purcell enhancement emerge naturally from the model's structure, explaining how confined electromagnetic fields qualitatively alter spontaneous emission and enable coherent quantum control. The model's simplicity belies its power: it has guided the design of cavity QED experiments across optical, microwave, and nanophotonic platforms, and remains central to quantum information science. The companion Python solver demonstrates how to diagonalize the Hamiltonian, extract dressed state properties, and simulate realistic dissipative dynamics, providing a computational foundation for designing and interpreting cavity QED experiments.