The Interaction Hamiltonian between light and matter§
We begin with the standard electric dipole interaction Hamiltonian:
H^I=−d^⋅E^(r)
Expanding the dipole operator d^=der(σ^++σ^−) and the quantized multi-mode electric field operator yields the full interaction expression:
H^I=−k,λ∑Ek(σ^++σ^−)der⋅(b^k,λek,λeik⋅r+b^k,λ†ek,λ∗e−ik⋅r)=−k,λ∑Ek(σ^+b^k,λder⋅ek,λeik⋅r+σ^+b^k,λ†der⋅ek,λ∗e−ik⋅r+σ^−b^k,λder⋅ek,λeik⋅r+σ^−b^k,λ†der⋅ek,λ∗e−ik⋅r)
Moving into a spinning basis (The Interaction Picture)§
Take our unperturbed baseline Hamiltonian as:
H^0=ℏωegσ^e+k,λ∑ℏωkb^k,λ†b^k,λ
To isolate the dynamics driven purely by the interaction, we transform our operators into the rotating frame via:
A^(t)=eiℏH^0tA^e−iℏH^0t
Atomic Excitation Operators§
By the Hadamard Lemma, the time-evolution of the operators depends entirely on their commutation relations with H^0. For the atomic system, we have:
[σ^e,σ^+]=σ^+
Using the identity eab^B^e−ab^=eaλB^ which holds when [b^,B^]=λB^, we evaluate the time-dependent raising operator:
σ^+(t)=σ^+eiωegt
Taking the Hermitian conjugate yields the lowering operator:
σ^−(t)=σ^−e−iωegt
Photon Field Operators§
Using the identical algebraic structure for the bosonic field modes where [b^†b^,b^]=−b^, we obtain:
b^k,λ(t)=b^k,λe−iωkt
b^k,λ†(t)=b^k,λ†eiωkt
Combined Operators in the Interaction Picture§
Substituting these time-dependent forms back into H^I reveals that the terms separate into two distinctly behaving categories based on their phase frequencies.
Invariants§
b^†b^=n^b^b^†=n^+1σ^+σ^−=σ^eσ^−σ^+=σ^g
Co-rotating Terms (Difference Frequency)§
These terms represent energy-conserving processes where an atomic transition is balanced by the creation or destruction of a photon:
b^k,λ†σ^−ei(ωk−ωeg)t
b^k,λσ^+e−i(ωk−ωeg)t
Counter-rotating Terms (Sum Frequency)§
These terms represent virtual fluctuations where a photon is created alongside an atomic excitation, or destroyed alongside a de-excitation:
b^k,λ†σ^+ei(ωk+ωeg)t
b^k,λσ^−e−i(ωk+ωeg)t
Rotating Wave Approximation§
When the coupling parameter g (which collects the dipole moments and field amplitudes) is small relative to the optical frequencies, the counter-rotating terms oscillate at extreme speeds (ωk+ωeg) and average out to zero over macroscopic timescales.
Assuming a near-resonant scenario where ωk≈ωeg, the phase clocks of the co-rotating terms slow down completely, dominating the physical evolution. Dropping the counter-rotating terms leaves us with the final RWA Interaction Hamiltonian:
H^int(t)=k,λ∑(ℏgk,λb^k,λσ^+e−i(ωk−ωeg)t+ℏgk,λ∗b^k,λ†σ^−ei(ωk−ωeg)t)
This will be established with rigour in the next article.