Quantum mechanical scattering of photon and nucleus explained

In pair production, a photon creates an electron positron pair. In the process of photons scattering in air (e.g. in lightning discharges), the most important interaction is the scattering of photons at the nuclei of atoms or molecules. The full quantum mechanical process of pair production can be described by the quadruply differential cross section given here:[1]

\begin{align} d4\sigma&=

2\alpha
Z
rm{fine
3c

2}{(2\pi)

2\hbar}|p
+||p
-| dE+
\omega3
d\Omega+d\Omega-d\Phi
|q|4

x \\ & x \left[-

2\sin
p
2\Theta
-
-
(E-|\cos\Theta
2
-)
--c|p
2-c
\left (4E
+

2q

2\right)\right.\\ &-
2\sin
p
2\Theta
+
+
(E+|\cos\Theta
2
+)
+-c|p
2-c
\left (4E
-

2q2\right)\\ &+2\hbar2\omega

2
2\sin
p
2\Theta
++p
2\sin
-
2\Theta
-
+
(E+-c|p+|\cos\Theta+)(E--c|p-|\cos\Theta-)
\\ &+2\left.
|p+||p-|\sin\Theta+\sin\Theta-\cos\Phi
(E+-c|p+|\cos\Theta+)(E--c|p-|\cos\Theta-)
2-c
\left(2E
-

2q2\right)\right].\\ \end{align}

with

\begin{align} d\Omega+&=\sin\Theta+ d\Theta+,\\ d\Omega-&=\sin\Theta- d\Theta-. \end{align}

This expression can be derived by using a quantum mechanical symmetry between pair production and Bremsstrahlung.

Z

is the atomic number,

\alphafine1/137

the fine structure constant,

\hbar

the reduced Planck constant and

c

the speed of light. The kinetic energies

Ekin,+/-

of the positron and electron relate to their total energies

E+,-

and momenta

p+,-

via

E+,-=Ekin,+/-+me

2
c
e
2
c
+,-

c2}.

Conservation of energy yields

\hbar\omega=E++E-.

The momentum

q

of the virtual photon between incident photon and nucleus is:
2-\left(\hbar
c
\begin{align} -q
-|
2+2|p
\omega\right)
+|\hbar
c

\omega\cos\Theta+

+2|p
-|\hbar
c

\omega\cos\Theta-\\ &-2|p+||p-|(\cos\Theta+\cos\Theta-+\sin\Theta+\sin\Theta-\cos\Phi), \end{align}

where the directions are given via:

\begin{align} \Theta+&=\sphericalangle(p+,k),\\ \Theta-&=\sphericalangle(p-,k),\\ \Phi&=Anglebetweentheplanes(p+,k)and(p-,k), \end{align}

where

k

is the momentum of the incident photon.

In order to analyse the relation between the photon energy

E+

and the emission angle

\Theta+

between photon and positron, Köhn and Ebert integrated [2] the quadruply differential cross section over

\Theta-

and

\Phi

. The double differential cross section is:
\begin{align} d2\sigma(E+,\omega,\Theta+)
dE+d\Omega+
6
= \sum\limits
j=1

Ij\end{align}

with
\begin{align} I
1&=2\piA
2\sin
\sqrt{(\Delta
2\Theta
+
-
} \\&\times\ln\left(\frac\right) \\&\times\left[-1-\frac{c\Delta^{(p)}_2}{p_-(E_+-cp_+\cos\Theta_+)}+\frac{p_+^2c^2\sin^2\Theta_+} {(E_+-cp_+\cos\Theta_+)^2}-\frac{2\hbar^2\omega^2p_-\Delta^{(p)}_2}{c(E_+-cp_+\cos \Theta_+)((\Delta^{(p)}_2)^2+4p_+^2p_-^2\sin^2\Theta_+)}\right], \\I_2&=\frac\ln\left(\frac\right), \\I_3&=\frac \\&\times\ln\Bigg(\Big((E_-+p_-c)(4p_+^2p_-^2\sin^2\Theta_+(E_--p_-c)+(\Delta^_1+\Delta^_2)((\Delta^_2E_-+\Delta^_1p_-c) \\&-\sqrt))\Big)\Big((E_--p_-c)(4p_+^2p_-^2\sin^2\Theta_+(-E_--p_-c) \\&+(\Delta^_1-\Delta^_2)((\Delta^_2E_-+\Delta^_1p_-c)-\sqrt))\Big)^\Bigg) \\&\times\left[\frac{c(\Delta^{(p)}_2E_-+\Delta^{(p)}_1p_-c)}{p_-(E_+-cp_+\cos\Theta_+)}\right.\\ &+\Big[((\Delta^{(p)}_2)^2+4p_+^2p_-^2\sin^2\Theta_+)(E_-^3+E_-p_-c)+p_-c(2 ((\Delta^{(p)}_1)^2-4p_+^2p_-^2\sin^2\Theta_+)E_-p_-c \\ &+\Delta^{(p)}_1\Delta^{(p)}_2(3E_-^2+p_-^2c^2))\Big]\Big[(\Delta^{(p)}_2E_-+\Delta^{(p)}_1p_-c)^2+4m^2c^4p_+^2p_-^2\sin^2\Theta_+\Big]^ \\&+\Big[-8p_+^2p_-^2m^2c^4\sin^2\Theta_+(E_+^2+E_-^2)-2\hbar^2\omega^2p_+^2\sin^2\Theta_+p_-c(\Delta^{(p)}_2E_-+\Delta^{(p)}_1p_-c) \\ &+2\hbar^2\omega^2p_- m^2c^3(\Delta^{(p)}_2E_-+\Delta^{(p)}_1p_-c)\Big]\Big[(E_+-cp_+\cos\Theta_+)((\Delta^{(p)}_2E_-+\Delta^{(p)}_1p_-c)^2+4m^2c^4p_+^2p_-^2\sin^2\Theta_+)\Big]^ \\&+\left.\frac\right], \\I_4&=\frac+\frac, \\I_5&=\frac \\&\times\left[\frac{\hbar^2\omega^2p_-^2}{E_+cp_+\cos\Theta_+} \Big[E_-[2(\Delta^{(p)}_2)^2((\Delta^{(p)}_2)^2-(\Delta^{(p)}_1)^2)+8p_+^2p_-^2\sin^2\Theta_+((\Delta^{(p)}_2)^2+(\Delta^{(p)}_1)^2)]\right.\\&+p_-c[2\Delta^{(p)}_1\Delta^{(p)}_2((\Delta^{(p)}_2)^2-(\Delta^{(p)}_1)^2)+16\Delta^{(p)}_1\Delta^{(p)}_2p_+^2p_-^2\sin^2\Theta_+]\Big]\Big[(\Delta^{(p)}_2)^2+4p_+^2p_-^2\sin^2\Theta_+\Big]^\\&+ \frac\\&-\Big[2E_+^2p_-^2\{2((\Delta^{(p)}_2)^2-(\Delta^{(p)}_1)^2)(\Delta^{(p)}_2E_-+\Delta^{(p)}_1p_-c)^2 +8p_+^2p_-^2\sin^2\Theta_+[((\Delta^{(p)}_1)^2+(\Delta^{(p)}_2)^2)(E_-^2+p_-^2c^2)\\ &+4\Delta^{(p)}_1\Delta^{(p)}_2E_-p_-c]\}\Big]\Big[(\Delta^{(p)}_2E_-+\Delta^{(p)}_1p_-c)^2+4m^2c^4p_+^2p_-^2\sin^2\Theta_+\Big]^\\&-\left.\frac\right], \\I_6&=-\frac \endand
\begin{align} A&=
3c
Z2
fine
(2\pi)2\hbar
|p+||p-|
\omega3
(p)
,\\ \Delta
1&:=-|p
2-\left(\hbar
c
-|

\omega\right) +2

\hbar
c

\omega|p+|\cos\Theta

(p)
2&:=2\hbar
c

\omega|pi|-2|p+||p-| \cos\Theta++2. \end{align}

This cross section can be applied in Monte Carlo simulations. An analysis of this expression shows that positrons are mainly emitted in the direction of the incident photon.

Notes and References

  1. Bethe, H.A., Heitler, W., 1934. On the stopping of fast particles and on the creation of positive electrons. Proc. Phys. Soc. Lond. 146, 83–112
  2. Koehn, C., Ebert, U., Angular distribution of Bremsstrahlung photons and of positrons for calculations of terrestrial gamma-ray flashes and positron beams, Atmos. Res. (2014), vol. 135-136, pp. 432-465