Unit 8: Nuclear Analytical Techniques: NAA, IDA & Radiocarbon Dating

Neutron Activation Analysis (NAA), developed by George de Hevesy and Hilde Levi in 1936, is an ultra-sensitive, non-destructive analytical technique for multi-elemental trace and ultra-trace determination (detection limits down to $10^{-9} - 10^{-12}\text{ g}$).

The Physical Principle of NAA

A sample containing an analyte target nucleus $^{A}_{Z}\text{X}$ is irradiated in a high thermal neutron flux $\Phi_{\text{th}}$ ($10^{12} - 10^{14}\text{ n/cm}^2\cdot\text{s}$) in a research reactor. The nucleus undergoes radiative capture $(n, \gamma)$:

The created radionuclide $^{A+1}_{Z}\text{X}$ decays with characteristic half-life $T_{1/2}$, emitting delayed gamma-ray photons with discrete energies that uniquely identify the isotope.

Mathematical Derivation of the Activation Saturation Equation

Let $N_0$ be the number of stable target nuclei in the sample, $\sigma_{\text{act}}$ the thermal neutron capture cross-section ($\text{cm}^2$), and $\Phi$ the neutron flux ($\text{n/cm}^2\cdot\text{s}$). Assuming target burnup is negligible ($\sigma \Phi \ll \lambda$):

Solving with initial condition $N^*(0) = 0$ via integrating factor $e^{\lambda t}$:

The induced activity $A(t_{\text{irr}}) = \lambda N^*(t_{\text{irr}})$ at the end of irradiation ($EOI$) is:

``` Induced Activity A(t) ▲ A_sat│─────────────────────────────────── Saturation Activity A_sat = N₀ σ Φ │ / │ / │ / A_sat/2 ────────────────────────/──────── (Irradiation Time t_irr = T₁/₂) │ / │_____________________/ └─────────────────────┴─────────────► Irradiation Time t_irr T₁/₂ ```

The term $(1 - e^{-\lambda t_{\text{irr}}})$ is the saturation factor:

• For $t_{\text{irr}} \ll T_{1/2}$: $A \approx N_0 \sigma \Phi (\lambda t_{\text{irr}})$, growing linearly.
• For $t_{\text{irr}} = T_{1/2}$: $A = 0.50 A_{\text{sat}}$ ($50\%$ of maximum achievable activity).
• For $t_{\text{irr}} = 3 T_{1/2}$: $A = 0.875 A_{\text{sat}}$ ($87.5\%$).
• For $t_{\text{irr}} \ge 6 - 7 T_{1/2}$: $A \to A_{\text{sat}} = N_0 \sigma_{\text{act}} \Phi$ (Saturation occurs; production rate equals decay rate, and further irradiation yields zero additional activity!).

Cooling and Counting Phases

Following irradiation, the sample is allowed to "cool" for decay time $t_d$ (to allow short-lived matrix interferences such as $^{28}\text{Al}$ [$2.24\text{ min}$] and $^{24}\text{Na}$ [$15.0\text{ h}$] to decay away):

The sample is then counted on an HPGe detector for counting time $t_c$. The total net gamma counts $C$ collected in the photopeak is:

where:

• $\epsilon_\gamma$ is the full-energy peak detection efficiency at that gamma energy.
• $I_\gamma$ is the absolute gamma emission probability (branching intensity).

Neutron Activation Analysis is implemented via two distinct protocols depending on matrix interferences:

1. Instrumental Neutron Activation Analysis (INAA)

In INAA, the sample is irradiated and counted purely instrumentally without chemical dissolution or separations:

• Advantages: Completely non-destructive (essential for precious archaeological artifacts, moon rocks, forensic evidence), zero reagent blank contamination, and handles hundreds of elements simultaneously.
• Protocol: Samples are typically irradiated twice:
• Short irradiation ($1 - 5\text{ min}$) for short-lived isotopes: $^{52}\text{V}$ ($3.75\text{ min}$), $^{28}\text{Al}$ ($2.24\text{ min}$), $^{56}\text{Mn}$ ($2.58\text{ h}$).
• Long irradiation ($10 - 50\text{ hours}$) followed by long cooling ($7 - 30\text{ days}$) for trace elements: $^{46}\text{Sc}$, $^{51}\text{Cr}$, $^{59}\text{Fe}$, $^{60}\text{Co}$, $^{75}\text{Se}$, $^{140}\text{La}$, $^{152}\text{Eu}$, $^{181}\text{Hf}$, $^{198}\text{Au}$.

2. Radiochemical Neutron Activation Analysis (RNAA)

When matrix radionuclides (such as $^{24}\text{Na}$, $^{38}\text{Cl}$, $^{82}\text{Br}$) emit overwhelming Compton backgrounds that drown out trace analyte peaks, post-irradiation chemical separation is performed:

• A non-radioactive stable carrier of the analyte element ($10 - 20\text{ mg}$) is added.
• The sample is digested in boiling mineral acids ($\text{HNO}_3 / \text{HF} / \text{HClO}_4$).
• The target element is separated via precipitation, solvent extraction, or ion-exchange chromatography.
• Yield corrections are determined gravimetrically using the added stable carrier.

The Standard Comparator Method

In practice, absolute quantification using the master NAA equation is difficult because neutron flux $\Phi$, cross-section $\sigma$, and absolute detector efficiency $\epsilon_\gamma$ have experimental uncertainties. Instead, analytical laboratories universally use the Standard Comparator Method: A certified reference standard containing known mass $m_{\text{std}}$ of the analyte is co-irradiated simultaneously in the exact same flux rabbit alongside the unknown sample of mass $m_{\text{unk}}$:

Because $\Phi, \sigma, t_{\text{irr}}$, and $\epsilon_\gamma$ are identical:

Correcting for cooling decay differences:

All complex nuclear parameters ($\Phi, \sigma, I_\gamma, \epsilon$) cancel completely, achieving analytical precisions better than $\pm 1\%$.

Isotope Dilution Analysis (IDA), pioneered by George de Hevesy, is a quantitative analytical method based on measuring the change in isotopic specific activity caused by mixing an unknown sample with an isotopically enriched spike.

The fundamental advantage of IDA is that quantitative (100%) chemical recovery is NOT required! Once isotopic equilibrium is achieved, any fraction of pure analyte isolated from the mixture retains the identical isotopic ratio.

1. Direct Isotope Dilution Analysis

Used to determine the unknown mass $m_x$ of a non-radioactive element in a complex sample matrix.

1. A known mass $m_1$ of the analyte labeled with a radioactive tracer having known activity $A_1$ (and specific activity $S_1 = A_1 / m_1$) is added to the sample.
2. The mixture is homogenized chemically to achieve complete isotopic exchange.
3. A small portion of pure analyte of mass $m_2$ is chemically isolated and its radioactivity $A_2$ is counted, giving diluted specific activity $S_2 = A_2 / m_2$.

By conservation of radioactivity:

Solving for the unknown mass $m_x$:

If the mass of the radioactive spike is negligible ($m_1 \ll m_x$):

``` Unknown Sample m_x (Stable) Spike m₁ with Activity A₁ (S₁ = A₁/m₁) │ │ └───────────────► MIX ◄─────────────┘ │ ▼ Complete Isotopic Equilibration Total Mass: m_x + m₁ Diluted Specific Activity: S₂ = A₁ / (m_x + m₁) │ ▼ Isolate Partial Pure Mass m₂ Measure S₂ = A₂ / m₂ ──► Compute m_x! ```

2. Reverse Isotope Dilution Analysis

Used when the substance to be determined is already radioactive ($A_x$, mass $m_x$, $S_x = A_x / m_x$). A known, large mass $m_1$ of pure non-radioactive carrier ($S_1 = 0$) is added. After equilibration and partial isolation of pure substance:

3. Substoichiometric Isotope Dilution Analysis

In classical direct IDA, measuring the diluted specific activity $S_2 = A_2 / m_2$ requires measuring both the radioactivity $A_2$ and the isolated mass $m_2$ (via gravimetry, spectrophotometry, or micro-balance). At sub-microgram trace levels, measuring $m_2$ accurately is impossible.

In 1958, Jaromír Růžička and Jiří Starý introduced Substoichiometric IDA:

1. Two solutions are prepared:

• Solution 1: Standard containing known mass $m_s$ with activity $A_s$.
• Solution 2: Sample containing unknown mass $m_x$ with identical spike activity $A_1$.

1. To both solutions, an exact, identical substoichiometric amount of chelating reagent (e.g., dithizone, cupferron, or EDTA) is added ($n_{\text{reagent}} < n_{\text{analyte}}$).
2. Because the reagent is limiting, it reacts with and extracts exactly the identical mass $m_{\text{sub}}$ of analyte from both solutions:

1. The activities of the extracted complexes ($a_x$ and $a_s$) are counted:

Because $m_{\text{sub}}$ cancels identically:

Mass determination is eliminated completely, extending detection limits to the picogram ($10^{-12}\text{ g}$) regime!

A radiometric titration is a volumetric analytical titration in which the endpoint is identified by monitoring the radioactivity of the solution or precipitate as a function of added titrant volume.

Operating Principles & Phase Separation

In radiometric titrimetry:

1. Either the analyte $A$, the titrant $B$, or an auxiliary indicator is labeled with a radioactive tracer.
2. The reaction must produce two separable physical phases (typically a solid precipitate and a liquid supernatant, or an organic solvent extraction phase).
3. After each increment of titrant, the phases are separated (via centrifugation, filtration, or settling) and the radioactivity of the supernatant phase is counted.

The Three Morphological Titration Curves

``` Case 1: Labeled Analyte Case 2: Labeled Titrant Case 3: Both Labeled Supernatant Activity Supernatant Activity Supernatant Activity ▲ ▲ ▲ A₀│\ │ / A₀│\ / │ \ │ / │ \ / │ \ │ / │ \ / │ \____________ │___________/ │ \_______/ └────┴───────────► Vol V └────┴──────────► Vol V └────┴───────► Vol V V_eq V_eq V_eq ```

1. Case 1: Only Analyte Labeled ($A^*$ + $B \to A^*B \downarrow$):

As non-radioactive titrant $B$ is added, the labeled analyte precipitates out. The activity of the supernatant decreases linearly until the equivalence point $V_{\text{eq}}$, where all analyte has precipitated. Beyond $V_{\text{eq}}$, the activity remains flat at the solubility product baseline. Example: Titration of radioactive $^{110m}\text{Ag}^+$ with non-radioactive $\text{Cl}^-$.

2. Case 2: Only Titrant Labeled ($A$ + $B^* \to AB^* \downarrow$):

Prior to equivalence, added radioactive titrant precipitates immediately with excess analyte; the supernatant activity remains near zero. Beyond $V_{\text{eq}}$, unreacted labeled titrant accumulates in the supernatant, causing activity to rise linearly. Example: Titration of non-radioactive $\text{SO}_4^{2-}$ with radioactive $^{133}\text{Ba}^{2+}$.

3. Case 3: Both Analyte and Titrant Labeled:

Activity falls linearly to equivalence and rises linearly thereafter, forming an inverted V-shaped curve.

Quantitative Advantages

• Unaffected by turbidity, intense color, or colloidal suspensions that make visual indicators useless.
• Operates at extreme dilutions ($10^{-5} - 10^{-7}\text{ M}$) where potentiometric electrodes lose Nernstian slope.
• Can be automated using continuous flow cells with scintillation counters.

Radiocarbon dating, developed by Willard F. Libby in 1949 (Nobel Prize in Chemistry, 1960), provides absolute chronometric dating of carbonaceous organic materials up to $\sim 50,000\text{ years}$ old.

Cosmogenic Production of Carbon-14

Primary galactic cosmic rays (predominantly relativistic protons) interact with upper atmospheric nuclei, creating spallation thermal neutrons. These thermal neutrons capture on atmospheric nitrogen-14 ($^{14}_{7}\text{N}$, $99.63\%$ natural abundance) via an exoergic $(n, p)$ reaction:

The production rate in the stratosphere and troposphere is approximately:

The Carbon Dynamic Reservoir & Secular Steady State

The newly formed $^{14}\text{C}$ atoms are rapidly oxidized to carbon monoxide and carbon dioxide:

The radioactive $^{14}\text{CO}_2$ mixes globally within weeks into the troposphere and dissolves into the global dynamic carbon reservoir (atmosphere, terrestrial biosphere, surface oceans, and deep oceans).

Carbon-14 decays via pure negative beta decay to stable nitrogen-14:

• Libby Half-Life (used by international convention for raw radiocarbon ages): $T_{1/2} = 5,568\pm 30\text{ years}$ ($\lambda_{\text{Libby}} = 1.2446 \times 10^{-4}\text{ yr}^{-1}$).
• Cambridge Half-Life (accurate physical value): $T_{1/2} = 5,730\pm 40\text{ years}$ ($\lambda = 1.2097 \times 10^{-4}\text{ yr}^{-1}$).

Because $^{14}\text{C}$ production has operated for millions of years ($t \gg T_{1/2}$), production and decay reached secular equilibrium:

In living equilibrium prior to the industrial revolution, the specific activity of carbon in all living organic tissue was:

Corresponding to an isotopic ratio:

The Post-Mortem Clock

Living organisms continuously assimilate $^{14}\text{CO}_2$ through photosynthesis (plants) or ingestion (animals), maintaining the steady-state isotopic ratio. Upon death, metabolic exchange halts immediately. The radioactive carbon-14 clock begins ticking as $^{14}\text{C}$ decays away:

``` Specific Activity A(t) [dpm/g C] 15.3 ▲ (Death: t = 0) │\ │ \ 7.65│--\-------------------------- T₁/₂ = 5,730 yr │ \ 3.83│----\------------------------ 2 T₁/₂ = 11,460 yr │ \ 1.91│------\---------------------- 3 T₁/₂ = 17,190 yr │ \_____________________ Background Detection Limit (~50,000 yr) └────────┴──────┴──────┴──────► Radiocarbon Age t ```

After 10 half-lives ($57,300\text{ years}$), the remaining $^{14}\text{C}$ activity is $0.015\text{ dpm/g C}$, reaching the analytical detection background.

Radiocarbon metrology underwent a revolution in the late 1970s with the invention of Accelerator Mass Spectrometry (AMS).

Radiometric Beta Counting Versus AMS

1. Beta Counting (Liquid Scintillation / Proportional Gas Counters):

Measures the radioactive decay rate $A = \lambda N$. Because the half-life of $^{14}\text{C}$ is long ($5,730\text{ years} \approx 3.0 \times 10^9\text{ minutes}$), only one atom out of every $4.3 \times 10^9$ decays each minute! To obtain 10,000 counts in 24 hours, traditional beta counting required $1.0 - 10.0\text{ grams}$ of pure elemental carbon (destroying significant portions of precious artifacts).

2. Accelerator Mass Spectrometry (AMS):

Directly counts individual $^{14}\text{C}$ atoms using an electrostatic tandem particle accelerator, bypassing radioactive decay entirely:

• Sample Mass: Requires only $0.2 - 1.0\text{ milligram}$ of carbon (a single fiber of wood, bone, or parchment)!
• Eliminating Isobaric Interference ($^{14}\text{N}$):

Stable $^{14}\text{N}$ has identical nominal mass ($A = 14$) and is $10^{12}$ times more abundant. In an AMS cesium sputter ion source, negative ions are produced. Nitrogen does not form a stable negative ion ($\text{N}^-$ is unbound), completely eliminating atomic nitrogen!

• Stripping Molecular Isobars ($^{12}\text{CH}_2^-$, $^{13}\text{CH}^-$):

The negative ions ($\text{C}^-$) accelerate into a terminal at $+2 - 3\text{ MV}$, passing through a carbon foil or argon gas stripper. Stripping strips $3 - 4$ electrons, transforming ions to $\text{C}^{3+}$ or $\text{C}^{4+}$, and completely dissociating all molecular ions via Coulomb explosion!

• High-resolution magnetic and electrostatic spectrometers count individual $^{14}\text{C}$ ions with precision better than $\pm 0.2\%$.

Dendrochronological Calibration Curves (IntCal)

Libby assumed that the atmospheric $^{14}\text{C} / ^{12}\text{C}$ production ratio was constant over geological time. In reality, it fluctuates due to:

1. Geomagnetic Field Modulation: Fluctuations in Earth's dipole magnetic moment deflect varying fractions of incoming cosmic rays.

2. Heliomagnetic Solar Cycles: Solar sunspot activity (Maunder Minimum, Spörer Minimum) modulates the solar wind magnetic shield.

3. The Suess Effect (Industrial Era): Burning vast quantities of fossil fuels (coal, oil, gas millions of years old, completely devoid of $^{14}\text{C}$) diluted atmospheric $^{14}\text{C}$ between 1850 and 1950.

4. Bomb Pulse (1950s-1960s): Atmospheric thermonuclear weapons testing nearly doubled atmospheric $^{14}\text{C}$ by 1963!

To convert conventional radiocarbon years into calibrated calendar dates ($\text{cal BC / cal AD}$), international calibration curves (e.g., IntCal20) calibrate raw radiocarbon determinations against independently dated tree rings (dendrochronology of bristlecone pines and Irish oaks extending back 14,000 years), speleothems, and varved lake sediments.

The physical identity of chemical properties between radioisotopes and stable isotopes of the same element enables the use of radiotracers across chemistry, biology, and materials science.

1. Organic and Inorganic Reaction Mechanism Elucidation

Radiotracers establish the exact bond cleavage and atom transfer pathways in chemical reactions:

• Ester Hydrolysis:

Hydrolyzing ethyl acetate with $^{18}\text{O}$-enriched water:

The heavy oxygen labels the acetic acid, proving conclusively that acyl-oxygen cleavage occurs rather than alkyl-oxygen cleavage.

• Photosynthetic Dark Reactions (Calvin Cycle):

Melvin Calvin fed unicellular green algae (Chlorella) with $^{14}\text{CO}_2$ for brief intervals ($2 - 5\text{ seconds}$), quenched the cells in boiling ethanol, and separated metabolites via 2D paper chromatography and autoradiography. He identified 3-phosphoglycerate as the initial product of photosynthetic carbon fixation (Nobel Prize in Chemistry, 1961).

2. Self-Diffusion in Solid-State Materials

Classical chemical gradients cannot measure self-diffusion (the diffusion of an element through its own crystal lattice, e.g., copper atoms diffusing through pure solid copper).

• A sub-micron layer of radioactive $^{64}\text{Cu}$ is electroplated onto a polished face of a pure copper cylinder.
• The specimen is annealed in a tube furnace at temperature $T$ for time $t$.
• Thin microtome slices of thickness $\Delta x$ are sectioned parallel to the interface and counted:

Plotting $\ln C(x)$ versus $x^2$ yields a slope of $-1/(4Dt)$, determining the self-diffusion coefficient $D$ with extreme accuracy. Repeating at multiple temperatures yields the activation energy $Q$ for vacancy migration via the Arrhenius equation:

3. Trace Equilibrium Constants and Solubility Products

Radionuclides determine solubility products ($K_{\text{sp}}$) of refractory precipitates far beyond the detection limits of atomic absorption or gravimetry. For lead sulfate ($\text{PbSO}_4$) labeled with $^{210}\text{Pb}$: Equilibrating labeled $\text{PbSO}_4$ with pure water, centrifuging, and counting the supernatant activity yields dissolved $[\text{Pb}^{2+}]$ directly down to $10^{-10}\text{ M}$.

Modern radiochemical analytical science has evolved advanced protocols that extend sensitivity and address national security metrology:

1. Epithermal Neutron Activation Analysis (ENAA)

In standard thermal NAA, high-abundance matrix elements with large thermal capture cross-sections (e.g., $^{23}\text{Na}, ^{45}\text{Sc}, ^{59}\text{Co}$) generate intense radioactivity that obscures trace elements. In ENAA, the sample is encapsulated inside a cadmium ($\text{Cd}$, thickness $1.0\text{ mm}$) or boron shield prior to irradiation. Cadmium possesses an enormous thermal capture cross-section ($\sigma_{\text{th}} \approx 20,000\text{ b}$ for $^{113}\text{Cd}$) that cuts off all neutrons below the cadmium cutoff energy ($E_{\text{Cd}} \approx 0.55\text{ eV}$), transmitting only epithermal and resonance neutrons. The Resonance Advantage Factor is:

where $I_0 = \int_{E_{\text{Cd}}}^\infty \sigma(E) \frac{dE}{E}$ is the resonance capture integral. For elements with massive resonance capture integrals (e.g., $\text{Au, U, Th, In, As, Sb, Mo}$), ENAA improves signal-to-noise ratios by factors of $50 - 500$!

2. The $k_0$-Standardization Method

In 1975, Frans De Corte developed the $k_0$-standardization method to eliminate the need for multi-element comparative standard solutions. Every reaction is calibrated relative to a single universal co-irradiated gold flux monitor ($^{197}\text{Au}$):

Because $k_0$ is a composite fundamental nuclear constant independent of experimental reactor parameters, measuring the gold monitor determines the concentration of all 65 detectable elements simultaneously from first principles.

3. Nuclear Forensics & Illicit Material Attribution

Nuclear forensics investigates interdicted nuclear materials (e.g., smuggled enriched uranium or plutonium) to trace their origin, reactor type, enrichment technology, and time since purification. Key signatures:

• Radiochronometry (Model Age): Determining the purification date by measuring daughter-to-parent decay ratios via High-Resolution ICP-MS:
• For Uranium: $^{230}\text{Th} / ^{234}\text{U}$ and $^{231}\text{Pa} / ^{235}\text{U}$.
• For Plutonium: $^{241}\text{Am} / ^{241}\text{Pu}$.
• Trace Elemental & Isotopic Fingerprinting: Minor uranium isotopes ($^{236}\text{U}$ indicates recycled reprocessed reactor fuel; $^{234}\text{U}$ depletion indicates gaseous diffusion vs centrifuge enrichment) and rare earth element distribution patterns.

Comparative Analytical Capabilities: Nuclear vs Atomic Spectrometric Techniques

| Technique | Analytical Principle | Detection Limit Range | Destructive? | Simultaneous Elements | Matrix Interferences | | :--- | :--- | :--- | :--- | :--- | :--- | | INAA | Thermal neutron $(n,\gamma)$ + HPGe | $10^{-9} - 10^{-12}\text{ g}$ | No (Non-destructive) | $40 - 60$ | Low (except high Na/Br/Cl) | | RNAA | Post-irradiation chemical separation | $10^{-12} - 10^{-14}\text{ g}$ | Yes (Destructive) | $1 - 10$ | None (Carrier chemistry) | | Direct IDA | Traced isotopic equilibration | $10^{-6} - 10^{-9}\text{ g}$ | Yes | Single element | None (Recovery independent) | | AMS | Direct atom counting via accelerator | $10^{-15} - 10^{-18}\text{ g}$ | Yes ($<1\text{ mg}$ sample) | Single isotope ($^{14}\text{C}, ^{10}\text{Be}$) | Zero isobaric background | | ICP-MS | Inductively coupled plasma ionization | $10^{-9} - 10^{-12}\text{ g}$ | Yes (Liquid digestion) | $70+$ | Spectral polyatomic isobars | | XRF | Inner-shell X-ray fluorescence | $10^{-3} - 10^{-6}\text{ g}$ | No | $30 - 40$ | Matrix absorption/enhancement |