Penelope Model
Total cross section
The total cross section of the Compton scattering process is determined from an analytical parameterization. For \(\gamma\) energy \(E\) greater than 5 MeV, the usual Klein-Nishina formula is used for \(\sigma(E)\). For a more accurate parameterization is used, which takes into account atomic binding effects and Doppler broadening [eal96]:
where:
Finally,
where \(J_{i0}\) is the value of the \(p_{z}\)-distribution profile \(J_{i}(p_{z})\) for the i-th atomic shell calculated in \(p_{z}=0\). The values of \(J_{i0}\) for the different shells of the different elements are tabulated from the Hartree-Fock atomic orbitals of Ref. [eal75].
The integration of Eq.(7) is performed numerically using the 20-point Gaussian method. For this reason, the initialization of the Penelope Compton process is somewhat slower than the Low Energy Livermore process.
Sampling of the final state
The polar deflection \(\cos\theta\) is sampled from the probability density function
(see Ref. [eal01] for details on the sampling algorithm). Once the direction of the emerging photon has been set, the active electron shell \(i\) is selected with relative probability equal to \(Z_{i} \Theta (E-U_{i}) n_{i} [p_{z}^{max} (E,\theta)]\).
A random value of \(p_{z}\) is generated from the analytical Compton profile [eal75]. The energy of the emerging photon is
where
The azimuthal scattering angle \(\phi\) of the photon is sampled uniformly in the interval \((0, 2\pi)\). It is assumed that the Compton electron is emitted with energy \(E_{e} = E-E'-U_{i}\), with polar angle \(\theta_{e}\) and azimuthal angle \(\phi_{e} = \phi + \pi\), relative to the direction of the incident photon. In this case \(\cos\theta_{e}\) is given by
Since the active electron shell is known, characteristic x-rays and electrons emitted in the de-excitation of the ionized atom can also be followed. The de-excitation is simulated as described in Atomic relaxation. For further details see [eal01].
Bibliography
- eal96
D. Brusa et al. Fast sampling algorithm for the simulation of photon compton scattering. NIM A, 379():167, 1996.
- eal01(1,2)
F. Salvat et al. Penelope - a code system for monte carlo simulation of electron and photon transport. Technical Report, Workshop Proceedings Issy-les-Moulineaux, France; AEN-NEA, 5-7 November 2001.
- eal75(1,2)
F.Biggs et al. Hartree-fock compton profiles for the elements. At. Data Nucl. Data Tables, 16():201, 1975.