# Lamb shift

In physics, the Lamb shift, named after Willis Lamb, is a difference in energy between two energy levels 2S1/2 and 2P1/2 (in term symbol notation) of the hydrogen atom which was not predicted by the Dirac equation, according to which these states should have the same energy.

Interaction between vacuum energy fluctuations and the hydrogen electron in these different orbitals is the cause of the Lamb shift, as was shown subsequent to its discovery. The Lamb shift has since played a significant role through vacuum energy fluctuations in theoretical prediction of Hawking radiation from black holes.

This effect was first measured in 1947 in the Lamb–Retherford experiment on the hydrogen microwave spectrum[1] and this measurement provided the stimulus for renormalization theory to handle the divergences. It was the harbinger of modern quantum electrodynamics developed by Julian Schwinger, Richard Feynman, Ernst Stueckelberg, Sin-Itiro Tomonaga and Freeman Dyson. Lamb won the Nobel Prize in Physics in 1955 for his discoveries related to the Lamb shift.

## Importance

On Lamb's 65th birthday, Freeman Dyson addressed him as follows: "Those years, when the Lamb shift was the central theme of physics, were golden years for all the physicists of my generation. You were the first to see that this tiny shift, so elusive and hard to measure, would clarify our thinking about particles and fields."[2]

## Derivation

This heuristic derivation of the electrodynamic level shift following Welton is from Quantum Optics.[3]

The fluctuation in the electric and magnetic fields associated with the QED vacuum perturbs the electric potential due to the atomic nucleus. This perturbation causes a fluctuation in the position of the electron, which explains the energy shift. The difference of potential energy is given by

${\displaystyle \Delta V=V({\vec {r}}+\delta {\vec {r}})-V({\vec {r}})=\delta {\vec {r}}\cdot \nabla V({\vec {r}})+{\frac {1}{2}}(\delta {\vec {r}}\cdot \nabla )^{2}V({\vec {r}})+\cdots }$

Since the fluctuations are isotropic,

${\displaystyle \langle \delta {\vec {r}}\rangle _{vac}=0,}$
${\displaystyle \langle (\delta {\vec {r}}\cdot \nabla )^{2}\rangle _{vac}={\frac {1}{3}}\langle (\delta {\vec {r}})^{2}\rangle _{vac}\nabla ^{2}.}$

So one can obtain

${\displaystyle \langle \Delta V\rangle ={\frac {1}{6}}\langle (\delta {\vec {r}})^{2}\rangle _{vac}\left\langle \nabla ^{2}\left({\frac {-e^{2}}{4\pi \epsilon _{0}r}}\right)\right\rangle _{at}.}$

The classical equation of motion for the electron displacement (δr)k induced by a single mode of the field of wave vector k and frequency ν is

${\displaystyle m{\frac {d^{2}}{dt^{2}}}(\delta r)_{\vec {k}}=-eE_{\vec {k}},}$

and this is valid only when the frequency ν is greater than ν0 in the Bohr orbit, ${\displaystyle \nu >\pi c/a_{0}}$. The electron is unable to respond to the fluctuating field if the fluctuations are smaller than the natural orbital frequency in the atom.

For the field oscillating at ν,

${\displaystyle \delta r(t)\cong \delta r(0)e^{-i\nu t}+c.c.,}$

therefore

${\displaystyle (\delta r)_{\vec {k}}\cong {\frac {e}{mc^{2}k^{2}}}E_{\vec {k}}={\frac {e}{mc^{2}k^{2}}}{\mathcal {E}}_{\vec {k}}\left(a_{\vec {k}}e^{-i\nu t+i{\vec {k}}\cdot {\vec {r}}}+h.c.\right)\qquad {\text{with}}\qquad {\mathcal {E}}_{\vec {k}}=\left({\frac {\hbar ck/2}{\epsilon _{0}\Omega }}\right)^{1/2},}$

where ${\displaystyle \Omega }$ is some large normalization volume (the volume of the hypothetical "box" containing the hydrogen atom). By the summation over all ${\displaystyle {\vec {k}},}$

{\displaystyle {\begin{aligned}\langle (\delta {\vec {r}})^{2}\rangle _{vac}&=\sum _{\vec {k}}\left({\frac {e}{mc^{2}k^{2}}}\right)^{2}\left\langle 0\left|(E_{\vec {k}})^{2}\right|0\right\rangle \\&=\sum _{\vec {k}}\left({\frac {e}{mc^{2}k^{2}}}\right)^{2}\left({\frac {\hbar ck}{2\epsilon _{0}\Omega }}\right)\\&=2{\frac {\Omega }{(2\pi )^{3}}}4\pi \int dkk^{2}\left({\frac {e}{mc^{2}k^{2}}}\right)^{2}\left({\frac {\hbar ck}{2\epsilon _{0}\Omega }}\right)&&{\text{since continuity of }}{\vec {k}}{\text{ implies }}\sum _{\vec {k}}\to 2{\frac {\Omega }{(2\pi )^{3}}}\int d^{3}k\\&={\frac {1}{2\epsilon _{0}\pi ^{2}}}\left({\frac {e^{2}}{\hbar c}}\right)\left({\frac {\hbar }{mc}}\right)^{2}\int {\frac {dk}{k}}\end{aligned}}}

This result diverges when no limits about the integral (at both large and small frequencies). As mentioned above, this method is expected to be valid only when ${\displaystyle \nu >\pi c/a_{0}}$, or equivalently ${\displaystyle k>\pi /a_{0}}$. It is also valid only for wavelengths longer than the Compton wavelength, or equivalently ${\displaystyle k. Therefore, one can choose the upper and lower limit of the integral and these limits make the result converge.

${\displaystyle \langle (\delta {\vec {r}})^{2}\rangle _{vac}\cong {\frac {1}{2\epsilon _{0}\pi ^{2}}}\left({\frac {e^{2}}{\hbar c}}\right)\left({\frac {\hbar }{mc}}\right)^{2}\ln {\frac {4\epsilon _{0}\hbar c}{e^{2}}}}$.

For the atomic orbital and the Coulomb potential,

${\displaystyle \left\langle \nabla ^{2}\left({\frac {-e^{2}}{4\pi \epsilon _{0}r}}\right)\right\rangle _{at}={\frac {-e^{2}}{4\pi \epsilon _{0}}}\int d{\vec {r}}\psi ^{*}({\vec {r}})\nabla ^{2}\left({\frac {1}{r}}\right)\psi ({\vec {r}})={\frac {e^{2}}{\epsilon _{0}}}|\psi (0)|^{2},}$

since it is known that

${\displaystyle \nabla ^{2}\left({\frac {1}{r}}\right)=-4\pi \delta ({\vec {r}}).}$

For p orbitals, the nonrelativistic wave function vanishes at the origin, so there is no energy shift. But for s orbitals there is some finite value at the origin,

${\displaystyle \psi _{2S}(0)={\frac {1}{(8\pi a_{0}^{3})^{1/2}}},}$

${\displaystyle a_{0}={\frac {4\pi \epsilon _{0}\hbar ^{2}}{me^{2}}}.}$

Therefore,

${\displaystyle \left\langle \nabla ^{2}\left({\frac {-e^{2}}{4\pi \epsilon _{0}r}}\right)\right\rangle _{at}={\frac {e^{2}}{\epsilon _{0}}}|\psi _{2S}(0)|^{2}={\frac {e^{2}}{8\pi \epsilon _{0}a_{0}^{3}}}}$.

Finally, the difference of the potential energy becomes:

${\displaystyle \langle \Delta V\rangle ={\frac {4}{3}}{\frac {e^{2}}{4\pi \epsilon _{0}}}{\frac {e^{2}}{4\pi \epsilon _{0}\hbar c}}\left({\frac {\hbar }{mc}}\right)^{2}{\frac {1}{8\pi a_{0}^{3}}}\ln {\frac {4\epsilon _{0}\hbar c}{e^{2}}}=\alpha ^{5}mc^{2}{\frac {1}{6\pi }}\ln {\frac {1}{\pi \alpha }},}$

where ${\displaystyle \alpha }$ is the fine-structure constant. This shift is about 1 GHz, very similar with the observed energy shift.

Welton's heuristic derivation of the Lamb shift is similar to, but distinct from, the calculation of the Darwin term using Zitterbewegung, a contribution to the fine structure that is of lower order in ${\displaystyle \alpha }$ than the Lamb shift.[4]:80–81

## Experimental work

In 1947 Willis Lamb and Robert Retherford carried out an experiment using microwave techniques to stimulate radio-frequency transitions between 2S1/2 and 2P1/2 levels of hydrogen.[5] By using lower frequencies than for optical transitions the Doppler broadening could be neglected (Doppler broadening is proportional to the frequency). The energy difference Lamb and Retherford found was a rise of about 1000 MHz of the 2S1/2 level above the 2P1/2 level.

This particular difference is a one-loop effect of quantum electrodynamics, and can be interpreted as the influence of virtual photons that have been emitted and re-absorbed by the atom. In quantum electrodynamics the electromagnetic field is quantized and, like the harmonic oscillator in quantum mechanics, its lowest state is not zero. Thus, there exist small zero-point oscillations that cause the electron to execute rapid oscillatory motions. The electron is "smeared out" and each radius value is changed from r to r + δr (a small but finite perturbation).

The Coulomb potential is therefore perturbed by a small amount and the degeneracy of the two energy levels is removed. The new potential can be approximated (using atomic units) as follows:

${\displaystyle \langle E_{\mathrm {pot} }\rangle =-{\frac {Ze^{2}}{4\pi \epsilon _{0}}}\left\langle {\frac {1}{r+\delta r}}\right\rangle .}$

The Lamb shift itself is given by

${\displaystyle \Delta E_{\mathrm {Lamb} }=\alpha ^{5}m_{e}c^{2}{\frac {k(n,0)}{4n^{3}}}\ \mathrm {for} \ \ell =0\,}$

with k(n, 0) around 13 varying slightly with n, and

${\displaystyle \Delta E_{\mathrm {Lamb} }=\alpha ^{5}m_{e}c^{2}{\frac {1}{4n^{3}}}\left[k(n,\ell )\pm {\frac {1}{\pi (j+{\frac {1}{2}})(\ell +{\frac {1}{2}})}}\right]\ \mathrm {for} \ \ell \neq 0\ \mathrm {and} \ j=\ell \pm {\frac {1}{2}},}$

with k(n,) a small number (< 0.05).

For a derivation of ΔELamb see for example:[6]

## In the hydrogen spectrum

In 1947, Hans Bethe was the first to explain the Lamb shift in the hydrogen spectrum, and he thus laid the foundation for the modern development of quantum electrodynamics. Bethe was able to derive the Lamb shift by implementing the idea of mass renormalization, which allowed him to calculate the observed energy shift as the difference between the shift of a bound electron and the shift of a free electron. [7] The Lamb shift currently provides a measurement of the fine-structure constant α to better than one part in a million, allowing a precision test of quantum electrodynamics.