Atomic Physics

Lamb Shift

The 1057 MHz gap that should not exist — where the empty vacuum reached into a hydrogen atom and nudged an electron

A tiny 1057 MHz splitting between hydrogen's 2s and 2p levels that Dirac theory predicts degenerate, caused by QED vacuum fluctuations.

  • 2s–2p splitting≈ 1057 MHz (≈ 4.4 µeV)
  • Dirac prediction2s₁⁄₂ and 2p₁⁄₂ exactly degenerate
  • CauseQED vacuum fluctuations (self-energy)
  • MeasuredLamb & Retherford, 1947
  • Nobel PrizeWillis Lamb, 1955
  • LegacyFirst confirmation of renormalization

Interactive visualization

Press play, or step through manually. The visualization is yours to drive — try it before reading on.

Open visualization fullscreen ↗

Watch the 60-second explainer

A condensed visual walkthrough — narrated, captioned, under a minute.

Definition

The Lamb shift is a small energy splitting between two hydrogen levels that the best pre-1947 theory said should be identical. Specifically, the 2s₁⁄₂ state lies about 1057 MHz higher than the 2p₁⁄₂ state, even though Dirac's relativistic equation puts them at exactly the same energy.

That "1057 MHz" is an energy expressed as a frequency through E = hf. In other units it is roughly 4.4 microelectronvolts — about a millionth of the energy that binds the electron in the n = 2 shell. It is a whisper. And explaining that whisper rebuilt physics.

The cause is the quantum electrodynamic (QED) vacuum. Empty space is not empty — it seethes with fluctuating electromagnetic fields and virtual particles. The electron interacts with this churn, and the interaction nudges its energy. Because an s-electron and a p-electron sample the vacuum differently near the nucleus, they are nudged by different amounts, and the degeneracy breaks.

How it works

Start with what Dirac theory says. Solve the Dirac equation for an electron in the hydrogen Coulomb potential, and the energy levels come out depending only on the principal quantum number n and the total angular momentum j — never on the orbital angular momentum l by itself. The 2s₁⁄₂ state has n = 2, l = 0, j = 1⁄2. The 2p₁⁄₂ state has n = 2, l = 1, j = 1⁄2. Same n, same j. Dirac says: same energy. Degenerate.

Now turn on the vacuum. In QED the electromagnetic field has a ground state that is not quiet — it has zero-point fluctuations in every mode. The bound electron is constantly buffeted by these fluctuating fields, so its position does a tiny random jitter, sometimes called Zitterbewegung-flavored smearing. Instead of sitting at a sharp point, the electron is effectively spread over a small fuzzy ball.

A smeared charge feels an averaged potential. Near the nucleus the Coulomb potential −e²/r is sharply peaked, so averaging it over the jitter weakens the binding right where it is strongest. The effect is largest exactly at the nucleus — and only s-states have a nonzero wavefunction amplitude there. The 2s electron, with finite density at r = 0, gets shifted upward; the 2p electron, whose density vanishes at the nucleus, barely notices. The levels split.

The full QED accounting has three named pieces:

  • Electron self-energy. The electron emits and reabsorbs virtual photons (it interacts with its own field). This is the dominant term and pushes the 2s level up by the bulk of the shift.
  • Vacuum polarization (Uehling term). A virtual electron–positron pair briefly screens the nucleus, modifying the potential at very short range. This contributes a small piece of the opposite sign for s-states.
  • Anomalous magnetic moment. The same vacuum coupling that gives the electron its g − 2 also feeds a small spin-orbit-like contribution.

Add them up with the right signs and the answer is ≈ 1057 MHz — exactly what Lamb and Retherford measured.

A worked example — Bethe's back-of-the-envelope

In June 1947, on a train home from the Shelter Island conference, Hans Bethe did the first estimate. Here is the spirit of it with concrete numbers.

The mean-square jitter of the electron driven by vacuum fluctuations scales as:

⟨(Δr)²⟩  ≈  (2α/π) · (ℏ/mc)² · ln(K_max / K_min)

where α ≈ 1/137 is the fine-structure constant, ℏ/mc ≈ 3.9 × 10⁻¹³ m is the electron's reduced Compton wavelength, and the logarithm runs over the range of photon energies that matter. The smearing changes the energy of an s-state by an amount proportional to ⟨(Δr)²⟩ times the curvature of the potential at the nucleus:

ΔE  ≈  (1/6) ⟨(Δr)²⟩ · ⟨∇²V⟩  =  (1/6) ⟨(Δr)²⟩ · 4π e² |ψ(0)|²

For the 2s state of hydrogen, |ψ(0)|² = 1/(8πa₀³) with a₀ the Bohr radius. The naive integral over photon energies diverges at the high end (K_max → ∞), which is the famous infinity. Bethe's move: cut it off at the electron rest energy mc² ≈ 511 keV, because above that the non-relativistic treatment is meaningless and the divergence is really part of the electron's observed mass. With that cutoff the logarithm becomes a finite number of order 7–8, and the arithmetic yields:

ΔE(2s − 2p)  ≈  1040 MHz   (Bethe, 1947, non-relativistic)
ΔE(2s − 2p)  ≈  1057 MHz   (full QED + experiment)

Getting 1040 out of 1057 from a train-ride estimate is one of the great calculations in physics. It proved the vacuum was the answer, and it showed that the infinity was tameable.

The numbers in context

Effect in hydrogenOriginTypical size (frequency)Described by
Gross structure (n levels)Coulomb binding~10¹⁵ Hz (eV scale)Bohr / Schrödinger
Fine structure (spin-orbit)Relativity + spin~10 GHz (α² of gross)Dirac equation
Lamb shift (2s–2p)QED vacuum fluctuations≈ 1057 MHzQED radiative correction
Hyperfine (21 cm line)Nuclear spin moment1420 MHzMagnetic dipole coupling
Electron g − 2QED vacuum (same family)shifts moment by ~0.1%QED radiative correction
2s natural linewidthMetastable lifetime ~0.12 s~1 Hz (very narrow)Two-photon decay

The Lamb shift sits in a beautiful gap: smaller than fine structure (which Dirac already explained), comparable to but distinct from hyperfine structure (which comes from the proton, not the vacuum), and produced by the same vacuum coupling that gives the electron its anomalous magnetic moment.

How Lamb and Retherford actually measured it

The genius of the 1947 experiment was avoiding optical spectroscopy entirely. Optically, the 2s–2p gap would be buried in the Doppler width of an emission line. Instead they used a microwave/radio-frequency resonance method born from wartime radar.

  • Produce a beam of hydrogen atoms and excite some to the metastable 2s₁⁄₂ state, which lives ~0.12 s — astronomically long for an excited atom, because it cannot decay to 1s by a single dipole photon.
  • Detect those metastable atoms downstream by the electrons they eject from a metal target.
  • Bathe the beam in tunable microwaves. When the frequency matches the 2s → 2p transition (the Lamb shift, ~1 GHz), atoms are driven into the short-lived 2p state, which decays to 1s almost instantly. The metastable beam is depleted, and the detector signal drops.
  • The microwave frequency at the dip is the Lamb shift. Reading it off gave ≈ 1000 MHz immediately, refined to ~1057 MHz.

Measuring an energy difference directly as a microwave frequency is exquisitely precise — far better than subtracting two optical lines. That precision is what made the small disagreement with Dirac theory undeniable.

Why this gave us renormalization

The self-energy integral diverges. Sum the electron's interaction with all vacuum photon modes and you get infinity, because there is no upper limit on photon energy. Before 1947 this was treated as an embarrassing flaw in quantum field theory. The Lamb shift turned it into a method.

Bethe's argument, sharpened by Kramers and then made rigorous by Schwinger, Feynman, Tomonaga, and Dyson, runs like this. A free electron also has a divergent self-energy — but we never see it separately, because that infinity is already baked into the mass we measure on the scale. The number m in our equations is the "bare" mass; the physical mass is bare mass plus the (infinite) self-energy. So when we compute the energy of a bound electron, we should subtract off the free-electron self-energy, keeping only the difference. The infinities cancel, and the leftover is finite — and it is the Lamb shift.

This subtraction — absorbing divergences into the definitions of measured quantities (mass, charge) — is renormalization. The Lamb shift was its first triumphant test: a finite prediction where there had been an infinity, agreeing with experiment. Everything in the Standard Model that we can calculate, from g − 2 to particle scattering at the LHC, rests on the renormalization program the Lamb shift validated.

Variants and regimes

  • Heavier hydrogen-like ions. The Lamb shift scales steeply with nuclear charge Z (roughly as Z⁴ for the leading self-energy, even faster for higher-order terms). In hydrogen-like uranium (Z = 92) the 1s Lamb shift reaches ~460 eV — a swing of about a hundred million in magnitude versus hydrogen — making heavy ions a stringent test of QED in ultra-strong fields.
  • Muonic hydrogen. Replace the electron with a muon (207× heavier) and the orbit shrinks ~200-fold, so the muon spends far more time inside the proton. The Lamb shift in muonic hydrogen is dominated by the proton's finite size, which is how the "proton radius puzzle" was uncovered around 2010.
  • Helium and few-electron systems. Lamb shifts in helium fine structure provide independent high-precision determinations of the fine-structure constant α.
  • Higher levels and circular Rydberg states. The shift falls off rapidly with n and with l, since both reduce the wavefunction density at the nucleus.
  • Bound-state QED today. Modern calculations include two-loop self-energy and recoil corrections; theory and experiment now agree to better than a part in 10⁵, limited mainly by knowledge of the proton charge radius.

Common pitfalls and misconceptions

  • "The Lamb shift is just fine structure." No. Fine structure splits different j and is contained in Dirac theory. The Lamb shift lifts a degeneracy that Dirac theory leaves intact (same n, same j, different l). It is a genuinely new, QED-only effect.
  • "It comes from the proton." That is hyperfine structure. The Lamb shift is overwhelmingly about the electron and the vacuum; the proton enters only as a small finite-size correction (large only for muonic atoms).
  • "Vacuum fluctuations are a metaphor." They produce a measured 1057 MHz number that matches calculation to five digits. The Casimir effect and g − 2 corroborate the same vacuum. It is as real as any energy in physics.
  • "The infinities mean QED is broken." The infinities are an artifact of using unobservable bare parameters. Renormalization removes them systematically, and the Lamb shift was the proof it works.
  • "The 2s level goes up because the electron gains energy from the photons." More precisely, the smearing weakens the binding near the nucleus, raising the energy (making it less bound). It is a change in how the electron samples the potential, not a simple energy deposit.
  • "You can see it in an ordinary spectrometer." Doppler broadening swamps a 1 GHz splitting on an optical line. It took a radio-frequency resonance method (and later Doppler-free laser spectroscopy) to resolve it.

Why it matters — applications and legacy

  • It launched modern QED. The Lamb shift, together with the electron g − 2, forced the development of renormalized quantum electrodynamics — still the most precisely tested theory in science.
  • It validated renormalization. The technique that tames infinities everywhere in particle physics earned its first experimental confirmation here.
  • Fundamental constants. Lamb-shift and related measurements feed the determination of α and the Rydberg constant to extreme precision.
  • The proton radius puzzle. Muonic-hydrogen Lamb-shift spectroscopy revealed a discrepancy in the proton's measured size, driving a decade of experiment and theory.
  • Strong-field QED. Lamb shifts in highly charged ions probe quantum electrodynamics in the most intense electric fields available in the lab.
  • Tests for new physics. Any deviation between measured and calculated Lamb shifts constrains exotic forces and particles, since QED's prediction is so sharp.

Derivation analysis — where the precision comes from

The leading self-energy shift for an s-state in hydrogen can be written schematically as:

ΔE_self  =  (4 α⁵ / 3π n³) · mc² · [ ln(1/α²) + small terms ]   (s-states)

Two features drive the answer. First, the α⁵ scaling: the shift is suppressed by five powers of the fine-structure constant relative to the rest energy, which is precisely why 13.6 eV of binding produces only a 4 µeV splitting. Second, the large logarithm ln(1/α²) ≈ ln(137²) ≈ 9.8 enhances the bare α⁵ estimate; this log is the finite remnant of Bethe's cutoff, the survivor of the cancelled infinity.

To reach modern accuracy you keep going: the Bethe logarithm (a numerically evaluated average excitation energy, ≈ 2.98 for 2s), two-loop QED terms of order α²(Zα)⁴, recoil corrections of order m/M (the proton is not infinitely heavy), and the finite nuclear size term proportional to the proton charge radius squared. The first three are pure theory and known to many digits; the last is the dominant uncertainty. That is why precise hydrogen and muonic-hydrogen Lamb-shift measurements double as the best probes of the proton's radius — the QED part is so clean that whatever is left over is the nucleus.

Frequently asked questions

What is the Lamb shift?

The Lamb shift is a tiny energy difference between the 2s₁⁄₂ and 2p₁⁄₂ levels of the hydrogen atom — about 1057 MHz, expressed as a frequency. Dirac's relativistic theory of the electron predicts these two states should have exactly the same energy (be degenerate). Experiment shows they don't: the 2s level sits slightly higher. The gap is the Lamb shift, and it forced physicists to take the quantum vacuum seriously.

Why does Dirac theory predict 2s and 2p to be degenerate?

In the Dirac equation for hydrogen, energy levels depend only on the principal quantum number n and the total angular momentum j, not separately on the orbital quantum number l. The 2s₁⁄₂ state (l=0) and the 2p₁⁄₂ state (l=1) both have n=2 and j=1⁄2, so Dirac predicts identical energies. The Lamb shift breaks this accidental degeneracy because the QED vacuum treats an s-electron (which has nonzero probability at the nucleus) differently from a p-electron (which vanishes there).

What causes the Lamb shift?

The dominant cause is the electron's self-energy: the electron continuously emits and reabsorbs virtual photons, and the quantum vacuum's fluctuating electromagnetic field jiggles the electron's position. This smears its charge over a small region, weakening how strongly it feels the nucleus's potential. Because the 2s electron spends time right at the nucleus and the 2p electron does not, the smearing raises the 2s energy more, splitting the levels. Vacuum polarization (the Uehling term) and the electron's anomalous magnetic moment contribute smaller amounts.

How big is the Lamb shift, in numbers?

The 2s₁⁄₂–2p₁⁄₂ splitting in hydrogen is about 1057 MHz — roughly 4.4 microelectronvolts, or about one part in a million of the 2p binding energy. As a fraction of the hydrogen ground-state energy (13.6 eV) it is around 3 in 10 million. Hans Bethe's first non-relativistic estimate in 1947 gave about 1040 MHz; the full QED calculation now agrees with measurement to better than a part in 100,000.

Who discovered the Lamb shift and won the Nobel Prize?

Willis Lamb and his graduate student Robert Retherford measured the splitting in 1947 at Columbia using microwave techniques developed during wartime radar research. Hans Bethe produced the first theoretical estimate within weeks. Lamb shared the 1955 Nobel Prize in Physics for the discovery (with Polykarp Kusch). The full theory by Feynman, Schwinger, Tomonaga and Dyson earned a later 1965 Nobel for QED itself.

How is the Lamb shift connected to renormalization?

A naive calculation of the self-energy gives infinity — the electron interacts with arbitrarily high-energy vacuum modes. Bethe's insight was that part of that infinity is already absorbed into the electron's observed mass; only the difference between a bound electron and a free electron is physical, and that difference is finite. Subtracting infinities by reassigning them to measured quantities is renormalization. The Lamb shift was the first place this trick produced a finite, correct, testable number.

Is the Lamb shift the same as fine structure or hyperfine structure?

No. Fine structure (spin-orbit coupling plus relativistic corrections) splits levels of different j and is fully described by the Dirac equation. Hyperfine structure comes from the nuclear spin's magnetic moment and is smaller still (the famous 1420 MHz 21 cm line). The Lamb shift is a QED radiative correction that lifts the degeneracy between same-j, different-l states that Dirac theory leaves untouched. It sits between fine and hyperfine structure in size, around 1057 MHz.

Why is the Lamb shift larger for s-states than p-states?

Only s-states (l=0) have a nonzero wavefunction amplitude at the nucleus. The vacuum fluctuations smear the electron's position, which matters most where the Coulomb potential is steepest — right at the nucleus. An s-electron samples that region; a p-electron's wavefunction vanishes there. So the self-energy correction shifts s-levels much more than p-levels, which is exactly why the 2s level rises above the 2p level instead of staying degenerate with it.