SEARCH · Search NASA
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.
The current state of neutron star structure calculations is reviewed. Uncertainties in the equation of state for matter at and above nuclear density remain. The role of the delta resonance, pion condensates, and quark matter is reviewed. It is found that recent models yield stable neutron star masses which are consistent with observational estimates.
A comprehensive 31 P nuclear magnetic resonance (NMR) study, combined with thermodynamic measurements and first-principles band-structure calculations, has been conducted to explore the ground state of the 𝑆=5/2 double trillium lattice antiferromagnet KSrFe 2 (PO 4 ) 3 . Our experimental results indicate that the magnetic ground state is neither a conventional three-dimensional (3D) long-range order (LRO) nor a pure gapless spin-liquid state, as conjectured previously [Boya et al., APL Mater. 10, 101103 (2022)]. Specifically, the observation of a nearly field-independent NMR linewidth below 𝑇* = (3.5 ± 0.4) K, and a significant enhancement of spin-spin relaxation rate 1/𝑇 2 below 2𝑇* (where 𝑇* is the characteristic temperature identified from the magnetic susceptibility), indicate a complex magnetic ground state where spin freezing coexists with persistent dynamics. Furthermore, we argue that the lack of magnetic LRO and the persistence of strong magnetic fluctuations in KSrFe 2 (PO 4 ) 3 are unlikely to originate from intersite K/Sr disorder, rather they arise due to intrinsic magnetic frustration. In conclusion, our findings position KSrFe 2 (PO 4 ) 3 into a broader family of geometrically frustrated magnets characterized by coexisting spin freezing and pronounced antiferromagnetic fluctuations, marking it as a promising platform for investigating exotic phenomena in 3D frustrated magnets.
The available world deep-inelastic scattering (DIS) data on proton and deuteron structure functions F 2 p , F 2 d , and their ratios are leveraged to extract the free neutron F 2 n structure function, the F 2 n / F 2 p ratio, and associated uncertainties using the latest nuclear effect calculations in the deuteron. Special attention is devoted to the normalization of the proton and deuteron experimental datasets and to the treatment of correlated systematic errors, as well as the quantification of procedural and theoretical uncertainties. The extracted F 2 n dataset is utilized to evaluate the Q 2 dependence of the Gottfried sum rule and the nonsinglet F 2 p − F 2 n moments. To facilitate replication of our study, as well as for general applications, we provide a comprehensive DIS database including all recent Jefferson Lab 6 GeV measurements, the extracted F n 2 , a modified CTEQ-JLab global parton distribution function fit named CJ15nlo_mod, and grids with calculated proton, neutron, and deuteron DIS structure functions. Published by the American Physical Society 2024
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.
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.
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.
Hartree-Fock-Slater calculation of internal conversion coefficients for magnetic multipoles for yttrium-87 with 0.05 or 0.15 mc-square gamma energy values
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.
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.
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.
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.
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.
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.
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.
Procedure for calculating rates of non-resonant nuclear reactions at stellar temperatures
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.
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.
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.