Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron Cross Section”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Measurement of U 238 ( n , n ' γ ) cross section data and their impact on reaction models

A better knowledge of (n,xn) reaction cross sections is important for both reaction modeling and energy applications. This article focuses on inelastic scattering of neutrons off 238 U for which improvements are needed to better constrain evaluations and solve inconsistencies in nuclear power reactor calculations. A new precise measurement of (n,xnγ) reaction cross sections on 238 U has been performed at the GELINA (Geel Electron LINear Accelerator) neutron facility operated by EC-JRC-Geel (Belgium) with the GRAPhEME (GeRmanium array for Actinides PrEcise MEasurements) setup. The prompt γ-ray spectroscopy method coupled to time-of-flight measurements is used to extract (n,xnγ) cross section values which can be further combined to infer the total neutron inelastic scattering cross section. Cross section data for 18 γ transitions (five never measured before) are presented and compared to the data in the literature. Here, emphasis is especially given to the uncertainty determination to produce partial cross section data as accurate as possible. Due to intrinsic limitations of the experimental method, the use of additional nuclear structure information coupled with theoretical modeling is required to determine the total (n, n') cross section over the whole neutron energy range. We have investigated modeling aspects of the 238 U(n,n'γ) cross sections related to the description of compound nucleus and preequilibirum mechanisms as well as the discrete part of nuclear structure. Through comparison between experimental and calculated (n,n'γ) cross sections, we pinpoint inaccuracies in the description of specific reaction mechanisms and challenge recently implemented models. This helps improving the whole modeling of the (n, n') reaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Development of Thermal Scattering Law and Cross Sections of Hydrogen in Paraffin Wax

Paraffin wax is frequently used as a neutron moderator and shielding material. The main component of paraffin wax is straight-chain alkanes (n-alkanes). The deposition of paraffin wax is primarily attributed to the crystallization of n-alkanes. It is important to gain a deeper understanding of the mechanisms underlying the behavior of paraffin wax, which would impact the thermal scattering Law (TSL) and cross sections and affect the analysis of neutronic and critical systems. In this work, a classical molecular dynamics (CMD) simulation model was used in LAMMPS to create the TSL and cross sections at room temperature and pressure. To generate the required velocity autocorrelation functions (VACF), previously published data were used to validate the approach and models by comparing them with key model parameters. The phonon density of state (DOS) was calculated using Fourier transformation of the normalized VACF. This DOS was used as the primary input to estimate the TSL (S($a$, $β$)) and cross sections of hydrogen in paraffin wax. The TSL and cross sections of hydrogen were estimated using the Full Law Analysis Scattering System Hub (FLASSH) code. The cross section of hydrogen in paraffin wax is consistent with other hydrocarbon materials such as polyethylene with deviations due to structure in the lowest energy region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Modeling and Design of Micro Flux Monitors for Neutron Fluence Measurements

Measuring neutron flux is critical to understanding the characteristics of nuclear reactors and can provide information on the energy spectrum, total neutron fluence, and operational history of a facility. Such measurements are typically performed by irradiating high-purity materials in an operating reactor for a known length of time. By measuring and identifying activation products in the irradiated materials and referencing the energy-dependent cross section of each isotope, the neutron energy spectrum and total fluence can be determined. Historical techniques use a variety of small, individually pressed foils or wires for dosimeter devices, which can result in excessively high levels of radioactivity and complications in handling after irradiation. The present work describes a new neutron dosimetry concept, named the micro flux monitor, that uses semiconductor device fabrication processes to deposit thin films of flux monitoring materials in precise geometries onto a SiC wafer. This approach allows multiple dosimetry materials to be located on a single device, simplifying irradiation and subsequent analysis. The geometry, and therefore mass, of each material can be custom-tailored for a specific application—considering the expected neutron flux, duration of irradiation, and device retrieval time of a facility. Additionally, hundreds of identical dosimetry monitors can be produced simultaneously, reducing fabrication costs. Finally, the extremely thin metal dimensions of these devices could improve the fidelity of neutron fluence measurements by reducing neutron self-shielding and gamma self-attenuation effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

An Analytic Benchmark for Neutron Boltzmann Transport with Downscattering—Part I: Flux and Eigenvalue Solutions

Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Thermal Neutron Scattering Law for Beryllium Hydride and Critical Mass Calculations

The thermal neutron scattering law (TSL) for crystalline beryllium hydride (BeH 2 ) is developed from first-principles ab initio lattice dynamics calculations and the impact of neutron thermalization in this material on critical mass is estimated. BeH 2 has a body-centered orthorhombic crystal structure with 12 molecules per unit cell and a theoretical density of 0.755 g/cm 3 . The vibrational (phonon) densities of states for H and Be bound in BeH 2 are determined using VASP density functional theory and PHONON lattice dynamics calculations. The TSLs for H bound in BeH 2 , H(BeH 2 ), and Be bound in BeH 2 , Be(BeH 2 ), are then evaluated in the incoherent approximation from the calculated H and Be partial phonon density of states using FLASSH. Finally, critical mass as a function of 235 U loading density for bare and reflected BeH 2 moderated spheres is predicted from MC21 Monte Carlo neutron transport calculations using ENDF/B-VIII.0 cross sections and the H(BeH 2 ) and Be(BeH 2 ) TSL evaluations. Comparisons are made to water (H 2 O), polyethylene (CH 2 ), and beryllium oxide (BeO) as moderators. These critical mass predictions are a refinement upon the prior work by Rao and Srinivasan that neglected thermal neutron scattering effects. The minimum critical mass of a BeH 2 moderated assembly is estimated to be 0.207 kg 235 U for a 0.20 m thick BeO reflected sphere and 0.178 kg 235 U for a 0.40 m thick BeO reflected sphere.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Potential neutron-induced γ-ray background on natural tellurium relevant to 130 Te 0 ν ββ decay searches at the CUORE and SNO+ detectors

Gamma-ray production cross-section data have been obtained for the inelastic neutron scattering reactions 126,128,130 Te(n,n'γ) at five neutron energies between 3.6 and 10 MeV. We report data for the γ-ray energy region relevant to 0ν ββ decay of 130 Te with a Q ββ value of 2527.515 keV, assuming natural-abundance tellurium, as used at CUORE and SNO+. The natural abundance of 130 Te, 128 Te, and 126 Te is 34%, 32%, and 19%, respectively. For CUORE the γ-ray cascades from the excited state in 130 Te at 2527.06 keV and in 126 Te at 2533.85 keV are of concern. For SNO+, accounting for its inferior energy resolution, an additional four levels are important in 130 Te, an additional nine levels in 128 Te, and an additional eight levels in 126 Te. Of these, we report neutron-induced γ-ray production cross sections for the strongest transitions: the 2581.15 keV level in 130 Te, the 2494.20, 2508.06, and 2630.14 keV levels in 128 Te, and the 2496.83 and 2585.46 keV level in 126 Te. Here, the largest cross-section values were found for cascade γ-ray transitions to the ground state, while direct transitions to the ground state are very weak and were not observed in the present work. Both the CUORE and SNO+ detectors, however, may not be able to distinguish between cascade transitions and direct transitions to the ground state, making the neutron-induced excitation of the 2527.06 keV state of 130 Te in particular a potential problem for 0νββ decay searches of 130Te, because it matches its Q ββ value.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Damage Pseudo Cross Section Generator

This software is a python script for generating group damage pseudo-cross-sections for specific elements, and tabulating them in forms useful to python or to office-open workbooks (xlsx). The damage pseudo-cross-section is needed for calculating how much neutron radiation damage a structural material receives in a reactor measured in displacements per atom (DPA). The process of generating these pseudo-cross-sections is laborious as it requires running the nuclear data processing code, NJOY, for hundreds of isotopes, and then post-processing the data. This also allows many opportunities for human error to be involved. This script eliminates these problems by fully automating the process.

Gale, MicahD.↗
DOE OSTI2022

Investigation of the impact of difference between FRENDY and NJOY2016 on neutronics calculations

A nuclear data library is used as a starting input for all subsequent neutronics calculations. NJOY has been used worldwide as a nuclear data processing code to create cross section libraries for a long time. For the verification of NJOY method and for providing an alternative nuclear data processing tool, JAEA has been developed the new nuclear library processing code FRENDY. In this paper, nuclear calculations were performed using the ACE files and the multigroup libraries created by both FRENDY and NJOY, and the impacts on the neutronics characteristics due to nuclear data processing were investigated using those libraries. MCNP was used to compare the ACE files by calculating many benchmark problems including ICSBEP and it was confirmed that the k-eff values generally agree with each other within the range of statistical errors. The multigroup cross sections are verified by the BWR design codes LANCR/AETNA through calculation of a commercial BWR5 equilibrium core loaded with 9*9 fuels. It was confirmed that fuel assembly and core characteristics are consistent with each other. From the above investigations, it was confirmed that FRENDY can provide comparable continuous/multi-group neutron cross sections with NJOY. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Design and commissioning of the Correlated Gamma-Neutron Array for sCattering (CoGNAC)

The Correlated Gamma-Neutron Array for sCattering (CoGNAC) is a dual $n$-$γ$ detector array fielded at the Weapons Neutron Research facility at the Los Alamos Neutron Science Center. Here, the detector array utilizes elpasolite detectors and liquid scintillators to perform measurements of neutron elastic and inelastic scattering cross sections and emitted particle angular distributions. This manuscript details the design and experimental setup of CoGNAC as well as some of the preliminary studies performed to select the elpasolite detectors deployed in the array. Examples of the analysis methods available to CoGNAC are given in the context of past CoGNAC measurements of neutron scattering reactions on 28 Si and 56 Fe.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Study of the deformed halo nucleus 31 Ne with Glauber model based on microscopic self-consistent structures

In this work, the exotic deformed nucleus 31 Ne is studied with an approach that combines self-consistent structure and reaction theories. The deformed relativistic Hartree-Bogoliubov theory in continuum (DRHBc) is utilized to demonstrate that deformation and pairing correlations give rise to a halo structure with a large-amplitude p-wave configuration in 31 Ne. Then the valence nucleon wave functions and angle-averaged density distributions of 30 Ne from this theory are used as input for the Glauber reaction model to study the observables of neutron-rich neon isotopes to search halo signatures. With NL1 effective interaction, our predictions of the reaction cross sections for these exotic neon isotopes on a carbon target can better reproduce the experimental data than those from the relativistic mean field model for a spherical shape with resonances and pairing correlation contributions, and are roughly 3.3% (~ 50 mb) larger than those with Gaussian function-fitted densities of the core nuclei. The calculated one-neutron removal cross section at 240 MeV/nucleon, and the inclusive longitudinal momentum distribution of the 30 Ne residues from the 31Ne breakup reaction are largely improved over previous theoretical predictions and agree well with data. These reaction evaluations indicate a dilute density distribution in coordinate space and are a canonical signature of a halo structure. Moreover, our predictions with the NL3 and PK1 effective interactions give slightly better descriptions of reaction observables for exotic neon isotopes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Energy-resolved fast-neutron radiography using an event-mode neutron imaging detector

Energy-resolved fast-neutron radiography is a powerful non-destructive technique that can be used to remotely measure the quantity and distribution of elements and isotopes in a sample. This is done by comparing the energy-dependent neutron transmission of a sample with the known cross-sections of individual isotopes. The reconstruction of the composition is possible due to the unique features (e.g. resonances) in the cross-sections of individual isotopes. At short-pulsed (≲ 1 ns) neutron sources, such information is accessible via time-of-flight neutron imaging in principle, but requires a detector with nanosecond temporal resolution. Conventional neutron detectors can meet this requirement only by heavily compromising spatial resolution or efficiency. Here, we present a unique approach on fast neutron resonance radiography using a scintillator-based event-mode imaging detector at a short-pulsed neutron source, including first results on spatially mapped resonance profiles using MeV neutrons. The event mode approach applied in the presented detector allows recording of individual neutron interactions with nanosecond precision in time and sub-mm resolution in space. As a result, the entire available neutron energy spectrum can be measured for each pulse. At the same time, the use of a thick scintillator screen and lenses to focus the produced light results in a highly flexible field of view and a high interaction probability in the sensitive volume of the detector.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant at or Near Delayed Critical [Abstract]

Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2020

IAEA Photonuclear Data Library 2019

Here, photo-induced reaction cross section data are of importance for a variety of current or emerging applications, such as radiation shielding design and radiation transport analyses, calculations of absorbed dose in the human body during radiotherapy, physics and technology of fission reactors (influence of photo-reactions on neutron balance) and fusion reactors (plasma diagnostics and shielding), activation analyses, safeguards and inspection technologies, nuclear waste transmutation, medical isotope production and astrophysical applications. To address these data needs the IAEA Photonuclear Data library was produced in 1999, containing evaluated photo-induced cross sections and neutron spectra for 164 nuclides which were deemed relevant for the applications. Since the release of the IAEA Photonuclear Data Library however, new experimental data as well as new methods to assess the reliability of experimental cross sections have become available. Theoretical models and input parameters used to evaluate photo-induced reactions have improved significantly over the years. In addition, new measurements of partial photoneutron cross sections using mono-energetic photon beams and advanced neutron detection systems have been performed allowing for the validation of the evaluations and assessments of the experimental data. Furthermore, technological advances have led to the construction of new and more powerful gamma-beam facilities, therefore new data needs are emerging. We report our coordinated efforts to address these data needs and present the results of the new up-to-date evaluations included in the new updated IAEA Photonuclear Data Library consisting of 219 nuclides. The new library includes 188 new evaluations produced by the CRP evaluators, and one evaluation taken from the JENDL/PD-2016 library, while 20 evaluations were retained from the previous 1999 IAEA Photonuclear Data Library. In most of the cases, the photon energy goes up to 200 MeV. A total of 55 nuclides are new in this library reflecting the progress in measurements but also the developing data needs. In this paper we discuss the new assessment method and make recommendations to the user community in cases where the experimental data are discrepant and the assessments disagree. In addition, in the absence of experimental data, we present model predictions for photo-induced reaction cross section on nuclides of potential interest to medical radioisotope production.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Uranium integral fission product yields for a spectrally-shaped 14.1 MeV neutron source at the National Ignition Facility

This paper describes the experimental results for an energy tuning assembly created to modify the National Ignition Facility deuterium–tritium fusion neutron source into a notional thermonuclear and prompt fission neutron spectrum, which has applications in integral measurements, nuclear data benchmarks, and radiation effects on microelectronics. The Monte Carlo neutron transport utilized MCNP5 to estimate the ETA-modified fluence using the ENDF-B/VIII.0 and IRDFF-II continuous energy nuclear data libraries, and SCALE Sampler was used to estimate the systematic nuclear data covariance using ENDF-B/VII.1 and IRDFF-II in a 252-group structure. The experiment fielded eight activation foils and a highly enriched uranium sample. This provided fifteen reaction channels that were used in a forward-fit comparison to the modeled results and to unfold the neutron spectrum using STAYSL. Gamma-ray spectrometry was performed on the activation and highly enriched uranium foils, and the reduced Χ 2 between the modeled and experimental values was 1.21. The results from the STAYSL unfold, reduced Χ 2 = 1.62, indicated that the modeled neutron spectrum was achieved and the systematic nuclear data uncertainty associated with the neutron transport and activation product cross sections was representative of the experiment. Integral cumulative fission product yield data were collected for 37 mass chains with a combination of gamma-ray spectrometry and radiochemical analysis. Fission product analysis was generally in agreement with two models using a semi-empirical fit and the General Observables of Fission code, with the exception of mass chains 88, 109, 111, 112, 113, 129, 139, 142, 144, 151, and 156.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Many-body Nuclear Dynamics

The principal goal of this work was to better understand the fusion of neutron-rich nuclei a topic relevant to the fields of both nuclear physics and nuclear astrophysics. The work provides insight into the structure and reactions of neutron-rich nuclei namely the extent of their neutron density distribution and its polarizability. By comparing the fusion excitation functions for a chain of isotopes with a common target nucleus changes in the attractive nuclear potential are assessed. This change in the attractive potential is related to changes in the neutron density distribution with increasing number of neutrons or changes in fusion dynamics with increasing neutron number. The impact of the pairing of valence neutrons and protons on the fusion cross-section is also examined. Measurement of an isotopic chain is a powerful tool to address this topic. The experimental program made use of several different accelerator facilities. The core of the experimental work involved experiments at the ReA3 accelerator and the Facility for Rare Isotope Beams (FRIB) situated at Michigan State University, at GANIL, the French national nuclear physics laboratory, and the University of Notre Dame. At ReA3 the degree to which α-clusters associated with fusion of 28,30,32 Si + 28 Si result from the collision dynamics or reflect an initial α-cluster structure was explored. Alpha clusters observed in fusion reactions exceed the predictions of the standard statistical model. This experiment provides a measure of how clusterization is impacted by the increasing neutron-richess of the system. In the experiment at GANIL we investigated fusion in 19 O + 12 C and 20 O + 12 C. This experiment utilized the recently developed active-target detector MuSIC@Indiana. Using this proven, highly efficient, active-target detector we measured the fusion excitation function for these neutron-rich systems. Grounded by these experimental measurements, theoretical models were used to examine the role of unpaired valence neutrons on the fusion cross-section at energies just above the fusion barrier.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Learning nuclear cross sections across the chart of nuclides with graph neural networks

We explore the use of deep learning techniques to learn how nuclear cross sections change as we add or remove protons and neutrons. As a proof of principle, we focus on the neutron-induced reactions in the fast energy regime. Our approach follows a two-stage learning framework. First, we apply representation learning to encode cross section data into a latent space using either variational autoencoders (VAEs) or implicit neural representations (INRs). Then, we train graph neural networks (GNNs) on the resulting embeddings to predict missing values across the nuclear chart by leveraging the topological structure of neighboring isotopes. We demonstrate accurate cross section predictions within a 9 × 9 block of missing nuclei. We also find that the optimal GNN training strategy depends on the type of latent representation used, with VAE embeddings performing best under end-to-end optimization in the original space, while INR embeddings achieve better results when the GNN is trained only in the latent space. Furthermore, using clustering algorithms, we map groups of latent vectors into regions of the nuclear chart and show that VAEs and INRs can discover some of the neutron magic numbers. These findings suggest that deep-learning models based on the representation encoding of cross sections combined with graph neural networks hold significant potential in augmenting nuclear theory models, e.g., by providing reliable estimates of covariances of cross sections, including cross-material covariances.

Machine learning↗