Unit 6: Radiation Interaction with Matter & Energy Loss Mechanisms

Heavy charged particles—such as alpha particles ($^{4}\text{He}^{2+}$), protons ($p$), deuterons ($d$), and fission fragments—interact with matter almost entirely through inelastic Coulomb collisions with atomic electrons of the absorbing medium.

Because the mass of a heavy charged particle is thousands of times greater than the electron mass ($m_\alpha \approx 7300 m_e$), the projectile transfers only a tiny fraction of its kinetic energy in any single collision:

Consequently, heavy charged particles travel in essentially straight-line trajectories, losing energy continuously through hundreds of thousands of microscopic electrostatic interactions that ionize and excite absorber atoms.

The Bethe-Bloch Stopping Power Equation

The linear stopping power $-dE/dx$ (energy loss per unit path length, $\text{MeV/cm}$) of a heavy charged particle with charge $z e$ and velocity $v = \beta c$ traversing a medium with atomic number $Z$ and number density $N$ ($\text{atoms/cm}^3$) is given by the relativistic Bethe-Bloch formula:

where in SI units:

``` Stopping Power -dE/dx ▲ │ \ Relativistic Rise │ \ 1/v² Low-Energy Region (Density Effect δ) │ \ / │ \ Minimum / │ \ Ionizing / │ \__________________________________(MIP: βγ ≈ 3-4) └────────────────────────────────────────────────────────► Particle Velocity βγ ```

Physical Analysis of the Bethe-Bloch Terms

1. $1/v^2$ (or $1/\beta^2$) Dependence:

At non-relativistic velocities ($T \ll M c^2$), the stopping power is inversely proportional to the square of the projectile velocity. A slower particle spends more time in the electrostatic vicinity of each absorber electron, exerting a larger impulse and transferring more energy.

2. $z^2$ Dependence:

Stopping power scales with the square of the projectile charge. An alpha particle ($z = 2$) experiences 4 times the stopping power of a proton ($z = 1$) moving at the identical velocity.

3. Electron Density ($n_e = N Z = \rho N_A Z / A$):

Stopping power is directly proportional to the electron density of the absorbing medium.

4. Mean Excitation Potential ($I$):

The average orbital ionization energy of the target atoms, empirically parameterized by Felix Bloch:

5. Corrections:

• Shell Correction ($C/Z$): Accounts for orbital electron velocity when projectile speed is comparable to inner atomic shell electron speeds.
• Density Effect ($\delta/2$, Fermi): Relativistic dielectric polarization of the medium shields distant electrons, truncating the relativistic logarithmic rise in condensed media.

As a heavy charged particle penetrates matter, its kinetic energy steadily decreases. Because the Bethe-Bloch stopping power scales inversely with velocity ($-dE/dx \propto 1/v^2$), the rate of energy loss accelerates dramatically toward the very end of its track.

The Bragg Curve

Plotting the specific ionization (ion pairs generated per millimeter of path length) or stopping power $-dE/dx$ as a function of penetration depth produces the Bragg Curve:

``` Specific Ionization (Ion Pairs / mm) ▲ │ BRAGG PEAK │ /\ │ / \ │ / \ │ Plateau Region / \ │ ────────────────────────────────/ \ Tail (Straggling) │ \_____ └─────────────────────────────────────────────┴────► Penetration Depth x 0 R (Mean Range) ```

1. Plateau Region: Over the initial portion of the track, the high-energy particle travels at high velocity, depositing a relatively low, uniform dose.

2. Bragg Peak: As the particle slows into the keV energy regime, $-dE/dx$ surges to a sharp maximum, depositing the vast majority of its total kinetic energy within the final millimeters or micrometers of travel.

3. Sharp Drop-off: Immediately past the peak, the particle captures orbital electrons (neutralizing its charge from $z \to 0$), and its stopping power drops to zero at the mean range $R$.

Hadrontherapy Application in Oncology

The Bragg peak is the physical foundation of proton and carbon-ion radiation therapy: Unlike conventional megavoltage X-rays (which deposit their maximum dose near the skin surface and irradiate healthy tissue all the way through the patient), proton beams can be tuned in energy so that the sharp Bragg peak lands directly within deep-seated tumors (e.g., ocular or pediatric brain tumors), depositing zero exit dose in healthy tissue behind the tumor!

Range-Energy Empirical Scaling

The mean range $R$ of a heavy charged particle is the integral of reciprocal stopping power:

For alpha particles in dry air at STP ($15^\circ\text{C}, 1\text{ atm}$), the empirical Geiger Rule provides an accurate estimate for $4.0\text{ MeV} \le E_\alpha \le 8.5\text{ MeV}$:

For any other absorbing medium of density $\rho$ and effective atomic mass $A$, the range can be scaled via the Bragg-Kleeman Rule:

Using air as reference ($\rho_{\text{air}} \approx 1.225 \times 10^{-3}\text{ g/cm}^3, A_{\text{air}} \approx 14.6$):

For biological soft tissue ($\rho \approx 1.0\text{ g/cm}^3, A \approx 11$):

A $5.5\text{ MeV}$ alpha particle with a range of $4.0\text{ cm}$ in air travels only $\approx 40\,\mu\text{m}$ in soft tissue—less than the thickness of the dead cellular stratum corneum of human skin.

Unlike heavy ions, electrons ($\beta^-$) and positrons ($\beta^+$) have extremely small mass ($m_e$). When fast electrons traverse matter, they undergo two distinct, competing energy loss mechanisms:

1. Collisional (Ionization) Stopping Power

Electrons lose energy through inelastic Coulomb collisions with atomic electrons, described by the relativistic Bethe-Bloch formula modified for identical particles (Møller scattering for electrons, Bhabha scattering for positrons):

2. Radiative Stopping Power (Bremsstrahlung)

Because electrons have such tiny mass, when they pass close to a heavy nucleus with charge $+Ze$, the intense Coulomb field exerts massive centripetal acceleration $\vec{a} = \vec{F}/m_e \propto Ze / m_e$. According to classical Larmor electrodynamics, an accelerating charge radiates electromagnetic power proportional to the square of its acceleration:

Notice that Bremsstrahlung radiation is inversely proportional to the square of particle mass ($m^2$):

Bremsstrahlung is completely negligible for protons and alpha particles, but represents a dominant energy loss mechanism for fast electrons!

The radiative stopping power scales linearly with electron energy $E$ and quadratically with absorber atomic number $Z$:

``` Stopping Power -dE/dx ▲ │ / Radiative Loss (Bremsstrahlung) ~ Z² E │ / │ / │ Collisional / │ Loss ~ Z ln(E)/ │ ─────────────/───────── (Critical Energy E_c) │ / │ / └─────────────┴──────────────────────► Electron Energy E ```

The Critical Energy ($E_c$)

The ratio of radiative stopping power to collisional stopping power is given by the empirical rule:

where $E$ is in $\text{MeV}$.

The Critical Energy $E_c$ is defined as the electron energy at which radiative loss equals collisional loss:

• In Lead ($Z = 82$): $E_c \approx \frac{800}{82} \approx 9.8\text{ MeV}$. Above $10\text{ MeV}$, electrons in lead lose energy primarily by emitting Bremsstrahlung X-rays!
• In Water/Tissue ($Z_{\text{eff}} \approx 7.4$): $E_c \approx \frac{800}{7.4} \approx 108\text{ MeV}$.

Golden Rule of Radiation Shielding for Beta Emitters

To shield pure beta-emitting radioisotopes (such as $^{90}\text{Sr}/^{90}\text{Y}$ or $^{32}\text{P}$), high-$Z$ materials like lead must never be used directly! Directly placing lead around a high-energy beta source converts fast electrons into penetrating Bremsstrahlung X-ray photons via $Z^2$ scaling. Correct protocol dictates:

1. Primary shield: Low-$Z$ material (Lucite acrylic plastic, Plexiglas, or aluminum) to absorb beta particles with minimal Bremsstrahlung.
2. Secondary shield: Outer layer of lead to attenuate any residual low-energy X-rays.

Unlike alpha particles (which are ejected with sharp, discrete monoenergetic kinetic energies), beta particles are emitted with a continuous kinetic energy spectrum extending from zero up to a precise maximum endpoint energy $E_{\max}$ ($Q_\beta$).

The Crisis of Conservation Laws (1920s)

In two-body nuclear decay ($A \to B + \beta$), conservation of energy and momentum requires the ejected beta particle to carry away a unique discrete energy:

The experimental continuous spectrum (Chadwick, 1914) meant that the average kinetic energy carried by the beta particle was barely $\sim 30 - 40\%$ of $Q_\beta$. The missing energy appeared to vanish, prompting Niels Bohr to suggest that conservation of energy might hold only statistically in quantum mechanics!

``` Number of Beta Particles N(E) ▲ │ Continuous Beta Spectrum │ ___ │ / \ │ / \ │ / \ │ / \ │ / \ │_______/ \________ └───────┴──────┴─────────────────┴► Kinetic Energy T_β 0 E_avg E_max (Q_β Endpoint) ```

Wolfgang Pauli's Neutrino Hypothesis (1930)

In December 1930, Wolfgang Pauli proposed a "desperate remedy": the nucleus emits an elusive, electrically neutral, spin-$1/2$ fermion of negligible or zero rest mass alongside the beta electron:

Because three bodies share the decay energy $Q_\beta$, the electron and neutrino share the energy stochastically:

• When $E_\nu \to 0$: $T_e = E_{\max} = Q_\beta$ (the spectrum endpoint).
• When $T_e \to 0$: The neutrino carries away the entire decay energy undetected.

Enrico Fermi's Theory of Beta Decay (1934)

Enrico Fermi developed the quantum field theory of beta decay using time-dependent perturbation theory (Fermi's Golden Rule):

The transition probability per unit time for an electron to be emitted with momentum $p$ in interval $dp$ is governed by the two-body phase space volume of the electron and neutrino:

where:

• $G_F \approx 1.436 \times 10^{-62}\text{ J}\cdot\text{m}^3$ is the Fermi weak coupling constant.
• $F(Z, E)$ is the Fermi function, correcting for Coulomb attraction ($\beta^-$) or repulsion ($\beta^+$) between the outgoing electron and daughter nucleus.

The Fermi-Kurie Plot

Linearizing Fermi's spectral distribution:

Plotting $\sqrt{N(p) / [p^2 F(Z, E)]}$ versus electron kinetic energy $T_e$ yields a straight line whose horizontal axis intercept determines the exact decay endpoint energy $Q_\beta$ and neutrino mass limit ($m_\nu \approx 0$).

Unlike charged particles (which lose energy continuously in small increments), gamma-ray photons ($h\nu$) are electrically uncharged and interact with matter through discrete, catastrophic interactions in which the photon is either completely absorbed or scattered out of the beam. The three primary photon interaction mechanisms are the Photoelectric Effect, Compton Scattering, and Pair Production.

The Photoelectric Absorption Mechanism

In the photoelectric effect, an incident gamma photon strikes a tightly bound inner-shell atomic electron (predominantly the $K$-shell, $>80\%$ of events). The photon disappears completely, and its entire energy $h\nu$ is transferred to the atomic electron, which is ejected into the continuum as a photoelectron:

where $B_K$ is the binding energy of the $K$-shell electron.

Conservation of linear momentum requires that a completely free electron cannot absorb a photon and conserve both energy and momentum simultaneously; the interaction requires a bound electron where the residual target nucleus acts as a third body to absorb momentum recoil.

``` Incident Photon hν ────────────────────► ● K-Shell Electron │ ▼ Ejection Photoelectron: T_e = hν - B_K ```

Energy and Atomic Number Dependence

The atomic photoelectric cross-section $\tau$ per atom depends strongly on photon energy $E_\gamma = h\nu$ and absorber atomic number $Z$:

• High-$Z$ absorbers: Materials like lead ($Z = 82$, $Z^5 \approx 3.7 \times 10^9$) have astronomical photoelectric cross-sections compared to aluminum ($Z = 13$, $Z^5 \approx 3.7 \times 10^5$)—a factor of $10,000$ greater!
• Dominance Regime: The photoelectric effect dominates at low photon energies ($E_\gamma < 100\text{ keV}$ in tissue, $E_\gamma < 500\text{ keV}$ in lead).

Absorption Edges ($K$-Edge, $L$-Edges)

As photon energy decreases, the cross-section climbs steeply as $1/E^{3.5}$. However, when photon energy drops just below the binding energy of an electron shell ($h\nu < B_K$), photons suddenly lack sufficient energy to ionize electrons in that shell. The cross-section drops discontinuously by a factor of 5 to 10:

``` Photoelectric Cross Section τ ▲ │ /│ │ / │ │ / │ K-EDGE (hν = B_K) │ / │ │ / │___ │ L-Edges│ \ │_______/│______│____\___________ └───────┴───────┴────────────────► Photon Energy hν ```

For lead ($Z = 82$), the $K$-edge occurs at $B_K = 88.00\text{ keV}$. Photons with $88.1\text{ keV}$ are absorbed violently, while photons with $87.9\text{ keV}$ penetrate much more deeply.

De-excitation: Characteristic X-Rays Versus Auger Electrons

Photoelectric ionization leaves an inner-shell vacancy. An outer-shell electron drops into the vacancy, releasing transition energy $\Delta E = B_K - B_L$. This energy is emitted as either:

1. Characteristic X-ray photon with energy $h\nu = B_K - B_L$.

2. Auger Electron: The energy is transferred non-radiatively to an outer-shell electron, which is ejected with kinetic energy $T_{\text{Auger}} = (B_K - B_L) - B_M$.

The fluorescent yield $\omega_K$ (probability of X-ray emission vs Auger) scales with atomic number: $\omega_K \approx \frac{Z^4}{Z^4 + 10^6}$. Low-$Z$ tissue predominantly emits Auger electrons (delivering localized nanometer damage), while high-$Z$ lead emits characteristic X-rays.

At intermediate gamma-ray energies ($0.5\text{ MeV} \le E_\gamma \le 5\text{ MeV}$), the dominant interaction mechanism is Compton scattering: the elastic scattering of a photon by a loosely bound or "free" atomic electron ($h\nu \gg B_e$).

Derivation of the Compton Scattering Formula

Consider an incident photon of energy $E = h\nu$ and momentum $p = h\nu/c$ striking a stationary electron of rest mass $m_e$ at rest ($E_0 = m_e c^2$). The photon scatters at angle $\theta$ with energy $E' = h\nu'$, while the electron recoils at angle $\phi$ with kinetic energy $T_e$.

``` scattered photon hν' ^ / θ (Scattering Angle) Incident hν / ──────────────► ● (e⁻ at rest) \ \ φ v Recoil Electron T_e ```

1. Conservation of Energy:

Squaring both sides:

2. Conservation of Linear Momentum:

Squaring:

Multiplying by $c^2$:

3. Equating Momentum Expressions:

Cancelling common terms:

Dividing both sides by $2 m_e c^2 (h\nu)(h\nu')$:

Multiplying by $h c$:

where $\lambda_C = \frac{h}{m_e c} \approx 2.4263 \times 10^{-12}\text{ m} = 0.02426\text{ \AA}$ is the Compton wavelength of the electron.

Expressing scattered photon energy $E'$:

The Compton Edge Energy ($E_C$)

The recoil electron kinetic energy is:

The maximum energy transfer to the electron occurs during a head-on collision where the photon backscatters directly at $\theta = 180^\circ$ ($\cos 180^\circ = -1 \implies 1 - \cos\theta = 2$):

This sharp cutoff $E_C$ is the Compton Edge in gamma spectroscopy. The minimum energy of the backscattered photon is:

For high-energy gammas ($E \gg m_e c^2$):

Regardless of incident gamma energy, a backscattered photon can never have more than $\approx 256\text{ keV}$!

The Klein-Nishina Differential Cross-Section

Quantum electrodynamics (Oskar Klein and Yoshio Nishina, 1928) gives the differential cross-section per electron:

where $r_e = \frac{e^2}{4\pi\varepsilon_0 m_e c^2} \approx 2.818\text{ fm}$ is the classical electron radius. Because Compton scattering occurs with individual atomic electrons, the atomic cross-section scales strictly with atomic number:

At high photon energies, a third interaction mechanism appears: electron-positron pair production.

Pair Production Kinematics & Energy Threshold

In the presence of the strong Coulomb electric field of an atomic nucleus (to absorb recoil momentum), a high-energy gamma photon can spontaneously materialize into an electron-positron pair:

The threshold photon energy required for pair production in the nuclear field is exactly the sum of the rest mass energies of the two leptons:

For photon energies $h\nu < 1.022\text{ MeV}$, pair production is strictly impossible.

Any excess photon energy above $1.022\text{ MeV}$ is partitioned into the kinetic energies of the created electron and positron:

``` e⁻ (Electron, T_e⁻) ▲ / hν > 1.022 MeV ────────────► ● Nucleus (Recoil) \ \ ▼ e⁺ (Positron, T_e⁺) ──► Slows ──► Annihilation (Two 511 keV γ) ```

The atomic pair production cross-section $\kappa$ scales with the square of the nuclear charge:

Pair production dominates at high energies ($h\nu > 5\text{ MeV}$ in lead; $h\nu > 20\text{ MeV}$ in tissue).

Positron Annihilation & Escape Peaks

Once the created positron slows to thermal energies, it annihilates with an atomic electron in the medium:

To conserve momentum, the two annihilation photons are emitted collinearly back-to-back ($180^\circ$). In gamma spectroscopy, this produces characteristic peaks:

1. Full Energy Peak (Photopeak): Both $511\text{ keV}$ photons are absorbed ($E$).

2. Single Escape Peak: One $511\text{ keV}$ photon escapes the detector ($E - 511\text{ keV}$).

3. Double Escape Peak: Both annihilation photons escape ($E - 1022\text{ keV}$).

Total Linear & Mass Attenuation Coefficients

The total probability of photon interaction per unit distance is the sum of all three independent mechanisms:

where $\mu$ is the linear attenuation coefficient ($\text{cm}^{-1}$).

For a narrow, collimated monoenergetic gamma beam traversing thickness $x$:

To remove dependence on physical density $\rho$, we use the mass attenuation coefficient $\mu / \rho$ ($\text{cm}^2/\text{g}$):

where $\rho x$ is the area mass density ($\text{g/cm}^2$).

Half-Value Layer (HVL) and Tenth-Value Layer (TVL)

• Half-Value Layer (HVL): The thickness of shielding required to attenuate beam intensity by $50\%$ ($I = I_0 / 2$):

• Tenth-Value Layer (TVL): The thickness required to attenuate beam intensity by $90\%$ ($I = I_0 / 10$):

| Absorber Material | Density $\rho$ ($\text{g/cm}^3$) | Linear Attenuation $\mu$ ($1\text{ MeV}$) | Half-Value Layer HVL ($1\text{ MeV}$) | Tenth-Value Layer TVL ($1\text{ MeV}$) | | :--- | :--- | :--- | :--- | :--- | | Lead ($\text{Pb}$) | $11.35$ | $0.771\text{ cm}^{-1}$ | $0.90\text{ cm}$ ($9.0\text{ mm}$) | $2.99\text{ cm}$ | | Steel / Iron ($\text{Fe}$) | $7.87$ | $0.470\text{ cm}^{-1}$ | $1.47\text{ cm}$ | $4.90\text{ cm}$ | | Standard Concrete | $2.35$ | $0.149\text{ cm}^{-1}$ | $4.65\text{ cm}$ | $15.45\text{ cm}$ | | Water / Soft Tissue | $1.00$ | $0.0706\text{ cm}^{-1}$ | $9.82\text{ cm}$ | $32.61\text{ cm}$ |

Three half-value layers attenuate a beam to $(1/2)^3 = 12.5\%$; seven half-value layers attenuate to $<1\%$; ten half-value layers attenuate to $<0.1\%$ ($1/1024$).

Macroscopic dosimetry defines absorbed dose as an average quantity: $D = \Delta E / \Delta m$. However, within microscopic cellular targets (cell nucleus diameter $\sim 5 - 10\,\mu\text{m}$, DNA chromatin fiber diameter $\sim 30\text{ nm}$, DNA double helix diameter $\approx 2.0\text{ nm}$), radiation deposits energy in discrete, highly localized stochastic clusters known as track structures.

Formalism of ICRU Microdosimetry

In 1983, the International Commission on Radiation Units and Measurements (ICRU Report 36) established microdosimetry to describe stochastic energy deposition in sub-cellular volumes:

1. Energy Imparted ($\epsilon$): The stochastic sum of all energy transfers within a microscopic volume $V$.

2. Lineal Energy ($y$): The stochastic analogue of linear energy transfer (LET), defined as the energy imparted by a single tracking event divided by the mean chord length $\bar{l}$ of the volume:

3. Specific Energy ($z$): The stochastic analogue of absorbed dose:

As the target mass $m \to \infty$ or the number of independent tracks $n \to \infty$, the expectation value converges to macroscopic dose: $\langle z \rangle = D$.

``` Low-LET Electron Track: Sparse Ionizations High-LET Alpha Track: Dense Column ● ●●●●●●●●●●●●●●●●●●●●●●● ● ●●●●●●●●●●●●●●●●●●●●●●● ●●●●●●●●●●●●●●●●●●●●●●● ● (Clustered Damage: >10 DSBs (Simple repairable SSB) across single chromatin loop!) ```

Nanodosimetry and Clustered DNA Lesions (Complex Damage)

Monte Carlo track structure simulations (e.g., GEANT4-DNA, PARTRAC) demonstrate that:

• For low-LET radiation ($0.2\text{ keV}/\mu\text{m}$), ionizations are isolated; over $85\%$ of DNA lesions are isolated single-strand breaks or base damages that cellular repair enzymes fix with $<0.1\%$ error.
• For high-LET radiation ($100\text{ keV}/\mu\text{m}$ alpha particles), a single track traversing a cell nucleus deposits hundreds of ionizations along a continuous cylinder, producing Clustered DNA Lesions (Multiple Damaged Sites, MDS): multiple double-strand breaks, base oxidations, and abasic sites all clustered within $1 - 2$ helical turns of DNA. Cellular repair machinery (NHEJ, HR) cannot resolve clustered breaks, resulting in chromosomal fragmentation, genomic instability, and mitotic death.

Linear and Mass Attenuation Data for Photons Across Materials (0.01 to 10 MeV)

| Material | Density $\rho$ ($\text{g/cm}^3$) | $\mu/\rho$ at $0.05\text{ MeV}$ | $\mu/\rho$ at $0.10\text{ MeV}$ | $\mu/\rho$ at $0.50\text{ MeV}$ | $\mu/\rho$ at $1.00\text{ MeV}$ | $\mu/\rho$ at $5.00\text{ MeV}$ | | :--- | :--- | :--- | :--- | :--- | :--- | :--- | | Air (STP) | $0.001205$ | $0.208\text{ cm}^2/\text{g}$ | $0.155\text{ cm}^2/\text{g}$ | $0.087\text{ cm}^2/\text{g}$ | $0.0636\text{ cm}^2/\text{g}$ | $0.0275\text{ cm}^2/\text{g}$ | | Water | $1.000$ | $0.227$ | $0.171$ | $0.0969$ | $0.0707$ | $0.0303$ | | Aluminum ($Z=13$) | $2.699$ | $0.368$ | $0.170$ | $0.0845$ | $0.0615$ | $0.0284$ | | Iron ($Z=26$) | $7.874$ | $1.96$ | $0.372$ | $0.0841$ | $0.0599$ | $0.0314$ | | Lead ($Z=82$) | $11.35$ | $8.04$ | $5.55$ | $0.161$ | $0.0710$ | $0.0426$ |

At $50\text{ keV}$, lead is 35 times more attenuating per gram than water due to the $Z^4 / E^{3.5}$ photoelectric effect. At $1\text{ MeV}$, Compton scattering dominates, and mass attenuation coefficients across all materials become virtually identical ($\sim 0.06 - 0.07\text{ cm}^2/\text{g}$) because electron density per gram ($Z/A \approx 0.4 - 0.5$) is nearly constant!