Skip to content

A harmonic perturbation resolves a drive into positive- and negative-frequency components. Their phases distinguish absorption from stimulated emission, while their matrix elements decide which transitions are allowed. The finite-time calculation below retains both components before making a near-resonance approximation.

Required background. First-Order Transition Probability supplies the amplitude integral. Selection Rules explains symmetry restrictions; Fermi’s Golden Rule treats the continuum 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. Take ω>0\omega>0. The operator WW multiplies the negative exponential and need not be Hermitian. In particular, Wfi†=Wif∗W^\dagger_{fi}=W_{if}^*, which is generally different from Wfi∗W_{fi}^*. For a real cosine drive Acos⁡ωtA\cos\omega t with Hermitian AA, W=A/2W=A/2.

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],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],

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}\bigg[ 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) \bigg].

where

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

For a given nonzero matrix element, a nearly stationary phase permits coherent accumulation. For short pulses, both frequency components can contribute appreciably and must be added as complex amplitudes before squaring.

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. Parity and Selection Rules in Transition Rates derive these zeros. A forbidden first-order matrix element need not forbid a higher-order multiphoton transition.

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. For linear polarization with peak field E0E_0, W=−E0d⋅ϵ/2W=-E_0\mathbf d\cdot\boldsymbol\epsilon/2. Thus the resonant two-level Rabi angular frequency is 2∣Wfi∣/ℏ=E0∣dfi⋅ϵ∣/ℏ2|W_{fi}|/\hbar=E_0|\mathbf d_{fi}\cdot\boldsymbol\epsilon|/\hbar. The factor of two comes from resolving a real cosine, not from a change of physical field strength.

For H0=ℏω0σz/2H_0=\hbar\omega_0\sigma_z/2, compare W=gσ+W=g\sigma_+ with W=gσxW=g\sigma_x, where gg is real. Which first-order frequency components can excite the lower state?

Solution

Since ⟨+∣σ+∣−⟩=1\langle+|\sigma_+|-\rangle=1 and ⟨+∣σ−∣−⟩=0\langle+|\sigma_-|-\rangle=0, the circular drive has only the excitation amplitude proportional to g∫ei(ω0−ω)tdtg\int e^{i(\omega_0-\omega)t}dt. For W=gσxW=g\sigma_x, both W+−W_{+-} and W+−†W^\dagger_{+-} equal gg, so both difference- and sum-frequency integrals occur. The latter Hamiltonian is 2gcos⁡ωt σx2g\cos\omega t\,\sigma_x; its peak transverse coefficient is twice gg.

For the linear drive above, take T≪1/(ω0+ω)T\ll1/(\omega_0+\omega). Compare the full first-order excitation probability with the result obtained by dropping the counter-rotating term.

Solution

Both integrals approach TT, so c+(1)≃−2igT/ℏc_+^{(1)}\simeq-2igT/\hbar and P≃4g2T2/ℏ2P\simeq4g^2T^2/\hbar^2. Keeping only one component gives g2T2/ℏ2g^2T^2/\hbar^2. There has been no time to average the rapid component; “weak drive” alone does not justify the rotating-wave approximation on this time scale.

  • L. Allen and J. H. Eberly, Optical Resonance and Two-Level Atoms, Dover, 1987.
  • 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.