4 Cross Sections

To be able to compare our calculated amplitudes to experimental data, we need to a bit more work. We currently have a way to calculating the scattering amplitude for a single phase space point but we need to integrate this over (regions of) phase space to be able to make real predictions. For a scattering process this is referred to cross section and for a decay process it is the decay width. These are obtained by integrating the matrix element squared over all possible momenta and multiply with the correct flux.

4.1 Phase space

We must integrate over all possible four-momenta of each particle involved. While doing so we must ensure that all particles are on-shell, i.e. pi2=mi2p_{i}^{2}=m_{i}^{2}, that their energies are positive pi0>0p_{i}^{0}>0 (this is still Lorentz invariant!), and that momentum is conserved ∑piin=∑piout\sum p_{i}^{\rm in}=\sum p_{i}^{\rm out}.

Let us begin with the phase space of a single particle. The requirements set out above can be implemented with a δ\delta function and a Θ\Theta function

d⁢ϕi=d4⁢pi(2⁢π)4⁢(2⁢π)⁢δ⁢(pi2−mi2)⁢Θ⁢(pi0)=d4⁢pi(2⁢π)4⁢(2⁢π)⁢δ⁢((pi0)2−p→i2−mi2)⁢Θ⁢(pi0).\displaystyle{\rm d}\phi_{i}=\frac{{\rm d}^{4}p_{i}}{(2\pi)^{4}}(2\pi)\delta(p% _{i}^{2}-m_{i}^{2})\Theta(p_{i}^{0})=\frac{{\rm d}^{4}p_{i}}{(2\pi)^{4}}(2\pi)% \delta\Big{(}(p_{i}^{0})^{2}-\vec{p}_{i}^{2}-m_{i}^{2}\Big{)}\Theta(p_{i}^{0})\,. (103)

The δ\delta function can be used to solve the pi0p_{i}^{0} integral – almost. We first need to use the following identity (convince yourself that this is true)

δ⁢(x2−α2)=1|2⁢α|⁢[δ⁢(x+α)+δ⁢(x−α)].\displaystyle\delta(x^{2}-\alpha^{2})=\frac{1}{|2\alpha|}\Big{[}\delta(x+% \alpha)+\delta(x-\alpha)\Big{]}\,. (104)

Therefore,

d⁢ϕi\displaystyle{\rm d}\phi_{i} =d4⁢pi(2⁢π)4⁢|2⁢pi0|⁢(2⁢π)⁢[δ⁢(pi0−p→i2+mi2)+δ⁢(pi0+p→i2+mi2)]⁢Θ⁢(pi0)\displaystyle=\frac{{\rm d}^{4}p_{i}}{(2\pi)^{4}|2p_{i}^{0}|}(2\pi)\bigg{[}% \delta\Big{(}p_{i}^{0}-\sqrt{\vec{p}_{i}^{2}+m_{i}^{2}}\Big{)}+\delta\Big{(}p_% {i}^{0}+\sqrt{\vec{p}_{i}^{2}+m_{i}^{2}}\Big{)}\bigg{]}\Theta(p_{i}^{0}) (105)
=d4⁢pi(2⁢π)4⁢2⁢pi0⁢(2⁢π)⁢δ⁢(pi0−p→i2+mi2)=d3⁢pi(2⁢π)3⁢2⁢pi0,\displaystyle=\frac{{\rm d}^{4}p_{i}}{(2\pi)^{4}2p_{i}^{0}}(2\pi)\delta\Big{(}% p_{i}^{0}-\sqrt{\vec{p}_{i}^{2}+m_{i}^{2}}\Big{)}=\frac{{\rm d}^{3}p_{i}}{(2% \pi)^{3}2p_{i}^{0}}\,, (106)

where we have used the Θ\Theta function requiring pi0>0p_{i}^{0}>0. The condition that the particle must be on-shell, i.e. pi0=p→i2+mi2p_{i}^{0}=\sqrt{\vec{p}_{i}^{2}+m_{i}^{2}} is now implicit and for the sake of readability keep pi0p_{i}^{0} in the denominator. This object is manifestly Lorentz invariant even if it does not look like it is. Therefore it is called the Lorentz-invariant phase space (LIPS). The total phase space is now

d⁢Φ=(∏id⁢ϕi)⁢(2⁢π)4⁢δ(4)⁢(pin−∑ipi)\displaystyle{\rm d}\Phi=\Big{(}\prod_{i}{\rm d}\phi_{i}\Big{)}(2\pi)^{4}% \delta^{(4)}\Big{(}p_{\rm in}-\sum_{i}p_{i}\Big{)} (107)

To calculate a cross section, we need to add a flux factor that accounts for the amount of incoming particles. For scattering, this is usually given as (see Appendix B for a short derivation)

ℱ=1(2⁢E1)⁢(2⁢E2)⁢|v→1−v→2|=14⁢(p1⋅p2)2−m12⁢m22.\displaystyle\mathcal{F}=\frac{1}{(2E_{1})(2E_{2})|\vec{v}_{1}-\vec{v}_{2}|}=% \frac{1}{4\sqrt{(p_{1}\cdot p_{2})^{2}-m_{1}^{2}m_{2}^{2}}}\,. (108)

Here, EiE_{i} and v→i\vec{v}_{i} are the energies and velocities of the incoming particles. The second version is equivalent but a bit easier to work with. Note that despite only depending on Lorentz invariant quantities like p1⋅p2p_{1}\cdot p_{2} and mi2m_{i}^{2}, ℱ\mathcal{F} is only invariant for boosts along the beam axis. For boosts along any other axis, it is not invariant which fits well with our intuition of cross-sectional areas.

The cross section is therefore

σ=∫ℱ×dΦ×∑|ℳ|2,\displaystyle\sigma=\int\mathcal{F}\times{\rm d}\Phi\times\sum|\mathcal{M}|^{2% }\,, (109)

where the sum averages over initial states and sums over final states. For a decay width, the flux is instead

ℱ=12⁢M,\displaystyle\mathcal{F}=\frac{1}{2M}\,, (110)

so that

Γ=12⁢M×d⁢Φ×∑|ℳ|2.\displaystyle\Gamma=\frac{1}{2M}\times{\rm d}\Phi\times\sum|\mathcal{M}|^{2}\,. (111)

It is useful to remember the phase space for a 1→21\to 2 or 2→22\to 2 process. With incoming momentum PP and outgoing momenta p1p_{1} and p2p_{2},

d⁢Φ=d⁢ϕ1⁢d⁢ϕ2⁢(2⁢π)4⁢δ(4)⁢(P−p1−p2)=d3⁢p1(2⁢π)3⁢2⁢p10⁢d3⁢p2(2⁢π)3⁢2⁢p20⁢(2⁢π)4⁢δ(4)⁢(P−p1−p2).\displaystyle{\rm d}\Phi={\rm d}\phi_{1}{\rm d}\phi_{2}\ (2\pi)^{4}\delta^{(4)% }(P-p_{1}-p_{2})=\frac{{\rm d}^{3}p_{1}}{(2\pi)^{3}2p_{1}^{0}}\frac{{\rm d}^{3% }p_{2}}{(2\pi)^{3}2p_{2}^{0}}\ (2\pi)^{4}\delta^{(4)}(P-p_{1}-p_{2})\,. (112)

We can use three of the four δ\delta functions to perform the integration over p2p_{2}

d⁢Φ=d3⁢p1(2⁢π)3⁢2⁢p10⁢12⁢p20⁢(2⁢π)⁢δ⁢(P0−p10−p20).\displaystyle{\rm d}\Phi=\frac{{\rm d}^{3}p_{1}}{(2\pi)^{3}2p_{1}^{0}}\frac{1}% {2p_{2}^{0}}\ (2\pi)\delta(P^{0}-p_{1}^{0}-p_{2}^{0})\,. (113)

We require here that pi0=p→i2+mi2p_{i}^{0}=\sqrt{\vec{p}_{i}^{2}+m_{i}^{2}} and that p→2=−p→1\vec{p}_{2}=-\vec{p}_{1} in the rest frame of PP (this is often called the centre-of-mass frame). Next, we use spherical coordinates to write d3⁢p1=|p→1|2⁢d⁢|p→1|⁢d⁢Ω{\rm d}^{3}p_{1}=|\vec{p}_{1}|^{2}{\rm d}|\vec{p}_{1}|\ {\rm d}\Omega.

d⁢Φ=|p→1|2⁢d⁢|p→1|⁢d⁢Ω(2⁢π)3⁢2⁢p10⁢12⁢p20⁢(2⁢π)⁢δ⁢(P0−p10−p20).\displaystyle{\rm d}\Phi=\frac{|\vec{p}_{1}|^{2}{\rm d}|\vec{p}_{1}|\ {\rm d}% \Omega}{(2\pi)^{3}2p_{1}^{0}}\frac{1}{2p_{2}^{0}}\ (2\pi)\delta(P^{0}-p_{1}^{0% }-p_{2}^{0})\,. (114)

After an annoying but not particularly difficult calculation, we arrive at

d⁢Φ=116⁢π2⁢|p→1|P0⁢d⁢Ω.\displaystyle{\rm d}\Phi=\frac{1}{16\pi^{2}}\frac{|\vec{p}_{1}|}{P^{0}}{\rm d}% \Omega\,. (115)

The momentum |p→i||\vec{p}_{i}| is

|p→1|=|p→2|=((P0)2−(m1+m2)2)⁢((P0)2−(m1−m2)2)2⁢P0.\displaystyle|\vec{p}_{1}|=|\vec{p}_{2}|=\frac{\sqrt{\big{(}(P^{0})^{2}-(m_{1}% +m_{2})^{2}\big{)}\big{(}(P^{0})^{2}-(m_{1}-m_{2})^{2}\big{)}}}{2P^{0}}\,. (116)

4.2 Example: electron-muon scattering

We can continue our discussion from Section 3.4 and calculate the cross section. Using the above result for the phase space and P0=sP^{0}=\sqrt{s}, we have

d⁢σd⁢Ω=ℱ⁢116⁢π2⁢|p→1|P0⁢2⁢e4t2⁢((s−m2−M2)2+(u−m2−M2)2+2⁢t⁢(m2+M2)).\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}=\mathcal{F}\frac{1}{16\pi^{2}% }\frac{|\vec{p}_{1}|}{P^{0}}\frac{2e^{4}}{t^{2}}\Big{(}(s-m^{2}-M^{2})^{2}+(u-% m^{2}-M^{2})^{2}+2t\ (m^{2}+M^{2})\Big{)}\,. (117)

For simplicity we will set m=M=0m=M=0. Then we can write

t=−s⁢sin2⁡θ2,u=−s⁢cos2⁡θ2.\displaystyle t=-s\ \sin^{2}\frac{\theta}{2}\,,\qquad u=-s\ \cos^{2}\frac{% \theta}{2}\,. (118)

And therefore, with α=e2/4⁢π\alpha=e^{2}/4\pi

d⁢σd⁢Ω=α22⁢s⁢1+cos4⁡(θ/2)sin4⁡(θ/2).\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}=\frac{\alpha^{2}}{2s}\frac{1+% \cos^{4}(\theta/2)}{\sin^{4}(\theta/2)}\,. (119)

4.3 Example: e+⁢e−e^{+}e^{-} annihilation

The calculation we have just performed is very similar to the one needed for

e+⁢(p1)⁢e−⁢(p2)→μ+⁢(p3)⁢μ−⁢(p4).\displaystyle e^{+}(p_{1})e^{-}(p_{2})\to\mu^{+}(p_{3})\mu^{-}(p_{4})\,. (120)

Although this now involves anti-particles, there is still one single diagram at leading-order and the trace algebra is very similar. Indeed we can re-interpret the incoming e+e^{+} as an outgoing e−e^{-} with momentum −p1-p_{1} and the outgoing μ+\mu^{+} as an incoming μ−\mu^{-} with momentum −p3-p_{3}. Then we find

|ℳ⁢(e+⁢(p1)⁢e−⁢(p2)→μ+⁢(p3)⁢μ−⁢(p4))|2=|ℳ⁢(e−⁢(p2)⁢μ−⁢(−p3)→e−⁢(−p1)⁢μ−⁢(p4))|2\displaystyle\Big{|}\mathcal{M}(e^{+}(p_{1})e^{-}(p_{2})\to\mu^{+}(p_{3})\mu^{% -}(p_{4}))\Big{|}^{2}=\Big{|}\mathcal{M}(e^{-}(p_{2})\mu^{-}(-p_{3})\to e^{-}(% -p_{1})\mu^{-}(p_{4}))\Big{|}^{2} (121)

This is an example of crossing symmetry. Note in general that there is an additional minus sign for each fermion which swaps from the initial to final state or vice versa. This is because, for example,

∑s1us1⁢(p1)⁢u¯s1⁢(p1)=γ⋅p1+m⟶∑s1vs1⁢(−p1)⁢v¯s1⁢(−p1)=−γ⋅p1−m=−(γ⋅p1+m).\displaystyle\sum_{s_{1}}u_{s_{1}}(p_{1})\bar{u}_{s_{1}}(p_{1})=\gamma\cdot p_% {1}+m\ \longrightarrow\ \sum_{s_{1}}v_{s_{1}}(-p_{1})\bar{v}_{s_{1}}(-p_{1})=-% \gamma\cdot p_{1}-m=-(\gamma\cdot p_{1}+m)\,. (122)

In this case there are two minus signs whose combined effect gives just one.

We can therefore recycle our old calculation and write for m=0m=0 and M>0M>0

∑si|ℳ|2\displaystyle\sum_{s_{i}}|\mathcal{M}|^{2} =8⁢e4s2⁢(2⁢M4−4⁢M2⁢t+s2+2⁢s⁢t+2⁢t2)=4⁢q4⁢(2−(1−cos2⁡θ)⁢β2),\displaystyle=\frac{8e^{4}}{s^{2}}\Big{(}2M^{4}-4M^{2}t+s^{2}+2st+2t^{2}\Big{)% }=4q^{4}\Big{(}2-(1-\cos^{2}\theta)\beta^{2}\Big{)}\,, (123)

where we have rewrite the Mandelstam variable tt as

t=(p1−p3)2=s2⁢(−1+β22+β⁢cos⁡θ)withβ=1−4⁢M2s\displaystyle t=(p_{1}-p_{3})^{2}=\frac{s}{2}\Big{(}-\frac{1+\beta^{2}}{2}+% \beta\ \cos\theta\Big{)}\quad\text{with}\quad\beta=\sqrt{1-\frac{4M^{2}}{s}} (124)

The cross section is

d⁢σd⁢Ω=α24⁢s⁢β⁢(2−(1−cos2⁡θ)⁢β2).\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}=\frac{\alpha^{2}}{4s}\beta% \Big{(}2-(1-\cos^{2}\theta)\beta^{2}\Big{)}\,. (125)

The cross section is

d⁢σd⁢Ω=α24⁢s⁢β⁢(2−(1−cos2⁡θ)⁢β2).\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}=\frac{\alpha^{2}}{4s}\beta% \Big{(}2-(1-\cos^{2}\theta)\beta^{2}\Big{)}\,. (126)

We can now convert the above result to a total cross section by performing the integral over the solid angle. This gives

σ⁢(e+⁢e−→μ+⁢μ−)=4⁢π⁢α23⁢s.\displaystyle\sigma(e^{+}e^{-}\to\mu^{+}\mu^{-})=\frac{4\pi\alpha^{2}}{3s}\,. (127)

Now, when an electron and positron annihilate, other fermions may be produced. If these are quarks, they are then observed in the detector as hadrons. The same calculation gives

σ⁢(e+⁢e−→hadrons)=4⁢π⁢α23⁢s⁢Nc⁢∑i=1nfQi2,\displaystyle\sigma(e^{+}e^{-}\to\ {\rm hadrons})=\frac{4\pi\alpha^{2}}{3s}N_{% c}\sum_{i=1}^{n_{f}}Q_{i}^{2}\,, (128)

plus higher-order corrections, where there are NcN_{c} colours in each of the nfn_{f} massless flavours of quarks with charge QiQ_{i}. Therefore the ratio

R=σ⁢(e+⁢e−→μ+⁢μ−)σ⁢(e+⁢e−→hadrons)R=\frac{\sigma(e^{+}e^{-}\to\mu^{+}\mu^{-})}{\sigma(e^{+}e^{-}\to\ {\rm hadrons% })} (129)

has been used to measure the number of colours to be Nc=3N_{c}=3.

4.4 Example: Compton scattering

Let us calculate another process in QED, namely Compton scattering

e⁢(p1)⁢γ⁢(p2)→e⁢(p3)⁢γ⁢(p4).\displaystyle e(p_{1})\gamma(p_{2})\to e(p_{3})\gamma(p_{4})\,. (130)

We can follow the same recipe as above and calculate the two diagrams

ℳ\displaystyle\mathcal{M} =

A diagram with an electron and photon coming in (momenta p1 and p2) and meeting. The electron continues horizontally with momentum p1+p2 and then emits a photon with momentum p4 before leaving with momentum p3.
+

An electron comes in with momentum p1 and emits a photon with momentum p4 going out. Its line goes straight down with momentum p1-p4 before interacting with an incoming photon of momentum p2. It leaves with momentum p3
\displaystyle=\begin{gathered}\includegraphics{bmlimages/notes-12.svg}% \bml@image@depth{12}\bmlDescription{ A diagram with an electron and photon coming in (momenta p1 and p2) and % meeting. The electron continues horizontally with momentum p1+p2 and then % emits a photon with momentum p4 before leaving with momentum p3. }\end{gathered}\quad+\quad\begin{gathered}\includegraphics{bmlimages/notes-13.% svg}\bml@image@depth{13}\bmlDescription{ An electron comes in with momentum p1 and emits a photon with momentum p4 % going out. Its line goes straight down with momentum p1-p4 before interacting % with an incoming photon of momentum p2. It leaves with momentum p3 }\end{gathered}
(133)
=−i⁢e2⁢ϵ∗μ⁢(p4)⁢ϵν⁢(p2)⁢u¯⁢(p3)⁢(γμ⁢γ⋅(p1+p2)+m(p1+p2)2−m2⁢γν+γν⁢γ⋅(p1−p4)+m(p1−p4)2−m2⁢γμ)⁢u⁢(p1).\displaystyle=-ie^{2}\epsilon^{*\mu}(p_{4})\,\epsilon^{\nu}(p_{2})\,\,\bar{u}(% p_{3})\Bigg{(}\gamma_{\mu}\frac{\gamma\cdot(p_{1}+p_{2})+m}{(p_{1}+p_{2})^{2}-% m^{2}}\gamma_{\nu}+\gamma_{\nu}\frac{\gamma\cdot(p_{1}-p_{4})+m}{(p_{1}-p_{4})% ^{2}-m^{2}}\gamma_{\mu}\Bigg{)}u(p_{1})\,. (134)

You can check explicitly that this fulfils the Ward identity by replacing ϵ⁢(k)\epsilon(k) with kk (see tutorial sheet). We can square the amplitude, summing over fermion spins and photon polarisations and find

∑|ℳ|2=8⁢e4⁢(p1⋅p2p1⋅p4+p1⋅p4p1⋅p2+2⁢m2⁢(1p1⋅p2−1p1⋅p4)+m4⁢(1p1⋅p2−1p1⋅p4)2).\displaystyle\sum|\mathcal{M}|^{2}=8e^{4}\Bigg{(}\frac{p_{1}\cdot p_{2}}{p_{1}% \cdot p_{4}}+\frac{p_{1}\cdot p_{4}}{p_{1}\cdot p_{2}}+2m^{2}\bigg{(}\frac{1}{% p_{1}\cdot p_{2}}-\frac{1}{p_{1}\cdot p_{4}}\bigg{)}+m^{4}\bigg{(}\frac{1}{p_{% 1}\cdot p_{2}}-\frac{1}{p_{1}\cdot p_{4}}\bigg{)}^{2}\Bigg{)}\,. (135)

Here we have used the identities

γμ⁢γμ=4,γμ⁢γν⁢γμ=−2⁢γν\displaystyle\gamma_{\mu}\gamma^{\mu}=4\,,\qquad\gamma_{\mu}\gamma^{\nu}\gamma% ^{\mu}=-2\gamma^{\nu} (136)

from the tutorial to simplify the algebra before completing the trace. We have also used the rules

p1⋅p3=m2+p1⋅p2−p1⋅p4,\displaystyle p_{1}\cdot p_{3}=m^{2}+p_{1}\cdot p_{2}-p_{1}\cdot p_{4}\,, p2⋅p3=p1⋅p4,\displaystyle p_{2}\cdot p_{3}=p_{1}\cdot p_{4}\,, (137a)
p2⋅p4=p1⋅p2−p1⋅p4,\displaystyle p_{2}\cdot p_{4}=p_{1}\cdot p_{2}-p_{1}\cdot p_{4}\,, p3⋅p4=p1⋅p2,\displaystyle p_{3}\cdot p_{4}=p_{1}\cdot p_{2}\,, (137b)

that follow from momentum conservation. For completeness, we also give the amplitude squared for e−⁢(p1)⁢e+⁢(p2)→γ⁢(p3)⁢γ⁢(p4)e^{-}(p_{1})e^{+}(p_{2})\to\gamma(p_{3})\gamma(p_{4}) as

ℳ\displaystyle\mathcal{M} =

An electron comes in with momentum p1 and emits a photon with momentum p3 going out. The electron line continues downwards where it joins with a positron coming in (momentum p2) that has emitted a photon with momentum p4.
+

An electron comes in with momentum p1 and emits a photon with momentum p4 going out. The electron line continues downwards where it joins with a positron coming in (momentum p2) that has emitted a photon with momentum p3. The outgoing photons are visibly swapped.
\displaystyle=\begin{gathered}\includegraphics{bmlimages/notes-14.svg}% \bml@image@depth{14}\bmlDescription{ An electron comes in with momentum p1 and emits a photon with momentum p3 % going out. The electron line continues downwards where it joins with a % positron coming in (momentum p2) that has emitted a photon with momentum p4. }\end{gathered}\quad+\quad\begin{gathered}\includegraphics{bmlimages/notes-15.% svg}\bml@image@depth{15}\bmlDescription{ An electron comes in with momentum p1 and emits a photon with momentum p4 % going out. The electron line continues downwards where it joins with a % positron coming in (momentum p2) that has emitted a photon with momentum p3. % The outgoing photons are visibly swapped. }\end{gathered}
(140)
=−i⁢e2⁢ϵ∗μ⁢(p3)⁢ϵ∗ν⁢(p4)⁢v¯⁢(p2)⁢(γν⁢γ⋅(p1−p3)+m(p1−p3)2−m2⁢γμ⁢γμ⁢γ⋅(p1−p4)+m(p1−p4)2−m2⁢γν)⁢u⁢(p1).\displaystyle=-ie^{2}\epsilon^{*\mu}(p_{3})\,\epsilon^{*\nu}(p_{4})\,\,\bar{v}% (p_{2})\Bigg{(}\gamma_{\nu}\frac{\gamma\cdot(p_{1}-p_{3})+m}{(p_{1}-p_{3})^{2}% -m^{2}}\gamma_{\mu}\gamma_{\mu}\frac{\gamma\cdot(p_{1}-p_{4})+m}{(p_{1}-p_{4})% ^{2}-m^{2}}\gamma_{\nu}\Bigg{)}u(p_{1})\,. (141)

And for |ℳ|2|\mathcal{M}|^{2}, we have

∑|ℳ|2=8⁢e4⁢(p1⋅p3p1⋅p4+p1⋅p4p1⋅p3+2⁢m2⁢(1p1⋅p3+1p1⋅p4)−m4⁢(1p1⋅p3+1p1⋅p4)2).\displaystyle\sum|\mathcal{M}|^{2}=8e^{4}\Bigg{(}\frac{p_{1}\cdot p_{3}}{p_{1}% \cdot p_{4}}+\frac{p_{1}\cdot p_{4}}{p_{1}\cdot p_{3}}+2m^{2}\bigg{(}\frac{1}{% p_{1}\cdot p_{3}}+\frac{1}{p_{1}\cdot p_{4}}\bigg{)}-m^{4}\bigg{(}\frac{1}{p_{% 1}\cdot p_{3}}+\frac{1}{p_{1}\cdot p_{4}}\bigg{)}^{2}\Bigg{)}\,. (142)

This result could also have been obtained by setting p3↔−p2p_{3}\leftrightarrow-p_{2} and swapping the sign in (134).

Working in the rest frame of the incoming electron, we have

p1=(m,0,0,0),p2=(ω,0,0,ω),p4=(ω′,⋯,ω′⁢cos⁡θ).\displaystyle p_{1}=(m,0,0,0)\,,\quad p_{2}=(\omega,0,0,\omega)\,,\quad p_{4}=% (\omega^{\prime},\cdots,\omega^{\prime}\cos\theta)\,. (143)

We can evaluate the two scalar products we need as p1⋅p2=m⁢ωp_{1}\cdot p_{2}=m\omega and p1⋅p4=m⁢ω′p_{1}\cdot p_{4}=m\omega^{\prime}. Using (137b) which arises from m2=p32=(p1+p2−p4)2m^{2}=p_{3}^{2}=(p_{1}+p_{2}-p_{4})^{2}, we can write

ω⁢ω′⁢(1−cos⁡θ)=p2⋅p4=p1⋅p2−p1⋅p4=m⁢ω−m⁢ω′.\displaystyle\omega\,\omega^{\prime}(1-\cos\theta)=p_{2}\cdot p_{4}=p_{1}\cdot p% _{2}-p_{1}\cdot p_{4}=m\omega-m\omega^{\prime}\,. (144)

And therefore find for ω′\omega^{\prime} as a function of θ\theta

ω′=ω1+(ω/m)⁢(1−cos⁡θ).\displaystyle\omega^{\prime}=\frac{\omega}{1+(\omega/m)(1-\cos\theta)}\,. (145)

With this, we can write |ℳ|2|\mathcal{M}|^{2} very simple as

∑|ℳ|2=8⁢e4⁢(ωω′+ω′ω−sin2⁡θ).\displaystyle\sum|\mathcal{M}|^{2}=8e^{4}\bigg{(}\frac{\omega}{\omega^{\prime}% }+\frac{\omega^{\prime}}{\omega}-\sin^{2}\theta\bigg{)}\,. (146)

Note that this still depends on mm via the relation between ω′\omega^{\prime} and θ\theta. Adding the flux factor and phase space, we arrive at

d⁢σd⁢Ω=α22⁢m2⁢(ω′ω)2⁢(ωω′+ω′ω−sin2⁡θ).\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}=\frac{\alpha^{2}}{2m^{2}}% \bigg{(}\frac{\omega^{\prime}}{\omega}\bigg{)}^{2}\bigg{(}\frac{\omega}{\omega% ^{\prime}}+\frac{\omega^{\prime}}{\omega}-\sin^{2}\theta\bigg{)}\,. (147)

Let us think about this result means by considering various limits

non-relativistic

where ω≪m\omega\ll m where ω′≈ω\omega^{\prime}\approx\omega and

d⁢σd⁢Ω|ω≪m=α22⁢m2⁢(2−sin2⁡θ)=α22⁢m2⁢(1+cos2⁡θ).\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}\Bigg{|}_{\omega\ll m}=\frac{% \alpha^{2}}{2m^{2}}\big{(}2-\sin^{2}\theta\big{)}=\frac{\alpha^{2}}{2m^{2}}% \big{(}1+\cos^{2}\theta\big{)}\,. (148)

This is the Thomson cross section for the scattering of classical electromagnetic radiation by a free electron, which is symmetrical in the scattering angle, i.e. the photon is just as likely to scatter backwards as forwards.

high-energy

where ω≫m\omega\gg m, we have ω′≈m/(1−cos⁡θ)\omega^{\prime}\approx m/(1-\cos\theta) and

d⁢σd⁢Ω|ω≫m=α22⁢m2⁢ω⁢11−cos⁡θ\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}\Bigg{|}_{\omega\gg m}=\frac{% \alpha^{2}}{2m^{2}\omega}\frac{1}{1-\cos\theta} (149)

and the cross section is strongly peaked for small angles. This leads to a logarithmic enhancement when you perform the angular integration. These collinear logarithms arise whenever massless particles are emitted; this will be discussed in more detail in the phenomenology course.

small angle at high energy

Since ω>ω′\omega>\omega^{\prime}, we have in the high-energy limit 1−cos⁡θ>m/ω1-\cos\theta>m/\omega. Therefore, if θ≪1\theta\ll 1, i.e. small angle scattering,

d⁢σd⁢Ω|ω≫m,θ≪1=α2m2.\displaystyle\frac{{\rm d}\sigma}{{\rm d}\Omega}\Bigg{|}_{\omega\gg m,\theta% \ll 1}=\frac{\alpha^{2}}{m^{2}}\,. (150)

The forward (small scattering angle) Compton scattering cross section is then a valuable method to measure the QED coupling α\alpha.