Nuclear Activation Method

• Nuclear activation method is an experimental technique that irradiates targets to measure reaction cross sections by quantifying radioactive nuclei through delayed decay signatures.
• It employs methodologies like stacked-foil irradiation, high-resolution gamma-ray spectrometry, and accelerator mass spectrometry to derive excitation functions and yield benchmarks.
• The method offers non-destructive, parallel radionuclide tracking with precise cross section measurements, though it requires meticulous calibration and accounts for target degradation.

The nuclear activation method is an experimental technique in which a target is irradiated with charged particles, neutrons, photons, or muons and the number of radioactive nuclei produced is determined from their delayed decay signatures or by direct atom counting. In nuclear astrophysics it provides a direct, angle-integrated measurement of total reaction cross sections through the off-line counting of radioactive nuclei synthesized in a target, while in charged-particle reaction studies it is commonly implemented with stacked-foil irradiation and high-resolution $\gamma$-ray spectrometry to derive excitation functions, integral yields, and benchmarks for reaction-model libraries (Gyürky et al., 2019, Tárkányi et al., 2018). Special cases include accelerator mass spectrometry (AMS), cyclic activation for short-lived positron emitters, in-beam activation with pulsed beams, and activation diagnostics that use the target itself as an internal spectrometer (Gyürky et al., 2022, Niikura et al., 2023, Matsubara et al., 9 May 2026).

1. Physical basis and governing equations

At the microscopic level, activation rests on the reaction-rate relation $R = N \sigma \phi$, or, over an energy spectrum, $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$ (Grigoriev et al., 2016). For a thin foil irradiated with constant flux $\Phi$, the standard activation equation for the activity at end of bombardment is

$$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$

with $\lambda=\ln 2/T_{1/2}$ and $N$ obtained from the target mass, molar mass, isotopic abundance, and Avogadro’s number (Tárkányi et al., 2018). After a cooling time, the activity decays according to

$A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$

and inversion gives the elemental cross section $\sigma(E)$ for the effective beam energy in the foil (Tárkányi et al., 2018).

The measured quantity is usually not $A_{\rm EOB}$ directly but the net full-energy peak counts $R = N \sigma \phi$0 under a selected $R = N \sigma \phi$1 line. For irradiation time $R = N \sigma \phi$2, cooling time $R = N \sigma \phi$3, counting time $R = N \sigma \phi$4, detector efficiency $R = N \sigma \phi$5, and emission probability $R = N \sigma \phi$6, one widely used form is

$R = N \sigma \phi$7

which yields

$R = N \sigma \phi$8

for a thin target and constant flux (Tárkányi et al., 2013). If the beam current is not perfectly constant, the irradiation can be subdivided into short intervals and the saturation factor replaced by the corresponding time-resolved sum over measured fluxes (Gyürky et al., 2019).

Activation is not restricted to independent production. When the nuclide of interest is fed by the decay of a parent, cumulative cross sections require explicit treatment of in-growth and decay of both parent and daughter through Bateman-type relations (Ditrói et al., 2013). This is one reason why cooling-time planning and repeated counting sequences are central to the method.

2. Irradiation architectures and target realization

The characteristic implementation for charged-particle work is the stacked-foil irradiation technique. A single irradiation contains many thin targets separated by degraders and monitor foils, so that successive foils experience progressively lower projectile energies. Representative arrangements include alternating CeO$R = N \sigma \phi$9 targets and monitor foils, repeated natYb/Al stacks, natRe alloy foils interleaved with Al or Ti monitors, and natLa foils combined with Al degraders and monitors (Tárkányi et al., 2018, Tárkányi et al., 2013, Ditrói et al., 2013, Tárkányi et al., 2018). The same logic is used with metal foils, oxide pellets on backings, and sedimented oxide layers.

Projectile-energy degradation through the stack is calculated from Andersen–Ziegler stopping powers and then refined by simultaneously measured monitor reactions (Tárkányi et al., 2018). Frequently used monitor channels are $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$0, $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$1, $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$2, and $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$3, compared with IAEA recommended excitation functions to fix both the energy scale and the absolute particle flux (Tárkányi et al., 2018, Tárkányi et al., 2013). The effective beam energy in a foil is then assigned as the mean, midpoint, or median energy in that foil segment, depending on the analysis convention used in a given study (Ditrói et al., 2013, Tárkányi et al., 2014).

Beam-current measurement is typically provided by a Faraday cup or a Faraday-cup-like holder with secondary-electron suppression, with monitor reactions used as an external consistency correction (Tárkányi et al., 2018, Ditrói et al., 2013). Thermal and chemical stability are treated as practical design constraints: water-cooled holders, He atmospheres, beamline vacuum, and careful control of foil adhesion are used to avoid overheating, oxidation, or flaking (Tárkányi et al., 2013, Tárkányi et al., 2014, Tárkányi et al., 2018). In nuclear astrophysics, the same activation logic extends beyond solid foils to thin evaporated or sputtered targets, anodized films, windowless gas targets, differential pumping cells, and sealed gas cells with thin entrance foils (Gyürky et al., 2019).

3. Assay of activation products

High-resolution $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$4-ray spectrometry with HPGe detectors is the dominant assay route. Energy calibration is performed with standard multi-$R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$5 sources such as $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$6Eu, $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$7Ba, $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$8Co, and $R = \int_0^\infty N(E)\,\sigma(E)\,\phi(E)\,dE$9Cs, and the full-energy peak efficiency $\Phi$0 is determined for the same geometry, often by a polynomial fit in log–log space (Tárkányi et al., 2013, Ditrói et al., 2013, Tárkányi et al., 2014). Counting is usually carried out in multiple time series at different detector-to-sample distances so that strong lines, weak lines, and nuclides with different half-lives can all be measured without excessive dead time (Tárkányi et al., 2018, Ditrói et al., 2013).

Spectrum analysis consists of assigning $\Phi$1 lines from NuDat or NUDAT 2.6, fitting the peaks with Gaussian-plus-background models, extracting net peak areas, and applying corrections for decay during counting, dead time, pile-up, and, where needed, coincidence summing (Tárkányi et al., 2018, Tárkányi et al., 2013, Tárkányi et al., 2014). Repeated spectra at different cooling times are also used to identify line interferences and to separate independent, cumulative, and isomeric contributions (Ditrói et al., 2013).

The activation method is not confined to delayed $\Phi$2 rays. For electron-capture products, characteristic X-rays can be counted with LEPS-type Ge detectors; $\Phi$3 emitters can be assayed through the 511 keV annihilation photons; $\Phi$4 decays can be measured with surface-barrier detectors or etched-track detectors; and conversion or $\Phi$5 electrons can be recorded with dedicated electron spectrometers (Gyürky et al., 2019). The $\Phi$6 study is a particularly clear example of positron-emitter activation: the total cross section was extracted from the annihilation radiation following the $\Phi$7 decay of $\Phi$8 using cyclic irradiation and counting windows tailored to its 122.24 s half-life (Gyürky et al., 2022).

AMS constitutes a special activation case in which the number of produced radionuclides is measured directly as an atomic ratio $\Phi$9, giving $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$0 for total fluence $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$1 (Gyürky et al., 2019). This route is especially valuable for very long-lived products or for cases with no convenient decay radiation, and it is explicitly independent of decay schemes and half-lives (Gyürky et al., 2019).

4. Variants and extended meanings of activation

Several distinct variants have emerged from the same underlying formalism.

Variant	Observable	Characteristic use
Off-line activation	Delayed $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$2, X-ray, $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$3, or $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$4 activity	Cross sections and yields in astrophysics and charged-particle studies (Gyürky et al., 2019, Tárkányi et al., 2018)
AMS activation	Atomic ratio $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$5	Very long-lived products or weak decay signatures (Gyürky et al., 2019)
In-beam activation	Delayed $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$6 rays during interpulse windows	Muon-capture branching ratios and short-lived products (Niikura et al., 2023)
Internal activation diagnostic	Absolute yields of $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$7C and $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$8Be produced in the target	Reconstruction of in-solid proton spectra in laser-driven fusion (Matsubara et al., 9 May 2026)

In-beam activation is a variant of the classic method in which the radioactive decay of nuclei produced on-line during beam irradiation is measured between beam pulses rather than only after end of bombardment. In the palladium muon-capture study, the 20 ms interpulse period at ISIS/RAL served as a low-background decay-spectroscopy window, enabling production branching ratios for residual nuclei spanning half-lives from milliseconds to many hours in a single irradiation run (Niikura et al., 2023). This variant retains the activation logic but relocates the counting interval into the irradiation sequence itself.

The internal activation diagnostic used in laser-driven proton–boron fusion is a different extension. There, the boron target becomes an internal spectrometer: the absolute yields of $$A_{\rm EOB}=N\,\sigma\,\Phi\,(1-e^{-\lambda t_{\rm irr}}),$$9 from $\lambda=\ln 2/T_{1/2}$0 and $\lambda=\ln 2/T_{1/2}$1 from $\lambda=\ln 2/T_{1/2}$2 provide two independent moments of the in-solid proton energy distribution, from which an exponential-equivalent spectrum can be reconstructed and propagated to the side-channel fusion yield (Matsubara et al., 9 May 2026). This suggests that activation can function not only as a cross-section method but also as an inverse diagnostic of particle transport inside dense matter.

A separate terminological point is that “nuclear activation” can also denote externally activated fission or fusion in subcritical systems. In that usage, an external source $\lambda=\ln 2/T_{1/2}$3 seeds reactions in a blanket or subcritical core, with the steady-state flux scaling as $\lambda=\ln 2/T_{1/2}$4 in a one-group approximation (Grigoriev et al., 2016). This is conceptually distinct from the measurement method used in activation cross-section studies, although both rely on the same reaction-rate formalism.

5. Data reduction, uncertainty propagation, and model benchmarking

The extraction of excitation functions is a foil-by-foil procedure. For each foil, the activity at the start of counting is inferred from the net peak area, efficiency, line intensity, and decay corrections; the number of target atoms $\lambda=\ln 2/T_{1/2}$5 is computed from mass, areal density, composition, or isotopic abundance; the particle flux $\lambda=\ln 2/T_{1/2}$6 is determined from monitor reactions; and the cross section is assigned to the effective energy in that foil (Tárkányi et al., 2018, Tárkányi et al., 2014). When multiple $\lambda=\ln 2/T_{1/2}$7 lines are available, weighted averages are formed; when overlapping irradiations exist, mutual agreement in the overlap region is checked (Tárkányi et al., 2014).

Uncertainty treatment is usually based on quadratic summation of independent relative contributions. Reported one-sigma components include beam current or flux determination of $\lambda=\ln 2/T_{1/2}$8–$\lambda=\ln 2/T_{1/2}$9, target thickness or target-nuclei determination of $N$0–$N$1, detector efficiency of $N$2–$N$3, $N$4-emission probabilities of $N$5–$N$6 or, in some datasets, decay-data contributions up to $N$7, and counting statistics from $N$8–$N$9 depending on peak intensity (Tárkányi et al., 2013, Ditrói et al., 2013, Tárkányi et al., 2014, Tárkányi et al., 2018, Gyürky et al., 2022). Corresponding total cross-section uncertainties reported in the cited studies range from about $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$0 to $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$1 (Tárkányi et al., 2018, Ditrói et al., 2013, Tárkányi et al., 2014).

Integral or thick-target yields are obtained by integrating the measured excitation function over the stopping trajectory,

$A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$2

or equivalent forms using stopping powers $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$3 and numerical spline interpolation (Tárkányi et al., 2013, Ditrói et al., 2013, Ditrói et al., 2013). In the rhenium deuteron study, the same cross-section information was further transformed into thin-layer activation curves $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$4 by combining the energy dependence of $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$5 with the depth–energy relation from stopping powers (Ditrói et al., 2013).

A major role of activation measurements is the benchmarking of nuclear-reaction models and evaluated libraries. Experimental excitation functions from activation studies have been compared with TALYS-based TENDL libraries, ALICE-IPPE, EMPIRE, PHITS, and intermediate- or high-energy transport frameworks combining MCNPX with CASCADE/INPE (Tárkányi et al., 2018, Tárkányi et al., 2013, Ditrói et al., 2013, Moskalenko et al., 2021, Niikura et al., 2023). The comparisons are often only moderately successful: low-energy peaks may be reproduced while high-energy tails are underestimated, isomeric ratios can differ significantly from experiment, and discrepancies by a factor of two or more are reported for some channels (Ditrói et al., 2013, Tárkányi et al., 2014). At intermediate and high energies, multi-criteria model selection based on $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$6, RMS-type measures, and composite scores has been used to choose models over broad residual-nucleus domains (Moskalenko et al., 2021).

6. Domains of use, strengths, and limitations

The method is established across several research areas. In medical-isotope and activation studies it has been used to measure proton- and deuteron-induced excitation functions on Ce, La, Nd, Yb, Mn, and Re, with applications to therapeutic radionuclide production, thin-layer activation, and production-yield assessment (Tárkányi et al., 2018, Tárkányi et al., 2018, Tárkányi et al., 2014, Tárkányi et al., 2013, Ditrói et al., 2013, Ditrói et al., 2013). In nuclear astrophysics it is described as a cornerstone technique because it directly measures total cross sections for reactions that produce radioactive final nuclei, and it is especially useful when in-beam particle or $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$7 detection is difficult (Gyürky et al., 2019). The $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$8 measurement further shows that activation can anchor R-matrix fits for extrapolation to astrophysical energies while remaining largely independent of branching-ratio and angular-distribution issues that affect prompt $A(t)=A_{\rm EOB}e^{-\lambda t_{\rm cool}},$9 methods (Gyürky et al., 2022).

Activation is also central to activation and transmutation inventory analysis. In the Fispact-II context, pathway-reduced analysis and the Pathways-Based Metric rank the reactions most important for final activity or inventory uncertainty, helping to identify nuclear data whose refinement would most reduce uncertainty in predicted inventories (Arter et al., 2015). This suggests a broader methodological role for activation data: not only measuring $\sigma(E)$0, but also structuring downstream sensitivity and transmutation analyses.

Its strengths are recurrent across the cited literature. The method is non-destructive, multiple radionuclides can be followed in parallel, the measured quantity is the total, angle-integrated production cross section, and AMS provides a route entirely independent of decay schemes and half-lives (Tárkányi et al., 2018, Gyürky et al., 2019). Repeated cooling-counting cycles allow separation of short- and long-lived species and of isomeric or parent–daughter contributions (Tárkányi et al., 2018, Ditrói et al., 2013). In-beam activation extends the accessible half-life range toward the millisecond regime when an appropriate pulsed beam is available (Niikura et al., 2023).

The limitations are equally specific. Classic activation requires either a radioactive reaction product with an observable decay signature or access to AMS; it does not by itself provide information on individual partial transitions when only the total activity is measured (Gyürky et al., 2019, Gyürky et al., 2022). Accurate results depend on reliable stopping powers, monitor reactions, detector calibration, and control of interferences, self-absorption, beam losses, and target degradation (Tárkányi et al., 2018, Ditrói et al., 2013, Gyürky et al., 2019). Very short-lived products can be inaccessible in off-line work, and in-beam activation is not directly applicable at DC sources and cannot observe products with $\sigma(E)$1 in the excluded recovery window (Niikura et al., 2023). These constraints explain why activation studies emphasize cooling-time planning, geometry control, repeated calibrations, and full uncertainty propagation as best practices (Tárkányi et al., 2018, Tárkányi et al., 2018).