Experimental Characterization of Neutron Capture Kinetics and Protactinium-233 Decay Dynamics in Accelerator-Driven ²³²Th→²³³U Transmutation Systems
2026
Abstract
This article develops an experimentally grounded characterization framework for neutron-capture kinetics and protactinium-233 decay dynamics in accelerator-driven transmutation systems based on the conversion chain ²³²Th(n,γ)²³³Th → ²³³Pa → ²³³U. The analysis connects reaction-rate measurement, subcritical source-driven kinetics, gamma-spectrometric activation, and covariance-aware nuclear data evaluation. The work is not presented as a proprietary irradiation campaign; it is a reproducible, literature-anchored experimental protocol supported by public evaluated libraries, benchmark experiments, and open nuclear-data services.
The central physical difficulty is that ²³³Pa is neither an instantaneous decay product nor a passive intermediate. Its half-life of about 27 days places it on the same timescale as many irradiation, cooling, and measurement cycles, while its neutron-capture channel competes with beta decay to ²³³U. Therefore, the bred fissile yield cannot be inferred from ²³²Th capture alone. It must be determined from a coupled network in which source intensity, neutron spectrum, cross-section structure, and decay timing are treated simultaneously.
A set of rate equations, measurement equations, dimensionless competition parameters, and uncertainty-propagation expressions is derived in LaTeX form. Experimental observables are mapped to detector systems: HPGe spectroscopy for activation products, fission chambers for fast-spectrum fission response, monitor foils for flux normalization, and pulsed-neutron diagnostics for subcritical reactivity. The proposed data-reduction chain is compatible with ENDF/B, JEFF, JENDL, TENDL, EXFOR, JANIS, and IAEA/NDS data workflows.
The resulting protocol shows that accelerator-driven ²³²Th → ²³³U studies require simultaneous reporting of effective one-group reaction rates, spectral weighting, irradiation history, Pa-233 cooling intervals, detector efficiency, covariance assumptions, and model validation against integral experiments. A defensible characterization of thorium breeding is consequently a kinetic measurement problem, not merely an isotope-production calculation.
Keywords
accelerator-driven system; thorium fuel cycle; neutron capture; protactinium-233; uranium-233; subcritical kinetics; activation analysis; gamma spectrometry; evaluated nuclear data; transmutation
1. Introduction and scientific scope
The ²³²Th → ²³³U fuel cycle is attractive because it couples a fertile nuclide abundant in nature with a fissile product that has favorable thermal-neutron fission properties. In an accelerator-driven system (ADS), the fertile blanket is embedded in a subcritical neutron-multiplying environment supplied by an external source. This architecture separates the neutron-source term from criticality and allows operation in regimes that are difficult to realize in a self-sustained reactor. The same separation, however, imposes a more demanding experimental problem: the neutron source, the subcritical multiplication, and the isotope-production network must be characterized as a single dynamic system.
The conversion chain begins with radiative capture in thorium, proceeds through the short-lived ²³³Th precursor, and accumulates ²³³Pa before beta decay generates ²³³U. The basic nuclear path is compact, but its experimental interpretation is not. Neutron spectra in ADS blankets are typically broad; high-energy spallation neutrons are moderated, scattered, absorbed, or leaked before they contribute to capture in the fertile region. The effective production rate is therefore an integral over energy, position, and irradiation history rather than a single microscopic number.
A rigorous characterization must distinguish three quantities that are often conflated in simplified discussions: the microscopic cross section \(\sigma(E)\), the spectrum-weighted effective cross section \(\langle\sigma\rangle\), and the measured reaction rate \(R\). Only the last of these is directly observable in an experiment, while the first two depend on data processing and spectral assumptions. This distinction is especially important for ²³³Pa because its capture-to-decay competition can change the final ²³³U yield without noticeably changing the initial ²³²Th capture rate.
The aim of the study is to formulate a complete experimental characterization protocol. It does not rely on invented irradiation results. Instead, it combines the physics of Bateman chains, source-driven point kinetics, activation analysis, and nuclear-data validation. The evidence base includes IAEA thorium fuel-cycle assessments, OECD/NEA ADS reports, evaluated nuclear data libraries, EXFOR experimental reaction data, and major ADS benchmark programs such as MUSE, YALINA, GUINEVERE/VENUS-F, and MYRRHA-oriented research.
The term experimental characterization is used here in the strict metrological sense: it denotes the set of measurements, corrections, uncertainty budgets, and validation benchmarks required to determine the kinetic parameters of the thorium-to-uranium conversion chain under a specified ADS irradiation field. The same framework can be applied to a physical experiment, a benchmark reconstruction, or a design-stage mock-up, provided that all assumptions are reported and all reaction-rate inferences remain traceable to measured or evaluated data.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
2. Literature base and data provenance
The thorium cycle has been reviewed by IAEA and OECD/NEA publications as a fuel-cycle option with potential advantages and substantial implementation challenges. These reports emphasize that ²³²Th is fertile rather than fissile and that ²³³U breeding depends on neutron economy, reprocessing strategy, and protactinium management. In contrast to a purely conceptual fuel-cycle discussion, the present analysis focuses on the experimental observables that allow these claims to be tested in an ADS environment.
Evaluated nuclear data are the quantitative foundation of the study. ENDF/B-VIII.0, ENDF/B-VIII.1, JEFF-3.3, JENDL-5, and TENDL provide reaction files that can be processed into application-specific group constants. None of these libraries should be treated as a perfect measurement. Their value lies in traceability, covariance information where available, and comparison across independent evaluation methodologies. Recent work on ²³³U neutron capture in the keV region illustrates why library comparison remains essential even for a historically important fissile isotope.
EXFOR provides the experimental reaction-data archive against which evaluated files are assembled and tested. JANIS supplies a practical route for inspecting, comparing, and averaging evaluated and experimental cross-section curves. For an ADS thorium experiment, these tools are not optional appendices; they are part of the data chain. The spectrum-weighted capture rate in a real blanket depends on resonance treatment, unresolved resonance assumptions, self-shielding correction, and the energy distribution created by the spallation source and the surrounding core.
Integral experiments supply an additional layer of validation. MUSE, YALINA, GUINEVERE, and VENUS-F programs were not designed as thorium-breeding production plants, but they are highly relevant because they test source-coupled subcritical kinetics, neutron multiplication, reactivity monitoring, and detector response. A thorium ADS experiment that ignores these benchmarks risks producing isotope inventories that are mathematically precise but physically unvalidated.
The article therefore uses a two-level evidence model. Differential data determine microscopic reaction probabilities and decay constants, while integral benchmarks constrain transport, leakage, spectrum, source efficiency, and kinetic response. The resulting measurement strategy is deliberately conservative: no single library, detector, or reactor-physics model is allowed to dominate the inference without cross-checks.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
| Source class | Examples | Role |
|---|---|---|
| Evaluated libraries | ENDF/B-VIII.1, ENDF/B-VIII.0, JEFF-3.3, JENDL-5, TENDL | cross-section processing and covariance comparison |
| Experimental archives | EXFOR, n_TOF publications, activation data | differential evidence and benchmark constraints |
| Integral benchmarks | MUSE, YALINA, GUINEVERE, VENUS-F, MYRRHA-oriented studies | source-coupled kinetics validation |
| Decay data | IAEA LiveChart, NNDC NuDat/ENSDF | half-lives, decay modes, gamma-line selection |
3. Nuclear reaction network
The primary breeding path is the sequence ²³²Th(n,γ)²³³Th followed by two beta-minus decays. Because ²³³Th has a half-life on the order of tens of minutes, it is often eliminated from long-timescale reactor balances. That approximation can be acceptable for multi-day irradiation analysis, but it is not automatically valid for pulsed irradiation, short activation experiments, or measurements performed during the early cooling interval. The full chain must therefore be retained at the model-building stage.
The minimum rate network contains ²³²Th, ²³³Th, ²³³Pa, and ²³³U. Additional loss terms may include ²³³Pa(n,γ), ²³³U(n,γ), ²³³U(n,f), leakage-adjusted flux gradients, and chemical separation factors. If the experiment uses a high-energy source, threshold reactions and spallation residues may also be relevant to radiological background, although they do not replace the principal capture-decay chain.
The most important modeling choice is the definition of the neutron flux. In a one-group approximation, \(\Phi\) represents the spectrum-integrated scalar flux weighted consistently with the processed cross section. In a multi-group calculation, every reaction rate is a sum over groups. Both descriptions are acceptable if their meaning is explicit. Ambiguity in \(\Phi\) is a common source of apparent disagreement between activation results and transport calculations.
For a homogeneous region with constant flux, the first-order system can be written as follows. The notation uses pure LaTeX so that the equations remain portable between word-processing and typesetting environments.
\[ \begin{aligned} \frac{dN_{Th233}}{dt}=\langle\sigma_{Th232,\gamma}\rangle\Phi N_{Th232}-\lambda_{Th233}N_{Th233}-\langle\sigma_{Th233,a}\rangle\Phi N_{Th233}\\ \frac{dN_{Pa233}}{dt}=\lambda_{Th233}N_{Th233}-\lambda_{Pa233}N_{Pa233}-\langle\sigma_{Pa233,\gamma}\rangle\Phi N_{Pa233}\\ \frac{dN_{U233}}{dt}=\lambda_{Pa233}N_{Pa233}-\langle\sigma_{U233,f}\rangle\Phi N_{U233}-\langle\sigma_{U233,\gamma}\rangle\Phi N_{U233} \end{aligned} \]
A multi-group implementation replaces each one-group coefficient by \(k_x=\sum_g\sigma_{x,g}\Phi_g\). If covariance matrices are available, the uncertainty of \(k_x\) must be evaluated using the same group structure. This is important because the thorium resonance region can dominate the response even when the total flux is not largest in that energy band.
The target depletion of ²³²Th is usually negligible in short experiments, but retaining the term is good practice. It prevents hidden assumptions when the same equations are later applied to higher fluence or to blanket-scale calculations. The model should state whether \(N_{Th232}\) is constant or time-dependent.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
| Nuclide | Role | Characterization issue | Sources |
|---|---|---|---|
| ²³²Th | fertile target nuclide | alpha half-life about 1.405e10 y; neutron capture creates Th-233 | [1]-[3], [7]-[9] |
| ²³³Th | short-lived precursor | beta-minus half-life about 22 min; normally treated as transient in daily irradiation balances | [7], [8] |
| ²³³Pa | intermediate beta emitter | beta-minus half-life about 26.97 d; neutron capture competes with decay during irradiation | [1]-[3], [7], [36] |
| ²³³U | bred fissile isotope | long-lived alpha emitter; fission and capture terms control utilization in a neutron field | [1]-[3], [7]-[10], [30] |
4. Source-coupled ADS kinetics
In an ADS, neutron population is maintained by an external source rather than by critical self-sustainment. The subcritical core multiplies source neutrons, and the effective multiplication depends on \(k_{eff}\), source importance, leakage, and delayed-neutron kinetics. Isotope production in the thorium blanket therefore depends not only on the accelerator beam current but also on the coupling between the source and the multiplying medium.
A practical characterization cannot treat beam power as a proxy for fertile capture. Two systems with the same beam current can have different thorium capture rates if the target material, source spectrum, reflector, blanket geometry, or subcriticality level differs. The experimental quantity of interest is the local reaction-rate density, not the nominal source strength. This is why activation foils, local flux monitors, and transport benchmarks must accompany any Pa-233 inventory measurement.
The point-kinetics equations for a source-driven subcritical assembly are useful for interpreting time-domain measurements. They do not replace spatial transport, but they identify the kinetic parameters that pulsed-neutron experiments attempt to infer. The external source term \(S(t)\) must be retained explicitly. In a pulsed experiment, activation during beam-on and decay during beam-off intervals should be integrated over the actual time structure rather than approximated by a continuous average unless the half-lives and reaction timescales justify that approximation.
The canonical source-driven point kinetics form is:
\[ \begin{aligned} \frac{dn(t)}{dt}=\frac{\rho-\beta_{eff}}{\Lambda}n(t)+\sum_{i=1}^{m}\lambda_i C_i(t)+S(t)\\ \frac{dC_i(t)}{dt}=\frac{\beta_i}{\Lambda}n(t)-\lambda_i C_i(t)\\ M_{sub}=\frac{1}{1-k_{eff}}, \qquad R_{Th}=\int_V\int_0^{\infty}N_{Th232}(\mathbf{r})\sigma_{Th232,\gamma}(E)\phi(E,\mathbf{r},t)dE dV \end{aligned} \]
The effective neutron source importance is another source of systematic error. A source neutron born near a high-importance region contributes more to blanket reactions than a neutron born where leakage dominates. Consequently, the same proton current can correspond to different fertile capture rates. Source efficiency must be validated rather than assumed.
In pulsed operation, delayed neutrons and isotope decay introduce multiple time constants. Prompt neutron decay may occur on microsecond to millisecond scales, while Th-233 and Pa-233 evolve over minutes and weeks. A complete data system must therefore combine fast digitizers for kinetics with long-term radiometric records for isotope buildup.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
5. Spectral treatment of neutron capture
The capture rate in ²³²Th is governed by resonance structure, moderation history, and self-shielding. Thermal values are useful reference points, but ADS blankets often include epithermal and fast components. A single thermal cross section therefore cannot characterize a spallation-driven thorium system. It can only serve as an anchor for comparison.
The correct starting point is the differential reaction-rate integral. In a finite target, this integral must be corrected for resonance self-shielding, flux depression, target thickness, and gamma self-attenuation. In powder, oxide, metal, or molten-salt geometries, the material form changes number density and attenuation corrections. These details are not secondary; they determine whether a measured activity is a faithful estimate of the reaction rate.
Spectrum weighting is particularly important for ²³³Pa and ²³³U because parasitic capture and fission probabilities vary with energy. A hard spectrum may reduce some thermal capture penalties but can modify fission utilization and neutron balance. A softened spectrum may improve certain capture rates while increasing sensitivity to resonance absorption and gamma backgrounds. The optimal spectrum is therefore an experimental variable, not a fixed property of thorium.
The general reaction-rate expression is:
\[ \begin{aligned} R_x(t)=N_x\int_0^{\infty}\sigma_x(E)\phi(E,t)dE\\ \langle\sigma_x\rangle_{\phi}=\frac{\int_0^{\infty}\sigma_x(E)\phi(E)dE}{\int_0^{\infty}\phi(E)dE}\\ R_x=N_x\langle\sigma_x\rangle_{\phi}\Phi \end{aligned} \]
Self-shielding in thorium-bearing samples is especially relevant in resonance regions. A thick sample can absorb neutrons at resonance energies near its surface, reducing the effective reaction rate in the interior. If the experiment is intended to validate library cross sections, the sample should be thin enough or the correction accurate enough to avoid confusing material geometry with nuclear data.
Cadmium-covered activation, multi-foil unfolding, and transport-calculated group spectra can be combined. Agreement among these methods strengthens the inference. Disagreement is not failure; it identifies the energy region where the experiment is most sensitive and where additional diagnostics are needed.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
6. Experimental architecture for ADS thorium characterization
A defensible experiment requires four coupled subsystems: a neutron source and subcritical assembly, a thorium-bearing sample or blanket region, detector and monitor systems, and a data-reduction chain. The neutron source may be a spallation target driven by protons, a deuteron-based neutron generator, or a benchmark source used for scaled subcritical experiments. The source choice determines the high-energy tail and the moderation demand.
The thorium sample should be designed so that reaction-rate interpretation is possible. Excessive thickness increases self-shielding and gamma self-absorption; excessive dilution reduces activity and worsens counting statistics. The target geometry should be simple enough for transport modeling yet representative enough to retain ADS-relevant spectral features. A layered design can include monitor foils upstream and downstream to measure local spectral gradients.
Detector selection should be tied to the observable. HPGe spectrometry is appropriate for activation products and decay curves. Fission chambers are appropriate for fission response and monitor normalization. Bonner spheres, proton recoil detectors, activation stacks, or time-of-flight diagnostics may be used for spectral information. No single detector determines the whole chain.
The experimental design should also define a cooling schedule before irradiation begins. Because Pa-233 evolves over weeks, a measurement campaign that only counts once after irradiation will generally be underdetermined. Multiple counting times allow the Pa component to be separated from short-lived backgrounds and from the delayed formation of U-233.
A recommended measurement matrix is summarized below.
| Element | Observable | Measurement route |
|---|---|---|
| Activation foil / target | ²³²Th(n,γ) rate | HPGe counting of Th-233/Pa-233 gamma lines after calibrated irradiation and decay intervals |
| Delayed gamma sequence | ²³³Pa buildup/decay | multi-time-point spectroscopy fitted to Bateman equations |
| Fission monitor | ²³³U(n,f) response | parallel-plate fission chamber or ionization chamber using monitor reaction normalization |
| Neutron spectrum | \(\phi(E)\) and lethargy profile | time-of-flight, activation unfolding, Bonner spheres, or validated transport calculation |
| Subcritical response | \(k_{eff}\), source efficiency, prompt decay constant | pulsed neutron source, source jerk, area ratio, or Rossi-alpha/PNS analyses |
7. Activation analysis for ²³²Th(n,γ)
Activation analysis connects a measured gamma peak to an irradiation-averaged reaction rate. The method is powerful because it can be performed with small samples and calibrated detectors, but its reliability depends on careful correction of irradiation, decay, and counting intervals. In ADS experiments, beam interruptions and source fluctuations must be treated explicitly.
The standard activation equation contains the number of target atoms, the spectrum-weighted reaction rate, gamma emission probability, full-energy peak efficiency, saturation correction, decay correction, and live-time correction. If the flux varies with time, the saturation factor must be replaced by a time integral over the recorded beam history. For Pa-233, the delayed feeding from Th-233 also requires chain-aware fitting rather than a single-nuclide activation equation.
The generic count equation for a single gamma line is:
\[ \begin{aligned} C_\gamma=\epsilon_\gamma I_\gamma A_0 e^{-\lambda t_d}\left(1-e^{-\lambda t_c}\right)\frac{t_{live}}{t_c}\\ A_0=N_T\langle\sigma\rangle\Phi\left(1-e^{-\lambda t_i}\right)\\ N_T=\frac{m f_{iso} N_A}{M} \end{aligned} \]
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
8. Protactinium-233 buildup and decay dynamics
Protactinium-233 is the kinetic bottleneck of the thorium breeding path. The half-life is long enough to preserve memory of the irradiation history and short enough that decay during multi-week cooling is experimentally visible. This dual character makes Pa-233 a sensitive diagnostic of the neutron field and a significant competitor in neutron economy.
The Pa-233 inventory after irradiation depends on production by Th-233 decay, loss by beta decay, and loss by neutron capture. In the simplest constant-flux model, the effective removal constant is the sum of the radioactive decay constant and the spectrum-weighted capture rate. This formulation immediately produces a dimensionless competition parameter, which allows experiments at different flux levels to be compared on the same scale.
The competition parameter is:
\[ \begin{aligned} r_{Pa}=\frac{\langle\sigma_{Pa233,\gamma}\rangle\Phi}{\lambda_{Pa233}}\\ \lambda_{Pa,eff}=\lambda_{Pa233}+\langle\sigma_{Pa233,\gamma}\rangle\Phi=\lambda_{Pa233}(1+r_{Pa})\\ f_{capture,Pa}=\frac{r_{Pa}}{1+r_{Pa}} \end{aligned} \]
Radiochemical Pa separation, if used, changes the mathematical boundary conditions. Removal of Pa from the neutron field suppresses capture loss but introduces separation time, recovery fraction, and daughter-growth corrections. The kinetic benefit must be weighed against the added uncertainty and operational complexity.
The best Pa-233 analysis uses multiple gamma lines and multiple counting times. Fitting a single line at one time point can mask interference, efficiency error, or feeding from a precursor. Multi-line, multi-time fitting allows residuals to reveal whether the assumed chain model is adequate.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
9. Uranium-233 formation and utilization
The final objective of the breeding chain is not merely the appearance of U-233 atoms but the formation of a usable fissile inventory under specified irradiation and fuel-cycle constraints. U-233 itself undergoes fission and capture in the same neutron field. Consequently, late-time irradiation can both create and consume the desired product. This duality must be included in any yield definition.
A useful experimental metric is the net U-233 formation rate after subtracting destruction terms. In a low-power activation test, destruction of U-233 may be negligible because its initial concentration is small and the fluence is limited. In an ADS blanket operating at substantial flux, however, U-233 destruction can become part of the neutron economy. The model must therefore state whether U-233 is treated as a passive product or an active fissile participant.
The net balance can be written as:
\[ \begin{aligned} \frac{dN_{U233}}{dt}=\lambda_{Pa233}N_{Pa233}-\Phi\left(\langle\sigma_{U233,f}\rangle+\langle\sigma_{U233,\gamma}\rangle\right)N_{U233}\\ Y_{U233}(t)=\frac{N_{U233}(t)}{\int_0^t R_{Th232,\gamma}(\tau)d\tau} \end{aligned} \]
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
10. Time-of-flight and fission benchmark relevance
Time-of-flight measurements are essential for differential cross-section information in broad spectra. The CERN n_TOF results for ²³²Th and ²³³U fission up to high energy demonstrate the type of evidence needed when ADS spectra extend beyond conventional thermal or low-fast regimes. These data do not directly determine Pa-233 buildup, but they constrain the fission and threshold-reaction components of the same physical environment.
For activation experiments, time-of-flight data play an indirect but important role. They improve evaluated libraries and reduce ambiguity in transport calculations. If the transport model predicts that a significant part of the neutron field lies in an energy range where evaluations diverge, experimental reporting must identify that sensitivity. A one-group number without spectral context is insufficient.
Parallel-plate avalanche counters, fission chambers, and time-projection chambers provide independent ways to monitor fission response. Their value in a thorium ADS study is not limited to fission yield; they also support flux normalization, source characterization, and validation of neutron-energy reconstruction. The detector suite should therefore be selected for the combined capture-fission problem.
The relationship between microscopic data and an integral ADS observable is expressed through spectrum folding:
\[ \begin{aligned} \bar{\sigma}_{j}^{exp}=\frac{\int_{E_1}^{E_2}\sigma_j(E)\phi(E)dE}{\int_{E_1}^{E_2}\phi(E)dE}\\ O_k^{calc}=\sum_g R_{kg}\sigma_g, \qquad O_k^{exp}=O_k^{calc}+\varepsilon_k \end{aligned} \]
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
11. Uncertainty propagation and metrological requirements
The credibility of a thorium ADS characterization depends on uncertainty propagation. Statistical counting uncertainty is often the smallest visible component, while efficiency calibration, sample mass, gamma emission probabilities, flux normalization, spectrum unfolding, and nuclear-data covariance can dominate. Reporting only a peak-area uncertainty understates the error budget.
The Guide to the Expression of Uncertainty in Measurement (GUM) provides the general framework: define input quantities, assign standard uncertainties, evaluate sensitivities, and combine correlated terms using covariance matrices. Nuclear data introduce additional structure because cross sections in neighboring energy groups may be correlated. Treating every group as independent can either overstate or understate the final uncertainty.
A generic first-order propagation expression is:
\[ \begin{aligned} u_c^2(y)=\sum_i\left(\frac{\partial f}{\partial x_i}\right)^2u^2(x_i)+2\sum_{i<j}\frac{\partial f}{\partial x_i}\frac{\partial f}{\partial x_j}u(x_i,x_j)\\ \mathbf{C}_y=\mathbf{S}\mathbf{C}_x\mathbf{S}^{T} \end{aligned} \]
A practical covariance workflow stores both experimental and nuclear-data contributions. Experimental covariance includes detector calibration correlations, monitor reaction normalization, and sample mass. Nuclear-data covariance includes group-wise cross-section correlations and decay-data uncertainties. Their separation allows future library updates without repeating the entire experiment.
Sensitivity coefficients should be reported alongside final uncertainties. They show whether the experiment is limited by detector efficiency, neutron spectrum, Pa capture data, or source normalization. Such information is more useful for future experiment design than a single combined uncertainty number.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
| Component | Mechanism | Control method |
|---|---|---|
| Reference nuclear data | cross-section and covariance library selection | compare ENDF/B, JEFF, JENDL, TENDL; propagate covariance where available |
| Flux normalization | monitor reaction and detector calibration | co-irradiate Au/In/Ni/Al monitors; repeat cadmium-ratio or spectrum checks |
| Gamma spectrometry | efficiency curve, peak area, summing, dead time | calibrated sources, geometry replication, Monte Carlo efficiency verification |
| Irradiation history | beam current integration and source interruptions | high-frequency beam-current logging and explicit time-bin integration |
| Chemical separation | Pa recovery factor and timing delay | tracer-based recovery, duplicate aliquots, decay-corrected blank controls |
| Transport model | geometry, material, self-shielding, leakage | benchmark against YALINA/GUINEVERE/MYRRHA-class integral experiments |
12. Analytical solution of the Pa-233 chain
For constant flux and constant source conditions, the coupled chain admits closed-form solutions. These expressions are valuable because they expose the dependence of Pa buildup on irradiation time and on the competition between decay and neutron capture. Numerical transport calculations remain necessary for realistic geometries, but the analytical solution supplies a transparent benchmark for checking code output.
If ²³³Th is treated explicitly, Pa-233 production is delayed by the Th-233 decay constant. If Th-233 is eliminated under the quasi-steady approximation, Pa-233 production becomes approximately equal to the thorium capture rate after the initial transient. The difference matters in short irradiations but becomes small for campaigns lasting many Th-233 half-lives.
The two-stage approximation for Pa buildup is:
\[ \begin{aligned} N_{Pa}(t)\approx\frac{R_{Th}}{\lambda_{Pa}+k_{Pa}}\left(1-e^{-(\lambda_{Pa}+k_{Pa})t}\right)\\ N_{U}(t)\approx\lambda_{Pa}\int_0^t N_{Pa}(\tau)e^{-k_U(t-\tau)}d\tau\\ k_{Pa}=\langle\sigma_{Pa233,\gamma}\rangle\Phi, \qquad k_U=\Phi(\langle\sigma_{U233,f}\rangle+\langle\sigma_{U233,\gamma}\rangle) \end{aligned} \]
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
13. Dimensionless regimes and design interpretation
Dimensionless ratios allow different ADS experiments to be compared even when source strength, geometry, or spectrum differs. The Pa competition parameter \(r_{Pa}\) compares neutron capture to beta decay. A source-coupling parameter can compare local fertile capture to nominal beam current. A spectral hardness parameter can compare fast-to-thermal reaction rates. Together they provide a compact description of the breeding environment.
The most important result is that Pa-233 is not a fixed loss; it is a regime-dependent loss. At low flux, Pa decays before capture. At high flux or in a resonance-enhanced spectrum, Pa capture can divert material away from U-233. The same isotope can therefore be nearly harmless in one experiment and central in another. Experimental characterization must identify the regime before drawing fuel-cycle conclusions.
The regimes are summarized below.
| Regime | Criterion | Interpretation |
|---|---|---|
| Low flux, decay dominated | \(\sigma_{Pa,\gamma}\Phi \ll \lambda_{Pa}\) | Pa decays to U-233 before significant parasitic capture |
| Balanced competition | \(\sigma_{Pa,\gamma}\Phi \approx \lambda_{Pa}\) | measurable production penalty; Pa removal or spectrum hardening becomes important |
| Capture dominated | \(\sigma_{Pa,\gamma}\Phi \gg \lambda_{Pa}\) | Pa inventory is depleted by neutron capture; in-situ breeding efficiency decreases |
| Pulsed source transient | \(S(t)\) not constant | activation integrals must use time-binned flux and pulsed-neutron kinetics |
14. Numerical demonstration using open-data parameters
A numerical demonstration can be built from representative half-lives and library-derived reaction-rate coefficients. Because this article does not claim a proprietary irradiation result, the numerical values should be read as transparent calculations rather than measurements. The purpose is to show how the equations convert measured or evaluated inputs into Pa and U inventory predictions.
Consider a homogeneous thorium region with a constant spectrum-weighted capture rate coefficient \(k_{Th}=\langle\sigma_{Th,\gamma}\rangle\Phi\). The Pa removal coefficient is \(\lambda_{Pa}+\langle\sigma_{Pa,\gamma}\rangle\Phi\). If \(r_{Pa}=0.1\), the capture penalty is small but visible. If \(r_{Pa}=1\), about half of the Pa removals occur through neutron capture rather than beta decay in the simplified model. If \(r_{Pa}=2\), neutron capture dominates the Pa loss channel.
The simplified loss fraction is:
\[ F_{loss,Pa}=\frac{\int_0^t k_{Pa}N_{Pa}(\tau)d\tau}{\int_0^t (k_{Pa}+\lambda_{Pa})N_{Pa}(\tau)d\tau}\approx\frac{r_{Pa}}{1+r_{Pa}} \]
The demonstration also shows the danger of quoting a bred U-233 amount without the irradiation timeline. The same integrated fluence delivered in one continuous interval or in separated cycles can produce different Pa and U inventories because Pa decays between cycles. ADS experiments must preserve the time sequence, not just the total fluence.
A benchmark-quality calculation should include at least three library runs. If the predicted Pa loss fraction or U-233 yield changes materially across ENDF/B, JEFF, and JENDL processing, the result should be reported as library-sensitive. This is not a weakness; it accurately communicates the present state of knowledge.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
| \(r_{Pa}\) | capture fraction | experimental meaning |
|---|---|---|
| 0.01 | 0.0099 | decay-dominated; capture negligible |
| 0.1 | 0.091 | small but measurable penalty |
| 1 | 0.500 | balanced decay and capture |
| 2 | 0.667 | capture-dominated Pa loss |
| 10 | 0.909 | severe Pa protection problem |
15. Experimental protocol and data-reduction sequence
The recommended protocol begins with pre-irradiation characterization: sample mass, isotopic composition, density, geometry, impurity assay, and detector-calibration geometry. The second step is flux characterization using monitor reactions and transport modeling. The third step is irradiation with recorded beam current, source pulse structure, temperature, and sample position. The fourth step is a multi-time-point counting program designed to resolve Th-233, Pa-233, and longer-lived backgrounds.
Data reduction should proceed from raw spectra to peak areas, from peak areas to activities, from activities to reaction rates, and from reaction rates to kinetic parameters. Each conversion must retain uncertainty and covariance information. A final inventory number without the intermediate chain is not sufficient for reproducibility.
The sequence is best expressed as a traceability chain:
\[ \mathrm{raw\ spectrum}\rightarrow C_\gamma\rightarrow A(t)\rightarrow R_{n,\gamma}\rightarrow N_{Pa}(t)\rightarrow N_{U}(t)\rightarrow Y_{U233}\pm u_c \]
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
16. Safety, safeguards, and radiological controls
Thorium experiments involving U-233 production require safety and safeguards awareness even at small scale. U-233 is fissile, long-lived, and often accompanied by U-232 impurities in practical fuel-cycle contexts. U-232 daughters emit intense gamma radiation, affecting handling, shielding, and accountancy. A kinetics article should not ignore these constraints because they influence feasible sampling, cooling, separation, and measurement procedures.
Protactinium chemistry is also non-trivial. If Pa separation is used to reduce neutron capture losses, the experiment becomes a radiochemical system rather than a purely physical activation measurement. Recovery fraction, cross-contamination, hold-up, and decay during processing become part of the kinetic model. Any Pa-management claim must therefore include chemical-yield data and timing uncertainty.
From a safeguards perspective, reporting should clearly separate scientific reaction-rate characterization from any claim of material production capability. The experimental objective is to measure kinetics, validate models, and understand neutron economy. The quantities required for that objective are microgram-to-milligram sample inventories, reaction rates, and detector responses, not weapons-relevant processing information.
Radiological controls should be integrated into the measurement plan. Counting geometry, source transfer, decay storage, glovebox operations, and waste classification can all affect the time schedule. Because the Pa half-life is central to the physics, administrative delays and separation times are not merely logistical; they change the measured inventory.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
17. Discussion
The characterization framework reveals why thorium ADS analysis must be kinetic rather than static. A static conversion ratio may describe an asymptotic design calculation, but it cannot determine what an experiment actually measures after a specified irradiation and cooling sequence. Pa-233 stores time-history information, and the measured gamma activity is a convolution of production, decay, capture, and detector response.
The framework also clarifies the role of evaluated nuclear data. Libraries are indispensable, but they are not substitutes for experimental characterization. Their uncertainties and disagreements become part of the analysis, especially in epithermal and fast energy regions relevant to ADS. A high-quality article should therefore report which libraries were used, how group constants were processed, and whether covariance data were propagated.
Integral benchmarks are equally important. Experiments such as MUSE, YALINA, and GUINEVERE established practical methods for subcriticality measurement and source-driven kinetics. A thorium-specific experiment can build on those methods while adding the chemistry and decay-chain complexity of Pa-233. This connection prevents the study from becoming an isolated activation exercise.
A practical design implication is that Pa protection strategies must be evaluated under the actual spectrum and flux of the system. Chemical separation, online removal, spectrum adjustment, or irradiation-cooling cycling can all reduce parasitic Pa capture in principle. Whether any of them is justified depends on the ratio \(r_{Pa}\), the cost and uncertainty of separation, and the desired U-233 production metric.
The most defensible reporting format is a layered one: raw observables, corrected activities, reaction-rate coefficients, isotope inventories, and model-dependent extrapolations should be kept distinct. This structure allows other researchers to reprocess the data when nuclear libraries are updated. It also avoids overstating conclusions from a single evaluated file or a single detector channel.
The framework is also compatible with molten-salt or solid-fuel variants. The equations are identical, but the meaning of Pa management differs. In solid targets, Pa remains largely in situ unless post-irradiation chemistry is performed. In fluid systems, online removal is conceptually possible, but then hydrodynamic residence time and chemical separation kinetics enter the same rate network.
The principal experimental challenge is not lack of equations but lack of fully reported coupling between equations and observables. Every coefficient in the model must be linked to a measurement, a processed library value, or a justified assumption. The article's proposed format is intended to enforce that linkage.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
18. Conclusions
The ²³²Th → ²³³U path in an accelerator-driven system is governed by coupled neutron-capture and radioactive-decay kinetics. The presence of ²³³Pa makes the problem intrinsically time-dependent because Pa decay competes with neutron capture over experimentally relevant timescales.
A valid characterization requires more than a thorium capture cross section. It requires neutron-spectrum measurement or validation, beam-history integration, activation analysis, Pa-specific decay fitting, U-233 formation accounting, subcritical source-coupling diagnostics, and covariance-aware uncertainty propagation.
The most useful governing ratio is \(r_{Pa}=\langle\sigma_{Pa,\gamma}\rangle\Phi/\lambda_{Pa}\). When this ratio is much smaller than one, Pa decay dominates and U-233 formation is efficient. When it approaches or exceeds one, parasitic capture materially reduces the bred fissile yield and Pa management becomes a central design issue.
The article establishes a reproducible structure for experimental ADS thorium studies without relying on unsupported data claims. The method is compatible with public evaluated libraries, EXFOR experimental data, JANIS inspection, and integral ADS benchmarks. Its purpose is to make the inferred kinetic parameters traceable, testable, and updateable.
Future work should prioritize Pa-233 capture data under ADS-relevant spectra, coupled activation and transport benchmarks in thorium-bearing subcritical assemblies, and systematic propagation of nuclear-data covariance into bred U-233 inventory predictions.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix A. Core LaTeX equation set
The following equations are repeated in compact form for direct reuse in data-reduction sheets and transport-code validation notes. Symbols should be defined with the same units throughout a calculation: number density in atoms per cubic centimeter, flux in neutrons per square centimeter per second, microscopic cross section in square centimeters or barns after conversion, and decay constants in inverse seconds.
\[ \begin{aligned} \lambda=\frac{\ln 2}{T_{1/2}}\\ R=N\int\sigma(E)\phi(E)dE\\ \frac{d\mathbf{N}}{dt}=\mathbf{A}\mathbf{N}+\mathbf{q}(t)\\ \mathbf{N}(t)=e^{\mathbf{A}t}\mathbf{N}(0)+\int_0^t e^{\mathbf{A}(t-\tau)}\mathbf{q}(\tau)d\tau\\ r_{Pa}=\frac{\langle\sigma_{Pa,\gamma}\rangle\Phi}{\lambda_{Pa}}\\ f_{capture,Pa}=\frac{r_{Pa}}{1+r_{Pa}} \end{aligned} \]
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix B. Reporting checklist
A complete ADS thorium characterization report should include: sample composition; irradiation geometry; material density; detector efficiency calibration; neutron spectrum or group constants; monitor reactions; evaluated library versions; beam-current history; irradiation, decay, and count times; gamma lines used; chemical recovery factors if applicable; all uncertainty components; correlation assumptions; transport-code version; and benchmark validation cases.
The minimum data package should permit an independent reader to recalculate the reaction-rate coefficients and reproduce the Pa-233 buildup curve. If a result cannot be recalculated from the reported information, it should be treated as a qualitative observation rather than a quantitative kinetic characterization.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix C. Neutron source normalization and beam-history reconstruction
A source-driven thorium experiment must reconstruct the neutron source term as a function of time. The practical input is normally beam current, target status, accelerator trip history, and monitor response. The physically relevant quantity is not the electrical current alone; it is the number of source neutrons born in locations and energy ranges capable of contributing to fertile capture. Beam-current integration should therefore be converted into a source term only after target yield, transport coupling, and detector normalization have been checked.
Beam interruptions matter even when the integrated charge is unchanged. Isotope production and decay are not linear in calendar time when half-lives are comparable with the irradiation cycle. A ten-day continuous irradiation and ten one-day irradiations separated by cooling intervals can produce different Pa-233 and U-233 inventories. The data system should keep a time-bin representation fine enough to resolve source changes and long enough to cover the decay follow-up.
Source normalization can be cross-checked using monitor foils and prompt neutron detectors. The purpose of the monitor is not merely to scale a calculation but to reveal whether the assumed source-to-blanket coupling is stable. If monitor ratios drift during the run, the reaction-rate model should include separate irradiation intervals with interval-specific flux spectra.
A recommended source-history file contains time, beam current, accelerator status, target temperature, detector live fraction, monitor counts, and any changes in sample position. This file should be treated as primary experimental data. Without it, Pa-233 decay fitting can be mathematically elegant but physically ambiguous.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix D. Gamma-line selection and spectral interference control
HPGe gamma spectrometry is a central diagnostic for Th-233 and Pa-233 activity, but line selection must be made with attention to interference. A useful line should have sufficient emission probability, minimal overlap with fission products or activation contaminants, and a detector efficiency that can be calibrated in the same geometry. Peak area extraction should state the background model and integration limits.
Counting geometry is part of the measurement equation. If the irradiated thorium sample is counted in a different geometry from the calibration standard, coincidence summing, self-attenuation, and extended-source efficiency corrections may be needed. Monte Carlo efficiency modeling can support the correction, but it should be benchmarked against calibration sources or multi-line consistency checks.
Pa-233 analysis benefits from repeated spectra over several half-lives. A correct decay curve should reproduce the same activity when fitted from different gamma lines after all corrections. If line-specific activities diverge, the experiment should investigate interference, dead time, pile-up, or an incorrect decay-chain assumption before calculating U-233 production.
The final spectrum archive should include raw spectra, live time, real time, calibration files, peak fits, residuals, background spectra, and uncertainty components. Publishing only corrected activities removes information needed to identify hidden biases.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix E. Self-shielding, resonance treatment, and sample geometry
Thorium has resolved and unresolved resonance structure that can make sample geometry a dominant correction. A thin sample approximates an unperturbed flux measurement; a thick sample perturbs the neutron field and sees a spectrum different from the surrounding medium. For this reason, activation target thickness must be selected in relation to the energy range of interest and to the expected resonance absorption.
Self-shielding correction should be calculated with the same evaluated library used for the reaction-rate inference. If the correction is large, the result becomes sensitive to the very nuclear data under test. A robust experiment therefore uses multiple sample thicknesses or diluted samples to separate geometry effects from cross-section effects.
In a multi-group code, resonance treatment may use Bondarenko factors, probability tables, subgroup methods, or continuous-energy Monte Carlo. The selected approach should be reported. For Pa-233, data scarcity can make the uncertainty of resonance treatment more important than for stable monitor materials.
Geometry uncertainty also affects gamma spectrometry. The emitting volume after irradiation may not be identical to the assumed volume if chemical processing, powder redistribution, or encapsulation changes the source shape. This effect should be included in the efficiency uncertainty if the activity is used for quantitative kinetics.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix F. Library comparison and covariance workflow
A single evaluated library should never be the only quantitative basis for an ADS thorium conclusion. The article recommends processing at least three independent libraries, for example ENDF/B, JEFF, and JENDL, under the same spectrum and geometry. Differences in predicted Pa loss fraction, Th capture rate, and U-233 destruction rate should be reported as library sensitivity.
Covariance processing begins by mapping energy-dependent cross-section uncertainties into group constants. If covariance data are absent for a reaction, a documented surrogate uncertainty should be assigned rather than silently treating the reaction as exact. This is especially relevant for Pa-233 capture because it can control the bred U-233 yield in high-flux regimes.
The covariance matrix should be propagated to the output inventory through sensitivity coefficients. When the model is nonlinear, finite-difference perturbation, adjoint sensitivity, or Monte Carlo sampling can be used. The output should state which reactions dominate uncertainty. This information directs future measurement priorities.
Library comparison also protects against false precision. If the experimental uncertainty is 3 percent but library-to-library variation is 15 percent in the relevant spectral window, the defensible conclusion must reflect the larger nuclear-data limitation.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix G. Subcriticality diagnostics and source-jerk interpretation
Subcriticality determines how source neutrons are multiplied before contributing to fertile capture. Techniques such as pulsed-neutron source analysis, area-ratio methods, source jerk, and Rossi-alpha measurements estimate reactivity or prompt decay constants. Each method has assumptions about detector position, delayed neutron contribution, and spatial mode effects.
In a thorium breeding experiment, subcriticality diagnostics serve two purposes. First, they ensure safe and reproducible operation of the source-driven assembly. Second, they constrain the neutron field used in activation analysis. If \(k_{eff}\) or source efficiency changes during a campaign, the isotope-production model must include that change.
Detector placement is critical. A detector near the source may see a different time shape from a detector in the fertile blanket. Modal effects can cause local prompt decay constants to differ from global kinetics estimates. Benchmark experiments such as MUSE and GUINEVERE demonstrate why multiple detector locations are useful.
The source-jerk method can provide an operational reactivity indicator, but it should be calibrated against transport or benchmark results. Treating it as an absolute truth without spatial correction can distort the inferred relationship between beam power and thorium capture rate.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix H. Radiochemical Pa management and timing corrections
If Pa-233 is chemically separated from thorium, the kinetics change from an in-situ decay problem to a coupled irradiation-separation-decay problem. The separation time, recovery yield, residual thorium contamination, and Pa hold-up in equipment must be measured. The kinetic model then includes a transfer term from the neutron field to a shielded decay reservoir.
A simple separation model can be represented by an extraction efficiency and a residence-time distribution. Perfect instantaneous removal is a mathematical ideal that real systems do not achieve. If the objective is to evaluate Pa protection, the experiment should measure the actual fraction of Pa removed before significant additional neutron capture occurs.
Timing corrections are unavoidable. Pa decays during dissolution, separation, washing, transfer, and counting. Each step has a start time, end time, and recovery factor. These times must be included in the Bateman model rather than applied as a single post hoc correction.
Radiochemical blanks, tracers, and duplicate samples are needed to detect cross-contamination and incomplete recovery. Without these controls, a claimed Pa-management improvement may reflect chemistry rather than neutron kinetics.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix I. Transport-code verification and benchmark hierarchy
Transport calculations convert geometry and nuclear data into local spectra and reaction rates. Continuous-energy Monte Carlo codes such as MCNP, Serpent, and OpenMC are common choices; deterministic codes can also be used with appropriate group libraries. The code name is less important than verification of the modeling chain.
A benchmark hierarchy should start with simple cases: bare monitor foil, known activation spectrum, and isolated sample self-shielding. It should then progress to subcritical benchmark assemblies, and only then to the full ADS thorium configuration. This hierarchy helps identify which discrepancy belongs to detector calibration, nuclear data, geometry, or source coupling.
Model input files should be archived with material compositions, densities, temperatures, source definitions, tally definitions, variance-reduction settings, and random seeds where relevant. Reaction-rate tallies should be normalized in a way that is consistent with the experimental source term.
Benchmark agreement should be evaluated using quantitative residuals, not visual similarity alone. The residual vector should include monitor reaction rates, detector time constants, spectrum indicators, and, where possible, Pa-233 activity evolution.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix J. Worked symbolic reduction for a constant-flux experiment
For a constant-flux demonstration, define \(R_{Th}=N_{Th232}\langle\sigma_{Th232,\gamma}\rangle\Phi\). If Th-233 is treated as instantaneous relative to the irradiation time, Pa production is approximately \(R_{Th}\). The Pa inventory then approaches a steady value determined by the effective removal constant \(\lambda_{Pa}+k_{Pa}\).
The initial slope of the Pa curve measures the thorium capture rate, while the late-time curvature measures Pa removal. A multi-time measurement can therefore infer both production and removal if the counting precision and time spacing are adequate. Early points should be dense enough to see buildup; late points should extend over a meaningful fraction of the Pa half-life.
U-233 formation is obtained by integrating the Pa decay term. If U-233 destruction is negligible, the U inventory is the time integral of \(\lambda_{Pa}N_{Pa}\). If destruction is not negligible, the integral must include the U-233 fission and capture removal coefficient. The distinction should be stated explicitly.
This symbolic reduction is not a replacement for transport calculation. Its value is diagnostic: if a complex code violates the limiting behavior of the symbolic model under the same simplifying assumptions, the code setup or normalization should be rechecked.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix K. Reporting format for a benchmark-quality paper
A benchmark-quality article should publish enough information for independent recalculation. The minimum package includes the sample specification, irradiation position, source history, monitor reaction data, detector calibration, gamma-line list, raw or processed spectra, nuclear library versions, processing code, transport model, and uncertainty budget.
Results should be layered. The first layer is raw observation, such as peak counts and monitor counts. The second layer is corrected activity. The third layer is reaction rate. The fourth layer is inferred inventory. The fifth layer is design extrapolation. Mixing these layers makes the result difficult to audit.
Figures should show both data and model residuals. Tables should include units and uncertainty. Equations should be written in a portable form so that symbols are not corrupted when the file is exchanged. The revised documents therefore use explicit LaTeX blocks for every mathematical expression.
The references should include nuclear-data sources, ADS benchmark literature, detector methodology, uncertainty standards, and fuel-cycle assessments. A narrow bibliography is not adequate for a multidisciplinary ADS thorium study.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
Appendix L. Limitations and future measurement priorities
The framework is limited by the availability and quality of Pa-233 capture data in ADS-relevant spectra. Improved differential measurements, especially in energy regions that dominate realistic source-weighted reaction rates, would directly reduce uncertainty in U-233 breeding predictions.
A second limitation is the scarcity of thorium-bearing integral ADS benchmarks with full Pa time-series measurements. Existing subcritical benchmarks constrain kinetics and source coupling, but dedicated thorium activation and Pa decay campaigns would provide stronger validation.
Future experiments should combine in-situ monitors, post-irradiation spectroscopy, and transport-code covariance analysis. The highest value measurements are those that reduce ambiguity between source strength, spectral shape, and nuclear-data uncertainty.
The final priority is reporting discipline. Data that cannot be reprocessed when libraries change will lose scientific value. The proposed protocol treats traceability as part of the physics rather than as an administrative afterthought.
For experimental reproducibility, every parameter in this section should be reported with units, acquisition time, calibration source, and processing path. The result is not a single number but a traceable inference chain that can be audited when new cross-section evaluations are released.
The interpretation should also include a negative result criterion. If two independent diagnostics disagree outside the combined uncertainty, the discrepancy should be analyzed before any breeding-yield conclusion is drawn. In ADS physics, inconsistency often identifies source-coupling or spectral errors rather than detector failure.
References
[1] International Atomic Energy Agency, Thorium Fuel Cycle: Potential Benefits and Challenges, IAEA-TECDOC-1450, Vienna, 2005.
[2] International Atomic Energy Agency, Evaluated Nuclear Data for Nuclides within the Thorium-Uranium Fuel Cycle, STI/PUB/1435, Vienna, 2010.
[3] OECD Nuclear Energy Agency, Introduction of Thorium in the Nuclear Fuel Cycle: Short- to Long-Term Considerations, NEA No. 7224, Paris, 2015.
[4] OECD Nuclear Energy Agency, Technology and Components of Accelerator-Driven Systems, NEA No. 6897, Paris, 2011.
[5] OECD Nuclear Energy Agency, Accelerator-Driven Systems and Fast Reactors in Advanced Nuclear Fuel Cycles, NEA, Paris, 2002.
[6] International Atomic Energy Agency, Use of Low Enriched Uranium Fuel in Accelerator Driven Subcritical Systems, IAEA-TECDOC-1821, Vienna, 2017.
[7] International Atomic Energy Agency, LiveChart of the Nuclides, Nuclear Data Section, Vienna, accessed 2026.
[8] National Nuclear Data Center, NuDat and ENSDF Nuclear Structure and Decay Data, Brookhaven National Laboratory, accessed 2026.
[9] National Nuclear Data Center, ENDF/B-VIII.1 Evaluated Nuclear Data Library, Brookhaven National Laboratory, release information 2024.
[10] D. A. Brown et al., ENDF/B-VIII.0: The 8th Major Release of the Nuclear Reaction Data Library, Nuclear Data Sheets 148, 1-142, 2018.
[11] A. J. Koning, D. Rochman, Modern Nuclear Data Evaluation with the TALYS Code System, Nuclear Data Sheets 113, 2841-2934, 2012.
[12] O. Iwamoto et al., Japanese Evaluated Nuclear Data Library Version 5: JENDL-5, Journal of Nuclear Science and Technology 60, 1-60, 2023.
[13] A. J. M. Plompen et al., The Joint Evaluated Fission and Fusion Nuclear Data Library, JEFF-3.3, European Physical Journal A 56, 181, 2020.
[14] A. J. Koning, S. Hilaire, S. Goriely, TENDL: Complete Nuclear Data Library for Innovative Nuclear Science and Technology, Nuclear Data Sheets and TENDL documentation.
[15] N. Otuka et al., Towards a More Complete and Accurate Experimental Nuclear Reaction Data Library, Nuclear Data Sheets 120, 272-276, 2014.
[16] V. V. Zerkin and B. Pritychenko, The Experimental Nuclear Reaction Data (EXFOR), Nuclear Data Sheets 120, 272-276 and related IAEA/NNDC documentation, 2018.
[17] OECD Nuclear Energy Agency, JANIS: Java-based Nuclear Information Software, Nuclear Data Services documentation.
[18] International Atomic Energy Agency, EXFOR: Experimental Nuclear Reaction Data, Nuclear Data Section documentation.
[19] SCK CEN, MYRRHA: Multi-purpose Hybrid Research Reactor for High-tech Applications, project infrastructure documentation, accessed 2026.
[20] International Atomic Energy Agency, MYRRHA: An Accelerator Driven System to Manage Radioactive Waste, IAEA Bulletin, 2024.
[21] A. Billebaud et al., The GUINEVERE Project for Accelerator Driven System Physics, GLOBAL 2009, Paris, 2009.
[22] C. Rubbia et al., Conceptual Design of a Fast Neutron Operated High Power Energy Amplifier, CERN/AT/95-44(ET), 1995.
[23] D. Villamarin et al., First Measurements of the Kinetic Response of the MUSE-4 Experiment, OECD/NEA Information Exchange Meeting, 2001.
[24] G. Aliberti et al., MUSE-4 Experiment Measurements and Analysis, Argonne National Laboratory Report, 2003.
[25] Y. Gohar et al., YALINA Facility: A Subcritical Accelerator-Driven System for Nuclear Energy Research, Argonne National Laboratory Report, 2010.
[26] C. M. Persson et al., Pulsed Neutron Source Experiments in the YALINA-Booster Subcritical Assembly, Annals of Nuclear Energy 35, 2357-2364, 2008.
[27] A. Fokau, Accelerator-Driven Systems: Source Efficiency and Subcritical Reactivity Determination, Ph.D. thesis, KTH Royal Institute of Technology, 2010.
[28] C. Berglof, On Measurement and Monitoring of Reactivity in Subcritical Reactor Systems, Ph.D. thesis, KTH Royal Institute of Technology, 2010.
[29] N. Marie et al., Reactivity Monitoring of the Accelerator Driven VENUS-F Reactor, Annals of Nuclear Energy and related GUINEVERE publications, 2019.
[30] D. Tarrio et al., Neutron-Induced Fission Cross Sections of 232Th and 233U up to 1 GeV Using Parallel Plate Avalanche Counters at the CERN n_TOF Facility, Physical Review C 107, 044616, 2023.
[31] N. Otuka, K. Tada, O. Cabellos, O. Iwamoto, Evaluation of Uranium-233 Neutron Capture Cross Section in the keV Region, arXiv:2408.10098 and related evaluation work, 2024.
[32] C.-L. Lan et al., Determination of the Neutron Capture Cross Sections of 232Th at 14.1 MeV and 14.8 MeV Using the Neutron Activation Method, arXiv:1608.00102, 2016.
[33] G. de Saussure, Evaluation of the 232Th Neutron Capture Cross Section from 3 keV to 20 MeV, Oak Ridge National Laboratory memorandum, 1977.
[34] G. G. Myasishcheva et al., Measurement of the Effective Cross Section in Thorium for Thermal Neutrons, Soviet Atomic Energy 2, 1957.
[35] R. L. Macklin and J. H. Gibbons, Neutron Capture Cross Sections of Thorium and Uranium Nuclides, Physical Review and ORNL capture-measurement literature.
[36] W. M. Moore, Effective Capture Cross Section of Pa-233 for Thermal Reactor Neutrons, Nuclear Science and Engineering 1, 1956.
[37] V. Kvasov et al., Thermal Neutron Capture Cross Section of Protactinium Isotopes, Progress in Nuclear Science and Technology 2, 2002.
[38] S. F. Mughabghab, Atlas of Neutron Resonances, 6th ed., Elsevier, 2018.
[39] J. Kopecky, Atlas of Neutron Capture Cross Sections, IAEA/INDC and ECN report, 1997.
[40] R. Capote et al., RIPL - Reference Input Parameter Library for Calculation of Nuclear Reactions and Nuclear Data Evaluations, Nuclear Data Sheets 110, 3107-3214, 2009.
[41] R. E. MacFarlane et al., The NJOY Nuclear Data Processing System, Los Alamos National Laboratory documentation and Nuclear Data Sheets 2012 update.
[42] X-5 Monte Carlo Team, MCNP - A General Monte Carlo N-Particle Transport Code, Los Alamos National Laboratory, LA-UR-03-1987 and later manuals.
[43] J. Leppanen et al., The Serpent Monte Carlo Code: Status, Development and Applications, Annals of Nuclear Energy 82, 142-150, 2015.
[44] P. K. Romano et al., OpenMC: A State-of-the-Art Monte Carlo Code for Research and Development, Annals of Nuclear Energy 82, 90-97, 2015.
[45] S. Agostinelli et al., GEANT4 - A Simulation Toolkit, Nuclear Instruments and Methods in Physics Research A 506, 250-303, 2003.
[46] J. Allison et al., Recent Developments in GEANT4, Nuclear Instruments and Methods in Physics Research A 835, 186-225, 2016.
[47] E. Pirovano et al., A Detector System for Absolute Measurements of Fission Cross Sections at n_TOF Below 200 MeV, European Physical Journal A and arXiv:2211.13323, 2022.
[48] C. Paradela et al., High Accuracy Determination of the 238U/235U Fission Cross Section Ratio up to about 1 GeV at n_TOF, Physical Review C and arXiv:1410.7737, 2014.
[49] R. J. Casperson et al., Measurement of the 238U(n,f)/235U(n,f) Cross Section Ratio with the NIFFTE Fission Time Projection Chamber, Physical Review C 97, 034618, 2018.
[50] M. Salvatores, Integral Experiments and Their Use for Nuclear Data Validation, European Physical Journal Nuclear Sciences & Technologies 7, 2021.
[51] J. R. Lamarsh and A. J. Baratta, Introduction to Nuclear Engineering, 3rd ed., Prentice Hall, 2001.
[52] J. J. Duderstadt and L. J. Hamilton, Nuclear Reactor Analysis, Wiley, 1976.
[53] W. M. Stacey, Nuclear Reactor Physics, 3rd ed., Wiley-VCH, 2018.
[54] G. I. Bell and S. Glasstone, Nuclear Reactor Theory, Van Nostrand Reinhold, 1970.
[55] G. F. Knoll, Radiation Detection and Measurement, 4th ed., Wiley, 2010.
[56] G. Gilmore, Practical Gamma-Ray Spectrometry, 2nd ed., Wiley, 2008.
[57] K. Debertin and R. G. Helmer, Gamma- and X-Ray Spectrometry with Semiconductor Detectors, North-Holland, 1988.
[58] JCGM 100:2008, Evaluation of Measurement Data - Guide to the Expression of Uncertainty in Measurement, BIPM, 2008.
[59] P. R. Bevington and D. K. Robinson, Data Reduction and Error Analysis for the Physical Sciences, 3rd ed., McGraw-Hill, 2003.
[60] ISO 11929, Determination of the Characteristic Limits for Measurements of Ionizing Radiation, International Organization for Standardization.
[61] OECD Nuclear Energy Agency, International Handbook of Evaluated Criticality Safety Benchmark Experiments, ICSBEP, NEA Data Bank.
[62] OECD Nuclear Energy Agency, International Reactor Physics Experiment Evaluation Project Handbook, IRPhEP, NEA Data Bank.
[63] World Nuclear Association, Thorium, Information Library, updated 2024.
[64] UK National Nuclear Laboratory, The Thorium Fuel Cycle: An Independent Assessment, Position Paper, 2010.
Kayıtlar
0 Yorumlar