Engineering topics

Hagen, Gaute

Publications and source records attributed to Hagen, Gaute.

At least 19 records

DOE OSTI2024

Electromagnetic observables of open-shell nuclei from coupled-cluster theory

Here we develop a new method to describe electromagnetic observables of open-shell nuclei with two nucleons outside a closed shell. This approach combines the equation-of-motion coupled-cluster method for such systems and the Lorentz integral transform technique, expanding the applicability of coupled-cluster theory for electromagnetic observables beyond closed-(sub)shell nuclei. To validate this new approach, we compute the non-energy-weighted dipole sum rule and the dipole polarizability of 16,24 ⁢O in both the closed-(sub)shell and the new equation-of-motion coupled-cluster frameworks, finding agreement within error bars. We then analyze the evolution of the dipole polarizability along the oxygen and calcium isotopic chains. Our predictions agree well with available experimental data and other available theoretical calculations for the closed-(sub)shell 16,22 ⁢O and the open-shell 18 O. In the calcium isotopes, we observe that our dipole polarizability predictions for open-shell nuclei are lower than those of closed-(sub)shell nuclei. We expect that our predictions for 24 O and 54,56 ⁢Ca will motivate future experimental studies at the drip line.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Ab initio computations of strongly deformed nuclei near Zr 80

Nuclei around N ≈ Z ≈ 40 are strongly deformed and exhibit coexistence of shapes. These phenomena have challenged nuclear models. Here, we perform ab initio coupled-cluster computations of low-lying collective states and electromagnetic quadrupole transitions of the even-even nuclei 72 Kr, 76,78 Sr, 78,80 Zr, and 84 Mo starting from chiral nucleon-nucleon and three-nucleon forces. Our calculations reproduce the coexistence of oblate and prolate shapes in these nuclei, yield rotational bands and strong electromagnetic transitions, but are not accurate for some observables and nuclei. Finally, these results highlight the advances and challenges of ab initio computations of heavy deformed nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

40 Ca transverse response function from coupled-cluster theory

Here, we present calculations of the 40 Ca transverse response function obtained from coupled-cluster theory used in conjunction with the Lorentz integral transform method. We employ nuclear forces derived at next-to-next-to leading order in chiral effective field theory with and without Δ degrees of freedom. We first benchmark this approach on the 4 He nucleus and compare both the transverse sum rule and the response function to earlier calculations based on different methods. As expected from the power counting of the chiral expansion of electromagnetic currents and from previous studies, our results retaining only one-body term underestimate the experimental data for 4 He by about 20%. However, when the method is applied to 40 Ca at the same order of the expansion, response functions do not lack strength and agree well with the world electron scattering data. We discuss various sources of theoretical uncertainties and comment on the comparison of our results with the available experiments.

39 ≤ A ≤ 58↗

DOE OSTI2023

Surprising Charge-Radius Kink in the Sc Isotopes at N = 20

For this work, charge radii of neutron deficient 40 Sc and 41 Sc nuclei were determined using collinear laser spectroscopy. With the new data, the chain of Sc charge radii extends below the neutron magic number N=20 and shows a pronounced kink, generally taken as a signature of a shell closure, but one notably absent in the neighboring Ca, K, and Ar isotopic chains. Theoretical models that explain the trend at N=20 for the Ca isotopes cannot reproduce this puzzling behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Coupled-cluster theory for strong entanglement in nuclei

Atomic nuclei can exhibit shape coexistence and multireference physics that enters in their ground states, and accurately capturing the ensuing correlations and entanglement is challenging. Here, we address this problem by applying single-reference coupled-cluster theory based on spherical and deformed reference states and the tailored coupled-cluster method. The latter combines configuration interactions to capture static correlations with coupled-cluster theory for dynamic correlations. We compute the atomic nuclei 12 C, 28 Si, and 56 Ni and find that the tailored coupled-cluster method and the single-reference approach based on a deformed Hartree-Fock state yield the most accurate results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Level structures of 56,58 Ca cast doubt on a doubly magic 60 Ca

Gamma decays were observed in 56 Ca and 58 Ca following quasi-free one-proton knockout reactions from 57,59 Sc beams at ≈200 MeV/nucleon. For 56 Ca, a γ ray transition was measured to be 1456(12) keV, while for 58 Ca an indication for a transition was observed at 1115(34) keV. Both transitions were tentatively assigned as the $2^{+}_{1}\rightarrow0^{+}_{gs}$ decays, and were compared to results from ab initio and conventional shell-model approaches. A shell-model calculation in a wide model space with a marginally modified effective nucleon-nucleon interaction depicts excellent agreement with experiment for $2^{+}_{1}$ level energies, two-neutron separation energies, and reaction cross sections, corroborating the formation of a new nuclear shell above the N = 34 shell. Its constituents, the 0$f$ 5/2 and 0$g$ 9/2 orbitals, are almost degenerate. This degeneracy precludes the possibility for a doubly magic 60 Ca and potentially drives the dripline of Ca isotopes to 70 Ca or even beyond.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Trends of Neutron Skins and Radii of Mirror Nuclei from First Principles

The neutron skin of atomic nuclei impacts the structure of neutron-rich nuclei, the equation of state of nucleonic matter, and the size of neutron stars. Here we predict the neutron skin of selected light- and medium-mass nuclei using coupled-cluster theory and the auxiliary field diffusion Monte Carlo method with two- and three-nucleon forces from chiral effective field theory. We find a linear correlation between the neutron skin and the isospin asymmetry in agreement with the liquid-drop model and compare with data. We also extract the linear relationship that describes the difference between neutron and proton radii of mirror nuclei and quantify the effect of charge symmetry breaking terms in the nuclear Hamiltonian. In conclusion, our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Towards precise and accurate calculations of neutrinoless double-beta decay

We present the results of a National Science Foundation Project Scoping Workshop, the purpose of which was to assess the current status of calculations for the nuclear matrix elements governing neutrinoless double-beta decay and determine if more work on them is required. After reviewing important recent progress in the application of effective field theory, lattice quantum chromodynamics, and ab initio nuclear-structure theory to double-beta decay, we discuss the state of the art in nuclear-physics uncertainty quantification and then construct a roadmap for work in all these areas to fully complement the increasingly sensitive experiments in operation and under development. The roadmap includes specific projects in theoretical and computational physics as well as the use of Bayesian methods to quantify both intra- and inter-model uncertainties. Here, the goal of this ambitious program is a set of accurate and precise matrix elements, in all nuclei of interest to experimentalists, delivered together with carefully assessed uncertainties. Such calculations will allow crisp conclusions from the observation or non-observation of neutrinoless double-beta decay, no matter what new physics is at play.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

How to renormalize coupled cluster theory

Coupled cluster theory is an attractive tool to solve the quantum many-body problem because its singles and doubles (CCSD) approximation is computationally affordable and yields about 90% of the correlation energy. Capturing the remaining 10%, e.g., via including triples, is numerically expensive. In this work, we assume that short-range three-body correlations dominate and—following Lepage (arXiv:nucl-th/9706029)—that their effects can be included within CCSD by renormalizing the three-body contact interaction. We renormalize this contact in 16 O and obtain systematically improved CCSD results for 24 O, 20–34 Ne, 40,48 Ca, 78 Ni, 90 Zr, and 100 Sn.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Proton Distribution Radii of 16–24 O: Signatures of New Shell Closures and Neutron Skin

The root mean square radii of the proton density distribution in 16–24 O derived from measurements of charge changing cross sections with a carbon target at ~900A MeV together with the matter radii portray thick neutron skin for 22–24 O despite 22,24 O being doubly magic. Imprints of the shell closures at N=14 and 16 are reflected in local minima of their proton radii that provide evidence for the tensor interaction causing them. The radii agree with ab initio calculations employing the chiral NNLO sat interaction, though skin thickness predictions are challenged. Furthermore, shell model predictions agree well with the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Two-pion exchange as a leading-order contribution in chiral effective field theory

Pion exchange is the central ingredient to nucleon-nucleon interactions used in nuclear structure calculations, and one-pion exchange (OPE) enters at leading order in chiral effective field theory. In the 2S+1 L J = 1 S 0 partial wave, however, OPE and a contact term needed for proper renormalization fail to produce the qualitative, and quantitative, features of the scattering phase shifts. Cutoff variation also revealed a surprisingly low breakdown momentum Λ b ≈ 330 MeV in this partial wave. Here we show that potentials consisting of OPE, two-pion exchange (TPE), and a single contact address these problems and yield accurate and renormalization group (RG) invariant phase shifts in the 1 S 0 partial wave. Here, we demonstrate that a leading-order potential with TPE can be systematically improved by adding a contact quadratic in momenta. For momentum cutoffs Λ ≲ 500 MeV, the removal of relevant physics from TPE loops needs to be compensated by additional contacts to keep RG invariance. Inclusion of the Δ isobar degree of freedom in the potential does not change the strong contributions of TPE.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Ab initio predictions link the neutron skin of 208 Pb to nuclear forces

Heavy atomic nuclei have an excess of neutrons over protons, which leads to the formation of a neutron skin whose thickness is sensitive to details of the nuclear force. This links atomic nuclei to properties of neutron stars, thereby relating objects that differ in size by orders of magnitude. The nucleus 208 Pb is of particular interest because it exhibits a simple structure and is experimentally accessible. However, computing such a heavy nucleus has been out of reach for ab initio theory. By combining advances in quantum many-body methods, statistical tools and emulator technology, we make quantitative predictions for the properties of 208 Pb starting from nuclear forces that are consistent with symmetries of low-energy quantum chromodynamics. We explore 10 9 different nuclear force parameterizations via history matching, confront them with data in select light nuclei and arrive at an importance-weighted ensemble of interactions. We accurately reproduce bulk properties of 208 Pb and determine the neutron skin thickness, which is smaller and more precise than a recent extraction from parity-violating electron scattering but in agreement with other experimental probes. This work demonstrates how realistic two- and three-nucleon forces act in a heavy nucleus and allows us to make quantitative predictions across the nuclear landscape.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Angular-momentum projection in coupled-cluster theory: Structure of 34 Mg

Single- reference coupled-cluster theory is an accurate and affordable computational method for the nuclear many-body problem. For open-shell nuclei, the reference state typically breaks rotational invariance and angular momentum must be restored as a good quantum number. We perform angular-momentum projection after variation and employ the disentangled coupled-cluster formalism and a Hermitian approach. We compare our results with benchmarks for 8 Be and 20 Ne using a two-nucleon interaction from chiral effective field theory and for pf-shell nuclei within the traditional shell model. We compute the rotational band in the exotic nucleus 34 Mg and find agreement with data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

Ab initio coupled-cluster calculations of ground and dipole excited states in He 8

We perform coupled-cluster calculations of ground- and dipole excited-state properties of the 8 He halo nucleus with nucleon-nucleon and three-nucleon interactions from chiral effective field theory, both with and without explicit delta degrees of freedom. By increasing the precision in our coupled-cluster calculations via the inclusion of leading-order three-particle three-hole excitations in the cluster operator, we obtain a ground-state energy and a charge radius that are consistent with experiment, albeit with a slight underbinding. We also investigate the excited states induced by the electric dipole operator and present a discussion on the Thomas-Reiche-Kuhn and cluster sum rules. Lastly, we compute the electric dipole polarizability, providing a theoretical benchmark for future experimental determinations that will study this exotic nucleus.

6 ≤ A ≤ 19↗

DOE OSTI2022

Universal trend of charge radii of even-even Ca–Zn nuclei

Radii of nuclear charge distributions carry information about the strong and electromagnetic forces acting inside the atomic nucleus. Whereas the global behavior of nuclear charge radii is governed by the bulk properties of nuclear matter, their local trends are affected by quantum motion of proton and neutron nuclear constituents. The measured differential charge radii δ$\langle$$r^{2}_{c}$$\rangle$ between neutron numbers N = 28 and N = 40 exhibit a universal pattern as a function of n = N – 28 that is independent of the atomic number. In this work, we analyze this remarkable behavior in even-even nuclei from calcium to zinc using two state-of-the-art theories based on quantified nuclear interactions: the ab initio coupled cluster theory and nuclear density functional theory. Both theories reproduce the smooth rise of differential charge radii and their weak dependence on the atomic number. By considering a large set of isotopic chains, we show that this trend can be captured by just two parameters: the slope and curvature of δ$\langle$$r^{2}_{c}$$\rangle$(n). We demonstrate that these parameters show appreciable model dependence, and the statistical analysis indicates that they are not correlated with any single model property, i.e., they are impacted by both bulk nuclear properties as well as shell structure.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Effective shell-model interaction for nuclei “southeast” of 100 Sn

We construct an effective shell-model interaction for the valence space spanned by single-particle neutron and single-hole proton states in 100 Sn. Starting from chiral nucleon-nucleon and three-nucleon forces and single-reference coupled-cluster theory for 100 Sn we apply a second similarity transformation that decouples the valence space. The particle-particle components of the resulting effective interaction can be used in shell model calculations for neutron deficient tin isotopes. The hole-hole interaction can be used to calculate the N=50 isotones south of 100 Sn, and the full particle-hole interaction describes nuclei in the region southeast of 100 Sn. We compute low-lying excited states in selected nuclei southeast of 100 Sn, and find reasonable agreement with data. The presented techniques can also be applied to construct effective shell-model interactions for other regions of the nuclear chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Proton inelastic scattering reveals deformation in 8 He

A measurement of proton inelastic scattering of 8He at 8.25A MeV at TRIUMF shows a resonance at 3.54(6) MeV with a width of 0.89(11) MeV. The energy of the state is in good agreement with coupled cluster and no-core shell model with continuum calculations, with the latter successfully describing the measured resonance width as well. Its differential cross section analyzed with phenomenological collective excitation form factor and microscopic coupled reaction channels framework consistently reveals a large deformation parameter = 0.40(3), consistent with no-core shell model predictions of a large neutron deformation. This deformed double-closed shell at the neutron drip-line opens a new paradigm.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Mass measurements of 99–101 In challenge ab initio nuclear theory of the nuclide 100 Sn

The tin isotope 100 Sn is of singular interest for nuclear structure due to its closed-shell proton and neutron configurations. It is also the heaviest nucleus comprising protons and neutrons in equal numbers—a feature that enhances the contribution of the short-range proton–neutron pairing interaction and strongly influences its decay via the weak interaction. Decay studies in the region of 100 Sn have attempted to prove its doubly magic character but few have studied it from an ab initio theoretical perspective, and none of these has addressed the odd-proton neighbours, which are inherently more difficult to describe but crucial for a complete test of nuclear forces. Here we present direct mass measurements of the exotic odd-proton nuclide 100 In, the beta-decay daughter of 100 Sn, and of 99 In, with one proton less than 100 Sn. We use advanced mass spectrometry techniques to measure 99 In, which is produced at a rate of only a few ions per second, and to resolve the ground and isomeric states in 101 In. The experimental results are compared with ab initio many-body calculations. The 100-fold improvement in precision of the 100 In mass value highlights a discrepancy in the atomic-mass values of 100 Sn deduced from recent beta-decay results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗