Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nuclear reactions and processes”

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 73 records · Page 4

Sensitivity of the 13 C(α,n) 16 O S factor to the uncertainty in the level parameters of the near-threshold state

The 13 C(α,n) 16 O reaction is the main source of neutrons for the s and i processes in asymptotic giant branch stars and carbon-enhanced metal-poor stars, respectively. The reaction rate over the relevant temperature range from 0.1 to 0.3 GK translates into a center-of-mass energy range of 150 to 540 keV. Current measurements extend down to 300 keV, still requiring an extrapolation of the cross section. At these low energies, the high-energy tail of a 1/2 + state near the reaction threshold makes a significant contribution to the cross section, but its amplitude is still highly uncertain. In this work the uncertainties associated with the low-energy cross-section extrapolation are investigated, in particular the sensitivity to the energy, width, and Coulomb renormalized asymptotic normalization coefficient of the near-threshold resonance. Furthermore, recently it has been suggested that the energy of the near-threshold level may have a large impact on the extrapolation, but this is not found to be the case compared with the other sources of uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Investigation of the Role of the Role of Nuclear Physics in Heavy Element Nucleosynthesis, through the Study of Key Reactions, and the Improvement of Theoretical Reaction Rates (Final Report)

Neutrino-driven winds in core-collapse supernovae have been identified as a possible site for the production of elements heavier than iron. Traditionally, these neutrino-driven winds have been proposed as the site of the main r-process. Recent simulations fail to reproduce the conditions required for the main r-process. while they remain a promising site for producing the lightest elements beyond iron, e.g., Sr, Y, and Zr through the νp process. The efficiency of the νp process depends on the hydrodynamical conditions, the electron fraction (which is related to the neutrino properties), and the nuclear reactions on many short-lived nuclei with limited (if at all) experimental information. The reaction rates on these nuclei are based on theoretical predictions using the Hauser-Feshbach model. Recent sensitivity studies have highlighted the importance of neutron-induced reactions on these nuclei along the νp process path. This work aimed to experimentally constrain reaction rates that are known to play a key role in the neutrino-p process nucleosynthesis. A secondary subsequently-added objective was to start the implementation of techniques that improve the description of nuclear properties in the Hauser-Feshbach model by extending the microscopic nuclear level density description offered via the shell model to high excitation energies without using experiment-based renormalizations. The main objective of this work was the experimental constraint of the 56 Ni(n,p) 56 Co reaction rate via a measurement of the inverse reaction 56 Co(p,n) 56 Ni at the National Superconducting Cyclotron Laboratory (NSCL) and later the Facility for Rare Isotope Beams (FRIB). This reaction is considered the key one for determining the yields possible by the neutrino-p process. A technique for this type of measurement in inverse kinematics at low energies did not exist before this work. The work also had two secondary objectives. First, to contribute to efforts to measure the same reaction in direct kinematics using a radioactive target at Los Alamos National Laboratory (LANL), and second, to advance work to implement shell-model-deduced microscopic level densities in Hauser-Feshbach calculations. The project has resulted in the development of the first technique to perform (p,n) cross-section measurements in relevant-for-astrophysics low energies in inverse kinematics using a magnetic spectrometer or separator, and neutron detectors for neutron tagging. It has also resulted in the precise measurement of the cross-section of the 40 Ar(p,n) 40 K reaction in a proof-of-principle experiment realized by using a beam-line quadrupole of the ReA3 accelerator of NSCL/FRIB. As part of this project the technique was successfully adapted to make use of the superior acceptance of the Separator for Capture Reactions (SECAR) at FRIB. In this project, the required experimental setup simulations and beam optics were developed and tested with the measurement of the 58 Fe(p,n) 58 Cu reaction cross-section. Additionally, this project contributed with simulation work to the development of a technique to measure (n,p) reactions with radioactive targets at LANL, and which resulted in the measurement of the key 56 Ni(n,p) 56 Co reaction cross-section at neutron energies above ≈1 MeV. Last, the project initiated work in the development of shell model based level densities using the moments method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Pathfinding quantum simulations of neutrinoless double- β decay

We present results from co-designed quantum simulations of the neutrinoless double- β decay of a simple nucleus in 1+1D quantum chromodynamics using IonQ’s Forte-generation trapped-ion quantum computers. Electrons, neutrinos, and up and down quarks are distributed across two lattice sites and mapped to 32 qubits, with an additional 4 qubits used for flag-based error mitigation. A four-fermion interaction is used to implement weak interactions, and lepton-number violation is induced by a neutrino Majorana mass. Quantum circuits that prepare the initial nucleus and time evolve with the Hamiltonian containing the strong and weak interactions are executed on IonQ Forte Enterprise. Enabled by tuned model parameters, lepton-number violation is observed in real time, providing a clear signal of neutrinoless double- β decay. This was made possible by co-designing the simulation to maximally utilize the all-to-all connectivity and native gate-set available on IonQ’s quantum computers. Quantum circuit compilation techniques and co-designed error-mitigation methods, informed from executing benchmarking circuits with up to 2,356 two-qubit gates, enabled observables to be extracted with high precision. We discuss the potential of future quantum simulations to provide yocto-second resolution of the reaction pathways in these, and other, nuclear processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Low Energy Analysis of NuclEar Reactions v0.0.1

This is a lightweight, easy to use Python analysis package for analyzing public low-energy nuclear reaction data relevant to understanding basic nuclear processes, often of relevance to Solar Fusion, to improve our theoretical understanding of them, and enable connections with input from lattice QCD. It performs the analysis in a Bayesian Framework and supports Bayesian Model Averaging to provide a robust uncertainty quantification.

Walker-Loud, André↗

Electromagnetic Dissociation Cross Sections for High LET Fragments

Nuclear interaction cross sections are used in space radiation transport codes to calculate the probability of fragment emission in high energy nucleus-nucleus collisions. Strong interactions usually dominate in these collisions, but electromagnetic (EM) interactions can also sometimes be important. Strong interactions typically occur when the projectile nucleus hits a target nucleus, with a small impact parameter. For impact parameters larger than the sum of the nuclear radii, EM reactions dominate and the process is called electromagnetic dissociation (EMD) if one of the nuclei undergo fragmentation. Previous models of EMD have been used to calculate single proton (p) production, single neutron (n) production or light ion production, where a light ion is defined as an isotope of hydrogen (H) or helium (He), such as a deuteron (2H), a triton (3H), a helion (3He) or an alpha particle (4He). A new model is described which can also account for multiple nucleon production, such as 2p, 2n, 1p1n, 2p1n, 2p2n, etc. in addition to light ion production. Such processes are important to include for the following reasons. Consider, for example, the EMD reaction 56Fe + Al --> 52Cr + X + Al, for a 56Fe projectile impacting Al, which produces the high linear energy transfer (LET) fragment 52Cr. In this reaction, the most probable particles representing X are either 2p2n or 4He. Therefore, production of the high LET fragment 52Cr, must include the multiple nucleon production of 2p2n in addition to the light ion production of 4He. Previous models, such as the NUCFRG3 model, could only account for the 4He production process in this reaction and could not account for 2p2n. The new EMD model presented in this work accounts for both the light ion and multiple nucleon processes, and is therefore able to correctly account for the production of high LET products such as 52Cr. The model will be described and calculations will be presented that show the importance of light ion and multiple nucleon production. The work will also show that EMD reactions contribute most to those fragments with the highest LET.

Norbury, John↗

r -process Nucleosynthesis from Hyperaccreting Neutron Stars in Common Envelopes

We investigate nuclear reactions and feedback in hyperaccreting neutron star environments, considering accretion rates in the range 0.3–3 × 10 4 M ⊙ yr −1 , typical of short-period compact-object binaries in common envelopes. Our models account for weak reactions, neutrino energy loss, nuclear energy release, pair production, degenerate equations of state, and general relativistic hydrodynamics. Depending on the accretion rates, these systems can develop both proton- and neutron-rich atmospheres with strong convective instabilities linking the neutrino sphere to the outgoing accretion shock inside the radiation trapping zone. Convection drives nucleons through multiple heating and cooling cycles, with photodisintegration dominating during the heating phase and heavy element synthesis during the cooling phase, ejecting material with abundances that depend on the accretion rate and depth of the final decompression trajectory. The turbulent nature of convective currents is conducive to creating a wide range of nuclear products through a variety of effects, including nuclear statistical equilibrium freeze-out and the r-, p-, and γ-processes. We also observe a novel multistep process in reheated trajectories, consisting of proton-capture and photodissociation reactions operating on r-process seeds, producing overall neutron-deficient isotopes. A significant amount of infalling gas experiences high entropy and short (millisecond) freeze-out timescales capable of making r-process elements with high overabundances through a disequilibrium effect between neutrons and α-particles that does not require an excess of neutrons.

79 ASTRONOMY AND ASTROPHYSICS↗

Isochronic Evolution and the Radioactive Decay of r-process Nuclei

Abstract We report on the creation and application of a novel decay network that uses the latest data from experiment and evaluation. We use the network to simulate the late-time phase of the rapid neutron capture ( r ) process. In this epoch, the bulk of nuclear reactions, such as radiative capture, have ceased, and nuclear decays are the dominant transmutation channels. We find that the decay from short-lived to long-lived species naturally leads to an isochronic evolution in which nuclei with similar half-lives are populated at the same time. We consider random perturbations along each isobaric chain to initial solar-like r -process compositions to demonstrate the isochronic nature of the late-time phase of the r -process. Our analysis shows that detailed knowledge of the final isotopic composition allows for the prediction of late-time evolution with a high degree of confidence despite uncertainties that exist in astrophysical conditions and the nuclear physics properties of the most neutron-rich nuclei. We provide the time-dependent nuclear composition in the Appendix as supplemental material.

79 ASTRONOMY AND ASTROPHYSICS↗

Neural network emulation of spontaneous fission

Large-scale computations of fission properties are an important ingredient for nuclear reaction network calculations simulating rapid neutron-capture process (the 𝑟 process) nucleosynthesis. Due to the large number of fissioning nuclei potentially contributing to the 𝑟 process, a microscopic description of fission based on nuclear density functional theory (DFT) is computationally challenging. Here, we explore the use of neural networks (NNs) to construct DFT emulators capable of predicting potential energy surfaces and collective inertia tensors across the whole nuclear chart, starting from a minimal set of DFT calculations. We use constrained Hartree-Fock-Bogoliubov (HFB) calculations to predict the potential energy and collective inertia tensor in the axial quadrupole and octupole collective coordinates, for a set of nuclei in the 𝑟-process region. We then employ NNs to emulate the HFB energy and collective inertia tensor across the considered region of the nuclear chart. Least-action pathways characterizing spontaneous fission half-lives and fragment yields are then obtained by means of the nudged elastic band method. The potential energy predicted by NNs agrees with the DFT value to within a root-mean-square error of 500 keV, and the collective inertia components agree to within an order of magnitude. These results are largely independent of the NN architecture. The exit points on the outer turning line are found to be well emulated. For the spontaneous fission half-lives the NN emulation provides values that are found to agree with the DFT predictions within a factor of 10 3 across more than 70 orders of magnitude. Neural networks are able to emulate the potential energy and collective inertia well enough to reasonably predict physical observables. Future directions of study, such as the inclusion of additional collective degrees of freedom and active learning, will improve the predictive power of microscopic theory and further enable large-scale fission studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Phenomenological R-Matrix parameterization of direct, doorway, and compound nuclear reactions [Slides]

The presented R-matrix formalism could improve evaluations of resolved resonance nuclear cross section data by parameterizing direct reaction processes via an orthonormal channel mixing matrix and doorway state processes. Together, they can parameterize direct-semidirect capture. A formal similarity among the R-, K-, T- matrix formalisms was leveraged to parameterize direct and doorway processes. The advantage of this R-matrix parameterization is that it is easily extensible due to its application of projection operators directly upon the (H-E) on p. 6.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Hydrogen density mapping in biomolecular crystals through dynamic nuclear polarization

Many fundamental biological processes, including those in photosynthetic reaction centers and enzyme active sites, involve charge and energy transfer, bond cleavage, protonation and hydrogen bonding. Because H atoms play such central roles in these reactions, accurately determining their positions is essential. Yet, conventional X-ray crystallography primarily resolves the heavy atoms in biological structures and provides limited insight into hydrogen, even at atomic resolution. Neutron macromolecular crystallography (NMC) overcomes this limitation by offering exceptional sensitivity to hydrogen and deuterium. Here, we present a theoretical framework for the development of dynamic nuclear polarization NMC (DNP-NMC) techniques, which exploit the alignment of neutron and proton nuclear spins to enhance and tune the hydrogen signal contribution. The DNP-NMC approach advances the resolution of H atoms within biomolecular crystals, whether bound to protein residues or present in solvent. The method establishes key relationships for the coherent structure factor of polarized neutron scattering from hydrogenous matter. It theoretically achieves full accuracy in phase reconstruction and offers a path to improve neutron structure determination, achieving accuracies exceeding ≳80% by incorporating titration states. Using a variant of the hybrid input/output phase-retrieval algorithm, it allows recovery of the hydrogen density with ≳90% phase accuracy. In conclusion, we further discuss sources of experimental uncertainty for the upcoming DNP-enabled, quasi-Laue IMAGINE-X experiment at Oak Ridge National Laboratory's High Flux Isotope Reactor.

dynamic nuclear polarization↗

YAHFC: A Code Framework to Model Nuclear Reactions and Estimate Correlated Uncertainties

Reaction modeling is a key ingredient in designing experiments and interpreting their results, and is an essential component in the process of evaluating nuclear data and assembling nuclear data libraries used in nuclear technology applications. Typically, experimental data are available only for a handful of reaction channels and theory models are used to fill in the gaps. In addition, theory is often called upon as the arbitrator between discrepant data. Most importantly, theory and modeling are required for an accurate determination of uncertainties in the evaluated data and the correlations between the multiple channels. A fast, accurate, and flexible modeling capability has been developed at LLNL with the code system YAHFC (Yet Another Hauser-Feshbach Code). YAHFC is a Monte Carlo, Hauser-Feshbach code framework, making full use of dynamic memory allocation, derived types, and parallel computing. YAHFC can generate events to simulate experiments and is guiding experiments designed to measure inelastic neutron scattering from actinide targets. YAHFC is also being used to analyze decays from surrogate experiments, thereby enabling the inference of reaction cross sections inaccessible by direct measurement. Finally, by modeling nuclear reactions with constraints from experimental data, YAHFC can deliver complete nuclear data libraries, with evaluated uncertainties, using the modernized Generalized Nuclear Data Structure (GNDS).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Fission fragment decay simulations with the CGMF code

The CGMF code implements the Hauser-Feshbach statistical nuclear reaction model to follow the de-excitation of fission fragments by successive emissions of prompt neutrons and γ rays. The Monte Carlo technique is used to facilitate the analysis of complex distributions and correlations among the prompt fission observables. Starting from initial configurations for the fission fragments in mass, charge, kinetic energy, excitation energy, spin, and parity, Y(A,Z,KE,U,J,π), CGMF samples neutron and γ-ray probability distributions at each stage of the decay process, conserving energy, spin and parity. Nuclear structure and reaction input data from the RIPL3 library are used to describe fission fragment properties and decay probabilities. Characteristics of prompt fission neutrons, prompt fission γ rays, and independent fission yields can be studied consistently. Correlations in energy, angle and multiplicity among the emitted neutrons and γ rays can be easily analyzed as a function of the emitting fragments. Finally, in its current version, CGMF is limited to certain spontaneous fission reactions on Pu and Cf isotopes, as well as neutron-induced fission reactions up to 20 MeV on U, Np, and Pu isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Dihadron azimuthal correlations in deep-inelastic scattering off nuclear targets

We measured the nuclear dependence of the di-pion azimuthal correlation function in deepinelastic scattering (DIS) using the CEBAF Large Acceptance Spectrometer (CLAS) and a 5 GeV electron beam. As the nuclear-target size increases, transitioning from deuterium to carbon, iron, and lead, the correlation function broadens monotonically. Its shape exhibits a significant dependence on kinematics, including the transverse momentum of the pions and the difference in their rapidity. None of the various Monte-Carlo event generators we evaluated could fully replicate the observed correlation functions and nuclear effects throughout the entire phase space. As the first study of its kind in DIS experiments, this research provides an important baseline for enhancing our understanding of the interplay between the nuclear medium and the hadronization process in these reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Photo- and Electro-Induced Hadron Production from Nuclei at Jefferson Laboratory

Understanding many-body knockout processes is crucial for nuclear physics, particularly in photo- and electro-induced reactions. In turn, understanding two- and three-body forces, including higher-order forces, is vital for a complete understanding of atoms. We present photo-induced many-proton knockout processes, with multiplicities from 1 to 6, using 12C, CH2, and C4H9OH targets in the g9a FROST dataset. Our analysis covers photon energies from 600 to 4500 MeV, significantly expanding current world data. Comparing our experimental data to the state-of-the-art GiBUU model offers a new challenge in the model’s theoretical description of many-body processes. GiBUU reasonably describes the data at lower photon energies but struggles at higher energies and missing masses, likely due to missing processes, such as initial 3-pion photoproduction. Our results will inform future developments in describing proton knockout processes, indicating GiBUU’s overall reasonable description of many-proton knockout data up to around 2.2 GeV. We also assess various electro-induced reactions using 2D, 12C, and 40Ar targets in the RGM dataset. Our results, obtained at electron beam energies of 2, 4, and 6 GeV, are compared in detail to GENIE and GiBUU, two widely used theory models in neutrino oscillation experiments. Discrepancies between model predictions and experimental data underscore the need for refining the two theoretical models. Despite discrepancies, GiBUU provides a more accurate modelling of electro-induced reactions, especially for 40Ar - crucial for future neutrino oscillation facilities such as DUNE. Understanding the fundamental nuclear physics involved in neutrino-nuclei interactions is essential for reducing the systematic uncertainties in extracting neutrino oscillation parameters. Many-body processes significantly contribute to the background processes observed in neutrino-nuclei interactions, hence the results from both analyses are crucial for developing the theoretical framework for the underlying nuclear physics.

Williams, Rhidian↗

Neutron diffraction: a primer

Because of the neutron’s special properties, neutron diffraction may be considered one of the most powerful techniques for structure determination of crystalline and related matter. Neutrons can be released from nuclear fission, from spallation processes, and also from low-energy nuclear reactions, and they can then be used in powder, time-of-flight, texture, single crystal, and other techniques, all of which are perfectly suited to clarify crystal and magnetic structures. With high neutron flux and sufficient brilliance, neutron diffraction also excels for diffuse scattering, for in situ and operando studies as well as for high-pressure experiments of today’s materials. For these, the wave-like neutron’s infinite advantage (isotope specific, magnetic) is crucial to answering important scientific questions, for example, on the structure and dynamics of light atoms in energy conversion and storage materials, magnetic matter, or protein structures. In this primer, we summarize the current state of neutron diffraction (and how it came to be), but also look at recent advances and new ideas, e.g., the design of new instruments, and what follows from that.

36 MATERIALS SCIENCE↗

Development of A Monte Carlo Radiation Transport Code System For HEDS: Status Update

Modifications of the Monte Carlo radiation transport code HETC are underway to extend the code to include transport of energetic heavy ions, such as are found in the galactic cosmic ray spectrum in space. The new HETC code will be available for use in radiation shielding applications associated with missions, such as the proposed manned mission to Mars. In this work the current status of code modification is described. Methods used to develop the required nuclear reaction models, including total, elastic and nuclear breakup processes, and their associated databases are also presented. Finally, plans for future work on the extended HETC code system and for its validation are described.

Townsend, Lawrence W.↗