Search NASA⌕ Search

SEARCH · Search NASA

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

Mixed mode stress-intensity-factors in mode-3 loaded middle crack tension specimen

A three dimensional stress analysis of a middle-crack tension specimen subjected to mode-3 type loading was performed using fracture mechanics based finite element code FRAC3D. Three-dimensional stress intensity factors were calculated for a range of specimen thicknesses that represent the structures used in aerospace and nuclear industries. Calculated SIF for very thick specimen (thickness-to-crack length b/a greater than or equal to 30) agreed very well with the antiplane solution in the literature. The K(sub II) stress field exists near the intersection of the crack front and free surface in a boundary-layer region covers the complete thickness of the plate and K(sub II) dominates all through the thickness. For very thin plates (b/a is less than .1), the average K(sub II) is larger than K(sub III) (about 25% for b/a = 0.1).

Shivakumar, Kunigal N.↗

Uncertainties in tellurium-based dark matter searches stemming from nuclear structure uncertainties

Using tellurium dioxide as a target, we calculate uncertainties on 90% upper confidence limits of Galilean effective field theory (Galilean EFT) couplings to a weakly interacting massive particle (WIMP) dark matter candidate due to uncertainties in nuclear shell models. We find that these uncertainties in naturally occurring tellurium isotopes are comparable across the different Galilean EFT couplings to uncertainties in xenon, with some reaching over 100%. We also consider the effect these nuclear uncertainties have on estimates of the annual modulation of dark matter from these searches, finding that the uncertainties in the modulation amplitude are proportional to the nonmodulating upper confidence limit uncertainties. We also show that the determination of the modulation phase is insensitive to changes in the nuclear model for a given isotope.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

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↗

First β -delayed neutron spectroscopy of 24 O

The β decay of 24 O was investigated at NSCL/MSU using a combination of neutron and γ-ray spectroscopy. For the first time, the β-delayed neutron energy spectrum of 24 O was measured, revealing three intensely populated, isolated neutron-unbound states in 24 F. This allowed for the extraction of the decay strength in 24 F up to 6.2 MeV. A comprehensive comparison of the experimental results with various nuclear theories, ranging from the empirical shell model to the most advanced ab initio calculations, was conducted. While most theoretical predictions align with the experimental data for low-lying states, discrepancies arise at higher excitation energies. Finally, in the transition from 24 O to 24 F, shell model calculations using the empirical USDB interaction predicted the structure of both nuclei without invoking the need for a stronger proton-neutron tensor force, which was postulated for the neighboring isotone 25 F.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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↗

Sierra/SD – Verification Test Manual – 5.22

Verification and validation (V&V) of scientific computing programs are important at Sandia National Labs due to the expanding role of computational simulation in managing the United States nuclear stockpile. The complexities of structural response calculations used to analyze physical problems, the varieties of codes applied to the calculations, and the importance of accurate predictions when assessing field conditions demand confidence in the consistency and accuracy of computer codes. Confidence in the accuracy of the predictions arising from computer simulations must ultimately be gained through verification and validation. The Sierra salinas structural dynamics analysis code, Sierra/SD, is used at the DOE Laboratories, and in several DOD projects. The roles of Sierra/SD in the qualification of weapon systems and components for normal and hostile environments throughout the Stockpile-to-Target Sequence include to, • Redesign weapon components. • Certify weapon components and systems for target environments such as hypersonic vehicles. • Certify that components will survive the thermal mechanical shock loads associated with hostile environments. • Evaluate current stockpile issues, including issues associated with uncertainty quantification. • Address many other problems that are encountered in stockpile management. The Sierra/SD verification plan is described, and an evolving set of key verification tests are described in detail. The verification tests ensure the correctness of the mathematics and numerical algorithms associated with functionality describing engineering phenomena. Development is in accordance with a set of tailored Software Quality Engineering (SQE) practices. SQE practices guide the overall verification and validation effort.

97 MATHEMATICS AND COMPUTING↗

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)↗

Absolute nuclear charge radius by Na-like spectral line separation in high-Z elements

Abstract We describe a novel technique to determine absolute nuclear radii of high- Z nuclides. Utilizing accurate theoretical atomic structure calculations together with precise measurements of extreme ultraviolet transitions in highly charged ions this method allows for precise determinations of absolute nuclear charge radii based upon the well-known nuclear radii of their neighboring elements. This method can work for elements without stable isotopes, and its accuracy may be competitive with current methods (electron scattering and muonic x-ray spectroscopy).

Hosier, A. (ORCID:0000000328366605)↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗