Thomson Scattering

The cross section for Thomson scattering illustrates the need for ``negative energy'' states in our calculations. Recall that we got the correct cross section from the non-relativistic calculation and that Thomson also got the correct result from classical E&M.

In the Dirac theory, we have only one term in the interaction Hamiltonian,

\begin{displaymath}\bgroup\color{black} H_{int}=ie\gamma_4\gamma_k A_k \egroup\end{displaymath}

Because it is linear in \bgroup\color{black}$A$\egroup it can create a photon or annihilate a photon. Photon scattering is therefore second order (and proportional to \bgroup\color{black}$e^2$\egroup). The quantized photon field is

\begin{displaymath}\bgroup\color{black} A_\mu(x)={1\over\sqrt{V}}\sum\limits_{k\...
...rho}+a_{k,\alpha}^\dagger(0)e^{-ik_\rho x_\rho}\right). \egroup\end{displaymath}

The initial and final states are definite momentum states, as are the intermediate electron states. We shall first do the calculation assuming no electrons from the ``negative energy'' sea participate, other than to exclude transitions to those ``negative energy'' states. The initial and final states are therefore the positive energy plane wave states \bgroup\color{black}$\psi^{(r)}_{\vec{p}}$\egroup for \bgroup\color{black}$r=1,2$\egroup. The intermediate states must also be positive energy states since the ``negative energy'' states are all filled.

The computation of the scattering cross section follows the same steps made in the development of the Krammers-Heisenberg formula for photon scattering. There is no \bgroup\color{black}$A^2$\egroup term so we are just computing the two second order terms.

\begin{eqnarray*}
c_{\vec{p'},r';\vec{k'},\hat{\epsilon'}}^{(2)}(t)
&=&{-e^2\ove...
...t\limits_0^t\;dt_2 e^{i(E'-E+\hbar\omega'-\hbar\omega)t_2/\hbar}
\end{eqnarray*}


As in the earlier calculation, the photon states have been eliminated from the equation since they give a factor of 1 with the initial state photon being annihilated and the final state photon being created in each term.

Now lets take a look at one of the matrix elements. Assume the initial state electron is at rest and that the photon momentum is small.

\begin{displaymath}\bgroup\color{black} \langle\vec{p''}r''\vert i\gamma_4\gamma_n\epsilon_ne^{i\vec{k}\cdot\vec{x}}\vert\vec{p}r\rangle \egroup\end{displaymath}

For \bgroup\color{black}$\vec{p}=0$\egroup and \bgroup\color{black}$\vec{k}=0$\egroup, a delta function requires that \bgroup\color{black}$\vec{p''}=0$\egroup. It turns out that \bgroup\color{black}$u^{(r'')\dagger}_0\gamma_4\gamma_n u^{(r)}_0=0$\egroup, so that the cross section is zero in this limit.

\begin{eqnarray*}
\gamma_i&=&\pmatrix{0 & -i\sigma_i \cr i\sigma_i & 0 \cr} \\
...
...amma_4\gamma_i&=&\pmatrix{0 & -i\sigma_i \cr -i\sigma_i & 0 \cr}
\end{eqnarray*}


This matrix only connects \bgroup\color{black}$r=1,2$\egroup spinors to \bgroup\color{black}$r=3,4$\egroup spinors because of its off diagonal nature. So, the calculation yields zero for a cross section in contradiction to the other two calculations. In fact, since the photon momentum is not quite zero, there is a small contribution, but far too small.

The above calculation misses some important terms due to the ``negative energy'' sea. There are additional terms if we consider the possibility that the photon can elevate a ``negative energy'' electron to have positive energy.

\epsfig{file=figs/Thompson.eps,height=3in}
In one term, the initial state photon is absorbed by a ``negative energy'' electron, then the initial state electron fills the hole in the ``negative energy'' sea emitting the final state photon. In the other term, even further from the mass shell, a ``negative energy'' electron emits the final state photon and moves to a positive energy state, then the initial state electron absorbs the initial photon and fills the hole left behind in the sea. These terms are larger because the \bgroup\color{black}$\gamma_4\gamma_i$\egroup matrix connects positive energy and ``negative energy'' states.

\begin{eqnarray*}
c_{\vec{p'},r';\vec{k'}\hat{\epsilon'}}^{(2)}(t)
&=&{ie^2c^2\o...
...t\limits_0^t\;dt_2 e^{i(E'-E+\hbar\omega'-\hbar\omega)t_2/\hbar}
\end{eqnarray*}


The matrix element is to be taken with the initial electron at rest, \bgroup\color{black}$\hbar k«mc$\egroup, the final electron (approximately) at rest, and hence the intermediate electron at rest, due to a delta function of momentum conservation that comes out of the spatial integral.

Let the positive energy spinors be written as

\begin{displaymath}\bgroup\color{black} u^{(r)}_0=\pmatrix{\chi^{(r)}\cr 0} \egroup\end{displaymath}

and the ``negative energy'' spinors as

\begin{displaymath}\bgroup\color{black} u^{(r'')}_0=\pmatrix{0\cr \chi^{(r'')}} \egroup\end{displaymath}

The matrix \bgroup\color{black}$\gamma_4\gamma_i$\egroup connect the positive and ``negative energy'' spinors so that the amplitude can be written in terms of two component spinors and Pauli matrices.

\begin{eqnarray*}
\gamma_4\gamma_i&=&\pmatrix{0 & -i\sigma_i \cr -i\sigma_i & 0 ...
...mits_0^t\;dt_2 e^{i(E'-E+\hbar\omega'-\hbar\omega)t_2/\hbar} \\
\end{eqnarray*}


\begin{eqnarray*}
c_{\vec{p'},r';\vec{k'}\hat{\epsilon'}}^{(2)}(t)
&=&{ie^2\over...
...
\vert\hat{\epsilon}\cdot\hat{\epsilon'}\vert^2\delta_{rr'} \\
\end{eqnarray*}


This agrees with the other calculations and with experiment. The ``negative energy'' sea is required to get the right answer in Dirac theory. There are alternatives to the ``negative energy'' sea. Effectively we are allowing for the creation of electron positron pairs through the use of the filled negative energy states. The same result could be obtained with the possibility of pair creation, again with the hypothesis that a positron is a ``negative energy'' electron going backward in time.

Jim Branson 2013-04-22