Ionization

Method

The $G4eIonisation$ class provides the continuous and discrete energy losses of electrons and positrons due to ionization in a material according to the approach described in Section 7.1. The value of the maximum energy transferable to a free electron $T_{max}$ is given by the following relation:

$$T_{max} = \left\{ \begin{array}{ll} E-mc^2 & {for \hspace{.... ...E-mc^2)/2 & {for \hspace{.2cm} e^- } \\ \end{array} \right .$$ (8.1)

where $mc^2$ is the electron mass. Above a given threshold energy the energy loss is simulated by the explicit production of delta rays by Möller scattering ($e^- e^-$), or Bhabha scattering ($e^+ e^-$). Below the threshold the soft electrons ejected are simulated as continuous energy loss by the incident ${e^{\pm}}$.

Continuous Energy Loss

The integration of 7.1 leads to the Berger-Seltzer formula [1]:

$$\left. \frac{dE}{dx} \right]_{T < T_{cut}} = 2 \pi r_e^2 mc... ...)} {(I/mc^2)^2}+ F^{\pm} (\tau , \tau_{up}) - \delta \right ]$$ (8.2)

with

$$\begin{array}{ll} r_e & \mbox{classical electron radius:} ... ...max}) \\ \delta & \mbox{density effect function} . \end{array}$$

In an elemental material the electron density is

$$n_{el} = Z \: n_{at} = Z \: \frac{\mathcal{N}_{av} \rho}{A} .$$

$\mathcal{N}_{av}$ is Avogadro's number, $\rho$ is the material density, and $A$ is the mass of a mole. In a compound material

$$n_{el} = \sum_i Z_i \: n_{ati} = \sum_i Z_i \: \frac{\mathcal{N}_{av} w_i \rho}{A_i} ,$$

where $w_i$ is the proportion by mass of the $i^{th}$ element, with molar mass $A_i$ .

The mean excitation energies $I$ for all elements are taken from [2].

The functions $F^{\pm}$ are given by :

$$F^+ (\tau,\tau_{up}) = \ln(\tau\tau_{up} )$$ (8.3)

$$-\frac{\tau_{up}^2}{\tau}\left[\tau + 2 \tau_{up} - \frac{3\tau_{... ...}{2} - \tau \frac{\tau_{up}^3}{3} + \frac{\tau_{up}^4 } {4} \right) y^3 \right]$$

$$F^- (\tau,\tau_{up} ) = -1 -\beta^2$$ (8.4)

$$+\ln \left [(\tau - \tau_{up}) \tau_{up} \right ] + \frac{\tau}{\... ...{up}^2}{2} + ( 2\tau +1) \ln \left (1- \frac{\tau_{up}}{\tau} \right ) \right ]$$

where $y = 1/(\gamma+1)$.

The density effect correction is calculated according to the formalism of Sternheimer [3]:
$x$ is a kinetic variable of the particle : $x = \log_{10}(\gamma \beta) = \ln(\gamma^{2} \beta^{2})/4.606$,
and $\delta(x)$ is defined by

$$\begin{array}{rll} \mbox{for} & x < x_0 : & \delta(x) = 0 \\... ...\ \mbox{for} & x > x_1 : & \delta(x) = 4.606 x - C \end{array}$$ (8.5)

where the matter-dependent constants are calculated as follows:

$$\begin{array}{lcl} h\nu_p & = & \mbox{ plasma energy of the ... ... = & 4.606(x_a - x_0)/(x_1 - x_0)^m \\ m & = & 3 . \end{array}$$ (8.6)

For condensed media

$$\begin{array}{ll} I < 100 \: \mbox{eV} & \left \{ \begin{arr... ... x_0 = 0.326 C - 1.5 & x_1 = 3 \end{array} \right . \end{array}$$

and for gaseous media

$$\begin{array}{rlll} \mbox{for} & C < 10. & x_0 = 1.6 & x_1 =... ...r} & C \geq 13.804 & x_0 = 0.326 C -2.5 & x_1 = 5 . \end{array}$$

Total Cross Section per Atom and Mean Free Path

The total cross section per atom for Möller scattering ($e^- e^-$) and Bhabha scattering ($e^+ e^-$) is obtained by integrating Eq. 7.2. In GEANT4 $T_{cut}$ is always 1 keV or larger. For delta ray energies much larger than the excitation energy of the material ($T \gg I$), the total cross section becomes [1] for Möller scattering,

$$\sigma ( Z,E,T_{cut} ) = \frac {2 \pi r_e^2 Z}{\beta^2(\gamma -1)} \times$$ (8.7)

$$\left[\frac{(\gamma-1)^2} {\gamma^2}\left(\frac{1}{2}-x\right) +\frac{1}{x}-\frac{1}{1-x}-\frac{2\gamma-1}{\gamma^2}\ln \frac{1-x}{x}\right] ,$$

and for Bhabha scattering ($e^+ e^-$),

$$\sigma ( Z,E,T_{cut} ) = \frac{ 2 \pi r_e^2 Z }{(\gamma -1)} \times$$ (8.8)

$$\left [\frac {1 }{\beta^2} \left(\frac{1}{x}-1\right) + B_1 \ln x + B_2 (1-x) - \frac {B_3 } {2} ( 1-x^2 ) +\frac{B_4}{3}(1-x^3)\right] .$$

Here

$$\begin{array}{lcllcl} \gamma & = & E/mc^2 & B_1 & = & 2-y^... ...\ y & = & 1/(\gamma + 1) & B_4 & = & (1-2y)^3 . \end{array}$$

The above formulae give the total cross section for scattering above the threshold energies

$$T_{\rm Moller}^{\rm thr} =2T_{cut} \mbox{\hspace{2cm}and\hspace{2cm}} T_{\rm Bhabha}^{\rm thr} = T_{cut} .$$ (8.9)

In a given material the mean free path is then

$$\begin{array}{lll} \lambda = (n_{at} \cdot \sigma)^{-1} & o... ...\left( \sum_i n_{ati} \cdot \sigma_i \right)^{-1} . \end{array}$$ (8.10)

Simulation of Delta-ray Production

Differential Cross Section

For $T \gg I$ the differential cross section per atom becomes [1] for Möller scattering,

$$\frac{d\sigma }{d \epsilon } = \frac {2 \pi r_e^2 Z}{\beta^2(\gamma -1)} \times$$ (8.11)

$$\left[ \frac{(\gamma -1 )^2} {\gamma^2 }+\frac{1}{\epsilon} \left... ...n}\left(\frac{1} {1- \epsilon} - \frac{2 \gamma - 1} {\gamma^2 }\right) \right]$$

and for Bhabha scattering,

$$\frac{d \sigma}{d \epsilon}=\frac{2 \pi r_e^2 Z}{(\gamma -1)... ...c{B_1}{\epsilon}+B_2 - B_3 \epsilon + B_4 \epsilon^2\right] .$$ (8.12)

Here $\epsilon = T/(E-mc^2)$. The kinematical limits of $\epsilon$ are

$$\epsilon_0 = \frac{T_{cut}}{E-mc^2} \leq \epsilon \leq \frac... ...c^2} \leq \epsilon \leq 1 \mbox{\hspace{.2cm} for $e^+ e^-$} .$$

Sampling

The delta ray energy is sampled according to methods discussed in Chapter 2. Apart from normalization, the cross section can be factorized as

$$\frac{d\sigma}{d\epsilon}=f(\epsilon) g(\epsilon) .$$ (8.13)

For $e^- e^-$ scattering

$$f(\epsilon) = \frac{1}{\epsilon^2} \frac{\epsilon_0 }{1- 2\epsilon_0}$$ (8.14)

$$g(\epsilon) = \frac{4}{9\gamma^2 - 10 \gamma + 5}\left[(\gamma -1)^2 \epsilon^2... ...a -1) \frac{\epsilon} {1- \epsilon }+ \frac{\gamma^2}{(1- \epsilon )^2 }\right]$$ (8.15)

and for $e^+ e^-$ scattering

$$f(\epsilon) = \frac{1}{\epsilon^2} \frac{\epsilon_0}{1- \epsilon_0 }$$ (8.16)

$$g(\epsilon) = \frac{B_0 -B_1 \epsilon +B_2 \epsilon^2 -B_3 \epsilon^3 +B_4 \eps... ... ^4}{B_ 0-B_1\epsilon_0 +B_2\epsilon_0^2 -B_3 \epsilon_0^3 +B_4 \epsilon_0^4} .$$ (8.17)

Here $B_0=\gamma^2/(\gamma^2-1)$ and all other quantities have been defined above.

To choose $\epsilon$, and hence the delta ray energy,

1. $\epsilon$ is sampled from $f(\epsilon)$
2. the rejection function $g(\epsilon)$ is calculated using the sampled value of $\epsilon$
3. $\epsilon$ is accepted with probability $g(\epsilon)$.

After the successful sampling of $\epsilon$, the direction of the ejected electron is generated with respect to the direction of the incident particle. The azimuthal angle $\phi$ is generated isotropically and the polar angle $\theta$ is calculated from energy-momentum conservation. This information is used to calculate the energy and momentum of both the scattered incident particle and the ejected electron, and to transform them to the global coordinate system.

Status of this document

9.10.98 created by L. Urbán.
29.07.01 revised by M.Maire.
13.12.01 minor cosmetic by M.Maire.
24.05.02 re-written by D.H. Wright.
01.12.03 revised by V. Ivanchenko.

Bibliography

1. H.Messel and D.F.Crawford. Pergamon Press,Oxford,1970.
2. ICRU Report No. 37 (1984)
3. R.M.Sternheimer. Phys.Rev. B3 (1971) 3681.