Photon Interactions and the Impulse Approximation
Modeling photon transport in High-Purity Germanium (HPGe) detectors requires a rigorous microscopic treatment of two primary competitive processes: the photoelectric effect and incoherent (Compton) scattering. While global cross-sections are standardized, tracking individual shell vacancies and secondary electron cascades introduces unique physical challenges.
The Concept of Evaluated Cross-Sections (EPICS2023)
Experimental cross-sections measured across different global accelerator laboratories frequently exhibit non-trivial systematic discrepancies due to target impurities and detector geometries. Conversely, pure quantum electrodynamics (QED) calculations often utilize boundary approximations that underrepresented solid-state crystal impacts.
Evaluated Cross-Sections (such as the EPICS2023 database provided by the IAEA) represent an international consensus standard. Expert committees collect all available experimental data, filter out systematic outliers, enforce strict conservation laws (sum rules), and perform global mathematical optimization. The resulting values reflect the absolute highest precision macroscopic baseline currently known to physics.
Precision Benchmarks for Germanium (\(Z=32\)):
- Photoelectric Effect (\(>11.1\) keV): High precision (~1–2%). The K-edge cutoff is exceptionally well-mapped.
- Photoelectric Effect (\(<1.4\) keV): Moderate precision (~5–10%) due to solid-state crystal field perturbations affecting outer shell orbitals.
- Incoherent Scattering (\(>100\) keV): Near-perfect precision (~1%), strictly converging to the relativistic Klein-Nishina formulation for free electrons.
- Incoherent Scattering (\(1 - 20\) keV): Limited precision (~5–10%). Experimental data is scarce, and the classic atomic scattering function approximation acts as a smooth collective bound.
Compton Scattering on Bound Electrons
Standard databases like EPDL provide only the integrated incoherent cross-section (\(\sigma_{\text{incoherent}}\)) for the whole atom. For a true Monte Carlo track, this macroscopic threshold is insufficient. The engine must stochastically resolve: 1. Which specific subshell (\(1s, 2s, 2p...\)) lost the electron, in order to initialize the correct relaxation cascade. 2. The exact kinetic energy spectrum of the Compton recoil electron, which undergoes Doppler broadening due to the pre-collision orbital momentum of the bound state.
The Step-Function and Impulse Approximations
In classic simulation architectures, the Incoherent Scattering Function (ISF) approach by Hubbell is utilized:
Where \(S(x, Z)\) represents the collective atomic scattering function. While highly accurate for global energy deposition tracks, ISF treats the electron shell as a single collective cloud and collapses the recoil profile back to a free, stationary electron baseline.
To bypass this limitation, GEIM utilizes a multi-shell пороговая модель (Step-Function model). It enforces an absolute energy conservation bound: scattering cannot occur on a specific subshell unless the energy transferred to the electron exceeds its binding threshold (\(E_{\text{bind}}\)):
By implementing this binding energy cutoff directly into the boundary limits of the Klein-Nishina integral, the simulation splits the atomic cross-section into independent, shell-specific Monte Carlo channels.
At low energies, this causes shells to sequentially "freeze out" as photon energy drops below their ionization potentials, matching the true quantum behavior of independent bound states.