Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Models & methods for nuclear reactions”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Fusion of 12 C+ 24 Mg far below the barrier: Evidence for the hindrance effect

Background: The phenomenon of fusion hindrance may have important consequences on the nuclear processes occurring in astrophysical scenarios, if it is a general behavior of heavy-ion fusion at extreme subbarrier energies, including reactions involving lighter systems, e.g., reactions in the carbon and oxygen burning stages of heavy stars. The hindrance is generally identified by the observation of a maximum of the S-factor vs energy. Whether there is an S-factor maximum at very low energies for systems with a positive fusion Q value is an experimentally challenging question. Purpose: Our aim has been to search for evidence of fusion hindrance in 12 C+ 24 Mg which is a medium-light m system with positive Q value for fusion, besides the heavier cases where hindrance is recognized to be a general phenomenon. 12 C+ 24 Mg is very close to the 16 O + 16 O and 12 C + 12 C systems that are important for the late evolution of heavy stars. Methods: The experiment has been performed in inverse kinematics using the 24 Mg beam from the XTU Tandem accelerator of LNL in the energy range 26-52 MeV with an intensity of 4-8 pnA. The targets were 12 C evaporations 50 mu g/cm 2 thick, isotopically enriched to 99.9%. The fusion-evaporation residues were detected at small angles by a E-Delta E-ToF detector telescope following an electrostatic beam deflector. Results: Previous measurements of fusion cross section for 12 C+ 24 Mg were limited to above-barrier energies. In the present experiment the excitation function has been extended down to similar or equal to 15 mu b and it appears that the S factor develops a clear maximum vs energy, indicating the presence of hindrance. This is the first convincing evidence of an S factor maximum in a medium-light system with a positive fusion Q value. These results have been fitted following a recently suggested method and a detailed analysis within the coupled-channels model that has been performed using a Woods-Saxon potential and including the ground state rotational band of 24 Mg. The coupled-channels calculations give a good account of the data near and above the barrier but overpredict the cross sections at very low energies. Conclusions: The hindrance phenomenon is clearly observed in 12 C+ 24 Mg, and its energy threshold is in reasonable agreement with the systematics observed for several medium-light systems. The fusion cross sections at the hindrance threshold show that the highest value (sigma(s) = 1.6 mb) is indeed found for this system. Therefore it may even be possible to extend the measurements further down in energy to better establish the position of the S-factor maximum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

225 Ac-MACROPATATE: A Novel α-Particle Peptide Receptor Radionuclide Therapy for Neuroendocrine Tumors

Neuroendocrine tumors (NETs) express somatostatin receptors (SSTRs) 2 and 5. Modified variants of somatostatin, the cognate ligand for SSTR2 and SSTR5, are used in treatment for metastatic and locoregional disease. Peptide receptor radionuclide therapy with 177 Lu-DOTATATE (DOTA-octreotate), a β-particle–emitting somatostatin derivative, has demonstrated survival benefit in patients with SSTR-positive NETs. Despite excellent results, a subset of patients has tumors that are resistant to treatment, and alternative agents are needed. Targeted α-particle therapy has been shown to kill tumors that are resistant to targeted β-particle therapy, suggesting that targeted α-particle therapy may offer a promising treatment option for patients with 177 Lu-DOTATATE–resistant disease. Although DOTATATE can chelate the clinically relevant α-particle–emitting radionuclide 225 Ac, the labeling reaction requires high temperatures, and the resulting radioconjugate has suboptimal stability. Methods: We designed and synthesized MACROPATATE (MACROPA-octreotate), a novel radioconjugate capable of chelating 225 Ac at room temperature, and assessed its in vitro and in vivo performance. Results: MACROPATATE demonstrated comparable affinity to DOTATATE (dissociation constant, 21 nM) in U2-OS-SSTR2, a SSTR2-positive transfected cell line. 225 Ac-MACROPATATE demonstrated superior serum stability at 37°C over time compared with 225 Ac-DOTATATE. Biodistribution studies demonstrated higher tumor uptake of 225 Ac-MACROPATATE than of 225 Ac-DOTATATE in mice engrafted with subcutaneous H69 NETs. Therapy studies showed that 225 Ac-MACROPATATE exhibits significant antitumor and survival benefit compared with saline control in mice engrafted with SSTR-positive tumors. However, the increased accumulation of 225 Ac-MACROPATATE in liver and kidneys and subsequent toxicity to these organs decreased its therapeutic index compared with 225 Ac-DOTATATE. Conclusion: 225 Ac-MACROPATATE and 225 Ac-DOTATATE exhibit favorable therapeutic efficacy in animal models. Because of elevated liver and kidney accumulation and lower administered activity for dose-limiting toxicity of 225 Ac-MACROPATATE, 225 Ac-DOTATATE was deemed the superior agent for targeted α-particle peptide receptor radionuclide therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2021

Neutron-unbound states in 31 Ne

Background: The Island of Inversion near the N = 20 shell gap is home to nuclei with reordered single-particle energy levels compared to the spherical shell model. Studies of 31 Ne have revealed that its ground state has a halo component, characterized by a valence neutron orbiting a deformed 30 Ne core. This lightly-bound nucleus with a separation energy of only S n = 170 keV is expected to have excited states that are neutron unbound. Purpose: The purpose of this experiment was to investigate the low-lying excited states in 31 Ne that decay by the emission of a single neutron. Methods: An 89 MeV/u 33 Mg beam impinged on a segmented Be reaction target. Neutron unbound states in 31 Ne were populated via a two-proton knockout reaction. The 30 Ne fragment and associated neutron from the decay of 31 Ne were detected by the MoNA-LISA-Sweeper experimental setup at the National Superconducting Cyclotron Laboratory. Invariant mass spectroscopy was used to reconstruct the two-body decay energy ( 30 Ne + n). Results: The two-body decay energy spectrum exhibits two features: a low-lying peak at 0.30 ± 0.17 MeV and a broad enhancement at 1.50 ± 0.33 MeV, each fit with an energy-dependent asymmetric Breit-Wigner line shape representing a resonance in the continuum. Accompanying shell model calculations using the FSU interaction within NuShellX, combined with cross-section cal culations using the eikonal reaction theory, indicate that these peaks in the decay energy spectrum are caused by multiple resonant states in 31 Ne. Conclusions: In this study, excited states in 31 Ne were observed for the first time. Transitions from calculated shell model final states in 31 Ne to bound states in 30 Ne are in good agreement with the measured decay energy spectrum. Cross-section calculations for the two-proton knockout populating 31 Ne states as well as spectroscopic factors pertaining to the decay of 31 Ne into 30 Ne are used to examine the results within the context of the shell model expectations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Model-independent determination of the dipole response of 66 Zn using quasimonoenergetic and linearly polarized photon beams

Background: Photon strength functions are an important ingredient in calculations relevant for the nucleosynthesis of heavy elements. The relation to the photoabsorption cross section allows to experimentally constrain photon strength functions by investigating the photoresponse of atomic nuclei. Purpose: Here, we determine the photoresponse of 66 Zn in the energy region of 5.6 MeV to 9.9 MeV and analyze the contribution of the 'elastic' decay channel back to the ground state. In addition, for the elastic channel electric and magnetic dipole transitions were separated. Methods: Nuclear resonance fluorescence experiments were performed using a linearly polarized quasi-monoenergetic photon beam at the High Intensity γ-ray Source. Photon beam energies from 5.6 to 9.9 MeV with an energy spread of about 3% were selected in steps of 200–300 keV. Two high purity germanium detectors were used for the subsequent γ-ray spectroscopy. Results: Full photoabsorption cross sections are extracted from the data making use of the monoenergetic character of the photon beam. For the ground-state decay channel, the average contribution of electric and magnetic dipole strengths is disentangled. The average branching ratio back to the ground state is determined as well. Conclusions: The new results indicate lower cross sections when compared to the values extracted from a former experiment using bremsstrahlung on 66 Zn. In the latter, the average branching ratio to the ground state is estimated from statistical-model calculations in order to analyze the data. Corresponding estimates from statistical-model calculations underestimate this branching ratio compared to the values extracted from the present analysis, which would partly explain the high cross sections determined from the bremsstrahlung data.

59 ≤ A ≤ 89↗
DOE OSTI2023

Nuclear rotation at the fission limit in Rf 254

A ground-state rotational band in the fissile nucleus Rf254 was observed for the first time. Levels up to spin 14 and excitation energy of 1.56 MeV were observed. The Rf254 nuclei were produced using the Pb206(Ti50,2n) fusion-evaporation reaction. It is the weakest reaction channel ever studied using in-beam γ-ray spectroscopic methods. The reaction products were separated from the beam in the Argonne gas-filled analyzer (AGFA). The Rf254 nuclei were implanted into a double-sided Si strip detector at the AGFA focal plane and tagged with subsequent ground-state spontaneous fission decays using temporal and spatial correlations. Prompt γ rays in coincidence with the Rf254 recoils were detected in the Gammasphere array of Ge detectors. In order to identify the ground-state rotational band in Rf254, a method for identifying rotational bands in low statistics γ-ray spectra was developed. The deduced Rf254 kinematic moment of inertia is smaller compared to neighboring even-even nuclei. This is most likely associated with a slightly lower quadrupole deformation and stronger pairing correlations in Rf254. The behavior of the moment of inertia as a function of rotational frequency is similar to that of the lighter N=150 isotones Fm250 and No252.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

CALPHAD modeling of uranium nitride (UN) fabrication routes enabled by first-principles calculations

The thermochemical details of fabricating uranium nitride (UN) by ammonolysis of uranium tetraflouride (UF 4 ) were determined using density functional theory (DFT) and CALculation of PHAse Diagrams (CALPHAD) computational methods. The thermochemical data of all binary, ternary, and quaternary U-H-N-F phases were computed using DFT, and the data for the phases that have not been measured experimentally, including UN 2 and NH 4 F(g), were combined with existing experimentally-determined data for CALPHAD modeling. The DFT data were benchmarked using experimental Gibbs energy of reaction and experimental thermochemical data for individual species. Phase diagrams relevant to the ammonolysis reaction are depicted, showing regions of stability for solid U-N, U-F and U-N-F phases. An unidentified phase produced in a previous experiment was identified as UN 0.95 F 1.2 (UNF) by comparing its X-ray diffraction spectrum to the experimental spectrum, and its formation during the fabrication of UN from UF 4 is supported by the simulated phase diagram. Here, it is calculated that UN 2 can be produced by the ammonolysis of UF 4 , but requires elevated temperatures, high NH 3 (g) partial pressure, and large amounts of flowing NH 3 (g) to avoid solid flu oride impurities in the uranium nitride. Likewise, U 2 N 3 can be produced instead at temperatures greater than 980 K. The use of silane (SiH 4 ) gas was investigated as a potential additive in the ammonolysis fabrication route to speed removal of fluorine. The addition of SiH 4 (g) offers little advantage to the removal of fluorine, and adds the complication of Si 3 N 4 formation. The use of DFT to fill in missing data to perform CALPHAD calculations demonstrated here allows for the determination of more comprehensive and trustworthy phase diagrams than the use of existing experimental data alone.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

1 H NMR Quantification of Aromatic Monomers from Reductive Catalytic Fractionation

Reductive catalytic fractionation (RCF) can produce high yields of aromatic monomers from lignin in native biomass. Quantification of these aromatic monomers is a well-known but demanding task, in part due to the lack of commercially available standards. Here, we demonstrate 1 H NMR spectroscopy as a complementary method to rapidly quantify aromatic monomer concentrations in RCF oils. The method exhibited good agreement with measurements from ultrahigh pressure liquid chromatography (UHPLC) for 96 RCF oils with varying monomer selectivity, with average absolute deviations of individual monomer yields between 0.5 and 1.1 wt % (relative 11–17%) and R 2 values above 0.9 compared to conventional UHPLC quantification. Quantification of S-type monomers, including for 4-ethylsyringol and 4-propenylsyringol, was generally reliable. The validity of G-type monomer quantifications depended on reaction selectivity due to overlap between peaks of 4-ethylguaiacol and 4-(3-hydroxypropyl)-guaiacol. The method could be applied on crude RCF oils without needing to perform the liquid–liquid extraction typically done for RCF reactions, thereby providing a convenient way to quantify lignin extraction and aromatic monomer yield. Overall, 1 H NMR spectroscopy can serve as a rapid primary quantification or secondary validation method for RCF monomer yield and selectivity measurements.

1H NMR spectroscopy↗
DOE OSTI2024

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2023

Mass measurement of 27 P to constrain type-I x-ray burst models and validate the isobaric multiplet mass equation for the $A = 27, T = \frac{3}{2}$ isospin quartet

Background: Light curves are the primary observable of type-I x-ray bursts. Computational x-ray burst models must match simulations to observed light curves. Most of the error in simulated curves comes from uncertainties in rp process reaction rates, which can be reduced via precision mass measurements of neutron-deficient isotopes in the rp process path. Purpose: Perform a precise atomic mass measurement of 27 P. Use this new measurement to calculate rp process reaction rates and input these rates into an x-ray burst model to reduce simulated light curve uncertainty. Use the mass measurement of 27 P to validate the isobaric multiplet mass equation (IMME) for the A = 27 T = 3/2 isospin quartet which 27 P belongs to. Method: High-precision Penning trap mass spectrometry utilizing the time-of-flight ion cyclotron resonance technique was used to determine the atomic mass of 27 P. The MESA code (Modules for Experiments in Stellar Astrophysics) was then used to simulate x-ray bursts using a one-dimensional multizone model to produce updated light curves. Results: The mass excess of 27 P was measured to be –670.7(6) keV, a 14-fold precision increase over the mass reported in the 2020 Atomic Mass Evaluation (AME2020). The 26 Si(p, $γ$) 27 P– 27 P($γ$, p) 26 Si rate equilibrium has been determined to a higher precision based on the precision mass measurement of 27 P. x-ray burst light curves were produced with the MESA code using the new reaction rates. Changes in the mass of 27 P seem to have minimal effect on light curves, even in burster systems tailored to maximize impact. Conclusion: The mass of 27 P does not play a significant role in x-ray burst light curves. It is important to understand that more advanced models do not just provide more precise results, but often qualitatively different ones. This result brings us a step closer to being able to extract stellar parameters from individual x-ray burst observations. In addition, the IMME has been validated for the A = 27, T = 3/2 quartet. The normal quadratic form of the IMME using the latest data yields a reduced $χ$ 2 of 2.9. Finally, the cubic term required to generate an exact fit to the latest data matches theoretical attempts to predict this term.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Discovery of the new isotope 251 Lr and the α-decay fine structure in 253 L

The products of the 203, 205 Tl ( 50 Ti, 2n) fusion-evaporation reactions were studied using the recently commissioned Argonne Gas-Filled Analyzer at Argonne National Laboratory. Two α-decay activities with energies of 9210(19) and 9246(19) keV and half-lives of $42^{+42}_{-14}$ and $24.4^{+7.0}_{-4.5}$ ms were observed which were followed by the known α decays of 247 Md and 243 Es . They are interpreted as originating from the 1/2 ⁻ [521] and 7/2 ⁻ [514] single-proton Nilsson states in the hitherto unknown isotope 251 Lr . From the measured Q α values the 1/2 ⁻ level was placed 117(27) keV above the 7/2 - level in 251 Lr in contrast to 255 Lr where the 1/2 ⁻ level is the lowest. Also, the α decay of 253 Lr was studied in more detail and a new α line at 8660(20) keV was found and a new half-life value of 2.46(32) s for an isomeric state in 253 Lr was measured. The 251, 253, 255 Lr Q α values were compared with predictions of various mass models. The relative energies of the 1/2 ⁻ [521] and 7/2 ⁻ [514] single-proton Nilsson states in 251, 253, 255 Lr isotopes were compared with results of the cranking shell model with pairing treated using the particle-number-conserving method. The level separation and, in particular, the level order change between 251 Lr and 255 Lr was reproduced only when the hexacontetrapole deformation ϵ 6 was included in the calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Production, isolation, and shipment of clinically relevant quantities of astatine-211: A simple and efficient approach to increasing supply

The alpha emitter astatine-211 ( 211 At) is a promising candidate for cancer treatment based on Targeted Alpha (α) Therapy (TAT). A small number of facilities, distributed across the United States, are capable of accelerating α-particle beams to produce 211 At. However, challenges remain regarding strategic methods for shipping 211 At in a form adaptable to advanced radiochemistry reactions and other uses of the radioisotope. Purpose: Our method allows shipment of 211 At in various quantities in a form convenient for further radiochemistry. Procedures: For this study, a 3-octanone impregnated Amberchrom® CG300M resin bed in a column cartridge was used to separate 211 At from the bismuth matrix on site at the production accelerator (Texas A&M) in preparation for shipping. Aliquots of 6 M HNO 3 containing up to ≈2.22 GBq of 211 At from the dissolved target were successfully loaded and retained on columns. Exempt packages (<370 MBq) were shipped to a destination radiochemistry facility, University of Texas MD Anderson Cancer Center, in the form of a convenient air-dried column. Type A packages have been shipped overnight to University of Alabama at Birmingham. Main findings: Air-dried column hold times of various lengths did not inhibit simple and efficient recovery of 211 At. Solution eluted from the column was sufficiently high in specific activity to successfully radiolabel a model compound, 4-(4,4,5,5-tetramethyl-1,3,2-dioxaborolan-2-yl)aniline (1), with 211 At. The method to prepare and ship 211 At described in this manuscript has also been used to ship larger quantities of 211 At a greater distance to University of Alabama at Birmingham. Principal conclusions: The successful proof of this method paves the way for the distribution of 211 At from Texas A&M University to research institutions and clinical oncology centers in Texas and elsewhere. Finally, use of this simple method at other facilities has the potential increase the overall availability of 211 At for preclinical and clinical studies.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Diagnosing hot-spot symmetry in surrogate ignition experiments via secondary DT-neutron spectroscopy at the NIF

The directional energy spectrum of neutrons generated from the in-flight fusion reaction of 1-MeV tritons contains information about the hot-spot symmetry. The National Ignition Facility (NIF) fields Symmetry Capsule (Symcap) implosions, which have historically measured the symmetry of the radiation, drive by measuring the hot-spot shape via x-ray self-emission. Symcaps are used to tune the hot-spot symmetry for ignition experiments at the NIF. This work shows the relationship between directional secondary DT-n spectra and x-ray imaging data for a large database of Symcap implosions. A correlation is observed between the relative widths of the DT-n spectra measured with nTOFs and the shape measured with x-ray imaging. A Monte Carlo model, which computes the directional secondary DT-n spectrum, is used to interpret the results. A comparison of the x-ray and secondary DT-n data with the Monte Carlo model indicates that 56% of the variance between the two datasets is explained by a P2 asymmetry. More advanced simulations using HYDRA suggest that the unaccounted variance is due to P1 and P4 asymmetries present in the hot spot. The comparison of secondary DT-n data and x-ray imaging data to the modeling shows the DT-n data contain important information that supplements current P2 measurements and contain new information about the P1 asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Basic Energy Sciences Roundtable: Foundational Science to Accelerate Nuclear Energy Innovation

Energy security, availability, and reliability are among the greatest challenges facing the nation and the planet. An abundant potential source of energy resides in the fundamental atomic building blocks of the universe in the form of nuclear fission and fusion reactions. In fact, energy from nuclear fission currently provides the majority of the world’s zero-carbon electricity, and future fusion energy systems offer great promise; carbon-free nuclear energy technologies can be key to the world’s decarbonized energy future. Although contemporary fission systems use well-established technologies to supply safe and efficient baseload power, they could be more fuel efficient and less costly. Moving beyond massive light-water fission reactors to a variety of advanced nuclear systems—which will vary in size and operate in extremes of temperature, corrosivity, and other parameters—will place stringent conditions on materials and chemical systems. New demands will be placed on the coolants and solvents, the materials, and the monitoring tools used in these reactors. Fusion-based nuclear energy will require superior materials to withstand extremely high temperatures, plasma exposure, radiation damage, and implanted gases. The advantages associated with these new fission and fusion technologies will be realized only through continued advancements in the fundamental science underpinning our knowledge of the physics and chemistry of nuclear systems gained via improved experimental and computational methods. In July 2022, the U.S. Department of Energy’s Office of Basic Energy Sciences—in coordination with the Offices of Nuclear Energy, Fusion Energy Sciences, and Advanced Scientific Computing Research—held a virtual roundtable titled “Foundational Science to Accelerate Nuclear Energy Innovation” to discuss the scientific and technical barriers for advanced nuclear energy systems. Five priority research opportunities were identified to address these scientific and technical challenges and to accelerate progress toward the realization of next-generation fusion and fission energy systems. The foundational science gaps inhibiting the advancement of nuclear energy technologies are identified and tackled in five priority research opportunities. These opportunities pave the way to accelerate the development and ultimately the adoption of new nuclear energy systems. They include the fundamental aspects of ion-electron interactions; novel properties of next-generation coolants and solvents; interfacial dynamics, not only in solids, but in other aspects of nuclear reactors; novel operando and in situ monitoring and sensing; and artificial intelligence to accelerate condensed phases discovery. Building on the foundation established by previous Basic Energy Sciences workshops, these opportunities encompass recent advances in fundamental knowledge and focus on the experimental and computational methods needed to resolve major technical challenges for nuclear energy technologies. Through developing fundamental scientific insight as well as pushing the frontiers of modeling complex systems and probing the operation of materials and chemical systems in extreme environments, research motivated by the priorities identified here will further develop the promise, potential, and utilization of nuclear energy for a clean energy future.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Reconstructing Masses of Merging Neutron Stars from Stellar r -process Abundance Signatures

We report that neutron star mergers (NSMs) are promising astrophysical sites for the rapid neutron-capture ("r") process, but can their integrated yields explain the majority of heavy-element material in the Galaxy? One method to address this question implements a forward approach that propagates NSM rates and yields along with stellar formation rates and compares those results with observed chemical abundances of r-process-rich, metal-poor stars. In this work, we take the inverse approach by utilizing r-process-element abundance ratios of metal-poor stars as input to reconstruct the properties -especially the masses - of their neutron star (NS) binary progenitors. This novel analysis provides an independent avenue for studying the population of the original NS binary systems that merged and produced the r-process material now incorporated in Galactic metal-poor halo stars. We use ratios of elements typically associated with the limited-r-process and the actinide region to those in the lanthanide region (i.e., Zr/Dy and Th/Dy) to probe the NS masses of the progenitor merger. We find that NSMs can account for all r-process material in metal-poor stars that display r-process signatures, while simultaneously reproducing the present-day distribution of double-NS systems. Notably, with our model assumptions and the studied stellar sample, we postulate that the most r-process enhanced stars (the r–II stars) on their own would require progenitor NSMs of asymmetric systems that are distinctly different from present ones in the Galaxy. We also explore variations to the model and find that the predicted degree of asymmetry is most sensitive to the electron fraction of the remnant disk wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Thermodynamic Modeling of Complex Solid Solutions in the Lu-H-N System via Graph Neural Network Accelerated Monte Carlo Simulations

Metal hydrides are important across diverse applications, such as hydrogen storage, batteries, gas sensors, nuclear reactions, and high-temperature superconductivity. Previous computational studies of metal hydrides under extreme pressures, e.g., 𝑂⁡(10 2 ) ⁢GPa, usually treat them as stoichiometric compounds without considering interstitial lattice disorder. As pressures become more moderate in the 𝑂⁡(10 0 ) ⁢GPa and below range, hydrogen disorder at interstitial lattice sites becomes prominent, e.g., in the N-doped Lu hydride that was recently claimed superconducting near 1 GPa. Further adding compositional complexity from alloying and/or multielement interstitial occupation makes elucidating pressure- and temperature-dependent observables intractable by first-principles calculations alone. We therefore propose a lattice graph neural-network surrogate modeling approach to predict configuration- and pressure-dependent equation-of-state properties. Their efficiency permits Monte Carlo simulations to calculate Gibbs energies and pressure-dependent phase diagrams, thereby revealing insights into the synthesis conditions required for achieving desired phase equilibria. We demonstrate this concept for the compositionally complex cubic Lu(H,N,Va) 3 system where three constituents (hydrogen, nitrogen and vacancy) have disordered multielement interstitial occupancies and insights into pressure-dependent phase equilibria are critically needed, e.g., N-doping levels can significantly lower dehydrogenation temperatures and provide a new strategy to optimize hydrogen-storage alloys. This work can improve the thermodynamic understanding of the Lu-H-N system and help rational synthesis of N-doped Lu hydrides, but more generally demonstrates an efficient approach to model pressure-dependent thermodynamics of multicomponent solid solutions.

Monte Carlo methods↗
DOE OSTI2020

Mass-asymmetric fission of 205,207,209 Bi at energies close to the fission barrier using proton bombardment of 204,206,208 Pb

Background: Recent observation of mass-asymmetric fission in neutron-deficient Hg and Pt nuclei has reignited interest in fission fragment mass distributions close to Pb. Investigations at energies close to the fission barrier, where mass-asymmetric fission is expected to be most obvious and the sensitivity to shell effects is maximized, are limited in this mass region. Purpose: To measure fission mass distributions for 205,207,209 Bi nuclei at the lowest possible excitation energies to determine how the mass distributions change with excitation energy and the neutron number of the compound nucleus. Method: Proton beams bombarding targets of 204,206,208 Pb were used to study the fission of 205,207,209 Bi at energies from just above to 10 MeV above their fission barriers. Fission fragments were measured using the CUBE fission spectrometer. Fission fragment mass distributions were determined using a newly developed time difference analysis method. Mass distributions were characterized by triple-Gaussian fits to determine the systematic trends across each isotope with excitation energy. Results: Measured mass distributions of all three Bi isotopes exhibit a component of mass-asymmetric fission at all energies studied. The probability of mass-asymmetric fission decreases significantly with increasing excitation energy, from ≈70 to ≈40% over a 10-MeV range. Comparisons between the three Bi isotopes hint at an increase in the mass-symmetric fission yield with increasing neutron number, which could be due to a decrease in the difference between the symmetric and asymmetric fission barriers. The centroids of the mass-asymmetric peaks suggest that several deformed shell gaps in the fission fragments could be contributing to the presence of the mass-asymmetric fission mode with Z light ≃ 38, Z heavy ≃ 45, and N light ≃ 56 all present in the fission fragments. Conclusions: Measurements of fission mass distributions at the lowest possible excitation energies above the fission barrier provide an excellent platform to investigate the origins of the mass-asymmetric fission mode. Finally, further systematic measurements at these energies offer an opportunity to rigorously test new models of fission in this mass region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗