The Lamb Shift

In 1947, Willis E. Lamb and R. C. Retherford used microwave techniques to determine the splitting between the $2S_{1\over 2}$ and $2P_{1\over 2}$ states in Hydrogen to have a frequency of 1.06 GHz, (a wavelength of about 30 cm). (The shift is now accurately measured to be 1057.864 MHz.) This is about the same size as the hyperfine splitting of the ground state.
\epsfig{file=figs/lamb.eps,height=3in}

The technique used was quite interesting. They made a beam of Hydrogen atoms in the $2S_{1\over 2}$ state, which has a very long lifetime because of selection rules. Microwave radiation with a (fixed) frequency of 2395 MHz was used to cause transitions to the \bgroup\color{black}$2P_{3\over 2}$\egroup state and a magnetic field was adjusted to shift the energy of the states until the rate was largest. The decay of the $2P_{3\over 2}$ state to the ground state was observed to determine the transition rate. From this, they were able to deduce the shift between the \bgroup\color{black}$2S_{1\over 2}$\egroup and \bgroup\color{black}$2P_{1\over 2}$\egroup states.

Hans Bethe used non-relativistic quantum mechanics to calculate the self-energy correction to account for this observation.

\epsfig{file=figs/Bethe.eps,height=3in}
We now can compute the correction the same way he did.

\begin{eqnarray*}
\Delta E_n^{(obs)}&=&\Delta E_n+C\langle n\vert p^2\vert n\ran...
...vert\langle n\vert\vec{p}\vert j\rangle\vert^2\right)d\omega \\
\end{eqnarray*}


It is now necessary to discuss approximations needed to complete this calculation. In particular, the electric dipole approximation will be of great help, however, it is certainly not warranted for large photon energies. For a good E1 approximation we need \bgroup\color{black}$E_\gamma«1973$\egroup eV. On the other hand, we want the cut-off for the calculation to be of order \bgroup\color{black}$w_{cut-off}\approx mc^2/\hbar$\egroup. We will use the E1 approximation and the high cut-off, as Bethe did, to get the right answer. At the end, the result from a relativistic calculation can be tacked on to show why it turns out to be the right answer. (We aren't aiming for the worlds best calculation anyway.)

\begin{eqnarray*}
\Delta E^{(obs)}_n&=&{2\alpha\hbar\over 3\pi m^2c^2}\int\limit...
...t-off}}\right)\vert\langle n\vert\vec{p}\vert j\rangle\vert^2\\
\end{eqnarray*}


The log term varies more slowly than does the rest of the terms in the sum. We can approximate it by an average. Bethe used numerical calculations to determine that the effective average of \bgroup\color{black}$\hbar\omega_{nj}$\egroup is \bgroup\color{black}$8.9\alpha^2mc^2$\egroup. We will do the same and pull the log term out as a constant.

\begin{eqnarray*}
\Delta E^{(obs)}_n={2\alpha\hbar\over 3\pi m^2c^2}\log\left({\...
...s_j\omega_{nj}\vert\langle n\vert\vec{p}\vert j\rangle\vert^2\\
\end{eqnarray*}


This sum can now be reduced further to a simple expression proportional to the \bgroup\color{black}$\vert\psi_n(0)\vert^2$\egroup using a typical clever quantum mechanics calculation. The basic Hamiltonian for the Hydrogen atom is \bgroup\color{black}$H_0={p^2\over 2m}+V(r)$\egroup.

\begin{eqnarray*}[\vec{p},H_0]&=&[\vec{p},V]={\hbar\over i}\vec{\nabla}V \\
\la...
...ec{p}\vert j\rangle\langle j\vert\vec{\nabla}V\vert n\rangle \\
\end{eqnarray*}


This must be a real number so we may use its complex conjugate.

\begin{eqnarray*}
\left(\sum\limits_j (E_i-E_n)\langle n\vert\vec{p}\vert j\rang...
...rt n\rangle \\
&=&-{e^2\hbar^2\over 2}\vert\psi_n(0)\vert^2 \\
\end{eqnarray*}


Only the s states will have a non-vanishing probability to be at the origin with \bgroup\color{black}$\vert\psi_{n00}(0)\vert^2={1\over\pi n^3a_0^3}$\egroup and \bgroup\color{black}$a_0={\hbar\over\alpha mc}$\egroup. Therefore, only the s states will shift in energy appreciably. The shift will be.

\begin{eqnarray*}
\Delta E^{(obs)}_n&=&-{2\alpha\hbar\over 3\pi m^2c^2}\log\left...
...pi^2\hbar c}\log\left({1\over 8.9\alpha^2}\right)=1.041 GHz \\
\end{eqnarray*}


This agrees far too well with the measurement, considering the approximations made and the dependence on the cut-off. There is, however, justification in the relativistic calculation. Typically, the full calculation was made by using this non-relativistic approach up to some energy of the order of \bgroup\color{black}$\alpha mc^2$\egroup, and using the relativistic calculation above that. The relativistic free electron self-energy correction diverges only logarithmically and a very high cutoff can be usedwithout a problem. The mass of the electron is renormalized as above. The Lamb shift does not depend on the cutoff and hence it is well calculated. We only need the non-relativistic part of the calculation up to photon energies for which the E1 approximation is OK. The relativistic part of the calculation down to \bgroup\color{black}$\omega_{min}$\egroup yields.

\begin{displaymath}\bgroup\color{black}\Delta E_n={4\alpha^5\over 3\pi n^3}\left...
...r\omega_{min}}\right)+{11\over 24}-{1\over 5}\right)mc^2\egroup\end{displaymath}

The non-relativistic calculation gave.

\begin{displaymath}\bgroup\color{black}\Delta E_n={4\alpha^5\over 3\pi n^3}\log\...
...omega_{min}\over\vert\bar{\omega}_{nj}\vert}\right)mc^2 \egroup\end{displaymath}

So the sum of the two gives.

\begin{displaymath}\bgroup\color{black}\Delta E_n^{(obs)}
={4\alpha^5\over 3\pi ...
...{\omega}_{nj}}\right)+{11\over 24}-{1\over 5}\right)mc^2\egroup\end{displaymath}

The dependence on $\omega_{min}$ cancels. In this calculation, the \bgroup\color{black}$mc^2$\egroup in the log is the outcome of the relativistic calculation, not the cutoff. The electric dipole approximation is even pretty good since we did not need to go up to large photon energies non-relativistically and no E1 approximation is needed for the relativistic part. That's how we (and Bethe) got about the right answer.

The Lamb shift splits the $2S_{1\over 2}$ and $2P_{1\over 2}$ states which are otherwise degenerate. Its origin is purely from field theory. The experimental measurement of the Lamb shift stimulated theorists to develop Quantum ElectroDynamics. The correction increases the energy of s states. One may think of the physical origin as the electron becoming less pointlike as virtual photons are emitted and reabsorbed. Spreading the electron out a bit decreases the effect of being in the deepest part of the potential, right at the origin. Based on the energy shift, I estimate that the electron in the 2s state is spread out by about 0.005 Angstroms, much more than the size of the nucleus.

The anomalous magnetic moment of the electron, \bgroup\color{black}$g-2$\egroup, which can also be calculated in field theory, makes a small contribution to the Lamb shift.

Jim Branson 2013-04-22