Rabi Oscillations

Point a clean, coherent laser at a two-level atom and turn up the amplitude, and the atom does not simply absorb a photon and stay lit. It flickers. The probability of finding it in the excited state climbs to one, falls back to zero, and climbs again — a rhythmic sloshing of population between ground and excited states that repeats for as long as the coherence survives. This is a Rabi oscillation, and it is the most literal picture we have of what "light driving an atom" actually means. It takes the transition dipole moment and the electric-dipole coupling −μ⋅E-\boldsymbol\mu\cdot\mathbf E from the two articles before it and asks the sharpest possible question: not how often does the atom jump, but how does its quantum amplitude evolve, instant by instant, while the field is on. The answer overturns the rate-equation intuition inherited from Einstein's AA and BB coefficients. Rate equations describe a one-way drift toward equilibrium; Rabi oscillations are reversible — during one cycle the atom borrows a photon's energy from the field and hands it straight back by stimulated emission. Only when something breaks the phase memory (spontaneous emission, a collision, laser jitter) does the exchange become one-way and the familiar absorption line reappear. Rabi oscillation is the parent phenomenon; incoherent absorption is its damped, forgetful child. Governing everything is one frequency, the Rabi frequency Ω=μE0/ℏ\Omega=\mu E_0/\hbar, set by the dipole matrix element and the field amplitude — proportional to I\sqrt{I}, not intensity, the fingerprint of a coherent amplitude effect. The full story with damping is the optical Bloch equations, and the resonant-absorption line shape, the natural linewidth, and every qubit gate downstream are its consequences.

Key takeaways
  • A two-level atom driven by a near-resonant coherent field does not absorb once and stay excited — it cycles population reversibly between ground and excited states at the Rabi frequency.
  • On resonance the excited probability is Pe(t)=sin⁡2(Ωt/2)P_e(t)=\sin^2(\Omega t/2), swinging fully between 0 and 1, with Ω=μE0/ℏ\Omega=\mu E_0/\hbar.
  • Off resonance the oscillation is faster but incomplete: it runs at ΩR=Ω2+δ2\Omega_R=\sqrt{\Omega^2+\delta^2} and peaks at only Ω2/ΩR2<1\Omega^2/\Omega_R^2<1.
  • A π\pi pulse fully inverts the atom; a π/2\pi/2 pulse builds an equal superposition — the primitives of qubit control and Ramsey spectroscopy.
  • Because Ω∝I\Omega\propto\sqrt I the Rabi rate is amplitude-linear, distinct from an intensity-linear absorption rate; decoherence damps the oscillation into the rate-equation steady state.
  • Key idea: a near-resonant coherent field turns a two-level atom into a driven two-state oscillator whose population is exchanged reversibly at ΩR=Ω2+δ2\Omega_R=\sqrt{\Omega^2+\delta^2}, with the coupling Ω=μE0/ℏ\Omega=\mu E_0/\hbar fixed by the transition dipole moment and the field amplitude.
The core insight

A near-resonant coherent field turns a two-level atom into a driven two-state oscillator: population is exchanged reversibly at the generalised Rabi frequency ΩR=Ω2+δ2\Omega_R=\sqrt{\Omega^2+\delta^2}, with the coupling strength Ω=μE0/ℏ\Omega=\mu E_0/\hbar set by the transition dipole moment and the field amplitude.

Where this lives today

The Rabi oscillation became a control knob. In a trapped-ion, neutral-atom, or superconducting processor, coherent control is implemented with calibrated drive pulses: an ideal resonant π\pi pulse implements a population inversion (a NOT operation in the appropriate qubit basis), a π/2\pi/2 pulse prepares an equal superposition, and the pulse phase selects the rotation axis. Gate fidelity depends in part on how firmly Ω≫Γ\Omega\gg\Gamma — how many coherent cycles fit before relaxation interrupts them — but also on detuning, pulse-shape errors, dephasing, and other control imperfections. Engineers watch live "Rabi flops" on a scope while tuning the drive, reading Ω\Omega off the flopping period. The same physics runs the world's clocks. Ramsey's two-π/2\pi/2-pulse sequence, with a long free-evolution gap, interrogates the atomic phase and gives fringes far narrower than a single Rabi resonance — the heart of caesium fountains and optical-lattice clocks stable to 10−1810^{-18}, with detuning entering through the generalised Rabi frequency of equation (4). Strong resonant driving splits levels by Ω\Omega into an Autler–Townes doublet; far-detuned driving shifts them by Ω2/4δ\Omega^2/4\delta, the light shift that traps atoms in optical tweezers and lattices; and a coherent 2π2\pi pulse can slide through a resonant medium without net absorption (self-induced transparency). From atomic clocks to superconducting qubits, the same two-level coherent-control physics remains a central building block of modern quantum technology.

The physical picture

The cleanest way to see a Rabi oscillation is to stop thinking about photons and picture a spinning arrow. Any two-level system — ground ∣g⟩|g\rangle, excited ∣e⟩|e\rangle — maps exactly onto a spin-1/2, and its state maps onto a unit vector on the Bloch sphere: south pole "definitely ground", north pole "definitely excited", every point between a coherent superposition with a definite phase. The near-resonant field, in the frame that rotates with it, acts like a static effective field — or equivalently an effective rotation vector — and the Bloch state vector precesses about it. That precession is the Rabi oscillation, and the whole phenomenon reduces to one question: where does the torque axis point? On resonance, the effective rotation axis Ωeff\boldsymbol{\Omega}_{\rm eff} lies in the equatorial plane. A vector starting at the south pole precesses in a great circle straight over the top — ground to excited and back, a complete sweep with nothing out of reach. That is why resonant driving gives complete inversion. Detune the field and the axis tilts up out of the plane; the effective field grows longer, ΩR=Ω2+δ2\Omega_R=\sqrt{\Omega^2+\delta^2}, so precession speeds up, but the state now circles a tilted axis on a cone that never reaches the north pole. Faster but incomplete — the two signatures of off-resonant driving in one picture.

Figure
The Bloch-sphere torque picture. On resonance δ = 0 the drive torque Ω lies in the equatorial plane and sweeps the state vector from the ground pole g up to the excited pole e — a complete Rabi cycle. Detuned by δ = Ω the effective field tilts 45° and the vector precesses on a cone that reaches only the equator, so P_e peaks at just 0.5
One geometric question: where does the drive torque point? resonant δ = 0 excited e ground g full sweep Ω P_e reaches 1 detuned δ = Ω excited e ground g Ω_R 45° reaches only the equator: P_e,max = 0.5

Two facts fall straight out. First, in the ideal closed two-level model, energy is conserved: as the Bloch vector climbs, the atom gains energy from the driving field, and as it falls it returns that energy to the field by stimulated emission — the energy exchange is reversible over a coherent Rabi cycle. Second, the coupling scales with field amplitude, Ω=μE0/ℏ∝I\Omega=\mu E_0/\hbar\propto\sqrt I: for a fixed mode and polarization, doubling the intensity makes Ω\Omega rise by only 2≈1.41\sqrt2\approx1.41 (41%), rather than by a factor of two as an intensity-linear rate would — a sharp distinction between coherent Rabi dynamics and incoherent rate-equation dynamics. The second figure plots the excited probability against the dimensionless clock Ωt\Omega t. Read off the two lessons: the resonant curve reaches the ceiling Pe=1P_e=1 (a π\pi pulse), while the detuned curve is faster and clipped at 0.50.5.

Figure
Excited-state probability P_e versus the dimensionless time Ωt. On resonance green it is sin²(Ωt/2), swinging fully between 0 and 1 with a π pulse at Ωt = π. Detuned by δ = Ω amber it runs faster at Ω_R = √2 Ω yet reaches only P_e,max = Ω²/(Ω²+δ²) = 0.5
Resonant driving reaches P_e = 1; detuning is faster but clipped at 0.5 1.0 0.5 0 0 π 2π 3π 4π dimensionless time Ωt excited probability P_e π pulse: P_e = 1 δ = Ω: peak 0.5 at Ω_R = √2 Ω

The core equations

Everything starts from a two-level atom coupled to a classical field through the electric-dipole interaction. The bare atom carries a splitting ℏω0\hbar\omega_0 between ∣g⟩|g\rangle and ∣e⟩|e\rangle; the field adds the dipole coupling −μ⋅E(t)-\boldsymbol\mu\cdot\mathbf E(t) recalled from the electric-dipole-approximation article. This Hamiltonian is the entire input — no relaxation yet, just a coherent drive.

H^=ℏω0 ∣e⟩⟨e∣−μ⋅E(t)(1)\hat H = \hbar\omega_0\,|e\rangle\langle e| - \boldsymbol{\mu}\cdot\mathbf{E}(t) \qquad (1)

where: ω0\omega_0 — atomic transition angular frequency; ∣e⟩,∣g⟩|e\rangle,|g\rangle — excited and ground states; μ\boldsymbol\mu — atomic electric-dipole operator; E(t)=E0cos⁡ωt e^\mathbf E(t)=E_0\cos\omega t\,\hat{\mathbf e} — classical driving field of amplitude E0E_0 and polarization e^\hat{\mathbf e}. The reading: a two-state energy ladder shaken by an oscillating dipole force. The strength of that shaking is one number, the Rabi frequency: the dipole matrix element along the field, times the amplitude, over ℏ\hbar. Being linear in E0E_0 it is linear in I\sqrt I, and it has units of angular frequency because it is the rate at which the driven atom cycles.

Ω≡μ E0ℏ,μ=⟨e∣ e^⋅μ ∣g⟩(2)\Omega \equiv \frac{\mu\, E_0}{\hbar}, \qquad \mu = \langle e|\,\hat{\mathbf e}\cdot\boldsymbol{\mu}\,|g\rangle \qquad (2)

where: Ω\Omega — on-resonance Rabi angular frequency; μ\mu — transition dipole matrix element along the field; E0E_0 — peak field amplitude; ℏ\hbar — reduced Planck constant. The reading: a bigger dipole or stronger field cycles the atom faster, and the dependence on E0E_0 (not E02E_0^2) is what makes Rabi coherent, not rate-like. Substituting ∣ψ⟩=cg∣g⟩+cee−iω0t∣e⟩|\psi\rangle=c_g|g\rangle+c_e e^{-i\omega_0 t}|e\rangle into the Schrödinger equation and dropping the fast counter-rotating term (the rotating-wave approximation, derived below) collapses the dynamics to two coupled amplitude equations — the engine of the whole phenomenon.

ic˙e=−Ω2 e−iδt cg,ic˙g=−Ω2 e+iδt ce(3)i\dot c_e = -\tfrac{\Omega}{2}\,e^{-i\delta t}\,c_g,\qquad i\dot c_g = -\tfrac{\Omega}{2}\,e^{+i\delta t}\,c_e \qquad (3)

where: ce,cgc_e,c_g — excited and ground amplitudes; δ=ω−ω0\delta=\omega-\omega_0 — detuning from resonance; Ω\Omega — Rabi frequency. The reading: each amplitude drives the other at rate Ω/2\Omega/2, and detuning appears only as a rotating phase — the seed of both the off-resonant speed-up and its incompleteness. Solving the coupled pair (Derivation C) gives an oscillation whose frequency is not Ω\Omega but the length of the effective-field vector (Ω,0,−δ)(\Omega,0,-\delta) — the generalised Rabi frequency, the key off-resonant quantity.

ΩR=Ω2+δ2(4)\Omega_R = \sqrt{\Omega^2 + \delta^2} \qquad (4)

where: ΩR\Omega_R — generalised Rabi frequency; Ω\Omega — on-resonance value; δ\delta — detuning. The reading: detuning adds in quadrature, so the atom always cycles at least as fast as Ω\Omega, faster the further you tune away. Frequency and amplitude together give the Rabi formula — the exact excited-state probability for a ground-state atom at any detuning, the master result of the article.

Pe(t)=Ω2Ω2+δ2 sin⁡2 ⁣(ΩR t2)(5)P_e(t) = \frac{\Omega^2}{\Omega^2+\delta^2}\,\sin^2\!\left(\frac{\Omega_R\,t}{2}\right) \qquad (5)

where: Pe(t)P_e(t) — excited occupation; prefactor Ω2/(Ω2+δ2)\Omega^2/(\Omega^2+\delta^2) — the maximum transfer, unity on resonance; ΩR\Omega_R — sets the period. The reading: a fixed-shape sin⁡2\sin^2 oscillation whose height is throttled by detuning and whose speed is boosted by it. For an ideal resonant pulse with a fixed phase, what matters is not the peak Rabi frequency but its time integral, the pulse area — the one scalar that determines the population transfer, so different pulse shapes with the same area give the same rotation angle.

θ=∫0τΩ(t) dt(6)\theta = \int_0^{\tau}\Omega(t)\,\mathrm{d}t \qquad (6)

where: θ\theta — pulse area; τ\tau — duration; Ω(t)\Omega(t) — instantaneous Rabi frequency. The reading: θ=π\theta=\pi fully inverts, θ=π/2\theta=\pi/2 makes an equal superposition, θ=2π\theta=2\pi returns to ground — the alphabet of coherent control. Real atoms decohere, and adding relaxation at rate Γ\Gamma damps the oscillation to a steady state. Balancing drive against decay gives the saturated population — the bridge from Rabi dynamics to the ordinary absorption line.

Pess=12 Ω2/2δ2+Γ2/4+Ω2/2=12 s1+s,s=Ω2/2δ2+Γ2/4(7)P_e^{ss} = \frac{1}{2}\,\frac{\Omega^2/2}{\delta^2 + \Gamma^2/4 + \Omega^2/2} = \frac{1}{2}\,\frac{s}{1+s},\quad s=\frac{\Omega^2/2}{\delta^2+\Gamma^2/4} \qquad (7)

where: Γ\Gamma — total relaxation rate; ss — saturation parameter; PessP_e^{ss} — steady excited population. The reading: as s→∞s\to\infty it saturates at exactly 1/21/2 — stimulated emission balances absorption, and the atom can never be more than half-excited on average once coherence is gone.

Derivation AFrom the Schrödinger equation to the rotating-wave amplitude equations

The goal is to turn the dipole Hamiltonian of equation (1) into the two coupled amplitude equations (3) that generate Rabi oscillations, and to expose the one approximation — the rotating-wave approximation — on which the whole clean picture rests. We move into the frame co-rotating with the bare atomic phases and identify which of the two driving terms can be discarded and why, leaving a pair simple enough to solve in closed form while keeping all the physics that matters at optical frequencies.

Step 1. With the fast atomic phase factored out, ∣ψ⟩=cg∣g⟩+cee−iω0t∣e⟩|\psi\rangle=c_g|g\rangle+c_e e^{-i\omega_0 t}|e\rangle, substitute into iℏ ∂t∣ψ⟩=H^∣ψ⟩i\hbar\,\partial_t|\psi\rangle=\hat H|\psi\rangle with H^int=−μE0cos⁡(ωt)(∣e⟩⟨g∣+∣g⟩⟨e∣)\hat H_{int}=-\mu E_0\cos(\omega t)(|e\rangle\langle g|+|g\rangle\langle e|). Projecting onto ⟨e∣\langle e| and writing cos⁡ωt=12(eiωt+e−iωt)\cos\omega t=\tfrac12(e^{i\omega t}+e^{-i\omega t}) gives ic˙e=−Ω2(ei(ω−ω0)t+e−i(ω+ω0)t)cgi\dot c_e=-\tfrac{\Omega}{2}(e^{i(\omega-\omega_0)t}+e^{-i(\omega+\omega_0)t})c_g. Why: factoring e−iω0te^{-i\omega_0 t} moves us to the interaction picture, where only the slow drive–atom mismatch survives explicitly.

Step 2. Split the two exponentials by frequency. The first, eiδte^{i\delta t}, oscillates slowly at the detuning; the second, e−i(ω+ω0)te^{-i(\omega+\omega_0)t}, oscillates at ≈2ω0\approx2\omega_0 (hundreds of terahertz) and averages to nothing over the slow timescale 1/Ω1/\Omega. Why: a term that flips sign a million times per Rabi cycle contributes zero net drive.

Step 3. Discard the fast term. This is the rotating-wave approximation, valid whenever Ω,∣δ∣≪ω0\Omega,|\delta|\ll\omega_0, and it leaves exactly the coupled pair of equation (3). The neglected term is not truly zero — it produces the small Bloch–Siegert shift of order Ω2/4ω0\Omega^2/4\omega_0, measurable in strongly driven microwave and circuit-QED systems. Why: the RWA is a controlled truncation, not a fudge, and it fails in the ultrastrong-coupling regime where Ω\Omega approaches ω0\omega_0.

Go deeper · Derivation B · 1 minThe on-resonance solution complete inversion and pulse areas
Derivation BThe on-resonance solution complete inversion and pulse areas

The goal is to solve the amplitude equations on exact resonance and read off the two results that make Rabi oscillations useful: complete, deterministic inversion, and the π\pi and π/2\pi/2 pulses that are the primitive gates of coherent control. Setting δ=0\delta=0 removes the rotating phases and reduces the coupled system to a single harmonic-oscillator equation whose conclusions run all the way to quantum computing.

Step 1. On resonance δ=0\delta=0, equation (3) becomes ic˙e=−Ω2cgi\dot c_e=-\tfrac{\Omega}{2}c_g and ic˙g=−Ω2cei\dot c_g=-\tfrac{\Omega}{2}c_e. Differentiate the first and substitute the second to decouple them: c¨e=−(Ω/2)2ce\ddot c_e=-(\Omega/2)^2 c_e. This is simple harmonic motion in the amplitude at frequency Ω/2\Omega/2. Why: eliminating cgc_g turns two coupled equations into one textbook oscillator, and the appearance of Ω/2\Omega/2 — half the Rabi frequency — in the amplitude is exactly why the probability, which squares it, oscillates at Ω\Omega.

Step 2. With ce(0)=0c_e(0)=0, cg(0)=1c_g(0)=1 the solution is ce(t)=isin⁡(Ωt/2)c_e(t)=i\sin(\Omega t/2), cg(t)=cos⁡(Ωt/2)c_g(t)=\cos(\Omega t/2), so Pe(t)=∣ce∣2=sin⁡2(Ωt/2)P_e(t)=|c_e|^2=\sin^2(\Omega t/2), and probability is conserved: ∣ce∣2+∣cg∣2=1|c_e|^2+|c_g|^2=1. Why: unitary evolution pins total probability to one, so population genuinely sloshes between the states rather than leaking away — the signature of coherent, lossless dynamics.

Step 3. Evaluate at special times. At Ωt=π\Omega t=\pi, Pe=1P_e=1: a π\pi pulse completely inverts the atom. At Ωt=π/2\Omega t=\pi/2, Pe=1/2P_e=1/2: a π/2\pi/2 pulse prepares the equal superposition 12(∣g⟩+i∣e⟩)\frac{1}{\sqrt2}(|g\rangle+i|e\rangle) used in Ramsey interferometry. At Ωt=2π\Omega t=2\pi, the atom returns to the ground state with the state vector acquiring an overall minus sign, reflecting the spin-1/21/2 nature of the quantum state under a 2π2\pi rotation. Why: the outcome depends only on the accumulated area θ=∫Ω dt\theta=\int\Omega\,dt, so calibrating the area gives robust gates independent of pulse shape.

Go deeper · Derivation C · 1 minThe generalised Rabi formula and the damped steady state
Derivation CThe generalised Rabi formula and the damped steady state

The goal is twofold: first to solve the amplitude equations at arbitrary detuning for the master Rabi formula (5), showing how detuning both speeds up and caps the oscillation; second to add relaxation and watch it decay into the incoherent rate-equation steady state — showing that coherent cycling and irreversible absorption are the same physics on two timescales.

Step 1. Off resonance, eliminate cgc_g from equation (3) to get c¨e−iδc˙e+Ω24ce=0\ddot c_e-i\delta\dot c_e+\tfrac{\Omega^2}{4}c_e=0, with roots iδ/2±iΩR/2i\delta/2\pm i\Omega_R/2 and ΩR=Ω2+δ2\Omega_R=\sqrt{\Omega^2+\delta^2}. With ce(0)=0c_e(0)=0 the solution is ce∝eiδt/2sin⁡(ΩRt/2)c_e\propto e^{i\delta t/2}\sin(\Omega_R t/2), and ∣ce∣2|c_e|^2 gives Pe=Ω2Ω2+δ2sin⁡2 ⁣(ΩRt2)P_e=\dfrac{\Omega^2}{\Omega^2+\delta^2}\sin^2\!\big(\tfrac{\Omega_R t}{2}\big). Why: the detuning enters as a phase eiδt/2e^{i\delta t/2} that drops out of ∣ce∣2|c_e|^2 but survives inside ΩR\Omega_R — hence faster oscillation, throttled amplitude.

Step 2. Read the prefactor. On resonance it is 11 and PeP_e reaches the ceiling; for ∣δ∣≫Ω|\delta|\gg\Omega it collapses to Ω2/δ2≪1\Omega^2/\delta^2\ll1 — tiny fast wobbles underlying the AC Stark (light) shift. The maximum transfer Pe,max⁡=Ω2/(Ω2+δ2)P_{e,\max}=\Omega^2/(\Omega^2+\delta^2) is unity only on resonance. Why: this is the quantitative content of the Bloch-sphere cone — a tilted axis cannot reach the excited pole, and the shortfall is exactly δ2/ΩR2\delta^2/\Omega_R^2.

Step 3. Add relaxation. Recasting the amplitudes as the Bloch vector (u,v,w)(u,v,w), w=2Pe−1w=2P_e-1, with coherences decaying at Γ/2\Gamma/2 and inversion at Γ\Gamma, gives u˙=δv−Γ2u\dot u=\delta v-\tfrac{\Gamma}{2}u, v˙=−δu+Ωw−Γ2v\dot v=-\delta u+\Omega w-\tfrac{\Gamma}{2}v, w˙=−Ωv−Γ(w+1)\dot w=-\Omega v-\Gamma(w+1). For weak damping Γ≪Ω\Gamma\ll\Omega on resonance, Pe(t)≈12[1−e−3Γt/4cos⁡Ωt]P_e(t)\approx\tfrac12[1-e^{-3\Gamma t/4}\cos\Omega t]: many cycles before settling. Why: the envelope is the death of coherence — count surviving cycles and you measure Ω/Γ\Omega/\Gamma, the qubit figure of merit.

Step 4. Set the derivatives to zero for the steady state: Pess=12 s/(1+s)P_e^{ss}=\tfrac12\,s/(1+s), s=(Ω2/2)/(δ2+Γ2/4)s=(\Omega^2/2)/(\delta^2+\Gamma^2/4), saturating at 1/21/2 as s→∞s\to\infty, with power-broadened width Γ2+2Ω2\sqrt{\Gamma^2+2\Omega^2}. Why: in the overdamped limit Γ≫Ω\Gamma\gg\Omega the coherences follow adiabatically and PeP_e obeys a single rate equation — the Einstein picture with its Lorentzian line re-emerges as the memoryless limit of Rabi dynamics.

Limiting cases

The behaviour is set by three rates — coupling Ω\Omega, detuning δ\delta, relaxation Γ\Gamma — and each corner of that space has a clean, nameable limit.

CaseResultMeaning
Exact resonance, no dampingPe(t)=sin⁡2 ⁣(Ωt2)P_e(t)=\sin^2\!\left(\tfrac{\Omega t}{2}\right)Complete, undamped inversion — the ideal two-state oscillator; a π\pi pulse reaches Pe=1P_e=1.
Large detuning ∣δ∣≫Ω\lvert\delta\rvert\gg\OmegaPe,max⁡=Ω2Ω2+δ2→Ω2δ2≪1P_{e,\max}=\dfrac{\Omega^2}{\Omega^2+\delta^2}\to\dfrac{\Omega^2}{\delta^2}\ll1Fast, tiny wobbles; the atom stays essentially in the ground state — the regime of the AC Stark / light shift and dipole traps.
Strong damping Γ≫Ω\Gamma\gg\OmegaPe(t)→PessP_e(t)\to P_e^{ss} (overdamped)Oscillations vanish; monotonic approach to the rate-equation steady state — the Einstein absorption regime.
Weak damping Γ≪Ω\Gamma\ll\OmegaPe(t)≈12 ⁣[1−e−3Γt/4cos⁡Ωt]P_e(t)\approx\tfrac12\!\left[1-e^{-3\Gamma t/4}\cos\Omega t\right]Underdamped Rabi oscillations, many cycles visible, relaxing to Pess=1/2P_e^{ss}=1/2.
Strong saturation s→∞s\to\inftyPess→12P_e^{ss}\to\tfrac12Population caps at one-half; stimulated emission balances absorption — an atom cannot be driven fully excited on average.

Worked examples

1 · Rabi frequency and π-pulse time for a real atom (confidence builder)

A near-resonant laser drives an alkali D-line transition with dipole matrix element μ=2.5×10−29 C m\mu=2.5\times10^{-29}\ \mathrm{C\,m} (about 3 ea03\,ea_0, a realistic value) at intensity I=100 W/m2I=100\ \mathrm{W/m^2}. Find the field amplitude, the Rabi frequency in rad/s and in ordinary frequency, and the duration of a resonant π\pi pulse.

Step 1. Intensity to field: for a plane wave I=12cε0E02I=\tfrac12 c\varepsilon_0 E_0^2, so E0=2I/(cε0)=2(100)/(2.998×108⋅8.854×10−12)=7.54×104=274 V/mE_0=\sqrt{2I/(c\varepsilon_0)}=\sqrt{2(100)/(2.998\times10^8\cdot8.854\times10^{-12})}=\sqrt{7.54\times10^{4}}=\mathbf{274\ V/m} — the peak field the dipole feels.

Step 2. Rabi frequency from equation (2): Ω=μE0/ℏ=(2.5×10−29)(274.7)/(1.055×10−34)=6.51×107 rad/s\Omega=\mu E_0/\hbar=(2.5\times10^{-29})(274.7)/(1.055\times10^{-34})=6.51\times10^7\ {\rm rad/s}, i.e. Ω/2π=10.36 MHz\Omega/2\pi=10.36\ {\rm MHz}. This sits far below the transition frequency (∼5×1014\sim5\times10^{14} Hz), so the RWA is superb.

Step 3. A π\pi pulse needs Ωτ=π\Omega\tau=\pi, so τπ=π/Ω=48 ns\tau_\pi=\pi/\Omega=\mathbf{48\ ns}; a π/2\pi/2 pulse takes 24 ns\mathbf{24\ ns}. In that window the atom is carried deterministically from ground to fully excited.

Step 4. Scaling check: doubling the intensity to 200 W/m2200\ \mathrm{W/m^2} raises E0E_0 and Ω\Omega by 2\sqrt2 — from 10.4 to 14.6 MHz\mathbf{14.6\ MHz}, a 41 % increase, not the 100 % an intensity-linear rate would show. Result: E0=274E_0=274 V/m, Ω=6.5×107\Omega=6.5\times10^7 rad/s (10.4 MHz), τπ=48\tau_\pi=48 ns, Ω∝I\Omega\propto\sqrt I.

2 · How detuning throttles the transfer

The same atom is now driven off resonance. Compute the generalised Rabi frequency and the maximum excited population for detunings δ=Ω\delta=\Omega and δ=2Ω\delta=2\Omega, and say what happens to the oscillation period.

Step 1. For δ=Ω\delta=\Omega, equation (4) gives ΩR=2 Ω=1.41 Ω\Omega_R=\sqrt2\,\Omega=\mathbf{1.41\,\Omega} (41 % faster), and equation (5) gives Pe,max⁡=Ω2/(Ω2+δ2)=0.50P_{e,\max}=\Omega^2/(\Omega^2+\delta^2)=\mathbf{0.50} — the oscillation tops out at the equator, the clipped amber curve above.

Step 2. For δ=2Ω\delta=2\Omega, ΩR=5 Ω=2.24 Ω\Omega_R=\sqrt5\,\Omega=\mathbf{2.24\,\Omega} and Pe,max⁡=1/5=0.20P_{e,\max}=1/5=\mathbf{0.20}: twice the detuning more than doubles the speed while cutting the reach to a fifth.

Step 3. The trend: detuning always shortens the period (ΩR>Ω\Omega_R>\Omega) and lowers the ceiling (prefactor <1<1), the two locked together through the vector (Ω,0,−δ)(\Omega,0,-\delta). Result: δ=Ω⇒ΩR=1.41 Ω, Pe,max⁡=0.50\delta=\Omega\Rightarrow\Omega_R=1.41\,\Omega,\ P_{e,\max}=0.50; δ=2Ω⇒ΩR=2.24 Ω, Pe,max⁡=0.20\delta=2\Omega\Rightarrow\Omega_R=2.24\,\Omega,\ P_{e,\max}=0.20.

3 · Saturating a damped transition

Drive the same atom on resonance (δ=0\delta=0) with relaxation rate Γ\Gamma. Find the steady-state excited population at (a) the saturation point Ω2/2=Γ2/4\Omega^2/2=\Gamma^2/4 and (b) hard drive Ω2/2=100 Γ2/4\Omega^2/2=100\,\Gamma^2/4.

Step 1. On resonance s=(Ω2/2)/(Γ2/4)s=(\Omega^2/2)/(\Gamma^2/4), so case (a) s=1s=1 (drive equal to decay) and case (b) s=100s=100. From equation (7), Pess=12 s/(1+s)P_e^{ss}=\tfrac12\,s/(1+s): case (a) 0.25\mathbf{0.25}; case (b) 12(100/101)=0.495\tfrac12(100/101)=\mathbf{0.495} — within 1 % of the ceiling.

Step 2. However hard you drive, Pess→1/2P_e^{ss}\to1/2: once coherence is lost, stimulated emission empties ∣e⟩|e\rangle as fast as absorption fills it. This incoherent ceiling is what the coherent π\pi pulse of Example 1 briefly beats, reaching Pe=1P_e=1 before decoherence catches up. Result: s=1⇒0.25s=1\Rightarrow0.25; s=100⇒0.495s=100\Rightarrow0.495; ceiling 1/21/2.

Physics of the portrait

Rabi oscillations carry their discoverer's name because he found the resonance formula not in light but in magnetism — and the two turned out to be one equation.

Origins · 1 minIsidor Isaac Rabi and how this was discovered
II
Isidor Isaac Rabi
1898–1988 · molecular-beam magnetic resonance and the two-level rotation formula

Isidor Isaac Rabi was born in Rymanów, in what is now Poland, and raised in New York; at Columbia he built the molecular-beam laboratory that defined mid-century atomic physics. In 1937 he analysed a spin in a static magnetic field driven by a small rotating field at the resonant frequency, and derived the formula — maximum flip probability Ω2/(Ω2+δ2)\Omega^2/(\Omega^2+\delta^2), oscillation at Ω2+δ2\sqrt{\Omega^2+\delta^2} — that now bears his name. The next year his group (Rabi, Millman, Kusch, Zacharias) turned it into the molecular-beam magnetic-resonance method, measuring nuclear magnetic moments by sweeping a radio-frequency field across a beam and watching the resonant dip; for this he received the 1944 Nobel Prize. The formula was written for spins, but the mathematics is identical for a two-level atom driven through its transition dipole — which is why the coherent light–matter cycling here carries his name. His student Norman Ramsey extended the method into the separated-oscillatory-fields technique behind atomic clocks, and the same resonance logic runs through NMR, MRI, and every qubit driven by a calibrated pulse today.

History
1937
— Rabi derives the resonant spin-flip formula for a field rotating at the transition frequency
1938
— Rabi, Millman, Kusch and Zacharias demonstrate molecular-beam magnetic resonance
1944
— Rabi is awarded the Nobel Prize in Physics for the resonance method
1950
— Ramsey introduces the separated-oscillatory-fields method, sharpening the Rabi resonance
1969
— McCall and Hahn report self-induced transparency, a macroscopic 2π-pulse effect in optics

Going further

Texts. Foot, Atomic Physics, Ch. 7 (the interaction of atoms with radiation) · Loudon, The Quantum Theory of Light, 3rd ed., Ch. 2 · Allen & Eberly, Optical Resonance and Two-Level Atoms (the definitive coherent-transient reference) · Saleh & Teich, Fundamentals of Photonics, 3rd ed., Ch. 13–14 · Demtröder, Laser Spectroscopy, Vol. 1 · Steck, Quantum and Atom Optics (open lecture notes) for the optical Bloch equations. The key results to keep in view are: the RWA amplitude equations and Pe=sin⁡2(Ωt/2)P_e=\sin^2(\Omega t/2) at δ=0\delta=0; the general Rabi formula (5) with ΩR=Ω2+δ2\Omega_R=\sqrt{\Omega^2+\delta^2} and Pe,max⁡=Ω2/ΩR2P_{e,\max}=\Omega^2/\Omega_R^2; the steady state Pess=12 s/(1+s)P_e^{ss}=\tfrac12\,s/(1+s), s=(Ω2/2)/(δ2+Γ2/4)s=(\Omega^2/2)/(\delta^2+\Gamma^2/4), with power-broadened width Γ2+2Ω2\sqrt{\Gamma^2+2\Omega^2}; the coupling Ω=μE0/ℏ∝I\Omega=\mu E_0/\hbar\propto\sqrt I; and the anchor numbers E0=274E_0=274 V/m, Ω/2π=10.4\Omega/2\pi=10.4 MHz, τπ=48\tau_\pi=48 ns for μ=2.5×10−29\mu=2.5\times10^{-29} C·m at I=100I=100 W/m².

Prerequisites. Electric Dipole Approximation · Transition Dipole Moment · Stimulated Emission. Leads to. Radiative Lifetime and Natural Linewidth · Lorentzian Line Shape · Resonant Absorption.

Ask the tutor: Want to see the Rabi flop trace turn into a smooth absorption line as you dial up the decoherence rate Γ? Ask me to walk through the optical Bloch equations step by step.

Why does the excited population oscillate at Ω but the amplitude at Ω/2?

Because probability is amplitude squared. On resonance the amplitudes evolve as cos⁡(Ωt/2)\cos(\Omega t/2) and sin⁡(Ωt/2)\sin(\Omega t/2), each cycling at half the Rabi frequency. But the observable is Pe=∣ce∣2=sin⁡2(Ωt/2)P_e=|c_e|^2=\sin^2(\Omega t/2), and sin⁡2x=12(1−cos⁡2x)\sin^2 x=\tfrac12(1-\cos2x) doubles the frequency, so PeP_e oscillates at the full Ω\Omega. It is the same factor of two that makes a spin-1/2 need a 4π4\pi rotation to return exactly to itself while its population repeats every 2π2\pi of Ωt\Omega t. Keep the two clocks separate: amplitudes tick at Ω/2\Omega/2, populations at Ω\Omega.

Is energy conserved during a Rabi oscillation?

Yes, completely. While PeP_e climbs the atom absorbs energy; while it falls the atom returns that energy by stimulated emission, coherently and in phase with the drive, so the net exchange over a cycle is zero. This is exactly what makes Rabi cycling reversible and unlike absorption. Net one-way transfer appears only when decoherence (spontaneous emission, collisions) interrupts the descent so the atom cannot give back what it borrowed; then the exchange is irreversible and you measure genuine absorption. Reversible cycling and irreversible absorption are the same interaction with and without a working phase memory.

When is the rotating-wave approximation safe to use?

Whenever the Rabi frequency and detuning are both far below the transition frequency, Ω,∣δ∣≪ω0\Omega,|\delta|\ll\omega_0. For an optical transition ω0∼5×1014\omega_0\sim5\times10^{14} Hz while a strong Ω\Omega is ∼107\sim10^7–101010^{10} Hz, so the RWA is excellent — the discarded term oscillates at ≈2ω0\approx2\omega_0 and averages away. It matters only under extreme driving: the neglected term gives the Bloch–Siegert shift ∼Ω2/4ω0\sim\Omega^2/4\omega_0, seen in strongly driven NMR and superconducting circuits, and in the ultrastrong-coupling regime where Ω∼ω0\Omega\sim\omega_0 the simple formula fails and the full quantum Rabi model is needed. For essentially all of atomic and optical physics it is nothing to worry about.

What is the difference between the Rabi frequency Ω and the generalised Rabi frequency Ω_R?

Ω=μE0/ℏ\Omega=\mu E_0/\hbar is the on-resonance coupling — the cycling frequency when the laser sits exactly on the transition, δ=0\delta=0. The generalised Rabi frequency ΩR=Ω2+δ2\Omega_R=\sqrt{\Omega^2+\delta^2} is the actual oscillation frequency once you detune; because detuning adds in quadrature, ΩR≥Ω\Omega_R\ge\Omega always, and the two coincide only on resonance. The distinction matters twice: ΩR\Omega_R sets how fast the population oscillates, while the ratio Ω2/ΩR2\Omega^2/\Omega_R^2 sets how high it reaches. A detuned atom cycles faster yet never fully inverts — both facts encoded in the effective-field vector (Ω,0,−δ)(\Omega,0,-\delta) of length ΩR\Omega_R.

Why can't the steady-state excited population exceed one-half?

Once coherence is destroyed the atom obeys rate equations in which the field drives absorption (g→eg\to e) and stimulated emission (e→ge\to g) at equal rates. At Pe=1/2P_e=1/2 these balance — as many atoms driven down as up — so the population stops climbing however hard you drive, and spontaneous emission only pulls it back further. Hence Pess=12 s/(1+s)→1/2P_e^{ss}=\tfrac12\,s/(1+s)\to1/2. A coherent π\pi pulse can momentarily beat this ceiling and reach Pe=1P_e=1, but only within the decoherence time; once phase memory is lost the half-excitation cap reasserts itself. Full inversion is a coherent privilege, not a steady-state option.

How is a Rabi oscillation actually observed in the lab?

You prepare the atom (or ion, or qubit) in the ground state, drive it resonantly for a controlled time tt, switch off, and read out the excited population — often by fluorescence, since ∣e⟩|e\rangle scatters photons. Plot PeP_e versus tt and you see the sinusoidal "Rabi flop" of period 2π/Ω2\pi/\Omega, giving Ω\Omega directly; scanning the amplitude confirms Ω∝I\Omega\propto\sqrt I and scanning the detuning maps out ΩR\Omega_R. In practice the flops decay: the envelope's rate measures decoherence, and the number of clean cycles before it dies is the coherence figure of merit qubit engineers work to maximise.

Before this

Where this leads