Larmor Precession and Magnetic Resonance

Put a spin-12\tfrac12 particle in a static magnetic field and nothing dramatic seems to happen: the two energy levels split by the Zeeman gap, and a spin already aligned or anti-aligned with the field just sits there. The interesting physics lives in the superposition. A spin tilted away from the field axis does not stay put — its expectation value ⟨S⟩\langle\mathbf{S}\rangle sweeps around the field direction at a single, field-set rate called the Larmor frequency ω0=γB0\omega_0=\gamma B_0, tracing a cone exactly as a classical gyroscope precesses under gravity. This is Larmor precession, and the remarkable feature is that the quantum expectation values obey the same equation of motion as a classical magnetic moment, a clean instance of Ehrenfest's theorem. The precession by itself moves nothing between the two energy levels; ⟨Sz⟩\langle S_z\rangle is a constant of the motion. To drive real transitions you add a weak oscillating field transverse to B0B_0 and tune its frequency near ω0\omega_0. On resonance the spin flips completely back and forth between up and down — Rabi oscillations — and off resonance the flip amplitude is sharply suppressed. The resonance condition ω=ω0\omega=\omega_0 and the flip probability are captured by the Rabi formula, derived here in the rotating-wave approximation. This single two-level problem is the engine behind an enormous body of technology: nuclear magnetic resonance and MRI (protons at 42.57742.577 MHz per tesla), electron spin resonance (28.028.0 GHz per tesla), atomic clocks, and the control pulses that rotate a qubit on the Bloch sphere. This article derives the precession, the resonance, and the numbers, and shows where each claim is proved and where it is merely quoted.

Key takeaways
  • A spin-12\tfrac12 magnetic moment in a static field B=B0z^\mathbf{B}=B_0\hat{\mathbf z} has H=−μ⋅B=ω0SzH=-\boldsymbol{\mu}\cdot\mathbf{B}=\omega_0 S_z, splitting the levels by ℏω0\hbar\omega_0 with Larmor frequency ω0=γB0\omega_0=\gamma B_0.
  • A tilted spin precesses: ⟨Sx⟩,⟨Sy⟩\langle S_x\rangle,\langle S_y\rangle rotate about z^\hat{\mathbf z} at ω0\omega_0 while ⟨Sz⟩\langle S_z\rangle stays constant — the quantum ⟨S⟩\langle\mathbf{S}\rangle obeys the classical torque equation ⟨S⟩˙=γ ⟨S⟩×B\dot{\langle\mathbf S\rangle}=\gamma\,\langle\mathbf S\rangle\times\mathbf B (Ehrenfest).
  • Precession alone drives no transitions; ⟨Sz⟩\langle S_z\rangle is conserved. Level-to-level flips require a transverse oscillating field.
  • Adding a rotating field of frequency ω\omega gives Rabi oscillations: P↑→↓(t)=ω12Ω2sin⁡2 ⁣(12Ωt)P_{\uparrow\to\downarrow}(t)=\dfrac{\omega_1^2}{\Omega^2}\sin^2\!\big(\tfrac12\Omega t\big) with generalized Rabi frequency Ω=ω12+Δ2\Omega=\sqrt{\omega_1^2+\Delta^2} and detuning Δ=ω−ω0\Delta=\omega-\omega_0.
  • On resonance (Δ=0\Delta=0) the spin flips completely, P=sin⁡2(12ω1t)P=\sin^2(\tfrac12\omega_1 t), giving π\pi-pulses and π/2\pi/2-pulses; off resonance the flip amplitude is capped at ω12/Ω2\omega_1^2/\Omega^2.
  • Anchors: proton NMR 42.57742.577 MHz/T (so 63.963.9 MHz in a 1.51.5 T MRI), electron ESR 28.028.0 GHz/T (9.89.8 GHz in 0.350.35 T).
  • Key idea: a magnetic field makes the transverse spin components rotate at ω0=γB0\omega_0=\gamma B_0 but cannot change ⟨Sz⟩\langle S_z\rangle; only a field oscillating near ω0\omega_0 resonantly pumps population between the levels, and that resonance is the whole basis of NMR, ESR, and qubit control.
The core insight

A static field sets a precession rate ω0=γB0\omega_0=\gamma B_0 for the transverse spin but conserves ⟨Sz⟩\langle S_z\rangle; a transverse field tuned to ω0\omega_0 resonantly drives full population transfer at the Rabi frequency. Precession and resonance are the two faces of the same two-level Hamiltonian.

Where this lives today

The two-level physics of this article is, without exaggeration, one of the most technologically consequential calculations in all of quantum mechanics. Nuclear magnetic resonance reads the Larmor frequencies of nuclei (protons, 13^{13}C, 15^{15}N) shifted by their chemical environment, making it the workhorse of molecular structure determination in chemistry and structural biology. Magnetic resonance imaging turns the same proton signal into spatial images by adding field gradients so that Larmor frequency encodes position — no ionizing radiation, unlike X-rays. Electron spin resonance probes unpaired electrons in radicals and transition-metal centers. Atomic clocks lock an oscillator to a hyperfine resonance measured by exactly this Rabi technique, defining the SI second. And in quantum computing, every single-qubit gate is a Rabi rotation: a resonant microwave π\pi- or π/2\pi/2-pulse rotates a superconducting or spin qubit on the Bloch sphere, with the pulse area ω1t\omega_1 t setting the rotation angle. The Bloch-sphere language physicists use for qubits is precisely the geometry of ⟨S⟩\langle\mathbf S\rangle derived here.

The physical picture

Think of the spin's magnetic moment μ=γS\boldsymbol{\mu}=\gamma\mathbf{S} as a tiny bar magnet carrying angular momentum. A static field exerts a torque τ=μ×B\boldsymbol{\tau}=\boldsymbol{\mu}\times\mathbf{B}, and because the moment is rigidly tied to angular momentum, the torque does not tip the magnet toward the field the way it would tip a compass needle with friction — instead it makes the angular momentum vector precess around the field, sweeping out a cone at fixed polar angle. This is identical to a spinning top under gravity. The precession rate is the Larmor frequency ω0=γB0\omega_0=\gamma B_0, independent of the tilt angle: a spin tipped 30∘30^\circ and a spin tipped 80∘80^\circ circle the field axis at exactly the same rate. The quantum content is that ⟨S⟩\langle\mathbf{S}\rangle obeys precisely this classical equation. Ehrenfest's theorem turns the Heisenberg equation for S\mathbf{S} into ⟨S⟩˙=γ ⟨S⟩×B\dot{\langle\mathbf S\rangle}=\gamma\,\langle\mathbf S\rangle\times\mathbf B, so the expectation value traces the cone even though the spin itself is in a superposition of ∣↑⟩|{\uparrow}\rangle and ∣↓⟩|{\downarrow}\rangle with fixed populations. Nothing moves between the levels: ⟨Sz⟩\langle S_z\rangle is a constant.

Figure
A tilted spin precesses about the field at the Larmor frequency
Larmor precession of ⟨S⟩ about B₀ B₀ ẑ ⟨S⟩ ω₀ = γB₀ α ⟨Sₓ⟩ = ½ sinα cos(ω₀t) ⟨Sᵧ⟩ = ½ sinα sin(ω₀t) ⟨S_z⟩ = ½ cosα (constant) Tilt angle α is frozen; only the azimuth advances at ω₀.

To actually flip a spin you must break the axial symmetry. Add a small field rotating in the xyxy-plane at frequency ω\omega. Viewed from a frame co-rotating with this field, the fast Larmor rotation is largely undone, and what remains is an effective static field tilted away from z^\hat{\mathbf z} by an amount set by the detuning Δ=ω−ω0\Delta=\omega-\omega_0. When ω\omega exactly matches ω0\omega_0 the effective field lies flat in the transverse plane, and the spin precesses about it straight from up to down and back — a complete flip. Off resonance the effective field has a residual zz-component, the precession cone about it never reaches the south pole, and the flip is only partial. This is the whole content of magnetic resonance.

Figure
Rabi flip probability: sharp on resonance, capped off resonance
Spin-flip probability P(t) at three detunings 1 0 time → Δ = 0 (resonance): full flips to 1 Δ = ω₁: capped at ½ Δ = 3ω₁: barely flips Amplitude cap = ω₁² / (ω₁² + Δ²); oscillation rate = Ω = √(ω₁² + Δ²).

The core equations

The two-level Hamiltonian for a magnetic moment μ=γS\boldsymbol{\mu}=\gamma\mathbf{S} in a static field along z^\hat{\mathbf z} is a Zeeman term. Writing it with the sign absorbed into the rate keeps the algebra clean.

H0=−μ⋅B=−γB0 Sz=ω0 Sz,ω0≡γB0(1)H_0 = -\boldsymbol{\mu}\cdot\mathbf{B} = -\gamma B_0\,S_z = \omega_0\,S_z, \qquad \omega_0 \equiv \gamma B_0 \qquad (1)

where: γ\gamma is the gyromagnetic ratio (γ=g q/2m\gamma=g\,q/2m for a particle of charge qq, mass mm, and gg-factor gg); B0B_0 is the static field magnitude; Sz=ℏ2σzS_z=\tfrac\hbar2\sigma_z is the spin-zz operator; and ω0=γB0\omega_0=\gamma B_0 is the Larmor frequency. We have folded the sign of γ\gamma into ω0\omega_0 so that H0=ω0SzH_0=\omega_0 S_z; the eigenvalues are ±12ℏω0\pm\tfrac12\hbar\omega_0, split by the Zeeman gap ℏω0\hbar\omega_0. (For the electron γ<0\gamma<0, which flips the precession sense but not the rate.) Time evolution of an arbitrary initial spinor follows from the propagator U(t)=e−iH0t/ℏU(t)=e^{-iH_0t/\hbar}. For a state tilted from z^\hat{\mathbf z} by polar angle α\alpha, the expectation values are

⟨Sx⟩=ℏ2sin⁡α cos⁡ω0t,⟨Sy⟩=ℏ2sin⁡α sin⁡ω0t,⟨Sz⟩=ℏ2cos⁡α(2)\langle S_x\rangle = \tfrac\hbar2\sin\alpha\,\cos\omega_0 t,\quad \langle S_y\rangle = \tfrac\hbar2\sin\alpha\,\sin\omega_0 t,\quad \langle S_z\rangle = \tfrac\hbar2\cos\alpha \qquad (2)

where: α\alpha is the (constant) tilt angle of the initial spin. The transverse components rotate at ω0\omega_0 with fixed magnitude ℏ2sin⁡α\tfrac\hbar2\sin\alpha, tracing a circle, while ⟨Sz⟩\langle S_z\rangle is time-independent. This is precession: the polar angle is frozen, only the azimuth advances. The sense of rotation (+sin⁡ω0t+\sin\omega_0 t for ⟨Sy⟩\langle S_y\rangle with our H0=+ω0SzH_0=+\omega_0 S_z) is fixed by the sign of γ\gamma. The same content in vector form is the classical torque equation, which Ehrenfest's theorem hands us exactly for the quantum expectation value.

d⟨S⟩dt=γ ⟨S⟩×B(3)\frac{d\langle\mathbf{S}\rangle}{dt} = \gamma\,\langle\mathbf{S}\rangle\times\mathbf{B} \qquad (3)

where: B\mathbf{B} is the total magnetic field. Equation (3) is a closed equation for ⟨S⟩\langle\mathbf S\rangle (no higher moments appear because HH is linear in S\mathbf S), so the quantum precession is literally the classical one; a solid magnetization M=nγℏ⟨σ⟩/2\mathbf M=n\gamma\hbar\langle\boldsymbol\sigma\rangle/2 obeys the identical equation, which is why bulk NMR magnetization precesses too. To drive transitions, add a field rotating in the xyxy-plane at drive frequency ω\omega. The full Hamiltonian, transformed into the frame co-rotating with the drive and dropping the fast counter-rotating term (the rotating-wave approximation, RWA), becomes time-independent:

Hrot=ℏ2(−Δ σz+ω1 σx),Δ≡ω−ω0,ω1≡γB1(4)H_{\rm rot} = \tfrac\hbar2\big(-\Delta\,\sigma_z + \omega_1\,\sigma_x\big),\qquad \Delta \equiv \omega-\omega_0,\quad \omega_1 \equiv \gamma B_1 \qquad (4)

where: Δ\Delta is the detuning of the drive from resonance; ω1=γB1\omega_1=\gamma B_1 is the Rabi frequency set by the rotating field amplitude B1B_1; and σx,z\sigma_{x,z} are Pauli matrices. In the rotating frame the spin sees a static effective field with a zz-component ∝Δ\propto\Delta and a transverse component ∝ω1\propto\omega_1; it precesses about that tilted effective field. Solving (4) for a spin starting in ∣↑⟩|{\uparrow}\rangle gives the probability of finding it flipped to ∣↓⟩|{\downarrow}\rangle — the Rabi formula.

P↑→↓(t)=ω12ω12+Δ2 sin⁡2 ⁣(12ω12+Δ2  t)=ω12Ω2sin⁡2 ⁣(12Ωt)(5)P_{\uparrow\to\downarrow}(t) = \frac{\omega_1^2}{\omega_1^2+\Delta^2}\,\sin^2\!\Big(\tfrac12\sqrt{\omega_1^2+\Delta^2}\;t\Big) = \frac{\omega_1^2}{\Omega^2}\sin^2\!\big(\tfrac12\Omega t\big) \qquad (5)

where: Ω≡ω12+Δ2\Omega\equiv\sqrt{\omega_1^2+\Delta^2} is the generalized Rabi frequency. The prefactor ω12/Ω2\omega_1^2/\Omega^2 caps the maximum flip probability, and Ω\Omega sets the oscillation rate. On resonance (Δ=0\Delta=0): P=sin⁡2(12ω1t)P=\sin^2(\tfrac12\omega_1 t), reaching 11 at t=π/ω1t=\pi/\omega_1 (a π\pi-pulse). The width of the resonance in Δ\Delta is ∼ω1\sim\omega_1, so a weaker drive gives a sharper, more selective resonance.

Derivation

Derivation ALarmor precession from the Schrödinger equation

The goal is to show that a spin-12\tfrac12 state tilted away from the field axis has transverse spin components that rotate at ω0=γB0\omega_0=\gamma B_0 while its axial component stays fixed, and that this reproduces the classical torque equation exactly. The route is to write the initial tilted spinor in the SzS_z eigenbasis, apply the diagonal propagator e−iH0t/ℏe^{-iH_0t/\hbar} to attach a relative phase between the up and down amplitudes, and then compute the three spin expectation values, watching how that single relative phase turns into a rotating transverse vector. The one thing to keep straight is that the propagator only changes the relative phase of the two amplitudes, never their magnitudes — which is exactly why ⟨Sz⟩\langle S_z\rangle cannot change.

Step 1. Write the initial state tilted by polar angle α\alpha (azimuth zero): ∣χ(0)⟩=cos⁡α2 ∣↑⟩+sin⁡α2 ∣↓⟩|\chi(0)\rangle=\cos\tfrac\alpha2\,|{\uparrow}\rangle+\sin\tfrac\alpha2\,|{\downarrow}\rangle. This has ⟨Sz⟩=ℏ2(cos⁡2α2−sin⁡2α2)=ℏ2cos⁡α\langle S_z\rangle=\tfrac\hbar2(\cos^2\tfrac\alpha2-\sin^2\tfrac\alpha2)=\tfrac\hbar2\cos\alpha, as required for a spin pointing at polar angle α\alpha.

Step 2. Apply the propagator. Since H0=ω0Sz=ℏ2ω0σzH_0=\omega_0 S_z=\tfrac\hbar2\omega_0\sigma_z is diagonal, e−iH0t/ℏe^{-iH_0t/\hbar} multiplies ∣↑⟩|{\uparrow}\rangle by e−iω0t/2e^{-i\omega_0 t/2} and ∣↓⟩|{\downarrow}\rangle by e+iω0t/2e^{+i\omega_0 t/2}: ∣χ(t)⟩=cos⁡α2 e−iω0t/2∣↑⟩+sin⁡α2 e+iω0t/2∣↓⟩|\chi(t)\rangle=\cos\tfrac\alpha2\,e^{-i\omega_0 t/2}|{\uparrow}\rangle+\sin\tfrac\alpha2\,e^{+i\omega_0 t/2}|{\downarrow}\rangle. The magnitudes of the two amplitudes are untouched.

Step 3. Compute ⟨Sz⟩\langle S_z\rangle. Because the propagator left ∣cos⁡α2∣2|\cos\tfrac\alpha2|^2 and ∣sin⁡α2∣2|\sin\tfrac\alpha2|^2 unchanged, ⟨Sz(t)⟩=ℏ2cos⁡α\langle S_z(t)\rangle=\tfrac\hbar2\cos\alpha for all tt. The axial component is a constant of the motion — no population moves between levels.

Step 4. Compute the transverse components. Using Sx=ℏ2(∣↑⟩⟨↓∣+∣↓⟩⟨↑∣)S_x=\tfrac\hbar2(|{\uparrow}\rangle\langle{\downarrow}|+|{\downarrow}\rangle\langle{\uparrow}|), the cross terms carry the relative phase e±iω0te^{\pm i\omega_0 t}: ⟨Sx(t)⟩=ℏ2⋅2cos⁡α2sin⁡α2cos⁡ω0t=ℏ2sin⁡αcos⁡ω0t\langle S_x(t)\rangle=\tfrac\hbar2\cdot 2\cos\tfrac\alpha2\sin\tfrac\alpha2\cos\omega_0 t=\tfrac\hbar2\sin\alpha\cos\omega_0 t, and similarly ⟨Sy(t)⟩=ℏ2sin⁡αsin⁡ω0t\langle S_y(t)\rangle=\tfrac\hbar2\sin\alpha\sin\omega_0 t. The transverse vector has fixed length ℏ2sin⁡α\tfrac\hbar2\sin\alpha and rotates at ω0\omega_0.

Step 5. Identify the classical equation. Differentiating (2) gives ⟨Sx⟩˙=−ω0⟨Sy⟩\dot{\langle S_x\rangle}=-\omega_0\langle S_y\rangle, ⟨Sy⟩˙=+ω0⟨Sx⟩\dot{\langle S_y\rangle}=+\omega_0\langle S_x\rangle, ⟨Sz⟩˙=0\dot{\langle S_z\rangle}=0, which is exactly ⟨S⟩˙=γ⟨S⟩×B\dot{\langle\mathbf S\rangle}=\gamma\langle\mathbf S\rangle\times\mathbf B with B=B0z^\mathbf B=B_0\hat{\mathbf z} and ω0=γB0\omega_0=\gamma B_0. The quantum expectation value precesses identically to a classical moment. Verified: the companion verifier integrates the spin-12\tfrac12 TDSE for a tilted spinor (α=0.7\alpha=0.7) and confirms ⟨Sz⟩=12cos⁡α\langle S_z\rangle=\tfrac12\cos\alpha to machine precision, ⟨Sx⟩=12sin⁡αcos⁡ω0t\langle S_x\rangle=\tfrac12\sin\alpha\cos\omega_0 t and ⟨Sy⟩=+12sin⁡αsin⁡ω0t\langle S_y\rangle=+\tfrac12\sin\alpha\sin\omega_0 t exactly, and an FFT of ⟨Sx(t)⟩\langle S_x(t)\rangle recovers ω0\omega_0 to better than 1%1\%.

Go deeper · Derivation B · 2 minThe Rabi formula in the rotating-wave approximation
Derivation BThe Rabi formula in the rotating-wave approximation

The aim here is to derive the two-level flip probability (5) — the quantitative statement of magnetic resonance — from the Hamiltonian of a spin in a static field plus a rotating transverse field. The route is to remove the fast Larmor rotation by transforming into a frame co-rotating with the drive, which converts the explicitly time-dependent Hamiltonian into a static effective Hamiltonian describing precession about a tilted effective field; then the flip probability is just the precession of a spin about that fixed axis, read off geometrically or by diagonalizing the 2×22\times2 matrix. The subtlety to watch is the rotating-wave approximation: a linearly oscillating field is a sum of two counter-rotating components, and we keep only the one co-rotating with the spin, discarding a term that oscillates at 2ω2\omega and averages away when ω1≪ω0\omega_1\ll\omega_0.

Step 1. Total field and Hamiltonian: B(t)=B0z^+B1(cos⁡ωt x^−sin⁡ωt y^)\mathbf B(t)=B_0\hat{\mathbf z}+B_1(\cos\omega t\,\hat{\mathbf x}-\sin\omega t\,\hat{\mathbf y}), so H(t)=ω0Sz+ω1(cos⁡ωt Sx−sin⁡ωt Sy)H(t)=\omega_0 S_z+\omega_1(\cos\omega t\,S_x-\sin\omega t\,S_y) with ω1=γB1\omega_1=\gamma B_1. The transverse field rotates in the plane at frequency ω\omega.

Step 2. Transform to the rotating frame via R(t)=e−iωt Sz/ℏR(t)=e^{-i\omega t\,S_z/\hbar}, defining ∣χ~⟩=R†∣χ⟩|\tilde\chi\rangle=R^\dagger|\chi\rangle. The transformed Hamiltonian is H~=R†HR−iℏR†R˙\tilde H=R^\dagger H R-i\hbar R^\dagger\dot R. The rotation cancels the explicit time dependence of the transverse field, and the −iℏR†R˙=−ωSz-i\hbar R^\dagger\dot R=-\omega S_z term shifts the axial part.

Step 3. Collect terms: H~=(ω0−ω)Sz+ω1Sx=−Δ Sz+ω1Sx\tilde H=(\omega_0-\omega)S_z+\omega_1 S_x=-\Delta\,S_z+\omega_1 S_x with Δ=ω−ω0\Delta=\omega-\omega_0, i.e. H~=ℏ2(−Δ σz+ω1σx)\tilde H=\tfrac\hbar2(-\Delta\,\sigma_z+\omega_1\sigma_x) — equation (4). This is time-independent: a static effective field Beff∝(ω1,0,−Δ)\mathbf B_{\rm eff}\propto(\omega_1,0,-\Delta) in the rotating frame.

Step 4. The spin precesses about Beff\mathbf B_{\rm eff} at the generalized Rabi frequency Ω=ω12+Δ2\Omega=\sqrt{\omega_1^2+\Delta^2} (the magnitude of the effective-field vector). Starting from ∣↑⟩|{\uparrow}\rangle, the probability to be found in ∣↓⟩|{\downarrow}\rangle is the squared projection of the precessing state onto the south pole, which geometry (or diagonalizing the 2×22\times2) gives as P=ω12Ω2sin⁡2(12Ωt)P=\dfrac{\omega_1^2}{\Omega^2}\sin^2(\tfrac12\Omega t).

Step 5. Read off the limits. On resonance Δ=0\Delta=0: Ω=ω1\Omega=\omega_1 and P=sin⁡2(12ω1t)P=\sin^2(\tfrac12\omega_1 t), a complete flip at t=π/ω1t=\pi/\omega_1. Far off resonance ∣Δ∣≫ω1|\Delta|\gg\omega_1: the prefactor ω12/Ω2→0\omega_1^2/\Omega^2\to0, so the spin barely flips — the resonance is a genuine peak in Δ\Delta of width ∼ω1\sim\omega_1. Verified: the verifier evolves the rotating-frame Hamiltonian (4) directly for detunings Δ∈{0,0.5,1.7,−1.3}ω1\Delta\in\{0,0.5,1.7,-1.3\}\omega_1 and confirms the analytic Rabi formula (5) matches the numerical ∣⟨↓∣ψ(t)⟩∣2|\langle{\downarrow}|\psi(t)\rangle|^2 to ∼10−15\sim10^{-15}, that on resonance P(π/ω1)=1P(\pi/\omega_1)=1, and that the far-detuned amplitude caps at ω12/Ω2\omega_1^2/\Omega^2.

Limiting cases

Each regime of the two-level problem is a check on the formulas.

CaseResultMeaning
Spin aligned with field (α=0\alpha=0)⟨S⟩\langle\mathbf S\rangle static, no precessionAn eigenstate of H0H_0; only a stationary phase accrues, nothing observable moves.
On resonance (Δ=0\Delta=0)P=sin⁡2(12ω1t)P=\sin^2(\tfrac12\omega_1 t), full flipsEffective field lies in the transverse plane; the spin precesses straight from up to down.
Far off resonance (∣Δ∣≫ω1|\Delta|\gg\omega_1)Pmax⁡=ω12/Δ2→0P_{\max}=\omega_1^2/\Delta^2\to0Effective field nearly along z^\hat{\mathbf z}; the precession cone never reaches the south pole.
Classical / large spinM˙=γM×B\dot{\mathbf M}=\gamma\mathbf M\times\mathbf B (Bloch)⟨S⟩\langle\mathbf S\rangle obeys the classical torque equation exactly (Ehrenfest); bulk magnetization precesses identically.
Strong static field (ω0→∞\omega_0\to\infty)RWA becomes exactThe counter-rotating term oscillates ever faster and averages away, so the rotating-wave approximation improves.

Worked examples

Example 1 — NMR and ESR frequencies. Find the proton NMR frequency in a 1.51.5 T clinical MRI magnet and the electron ESR frequency in a 0.350.35 T laboratory magnet.

Step 1. Proton NMR. The proton gyromagnetic ratio is γp/2π=42.577\gamma_p/2\pi=42.577 MHz/T (this already includes the proton gg-factor gp≈5.586g_p\approx5.586). The Larmor frequency is f=γpB0/2π=(42.577 MHz/T)×(1.5 T)f=\gamma_p B_0/2\pi=(42.577\ \text{MHz/T})\times(1.5\ \text{T}).

Step 2. Evaluate: f=63.87 MHz≈63.9 MHzf=63.87\ \text{MHz}\approx\mathbf{63.9\ MHz}. This is the radio-frequency at which a 1.51.5 T MRI scanner transmits and receives — squarely in the FM-radio band, which is why MRI hardware uses RF coils.

Step 3. Electron ESR. The electron's effective ratio is ∣γe∣/2π=geμB/h=28.025|\gamma_e|/2\pi=g_e\mu_B/h=28.025 GHz/T (with ge=2.0023g_e=2.0023). In 0.350.35 T: f=(28.025 GHz/T)×(0.35 T)=9.81 GHzf=(28.025\ \text{GHz/T})\times(0.35\ \text{T})=\mathbf{9.81\ GHz} — the microwave X-band, and the reason ESR spectrometers use waveguides and cavities rather than RF coils. The ratio of the two, ≈28.0 GHz/42.6 MHz≈660\approx28.0\ \text{GHz}/42.6\ \text{MHz}\approx660, is essentially the proton-to-electron effective moment ratio — the electron's moment is far larger, so it resonates at far higher frequency in the same field. Takeaways: 63.963.9 MHz proton NMR at 1.51.5 T, 9.89.8 GHz electron ESR at 0.350.35 T.

Example 2 — Designing a π/2\pi/2-pulse. An NMR spectrometer applies an on-resonance rotating field giving a Rabi frequency ω1/2π=25\omega_1/2\pi=25 kHz. How long a pulse rotates the magnetization from z^\hat{\mathbf z} into the transverse plane (a π/2\pi/2-pulse), and how long fully inverts it (a π\pi-pulse)?

Step 1. On resonance the polar angle swept is θ(t)=ω1t\theta(t)=\omega_1 t. A π/2\pi/2-pulse needs θ=π/2\theta=\pi/2, so tπ/2=π/2ω1=14⋅2πω1=14⋅125 kHz=10 μst_{\pi/2}=\dfrac{\pi/2}{\omega_1}=\dfrac{1}{4}\cdot\dfrac{2\pi}{\omega_1}=\dfrac{1}{4}\cdot\dfrac{1}{25\,\text{kHz}}=10\ \mu\text{s}.

Step 2. A π\pi-pulse (full inversion, P=1P=1) needs θ=π\theta=\pi: tπ=π/ω1=2 tπ/2=20 μst_\pi=\pi/\omega_1=2\,t_{\pi/2}=20\ \mu\text{s}.

Step 3. Sanity check with the Rabi formula on resonance: P(tπ)=sin⁡2(12ω1tπ)=sin⁡2(π/2)=1P(t_\pi)=\sin^2(\tfrac12\omega_1 t_\pi)=\sin^2(\pi/2)=1, a complete flip; P(tπ/2)=sin⁡2(π/4)=12P(t_{\pi/2})=\sin^2(\pi/4)=\tfrac12, an equal superposition with the moment in the transverse plane. A 2525 kHz Rabi drive gives a 10 μ10\ \mus π/2\pi/2-pulse and a 20 μ20\ \mus π\pi-pulse.

Origins · 1 minlarmor precession and how this was discovered
LP
larmor precession
classical + quantum torque on a magnetic moment · The rotation of a magnetic moment about an applied field

Larmor precession is named for Joseph Larmor, who in 1897 derived that an orbiting charged particle's angular momentum precesses about an applied magnetic field at the frequency ω0=γB\omega_0=\gamma B, one of the earliest results linking magnetism to rotational dynamics. Quantum-mechanically the concept survives intact: the spin expectation value ⟨S⟩\langle\mathbf S\rangle precesses at exactly the Larmor frequency, a consequence of Ehrenfest's theorem that the two-level Hamiltonian is linear in the spin. The concept is the bridge between the classical picture of a spinning magnetized top and the quantum two-level system, and it sets the carrier frequency of every magnetic-resonance technique from NMR to qubit control.

History
1897
— Larmor derives the precession of an orbital magnetic moment about a field, ω=γB\omega=\gamma B
1922
— Stern and Gerlach demonstrate space quantization of the silver-atom magnetic moment
1938
— Rabi and coworkers observe molecular-beam magnetic resonance (Phys. Rev. 53, 318), the first resonant spin flips
1944
— Rabi awarded the Nobel Prize in Physics for the resonance method
1946
— Bloch (Stanford) and Purcell (Harvard) independently detect nuclear magnetic resonance in bulk matter; Bloch writes the phenomenological Bloch equations
1952
— Bloch and Purcell share the Nobel Prize in Physics for nuclear magnetic resonance
1973
— Lauterbur and Mansfield develop magnetic-resonance imaging from the NMR signal

Going further

References. Griffiths & Schroeter, Introduction to Quantum Mechanics, 3rd ed., §4.4 (spin in a magnetic field) · Sakurai & Napolitano, Modern Quantum Mechanics, 2nd ed., §2.1 (spin precession) and §5.5 (magnetic resonance) · Cohen-Tannoudji, Diu & Laloë, Quantum Mechanics, Ch. IV complement F · Slichter, Principles of Magnetic Resonance.

Provenance. A verifier should reproduce: (1) the precession law (2) by integrating the spin-12\tfrac12 TDSE for a tilted spinor — ⟨Sz⟩=12cos⁡α\langle S_z\rangle=\tfrac12\cos\alpha constant, transverse components rotating at ω0\omega_0; (2) the Rabi formula (5) against direct rotating-frame evolution of (4) for several detunings, matching to ∼10−15\sim10^{-15}; (3) the anchors proton 42.57742.577 MHz/T, 63.963.9 MHz at 1.51.5 T, electron 28.028.0 GHz/T, 9.89.8 GHz at 0.350.35 T, and the π\pi-pulse time tπ=π/ω1t_\pi=\pi/\omega_1.

Builds on: Spin and Spin Angular Momentum · The Stern-Gerlach Experiment · Time-Dependent Schrödinger Equation · Ehrenfest's Theorem. Leads to: Time-Dependent Perturbation Theory and Rabi Oscillations · Nuclear Magnetic Resonance and MRI · Atomic Clocks · Single-Qubit Gates and the Bloch Sphere.

Why does the spin precess instead of just aligning with the field like a compass needle?

A compass needle aligns because friction dissipates its energy until it settles at the minimum-energy orientation. An isolated spin has no such dissipation, and crucially its magnetic moment is rigidly tied to its angular momentum. The field exerts a torque μ×B\boldsymbol{\mu}\times\mathbf B perpendicular to both μ\boldsymbol{\mu} and B\mathbf B, and since torque equals the rate of change of angular momentum, the angular momentum vector rotates around the field rather than toward it — exactly like a spinning top under gravity precesses instead of falling over. Only when you add a relaxation mechanism (spin-lattice coupling, the T1T_1 process in NMR) does the spin actually relax toward alignment, and that happens on a much slower timescale than the precession.

Is Larmor precession a genuinely quantum effect or just classical?

The precession equation (3) is identical to the classical torque equation, so in that sense it is not distinctively quantum — a classical magnetized gyroscope precesses the same way, and this is a clean example of the classical limit emerging via Ehrenfest's theorem because the Hamiltonian is linear in the spin. What is quantum is that the state precessing is a superposition of just two discrete levels, that the energy gap ℏω0\hbar\omega_0 is quantized, and that a measurement of SzS_z yields only ±ℏ2\pm\tfrac\hbar2. The expectation value follows the classical cone, but any single measurement collapses the spin. So the trajectory of ⟨S⟩\langle\mathbf S\rangle is classical; the measurement statistics are quantum.

What exactly is the rotating-wave approximation, and when does it fail?

A linearly oscillating field 2B1cos⁡ωt x^2B_1\cos\omega t\,\hat{\mathbf x} decomposes into two circularly polarized components rotating in opposite senses at ±ω\pm\omega. Only the component co-rotating with the spin's Larmor precession stays nearly in step and drives the transition; the counter-rotating component runs at ≈2ω\approx2\omega relative to the spin and averages to nearly zero. The RWA keeps the co-rotating term and drops the counter-rotating one, which is excellent when ω1≪ω0\omega_1\ll\omega_0 (weak drive, near resonance). It fails for very strong driving (ω1∼ω0\omega_1\sim\omega_0), where the dropped term produces a small frequency shift (the Bloch-Siegert shift) and additional structure. In NMR and most qubit control ω1/ω0\omega_1/\omega_0 is tiny, so the RWA is superb.

Does the sign of the gyromagnetic ratio matter?

It flips the sense of precession but not the rate. For nuclei like the proton γ>0\gamma>0 and ⟨S⟩\langle\mathbf S\rangle precesses one way; for the electron γ<0\gamma<0 (its moment is antiparallel to its spin), so it precesses the opposite way at the same ∣γ∣B0|\gamma|B_0. The Larmor frequency magnitude ω0=∣γ∣B0\omega_0=|\gamma|B_0 and all resonance conditions are unchanged; only the direction of circulation on the Bloch sphere reverses. In most experiments only the magnitude matters, which is why tables quote ∣γ∣/2π|\gamma|/2\pi.

Why is a weaker drive better for resolving two nearby resonances?

The Rabi resonance in (5) has a width in detuning of order ω1\omega_1: the flip probability is appreciable only for ∣Δ∣≲ω1|\Delta|\lesssim\omega_1. So a weaker rotating field (smaller ω1\omega_1) produces a narrower resonance, letting you distinguish two nuclei whose Larmor frequencies differ only slightly (as in chemical-shift NMR). The tradeoff is speed: a smaller ω1\omega_1 means a longer π\pi-pulse (tπ=π/ω1t_\pi=\pi/\omega_1) and a slower experiment. This is the frequency-resolution-versus-speed tradeoff at the heart of pulse design, and it is a time-energy uncertainty relation in disguise: a longer pulse has a narrower frequency bandwidth.

How does this connect to the Bloch sphere used in quantum computing?

The Bloch sphere is exactly the space of directions of ⟨S⟩\langle\mathbf S\rangle (equivalently ⟨σ⟩\langle\boldsymbol\sigma\rangle) for a two-level system. Free precession about B0B_0 is rotation about the sphere's polar (zz) axis at ω0\omega_0; a resonant drive is rotation about a transverse axis at ω1\omega_1. A π\pi-pulse is a 180∘180^\circ rotation (bit flip, the XX gate up to phase), a π/2\pi/2-pulse is a 90∘90^\circ rotation that creates superpositions (a Hadamard-like operation up to phase), and the rotation axis in the transverse plane is set by the drive's phase. Every single-qubit gate is one of these Rabi rotations, so the physics derived here is literally the instruction set of a spin or superconducting qubit.

If precession conserves energy, where does the NMR signal come from?

The detected NMR signal comes from the transverse magnetization precessing in the coil, which by Faraday induction generates an oscillating EMF at the Larmor frequency — this is the free-induction decay. To get transverse magnetization in the first place you apply a π/2\pi/2-pulse (a driven, resonant process that does exchange energy with the RF field), tipping the equilibrium axial magnetization into the transverse plane. After the pulse ends, the transverse magnetization precesses freely (conserving energy) and is what the coil picks up, decaying over the T2T_2 time as spins dephase. So energy is exchanged during the pulse; the readout is free precession afterward.

Before this