# $e^+ e^- \to \mu^+ \mu^-$ cross section

According to Schwartz (Quantum Field Theory and the Standard Model, section 5.3), the matrix element needed for the calculation of the cross section for the reaction $e^+ e^- \to \mu^+\mu^-$ is given by $$M = \sum_\epsilon \left<s_1s_2|\epsilon\right>\left<\epsilon|s_3s_4\right>,$$ where $\epsilon$ is the polarization of the intermediate photon, $s_1$ and $s_2$ are the spins of the incoming particles and $s_3$ and $s_4$ the spins of the outgoing ones.

There are four possible initial states:

$$\left. |s_1s_2\right> = \left. |\leftrightarrow\leftrightarrow\right>, \left.|\updownarrow\updownarrow\right>, \left.|\updownarrow\leftrightarrow\right>, \left.|\leftrightarrow\updownarrow\right>,$$ meaning that both the electron and the positron are polarized along the $x$-axis or the $y$-axis or one is polarized along the $x$-axis and the other one along the $y$-axis (when moving into the $z$-direction). Since the photon is spin 1, only the first two states could annihilate into a photon which has the two possible polarizations $$\epsilon_1 = (0,1,0,0)\qquad\text{and}\qquad\epsilon_2 = (0,0,1,0).$$ Similarly, the possible polarizations of the photon looking at the $\mu^+\mu^-$ pair can be $$\bar{\epsilon}_1 = (0,1,0,0)\qquad\text{and}\qquad \bar{\epsilon}_2 = (0,0,\cos\theta, -\sin\theta)$$ with $\theta$ the angle to the $e^+e^-$-axis.

My question is how to explicitly calculate the matrix elements needed for the cross section. According to Schwartz, the only two non-vanishing matrix elements are $$M_1 = \epsilon^1\bar{\epsilon}_1+\epsilon^1\bar{\epsilon}_2 = -1$$ and $$M_2 = \epsilon^2\bar{\epsilon}_1+\epsilon^2\bar{\epsilon}_2 = -\cos\theta,$$ but I don't know how to arrive at this using the formula from above. Hopefully somebody can help.

In a typical collider experiment the momentum vectors are $$p_1=\begin{pmatrix}E\\0\\0\\p\end{pmatrix}\qquad p_2=\begin{pmatrix}E\\0\\0\\-p\end{pmatrix}\qquad p_3=\begin{pmatrix} E\\ \rho\sin\theta\cos\phi\\ \rho\sin\theta\sin\phi\\ \rho\cos\theta \end{pmatrix} \qquad p_4=\begin{pmatrix} E\\ -\rho\sin\theta\cos\phi\\ -\rho\sin\theta\sin\phi\\ -\rho\cos\theta \end{pmatrix}$$

where $$E$$ is beam energy, $$p=\sqrt{E^2-m^2}$$, $$\rho=\sqrt{E^2-M^2}$$, $$m$$ is electron mass $$0.51\,\text{MeV}$$, and $$M$$ is muon mass $$106\,\text{MeV}$$. The spinors are $$\begin{gather*} u_{11}=\begin{pmatrix}E+m\\0\\p\\0\end{pmatrix}\quad v_{21}=\begin{pmatrix}-p\\0\\E+m\\0\end{pmatrix}\quad u_{31}=\begin{pmatrix}E+M\\0\\p_3^z\\p_3^x+ip_3^y\end{pmatrix}\quad v_{41}=\begin{pmatrix}p_4^z\\p_4^x+ip_4^y\\E+M\\0\end{pmatrix} \\ u_{12}=\begin{pmatrix}0\\E+m\\0\\-p\end{pmatrix}\quad v_{22}=\begin{pmatrix}0\\p\\0\\E+m\end{pmatrix}\quad u_{32}=\begin{pmatrix}0\\E+M\\p_3^x-ip_3^y\\-p_3^z\end{pmatrix}\quad v_{42}=\begin{pmatrix}p_4^x-ip_4^y\\-p_4^z\\0\\E+M\end{pmatrix} \end{gather*}$$

The last digit in a spinor subscript is 1 for spin up and 2 for spin down. Note that the spinors are not individually normalized. Instead, a combined spinor normalization constant $$N=(E+m)^2(E+M)^2$$ will be used where needed.

This is the probability density for muon production. Symbol $$s=(p_1+p_2)^2=4E^2$$, symbol $$s_j$$ selects the spin of spinor $$j$$, and $$e$$ is electron charge. $$\begin{equation*} |\mathcal{M}(s_1,s_2,s_3,s_4)|^2 =\frac{e^4}{s^2}\frac{1}{N}\left|(\bar{u}_3\gamma_\mu v_4)(\bar{v}_2\gamma^\mu u_1)\right|^2 \end{equation*}$$

The expected probability density $$\langle|\mathcal{M}|^2\rangle$$ is computed by summing $$|\mathcal{M}|^2$$ over all spin states and dividing by the number of inbound states. There are four inbound states. \begin{align*} \langle|\mathcal{M}|^2\rangle &=\frac{1}{4}\sum_{s_1=1}^2\sum_{s_2=1}^2\sum_{s_3=1}^2\sum_{s_4=1}^2|\mathcal{M}(s_1,s_2,s_3,s_4)|^2 \\ &=\frac{e^4}{4s^2}\sum_{s_1=1}^2\sum_{s_2=1}^2\sum_{s_3=1}^2\sum_{s_4=1}^2 \frac{1}{N}\left|(\bar{u}_3\gamma_\mu v_4)(\bar{v}_2\gamma^\mu u_1)\right|^2 \end{align*}

Another way to compute $$\langle|\mathcal{M}|^2\rangle$$ is to use the Casimir trick. $$\begin{equation*} \langle|\mathcal{M}|^2\rangle =\frac{e^4}{4s^2} \mathop{\rm Tr}\left((\not p_3+M)\gamma^\mu(\not p_4-M)\gamma^\nu\right) \mathop{\rm Tr}\left((\not p_2-m)\gamma_\mu(\not p_1+m)\gamma_\nu\right) \end{equation*}$$

Here is a third way to compute $$\langle|\mathcal{M}|^2\rangle$$. $$\begin{equation*} \langle|\mathcal{M}|^2\rangle =\frac{e^4}{4s^2} \left( 32 (p_1\cdot p_3) (p_2\cdot p_4) + 32 (p_1\cdot p_4) (p_2\cdot p_3) + 32 m^2 (p_3\cdot p_4) + 32 M^2 (p_1\cdot p_2) + 64 m^2 M^2 \right) \end{equation*}$$

For the momentum vectors given above the result is $$\begin{equation*} \langle|\mathcal{M}|^2\rangle =e^4\left(1+\cos^2\theta+\frac{m^2+M^2}{E^2}\sin^2\theta+\frac{m^2M^2}{E^4}\cos^2\theta\right) \end{equation*}$$

The Stanford Linear Collider had a collision energy of $$2E=91$$ GeV. For beam energies such as SLC where $$E\gg M$$ the above equation can be approximated as $$\langle|\mathcal{M}|^2\rangle=e^4(1+\cos^2\theta)$$

The differential cross section is $$\frac{d\sigma}{d\Omega} =\frac{\langle|\mathcal{M}|^2\rangle}{64\pi^2s} =\frac{e^4}{256\pi^2E^2}(1+\cos^2\theta)$$

Recall that $$e^2=4\pi\alpha$$ hence $$\begin{equation*} \frac{d\sigma}{d\Omega}=\frac{\alpha^2}{16E^2}(1+\cos^2\theta) \end{equation*}$$

The total cross section calculation requires the following definite integral. $$\int_\Omega(1+\cos^2\theta)\,d\Omega =\int_0^{2\pi}\int_0^\pi(1+\cos^2\theta)\sin\theta\,d\theta\,d\phi =\frac{8}{3}\int_0^{2\pi}d\phi =\frac{16\pi}{3}$$

Hence the total cross section is $$\sigma =\int_\Omega d\sigma =\int_\Omega\frac{\alpha^2}{16E^2}(1+\cos^2\theta)\,d\Omega =\frac{\alpha^2}{16E^2}\frac{16\pi}{3} =\frac{\pi\alpha^2}{3E^2}$$