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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Exploiting 20 Ne Isotopes for Precision Characterizations of Collectivity in Small Systems

Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of 16 O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing 16 O + 16 O data with collisions of 20 Ne ions. We couple both nuclear lattice effective field theory (NLEFT) and projected generator coordinate method (PGCM) ab initio descriptions of the structure of 20 Ne and 16 O to hydrodynamic simulations of 16 O + 16 O and 20 Ne + 20 Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of 20 Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of 20 Ne + 20 Ne collisions is enhanced by as much as 1.174⁢(8) stat ⁢(31) syst for NLEFT and 1.139⁢(6) stat ⁢(39) syst for PGCM relative to 16 O + 16 O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Toward scalable bound-to-resonance extrapolations for few- and many-body systems

In open quantum many-body systems, the theoretical description of resonant states of many particles strongly coupled to the continuum can be challenging. Such states are commonplace in, for example, exotic nuclei and hadrons, and can reveal important information about the underlying forces at play in these systems. In this work, we demonstrate that the complex-augmented eigenvector continuation (CA-EC) method, originally formulated for the two-body problem with uniform complex scaling, can reliably perform bound-to-resonance extrapolations for genuine three-body resonances having no bound subsystems. Here, we first establish that three-body bound-to-resonance extrapolations are possible by benchmarking different few-body approaches, and we provide arguments to explain how the extrapolation works in the many-body case. We furthermore pave the way towards scalable resonance extrapolations in many-body systems by showing that the CA-EC method also works in the Berggren basis, studying a realistic application using the Gamow shell model.

Ab initio calculations↗
DOE OSTI2024

Occurrence of Gravitational Collapse in the Accreting Neutron Stars of Binary-driven Hypernovae

The binary-driven hypernova (BdHN) model proposes long gamma-ray bursts (GRBs) originate in binaries composed of a carbon–oxygen (CO) star and a neutron star (NS) companion. The CO core collapse generates a newborn NS and a supernova that triggers the GRB by accreting onto the NSs, rapidly transferring mass and angular momentum to them. This article aims to determine the conditions under which a black hole (BH) forms from NS collapse induced by the accretion and the impact on the GRB's observational properties and taxonomy. We perform three-dimensional, smoothed particle hydrodynamics simulations of BdHNe using up-to-date NS nuclear equations of state, with and without hyperons, and calculate the structure evolution in full general relativity. We assess the binary parameters leading either NS in the binary to the critical mass for gravitational collapse into a BH and its occurrence time, t col . We include a nonzero angular momentum of the NSs and find that t col ranges from a few tens of seconds to hours for decreasing NS initial angular momentum values. BdHNe I are the most compact (about 5 minute orbital period), promptly form a BH, and release ≳10 52 erg of energy. They form NS–BH binaries with tens of kiloyears merger timescales by gravitational-wave emission. BdHNe II and III do not form BHs, and release ~10 50 –10 52 erg and ≲10 50 erg of energy, respectively. They form NS–NS binaries with a range of merger timescales larger than for NS–BH binaries. In some compact BdHNe II, either NS can become supramassive, i.e., above the critical mass of a nonrotating NS. Magnetic braking by a 10 13 G field can delay BH formation, leading to BH–BH or NS–BH with tens of kiloyears merger timescales.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2016

Collisional Dissociation of CO: ab initio Potential Energy Surfaces and Quasiclassical Trajectory Rate Coefficients

We have generated accurate global potential energy surfaces for CO+Ar and CO+O that correlate with atom-diatom pairs in their ground electronic states based on extensive ab initio electronic structure calculations and used these potentials in quasi-classical trajectory nuclear dynamics calculations to predict the thermal dissociation rate coefficients over 5000- 35000 K. Our results are not compatible with the 20-45 year old experimental results. For CO + Ar we obtain fairly good agreement with the experimental rate coefficients of Appleton et al. (1970) and Mick and Roth (1993), but our computed rate coefficients exhibit a stronger temperature dependence. For CO + O our dissociation rate coefficient is in close agreement with the value from the Park model, which is an empirical adjustment of older experimental results. However, we find the rate coefficient for CO + O is only 1.5 to 3.3 times larger than CO + Ar over the temperature range of the shock tube experiments (8000-15,000 K). The previously accepted value for this rate coefficient ratio is 15, independent of temperature. We also computed the rate coefficient for the CO + O ex- change reaction which forms C + O2. We find this reaction is much faster than previously believed and is the dominant process in the removal of CO at temperatures up to 16,000 K. As a result, the dissociation of CO is accomplished in two steps (react to form C+O2 and then O2 dissociates) that are endothermic by 6.1 and 5.1 eV, instead of one step that requires 11.2 eV to break the CO bond.

exchange rate coefficient↗
DOE OSTI2025

Hydrogen and water interactions with CrMnFeCoNi alloy from density functional theory calculations

High entropy alloys (HEAs) are a promising class of materials with remarkable mechanical and catalytic properties. Among these, the quinary CrMnFeCoNi alloy (also called “Cantor alloy”) has attracted considerable attention given its thermodynamic stability and remarkable mechanical properties under different temperatures. Given that various degradation mechanisms involve multiple contaminants, such as hydrogen and water in hydrogen embrittlement and surface poisoning, respectively, understanding their interactions with the Cantor alloy is critical for its practical applications as structural, nuclear, or hydrogen storage material. In this work, we perform first-principles calculations based on Density Functional Theory (DFT) to investigate such interactions when considering various microstructures, including bulk materials and those containing certain defects, such as grain boundaries, stacking faults, and vacancies. We also employ Global Sensitivity Analysis to identify the importance of different factors in the stability of the impurities. We find that the accuracy of the H formation energy is significantly affected by spin polarization and chemical short-range order. The study also identifies a strong tendency for hydrogen interstitials to segregate to Σ5(210)/[001] symmetric tilt grain boundary, even when H concentrations are high, suggesting that a certain type of grain boundaries acts as H sinks within the alloy. Further, this result is reinforced by the low formation energy of vacancy-hydrogen complexes, which can contain multiple hydrogen atoms. Finally, the surface reactivity analysis reveals that the adsorption energy of oxygen and hydroxyl groups is highly sensitive to the specific metal atom involved in the binding, with a clear preference for chromium atoms, which could have implications for the alloy’s oxidation and corrosion behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB + QRPA calculations with the Gogny interaction

The zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of-the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98 Zr . This region includes a benchmarking case and offers insights into nuclear deformation phenomena. Here, to investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the quasiparticle random phase approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98 Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the 𝐾 𝜋 = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the 𝐾 𝜋 = 0 − and 𝐾 𝜋 = 1 − components. Intrinsic transition densities and associated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e., we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Enabling event-by-event precision in γ-ray cascades for neutron-induced reactions

Neutron-induced γ-ray spectra provide key inputs for modern active interrogation applications. A precise modeling of the nuclear reaction and subsequent emission of γ rays is challenging and often impossible due to limitations on evaluated data file formats and nuclear transport simulation codes. We present a framework that addresses these challenges by combining experimental data and reaction-model calculation outputs into an extended candidate version of the Generalized Nuclear Data Structure (GNDS) file, the successor format for the legacy Evaluated Nuclear Data File (ENDF-6). This proposed GNDS hierarchical format contains all the necessary ingredients for inline γ-ray cascade reproduction with event-by-event precision, including continuum–continuum and continuum–discrete transitions following neutron-capture and inelastic neutron scattering reactions. Cascade-event generation based on our approach demonstrates improved energy conservation on an event-by-event basis and permits the use of γ-γ coincidences in applications. This work offers, for the first time, a method to generate neutron-capture and inelastic neutron-scattering γ-ray cascades where energy conservation, correlations, and experimental primaries are fully accounted for.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Specimen sizing and remaining life sample calculations for high temperature reactor material surveillance

Advanced high-temperature nuclear reactors rely on structural components that will operate for decades under combined mechanical, thermal, and environmental loading. Materials surveillance programs are a promising strategy for managing the resulting uncertainty in long-term structural integrity by monitoring degradation in service using passively actuated mechanical test articles. Previous reports have developed a simplified, spreadsheetimplementable framework for sizing these test articles and for inferring accumulated creep damage and remaining life from ex-situ test data. This report advances that work toward practical deployment by providing a sample-problem book: a collection of worked, end-toend examples in which the ASME Section III, Division 5 design analysis of a representative high-temperature reactor component is carried through to a sized surveillance article, verified with detailed finite-element analysis, and concluded with a remaining-life assessment based on an assumed ex-situ creep-rate measurement on the retrieved specimen. The report also summarizes ongoing ANL engagement with ASTM Committee E10 on Nuclear Technology and Applications toward drafting a standard covering surveillance procedures for advanced reactors.

Barua, Bipul (ORCID:0000000247184113)↗
DOE OSTI2025

Coupling between collective modes in the deformed 98 Zr nucleus: Insights from consistent HFB+QRPA calculations with the Gogny interaction

The Zirconium isotopes exhibit structural properties that present multiple challenges to nuclear theory. Investigations of the coupling present within isoscalar modes and within isovector modes are scarce but important for advancing our understanding of the microscopic picture of nuclei. To explore some of these underlying coupling features, and to test the predictive power of a state-of the-art nuclear structure approach, we provide a detailed analysis of the properties of 90,96,98Zr. This region includes a benchmarking case and offers insights into nuclear deformation phenomena. To investigate the coupling between collective modes in deformed nuclei, we focused our analysis on the ground and excited-state properties of these isotopes, employing a consistent approach with the axially-symmetric deformed Hartree-Fock-Bogoliubov (HFB) and the Quasiparticle Random Phase Approximation (QRPA) framework, both using the Gogny D1M force. This approach effectively describes both low-lying and giant-resonance states. We devoted special attention to the deformed 98Zr nucleus, where we confirm the existence of coupling between monopole and quadrupole excitations through the K π = 0 + QRPA components and demonstrate an analogous dipole-octupole coupling through the K π = 0 − and K π = 1 − components. Intrinsic transition densities and associ ated radial projections illustrate the coupling. Our work complements and extends earlier studies carried out using density-functional-based methods and notably, we included the complete Coulomb interaction also in the pairing fields, i.e. we treat terms exactly that are approximated in typical calculations that use the Gogny D1 and D2 interaction families.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Robust ab initio predictions for dimensionless ratios of 𝐸⁢2 and radius observables. II. Estimation of 𝐸⁢2 transition strengths by calibration to the charge radius

Converged results for 𝐸⁢2 observables are notoriously challenging to obtain in ab initio no-core configuration interaction approaches. Matrix elements of the 𝐸⁢2 operator are sensitive to the large-distance tails of the nuclear wave function, which converge slowly in an oscillator basis expansion. Similar convergence challenges beset ab initio prediction of the nuclear charge radius. However, we exploit systematic correlations between the calculated 𝐸⁢2 and radius observables to yield meaningful predictions for relations among these observables. In particular, we examine ab initio predictions for dimensionless ratios of the form 𝐵⁡(𝐸⁢2)/(𝑒 2 ⁢𝑟 4 ) for nuclei throughout the 𝑝 shell. Finally, meaningful predictions for 𝐸⁢2 transition strengths may then be made by calibrating to the ground-state charge radius if experimentally known.

ab initio calculations↗
DOE OSTI2024

Magnetic properties of the double trillium lattice antiferromagnet KBaCr 2 ⁢(PO 4 ) 3

Trillium lattices formed by corner-shared triangular units are the platform for magnetic frustration in three dimensions. Herein, we report structural and magnetic properties of the Cr-based (S = 3/2) double trillium lattice material KBaCr 2 ⁢(PO 4 ) 3 studied by x-ray diffraction, dc and ac magnetization, heat capacity, thermal conductivity, and 31 P nuclear magnetic resonance (NMR) measurements complemented by density-functional band-structure calculations. Magnetization, heat capacity, and 31 P NMR measurements reveal the magnetic transition at T N⁢1 ≃ 12.5 K in zero field followed by another transition (T N⁢2 ) at low temperatures in weak applied fields. The 31 P spin-lattice relaxation rate in the ordered state follows the T 3 behavior indicative of the two-magnon Raman process. Finally, the spin lattice of KBaCr 2 ⁢(PO 4 ) 3 comprises two crystallographically nonequivalent ferromagnetic sublattices that are coupled antiferromagnetically, thus eliminating frustration in this trillium network and leading to weak ferromagnetism in the vicinity of T N⁢1 .

36 MATERIALS SCIENCE↗
DOE OSTI2024

Magnetic structure of A ≤ 10 nuclei using the Norfolk nuclear models with quantum Monte Carlo methods

Here we present quantum Monte Carlo calculations of magnetic moments, form factors, and densities of A ≤ 10 nuclei within a chiral effective field theory approach. We use the Norfolk two- and three-body chiral potentials and their consistent electromagnetic one- and two-nucleon current operators. We find that two-body contributions to the magnetic moment can be large (up to ≈ 33% in A = 9 systems). We study the model dependence of these observables and place particular emphasis on investigating their sensitivity to using different cutoffs to regulate the many-nucleon operators. Calculations of elastic magnetic form factors for A ≤ 10 nuclei show excellent agreement with the data out to momentum transfers q ≈ 3 fm -1 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Confinement of quasi-atomic structures in Ti 2 N and Ti 3 N 2 MXene electrides

Metal carbides, nitrides, or carbonitrides of early transition metals, better known as MXenes, possess notable structural, electrical, and magnetic properties. Analyzing electronic structures by calculating structural stability, band structure, density of states, Bader charge transfer, and work functions utilizing first principle calculations, we revealed that titanium nitride MXenes, namely Ti 2 N and Ti 3 N 2 , have excess anionic electrons in their lattice voids, making them MXene electrides. Bulk Ti 3 N 2 has competing antiferromagnetic (AFM) and ferromagnetic (FM) configurations with slightly more stable AFM configuration, while the Ti 2 N MXene is nonmagnetic. Although Ti 3 N 2 favors AFM configuration with hexagonal crystal systems having 6/ mmm point group symmetry, Ti 3 N 2 does not support altermagnetism. The monolayer of the Ti 3 N 2 MXene is a ferromagnetic electride. These unique properties of having non-nuclear interstitial anionic electrons in the electronic structure of titanium nitride MXene have not yet been reported in the literature. Density functional theory calculations show TiN is neither an electride, MXene, or magnetic.

Anionic electrons↗
DOE OSTI2025

Building bridges between nuclear structure and reactions

The past couple of decades have seen tremendous advances in nuclear structure and reaction theory. Innovative theory frameworks for describing the nuclear many-body system, increasingly powerful computers, and opportunities for confronting theory predictions with data on unstable nuclei, have been driving the field. An important goal is to move from phenomenological ingredients in reaction calculations to predictive theories based on microscopic frameworks. We discuss ongoing efforts aimed at integrating microscopic descriptions of nuclear structure into reaction predictions for medium-mass and heavy nuclei. This contribution highlights areas where Eric Bauge, a champion for building bridges, has made important contributions by encouraging and enabling collaborations between communities with complementary expertise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Single nuclear spin detection and control in a van der Waals material

Optically active spin defects in solids are leading candidates for quantum sensing and quantum networking. Recently, single spin defects were discovered in hexagonal boron nitride (hBN), a layered van der Waals (vdW) material. Owing to its two-dimensional structure, hBN allows spin defects to be positioned closer to target samples than in three-dimensional crystals, making it ideal for atomic-scale quantum sensing, including nuclear magnetic resonance (NMR) of single molecules. However, the chemical structures of these defects remain unknown and detecting a single nuclear spin with a hBN spin defect has been elusive. Here we report the creation of single spin defects in hBN using 13 C ion implantation and the identification of three distinct defect types based on hyperfine interactions. We observed both S = 1/2 and S = 1 spin states within a single hBN spin defect. We demonstrated atomic-scale NMR and coherent control of individual nuclear spins in a vdW material, with a π-gate fidelity up to 99.75% at room temperature. By comparing experimental results with density functional theory (DFT) calculations, we propose chemical structures for these spin defects. Our work advances the understanding of single spin defects in hBN and provides a pathway to enhance quantum sensing using hBN spin defects with nuclear spins as quantum memories.

Quantum metrology↗
NASA NTRS1981

Neutron star models

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.

Canuto, V.↗
DOE OSTI2025

Short-range spin freezing in the double trillium lattice spin-liquid candidate KSrFe 2 ⁢(PO 4 ) 3 revealed via 31 P NMR

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.

Exotic phases of matter↗
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↗