Chemistry / Nuclear & Radiochemistry Decay Kinetics, Fission, Fusion & Radiometrics 100% Free Open Access
Chapter 7 β€’ Theory & Derivations

Unit 7: Radiation Detection Systems: Gas, Scintillation & Semiconductors

Comprehensive physical and electronic treatise on radiation detection instrumentation: Townsend avalanche kinetics and the six gas detector operating regimes; ionization chambers, proportional counters, and Geiger-MΓΌller tubes; organic and inorganic scintillation mechanisms; photomultiplier tubes and silicon photomultipliers; solid-state semiconductor bandgap physics; High-Purity Germanium (HPGe) gamma spectroscopy; and multichannel pulse height analysis.

Β§7.1 Gas-Filled Radiation Detectors: The Six Operating Voltage Regimes & Charge Multiplication

Gas-filled radiation detectors represent the oldest and most versatile class of radiation detection instrumentation. A gas detector consists of a gas-filled chamber containing two electrodes across which an external high-voltage electric field is applied: an outer cylindrical cathode and a thin central axial anode wire.

When ionizing radiation enters the active gas volume, it creates a quantity of primary electron-ion pairs proportional to the absorbed energy:

$$n_0 = \frac{E_{\text{dep}}}{W}$$

where $W$ is the average energy required to produce one electron-ion pair in the gas ($W \approx 26\text{ eV}$ in argon, $34\text{ eV}$ in air, $41\text{ eV}$ in helium).

The Characteristic Pulse Height Versus Voltage Curve

If an ionization event creates $n_0$ primary ion pairs, the total charge collected at the anode depends fundamentally on the applied voltage $V$:

``` Total Charge Collected Q β–² Continuous Discharge (VI) β”‚ / β”‚ Geiger-MΓΌller / β”‚ Plateau (V) / β”‚ /────────────/ β”‚ Limited / β”‚ Proportionality / β”‚ (IV) / β”‚ Proportional / β”‚ Region (III) / β”‚ /───────────────/ β”‚ Ionization / β”‚ Chamber (II)/ β”‚ β”Œβ”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”˜ β”‚/ Recombination (I) └──────────────────────────────────────────────────────► Applied Voltage V ```

1. Region I: Recombination Region:

The applied electric field is too weak ($E < 100\text{ V/cm}$) to overcome electrostatic attraction. Electrons and positive ions recombine before reaching the electrodes; charge collection is incomplete and voltage-dependent.

2. Region II: Ionization Chamber Region:

The electric field is strong enough to achieve complete charge collection ($100\%$ collection efficiency) before recombination occurs. The pulse height forms a flat saturation plateau where collected charge equals the initial ionization charge:

$$Q = n_0 e \quad (\text{Multiplication Factor } M = 1)$$

Pulse height is strictly proportional to deposited energy, but pulse amplitudes are extremely small ($\sim 10^{-14}\text{ C} \sim \text{microvolts}$), requiring sensitive electrometer amplification.

3. Region III: Proportional Region:

Near the thin central anode wire of radius $a$, the radial electric field escalates:

$$E(r) = \frac{V}{r \ln(b/a)}$$

When $E(r) > 10^4\text{ V/cm}$, free electrons gain enough kinetic energy between mean free paths to cause secondary impact ionization of gas atoms. This initiates a Townsend electron avalanche. The total collected charge is:

$$Q = M \cdot n_0 e$$

where $M$ is the gas multiplication factor ($M \sim 10^3 - 10^5$). Crucially, $M$ is independent of the initial ionization $n_0$; the output pulse height remains strictly proportional to the energy deposited by the incident radiation, enabling energy spectroscopy of alpha and beta particles.

4. Region IV: Region of Limited Proportionality:

The electron avalanche becomes so massive that the dense sheath of slow positive ions around the anode wire shields the electric field (space charge effect). Proportionality is lost.

5. Region V: Geiger-MΓΌller (GM) Region:

The electric field is so high that ultraviolet photons emitted by excited gas atoms trigger secondary avalanches along the entire length of the anode wire. A single primary electron triggers a full discharge ($M \sim 10^8 - 10^{10}$). Output pulse heights are huge ($\sim 1 - 2\text{ Volts}$) and completely independent of the energy or type of incident radiation.

6. Region VI: Continuous Discharge Region:

The electric field exceeds the dielectric breakdown strength of the gas; a continuous, damaging glow discharge occurs without radiation.

Β§7.2 Ionization Chambers & Proportional Counters: Cavity Theory & Neutron Proportional Tubes

Ionization chambers and proportional counters operate in distinct voltage regimes, serving complementary roles in radiation metrology.

Ionization Chamber Instrumentation and Cavity Theory

Ionization chambers operate with $M = 1$, measuring either individual ionization pulses (pulse mode) or the integrated steady-state saturation current (current mode):

$$I_{\text{sat}} = \left(\frac{dE}{dt}\right) \frac{e}{W} = \dot{D} \cdot m \cdot \left(\frac{e}{W}\right)$$

Ion chambers are the gold standard for reference dosimetry calibrations because their response is absolute and directly relates to the definition of radiation dose.

According to the Bragg-Gray Cavity Principle, the absorbed dose $D_{\text{med}}$ in an absorbing medium surrounding a small gas cavity is related to the dose $D_{\text{gas}}$ measured in the cavity gas:

$$D_{\text{med}} = D_{\text{gas}} \cdot \bar{s}_{\text{med, gas}} = \left(\frac{Q}{m_{\text{gas}}} \frac{W_{\text{gas}}}{e}\right) \bar{s}_{\text{med, gas}}$$

where $\bar{s}_{\text{med, gas}}$ is the ratio of mass collisional stopping powers of the medium to the gas averaged over the electron spectrum.

Proportional Counters and Counting Gases

Proportional counters typically utilize a gas mixture known as P-10:

  • $90\%$ Argon ($\text{Ar}$): Noble gas providing high ionization density and low excitation threshold.
  • $10\%$ Methane ($\text{CH}_4$): Polyatomic "quench gas" whose rotational and vibrational energy levels absorb ultraviolet photons, preventing spurious photoemission from the cathode.

Thermal Neutron Detection via Proportional Tubes

Because neutrons carry no electrical charge, they cannot ionize gases directly. Thermal neutrons are detected in proportional counters by filling the tube with gases that undergo exoergic neutron reactions yielding high-energy charged particles:

1. Boron Trifluoride ($\text{BF}_3$) Proportional Tubes:

Utilizes enriched $^{10}\text{B}$ ($96\%$ enrichment):

$$^{10}_5\text{B} + ^1_0n_{\text{th}} \longrightarrow \begin{cases} ^7_3\text{Li} + ^4_2\alpha + 2.792\text{ MeV} & (6\% \text{ to ground state}) \\ ^7_3\text{Li}^* + ^4_2\alpha + 2.310\text{ MeV} \quad (\gamma = 478\text{ keV}) & (94\% \text{ to excited state}) \end{cases}$$

The thermal capture cross-section is $\sigma_{\text{th}} = 3,840\text{ barns}$. The alpha particle ($1.47\text{ MeV}$) and lithium recoil nucleus ($0.84\text{ MeV}$) create massive ionization pulses ($\approx 80,000$ ion pairs) that dwarf background gamma pulses, providing excellent gamma-neutron discrimination via simple discriminator thresholds.

2. Helium-3 ($^3\text{He}$) Proportional Tubes:

$$^3_2\text{He} + ^1_0n_{\text{th}} \longrightarrow ^3_1\text{H} + ^1_1p + 0.764\text{ MeV} \quad (\sigma_{\text{th}} = 5,330\text{ barns})$$

Provides superior neutron detection efficiency and chemical non-toxicity, making $^3\text{He}$ tubes the global standard for neutron border monitors and nuclear safeguards.

Β§7.3 Geiger-MΓΌller Counters: Townsend Avalanches, Halogen Quenching & Dead Time Kinetics

The Geiger-MΓΌller (GM) counter is the most widely deployed portable radiation survey instrument due to its exceptional sensitivity, ruggedness, and simplicity.

The GM Avalanche Propagation Mechanism

In the GM region ($V \approx 900 - 1400\text{ V}$), the electric field near the anode wire is so high that electrons in an avalanche excite argon atoms to radiative states. Within picoseconds, these atoms de-excite by emitting UV photons ($h\nu \sim 11 - 15\text{ eV}$). Because argon is transparent to its own UV emission, these photons travel freely through the gas volume and strike other gas molecules or the cathode wall, ejecting photoelectrons. Each photoelectron initiates an independent Townsend avalanche at a new location along the wire:

``` Primary Ionization ──► Initial Avalanche ──► UV Photon Emission β”‚ β–Ό Cathode Photoemission ◄── UV Propagation ──► Secondary Avalanches β”‚ β”‚ └─────────────────► FULL WIRE DISCHARGE β—„β”˜ ```

The discharge propagates axially until a dense sheath of slow positive argon ions envelops the entire central anode wire. Because positive ions move $\sim 1000$ times slower than electrons, this positive space-charge sheath reduces the effective electric field below the threshold needed for multiplication, extinguishing the discharge.

Quenching Mechanisms: Organic Versus Halogen

When the sheath of positive ions reaches the cathode wall, they neutralize by capturing electrons:

$$\text{Ar}^+ + e^- \longrightarrow \text{Ar}^* \longrightarrow \text{Ar} + h\nu$$

The energy released equals the ionization potential of argon ($I_{\text{Ar}} = 15.76\text{ eV}$). This exceeds the work function of the cathode metal ($\Phi \approx 4 - 5\text{ eV}$), liberating secondary electrons that would reignite a spurious second pulse, resulting in continuous pulsing.

To quench this, a small percentage ($0.1 - 1\%$) of a quench gas is added:

  • Halogen Quenching: A halogen vapor (bromine $\text{Br}_2$ or chlorine $\text{Cl}_2$) is introduced. Because the ionization potential of $\text{Br}_2$ ($10.5\text{ eV}$) is lower than argon ($15.8\text{ eV}$), charge exchange occurs:
$$\text{Ar}^+ + \text{Br}_2 \longrightarrow \text{Ar} + \text{Br}_2^+$$

All positive ions arriving at the cathode are $\text{Br}_2^+$ molecules. Upon electron capture at the wall, the molecule neutralizes by dissociating into neutral atoms:

$$\text{Br}_2^+ + e^- \longrightarrow \text{Br} + \text{Br} \quad (\text{Non-radiative dissociation})$$

Unlike organic quenchers (which permanently degrade after $\sim 10^8$ counts), halogen atoms spontaneously recombine ($\text{Br} + \text{Br} \to \text{Br}_2$), giving halogen-quenched GM tubes an infinite operational lifespan!

Dead Time ($\tau$) Kinetics: Paralyzable Versus Non-Paralyzable Models

While the positive ion sheath is drifting away from the anode wire, the detector is insensitive to new incoming ionizing particles. This duration is the dead time $\tau$ (typically $\tau \approx 50 - 200\,\mu\text{s}$ in GM tubes).

Let $n$ be the true particle interaction rate and $m$ be the recorded count rate:

1. Non-Paralyzable Model:

The detector is dead for a fixed dead time $\tau$ following each recorded event. Any radiation arriving during $\tau$ is ignored and does not extend the dead time. The fraction of dead time per unit time is $m \tau$. The true count rate is:

$$m = n (1 - m \tau) \implies n = \frac{m}{1 - m \tau}$$

2. Paralyzable Model:

Each interaction extends the dead time by another period $\tau$, even if the event was not recorded. By Poisson statistics, the probability of zero events occurring in interval $\tau$ is $e^{-n \tau}$:

$$m = n e^{-n \tau}$$

At extremely high count rates ($n \to \infty$), the recorded count rate in a paralyzable detector drops toward zeroβ€”a dangerous failure mode where a catastrophic radiation field reads zero on an uncompensated meter!

Β§7.4 Scintillation Detectors: Inorganic Phosphors, Activator Luminescence & Organic Scintillators

Scintillation detectors convert the kinetic energy of ionizing radiation into a flash of optical or ultraviolet photons through luminescent excitation of transparent scintillator materials.

1. Inorganic Crystal Scintillators (Bandgap Luminescence)

Inorganic scintillators are wide-bandgap crystalline dielectrics doped with trace impurity activators. The premier gamma detection scintillator is Thallium-activated Sodium Iodide, $\text{NaI(Tl)}$:

  • High Stopping Power: High density ($\rho = 3.67\text{ g/cm}^3$) and high effective atomic number ($Z_{\text{I}} = 53$) provide high photoelectric gamma absorption.
  • Scintillation Mechanism:

Ionizing radiation promotes electrons from the crystal valence band to the conduction band, creating electron-hole pairs and free excitons. In a pure crystal, radiative de-excitation back to the valence band emits photons whose energy equals the bandgap ($E_g \approx 6\text{ eV}$ in UV), which are immediately reabsorbed by the crystal (self-absorption). Doping with $\approx 0.1\%$ thallium ($\text{Tl}^+$) creates localized energy states within the forbidden bandgap:

``` CONDUCTION BAND ══════════════════════════════════════════════ β”‚ β–² β–Ό Ionization β”‚ Exciton Migration ─────────────────────────┴──────────────────── ACTIVATOR STATES (Tl⁺) ─── Excited Level (Β³P₁) β”‚ β–Ό Optical Emission (Ξ» β‰ˆ 415 nm, Visible Blue) ─── Ground Level (ΒΉSβ‚€) ══════════════════════════════════════════════ VALENCE BAND ```

De-excitation through the thallium activator levels emits visible blue photons ($\lambda_{\max} \approx 415\text{ nm}$, $h\nu \approx 3.0\text{ eV}$). Because $3.0\text{ eV} < E_g$, the crystal is completely transparent to its own scintillation light!

  • Scintillation Light Yield: $\approx 38,000\text{ optical photons per MeV}$ of absorbed gamma energy.
  • Decay Time: Exponential decay with primary time constant $\tau \approx 230\text{ ns}$.

Other key inorganic scintillators:

  • $\text{CsI(Tl)}$: Higher stopping power ($\rho = 4.51\text{ g/cm}^3$), emission at $550\text{ nm}$ matching silicon photodiodes.
  • $\text{BGO}$ ($\text{Bi}_4\text{Ge}_3\text{O}_{12}$): Bismuth ($Z = 83$), density $7.13\text{ g/cm}^3$, used in Positron Emission Tomography (PET).
  • $\text{LaBr}_3\text{:Ce}$: Ultrafast decay ($\tau \approx 16\text{ ns}$) and exceptional energy resolution ($2.6\%$ at $662\text{ keV}$).

2. Organic Scintillators (Molecular Transitions)

Organic scintillators include aromatic hydrocarbon crystals (anthracene, stilbene), plastic polymers (PVT doped with PPO/POPOP), and liquid scintillation cocktails:

  • Luminescence Mechanism: Driven by transitions between $\pi$-electron molecular orbitals ($S_1 \to S_0$) of individual benzene rings, entirely independent of physical crystalline state.
  • Fast Decay Time: Ultrafast fluorescence lifetimes ($\tau \approx 1 - 3\text{ ns}$), ideal for sub-nanosecond timing coincidence and time-of-flight measurements.
  • Low-$Z$: Composed of carbon and hydrogen ($Z \approx 6$), minimizing photoelectric absorption; useful for beta counting and fast neutron detection via recoil protons.
  • Liquid Scintillation Counting (LSC): The radioactive analyte (e.g., low-energy beta emitters $^3\text{H}$ [$18.6\text{ keV}$] or $^{14}\text{C}$ [$156\text{ keV}$]) is dissolved directly in an aromatic liquid solvent cocktail containing primary and secondary fluor solutes, achieving $4\pi$ geometry with zero self-absorption.

Β§7.5 Photomultiplier Tubes & Optical Readout: Dynode Cascades, Quantum Efficiency & Modern SiPMs

The optical scintillation flash produced in a scintillator crystal is too faint for direct electronic measurement. The Photomultiplier Tube (PMT) converts these optical photons into a measurable electronic charge pulse with a gain of $10^6 - 10^7$.

The Photomultiplier Tube Architecture

A PMT consists of an evacuated glass envelope housing:

1. Photocathode: A thin, semitransparent layer of low-work-function photoemissive material (bialkali, e.g., $\text{Sb-Rb-Cs}$ or $\text{Sb-K-Cs}$) deposited on the interior entrance window.

Optical photons strike the photocathode and eject electrons into the vacuum via the photoelectric effect:

$$\text{Quantum Efficiency (QE)} \equiv \frac{\text{Photoelectrons Ejected}}{\text{Incident Scintillation Photons}} \approx 25 - 35\%$$

2. Focusing Electrode: Electrostatic lens directing photoelectrons toward the first dynode.

3. Dynode Electron Multiplication Cascade:

A series of $n$ curved electrodes (typically $n = 10 - 14$ dynodes) held at progressively higher positive potentials via a resistive voltage divider chain ($\Delta V \approx 100\text{ V}$ per stage). When an electron strikes dynode $i$, it liberates $\delta$ secondary electrons (secondary emission ratio $\delta \approx 4 - 6$):

$$\delta \propto (\Delta V)^k \quad (k \approx 0.7 - 0.8)$$

``` Scintillation Crystal ═════════════════════ β”‚ β”‚ β”‚ (Optical Photons hΞ½) β–Ό β–Ό β–Ό β”Œβ”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β”€β” Photocathode (QE β‰ˆ 30%) β””β”€β”€β”€β”€β”€β”€β”€β”€β”€β”¬β”€β”€β”€β”€β”€β”€β”€β”€β”€β”˜ β”‚ (Primary Photoelectrons) β–Ό ( D₁ ) ──► Ξ΄ Secondary Electrons β”‚ β–Ό ( Dβ‚‚ ) ──► δ² Electrons β”‚ β–Ό ( D₃ ) ──► δ³ Electrons ... β–Ό [ ANODE ] ──► Total Charge Pulse Q = e N_pe δⁿ ```

The total current gain $G$ of a PMT with $n$ dynode stages is:

$$G = \delta^n$$

For a 10-stage PMT with $\delta = 4.5$:

$$G = (4.5)^{10} \approx 3.4 \times 10^6$$

A single photoelectron produces an anode charge packet containing several million electrons within a rise time of $\approx 1 - 2\text{ ns}$.

Silicon Photomultipliers (SiPMs)

Modern radiation instrumentation increasingly replaces bulky, fragile vacuum PMTs with solid-state Silicon Photomultipliers (SiPMs):

  • Consists of a high-density array of thousands of micro-pixel avalanche photodiodes (APDs) connected in parallel on a common silicon substrate.
  • Operates in Geiger mode above breakdown voltage ($V_{\text{bias}} \approx 30 - 60\text{ V}$).
  • Advantages: Immune to strong magnetic fields (crucial for simultaneous PET-MRI imaging), compact footprint ($3 \times 3\text{ mm}$), low operating voltage, and high Photon Detection Efficiency (PDE $>50\%$).

Β§7.6 Semiconductor Radiation Detectors: Solid-State Band Theory, Fano Factor & Charge Collection

Semiconductor detectors function as solid-state ionization chambers. Instead of creating electron-ion pairs in a gas, ionizing radiation creates electron-hole ($e^--h^+$) pairs in a crystalline semiconductor lattice.

The Fundamental Advantage: Energy Resolution

In a gas detector, the average energy required to create one ion pair is $W \approx 30\text{ eV}$. In a scintillator-PMT combination, creating one photoelectron at the photocathode requires $\approx 100 - 300\text{ eV}$ of absorbed gamma energy. In semiconductor crystals, the bandgap $E_g$ between the valence and conduction bands is narrow:

  • Silicon ($\text{Si}$): $E_g = 1.12\text{ eV} \implies \epsilon = 3.62\text{ eV}$ per $e^--h^+$ pair.
  • Germanium ($\text{Ge}$): $E_g = 0.67\text{ eV} \implies \epsilon = 2.96\text{ eV}$ per $e^--h^+$ pair.

For an identical absorbed energy $E_{\text{dep}} = 1.0\text{ MeV}$:

  • $\text{NaI(Tl)}$ creates $\approx 4,000$ photoelectrons at the PMT photocathode.
  • Germanium creates:
$$N = \frac{1,000,000\text{ eV}}{2.96\text{ eV}} \approx 338,000\text{ electron-hole pairs}$$

Because Poisson statistical uncertainty scales as $\sigma_N / N = 1/\sqrt{N}$, having $85$ times more charge carriers reduces statistical variance dramatically, resulting in an unprecedented improvement in energy resolution!

The Fano Factor ($F$)

Because energy loss in a crystal lattice occurs through two competing channelsβ€”discrete ionization of valence electrons and excitation of non-ionizing acoustic lattice phonons (heat)β€”the individual ionization events are not statistically independent. The variance in the number of created charge carriers is reduced below the classical Poisson limit by the Fano Factor $F < 1$:

$$\sigma_N^2 = F \cdot N = F \left(\frac{E}{\epsilon}\right)$$

For germanium: $F \approx 0.08 - 0.10$; for silicon: $F \approx 0.11$.

The intrinsic statistical Full Width at Half Maximum (FWHM) energy resolution is:

$$\text{FWHM}_{\text{stat}} = 2.355 \cdot \sigma_E = 2.355 \sqrt{F \cdot \epsilon \cdot E}$$

For a $1.332\text{ MeV}$ gamma ray in germanium:

$$\text{FWHM}_{\text{stat}} = 2.355 \sqrt{(0.08)(2.96\text{ eV})(1.332 \times 10^6\text{ eV})} = 2.355 \sqrt{315,417} \approx 1.32\text{ keV}$$
$$\frac{\text{FWHM}}{E} = \frac{1.32\text{ keV}}{1332\text{ keV}} \approx 0.10\%$$

Compared to $\sim 6.0\%$ for $\text{NaI(Tl)}$, semiconductor detectors provide 60 times sharper spectral peaks, enabling the resolution of closely spaced gamma multiplets that appear as single blobs in scintillation detectors.

Β§7.7 High-Purity Germanium (HPGe) Spectrometry & Pulse Height Multichannel Analysis (MCA)

Germanium possesses a significantly higher atomic number ($Z = 32$) and density ($\rho = 5.32\text{ g/cm}^3$) than silicon ($Z = 14, \rho = 2.33\text{ g/cm}^3$), making it the premier semiconductor material for gamma-ray spectroscopy ($Z^4$ photoelectric scaling).

High-Purity Germanium (HPGe) Crystal Technology

In intrinsic semiconductors at room temperature ($T = 300\text{ K}$), the narrow bandgap ($E_g = 0.67\text{ eV}$) allows thermal energy ($k_B T \approx 0.026\text{ eV}$) to excite electrons across the gap, generating an enormous thermal leakage current that swamps radiation signals. To function as a radiation detector:

1. Ultra-Purification: The crystal must be zone-refined to unprecedented impurity concentrations ($|N_A - N_D| \le 10^{10}\text{ atoms/cm}^3$β€”less than one impurity atom per trillion germanium atoms!). This enables creation of wide depletion depths ($W > 3 - 5\text{ cm}$) at reverse bias voltages of $V \approx 2000 - 5000\text{ V}$:

$$W = \sqrt{\frac{2\varepsilon V}{e |N_A - N_D|}}$$

2. Cryogenic Cooling: HPGe detectors must be operated at liquid nitrogen temperatures ($77\text{ K}$, $-196^\circ\text{C}$) via a dewar cryostat or closed-cycle mechanical Stirling cooler to freeze out thermal charge carriers.

``` TYPICAL HPGe GAMMA-RAY PULSE HEIGHT SPECTRUM Counts β–² β”‚ PHOTOPEAK (Full Energy Absorption) β”‚ β”‚ β”‚ β–Ό β”‚ | | β”‚ /| |\ β”‚ / | | \ Single Escape Peak (E - 511 keV) β”‚ COMPTON / | | \ β”‚ β”‚ EDGE / | | \ β–Ό β”‚ | / | | \ | | β”‚ / \___/ | | \| | Backscatter Peak (~180-250 keV) β”‚ / Compton | | | | β”‚ β”‚ / Continuum | | | | β–Ό β”‚ / | | | | | | └────────────────┴─┴──────┴─┴──┴─┴────────────────► Channel (Energy E) ```

Multichannel Analyzer (MCA) Spectral Morphology

The preamplifier pulse is shaped by a spectroscopic amplifier into a semi-Gaussian voltage pulse whose peak amplitude $V_{\text{peak}} \propto E_{\text{dep}}$. The Multichannel Analyzer (MCA) digitizes $V_{\text{peak}}$ via an Analog-to-Digital Converter (ADC, typically $4096 - 16384$ channels) and increments the corresponding memory register.

The resulting gamma spectrum displays distinct physical features:

1. Photopeak (Full Energy Peak): Full absorption of $E_\gamma$ through photoelectric effect or multiple Compton events followed by photoelectric absorption.

2. Compton Continuum: Broad plateau below the photopeak corresponding to Compton scattering where the scattered photon escapes the crystal.

3. Compton Edge: Sharp drop at $E_C = \frac{2E_\gamma^2}{m_e c^2 + 2E_\gamma}$.

4. Backscatter Peak: Gammas scattering backward ($180^\circ$) from shielding/cryostat walls into the crystal ($E \approx 180 - 250\text{ keV}$).

5. Escape Peaks: For $E_\gamma > 1.022\text{ MeV}$, pair production generates two $511\text{ keV}$ annihilation photons. If one escapes: Single Escape Peak ($E - 511\text{ keV}$). If both escape: Double Escape Peak ($E - 1022\text{ keV}$).

Β§7.8 Digital Signal Processing & Pulse Shape Discrimination in Radiation Spectrometry

Modern radiation spectroscopy has transitioned from analog shaping amplifiers to high-speed Digital Signal Processing (DSP) and digital pulse processors (DPPs).

Flash ADC Digitization and Digital Trapezoidal Filtering

In digital spectrometers, the continuous exponential charge pulse emerging from a preamplifier ($V(t) = V_0 e^{-t/\tau}$) is digitized directly by a high-speed Flash ADC ($14 - 16\text{ bits}$ at $100 - 500\text{ MSamples/s}$). The digitized waveform $x[n]$ is shaped in real time using a Digital Trapezoidal Filter:

$$y[n] = y[n-1] + (x[n] - x[n-k] - x[n-l] + x[n-k-l]) + M_c \sum_{i=n-k}^{n-1} x[i]$$

where $k$ is the peaking rise time, $l$ is the flat-top duration, and $M_c$ is the pole-zero cancellation constant. The flat top eliminates pulse ballistic deficit (caused by variations in charge collection time across large HPGe crystals), while the sharp symmetric sides optimize signal-to-noise ratio at count rates exceeding $500,000\text{ cps}$ without peak shift or baseline distortion.

Pulse Shape Discrimination (PSD) for Neutron-Gamma Separation

Certain scintillation detectors (e.g., stilbene, organic liquid scintillators BC-501A/EJ-301, and dual-mode CLYC) exhibit scintillation decay profiles that depend on the ionization density of the particle. Scintillation light consists of two components:

  • Prompt Fluorescence ($\tau \sim 2 - 5\text{ ns}$): Arises from radiative decay of excited singlet states ($S_1 \to S_0$).
  • Delayed Phosphorescence ($\tau \sim 100 - 500\text{ ns}$): Arises from bimolecular annihilation of long-lived triplet states ($T_1 + T_1 \to S_1^* + S_0$).

High-LET particles (such as recoil protons from fast neutron elastic scattering) produce intense ionization density, promoting extensive triplet-triplet collisions and generating a much larger delayed light component than low-LET Compton electrons from gamma rays.

Digital spectrometers calculate the Charge Comparison Discrimination Parameter:

$$\text{PSD} = \frac{Q_{\text{tail}}}{Q_{\text{total}}} = \frac{\int_{t_{\text{gate}}}^{t_{\text{end}}} V(t) dt}{\int_0^{t_{\text{end}}} V(t) dt}$$

Plotting $\text{PSD}$ versus total energy yields two completely separated 2D branches, achieving real-time discrimination of fast neutrons from gamma backgrounds with Figures of Merit ($\text{FOM} = \Delta \text{Peak} / [\text{FWHM}_1 + \text{FWHM}_2]$) exceeding $2.5$.

Comprehensive Performance Specifications of Radiation Detectors

| Detector Category | Active Material | Average Energy per Ion/Carrier | FWHM Resolution @ $662\text{ keV}$ | Operational Temperature | Key Advantages & Applications | | :--- | :--- | :--- | :--- | :--- | :--- | | Gas Ion Chamber | Air / Ar ($M=1$) | $34.0\text{ eV}$ | Not applicable | Room temp ($300\text{ K}$) | Absolute reference dosimetry | | Proportional Tube | P-10 ($M=10^4$) | $26.0\text{ eV}$ | $\approx 12 - 15\%$ | Room temp | X-ray spectroscopy, neutron detection ($^3\text{He}$) | | GM Counter | $\text{Ar} + \text{Br}_2$ ($M=10^9$) | Full discharge | No energy resolution | Room temp | High-sensitivity survey meters | | Inorganic Scintillator | $\text{NaI(Tl)} + \text{PMT}$ | $\sim 100\text{ eV/pe}$ | $6.5 - 7.5\%$ | Room temp | Field gamma spectrometry, borehole logging | | Ultra-Dense Scintillator | $\text{BGO} (\text{Bi}_4\text{Ge}_3\text{O}_{12})$ | $\sim 300\text{ eV/pe}$ | $10 - 12\%$ | Room temp | Positron Emission Tomography (PET) | | Fast Lanthanide Scint | $\text{LaBr}_3\text{:Ce}$ | $\sim 60\text{ eV/pe}$ | $2.6 - 2.9\%$ | Room temp | Fast timing coincidence, homeland security | | HPGe Semiconductor | High-Purity Ge | $2.96\text{ eV}$ | $0.15 - 0.20\%$ ($1.2\text{ keV}$) | Cryogenic ($77\text{ K}$) | Gold standard gamma spectroscopy | | CZT Semiconductor | $\text{Cd}_{0.9}\text{Zn}_{0.1}\text{Te}$ | $4.64\text{ eV}$ | $1.5 - 2.5\%$ | Room temp | Handheld isotope identification devices |

Intermediate Example 7.1: HPGe Semiconductor FWHM Energy Resolution from Fano Factor Statistics

A High-Purity Germanium (HPGe) spectrometer has an effective Fano factor $F = 0.085$ and requires an average energy $\epsilon = 2.96\text{ eV}$ to create an electron-hole pair at $77\text{ K}$. Electronic noise contributes an independent electronic FWHM of $\text{FWHM}_{\text{noise}} = 0.85\text{ keV}$.

  1. For the $1.3325\text{ MeV}$ gamma ray of Cobalt-60 ($^{60}\text{Co}$), calculate the mean number of charge carriers $N$ generated.
  2. Determine the intrinsic statistical energy resolution $\text{FWHM}_{\text{stat}}$ in $\text{keV}$.
  3. Calculate the total overall energy resolution $\text{FWHM}_{\text{total}}$ in $\text{keV}$ and as a percentage of peak energy.

Step 1: Mean Number of Charge Carriers $N$

$$N = \frac{E}{\epsilon} = \frac{1,332,500\text{ eV}}{2.96\text{ eV}} \approx 450,169\text{ electron-hole pairs}$$

Step 2: Intrinsic Statistical Resolution

The standard deviation of carrier count:

$$\sigma_N = \sqrt{F \cdot N} = \sqrt{0.085 \times 450,169} = \sqrt{38,264} \approx 195.61\text{ pairs}$$

Energy standard deviation:

$$\sigma_E = \sigma_N \cdot \epsilon = 195.61 \times 2.96\text{ eV} \approx 579.0\text{ eV} = 0.579\text{ keV}$$

Statistical FWHM:

$$\text{FWHM}_{\text{stat}} = 2.35482 \cdot \sigma_E = 2.35482 \times 0.5790\text{ keV} \approx 1.363\text{ keV}$$

Step 3: Total Overall Energy Resolution

Electronic noise and statistical fluctuations add in quadrature:

$$\text{FWHM}_{\text{total}} = \sqrt{\text{FWHM}_{\text{stat}}^2 + \text{FWHM}_{\text{noise}}^2} = \sqrt{(1.363\text{ keV})^2 + (0.850\text{ keV})^2}$$
$$\text{FWHM}_{\text{total}} = \sqrt{1.8578 + 0.7225} = \sqrt{2.5803} \approx 1.606\text{ keV}$$

Percentage resolution:

$$\% \text{ Resolution} = \frac{\text{FWHM}_{\text{total}}}{E} \times 100\% = \frac{1.606\text{ keV}}{1332.5\text{ keV}} \times 100\% \approx 0.120\%$$

The peak width is only $1.61\text{ keV}$ ($0.12\%$), resolving narrow gamma signatures.

Easy Example 7.2: Geiger-MΓΌller Dead Time Correction for High Count Rate Radiation Survey

A non-paralyzable Geiger-MΓΌller survey meter with dead time $\tau = 120.0\,\mu\text{s}$ records a count rate of $m = 35,000\text{ counts per minute (cpm)}$ in an industrial radiography enclosure.

  1. Convert the observed count rate to counts per second ($\text{cps}$).
  2. Calculate the fraction of time the detector is dead.
  3. Compute the true incident interaction rate $n$ in $\text{cpm}$ and determine the percentage counting loss.

Step 1: Convert to Counts per Second

$$m = \frac{35,000\text{ cpm}}{60\text{ s/min}} \approx 583.33\text{ cps}$$

Step 2: Fractional Dead Time

The dead time is $\tau = 120.0\,\mu\text{s} = 1.20 \times 10^{-4}\text{ s}$.

$$m \tau = (583.33\text{ s}^{-1})(1.20 \times 10^{-4}\text{ s}) = 0.0700 \quad (7.00\%)$$

The detector is dead $7\%$ of the time.

Step 3: True Count Rate and Counting Loss

Using the non-paralyzable dead-time equation:

$$n = \frac{m}{1 - m \tau} = \frac{583.33\text{ cps}}{1 - 0.0700} = \frac{583.33}{0.9300} \approx 627.24\text{ cps}$$

In cpm:

$$n_{\text{cpm}} = 627.24 \times 60 \approx 37,634\text{ cpm}$$

Counting loss:

$$\% \text{ Loss} = \frac{n - m}{n} \times 100\% = \frac{37,634 - 35,000}{37,634} \times 100\% = \frac{2,634}{37,634} \times 100\% \approx 7.00\%$$

Uncorrected, the survey meter underestimates the radiation field by $2,634\text{ cpm}$ ($7\%$).

Intermediate Example 7.3: Two-Source Method for Experimental Determination of Detector Dead Time

In the experimental two-source method, the dead time $\tau$ of a counter is determined using two sealed sources ($S_1$ and $S_2$):

  • Background alone: $R_b = 25\text{ cpm}$
  • Source 1 alone: $R_1 = 12,450\text{ cpm}$
  • Source 2 alone: $R_2 = 15,380\text{ cpm}$
  • Sources 1 and 2 together: $R_{12} = 26,920\text{ cpm}$
  1. Convert all rates to counts per second ($\text{cps}$).
  2. Using the standard two-source formula:
$$\tau \approx \frac{R_1 + R_2 - R_{12} - R_b}{R_{12}^2 - R_1^2 - R_2^2}$$

compute the experimental dead time $\tau$ in microseconds ($\mu\text{s}$).

Step 1: Rates in Counts per Second

$$R_b = \frac{25}{60} \approx 0.417\text{ cps}$$
$$R_1 = \frac{12,450}{60} = 207.50\text{ cps}$$
$$R_2 = \frac{15,380}{60} = 256.333\text{ cps}$$
$$R_{12} = \frac{26,920}{60} = 448.667\text{ cps}$$

Step 2: Compute Dead Time $\tau$

Numerator:

$$\Delta R = R_1 + R_2 - R_{12} - R_b = 207.50 + 256.333 - 448.667 - 0.417 = 463.833 - 449.084 = 14.749\text{ cps}$$

Denominator:

$$R_{12}^2 - R_1^2 - R_2^2 = (448.667)^2 - (207.50)^2 - (256.333)^2$$
$$(448.667)^2 \approx 201,302$$
$$(207.50)^2 \approx 43,056$$
$$(256.333)^2 \approx 65,707$$
$$\text{Denominator} = 201,302 - 43,056 - 65,707 = 201,302 - 108,763 = 92,539\text{ cps}^2$$

Compute $\tau$:

$$\tau = \frac{14.749\text{ s}^{-1}}{92,539\text{ s}^{-2}} \approx 1.5938 \times 10^{-4}\text{ s} = 159.4\,\mu\text{s}$$

The detector has an experimental dead time of $\tau \approx 159\,\mu\text{s}$.

Easy Example 7.4: Photomultiplier Tube Gain and Anode Charge Pulse Amplitude

A 10-stage photomultiplier tube (PMT) has an average secondary emission ratio $\delta = 4.80$ per dynode. A scintillation event in a $\text{NaI(Tl)}$ crystal ejects $N_{pe} = 2,500\text{ photoelectrons}$ from the photocathode into the first dynode.

  1. Calculate the overall electron current gain $G = \delta^n$ of the PMT.
  2. Determine the total electrical charge $Q$ collected at the anode in picocoulombs ($\text{pC}$).
  3. If the anode load circuit has a capacitance $C = 50.0\text{ pF}$, calculate the peak voltage amplitude $V_{\text{peak}}$ of the output pulse.

Step 1: PMT Current Gain $G$

$$G = \delta^n = (4.80)^{10} \approx 6.4925 \times 10^6$$

The multiplication gain is $\approx 6.49 \times 10^6$.

Step 2: Total Anode Charge $Q$

$$Q = N_{pe} \cdot G \cdot e$$
$$Q = (2500)(6.4925 \times 10^6)(1.60218 \times 10^{-19}\text{ C}) = 1.6231 \times 10^{10} \times 1.60218 \times 10^{-19}\text{ C} \approx 2.6005 \times 10^{-9}\text{ C} = 2,600\text{ pC} = 2.60\text{ nC}$$

Step 3: Peak Voltage Amplitude

$$V_{\text{peak}} = \frac{Q}{C} = \frac{2.6005 \times 10^{-9}\text{ C}}{50.0 \times 10^{-12}\text{ F}} \approx 52.01\text{ Volts}$$

With a $50\text{ pF}$ load, the output pulse is a massive $52.0\text{ Volts}$ (in practice, smaller anode loads or preamplifiers shape this into a $0.5 - 2\text{ V}$ pulse).

Intermediate Example 7.5: Proportional Counter Electric Field Gradient and Gas Avalanche Threshold

A cylindrical proportional counter has an inner cathode radius $b = 1.50\text{ cm}$ and a central tungsten anode wire of radius $a = 25.0\,\mu\text{m}$ ($0.00250\text{ cm}$). The tube is filled with P-10 gas at 1 atmosphere, for which the threshold electric field for Townsend secondary ionization is $E_{\text{crit}} = 1.00 \times 10^4\text{ V/cm}$. The applied voltage between anode and cathode is $V_0 = 1,800\text{ V}$.

  1. Calculate the electric field at the cathode surface ($r = b$) and at the anode surface ($r = a$).
  2. Determine the critical avalanche radius $r_{\text{crit}}$ inside which secondary multiplication occurs.
  3. Express the avalanche volume as a percentage of the total detector gas volume.

Step 1: Electric Field Formula

$$E(r) = \frac{V_0}{r \ln(b/a)}$$

Calculate logarithmic term:

$$\frac{b}{a} = \frac{1.50\text{ cm}}{0.00250\text{ cm}} = 600.0$$
$$\ln\left(\frac{b}{a}\right) = \ln(600) \approx 6.3969$$

At anode surface ($r = a = 0.00250\text{ cm}$):

$$E(a) = \frac{1800\text{ V}}{(0.00250\text{ cm})(6.3969)} = \frac{1800}{0.015992} \approx 1.1255 \times 10^5\text{ V/cm}$$

At cathode surface ($r = b = 1.50\text{ cm}$):

$$E(b) = \frac{1800\text{ V}}{(1.50\text{ cm})(6.3969)} = \frac{1800}{9.5954} \approx 187.6\text{ V/cm}$$

Step 2: Critical Avalanche Radius $r_{\text{crit}}$

Set $E(r_{\text{crit}}) = E_{\text{crit}} = 1.00 \times 10^4\text{ V/cm}$:

$$r_{\text{crit}} = \frac{V_0}{E_{\text{crit}} \ln(b/a)} = \frac{1800\text{ V}}{(1.00 \times 10^4\text{ V/cm})(6.3969)} = \frac{1800}{63,969} \approx 0.02814\text{ cm} = 281.4\,\mu\text{m}$$

Step 3: Avalanche Volume Percentage

Volume scales as radius squared:

$$\frac{V_{\text{avalanche}}}{V_{\text{total}}} = \frac{\pi (r_{\text{crit}}^2 - a^2) L}{\pi (b^2 - a^2) L} \approx \left(\frac{r_{\text{crit}}}{b}\right)^2 = \left(\frac{0.02814\text{ cm}}{1.50\text{ cm}}\right)^2 = (0.01876)^2 \approx 0.000352 \quad (0.035\%)$$

Secondary multiplication is confined strictly to a microscopic sheath extending only $256\,\mu\text{m}$ from the wire, occupying less than $0.04\%$ of the tube volume! This confinement ensures all avalanches experience identical gas gain.

Intermediate Example 7.6: Scintillation Pulse Light Output and Preamplifier Statistics in NaI(Tl)

A $2.0\text{ inch} \times 2.0\text{ inch}$ $\text{NaI(Tl)}$ detector detects a $661.7\text{ keV}$ gamma ray from $^{137}\text{Cs}$.

  • Scintillation efficiency: $38.0\text{ optical photons/keV}$
  • Light collection efficiency onto photocathode: $\eta_{\text{coll}} = 75\%$
  • Photocathode quantum efficiency: $\text{QE} = 28\%$
  1. Calculate the total number of optical photons $N_{\text{phot}}$ generated in the crystal.
  2. Determine the number of photoelectrons $N_{pe}$ reaching the first dynode.
  3. Assuming Poisson statistics for photoelectron production, calculate the theoretical statistical limit of energy resolution $\text{FWHM} / E$ in percent.

Step 1: Scintillation Photons Generated

$$N_{\text{phot}} = 661.7\text{ keV} \times 38.0\text{ photons/keV} \approx 25,145\text{ optical photons}$$

Step 2: Photoelectrons Ejected

$$N_{pe} = N_{\text{phot}} \cdot \eta_{\text{coll}} \cdot \text{QE} = 25,145 \times 0.75 \times 0.28 \approx 5,280\text{ photoelectrons}$$

Step 3: Statistical Resolution Limit

The relative statistical standard deviation of photoelectron count is:

$$\frac{\sigma}{N_{pe}} = \frac{1}{\sqrt{N_{pe}}} = \frac{1}{\sqrt{5280}} = \frac{1}{72.66} \approx 0.01376$$

Statistical Full Width at Half Maximum (FWHM):

$$\frac{\text{FWHM}_{\text{stat}}}{E} = 2.355 \left(\frac{\sigma}{N_{pe}}\right) = 2.355 \times 0.01376 \approx 0.0324 \quad (3.24\%)$$

Accounting for dynode multiplication variance ($1 + 1/\delta$) and non-proportional crystal response, the actual observed FWHM of $\text{NaI(Tl)}$ at $662\text{ keV}$ is $\approx 6.5 - 7.0\%$, consistent with theoretical limits.

Advanced Example 7.7: HPGe Cryogenic Depletion Depth and Electric Bias Field

A planar p-type High-Purity Germanium detector has a net acceptor concentration $|N_A - N_D| = 1.20 \times 10^{10}\text{ cm}^{-3}$. The dielectric constant of germanium is $\varepsilon_r = 16.0$ ($\varepsilon = \varepsilon_r \varepsilon_0 = 16.0 \times 8.854 \times 10^{-14}\text{ F/cm}$).

  1. Calculate the reverse bias voltage $V$ required to achieve a full depletion depth of $W = 3.00\text{ cm}$.
  2. Determine the maximum electric field $E_{\max}$ in the crystal at this bias voltage.

Step 1: Required Reverse Bias Voltage

From semiconductor junction theory:

$$W = \sqrt{\frac{2 \varepsilon V}{e |N_A - N_D|}} \implies W^2 = \frac{2 \varepsilon V}{e |N_A - N_D|}$$
$$V = \frac{e |N_A - N_D| W^2}{2 \varepsilon}$$

Given:

  • $e = 1.60218 \times 10^{-19}\text{ C}$
  • $|N_A - N_D| = 1.20 \times 10^{10}\text{ cm}^{-3}$
  • $W = 3.00\text{ cm} \implies W^2 = 9.00\text{ cm}^2$
  • $\varepsilon = 16.0 \times 8.854 \times 10^{-14}\text{ F/cm} = 1.4166 \times 10^{-12}\text{ F/cm}$

Numerator:

$$e |N_A - N_D| W^2 = (1.60218 \times 10^{-19})(1.20 \times 10^{10})(9.00) = 1.73035 \times 10^{-8}\text{ C/cm}$$

Denominator:

$$2 \varepsilon = 2 \times 1.4166 \times 10^{-12}\text{ F/cm} = 2.8333 \times 10^{-12}\text{ F/cm}$$
$$V = \frac{1.73035 \times 10^{-8}}{2.8333 \times 10^{-12}} \approx 6,107\text{ Volts}$$

A bias of $\approx 6,100\text{ Volts}$ is required to fully deplete a $3\text{ cm}$ thick planar HPGe crystal.

Step 2: Maximum Electric Field $E_{\max}$

In a planar one-sided junction, the electric field increases linearly to a maximum at the junction interface:

$$E_{\max} = \frac{2 V}{W} = \frac{2(6107\text{ V})}{3.00\text{ cm}} \approx 4,071\text{ V/cm}$$

This field ($>1000\text{ V/cm}$) exceeds the saturation drift velocity threshold for both electrons and holes ($\approx 10^7\text{ cm/s}$), ensuring fast, complete charge collection.

Intermediate Example 7.8: Digital Pulse Processor Ballistic Deficit Elimination and Trapezoidal Peaking

In a coaxial High-Purity Germanium detector, charge collection times vary from $t_{\text{coll}} = 150\text{ ns}$ to $400\text{ ns}$ depending on whether ionizing events occur near the central core contact or the outer circumference. An analog RC-(CR) filter with shaping time $\tau = 1.0\,\mu\text{s}$ experiences a peak amplitude loss of $8.5\%$ for the slowest pulses (ballistic deficit).

  1. Explain how a digital trapezoidal filter with flat-top duration $L_{\text{flat}}$ completely eliminates ballistic deficit.
  2. Determine the minimum flat-top duration $L_{\text{flat}}$ required in microseconds.

Step 1: Mechanism of Ballistic Deficit Elimination

In analog semi-Gaussian shaping, the pulse peak occurs at a single point in time ($t_{\text{peak}} \approx 2\tau$). If charge collection is prolonged, some charge has not yet arrived when the shaping network peaks, causing a deficit in pulse height that broadens the spectral line.

A Digital Trapezoidal Filter convolves the digitized step pulse with a finite impulse response (FIR) filter having a flat plateau (flat top) of duration $L_{\text{flat}}$. The height of the flat top represents the true total integrated charge collected at the electrode, independent of when individual charge carriers arrived.

Step 2: Minimum Flat-Top Duration

To ensure complete, invariant charge collection:

$$L_{\text{flat}} \ge t_{\text{coll,\max}} - t_{\text{coll,\min}} = 400\text{ ns} - 150\text{ ns} = 250\text{ ns}$$

Adding an electronic safety margin of $100\text{ ns}$:

$$L_{\text{flat}} \ge 0.35 - 0.50\,\mu\text{s} \quad (350 - 500\text{ ns})$$

Setting the flat top to $0.50\,\mu\text{s}$ eliminates ballistic deficit completely, restoring intrinsic HPGe energy resolution even in large $150\%$ relative efficiency crystals!

Intermediate Example 7.9: Fast Neutron-Gamma Pulse Shape Discrimination Figure of Merit (FOM)

In an organic liquid scintillation detector (EJ-301), pulse shape discrimination (PSD) separates fast neutrons (recoil protons) from gamma rays (Compton electrons). The distribution of the charge ratio $\text{PSD} = Q_{\text{tail}} / Q_{\text{total}}$ reveals two Gaussian peaks:

  • Gamma peak: centroid $P_\gamma = 0.185$, $\text{FWHM}_\gamma = 0.035$
  • Neutron peak: centroid $P_n = 0.310$, $\text{FWHM}_n = 0.045$
  1. Formulate the standard Figure of Merit ($\text{FOM}$) for pulse shape discrimination:
$$\text{FOM} = \frac{P_n - P_\gamma}{\text{FWHM}_n + \text{FWHM}_\gamma}$$
  1. Compute the $\text{FOM}$ and determine whether good separation ($\text{FOM} \ge 1.25$) is achieved.
  2. Compute the peak separation distance in units of combined standard deviations ($\sigma_\gamma + \sigma_n$).

Step 1: Figure of Merit Formulation

$$\text{FOM} = \frac{\Delta \text{Peak}}{\text{FWHM}_n + \text{FWHM}_\gamma} = \frac{P_n - P_\gamma}{\text{FWHM}_n + \text{FWHM}_\gamma}$$

Step 2: Numerical Calculation

Peak separation:

$$\Delta \text{Peak} = P_n - P_\gamma = 0.310 - 0.185 = 0.125$$

Sum of FWHMs:

$$\text{FWHM}_n + \text{FWHM}_\gamma = 0.045 + 0.035 = 0.080$$
$$\text{FOM} = \frac{0.125}{0.080} = 1.5625$$

Because $\text{FOM} = 1.56 > 1.25$, the detector achieves clean, excellent discrimination with less than $0.01\%$ mutual misclassification!

Step 3: Separation in Standard Deviations

Convert FWHM to $\sigma$ ($FWHM = 2.355 \sigma$):

$$\sigma_n = \frac{0.045}{2.355} \approx 0.01911$$
$$\sigma_\gamma = \frac{0.035}{2.355} \approx 0.01486$$
$$\sigma_n + \sigma_\gamma = 0.03397$$

Number of standard deviations:

$$\text{Separation} = \frac{0.125}{0.03397} \approx 3.68 \sigma$$

The peaks are separated by nearly $3.7$ standard deviations, ensuring reliable real-time neutron detection.

Solved Honors Problems & Derivations

Step-by-step rigorous solutions with full physical, thermodynamic, and nuclear kinematic validation.