Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Neutron Capture”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Porosity in nuclear graphite and its impact on nuclear reactor science and criticality safety applications

Porosity in nuclear-grade graphite significantly influences its low-energy neutron scattering, yet its effect on underlying phonon properties remains debated. Here, this work integrates inelastic and small-angle neutron scattering (INS/SANS) experiments, advanced atomistic simulations with a novel machine-learned potential (DeepMD), total cross-section measurements, and neutronics calculations (SCALE, MCNP, OpenMC) to investigate porosity’s impact on neutron thermalization. INS measurements on diverse graphite grades reveal no discernible porosity effect on phonon spectra, which align with crystalline graphite. Conversely, total cross-section data below ≈10 meV show increased scattering attributable to SANS. Our DeepMD simulations demonstrate that realistic micropores do not distort phonon spectra, challenging the assumptions in current ENDF/B-VIII.1 porosity thermal scattering laws (TSLs). These TSLs, based on random atom removal, produce unphysical phonon spectra and inflate inelastic cross-sections. Augmenting a crystalline TSL with an SANS component accurately captures experimental total cross-sections. Neutronics benchmarks (ICSBEP/IRPhE) show ENDF porosity TSLs unphysically increase neutron multiplication factor, keff. Crucially, incorporating SANS physics (NCrystal/OpenMC) indicates accurately modeled porosity negligibly affects keff, reactor physics, or criticality safety.

Critical benchmarks↗
DOE OSTI2021

Neutron Cross Section Measurement In The Protodune-SP Experiment

Understanding the detector response to neutrons will be critical for performing neutrino oscillation analyses in the next-generation Deep Underground Neutrino Experiment (DUNE). The DUNE physics program is centered around measuring the neutrino flavor composition as a function of their energy both at the near and the far detector. Neutrinos in the DUNE beam will have energies ranging between 100~MeV and 10~GeV, which is significant, because individual neutrino energies will not be known beforehand and will have to be reconstructed. Neutrino interactions in DUNE will produce leptons and hadrons -- including protons, pions, and neutrons. Neutrons can transport energy away from their origin and sometimes go undetected. In addition to the primary neutrons produced by the neutrino, subsequent interactions of the charged hadrons can result in secondary neutrons. Neutrons are a source of missing energy and will bias the neutrino energy measurement. Currently, there is also a 20\% energy scale uncertainty and a 40\% uncertainty on the energy resolution for neutrons in DUNE, which must be addressed. ProtoDUNE Single-Phase (ProtoDUNE-SP) is a 770-ton prototype for the DUNE far detector and was designed to both validate the technology that will be employed in DUNE and to measure cross sections for the charged hadrons (pions, protons, and kaons) at the relevant energies for DUNE. The ProtoDUNE-SP experiment, therefore, is in a unique position to characterize the secondary neutron component for DUNE. This is achieved by searching for candidate neutron interactions in ProtoDUNE-SP events and using these to facilitate a measurement of the neutron inelastic cross section as well as an estimate of the neutron energy and number. The cross section measurement presented here is based on neutrons produced in 1~GeV, $\pi^+$ events captured in 2018 by ProtoDUNE-SP in accordance with the production and cross section models in the GEANT4 simulation toolkit, version 4.1 0.6p1. T he best-fit neutron inelastic cross section, in the kinetic energy range of 114 to 314~MeV, is $1.24_{-0.08}^{+0.10}$~(stat.~$\oplus$~syst.)~barns.

Rivera, David Orlando↗
DOE OSTI2026

Pulsed-Neutron Experiments at the Inherently Safe Subcritical Assembly

The pulsed neutron technique is a powerful, dynamic method to assay the reactivity of a multiplying system. This work presents the novel application of the pulsed neutron technique to the Inherently Safe Subcritical Assembly, an experimental configuration accepted by the International Criticality Safety Benchmark Evaluation Project Handbook. The experiments were replicated with COG11.3, a continuous-energy Monte Carlo code. The pulsed neutron data were analyzed using the Sjöstrand and Gozani area-ratio methods and by extracting the prompt neutron decay constant. Subsequent static k-eigenvalue and 𝛼-eigenvalue simulations were also performed for the same configurations. Neutron detector dead-time effects from the experiments were corrected using the Backwards Extrapolation Method and shown to have a negligible impact on the estimated reactivities. The results highlight that capturing time-dependent effects like delayed neutron precursor buildup are essential to accurately reproduce experimental results. They also highlight the importance of shielding the detectors from generator source neutrons in deeply subcritical configurations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Test for the Origin of Solar Mass Black Holes

Solar-mass black holes with masses in the range of ~1–2.5 M ⊙ are not expected from conventional stellar evolution, but can be produced naturally via neutron star (NS) implosions induced by capture of small primordial black holes (PBHs) or from accumulation of some varieties of particle dark matter. We argue that a unique signature of such “transmuted” solar-mass BHs is that their mass distribution would follow that of the NSs. This would be distinct from the mass function of black holes in the solar-mass range predicted either by conventional stellar evolution or early Universe PBH production. We propose that analysis of the solar-mass BH population mass distribution in a narrow mass window of ~1–2.5 M ⊙ can provide a simple yet powerful test of the origin of these BHs. Recent LIGO/VIRGO gravitational wave (GW) observations of the binary merger events GW190425 and GW190814 are consistent with a BH mass in the range ~1.5–2.6 M ⊙ . Though these results have fueled speculation on dark matter-transmuted solar-mass BHs, we demonstrate that it is unlikely that the origin of these particular events stems from NS implosions. Data from upcoming GW observations will be able to distinguish between solar-mass BHs and NSs with high confidence. This capability will facilitate and enhance the efficacy of our proposed test.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Celestial-body focused dark matter annihilation throughout the Galaxy

Indirect detection experiments typically measure the flux of annihilating dark matter (DM) particles propagating freely through galactic halos. We consider a new scenario where celestial bodies “focus” DM annihilation events, increasing the efficiency of halo annihilation. In this setup, DM is first captured by celestial bodies, such as neutron stars or brown dwarfs, and then annihilates within them. If DM annihilates to sufficiently long-lived particles, they can escape and subsequently decay into detectable radiation. This produces a distinctive annihilation morphology, which scales as the product of the DM and celestial body densities, rather than as DM density squared. We show that this signal can dominate over the halo annihilation rate in γ -ray observations in both the Milky Way Galactic center and globular clusters. We use Fermi and H.E.S.S. data to constrain the DM-nucleon scattering cross section, setting powerful new limits down to ~10 -39 cm 2 for sub-GeV DM using brown dwarfs, which is up to 9 orders of magnitude stronger than existing limits. We demonstrate that neutron stars can set limits for TeV-scale DM down to about 10 -47 cm 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Cascade Models in Simulation of Extended Heavy Targets Irradiated by Accelerated Proton and Deuteron Beams

The paper presents a survey of the main numerical models used for simulation of interaction of accelerated particle beams with target nuclei. These models form the core of the software for simulation of various experiments and experimental facilities both for scientific and applied purposes. The beam and target parameters considered in detail in this study (protons and deuterons with energies from 0.66 to 4 AGeV and bulk U targets) cover the range of interest in development of new concepts of nuclear power production aided by accelerated particle beams.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

The Evolution of Inclined Binary Black Holes in the Disks of Active Galactic Nuclei

The accretion disks that fuel active galactic nuclei (AGNs) may house numerous stars and compact objects, formed in situ or captured from nearby star clusters. Embedded neutron stars and black holes may form binaries and eventually merge, emitting gravitational waves detectable by LIGO/VIRGO. AGN disks are a particularly promising environment for the production of high-mass gravitational-wave events involving black holes in the pair-instability mass gap, and may facilitate electromagnetic counterparts to black hole binary mergers. However, many orders of magnitude separate the typical length scales of binary formation and those on which gravitational waves can drive binary inspirals, making binary mergers inside the disk uncertain. Previous hydrodynamical simulations of binaries have either been restricted to two dimensions entirely, or focused on binaries aligned with the midplane of the disk. Herein we present the first three-dimensional, high-resolution, local-shearing-box, inviscid hydrodynamical simulations of disk-embedded binaries over a range of orbital inclinations. We find that retrograde binaries can shrink up to 4 times as quickly as prograde binaries, and that all binaries not perfectly aligned (or anti-aligned) with the AGN disk are driven into alignment. An important consequence of this is that initially retrograde binaries will traverse the inclinations where von Zeipel–Lidov–Kozai oscillations can drive binary eccentricities to large values, potentially facilitating mergers. We also find that interactions with the AGN disk may excite eccentricities in retrograde binaries and cause the orbits of embedded binaries to precess.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

NCERC-CL TREAT Free Field Characterization Foil Reaction Rate Results

This report summarizes the analysis and preliminary MCNP reaction rate modeling of free field characterization measurements performed in the Transient Reaction Test Facility (TREAT) at Idaho National Laboratory (INL) between 12/18/2023 and 2/22/2024. A series of seven irradiations were performed using the Big-BUSTER (Broad Use Specimen Transient Experiment Rig) core configuration. Foil sets including 19.5% LEU-Zr alloy wire, S, Au, Fe, Ni, Co, Ti, and Zr were deployed in each irradiation. Identical foil sets were supplied to Los Alamos National Laboratory (LANL), Lawrence Livermore National Laboratory (LLNL), and INL. All foils were supplied by INL. Additional details on the experiment can be found in [2]. The foil sets were shipped to the National Criticality Experiments Research Center Counting Laboratory (NCERC-CL) branch at LANL (NCERC-CL NISC) and were received on 3/25/2024. Additional details on receipt and counting are reported in [3]. All reactor dosimetry measurements were performed in adherence to the ASTM standards applicable to reactor dosimetry. Measured reaction rates from Au, Co, Fe, Ni, and Ti are reported. Select reaction rate ratios and simulated reaction rates and ratios are discussed. Reaction rates for all measured reaction products were modeled with MCNP 6.3.1 and a TREAT input deck supplied by Edward Lum. Neutron emission estimates calculated using the measured and simulated reaction rates show excellent agreement between the three fast-threshold n-p reactions with the percent differences being < 9% between the three reactions. The capture reactions had poorer, but still reasonable, agreement with < 20% percent differences between the neutron emission estimates of the three capture reactions. The percent difference between the fast-threshold and capture reactions is ∼ 100% for each reaction. This supports the TREAT model is not accurately modeling the neutron spectrum and additional measurements are required to characterize TREAT and match measurements to simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

R-Matrix Analysis of the Neutron-Induced Cross Sections on 143 Nd Measured at LANSC [Slides]

An established use of R-Matrix codes is being implemented at LANSCE. There are unique features accessible with detectors at LANSCE that can be used to perform more complete R-Matrix analysis. Further, Spin separation is a powerful technique that can be applied in DANCE measurements. Additionally, a unique method to measure transmission with DICER is being undertaken. Transmission and capture data were analyzed in this work. The study can be extended to other channels. The results of the 143 Nd data are presented in this work. New resonances found in this work: 21 compared to ENDF/B-VIII.0 and JEFF-3.3, and 24 compared to JENDL-5. Wrong spin assignment in 23 resonances in ENDF/B-VIII.0, 25 in JEFF-3.3 and 26 in JENDL-5. The 143 Nd and 147,149 Sm have been analyzed. Exciting future for new measurements and analysis!

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Connecting relativistic density functional theory to microscopic calculations

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Neutrino and gamma-ray signatures of inelastic dark matter annihilating outside neutron stars

We present a new inelastic dark matter search: neutron stars in dark matter-rich environments capture inelastic dark matter which, for interstate mass splittings between about 45 - 285 MeV, will annihilate away before becoming fully trapped inside the object. This means a sizable fraction of the dark matter particles can annihilate while being outsidethe neutron star, producing neutron star-focused gamma-rays and neutrinos. We analyze this effect for the first time and target the neutron star population in the Galactic Center, where the large dark matter and neutron star content makes this signal most significant. Depending on the assumed neutron star and dark matter distributions, we set constraints on the dark matter-nucleon inelastic cross-section using existing H.E.S.S. observations. We also forecast the sensitivity of upcoming gamma-ray and neutrino telescopes to this signal, which can reach inelastic cross-sections as low as ∼ 2 × 10 -47 cm 2 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Uncertainty Analysis for a Criticality Benchmark [Slides]

This presentation discusses how KRUSTY was analyzed with LANL’s MCNP neutron transport code, Monte Carlo N-Particle. In the Monte Carlo method, neutrons are born in fission events and fly around the geometry having events (e.g., scatter, capture, fission, escape) according to the natural probabilities. “Natural probabilities” are given as neutron cross sections , which are measured and evaluated nuclear data. One of the goals of the benchmark project is to improve the nuclear cross sections (and the codes that use them). In this work, the response of interest is called k eff , which describes the neutron multiplication in a near critical system. Like most quantities, k eff is not measured directly; it is only inferred. In summation, MCNP results for k eff matched the measurements extremely well. The evaluated uncertainties are asymmetric and smaller than 0.00100 (e.g., + 0.00080 / – 0.00062). The largest uncertainties are due to alignment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Uncertainty Analysis for a Criticality Benchmark [Slides]

This presentation discusses how KRUSTY was analyzed with LANL’s MCNP neutron transport code, Monte Carlo N-Particle. In the Monte Carlo method, neutrons are born in fission events and fly around the geometry having events (e.g., scatter, capture, fission, escape) according to the natural probabilities. “Natural probabilities” are given as neutron cross sections , which are measured and evaluated nuclear data. One of the goals of the benchmark project is to improve the nuclear cross sections (and the codes that use them). In this work, the response of interest is called k eff , which describes the neutron multiplication in a near critical system. Like most quantities, k eff is not measured directly; it is only inferred. In summation, MCNP results for k eff matched the measurements extremely well. The evaluated uncertainties are asymmetric and smaller than 0.00100 (e.g., + 0.00080 / – 0.00062). The largest uncertainties are due to alignment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Selected Uses of TSUNAMI in Critical Experiment Design and Analysis

Validation in criticality safety is performed by comparing the results of critical experiments with the calculated results from models of the experiments using the computational method to be validated. Laboratory critical experiments are controlled systems that achieve a k eff of approximately 1 and enable investigation of the parameters at which such a critical condition is achieved. For the critical experiments used in a validation to capture the biases of the materials and neutron energy spectra of interest, those materials must be included in the experiment such that they influence k eff or another observable parameter with statistical significance. This paper discusses the use of sensitivity uncertainty (S/U) methods to develop critical experiments for various purposes. S/U techniques are useful for understanding the underlying components of nuclear data which affect the k eff or another parameter of a given configuration. S/U calculations are most commonly used to compare existing experiments to applications of interest; however, S/U techniques can also be used to identify, optimize, or assess features of proposed experiments so that they can better test specific portions of nuclear data or match an application of interest. The S/U techniques discussed here are from the TSUNAMI code system. The two primary codes discussed in this work are TSUNAMI-3D, which implements the KENO criticality code to calculate the sensitivity of k eff to nuclear data, and TSAR, which calculates the sensitivity of a reactivity difference between two configurations based on their TSUNAMI-3D generated sensitivity profiles. The methods used in these tools are discussed in more detail in the SCALE manual. This paper is one of a series on the development and use of TSUNAMI tools. The other papers address development of TSUNAMI methods and a review of TSUNAMI applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Criticality Experiments to Reduce Compensating Errors in Plutonium Nuclear Data

Compensating errors between nuclear data observables in a library can adversely impact application simulations. The primary goal of the EUCLID project (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) is to reduce compensating errors in nuclear data. A new criticality experiment, described in this work, was designed with the specific target nuclear data of 239 Pu fission, inelastic scattering, elastic scattering, capture, nu-bar, and prompt fission neutron spectrum (PFNS). This work will focus on the design and execution of the EUCLID experiment, performed on the Planet vertical lift critical assembly machine at the National Criticality Experiments Research Center (NCERC). The criticality experiment includes two different configurations with very different geometries: one is cube-like to minimize neutron leakage while the other is slab-like to maximize leakage. Having these two widely varying configurations allows the scattering sensitivities of 239 Pu to the neutron multiplication factor to be greatly changed while minimally impacting the other cross section sensitivities. Both configurations utilize the Pu ZPPR (Zero Power Physics Reactor) plates as fuel. The experiments were designed using a D-Optimality criteria, which is an optimization method minimizing the log-determinant of the adjusted nuclear data covariance for the target reactions. These experiments include not only inference of k eff , as done in all critical benchmark experiments, but several other responses as well, such as neutron multiplication measurements and reaction rate ratios. After analysis of the measured data is complete, adjustment of nuclear data will be performed to assess whether the new experimental data successfully reduced compensating errors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Direct NeTS sampling of nuclear graphite $S(α, β, T)$ in Serpent

For advanced reactor applications, Neural Thermal Scattering (NeTS) modules were developed to predict the thermal scattering law (TSL or $S(α, β, T)$) of a nuclear graphite neutron moderator. NeTS are multi-layer, feedforward artificial neural networks, which act as universal function approximators designed for TSL datasets. In this case, a 4-layer neural network with 164 neurons per layer is trained using FLASSH evaluated data in PyTorch and serialized as a torchscript dictionary to predict $S(α, β, T)$ on-the-fly. Relative, absolute and maximum percent deviations of NeTS from File 7 data generated using the FLASSH code are on the order of 0.01%, 0.1% and 1%, respectively, with low inference latencies of 0.000172 s per $S(α, β, T)$ at a given temperature. Capturing the full dimensionality of possible inelastic neutron-lattice interactions, NeTS functionality is embedded in the Serpent Monte Carlo code, where $S(α, β, T)_{NeTS}$ sampling is conducted on-the-fly and compared to ACE look-up-tables for predicting TREAT criticality. k-eff differences between sampling algorithms of 6 pcm are observed and are within the order of Monte Carlo uncertainty. Compared to discrete and continuous-energy ACE files (30 MB and 131 MB per temperature), the NeTS format is on the order of 200–300 kB for a continuous-temperature, interpolation-free representation of $S(α, β, T)$ and cross sections. NeTS-in-Serpent runtimes comparable with ACE look-up tables are achieved by scaling NeTS for high performance computing architectures with hybrid OpenMP + MPI parallelization. This work validates a novel, self-contained reactor physics framework for predictive cross sections, and demonstrates a general methodology for embedding modern machine learning libraries within existing neutronic analysis frameworks.

Nuclear Criticality Safety Program (NCSP)↗