Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nuclear structure calculations”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Hyperfine spectroscopy and laser cooling of the fermionic isotopes 47 Ti and 49 Ti

Here, we report on magneto-optical trapping of the two fermionic isotopes of atomic titanium, 47 Ti and 49 Ti . Unlike the even mass-number isotopes, which were recently laser cooled, 47 Ti and 49 Ti have nonzero nuclear spins and, consequently, their atomic levels are split by hyperfine structure. Combining and comparing theoretical calculations and atomic beam-spectroscopy measurements, we determine the hyperfine structures and isotope shifts of the 3⁢𝑑 2 ⁢4⁢𝑠 2 𝑎 3 ⁢𝐹 4 → 3⁢𝑑 2 ⁢( 3 𝑃)⁢4⁢𝑠⁢4⁢𝑝⁢( 3 𝑃 𝑜 ) 𝑦 5 ⁢𝐷$^o_4$ optical-pumping transition at optical wavelength 391 nm and the 3⁢𝑑 3 ⁢(4𝐹)⁢4⁢𝑠 𝑎 5 ⁢𝐹 5 → 3⁢𝑑 3 ⁢( 4 𝐹)⁢4⁢𝑝 𝑦 5 ⁢𝐺$^𝑜_6$ laser-cooling transition at wavelength 498 nm. With this information, we produce magneto-optical traps of both 47 Ti and 49 Ti by applying two additional tones of light to repump atoms to the maximum-spin states on the laser-cooling transition. Directly loading from the atomic flux of a titanium sublimation pump, we produce 47 Ti and 49 Ti traps with 731(190) and 1142(240) atoms, and with lifetimes of 330⁢(15) and 310⁢(8)⁢ ms, respectively.

atomic spectra↗
DOE OSTI2026

Neutron skins probed in proton knockout from neutron-rich nuclei

Proton-induced quasifree knockout reactions provide a powerful probe of nuclear single-particle structure and reaction dynamics in both stable and neutron-rich nuclei. Here, in this work, we develop a unified theoretical framework for the calculation of inclusive (𝑝, 2⁢𝑝) and sequential (𝑝, 3⁢𝑝) reaction cross sections and fragment momentum distributions at intermediate and relativistic energies. The approach is based on a probabilistic extension of Glauber multiple-scattering theory combined with microscopic nuclear densities obtained from Hartree-Fock-Bogoliubov calculations using Skyrme energy-density functionals. We focus in particular on the sensitivity of total cross sections and longitudinal momentum dispersions to neutron-skin thickness along isotopic chains. Our results indicate that both (𝑝, 2⁢𝑝) and (𝑝, 3⁢𝑝) reactions exhibit a systematic decrease of cross section and momentum width with increasing neutron excess, reflecting enhanced attenuation and surface bias induced by neutron skins. The effect is significantly stronger for two-proton removal, suggesting that (𝑝, 3⁢𝑝) reactions may offer enhanced sensitivity to isovector nuclear structure. These findings establish proton-induced knockout reactions as complementary hadronic probes of neutron skins and the density dependence of the nuclear symmetry energy.

direct reactions↗
DOE OSTI2026

Low-energy 17 O(𝑛,𝛾)⁢ 18 O reaction within the microscopic potential model and its role for the weak 𝑟 process

The neutron radiative capture reaction 17 O ⁡(𝑛,𝛾) ⁢18 O plays a pivotal role in both nuclear structure studies and astrophysical nucleosynthesis, particularly in the formation of elements during hydrostatic and explosive stellar environments. We calculated the 17 O ⁡(𝑛,𝛾) ⁢18 O cross section within the Skyrme Hartree-Fock potential model and analyzed electric dipole 𝐸⁢1 transitions to both positive- and negative-parity states below the α-decay threshold in 18 O. Our cross sections are significantly different from the data available in commonly used libraries. We further investigate the impact of the new calculated cross section on weak 𝑟-process nucleosynthesis using large-scale reaction network calculations across a wide range of electron fractions and entropies. Our results show that the 17 O ⁡(𝑛,𝛾) ⁢18 O reaction rate significantly influences the production of first 𝑟-process peak elements, such as strontium, under specific astrophysical conditions. This study highlights the importance of accurate nuclear data for light isotopes in modeling heavy-element synthesis and provides updated reaction rates for future nucleosynthesis simulations.

6 ≤ A ≤ 19↗
NASA NTRS1985

Photo nuclear energy loss term for muon-nucleus interactions based on xi scaling model of QCD

Extensive air showers (EMC) experiments discovered a significant deviation of the ratio of structure functions of iron and deuteron from unity. It was established that the quark parton distribution in nuclei are different from the corresponding distribution in the nucleus. It was examined whether these results have an effect on the calculation of photo nucleus energy loss term for muon-nucleus nuclear interaction. Though the EMC and SLAC data were restricted to rather large q sq region it is expected that the derivation would persist even in the low q sq domain. For the ratio of iron and deuteron structure function a rather naive least square fit of the form R(x) = a + bx was taken and it is assumed that the formula is valid for the whole q sq region the absence of any knowledge of R(x) for small q sq.

Roychoudhury, R.↗
DOE OSTI2025

Structure and dynamics of open-shell nuclei from spherical coupled-cluster theory

We extend the spherical coupled-cluster ab initio method for open-shell nuclei where two nucleons are removed from a shell subclosure. Following the recent implementation of the two-particle-attached approach [Phys. Rev. C 110, 044306 (2024)], we focus on the two-particle-removed method. Using the equations-of-motion framework, we address both nuclear structure and dipole response functions by coupling coupled-cluster theory with the Lorentz integral transform technique. We perform calculations using chiral interactions, including three-nucleon forces, and estimate many-body uncertainties by comparing different coupled-cluster truncation schemes. Here, we validate our approach by studying ground-state energies, excited states, and electric dipole polarizabilities in the oxygen and calcium isotopic chains. For binding energies and selected low-lying excited states, we achieve an accuracy comparable to that of the established closed-shell coupled-cluster theory and generally agree with experiment. Finally, we underestimate experimental data for electric dipole polarizabilities, particularly in calcium isotopes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Imaging nuclear shape through anisotropic and radial flow in high-energy heavy-ion collisions

Most atomic nuclei exhibit ellipsoidal shapes characterized by quadrupole deformation β 2 and triaxiality γ, and sometimes even a pear-like octupole deformation β 3 . The STAR experiment introduced a new ‘imaging-by-smashing’ technique to image the nuclear global shape by colliding nuclei at ultra-relativistic speeds and analyzing outgoing debris. Features of nuclear shape manifest in collective observables like anisotropic flow vn and radial flow via mean transverse momentum [p T ]. We present new measurements of the variances of v n (n = 2, 3, and 4) and [p T ], and the covariance of $v$$^{2}_{n}$ with [p T ], in collisions of highly deformed 238 U and nearly spherical 197 Au. Ratios of these observables between the two systems effectively suppress common final-state effects, isolating the strong impact of uranium’s deformation. By comparing results with state-of-the-art hydrodynamic model calculations, we extract β 2U and γ U values consistent with those deduced from low-energy nuclear structure measurements. Measurements of v 3 and its correlation with [p T ] also provide the first experimental suggestion of a possible octupole deformation for 238 U. These findings provide significant support for using high-energy collisions to explore nuclear shapes on femtosecond timescales, with implications for both nuclear structure and quark-gluon plasma studies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Impact of T - and ρ -dependent decay rates and new (n, γ ) cross-sections on the s process in low-mass asymptotic giant branch stars

Aims. We study the impact of nuclear input related to weak-decay rates and neutron-capture reactions on predictions for the slow neutron-capture process (s process) in asymptotic giant branch (AGB) stars. We provide the first database of surface abundances and stellar yields of the isotopes heavier than iron from the Monash models. Methods. We ran nucleosynthesis calculations with the Monash post-processing code for seven stellar structure evolution models of low-mass AGB stars with three different sets of nuclear inputs. The reference set has constant decay rates and represents the set used in the previous Monash publications. The second set contains the temperature and density dependence of β decays and electron captures based on the default rates of nuclear NETwork GENerator (NETGEN). In the third set, we further update 92 neutron-capture rates based on re-evaluated experimental cross sections from the ASTrophysical Rate and rAw data Library. We compare and discuss the predictions of the sets relative to each other in terms of isotopic surface abundances and total stellar yields. We also compare the results to isotopic ratios measured in presolar stardust silicon carbide (SiC) grains from AGB stars. Results. The new sets of models result in a ∼66% solar s-process contribution to the p-nucleus 152 Gd, confirming that this isotope is predominantly made by the s process. The nuclear input updates result in predictions for the 80 Kr/ 82 Kr ratio in the He intershell and surface 64 Ni/ 58 Ni, 94 Mo/ 96 Mo, and 137 Ba/ 136 Ba ratios that are more consistent with the corresponding ratios measured in stardust; however, the new predicted 138 Ba/ 136 Ba ratios are higher than the typical values of the SiC grains. The W isotopic anomalies are in agreement with data from the analyses of other meteoritic inclusions. We confirm that the production of 176 Lu and 205 Pb is affected by too large uncertainties in their decay rates from NETGEN.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Weighted active space protocol for multireference machine-learned potentials

Multireference methods such as multiconfiguration pair-density functional theory accurately capture electronic correlation in systems with strong multiconfigurational character, but their cost precludes direct use in molecular dynamics. Combining these methods with machine-learned interatomic potentials (MLPs) can extend their reach. However, the sensitivity of multireference calculations to the choice of the active space complicates the consistent evaluation of energies and gradients across structurally diverse nuclear configurations. To overcome this limitation, we introduce the weighted active space protocol (WASP), a systematic approach to assign a consistent active space for a given system across uncorrelated configurations. By integrating WASP with MLPs and enhanced sampling techniques, we propose a data-efficient active learning cycle that enables the training of an MLP on multireference data. We demonstrated the approach on the TiC + -catalyzed C–H activation of methane, a reaction that poses challenges for Kohn–Sham density functional theory due to its significant multireference character. This framework enables accurate and efficient modeling of catalytic dynamics, establishing a paradigm for simulating complex reactive processes beyond the limits of conventional electronic-structure methods.

enhanced sampling↗
DOE OSTI2026

Halo Nuclei from Ab Initio Nuclear Theory

A realistic description of halo nuclei, characterized by low-lying breakup thresholds, requires a proper treatment of continuum effects. We have developed an ab initio approach, the No-Core Shell Model with Continuum (NCSMC), capable of describing both bound and unbound states in light nuclei in a unified way. With chiral two- and three-nucleon interactions as the only input, we can predict the structure and dynamics of halo and other light nuclei and, by comparing to available experimental data, test the quality of chiral nuclear forces. We review NCSMC calculations of weakly bound states and resonances of the exotic halo nuclei 6He, 8B, 11Be, and 15C. For the latter, we discuss its production in the capture reaction 14C(n,𝛾 )15C. We highlight the challenges of a description of 6He as a Borromean n-n-4He system. Finally, we present our calculations of excited states in 10Be exhibiting a one-neutron halo structure and a large scale No-Core Shell Model investigation of 11Li as a precursor of a full n-n-9Li NCSMC study.

Navrátil, Petr↗
DOE OSTI2024

Multipole responses in fissioning nuclei and their uncertainties

Electromagnetic multipole responses are key inputs to model the structure, decay, and reactions of atomic nuclei. With the introduction of the finite amplitude method (FAM), large-scale calculations of the nuclear linear response in heavy deformed nuclei have become possible. This work provides a detailed study of multipole responses in actinide nuclei with Skyrme energy density functionals. We quantify both systematic and statistical uncertainties induced by the functional parametrization in FAM calculations. We also extend the FAM formalism to perform blocking calculations with the equal filling approximation for odd-mass and odd-odd nuclei, and analyze the impact of blocking configurations on the response. In conclusion, by examining the entire plutonium isotopic chain from the proton to the neutron drip line, we find a large variability of the response with the neutron number and study how it correlates with the deformation of the nuclear ground state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Pre-Equilibrium De-Excitations for Neutrino-Nucleus Interactions

The Deep Underground Neutrino Experiment (DUNE) is sensitive to MeV-scale energy depositions from low-energy astrophysical neutrinos, including those from core-collapse supernovae. Interpreting these detector signals requires accurate modeling of the nuclear de-excitation that follows from the neutrino-nucleus interaction. The MARLEY (Model of Argon Reaction Low-Energy Yields) event generator specializes in the low-energy regime. MARLEY currently assumes the residual nucleus equilibrates immediately after the primary interaction. This omits the intermediate pre-equilibrium stage in which energy redistributes among nucleons until statistical equilibrium is reached. While pre-equilibrium effects are well established for nucleon-induced reactions, they have not previously been studied for neutrino-nucleus interactions. This work addresses that gap by implementing a two-component exciton model, which is the first dedicated treatment of pre-equilibrium de-excitation for neutrino-nucleus interactions, restructured around MARLEY's existing class hierarchy to prepare for direct integration, including particle-hole state densities, internal transition rates, and pre-equilibrium particle emission. The calculations show encouraging agreement with the TALYS-2.2 nuclear reaction code for neutron-nucleus interactions. We further propose a concrete integration path into the full MARLEY event generator, including derived class structure and an extended event record for pre-equilibrium vertices in support of future reweighting. Remaining work focuses on refining the emission width calculation, adding $\gamma$-ray emission, and completing this integration to quantify the impact of pre-equilibrium effects on the expected low-energy neutrino signals in DUNE and similar experiments.

Visser, Erin [Michigan State U., East Lansing (mai↗
DOE OSTI2025

CONFLUX: A standardized framework to calculate reactor antineutrino flux

Nuclear fission reactors are abundant sources of antineutrinos for neutrino physics experiments. The flux and spectrum of antineutrinos emitted by a reactor can indicate its activity and composition, suggesting potential applications of neutrino measurements beyond fundamental scientific studies that may be valuable to society. The utility of reactor antineutrinos for applications and fundamental science is dependent on the availability of precise predictions of these emissions. For example, in the last decade, disagreements between reactor antineutrino measurements and models have inspired revision of reactor antineutrino calculations and standard nuclear databases as well as searches for new fundamental particles not predicted by the Standard Model of particle physics. Past predictions and descriptions of the methods used to generate them are documented to varying degrees in the literature, with different modeling teams incorporating a range of methods, input data, and assumptions. The resulting difficulty in accessing or reproducing past models and reconciling results from differing approaches complicates the future study and application of reactor antineutrinos. The CONFLUX (Calculation Of Neutrino FLUX) software framework is a neutrino prediction tool built with the goal of simplifying, standardizing, and democratizing the process of reactor antineutrino flux calculations. CONFLUX includes three primary methods for calculating the antineutrino emissions of nuclear reactors or individual beta decays that incorporate common nuclear data and beta decay theory. The software is prepackaged with the current nuclear databases, including ENDF.B/VIII, JEFF-3.3, and ENSDF, and it includes the capability to predict time-dependent reactor emissions, adjust nuclear database or beta decay inputs/assumptions, and propagate related sources of uncertainty. Here, this paper describes the CONFLUX software structure, details the methods used for flux and spectrum calculations, and provides examples of potential use cases.

Zhang, Xianyi [Lawrence Livermore National Laborat↗
DOE OSTI2024

Structure Factors for Hot Neutron Matter from Ab Initio Lattice Simulations with High-Fidelity Chiral Interactions

We present the first ab initio lattice calculations of spin and density correlations in hot neutron matter using high-fidelity interactions at next-to-next-to-next-to-leading order in chiral effective field theory. These correlations have a large impact on neutrino heating and shock revival in core-collapse supernovae and are encapsulated in functions called structure factors. Unfortunately, calculations of structure factors using high-fidelity chiral interactions were well out of reach using existing computational methods. In this Letter, we solve the problem using a computational approach called the rank-one operator (RO) method. The RO method is a general technique with broad applications to simulations of fermionic many-body systems. It solves the problem of exponential scaling of computational effort when using perturbation theory for higher-body operators and higher-order corrections. Using the RO method, we compute the vector and axial static structure factors for hot neutron matter as a function of temperature and density. Here, the ab initio lattice results are in good agreement with virial expansion calculations at low densities but are more reliable at higher densities. Random phase approximation codes used to estimate neutrino opacity in core-collapse supernovae simulations can now be calibrated with ab initio lattice calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

𝑙-forbidden 𝑀⁢1 strengths near 100 Sn from knockout reactions in Cd and Sn

Neutron knockout reactions on beams of 104,102 Cd and 104 Sn are presented. States in the residual 103,101 Cd and 103 Sn nuclei are populated, including low-lying 7/2 + states of 𝜈⁢𝑔 7/2 character. These states have half-lives ≈ 400 ps due to their low energy and hindered 𝐵⁡(𝑀⁢1;7/2 + →5/2 + ) strengths. The excited-state half-lives were measured using their Doppler-shifted line shapes, and the resulting 𝐵⁡(𝑀⁢1) strengths are compared to valence space in medium similarity renormalization group (VS-IMSRG) calculations. Lastly, the VS-IMSRG calculations underpredict the 𝑙-forbidden 𝑀⁢1 strengths in the 100 Sn region, as well as in other regions of the nuclear chart near 40 Ca and 208 Pb.

90 ≤ A ≤ 149↗
NASA NTRS2009

Using FLUKA to Calculate Spacecraft: Single Event Environments: A Practical Approach

The FLUKA nuclear transport and reaction code can be developed into a practical tool for calculation of spacecraft and planetary surface asset SEE and TID environments. Nuclear reactions and secondary particle shower effects can be estimated with acceptable accuracy both in-flight and in test. More detailed electronic device and/or spacecraft geometries than are reported here are possible using standard FLUKA geometry utilities. Spacecraft structure and shielding mass. Effects of high Z elements in microelectronic structure as reported previously. Median shielding mass in a generic slab or concentric sphere target geometry are at least approximately applicable to more complex spacecraft shapes. Need the spacecraft shielding mass distribution function applicable to the microelectronic system of interest. SEE environment effects can be calculated for a wide range of spacecraft and microelectronic materials with complete nuclear physics. Evaluate benefits of low Z shielding mass can be evaluated relative to aluminum. Evaluate effects of high Z elements as constituents of microelectronic devices. The principal limitation on the accuracy of the FLUKA based method reported here are found in the limited accuracy and incomplete character of affordable heavy ion test data. To support accurate rate estimates with any calculation method, the aspect ratio of the sensitive volume(s) and the dependence must be better characterized.

Koontz, Steve↗
NASA NTRS1987

An algorithm for a semiempirical nuclear fragmentation model

In the present paper, an algorithm for HZE (High Charge and Energy) fragmentation based upon a combination of a two step abrasion/ablation model and electromagnetic dissociation is presented. Development of the model and detailed comparison with available experimental data are given elsewhere. The abrasion process accounts for the removal of nuclear matter in the overlap region of the colliding ions. An average transmission factor is used for the projectile and target nuclei at a given impact parameter to account for the finite mean free path in nuclear matter. The ions are treated otherwise on a geometric basis assuming uniform spheres. The ablation process is treated as a single nucleon-evaporation for every 10 MeV of excitation energy as used by Bowman in the original form of the model. The charge distribution of final fragments are calculated according to the Rudstam formula, except for some correction in mass 5, 8 and 9 fragments which show strong structure effects and correspondingly significant deviation from Rudstam's values. The nuclear electromagnetic dissociation is based on the Weizsacker-Williams (WW) method of virtual quanta where due to its simplicity, the virtual photon spectrum for individual multipoles, and finite extent of the charge distribution are not included. Comparisons of the model are made with the available experimental data here and more extensively elsewhere.

NASA Program Space Medicine↗