SEARCH · Engineering Papers
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.
Here, we study (nuclear) fusion reactions in beryllium. We argue that some are nuclear long-distance molecular processes of interest for low-energy nuclear reactions. For these, we develop a novel reaction model. We note a chain of reactions that can naturally arise in beryllium target created by nonequilibrium proton reactions.
We present calculations of neutron thermal cross sections, Westcott factors, resonance integrals, Maxwellian-averaged cross sections, astrophysical reaction rates, and solar system r-process abundances using the latest data from the major evaluated nuclear libraries for 849 ENDF target materials. The recent release of ENDF/B-VIII.1 library, progress in 252 Cf(SF) evaluation, extensive analysis of newly-evaluated neutron reaction cross sections, neutron covariances, and improvements in data processing techniques motivated us to calculate the nuclear industry and neutron physics parameters, produce s-process Maxwellian-averaged cross sections and astrophysical reaction rates, extract r-process abundances, systematically calculate uncertainties, and provide additional insights on currently available neutron-induced reaction data.
The ground state half-lives of 69 Ge, 73 Se, 83 Sr, 63 Zn, and the half-life of the 1/2 - isomer in 85 Sr have been measured with high precision using the photoactivation technique at an unconventional bremsstrahlung facility that features a repurposed medical electron linear accelerator. The γ -ray activity was counted over about 6 half-lives with a high-purity germanium detector, enclosed into an ultra low-background lead shield. The measured half-lives are: T 1/2 ( 69 Ge) = 38.82 ± 0.07 (stat) ± 0.06 (sys) h; T 1/2 ( 73 Se) = 7.18 ± 0.02 (stat) ± 0.004 (sys) h; T 1/2 ( 83 Sr) = 31.87 ± 1.16 (stat) ± 0.42 (sys) h; T 1/2 ( 85m Sr) = 68.24 ± 0.84 (stat) ± 0.11 (sys) min; T 1/2 ( 63 Zn) = 38.71 ± 0.25 (stat) ± 0.10 (sys) min. These high-precision half-life measurements will contribute to a more accurate determination of corresponding ground-state photoneutron reaction rates, which are part of a broader effort of constraining statistical nuclear models needed to calculate stellar nuclear reaction rates relevant for the astrophysical p -process nucleosynthesis.
Type Ia supernovae are bright thermonuclear explosions that are important to numerous areas of astronomy. However, the origins of these events are poorly understood. One proposed setting is that of a near Chandrasekhar mass white dwarf that undergoes runaway carbon burning in the core. During the thousand years leading up to the explosion, the white dwarf undergoes a simmering phase where slow carbon burning heats the core and drives convection. A poorly understood aspect of this phase is the convective Urca process, which links convection with weak nuclear reactions. We use the low Mach number code MAESTROeX to perform full 3D simulations as is required to accurately capture the turbulent convection. We present simulations with and without the A=23 convective Urca process, which have relaxed to a steady state. We characterize the effects of the convective Urca process on the neutrino losses, the nuclear energy generation, and the convective boundary. We find that the size of the convection zone is substantially reduced by the convective Urca process, though convection still extends past the Urca shell. Our findings on the structure of the convective zone and the compositional changes can be used to inform 1D stellar models that track the longer-timescale evolution.
We address the calculation of microcanonical reaction rates for processes involving significant nuclear quantum effects using ring-polymer molecular dynamics (RPMD), both with and without electronically non-adiabatic transitions. After illustrating the shortcomings of the naive free-particle direct-shooting method, in which the temperature of the internal ring-polymer modes is set to the translational energy scale, we investigate alternative strategies based on the expression for the microcanonical rate in terms of the inverse Laplace transform of the thermal reaction rate. It is shown that simple application of the stationary-phase approximation (SPA) dramatically improves the performance of the microcanonical rates using RPMD, particularly in the low-energy region where tunneling dominates. Using the SPA as a Bayesian prior, numerically exact RPMD microcanonical rates are then obtained using maximum entropy inversion of the thermal reaction rates for both electronically adiabatic and non-adiabatic model systems. Finally, the direct-shooting method is revisited using the SPA-determined temperature for the internal ring-polymer modes, leading to a simple, direct-simulation method with improved accuracy in the tunneling regime. This work suggests a general strategy for the extraction of microcanonical dynamical quantities from RPMD (or other approximate thermal) simulations.
Advanced reprocessing schemes such as the Zircex process are important for continued development of the nuclear fuel cycle. The hydrochlorination reaction of the Zircex process flow sheet for metallic fuels has a known off-gas stream of H 2 (g) and HCl. The flow sheet can be simplified through HCl recycling, which requires detection and quantification of any residual H 2 (g) for safety and purity. In this work, commercially available Raman spectroscopy was applied for low-level, on-line detection of H 2 (g) in process gas streams, and parameters were adjusted to optimize the H 2 (g) Raman signal. Increasing the number of scans and exposure time in the spectrometer increased the signal-to-noise ratio of the H 2 (g) Raman signal, while the gas flow rate was optimized at 2.3 L/min. The limit of detection for H 2 (g) was estimated to be 2,100 ppm H 2 (g) in an N 2 (g) background. This scoping study for the application of Raman spectroscopy for on-line measurements of H 2 (g) shows that H 2 (g) can be detected at low levels (approximately 5% of the lower explosive limit (LEL)) in flowing gas streams.
R-matrix formalism is extended beyond compound nuclear (CN) resonant reactions to include parameterization of direct as well as doorway processes. Direct processes in the R-matrix exterior are parameterized by a unitary matrix that introduces mixing among wave function coefficients of the incoming and outgoing wave function components at the R-matrix channel surface. Doorway processes are parameterized by separating the Hilbert space of the interior R-matrix region into its doorway and CN subspaces, from which doorway state eigenenergies, reduced width amplitudes, and the strengths of their coupling to CN levels appear as new R-matrix parameters. Parameterization of generalized as well as the conventional Reich–Moore approximation for eliminated capture channels in the presence of direct, doorway, and CN processes is presented along with a complex-valued scattering length with contributions from direct, doorway, and CN capture processes. Derivation of Brune’s alternative R-matrix parameters is extended to include doorway states. This work suggests how R-matrix formalism could be extended further by adopting the concepts from related reaction formalisms.
Abstract Nuclear reactions heat and cool the crust of accreting neutron stars and need to be understood to interpret observations of X-ray bursts and long-term cooling in transiently accreting systems. It was recently suggested that previously ignored neutron transfer reactions may play a significant role in the nuclear processes. We present results from full nuclear network calculations that now include these reactions and determine their impact on crust composition, crust impurity, heating, and cooling. We find that a large number of neutron transfer reactions indeed occur and impact crust models. In particular, we identify a new type of reaction cycle that brings a pair of nuclei across the nuclear chart into equilibrium via alternating neutron capture and neutron release, interspersed with a neutron transfer. While neutron transfer reactions lead to changes in crust model predictions and need to be considered in future studies, previous conclusions concerning heating, cooling, and compositional evolution are remarkably robust.
Solar flares studied in the gamma ray region provide essential information on accelerated nuclei that can be obtained in no other way. A multitude of physical processes, such as particle acceleration, nuclear reactions, positron and neutron physics, and kinematical line broadening, come into consideration at gamma ray energies. Gamma ray observations are complementary to hard X ray observations, since both provide information on accelerated particles. It appears that only in the gamma ray region do these particles produce distinct spectral lines.
New astronomical observations point to a nucleosynthesis picture that goes beyond what was accepted until recently. The intermediate “i” process was proposed as a plausible scenario to explain some of the unusual abundance patterns observed in metal-poor stars. The most important nuclear physics properties entering i-process calculations are the neutron-capture cross sections and they are almost exclusively not known experimentally. In this report we provide the first experimental constraints on the 139 Ba(n,γ) 140 Ba reaction rate, which is the dominant source of uncertainty for the production of lanthanum, a key indicator of i-process conditions. This is an important step towards identifying the exact astrophysical site of stars carrying the i-process signature.
This study introduces advancements to the numerical solutions employed in the processing of nuclear data for fusion applications. It leverages the convolution theorem and Fourier transform techniques to enhance computational efficiency and broaden applicability. Building upon a previously reported discrete Hankel transform approach for Doppler broadening, this work refines the solution of convolution integrals central to these applications. The methodology provides a general and unified framework for evaluating any convolution operation, regardless of whether the underlying problem involves temperature effects in nuclear reactions. The applicability to the nuclear data processing for fusion is demonstrated by deriving the convolution integrals for some of the fusion-related quantities. As before, the convolution operation utilizes a Gaussian-based kernel; however, the discrete Hankel transform of order $𝛼$ = $\frac{1}{2}$ is now applied to the forward Fourier transform of the nonkernel argument, rather than the inverse Fourier transform. This modification eliminates the need for the integration of the nonkernel, cross section–based function, which is a step that posed challenges for certain pointwise cross-section representations. It also removes the requirement for cross-section linearization. Optimized for graphics processing unit architectures, the approach significantly improves computational performance. These advancements are currently under evaluation as the foundation for the next-generation thermonuclear data file processing codes being developed at Lawrence Livermore National Laboratory.
Abstract X-ray bursts are among the brightest stellar objects frequently observed in the sky by space-based telescopes. A type-I X-ray burst is understood as a violent thermonuclear explosion on the surface of a neutron star, accreting matter from a companion star in a binary system. The bursts are powered by a nuclear reaction sequence known as the rapid proton capture process (rp process), which involves hundreds of exotic neutron-deficient nuclides. At so-called waiting-point nuclides, the process stalls until a slower β + decay enables a bypass. One of the handful of rp process waiting-point nuclides is 64 Ge, which plays a decisive role in matter flow and therefore the produced X-ray flux. Here we report precision measurements of the masses of 63 Ge, 64,65 As and 66,67 Se—the relevant nuclear masses around the waiting-point 64 Ge—and use them as inputs for X-ray burst model calculations. We obtain the X-ray burst light curve to constrain the neutron-star compactness, and suggest that the distance to the X-ray burster GS 1826–24 needs to be increased by about 6.5% to match astronomical observations. The nucleosynthesis results affect the thermal structure of accreting neutron stars, which will subsequently modify the calculations of associated observables.
The competing 22 Ne(α, γ) 26 Mg and 22 Ne(α,n) 25 Mg reactions control the production of neutrons for the weak s process in massive and asymptotic giant branch (AGB) stars. In both systems, the ratio between the corresponding reaction rates strongly impacts the total neutron budget and strongly influences the final nucleosynthesis. A number of experimental studies have been performed over recent years which necessitate the reevaluation of the 26 Mg and 22 Ne(α, γ) 25 Mg reactions control the production 26 Mg and 22 Ne(α, n) 25 Mg reaction rates. Evaluations of the reaction rates following the collection of new nuclear data presently show differences of up to a factor of 500, resulting in considerable uncertainty in the resulting nucleosynthesis. To reevaluate the 22 Ne(α, γ) 26 Mg and 22 Ne(α,n) 25 Mg reaction rates using updated nuclear data from a number of sources including updating spin and parity assignments. With updated spin and parity assignments, the levels which can contribute to the reaction rates are identified. The reaction rates are computed using a Monte Carlo method which has been used for previous evaluations of the reaction rates in order to focus solely on the changes due to modified nuclear data. The evaluated 22 Ne(α, γ) 26 Mg reaction rate remains substantially similar to that of Longland et al. but, including recent results from Texas A&M, the 22 Ne(α, n) 25 Mg reaction rate is lower at a range of astrophysically important temperatures. Stellar models computed with newton and mesa predict decreased production of the weak branch s process due to the decreased efficiency of 22 Ne as a neutron source. Using the new reaction rates in the mesa model results in 96 Zr / 94 Zr and 135 Ba / 136 Ba ratios in much better agreement with the measured ratios from presolar SiC grains. The 22 Ne + α reaction rates 22 Ne(α, γ) 26 Mg and 22 Ne(α, n) 25 Mg have been recalculated based on more recent nuclear data. The 22 Ne(α, γ) 26 Mg reaction rate remains substantially unchanged since the previous evaluation but the 22 Ne(α, γ) 25 Mg reaction rate is substantially decreased due to updated nuclear data. This results in significant changes to the nucleosynthesis in the weak branch of the s process.
A discussion is given of the various roles that nuclear reactions may play in production of anomalous abundances of elements in peculiar A stars. The effects which may be expected to occur both in the surface nuclear reactions and in some possible internal reactions that can occur in advanced stages of stellar evolution are considered. It is suggested that various features of peculiar A stars may require simultaneous operation of two or more of the processes of surface diffusion, surface nuclear reactions, and internal nuclear reactions.
ABSTRACT Massive stars are crucial to galactic chemical evolution for elements heavier than iron. Their contribution at early times in the evolution of the Universe, however, is unclear due to poorly constrained nuclear reaction rates. The competing 17O(α, γ)21Ne and 17O(α, n)20Ne reactions strongly impact weak s-process yields from rotating massive stars at low metallicities. Abundant 16O absorbs neutrons, removing flux from the s-process, and producing 17O. The 17O(α, n)20Ne reaction releases neutrons, allowing continued s-process nucleosynthesis, if the 17O(α, γ)21Ne reaction is sufficiently weak. While published rates are available, they are based on limited indirect experimental data for the relevant temperatures and, more importantly, no uncertainties are provided. The available nuclear physics has been evaluated, and combined with data from a new study of astrophysically relevant 21Ne states using the 20Ne(d, p)21Ne reaction. Constraints are placed on the ratio of the (α, n)/(α, γ) reaction rates with uncertainties on the rates provided for the first time. The new rates favour the (α, n) reaction and suggest that the weak s-process in rotating low-metallicity stars is likely to continue up to barium and, within the computed uncertainties, even to lead.
We review the ab initio symmetry-adapted (SA) framework for determining the structure of stable and unstable nuclei, along with related electroweak, decay, and reaction processes. This framework utilizes the dominant symmetry of nuclear dynamics, the shape-related symplectic [Formula: see text] symmetry, which has been shown to emerge from first principles and to expose dominant degrees of freedom that are collective in nature, even in the lightest species or seemingly spherical states. This feature is illustrated for a broad scope of nuclei ranging from helium to titanium isotopes, enabled by recent developments of the ab initio SA no-core shell model expanded to the continuum through the use of the SA basis and that of the resonating group method. The review focuses on energies, electromagnetic transitions, quadrupole and magnetic moments, radii, form factors, and response function moments for ground-state rotational bands and giant resonances. The method also determines the structure of reaction fragments that is used to calculate decay widths and α-capture reactions for simulated X-ray burst abundance patterns, as well as nucleon–nucleus interactions for cross sections and other reaction observables.
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.
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.