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

Connecting relativistic density functional theory to microscopic calculations

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

A STUDY OF SHORT-RANGE CORRELATED PAIR FORMATION MECHANISMS

Short-Range Correlations (SRCs) refers to pairs of nucleons that are temporary high density fluctuations with high relative momenta and lower center-of-mass momenta com pared to the nuclear Fermi momentum (kF). SRCs account for 20–25% of the nucleons in medium to heavy nuclei, make up essentially all nucleons with momentum greater than kF, and contribute most of the kinetic energy carried by nucleons in nuclei. The existing semi-inclusive and exclusive measurements only cover a handful of light nu clei or heavy elements. This does not allow for a systematic study of the dependence of SRC pairs on nuclear mass and proton-neutron asymmetry. It also does not allow for insights into SRC pairing mechanisms. Therefore, we systematically studied the individual probabilities for finding SRC protons in symmetric and neutron-rich asymmetric nuclei d, 9Be, 10B, 11B, 12C, 40Ca, 48Ca, 54Fe, and 197Au. We measured the (e,e'p) reaction in kinematics dominated by scattering off mean-field nucleons (k = kF) and nucleons in SRC pairs (k = kF) at the Thomas Jefferson National Accelerator Facility (JLab) in Hall C of the Continuous Electron Beam Accelerator Facility (CEBAF) in the Fall of 2022. The measured results were used to determine the SRC pairing probabilities for protons to examine how pairing depends on nuclear mass, proton-neutron asymmetry, and nuclear shell structure. The extracted cross-section ratios were also compared to theoretical calculations. We found that SRC pair formation depends more on the nuclear shell structure with sharp increases locally within the general trend of a slower increase with larger A. We also found that intra-shell pairing has a much larger influence than inter-shell pairing. Comparisons to theory suggest that angular momentum selection rules are important to SRC pair formation and can provide new constraints for new theoretical models.

Swan, Noah [Old Dominion Univ., Norfolk, VA (Unite↗
DOE OSTI2025

Constrained nuclear–electronic orbital method for periodic density functional theory: Application to H 2 chemisorption on Si(001) surfaces

The nuclear–electronic orbital (NEO) method provides a powerful computational framework for incorporating nuclear quantum effects (NQE) in electronic structure calculations beyond the Born–Oppenheimer approximation. By incorporating additional constraints to the position operator on quantum particles like protons, the NEO method enables calculation of effective potential that accounts for NQE. Here, in this work, we present a new constrained NEO (cNEO) formulation for density functional theory (cNEO-DFT) calculations in the context of extended periodic systems. Using the nudged elastic band method, we discuss an application of the cNEO-DFT approach to studying the adsorption of a hydrogen molecule on the Si(001) surfaces. The calculation shows how NQE impacts the reaction energetics. The proton density changes are computed along the reaction pathways. This work demonstrates the capability of the new cNEO-DFT method to study a wide range of chemical processes, such as surface reactions where the quantum nature of light atoms like protons is non-negligible.

Chemical processes↗
DOE OSTI2025

Radioisotope production at the Spallation Neutron Source: Design concept of isotope production target

Upon completion of the Second Target Station (STS) Project in the mid 2030s, the Spallation Neutron Source accelerator at Oak Ridge National Laboratory will deliver a 2.7 MW proton beam to the neutron production targets. In the post-STS phase, the accelerator will have a reserve beam power capacity of at least 100 kW beyond what the two neutron production targets will receive, which could potentially be ramped up to 300 kW. In this paper, a design concept for a radioisotope production target that could utilize 250 kW of the reserve beam power capacity is presented. The target consists of thorium discs encapsulated in 316L austenitic steel shells that are cooled by water. The estimated post-irradiation activity of Ac-225 and Ra-225, critical medical radioisotopes used in targeted alpha therapy cancer treatment, is calculated at the end of bombardment after a 14 day long irradiation time. Thermal and structural analyses are performed on the basis of calculated nuclear heating data. The technical feasibility of a high-power target under a 250-kW beam load with an extremely low duty factor of $3.5\cdot 10^{-6}$ is presented from thermal, structural and fatigue lifetime perspectives.

Lee, Yong Joong [ORNL] (ORCID:0000000298381723)↗
DOE OSTI2025

Electron correlation and relativistic effects in the excited states of radium monofluoride

Highly accurate and precise electronic structure calculations of heavy radioactive atoms and their molecules are important for several research areas, including chemical, nuclear, and particle physics. Ab initio quantum chemistry can elucidate structural details in these systems that emerge from the interplay of relativistic and electron correlation effects, but the large number of electrons complicates the calculations, and the scarcity of experiments prevents insightful theory-experiment comparisons. Here we report the spectroscopy of the 14 lowest excited electronic states in the radioactive molecule radium monofluoride (RaF), which is proposed as a sensitive probe for searches of new physics. The observed excitation energies are compared with state-of-the-art relativistic Fock-space coupled cluster calculations, which achieve an agreement of ≥99.64% (within ~12 meV) with experiment for all states. Guided by theory, a firm assignment of the angular momentum and term symbol is made for 10 states and a tentative assignment for 4 states. The role of high-order electron correlation and quantum electrodynamics effects in the excitation energies is studied and found to be important for all states.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Neutrinoless Double Beta Decay of 136 Xe and Related Nuclear Structure Studies (Final Scientific Report)

The search for neutrinoless double-beta decay (0nbb) explores new physics by directly probing the unknown mass scale and possible Majorana nature of the neutrino. The nEXO experiment, one of the two leading proposed ton-scale projects in the US, has a projected sensitivity to the 136 Xe 0nbb half-life of 10 28 years. The interpretation of possible signals in this next generation of experiments would be complicated by significant variations in theoretical calculations of the nuclear matrix elements (NMEs) of the decay, and nuclear structure measurements testing those theories can help address that uncertainty. We have 136 Xe(p,n) 136 Cs reaction at TUNL to deduce the level-scheme of states in 136 Cs through which the lowest-lying 1+ state decays. We find that over 99% of such decays will proceed through at least one isomeric state with a lifetime of order 100ns, which would enable large xenon detectors to employ a delayed-coincidence technique to search for charge-exchange processes including solar neutrino interactions. We have carried out measurements at TUNL of the 134 Xe( 3 He,n) 136 Ba and 136 Xe( 3 He,n) 138 Ba reactions to low-lying 0 + states in the residuals to probe the BCS assumption for QRPA NME calculations for the initial and final nuclei in 136 Xe 0nbb. While the analyses of these reactions is not yet complete and is ongoing, our initial results indicate tension with the BCS assumption. We have assembled a thermosyphon cooling R&D system for nEXO which has demonstrated up to 500 W of cooling in vacuum. We have developed a simple thermosyphon cooling simulation for nEXO and planned laboratory tests with the R&D system to benchmark it.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Improved structure of calcium isotopes from ab initio calculations

The in-medium similarity renormalization group (IMSRG) is a powerful and flexible many-body method to compute the structure of nuclei starting from nuclear forces. Recent developments have extended the IMSRG from its standard truncation at the normal-ordered two-body level, the IMSRG(2), to a precision approximation including normal-ordered three-body operators, the IMSRG(3)-N 7 . This improvement provides a more precise solution to the many-body problem and makes it possible to quantify many-body uncertainties in IMSRG calculations. We explore the structure of 44,48,52 Ca using the IMSRG(3)-N 7 , focusing on understanding existing discrepancies of the IMSRG(2) to experimental results. We find a significantly better description of the first 2 + excitation energy of 48 Ca, improving the description of the shell closure at N=28. At the same time, we find that the IMSRG(3)-N 7 corrections to charge radii do not resolve the systematic underprediction of the puzzling large charge radius difference between 52 Ca and 48 Ca. We present estimates of many-body uncertainties of IMSRG(2) calculations applicable also to other systems based on the size extensivity of the method.

39 ≤ A ≤ 58↗
DOE OSTI2025

Robust ab initio predictions for dimensionless ratios of 𝐸⁢2 and radius observables. I. Electric quadrupole moments and deformation

We report that 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 for nuclei throughout the 𝑝 shell. Meaningful predictions for electric quadrupole moments may then be made by calibrating to the ground-state charge radius, if experimentally known, or vice versa. Moreover, these dimensionless ratios provide ab initio insight into the nuclear quadrupole deformation.

ab initio calculations↗
DOE OSTI2025

Photonuclear tomography in ultraperipheral heavy-ion collisions

We present a theoretical investigation of photonuclear tomography as a novel technique for probing the internal structure of nuclei. In this approach, ultraperipheral heavy-ion collisions (UPCs) serve as a source of intense fluxes of virtual photons, which induce coherent production of vector mesons. By analyzing the probabilities and cross sections of these photon-induced processes, we propose a methodology for reconstructing the spatial distribution of nucleons within the nucleus. Our framework provides a systematic way to access information on the nuclear geometry probed in UPCs, offering new opportunities for studies of nuclear structure using particle production as a probe. Numerical calculations for selected examples illustrate the feasibility and potential of this method.

particle & resonance production↗
DOE OSTI2025

Vector meson production in ultraperipheral heavy ion collisions

We review model calculations of exclusive vector meson production in ultraperipheral heavy ion collisions. We highlight differences and similarities between different dipole models and leading twist shadowing calculations. Recent color glass condensate calculations are presented with focus on effects from nuclear structure and azimuthal anisotropies driven by interference effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Comprehensive review of 2 β decay half-lives

Here, the double-beta (2 β )-decay is the rarest nuclear physics process, and its experimental half-lives (T 1/2 ) exceed the age of the Universe from nine to fourteen orders of magnitude. Double-beta decay was observed, and its half-life was measured in 14 parent nuclei using direct, radiochemical, and geochemical methods. The decay observables are analyzed using the Evaluated Nuclear Structure Data File (ENSDF) procedures, and the recommended T 1/2 were deduced. Using the calculated values of phase factors, the effective nuclear matrix elements were extracted and compared with available data. Thousands of theoretical and experimental works have been dedicated to these topics in the last 85 years, and we present two data sets of recommended values to encapsulate the results.

2β-decay↗
DOE OSTI2024

Comprehensive Review of 2$β$ Decay Half-Lives

The double-beta (2β)-decay is the rarest nuclear physics process, and its experimental half-lives (T 1/2 ) exceed the age of the Universe from nine to fourteen orders of magnitude. Double-beta decay was observed, and its half-life was measured in 14 parent nuclei using direct, radiochemical, and geochemical methods. The decay observables are analyzed using the Evaluated Nuclear Structure Data File (ENSDF) procedures, and the recommended T 1/2 were deduced. Using the calculated values of phase factors, the effective nuclear matrix elements were extracted and compared with available data. Thousands of theoretical and experimental works have been dedicated to these topics in the last 85 years, and we present two data sets of recommended values to encapsulate the results.

2β-decay↗
DOE OSTI2025

Determination of nuclear charge radius by extreme-ultraviolet spectroscopy of Na-like ions

We report on a method for determining the absolute nuclear charge radius of high-𝑍 elements using extreme-ultraviolet spectroscopy of highly charged Na-like ions in tandem with highly accurate atomic structure calculations of transition energy differences. The application of this method has reduced the nuclear charge radius uncertainty of 191 Ir by a factor of 8 from the currently accepted literature value, with a recently reported charge radius of 5.435(12) fm. The result reduces the charge radius uncertainty along the full Ir isotopic chain when combined with prior optical isotope shift measurements. The technique utilizes only a few million ions stored in an ion trap, which should apply to measurements with small quantities of radioactive nuclei.

Atom & ion trapping & guiding↗
DOE OSTI2026

Ab initio many-fermion structure calculations on a quantum computer

To overcome the limitations of existing algorithms for solving self-bound quantum many-body problems—such as those encountered in nuclear and particle physics—that access only a restricted subset of energy levels and provide limited structural information, we introduce and demonstrate a novel quantum-classical approach capable of resolving the complete bound-state spectrum. This method also provides the total angular momentum 𝐽 associated with each eigenstate. Here, our approach is based on expressing the Hamiltonian in second-quantized form within a novel input model combined with a scan scheme, enabling broad applicability to configuration-interaction calculations across diverse fields. We apply this hybrid method to compute, for the first time, the bound-state spectrum together with corresponding 𝐽 values of 20 O using a realistic strong-interaction Hamiltonian. Our approach applies to hadron spectra and 𝐽 values solved in the relativistic basis light-front quantization approach.

Du, Weijie [Chinese Academy of Sciences (CAS), Lan↗
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↗