Past exam of the atomic and laser physics course of the University of Oxford / 2022 / 3

Choose the relative phases of the two energy eigenstates so that $d=⟨{\psi }_2|{\hat{D}}_0|{\psi }_1⟩$ is real and nonnegative. Define the resonant Rabi angular frequency $\Omega =d/\hslash$, the Rabi detuning $\delta ={\omega }_0-\omega$, and $\Sigma ={\omega }_0+\omega$. Substitution into the Schrödinger equation, followed by projection onto the two eigenstates, gives the exact interaction picture equationsKeep both probability amplitudes dynamical: replacing $c_1$ by one would incorrectly impose weak excitation and miss complete population transfer.

For $|\delta |,\Omega \ll {\omega }_0$, apply the rotating-wave approximation to discard the rapidly oscillating terms at $\Sigma$. Differentiating the second remaining equation and using the first givesThe initial conditions are $c_2(0)=0$ and ${\dot{c}}_2(0)=-i\Omega /2$. With ${\Omega }_R=\sqrt{{\Omega }^2+{\delta }^2}$, their solution isConsequently,This result permits order-one excitation, including full population transfer on resonance, even though the drive is small compared with the optical angular frequency.

For an average from switch-on to time $\tau$, integrate the squared sine:On resonance this limiting average is $1/2$, although the instantaneous probability repeatedly reaches one. The zero-coupling, zero-Rabi detuning case is understood separately as $P_2=0$.

Near resonance alone does not make this an exact result for arbitrarily strong driving. The stated sinusoidal Hamiltonian also needs the rotating-wave approximation. A concrete counterexample is $\omega ={\omega }_0$ and $\Omega ={\omega }_0$: the exact equations give $P_2(t)={\Omega }^2{t}^2+O(t^4)$ at switch-on, while the claimed rotating-wave sinusoid gives ${\Omega }^2{t}^2/4+O(t^4)$. Thus “without assuming weak radiation” can mean without assuming small excitation, but cannot remove the drive-to-carrier scale separation. The infinite-time average above belongs to the rotating-wave model; omitted corrections can accumulate over sufficiently long times.

Rabi oscillations require a coherent drive with a stable phase over several population cycles and an isolated, appreciably coupled pair of quantum states. For large-amplitude oscillations, the Rabi detuning must be comparable to or smaller than the Rabi angular frequency. The population lifetime, optical coherence time, transit time and drive coherence time should all exceed $2\pi /{\Omega }_R$ by enough to resolve the cycles. An ensemble also needs sufficiently uniform Rabi frequencies and Rabi detunings that phase averaging does not wash the oscillations out. The usual sinusoidal formula additionally assumes the rotating-wave approximation.

The Einstein coefficients describe incoherent transition rates, rather than the evolving optical probability amplitudes. Their simplest use is appropriate for broadband, randomly phased radiation, with field correlation time short relative to population evolution, or for rapidly dephased chemical atoms whose optical coherences can be eliminated. The radiation energy density per unit frequency must be smooth over the relevant spectral line. Stimulated upward and downward transitions and spontaneous emission then enter rate equations, without predicting coherent population return. A short observation time can give a small excitation probability even for a coherent field; that condition alone does not justify an Einstein coefficients rate description. Coherent laser irradiation can also admit steady rate equations after optical coherence has relaxed, with the appropriate saturation response.

Let the upper-state population lifetime be $T_1=1/\Gamma$, and take the lower energy level to be stable. Spontaneous emission gives survival probability $e^{-\Gamma t}$, hence a field or excited-state probability amplitude proportional to $e^{-\Gamma t/2}{e}^{-i{\omega }_{0}t}$ for $t\ge 0$. Its Fourier transform at angular frequency $\omega$ is proportional toThe squared modulus therefore has the Lorentzian line shapeHere $L$ is normalized over the full mathematical frequency axis. The full width at half maximum in ordinary frequency is $1/(2\pi T_1)$. If both energy levels decay, their population decay rates add in the natural linewidth; pure dephasing adds to the amplitude-decay rate instead.

Impact-regime pressure broadening provides another example. Independent phase-randomizing collisions with rate ${\gamma }_c$ give optical correlation proportional to $e^{-{\gamma }_c|t|}$, whose Fourier transform is a Lorentzian line shape of full width at half maximum $2{\gamma }_c$. This assumes short, uncorrelated impacts: a general static distribution of perturbing fields need not give a Lorentzian line shape.

Use first-order time-dependent perturbation theory, setting $c_1=1$ in the equation for $c_2$ and retaining the generally complex matrix element $d=⟨{\psi }_2|{\hat{D}}_0|{\psi }_1⟩$. With $\delta ={\omega }_0-\omega$ and $\Sigma ={\omega }_0+\omega$, the probability amplitude isFor $t>0$, split the integral at the pulse maximum. For either $x=\delta$ or $x=\Sigma$,Thus the full first-order answer, without discarding the counter-rotating contribution, isAt long times it givesThis full expression is not exactly one squared Lorentzian line shape. To obtain the intended near-resonance spectral line profile, additionally assume $b,|\delta |\ll {\omega }_0$, so the pulse has a narrow positive-frequency lobe and $F(\Sigma ,t)$ is negligible. ThenThe weak-excitation condition near resonance is $|d|/(\hslash b)\ll 1$. The squared Lorentzian line shape occurs because the symmetric exponential pulse has a Lorentzian line shape amplitude Fourier transform, which is then squared to obtain probability.

At half the maximum, $[1+(\delta /b{)}^2{]}^2=2$, soThe full width at half maximum refers to varying the drive frequency through the positive-frequency resonance. For a concrete failure of the unqualified squared-Lorentzian line shape claim, set $b=\omega ={\omega }_0$ while keeping $|d|/(\hslash b)\ll 1$. The full first-order resonant probability is $36/25$ times the one-lobe prediction. Weak radiation alone therefore does not justify dropping the second Fourier transform lobe.