Skip to content

Harmonic Perturbations

A harmonic perturbation is a weak drive with a single angular frequency. It is the basic model for absorption, stimulated emission, driven transitions, spectroscopy, and the weak-drive limit of two-level dynamics. Resonant Driving analyzes detuning, finite-pulse linewidth, and the long-time breakdown of first order. Linear Response Preview reorganizes the same transition frequencies and matrix elements into a susceptibility. The coherent two-level version is introduced in Rabi Oscillations: First Encounter, with the explicit perturbative match derived in Rabi Formula in the Weak-Drive Limit.

Write the perturbation as

V(t)=We−iωt+W†eiωt,V(t) = W e^{-i\omega t} + W^\dagger e^{i\omega t},

so that V(t)V(t) is Hermitian. The operator WW carries the part of the drive that oscillates with frequency ω\omega.

Let

ωfi=Ef−Eiℏ.\omega_{fi} = \frac{E_f-E_i}{\hbar}.

The first-order transition amplitude is

cf(1)(T)=−iℏ∫0Tdt eiωfit⟨f∣V(t)∣i⟩.c_f^{(1)}(T) = - \frac{i}{\hbar} \int_0^T dt\, e^{i\omega_{fi}t} \langle f\vert V(t)\vert i\rangle.

Substituting the harmonic perturbation gives

cf(1)(T)=−iℏ[Wfi∫0Tdt ei(ωfi−ω)t+Wfi†∫0Tdt ei(ωfi+ω)t],\begin{aligned} c_f^{(1)}(T) = - \frac{i}{\hbar} \bigg[ &W_{fi} \int_0^T dt\,e^{i(\omega_{fi}-\omega)t} \\ &+ W^\dagger_{fi} \int_0^T dt\,e^{i(\omega_{fi}+\omega)t} \bigg], \end{aligned}

where

Wfi=⟨f∣W∣i⟩,Wfi†=⟨f∣W†∣i⟩.W_{fi}=\langle f\vert W\vert i\rangle, \qquad W^\dagger_{fi}=\langle f\vert W^\dagger\vert i\rangle.

For any detuning Δ\Delta,

∫0Tdt eiΔt=eiΔT/2Tsin⁡(ΔT/2)ΔT/2.\int_0^T dt\,e^{i\Delta t} = e^{i\Delta T/2} T \frac{\sin(\Delta T/2)}{\Delta T/2}.

Thus

cf(1)(T)=−iTℏ[Wfiei(ωfi−ω)T/2sinc⁡((ωfi−ω)T2)+Wfi†ei(ωfi+ω)T/2sinc⁡((ωfi+ω)T2)],c_f^{(1)}(T) = - \frac{iT}{\hbar} \left[ W_{fi} e^{i(\omega_{fi}-\omega)T/2} \operatorname{sinc} \left( \frac{(\omega_{fi}-\omega)T}{2} \right) + W^\dagger_{fi} e^{i(\omega_{fi}+\omega)T/2} \operatorname{sinc} \left( \frac{(\omega_{fi}+\omega)T}{2} \right) \right],

where

sinc⁡x=sin⁡xx.\operatorname{sinc}x=\frac{\sin x}{x}.

The amplitude is largest when one of the oscillatory phases is nearly stationary.

If Ef>EiE_f>E_i, then ωfi>0\omega_{fi}>0. The term with denominator ωfi−ω\omega_{fi}-\omega is resonant when

ω≈ωfi.\omega\approx\omega_{fi}.

This is the absorption condition: the drive supplies energy ℏω\hbar\omega.

If Ef<EiE_f<E_i, then ωfi<0\omega_{fi}<0. The term with denominator ωfi+ω\omega_{fi}+\omega is resonant when

ω≈−ωfi.\omega\approx-\omega_{fi}.

This is the stimulated-emission condition: the system loses energy ℏω\hbar\omega to the drive.

The labels “absorption” and “emission” depend on which transition is being considered and how the external field is modeled. The mathematical test is always the phase-matching condition in the time integral.

Rotating and Counter-Rotating Contributions

Section titled “Rotating and Counter-Rotating Contributions”

Near a positive-frequency transition ωfi≈ω\omega_{fi}\approx\omega, the We−iωtW e^{-i\omega t} term varies slowly in the interaction picture. The W†eiωtW^\dagger e^{i\omega t} term varies with frequency approximately 2ω2\omega and often averages out.

Keeping only the slowly rotating term is the rotating-wave approximation. It is not part of first-order perturbation theory itself; it is an additional approximation based on near resonance and weak drive. It fails for strong drives, very short pulses, or regimes where counter-rotating effects accumulate.

The first-order probability is

Pi→f(T)≈∣cf(1)(T)∣2.P_{i\to f}(T) \approx \lvert c_f^{(1)}(T)\rvert^2.

Near a single resonance, this becomes approximately

Pi→f(T)≈∣Wfi∣2T2ℏ2sinc⁡2((ωfi−ω)T2).P_{i\to f}(T) \approx \frac{\lvert W_{fi}\rvert^2T^2}{\hbar^2} \operatorname{sinc}^2 \left( \frac{(\omega_{fi}-\omega)T}{2} \right).

The finite-time line shape has width of order 1/T1/T. In the long-time continuum limit, this finite-time peak becomes the delta-function structure behind Fermi’s golden rule.

The drive can be perfectly resonant and still produce no transition if the matrix element vanishes:

Wfi=0.W_{fi}=0.

Such zeros may come from parity, angular momentum, spin, translation symmetry, or other quantum numbers. The symmetry logic belongs in Why Symmetry Matters and Parity; this page uses those zeros inside the transition amplitude.

In the electric dipole approximation, a classical monochromatic electric field gives an interaction of the schematic form

V(t)=−d⋅E0cos⁡ωt,V(t) = - \mathbf d\cdot\mathbf E_0\cos\omega t,

where d\mathbf d is the dipole operator. Writing the cosine as positive- and negative-frequency parts puts the perturbation in harmonic form. Transition strengths are controlled by dipole matrix elements such as

⟨f∣d⋅ϵ∣i⟩,\langle f\vert\mathbf d\cdot\boldsymbol\epsilon\vert i\rangle,

where ϵ\boldsymbol\epsilon is the polarization direction.

  • Calling a transition resonant while its matrix element is zero by symmetry.
  • Forgetting the counter-rotating term before justifying a rotating-wave approximation.
  • Treating the finite-time sinc peak as an exact delta function at finite TT.
  • Letting first-order probabilities grow beyond the perturbative regime.
  • Confusing angular frequency ω\omega with ordinary frequency.
  1. Derive the finite-time integral for a detuning Δ\Delta.
Solution

Compute

∫0Tdt eiΔt=eiΔT−1iΔ.\int_0^T dt\,e^{i\Delta t} = \frac{e^{i\Delta T}-1}{i\Delta}.

Factor out eiΔT/2e^{i\Delta T/2}:

eiΔT−1=eiΔT/2(eiΔT/2−e−iΔT/2)=2ieiΔT/2sin⁡ΔT2.e^{i\Delta T}-1 = e^{i\Delta T/2} \left( e^{i\Delta T/2}-e^{-i\Delta T/2} \right) = 2i e^{i\Delta T/2}\sin\frac{\Delta T}{2}.

Therefore

∫0Tdt eiΔt=eiΔT/2Tsin⁡(ΔT/2)ΔT/2.\int_0^T dt\,e^{i\Delta t} = e^{i\Delta T/2} T \frac{\sin(\Delta T/2)}{\Delta T/2}.
  1. For Ef>EiE_f>E_i, which term is resonant when ω>0\omega>0?
Solution

If Ef>EiE_f>E_i, then ωfi>0\omega_{fi}>0. The denominator ωfi−ω\omega_{fi}-\omega can vanish when ω=ωfi\omega=\omega_{fi}. The denominator ωfi+ω\omega_{fi}+\omega is then far from zero for positive ω\omega. Thus the term proportional to WfiW_{fi} is resonant.

  1. Why does a zero matrix element suppress a resonant transition?
Solution

The resonant enhancement multiplies the matrix element. If Wfi=0W_{fi}=0 by symmetry, the resonant term in the amplitude vanishes. Resonance can enhance an allowed transition, but it cannot make a forbidden matrix element nonzero unless the symmetry assumptions are broken.

  • C. Cohen-Tannoudji, B. Diu, and F. Laloë, Quantum Mechanics, Wiley, 1977.
  • J. J. Sakurai and J. Napolitano, Modern Quantum Mechanics, 3rd ed., Cambridge University Press, 2020.
  • R. Shankar, Principles of Quantum Mechanics, 2nd ed., Springer, 1994.
  • L. Allen and J. H. Eberly, Optical Resonance and Two-Level Atoms, Dover, 1987.