Radiation Interaction with Matter & Detection Systems

1. Linear Stopping Power & The Relativistic Bethe-Bloch Formula

As a fast heavy charged particle (proton, alpha, fission fragment) traverses matter, it interacts primarily via the long-range Coulomb force with atomic electrons, transferring energy through excitation and ionization. The average linear rate of energy loss per unit path length, $-dE/dx$, is governed by the quantum mechanical Bethe-Bloch formula:

where:

• $z e$ and $v = \beta c$ are the charge and velocity of the incident projectile ($\gamma = 1 / \sqrt{1 - \beta^2}$).
• $Z, A, \rho$ are the atomic number, mass number, and mass density of the absorbing medium.
• $r_e = \frac{e^2}{4\pi\varepsilon_0 m_e c^2} \approx 2.818\text{ fm}$ is the classical electron radius.
• $I \approx 10 Z\text{ eV}$ is the mean excitation energy of the absorber atoms.
• $T_{\text{max}} = \frac{2 m_e c^2 \beta^2 \gamma^2}{1 + 2\gamma (m_e / M) + (m_e / M)^2} \approx 2 m_e c^2 \beta^2 \gamma^2$ is the maximum kinetic energy transferable to an atomic electron in a single head-on collision.
• $\delta(\beta)$ is the high-velocity density effect correction (dielectric polarization of the medium).
• $C/Z$ is the inner-shell correction (for projectile velocities comparable to inner orbital electron speeds).

2. Physical Regimes & The Bragg Peak

The Bethe-Bloch formula reveals essential physical properties of charged particle stopping:

1. Low Velocity Regime ($v \ll c$): Stopping power varies inversely with the square of velocity:
   $$-\frac{dE}{dx} \propto \frac{z^2}{v^2} \propto \frac{z^2}{E}$$
   As the charged particle penetrates matter and slows down, its rate of energy loss does not decrease—it increases dramatically, reaching a sharp maximum immediately prior to coming to rest. This pronounced ionization spike at the end of the particle range is the Bragg Peak.
2. Hadron Cancer Therapy: Protons ($150 - 250\text{ MeV}$) and Carbon ions ($400\text{ MeV/u}$) deposit minimal radiation dose in healthy surface tissue, concentrating their massive destructive Bragg peak precisely within deep-seated tumors.
3. Minimum Ionizing Particles (MIP): At $\beta\gamma \approx 3 - 4$ ($\beta \approx 0.96$), $-dE/dx$ reaches a broad global minimum of approximately $1 - 2\text{ MeV}\cdot\text{cm}^2/\text{g}$ for all singly charged particles ($z=1$).

3. Particle Range & Range Straggling

The total mean path length or range $R$ of a charged particle of initial kinetic energy $E_0$ is evaluated under the Continuous Slowing Down Approximation (CSDA):

Because energy loss arises from a large number of discrete statistical collisions, individual particle ranges fluctuate around the mean $R_0$ with a Gaussian distribution characterized by the range straggling parameter $\sigma_R / R_0 \approx 1 - 2\%$.

1. Collisional Energy Loss of Fast Electrons

Electrons ($e^-$) and positrons ($e^+$) differ fundamentally from heavy ions due to their minuscule mass ($m_e$ equals target electron mass) and quantum mechanical indistinguishability. The collisional stopping power is given by the relativistic Bethe-Møller formula:

Because $M_{\text{proj}} = M_{\text{target}}$, electrons undergo catastrophic wide-angle Coulomb scatterings, following tortuous, zigzag trajectories through the absorber. Consequently, the actual penetration depth (projected range) is typically $30 - 50\%$ shorter than the total integrated path length.

2. Radiative Energy Loss: Bremsstrahlung

When a fast electron experiences intense transverse acceleration in the strong electric field of an atomic nucleus, classical and quantum electrodynamics dictate that it radiates electromagnetic energy as Bremsstrahlung ("braking radiation"):

where $X_0$ is the radiation length of the absorbing material (in $\text{g/cm}^2$):

Bremsstrahlung losses scale proportional to the square of the target atomic number $Z^2$ and linearly with electron energy $E$.

3. The Critical Energy ($E_c$)

The total stopping power of electrons is the sum of collisional and radiative components:

The ratio of radiative to collisional losses is parameterized empirically as:

The energy at which radiative and collisional losses are exactly equal is the Critical Energy $E_c$:

• In Lead ($Z = 82$): $E_c \approx 9.5\text{ MeV}$. For $E > 10\text{ MeV}$, Bremsstrahlung completely dominates.
• In Air / Water ($Z_{\text{eff}} \approx 7.5$): $E_c \approx 85\text{ MeV}$. Collisional ionization dominates up to very high energies.
• Shielding Precaution: Energetic beta sources ($^{90}\text{Sr}$-$^{90}\text{Y}$) must never be shielded directly with Lead, as intense secondary Bremsstrahlung X-rays would be generated. Instead, low-$Z$ materials (Lucite, aluminum) are employed first to slow the electrons with minimal radiation, followed by high-$Z$ lead to attenuate residual gammas.

1. Neutron Interactions and Kinematics of Elastic Scattering

Because the neutron carries zero electric charge, it experiences no Coulomb forces and interacts exclusively with atomic nuclei via the strong force (or magnetic dipole interaction). In an elastic collision of a neutron of mass $m \approx 1$ and laboratory kinetic energy $E$ with a target nucleus of mass number $A$ initially at rest, the ratio of scattered neutron energy $E'$ to initial energy $E$ as a function of the center-of-mass scattering angle $\theta_{\text{cm}}$ is:

where the collision collision parameter $\alpha$ is:

The minimum energy after a single head-on collision ($\theta_{\text{cm}} = 180^\circ$) is $E'_{\text{min}} = \alpha E$:

• For Hydrogen ($A = 1$): $\alpha = 0$. In a head-on collision with a proton, the neutron can transfer $100\%$ of its kinetic energy in a single collision ($E'_{\text{min}} = 0$). Light water ($H_2O$) is therefore an extraordinarily compact moderator.
• For Carbon ($A = 12$): $\alpha = (11/13)^2 \approx 0.716$. The maximum energy loss per collision is only $28.4\%$.
• For Uranium ($A = 238$): $\alpha \approx 0.983$. The neutron loses at most $1.7\%$ of its energy per collision, making heavy nuclei useless as moderators.

2. Average Logarithmic Energy Decrement ($\xi$)

Assuming isotropic scattering in the center-of-mass frame (valid for s-wave scattering at $E_n < 10\text{ MeV}$), the probability distribution of $E'$ is uniform between $\alpha E$ and $E$. The mean loss of the logarithm of energy per collision, denoted $\xi$, is an invariant constant independent of initial neutron energy:

For $A > 1$, a highly accurate Taylor series approximation is $\xi \approx \frac{2}{A + 2/3}$.

• For Hydrogen ($A=1$): $\xi = 1.000$.
• For Deuterium ($A=2$): $\xi = 0.725$.
• For Carbon ($A=12$): $\xi \approx 0.158$.

The average number of collisions $N$ required to slow down a fast fission neutron from $E_0 = 2\text{ MeV}$ to thermal energy $E_{\text{th}} = 0.025\text{ eV}$ is:

In Hydrogen, $N \approx 18$ collisions; in Deuterium, $N \approx 25$; in Graphite, $N \approx 115$ collisions.

3. Moderating Ratio

An ideal moderator must possess not only high slowing-down power $\xi \Sigma_s$, but also an extremely small thermal neutron absorption cross-section $\Sigma_a$. The overall figure of merit is the Moderating Ratio:

Heavy water ($D_2O$) boasts the world's highest moderating ratio ($\text{MR} \approx 5800$) due to the near-zero neutron absorption of deuterium, enabling CANDU reactors to achieve criticality using unenriched natural uranium.

1. Gas Ionization Regimes vs Applied Voltage

Gas-filled detectors consist of a cylindrical conductive chamber filled with gas (e.g., Argon + quench gas) with a central thin anode wire maintained at positive potential $V$. Ionizing radiation creates electron-ion pairs along its track (average energy required to produce one ion pair in gas is $W \approx 30 - 35\text{ eV}$). Plotting the collected pulse charge $Q$ against applied cathode-to-anode voltage reveals five characteristic operational regimes:

1. Recombination Region: Low electric field; positive ions and electrons recombine before reaching the electrodes. No steady signal.
2. Ionization Chamber Region (Regime II): Electric field is sufficient to sweep all primary ion pairs to electrodes with zero recombination. Gas multiplication factor $M = 1$. Pulse amplitude is directly proportional to deposited energy, but signals are minuscule ($10^{-15}\text{ C}$), requiring ultra-low-noise electrometers. Used for beam monitoring and gamma dosimeters.
3. Proportional Counter Region (Regime III): Near the thin anode wire ($r_a \sim 25\text{ \mu m}$), the radial electric field $E(r) = \frac{V}{r \ln(b/a)}$ exceeds $10^6\text{ V/m}$. Primary electrons gain sufficient kinetic energy between collisions to ionize gas atoms, triggering a localized Townsend avalanche. Gas gain is linear: $M \approx 10^3 - 10^5$. Pulse height is strictly proportional to initial particle energy, allowing energy spectroscopy of alpha and beta particles.
4. Limited Proportionality (Regime IV): Positive ion space charge shields the anode wire, causing non-linear saturation.
5. Geiger-Müller (GM) Region (Regime V): High electric field causes UV photons emitted in the avalanche to induce secondary avalanches throughout the entire length of the anode wire. The discharge terminates only when a sheath of slow positive ions encompasses the wire, collapsing the electric field. Gas gain reaches $M \sim 10^8$. Every ionizing event produces an identical, massive pulse ($\sim 1\text{ V}$), destroying all energy information but enabling sensitive count detection.

2. Geiger Tube Quenching & Dead Time

As the positive ion sheath drifts to the cathode wall, neutralizing ions could extract secondary electrons and trigger spurious repeat discharges. To prevent this, a quench gas is added:

• Organic Quenching: Ethanol or ethyl formate ($10\%$). Quench molecules have lower ionization potential than Argon, absorbing positive charges and dissipating energy via harmless molecular dissociation (tubes have a finite lifespan of $\sim 10^9$ counts).
• Halogen Quenching: Bromine ($Br_2$) or Chlorine ($Cl_2$). Dissociated halogen atoms spontaneously recombine, providing infinite tube operational lifetime.

The dead time $\tau$ of a GM tube is the time window ($\sim 50 - 200\text{ \mu s}$) during which a subsequent incoming particle cannot produce a pulse. For a measured count rate $R_m$, the true count rate $R_{\text{true}}$ in the non-paralyzable model is:

1. Mechanism of Inorganic and Organic Scintillators

Scintillation detectors convert the energy deposited by ionizing radiation into a burst of visible or ultraviolet fluorescence photons:

• Inorganic Scintillators (e.g., $\text{NaI(Tl)}$, $\text{CsI(Tl)}$, $\text{BGO}$, $\text{LaBr}_3\text{:Ce}$): Crystalline insulators doped with activator impurities (Thallium). Ionizing radiation excites electrons into the conduction band, creating electron-hole pairs that migrate to activator luminescence centers $\text{Tl}^+$. Transition to the ground state emits optical photons ($\lambda \approx 415\text{ nm}$ for NaI:Tl) with high light yield ($\sim 38,000\text{ photons/MeV}$) and high density ($\rho = 3.67\text{ g/cm}^3$) / high atomic number ($Z=53$), making NaI(Tl) the premier standard for gamma-ray detection.
• Organic Scintillators (Anthracene, Stilbene, Plastic Scintillators): Fluorescence arises from transitions between $\pi$-electron molecular energy levels. Possess sub-nanosecond decay times ($\tau \sim 1 - 3\text{ ns}$), making them ideal for high-speed coincidence timing and neutron detection.

2. The Photomultiplier Tube (PMT) Operation

The scintillation crystal is optically coupled to a Photomultiplier Tube (PMT) that converts the weak optical burst into a measurable electrical charge pulse:

1. Photocathode: Optical photons strike a thin photosensitive layer (e.g., bialkali Sb-K-Cs), releasing photoelectrons via the external photoelectric effect with quantum efficiency $\eta \approx 20 - 30\%$.
2. Electron Focusing: Electrostatic focusing electrodes direct photoelectrons toward the first dynode.
3. Dynode Multiplication Chain: A sequence of $N \approx 10 - 14$ dynodes at escalating positive potentials ($\Delta V \approx 100\text{ V}$ per stage). Each incoming electron releases $\delta \approx 3 - 5$ secondary electrons upon striking a dynode surface. Total electron multiplication gain is:
   $$G = \delta^N \approx (4)^{10} \approx 10^6 - 10^7$$
4. Anode Collection: The resulting packet of $\sim 10^7$ electrons is collected at the anode, producing a voltage pulse across load resistor $R_L$:
   $$V(t) = \frac{Q}{C} e^{-t / RC} = \frac{N_{pe} e G}{C} e^{-t / RC}$$
   The peak voltage amplitude is strictly proportional to the energy deposited by the incident gamma ray in the crystal.

1. Semiconductor Diode Principles (HPGe & Si Surface Barrier)

Semiconductor radiation detectors function as solid-state ionization chambers. An incoming charged particle or photon creates electron-hole pairs across the semiconductor band gap $E_g$ ($E_g = 1.12\text{ eV}$ for Silicon, $E_g = 0.67\text{ eV}$ for Germanium). The average ionization energy required to create one electron-hole pair is:

Because $w$ in semiconductors is an order of magnitude smaller than in gas ($W \approx 30\text{ eV}$) and two orders smaller than in scintillators ($W_{\text{scint}} \approx 100 - 300\text{ eV}$ per photoelectron), a given energy deposition $E$ creates a vastly larger number of charge carriers $N = E / w$, yielding unparalleled energy resolution.

• High-Purity Germanium (HPGe): Germanium crystals purified to residual impurity densities $< 10^{10}\text{ atoms/cm}^3$ enable planar and coaxial depletion depths of several centimeters. Due to small band gap ($0.67\text{ eV}$), HPGe must be operated at liquid nitrogen temperatures ($77\text{ K}$) to suppress thermal leakage currents. HPGe is the absolute gold standard for high-resolution gamma spectroscopy (FWHM $< 0.15\%$ at $1.33\text{ MeV}$).
• Silicon Detectors (Passivated Implanted Planar Silicon - PIPS): Silicon has a larger band gap, allowing room temperature operation. Used for charged particle spectroscopy (alpha, proton) and X-ray fluorescence (SDD).

2. Energy Resolution & The Fano Factor

The theoretical statistical variance in the number of created charge carriers $N$ is reduced below Poisson statistics because carrier generation events are not mutually independent (constrained by total energy conservation). This reduction is quantified by the Fano Factor ($F \approx 0.10 - 0.12$ for HPGe/Si):

The statistical Full Width at Half Maximum (FWHM) of an energy peak is:

For a $1.332\text{ MeV}$ gamma ray in HPGe ($F = 0.11, w = 2.96\text{ eV}$):

In contrast, an $\text{NaI(Tl)}$ scintillator achieves an FWHM of $\approx 80\text{ keV}$ ($6\%$) at the same energy—making HPGe over 50 times sharper.

3. Counting Statistics and Error Propagation

Radioactive disintegrations obey the Poisson distribution. For a total observed count $N$ recorded over duration $t$, the standard deviation is $\sigma_N = \sqrt{N}$, and the relative fractional uncertainty is:

When subtracting a background count $N_b$ (measured over $t_b$) from a gross sample count $N_g$ (measured over $t_g$), the net count rate $R_{\text{net}} = \frac{N_g}{t_g} - \frac{N_b}{t_b}$ has an uncertainty: