Part IV, Chapter 4 | Page 1 of 4

e⁺e⁻ → μ⁺μ⁻ Annihilation

The complete tree-level calculation: from Feynman diagram to total cross section

4.1 The Process and Its Feynman Diagram

We compute the simplest QED scattering process: an electron and positron annihilate into a muon-antimuon pair. At tree level, there is exactly one Feynman diagram — the s-channel, where the $e^+e^-$ pair annihilates into a virtual photon which then creates the $\mu^+\mu^-$ pair.

We label the momenta as: $e^-(p_1) + e^+(p_2) \to \mu^-(p_3) + \mu^+(p_4)$. Since $m_\mu \gg m_e$, we work in the limit $m_e = 0$ but keep $m_\mu$ general initially, then take the high-energy limit $s \gg m_\mu^2$ at the end.

💡Why Only One Diagram?

Unlike e⁺e⁻ → e⁺e⁻ (Bhabha scattering), there is no t-channel diagram here because the electron cannot emit a photon and turn into a muon — lepton flavor is conserved at each QED vertex. The photon must carry the full center-of-mass energy, making this a clean probe of the electromagnetic interaction.

4.2 Writing the Amplitude from Feynman Rules

Reading the diagram from right to left, the QED Feynman rules give us:

Electron vertex: The incoming electron ($u(p_1)$) and positron ($\bar{v}(p_2)$) couple to the photon via $-ie\gamma^\mu$. The electron current is:

$$j^\mu_e = \bar{v}(p_2)(-ie\gamma^\mu)u(p_1)$$

Photon propagator: The virtual photon with momentum $q = p_1 + p_2$ contributes (in Feynman gauge):

$$\frac{-ig_{\mu\nu}}{q^2} = \frac{-ig_{\mu\nu}}{(p_1 + p_2)^2} = \frac{-ig_{\mu\nu}}{s}$$

Muon vertex: The outgoing muon ($\bar{u}(p_3)$) and antimuon ($v(p_4)$) give:

$$j^\nu_\mu = \bar{u}(p_3)(-ie\gamma^\nu)v(p_4)$$

Assembling the amplitude: Contracting everything together:

$$\boxed{i\mathcal{M} = \frac{(-ie)^2}{s}\left[\bar{v}(p_2)\gamma^\mu u(p_1)\right]\left[\bar{u}(p_3)\gamma_\mu v(p_4)\right]}$$

or equivalently:

$$\mathcal{M} = \frac{-e^2}{s}\left[\bar{v}(p_2)\gamma^\mu u(p_1)\right]\left[\bar{u}(p_3)\gamma_\mu v(p_4)\right]$$

4.3 Squaring the Amplitude and Spin Sums

For unpolarized beams and unobserved final-state spins, we compute the spin-averaged squared amplitude:

$$\overline{|\mathcal{M}|^2} = \frac{1}{4}\sum_{\text{spins}}|\mathcal{M}|^2$$

The factor 1/4 comes from averaging over the two spin states of each initial particle (electron and positron). Writing out $|\mathcal{M}|^2 = \mathcal{M}\mathcal{M}^*$:

$$|\mathcal{M}|^2 = \frac{e^4}{s^2}\left[\bar{v}(p_2)\gamma^\mu u(p_1)\right]\left[\bar{u}(p_3)\gamma_\mu v(p_4)\right]\left[\bar{v}(p_2)\gamma^\nu u(p_1)\right]^*\left[\bar{u}(p_3)\gamma_\nu v(p_4)\right]^*$$

Using the identity $[\bar{\psi}_1 \Gamma \psi_2]^* = \bar{\psi}_2 \gamma^0 \Gamma^\dagger \gamma^0 \psi_1 = \bar{\psi}_2 \bar{\Gamma} \psi_1$ where$\bar{\Gamma} = \gamma^0 \Gamma^\dagger \gamma^0$ (and $\overline{\gamma^\nu} = \gamma^\nu$), the spin sums factorize into two independent traces using the completeness relations:

$$\sum_s u^s(p)\bar{u}^s(p) = \not{p} + m, \qquad \sum_s v^s(p)\bar{v}^s(p) = \not{p} - m$$

The electron-side spin sum becomes a trace (setting $m_e = 0$):

$$L_e^{\mu\nu} = \sum_{\text{spins}} [\bar{v}\gamma^\mu u][\bar{v}\gamma^\nu u]^* = \text{Tr}\left[\not{p}_2 \gamma^\mu \not{p}_1 \gamma^\nu\right]$$

The muon-side spin sum (keeping $m_\mu$):

$$L_\mu^{\mu\nu} = \sum_{\text{spins}} [\bar{u}\gamma_\mu v][\bar{u}\gamma_\nu v]^* = \text{Tr}\left[(\not{p}_3 + m_\mu)\gamma^\mu (\not{p}_4 - m_\mu)\gamma^\nu\right]$$

The full spin-summed result is:

$$\sum_{\text{spins}}|\mathcal{M}|^2 = \frac{e^4}{s^2}\, L_{e\,\mu\nu}\, L_\mu^{\mu\nu}$$

4.4 Evaluating the Traces

Electron trace: Using $\text{Tr}[\gamma^\alpha\gamma^\beta\gamma^\mu\gamma^\nu] = 4(g^{\alpha\beta}g^{\mu\nu} - g^{\alpha\mu}g^{\beta\nu} + g^{\alpha\nu}g^{\beta\mu})$:

$$L_e^{\mu\nu} = \text{Tr}[\not{p}_2\gamma^\mu\not{p}_1\gamma^\nu] = 4\left[p_2^\mu p_1^\nu + p_2^\nu p_1^\mu - (p_1 \cdot p_2)g^{\mu\nu}\right]$$

Muon trace: Expanding $(\not{p}_3 + m_\mu)\gamma^\mu(\not{p}_4 - m_\mu)\gamma^\nu$, the terms with odd numbers of gamma matrices vanish in the trace. The $m_\mu^2$ term gives $-m_\mu^2 \text{Tr}[\gamma^\mu\gamma^\nu] = -4m_\mu^2 g^{\mu\nu}$. Thus:

$$L_\mu^{\mu\nu} = 4\left[p_3^\mu p_4^\nu + p_3^\nu p_4^\mu - (p_3 \cdot p_4 + m_\mu^2)g^{\mu\nu}\right]$$

Contracting the tensors: We compute $L_{e\,\mu\nu} L_\mu^{\mu\nu}$ by contracting term by term. Using$g_{\mu\nu}g^{\mu\nu} = 4$ in 4 dimensions:

\begin{align*} L_{e\,\mu\nu}L_\mu^{\mu\nu} &= 16\Big[ (p_2 \cdot p_3)(p_1 \cdot p_4) + (p_2 \cdot p_4)(p_1 \cdot p_3) \\ &\quad - (p_1 \cdot p_2)(p_3 \cdot p_4 + m_\mu^2) \\ &\quad - (p_3 \cdot p_4 + m_\mu^2)(p_1 \cdot p_2) \\ &\quad + (p_1 \cdot p_2)(p_3 \cdot p_4 + m_\mu^2) \cdot 4 / 4 \cdot \text{(from } g\text{-terms)} \Big] \end{align*}

After careful bookkeeping, the contraction simplifies to:

$$L_{e\,\mu\nu}L_\mu^{\mu\nu} = 32\left[(p_1 \cdot p_3)(p_2 \cdot p_4) + (p_1 \cdot p_4)(p_2 \cdot p_3) + m_\mu^2(p_1 \cdot p_2)\right]$$

💡Trace Contraction Strategy

The contraction of two rank-2 tensors involves 9 pairs of terms. But symmetry cuts the work: $L_e^{\mu\nu}$ and $L_\mu^{\mu\nu}$ share the same tensor structure (symmetric, traceless part plus trace), so many cross terms simplify. In practice, experienced physicists use Mandelstam variables to bypass the explicit contraction.

4.5 Result in Mandelstam Variables

The Mandelstam variables for this process are:

$$s = (p_1+p_2)^2, \quad t = (p_1-p_3)^2, \quad u = (p_1-p_4)^2$$

with the constraint $s + t + u = 2m_\mu^2$ (setting $m_e = 0$). The dot products become:

\begin{align*} p_1 \cdot p_2 &= s/2 \\ p_1 \cdot p_3 &= -t/2 + m_\mu^2/2 \quad \Rightarrow \quad -(t - m_\mu^2)/2 \\ p_1 \cdot p_4 &= -u/2 + m_\mu^2/2 \quad \Rightarrow \quad -(u - m_\mu^2)/2 \end{align*}

By momentum conservation, $p_2 \cdot p_3 = p_1 \cdot p_4$ and $p_2 \cdot p_4 = p_1 \cdot p_3$ (in the massless electron limit). Substituting into our trace result:

$$\overline{|\mathcal{M}|^2} = \frac{1}{4}\cdot\frac{e^4}{s^2}\cdot 32\left[\frac{(t-m_\mu^2)^2}{4} + \frac{(u-m_\mu^2)^2}{4} + \frac{m_\mu^2 s}{2}\right]$$

In the high-energy limit $s \gg m_\mu^2$, we set $m_\mu = 0$ everywhere and use $s + t + u = 0$:

$$\boxed{\overline{|\mathcal{M}|^2} = \frac{2e^4}{s^2}\left(t^2 + u^2\right)} \qquad (m_e = m_\mu = 0)$$

4.6 Differential Cross Section in the CM Frame

In the center-of-mass frame with all particles massless, $|\vec{p}_f| = |\vec{p}_i|$ and the master formula gives:

$$\frac{d\sigma}{d\Omega} = \frac{|\mathcal{M}|^2}{64\pi^2 s}$$

We express $t$ and $u$ in terms of the scattering angle $\theta$. With $E = \sqrt{s}/2$ for each particle:

$$t = -2E^2(1 - \cos\theta) = -\frac{s}{2}(1 - \cos\theta)$$
$$u = -2E^2(1 + \cos\theta) = -\frac{s}{2}(1 + \cos\theta)$$

Therefore:

$$t^2 + u^2 = \frac{s^2}{4}\left[(1-\cos\theta)^2 + (1+\cos\theta)^2\right] = \frac{s^2}{2}(1 + \cos^2\theta)$$

Putting it all together and writing $e^2 = 4\pi\alpha$:

$$\boxed{\frac{d\sigma}{d\Omega} = \frac{\alpha^2}{4s}(1 + \cos^2\theta)}$$

This is one of the most celebrated results in QED. The $(1 + \cos^2\theta)$ angular dependence is a direct signature of spin-1/2 particles coupling to a spin-1 mediator.

4.7 Total Cross Section

Integrating over the full solid angle $d\Omega = \sin\theta\, d\theta\, d\phi$:

$$\sigma = \int \frac{d\sigma}{d\Omega}\, d\Omega = \frac{\alpha^2}{4s}\int_0^{2\pi}d\phi \int_0^\pi (1+\cos^2\theta)\sin\theta\, d\theta$$

The $\phi$ integral gives $2\pi$. For the $\theta$ integral, let $x = \cos\theta$:

$$\int_{-1}^{1}(1 + x^2)\, dx = \left[x + \frac{x^3}{3}\right]_{-1}^{1} = \left(1 + \frac{1}{3}\right) - \left(-1 - \frac{1}{3}\right) = \frac{8}{3}$$

Therefore:

$$\boxed{\sigma(e^+e^- \to \mu^+\mu^-) = \frac{4\pi\alpha^2}{3s}}$$

💡The Muon Point Cross Section

This result, $\sigma = 4\pi\alpha^2/(3s)$, is so fundamental that it serves as the "standard candle" of particle physics. At $\sqrt{s} = 10$ GeV, it gives$\sigma \approx 0.87$ nb. All other $e^+e^-$ cross sections are compared to this as a ratio (the R-ratio), which we explore on Page 4. The $1/s$ behavior is universal for point-like fermion pair production.

Key Concepts (Page 1)

  • • e⁺e⁻ → μ⁺μ⁻ has a single s-channel diagram (lepton flavor conservation)
  • • Feynman rules give M = (-e²/s)[v̄γ¹u][ūγ⊂μv]
  • • Spin averaging: factor of 1/4 for two spin-1/2 initial particles
  • • Completeness relations convert spin sums into traces over gamma matrices
  • • Trace identities: Tr[γ𝛼γ𝛽γ𝛾γ𝛿] = 4(g𝛼𝛽g𝛾𝛿 - g𝛼𝛾g𝛽𝛿 + g𝛼𝛿g𝛽𝛾)
  • • Result: dσ/dΩ = α²(1 + cos²θ)/(4s) and σ = 4πα²/(3s)
← Cross Sections
Page 1 of 4
Chapter 4: Elementary QED Processes I
Compton Scattering →