BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

G. Hagen

G. Hagen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
nucl-thnucl-exphysics.comp-phcond-mat.quant-gashep-excond-mat.mes-hallhep-phphysics.atom-phquant-ph

Catalog footprint

What is connected

51works
9topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2023arXiv

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. Our results for the mirror-difference charge radii and binding energies per nucleon agree with existing data.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Systematics of E2 strength in the sd-shell with the valence-space in-medium similarity renormalization group

Background: Recent developments in {\it ab initio} nuclear theory demonstrate promising results in medium- to heavy-mass nuclei. A particular challenge for many of the many-body methodologies, however, is an accurate treatment of the electric-quadrupole, $E2$, strength associated with collectivity. Purpose: The valence-space in-medium similarity renormalization group (VS-IMSRG) is a particularly powerful method for accessing medium- and high-mass nuclei but has been found to underpredict $E2$ strengths. The purpose of this work is to evaluate the isospin dependence of this underprediction. Methods: We perform a systematic comparison of valence-space in-medium similarity renormalization group (VS-IMSRG) calculations with available literature. We make use of isoscalar and isovector contributions to the $E2$ matrix elements to assess isoscalar and isovector contributions to the missing strength. Results: It is found that the $E2$ strength is consistent throughout $T_z=\left|\frac{1}{2}\right|$, $T_z=\left|1\right|$, $T_z=\left|\frac{3}{2}\right|$ and $T_z=2$ pairs within the $sd$-shell. Furthermore, no isovector contribution to the deficiency is identified. Conclusions: A comparison with toy-models and coupled-cluster calculations is used to discuss potential origins of the missing strength, which arises from missing many-particle, many-hole excitations out of the model space. The absence of any significant isovector contribution to the missing $E2$ strength indicates that the $E2$ strength discrepancy, and therefore any correction, is largely independent of the isospin of the nuclei in question.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nuclear charge radii of Na isotopes: A tale of two theories

The accuracy of atomic theory calculations limits the extraction of nuclear charge radii from isotope shift measurements of odd-proton nuclei. For Na isotopes, though precise spectroscopic measurements have existed since more than half a century, calculations by different methods offer a wide range of values. Here, we present accurate atomic calculations to reliably extract the Na charge radii. By combining experimental matter radii with nuclear coupled-cluster calculations based on nucleon-nucleon and three-nucleon forces, we constrain the parameters obtained from the atomic calculations. Therefore, this study guides atomic theory and highlights the importance of using accurate atomic and nuclear computations in our understanding of the size of light nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Nuclear Charge Radii of the Nickel Isotopes $^{58-68,70}$Ni

Collinear laser spectroscopy is performed on the nickel isotopes $^{58-68,70}$Ni, using a time-resolved photon counting system. From the measured isotope shifts, nuclear charge radii $R_c$ are extracted and compared to theoretical results. Three ab initio approaches all employ, among others, the chiral interaction NNLO$_{\rm sat}$, which allows an assessment of their accuracy. We find agreement with experiment in differential radii $δ\left\langle r_\mathrm{c}^2 \right\rangle$ for all employed ab initio methods and interactions, while the absolute radii are consistent with data only for NNLO$_{\rm sat}$. Within nuclear density functional theory, the Skyrme functional SV-min matches experiment more closely than the Fayans functional Fy($Δr$,HFB).

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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={^1S_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 of about 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 $^1S_0$ partial wave. We demonstrate that a leading-order potential with TPE can be systematically improved by adding a contact quadratic in momenta. For momentum cutoffs below about 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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A two-neutron halo is unveiled in $^{29}$F

We report the measurement of reaction cross sections ($σ_R^{\rm ex}$) of $^{27,29}$F with a carbon target at RIKEN. The unexpectedly large $σ_R^{\rm ex}$ and derived matter radius identify $^{29}$F as the heaviest two-neutron Borromean halo to date. The halo is attributed to neutrons occupying the $2p_{3/2}$ orbital, thereby vanishing the shell closure associated with the neutron number $N = 20$. The results are explained by state-of-the-art shell model calculations. Coupled-cluster computations based on effective field theories of the strong nuclear force describe the matter radius of $^{27}$F but are challenged for $^{29}$F.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Charge radius of the short-lived $^{68}$Ni and correlation with the dipole polarizability

We present the first laser spectroscopic measurement of the neutron-rich nucleus $^{68}$Ni at the \mbox{$N=40$} subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α_{\rm D}$ has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus $^{48}$Ca. Three-particle--three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of $^{68}$Ni.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Coulomb sum rule for $^4$He and $^{16}$O from coupled-cluster theory

We demonstrate the capability of coupled-cluster theory to compute the Coulomb sum rule for the $^4$He and $^{16}$O nuclei using interactions from chiral effective field theory. We perform several checks, including a few-body benchmark for $^4$He. We provide an analysis of the center-of-mass contaminations, which we are able to safely remove. We then compare with other theoretical results and experimental data available in the literature, obtaining a fair agreement. This is a first and necessary step towards initiating a program for computing neutrino-nucleus interactions from first principles and supporting the experimental long-baseline neutrino program with a state-of-the-art theory that can reach medium-mass nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Coherent elastic neutrino-nucleus scattering on 40Ar from first principles

Coherent elastic neutrino scattering on the 40Ar nucleus is computed with coupled-cluster theory based on nuclear Hamiltonians inspired by effective field theories of quantum chromodynamics. Our approach is validated by calculating the charge form factor and comparing it to data from electron scattering. We make predictions for the weak form factor, the neutron radius, and the neutron skin, and estimate systematic uncertainties. The neutron-skin thickness of 40Ar40 is consistent with results from density functional theory. Precision measurements from coherent elastic neutrino-nucleus scattering could potentially be used to extract these observables and help to constrain nuclear models.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effective field theory in the harmonic oscillator basis

We develop interactions from chiral effective field theory (EFT) that are tailored to the harmonic oscillator basis. As a consequence, ultraviolet convergence with respect to the model space is implemented by construction and infrared convergence can be achieved by enlarging the model space for the kinetic energy. In oscillator EFT, matrix elements of EFTs formulated for continuous momenta are evaluated at the discrete momenta that stem from the diagonalization of the kinetic energy in the finite oscillator space. By fitting to realistic phase shifts and deuteron data we construct an effective interaction from chiral EFT at next-to-leading order. Many-body coupled-cluster calculations of nuclei up to 132Sn exhibit a fast convergence of ground-state energies and radii in feasible model spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Electric dipole polarizability from first principles calculations

The electric dipole polarizability quantifies the low-energy behaviour of the dipole strength and is related to critical observables such as the radii of the proton and neutron distributions. Its computation is challenging because most of the dipole strength lies in the scattering continuum. In this paper we combine integral transforms with the coupled-cluster method and compute the dipole polarizability using bound-state techniques. Employing different interactions from chiral effective field theory, we confirm the strong correlation between the dipole polarizability and the charge radius, and study its dependence on three-nucleon forces. We find good agreement with data for the 4He, 40Ca, and 16O nuclei, and predict the dipole polarizability for the rare nucleus 22O.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Emergent properties of nuclei from ab initio coupled-cluster calculations

Emergent properties such as nuclear saturation and deformation, and the effects on shell structure due to the proximity of the scattering continuum and particle decay channels are fascinating phenomena in atomic nuclei. In recent years, ab initio approaches to nuclei have taken the first steps towards tackling the computational challenge of describing these phenomena from Hamiltonians with microscopic degrees of freedom. This endeavor is now possible due to ideas from effective field theories, novel optimization strategies for nuclear interactions, ab initio methods exhibiting a soft scaling with mass number, and ever-increasing computational power. This paper reviews some of the recent accomplishments. We also present new results. The recently optimized chiral interaction NNLO$_{\rm sat}$ is shown to provide an accurate description of both charge radii and binding energies in selected light- and medium-mass nuclei up to $^{56}$Ni. We derive an efficient scheme for including continuum effects in coupled-cluster computations of nuclei based on chiral nucleon-nucleon and three-nucleon forces, and present new results for unbound states in the neutron-rich isotopes of oxygen and calcium. The coupling to the continuum impacts the energies of the $J^π= {1/2}^-,{3/2}^-,{7/2}^-,{3/2}^+$ states in $^{17,23,25}$O, and - contrary to naive shell-model expectations - the level ordering of the $J^π= {3/2}^+,{5/2}^+,{9/2}^+$ states in $^{53,55,61}$Ca.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Open $sd$-shell nuclei from first principles

We extend the ab initio coupled-cluster effective interaction (CCEI) method to deformed open-shell nuclei with protons and neutrons in the valence space, and compute binding energies and excited states of isotopes of neon and magnesium. We employ a nucleon-nucleon and three-nucleon interaction from chiral effective field theory evolved to a lower cutoff via a similarity renormalization group transformation. We find good agreement with experiment for binding energies and spectra, while charge radii of neon isotopes are underestimated. For the deformed nuclei $^{20}$Ne and $^{24}$Mg we reproduce rotational bands and electric quadrupole transitions within uncertainties estimated from an effective field theory for deformed nuclei, thereby demonstrating that collective phenomena in $sd$-shell nuclei emerge from complex ab initio calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Proton distribution radii of $^{12-19}$C illuminate features of neutron halos

Proton radii of $^{12-19}$C densities derived from first accurate charge changing cross section measurements at 900$A$ MeV with a carbon target are reported. A thick neutron surface evolves from $\sim$ 0.5 fm in $^{15}$C to $\sim$ 1 fm in $^{19}$C. The halo radius in $^{19}$C is found to be 6.4$\pm$0.7 fm as large as $^{11}$Li. Ab initio calculations based on chiral nucleon-nucleon and three-nucleon forces reproduce well the radii.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Structure of $^{78}$Ni from first principles computations

Doubly magic nuclei have a simple structure and are the cornerstones for entire regions of the nuclear chart. Theoretical insights into the supposedly doubly magic $^{78}$Ni and its neighbors are challenging because of the extreme neutron-to-proton ratio and the proximity of the continuum. We predict the $J^π=2_1^+$ state in $^{78}$Ni from a correlation with the $J^π=2_1^+$ state in $^{48}$Ca using chiral nucleon-nucleon and three-nucleon interactions. Our results confirm that $^{78}$Ni is doubly magic, and the predicted low-lying states of $^{79,80}$Ni open the way for shell-model studies of many more rare isotopes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unexpectedly large charge radii of neutron-rich calcium isotopes

Despite being a complex many-body system, the atomic nucleus exhibits simple structures for certain "magic" numbers of protons and neutrons. The calcium chain in particular is both unique and puzzling: evidence of doubly-magic features are known in 40,48Ca, and recently suggested in two radioactive isotopes, 52,54Ca. Although many properties of experimentally known Ca isotopes have been successfully described by nuclear theory, it is still a challenge to predict their charge radii evolution. Here we present the first measurements of the charge radii of 49,51,52Ca, obtained from laser spectroscopy experiments at ISOLDE, CERN. The experimental results are complemented by state-of-the-art theoretical calculations. The large and unexpected increase of the size of the neutron-rich calcium isotopes beyond N = 28 challenges the doubly-magic nature of 52Ca and opens new intriguing questions on the evolution of nuclear sizes away from stability, which are of importance for our understanding of neutron-rich atomic nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate nuclear radii and binding energies from a chiral interaction

With the goal of developing predictive ab-initio capability for light and medium-mass nuclei, two-nucleon and three-nucleon forces from chiral effective field theory are optimized simultaneously to low-energy nucleon-nucleon scattering data, as well as binding energies and radii of few-nucleon systems and selected isotopes of carbon and oxygen. Coupled-cluster calculations based on this interaction, named NNLO$_{\rm sat}$, yield accurate binding energies and radii of nuclei up to $^{40}$Ca, and are consistent with the empirical saturation point of symmetric nuclear matter. In addition, the low-lying collective $J^π=3^-$ states in $^{16}$O and $^{40}$Ca are described accurately, while spectra for selected $p$- and $sd$-shell nuclei are in reasonable agreement with experiment.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Charge, neutron, and weak size of the atomic nucleus

What is the size of the atomic nucleus? This deceivably simple question is difficult to answer. While the electric charge distributions in atomic nuclei were measured accurately already half a century ago, our knowledge of the distribution of neutrons is still deficient. In addition to constraining the size of atomic nuclei, the neutron distribution also impacts the number of nuclei that can exist and the size of neutron stars. We present an ab initio calculation of the neutron distribution of the neutron-rich nucleus $^{48}$Ca. We show that the neutron skin (difference between radii of neutron and proton distributions) is significantly smaller than previously thought. We also make predictions for the electric dipole polarizability and the weak form factor; both quantities are currently targeted by precision measurements. Based on ab initio results for $^{48}$Ca, we provide a constraint on the size of a neutron star.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Giant and pigmy dipole resonances in 4He, 16,22O, and 40Ca from chiral nucleon-nucleon interactions

We combine the coupled-cluster method and the Lorentz integral transform for the computation of inelastic reactions into the continuum. We show that the bound-state-like equation characterizing the Lorentz integral transform method can be reformulated based on extensions of the coupled-cluster equation-of-motion method, and we discuss strategies for viable numerical solutions. Starting from a chiral nucleon-nucleon interaction at next-to-next-to-next-to-leading order, we compute the giant dipole resonances of 4He, 16,22O and 40Ca, truncating the coupled-cluster equation-of-motion method at the two-particle-two-hole excitation level. Within this scheme, we find a low-lying E1 strength in the neutron-rich 22O nucleus, which compares fairly well with data from [Leistenschneider et al. Phys. Rev. Lett. 86, 5442 (2001)]. We also compute the electric dipole polariziability in 40Ca. Deficiencies of the employed Hamiltonian lead to overbinding, too small charge radii and a too small electric dipole polarizability in 40Ca.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Ab-initio coupled-cluster effective interactions for the shell model: Application to neutron-rich oxygen and carbon isotopes

We derive and compute effective valence-space shell-model interactions from ab-initio coupled-cluster theory and apply them to open-shell and neutron-rich oxygen and carbon isotopes. Our shell-model interactions are based on nucleon-nucleon and three-nucleon forces from chiral effective-field theory. We compute the energies of ground and low-lying states, and find good agreement with experiment. In particular our calculations are consistent with the N=14, 16 shell closures in oxygen-22 and oxygen-24, while for carbon-20 the corresponding N=14 closure is weaker. We find good agreement between our coupled-cluster effective-interaction results with those obtained from standard single-reference coupled-cluster calculations for up to eight valence neutrons.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Effects of three-nucleon forces and two-body currents on Gamow-Teller strengths

We optimize chiral interactions at next-to-next-to leading order to observables in two- and three-nucleon systems, and compute Gamow-Teller transitions in carbon-14, oxygen-22 and oxygen-24 using consistent two-body currents. We compute spectra of the daughter nuclei nitrogen-14, fluorine-22 and fluorine-24 via an isospin-breaking coupled-cluster technique, with several predictions. The two-body currents reduce the Ikeda sum rule, corresponding to a quenching factor q^2 ~ 0.84-0.92 of the axial-vector coupling. The half life of carbon-14 depends on the energy of the first excited 1+ state, the three-nucleon force, and the two-body current.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Efimov Physics around the neutron rich Calcium-60 isotope

We calculate the neutron-Calcium-60 S-wave scattering phase shifts using state of the art coupled-cluster theory combined with modern ab initio interactions derived from chiral effective theory. Effects of three-nucleon forces are included schematically as density dependent nucleon-nucleon interactions. This information is combined with halo effective field theory in order to investigate the Calcium-60-neutron-neutron system. We predict correlations between different three-body observables and the two-neutron separation energy of Calcium-62. This provides evidence of Efimov physics along the Calcium isotope chain. Experimental key observables that facilitate a test of our findings are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Excited states in the neutron-rich nucleus 25F

The structure of the nucleus 25F was investigated through in-beam γ-ray spectroscopy of the fragmentation of 26Ne and 27,28Na ion beams. Based on the particle-γ and particle-γγ coincidence data, a level scheme was constructed and compared with shell model and coupled-cluster calculations. Some of the observed states were interpreted as quasi single-particle states built on top of the closed-shell nucleus 24O, while the others were described as states arising from coupling of a single proton to the 2+ core excitation of 24O.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Infrared extrapolations for atomic nuclei

Harmonic oscillator model-space truncations introduce systematic errors to the calculation of binding energies and other observables. We identify the relevant infrared scaling variable and give values for this nucleus-dependent quantity. We consider isotopes of oxygen computed with the coupled-cluster method from chiral nucleon-nucleon interactions at next-to-next-to-leading order and show that the infrared component of the error is sufficiently understood to permit controlled extrapolations. By employing oscillator spaces with relatively large frequencies, well above the energy minimum, the ultraviolet corrections can be suppressed while infrared extrapolations over tens of MeVs are accurate for ground-state energies. However, robust uncertainty quantification for extrapolated quantities that fully accounts for systematic errors is not yet developed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

An optimized chiral nucleon-nucleon interaction at next-to-next-to-leading order

We optimize the nucleon-nucleon interaction from chiral effective field theory at next-to-next- to-leading order. The resulting new chiral force NNLOopt yields χ^2 \approx 1 per degree of freedom for laboratory energies below approximately 125 MeV. In the A = 3, 4 nucleon systems, the contributions of three-nucleon forces are smaller than for previous parametrizations of chiral interactions. We use NNLOopt to study properties of key nuclei and neutron matter, and demonstrate that many aspects of nuclear structure can be understood in terms of this nucleon-nucleon interaction, without explicitly invoking three-nucleon forces.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coupled-cluster calculations of nucleonic matter

Background: The equation of state (EoS) of nucleonic matter is central for the understanding of bulk nuclear properties, the physics of neutron star crusts, and the energy release in supernova explosions. Purpose: This work presents coupled-cluster calculations of infinite nucleonic matter using modern interactions from chiral effective field theory (EFT). It assesses the role of correlations beyond particle-particle and hole-hole ladders, and the role of three-nucleon-forces (3NFs) in nuclear matter calculations with chiral interactions. Methods: This work employs the optimized nucleon-nucleon NN potential NNLOopt at next-to-next-to leading-order, and presents coupled-cluster computations of the EoS for symmetric nuclear matter and neutron matter. The coupled-cluster method employs up to selected triples clusters and the single-particle space consists of a momentum-space lattice. We compare our results with benchmark calculations and control finite-size effects and shell oscillations via twist-averaged boundary conditions. Results: We provide several benchmarks to validate the formalism and show that our results exhibit a good convergence toward the thermodynamic limit. While neutron matter with interactions from chiral EFT is perturbative, symmetric nuclear matter requires nonperturbative approaches. Correlations beyond the standard particle-particle ladder approximation yield non-negligible contributions. The saturation point of symmetric nuclear matter is sensitive to the employed 3NFs and the employed regularization scheme. 3NFs with nonlocal cutoffs exhibit a considerably improved convergence than their local cousins. We are unable to find values for the parameters of the short-range part of the local 3NF that simultaneously yield acceptable values for the saturation point in symmetric nuclear matter and the binding energies of light nuclei.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coupled-cluster computations of atomic nuclei

In the past decade, coupled-cluster theory has seen a renaissance in nuclear physics, with computations of neutron-rich and medium-mass nuclei. The method is efficient for nuclei with product-state references, and it describes many aspects of weakly bound and unbound nuclei. This report reviews the technical and conceptual developments of this method in nuclear physics, and the results of coupled-cluster calculations for nucleonic matter, and for exotic isotopes of helium, oxygen, calcium, and some of their neighbors.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Coupled-cluster studies of infinite nuclear matter

The aim of this work is to develop the relevant formalism for performing coupled-cluster (CC) calculations in nuclear matter and neutron star matter, including thereby important correlations to infinite order in the interaction and testing modern nuclear forces based on chiral effective field theory. Our formalism includes the exact treatment of the so-called Pauli operator in a partial wave expansion of the equation of state. Nuclear and neutron matter calculations are done using a coupled particle-particle and hole-hole ladder approximation. The coupled ladder equations are derived as an approximation of CC theory, leaving out particle-hole and non-linear diagrams from the CC doubles amplitude equation. This study is a first step toward CC calculations for nuclear and neutron matter. We present results for both symmetric nuclear matter and pure neutron matter employing state-of-the-art nucleon-nucleon interactions based on chiral effective field theory. We employ also the newly optimized chiral interaction [A. Ekström et al., Phys. Rev. Lett. 110, 192502 (2013)] to study infinite nuclear matter. The ladder approximation method and corresponding results are compared with conventional Brueckner-Hartree-Fock theory.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spectroscopy of $^{26}$F to probe proton-neutron forces close to the drip line

A long-lived $J^π=4_1^+$ isomer, $T_{1/2}=2.2(1)$ms, has been discovered at 643.4(1) keV in the weakly-bound $^{26}_{9}$F nucleus. It was populated at GANIL in the fragmentation of a $^{36}$S beam. It decays by an internal transition to the $J^π=1_1^+$ ground state (82(14)%), by $β$-decay to $^{26}$Ne, or beta-delayed neutron emission to $^{25}$Ne. From the beta-decay studies of the $J^π=1_1^+$ and $J^π=4_1^+$ states, new excited states have been discovered in $^{25,26}$Ne. Gathering the measured binding energies of the $J^π=1_1^+-4_1^+$ multiplet in $^{26}_{9}$F, we find that the proton-neutron $π0d_{5/2} ν0d_{3/2}$ effective force used in shell-model calculations should be reduced to properly account for the weak binding of $^{26}_{9}$F. Microscopic coupled cluster theory calculations using interactions derived from chiral effective field theory are in very good agreement with the energy of the low-lying $1_1^+,2_1^+,4_1^+$ states in $^{26}$F. Including three-body forces and coupling to the continuum effects improve the agreement between experiment and theory as compared to the use of two-body forces only.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Universal properties of infrared oscillator basis extrapolations

Recent work has shown that a finite harmonic oscillator basis in nuclear many-body calculations effectively imposes a hard-wall boundary condition in coordinate space, motivating infrared extrapolation formulas for the energy and other observables. Here we further refine these formulas by studying two-body models and the deuteron. We accurately determine the box size as a function of the model space parameters, and compute scattering phase shifts in the harmonic oscillator basis. We show that the energy shift can be well approximated in terms of the asymptotic normalization coefficient and the bound-state momentum, discuss higher-order corrections for weakly bound systems, and illustrate this universal property using unitarily equivalent calculations of the deuteron.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ab-initio approach to effective single-particle energies in doubly closed shell nuclei

The present work discusses, from an ab-initio standpoint, the definition, the meaning, and the usefulness of effective single-particle energies (ESPEs) in doubly closed shell nuclei. We perform coupled-cluster calculations to quantify to what extent selected closed-shell nuclei in the oxygen and calcium isotopic chains can effectively be mapped onto an effective independent-particle picture. To do so, we revisit in detail the notion of ESPEs in the context of strongly correlated many-nucleon systems and illustrate the necessity to extract ESPEs through the diagonalization of the centroid {\it matrix}, as originally argued by Baranger. For the purpose of illustration, we analyse the impact of correlations on observable one-nucleon separation energies and non-observable ESPEs in selected closed-shell oxygen and calcium isotopes. We then state and illustrate the non-observability of ESPEs. Similarly to spectroscopic factors, ESPEs can indeed be modified by a redefinition of inaccessible quantities while leaving actual observables unchanged. This leads to the absolute necessity to employ consistent structure and reaction models based on the same nuclear Hamiltonian to extract the shell structure in a meaningful fashion from experimental data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Continuum effects and three-nucleon forces in neutron-rich oxygen isotopes

We employ interactions from chiral effective field theory and compute binding energies, excited states, and radii for isotopes of oxygen with the coupled-cluster method. Our calculation includes the effects of three-nucleon forces and of the particle continuum, both of which are important for the description of neutron-rich isotopes in the vicinity of the nucleus O-24. Our main results are the placement of the neutron drip-line at O-24, the assignment of spins, parities and resonance widths for several low-lying states of the drip-line nucleus, and an efficient approximation that incorporates the effects of three-body interactions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Elastic proton scattering of medium mass nuclei from coupled-cluster theory

Using coupled-cluster theory and interactions from chiral effective field theory, we compute overlap functions for transfer and scattering of low-energy protons on the target nucleus 40-Ca. Effects of three-nucleon forces are included phenomenologically as in-medium two-nucleon interactions. Using known asymptotic forms for one-nucleon overlap functions we derive a simple and intuitive way of computing scattering observables such as elastic scattering phase shifts and cross sections. As a first application and proof-of-principle, we compute phase shifts and differential interaction cross sections at energies 9.6 MeV and 12.44 MeV and compare with experimental data. Our computed diffraction minima are in good agreement with experiment, while we tend to overestimate the cross sections at large scattering angles.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Evolution of shell structure in neutron-rich calcium isotopes

We employ interactions from chiral effective field theory and compute the binding energies and low-lying excitations of calcium isotopes with the coupled-cluster method. Effects of three-nucleon forces are included phenomenologically as in-medium two-nucleon interactions, and the coupling to the particle continuum is taken into account using a Berggren basis. The computed ground-state energies and the low-lying 2+ states for the isotopes 42,48,50,52Ca are in good agreement with data, and we predict the excitation energy of the first 2+ state in 54Ca at 1.9 MeV, displaying only a weak sub-shell closure. In the odd-mass nuclei 53,55,61Ca we find that the positive parity states deviate strongly from the naive shell model.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Living on the edge of stability, the limits of the nuclear landscape

A first-principles description of nuclear systems along the drip lines presents a substantial theoretical and computational challenge. In this paper, we discuss the nuclear theory roadmap, some of the key theoretical approaches, and present selected results with a focus on long isotopic chains. An important conclusion, which consistently emerges from these theoretical analyses, is that three-nucleon forces are crucial for both global nuclear properties and detailed nuclear structure, and that many-body correlations due to the coupling to the particle continuum are essential as one approaches particle drip lines. In the quest for a comprehensive nuclear theory, high performance computing plays a key role.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Time-dependent coupled-cluster method for atomic nuclei

We study time-dependent coupled-cluster theory in the framework of nuclear physics. Based on Kvaal's bi-variational formulation of this method [S. Kvaal, arXiv:1201.5548], we explicitly demonstrate that observables that commute with the Hamiltonian are conserved under time evolution. We explore the role of the energy and of the similarity-transformed Hamiltonian under real and imaginary time evolution and relate the latter to similarity renormalization group transformations. Proof-of-principle computations of He-4 and O-16 in small model spaces, and computations of the Lipkin model illustrate the capabilities of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

UNEDF: Advanced Scientific Computing Collaboration Transforms the Low-Energy Nuclear Many-Body Problem

The demands of cutting-edge science are driving the need for larger and faster computing resources. With the rapidly growing scale of computing systems and the prospect of technologically disruptive architectures to meet these needs, scientists face the challenge of effectively using complex computational resources to advance scientific discovery. Multidisciplinary collaborating networks of researchers with diverse scientific backgrounds are needed to address these complex challenges. The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper describes UNEDF and identifies attributes that classify it as a successful computational collaboration. We illustrate significant milestones accomplished by UNEDF through integrative solutions using the most reliable theoretical approaches, most advanced algorithms, and leadership-class computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Ab initio coupled-cluster theory for open-shell nuclei

We develop a new method to describe properties of truly open-shell nuclei. This method is based on single-reference coupled-cluster theory and the equation-of-motion method with extensions to nuclei with $A\pm 2$ nucleons outside a closed shell. We perform proof-of-principle calculations for the ground states of the helium isotopes $^{3-6}$He and the first excited $2^+$ state in $^6$He. The comparison with exact results from matrix diagonalization in small model spaces demonstrates the accuracy of the coupled-cluster methods. Three-particle--one-hole excitations of $^4$He play an important role for the accurate description of $^6$He. For the open-shell nucleus $^6$He, the computational cost of the method is comparable with the coupled-cluster singles-and-doubles approximation while its accuracy is similar to coupled-cluster with singles, doubles and triples excitations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Exploring the anomaly in the interaction cross section and matter radius of 23O

New measurements of the interaction cross sections of 22,23O at 900A MeV performed at the GSI, Darmstadt are reported that address the unsolved puzzle of the large cross section previously observed for 23O. The matter radii for these oxygen isotopes extracted through a Glauber model analysis are in good agreement with the new predictions of the ab initio coupled-cluster theory reported here. They are consistent with a 22O+neutron description of 23O as well.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quenching of spectroscopic factors for proton removal in oxygen isotopes

We present microscopic coupled-cluster calculations of the spectroscopic factors for proton removal from the closed-shell oxygen isotopes $^{14,16,22,24,28}$O with the chiral nucleon-nucleon interaction at next-to-next-to-next-to-leading order. We include coupling-to-continuum degrees of freedom by using a Hartree-Fock basis built from a Woods-Saxon single-particle basis. This basis treats bound and continuum states on an equal footing. We find a significant quenching of spectroscopic factors in the neutron-rich oxygen isotopes, pointing to enhanced many-body correlations induced by strong coupling to the scattering continuum above the neutron emission thresholds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

UNEDF: Advanced Scientific Computing Transforms the Low-Energy Nuclear Many-Body Problem

The UNEDF SciDAC collaboration of nuclear theorists, applied mathematicians, and computer scientists is developing a comprehensive description of nuclei and their reactions that delivers maximum predictive power with quantified uncertainties. This paper illustrates significant milestones accomplished by UNEDF through integration of the theoretical approaches, advanced numerical algorithms, and leadership class computational resources.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio computation of circular quantum dots

We perform coupled-cluster and diffusion Monte Carlo calculations of the energies of circular quantum dots up to 20 electrons. The coupled-cluster calculations include triples corrections and a renormalized Coulomb interaction defined for a given number of low-lying oscillator shells. Using such a renormalized Coulomb interaction brings the coupled-cluster calculations with triples correlations in excellent agreement with the diffusion Monte Carlo calculations. This opens up perspectives for doing ab initio calculations for much larger systems of electrons.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio computation of the 17F proton-halo state and resonances in A = 17 nuclei

We perform coupled-cluster calculations of the energies and lifetimes of single-particle states around the doubly magic nucleus 16O based on chiral nucleon-nucleon interactions at next-to-next- to-next-to-leading order. To incorporate effects from the scattering continuum, we solve the coupled- cluster equations with a Gamow-Hartree-Fock basis. Our calculations for the 1/2+ proton-halo state in 17F and the 1/2+ state in 17O agree well with experiment, while the calculated spin-orbit splitting between d5/2 and d3/2 states is too small due to the lack of three-nucleon forces. We find that continuum effects yield a significant amount of additional binding energy for the 1/2+ and 3/2+ states in 17O and 17F.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ab initio coupled-cluster approach to nuclear structure with modern nucleon-nucleon interactions

We perform coupled-cluster calculations for the doubly magic nuclei 4He, 16O, 40Ca and 48Ca, for neutron-rich isotopes of oxygen and fluorine, and employ "bare" and secondary renormalized nucleon-nucleon interactions. For the nucleon-nucleon interaction from chiral effective field theory at order next-to-next-to-next-to leading order, we find that the coupled-cluster approximation including triples corrections binds nuclei within 0.4 MeV per nucleon compared to data. We employ interactions from a resolution-scale dependent similarity renormalization group transformations and assess the validity of power counting estimates in medium-mass nuclei. We find that the missing contributions due to three-nucleon forces are consistent with these estimates. For the unitary correlator model potential, we find a slow convergence with respect to increasing the size of the model space. For the G-matrix approach, we find a weak dependence of ground-state energies on the starting energy combined with a rather slow convergence with respect to increasing model spaces. We also analyze the center-of-mass problem and present a practical and efficient solution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Closed-shell properties of $^{24}$O with {\em ab initio} coupled-cluster theory

We present an \emph{ab initio} calculation of spectroscopic factors for neutron and proton removal from $^{24}$O using the coupled-cluster method and a state-of-the-art chiral nucleon-nucleon interaction at next-to-next-to-next-to-leading order. In order to account for the coupling to the scattering continuum we use a Berggren single-particle basis that treats bound, resonant, and continuum states on an equal footing. We report neutron removal spectroscopic factors for the $^{23}$O states $J^π = 1/2^+$, $5/2^+$, $3/2^-$ and $1/2^-$, and proton removal spectroscopic factors for the $^{23}$N states $1/2^-$ and $3/2^-$. Our calculations support the accumulated experimental evidence that $^{24}$O is a closed-shell nucleus.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Computation of spectroscopic factors with the coupled-cluster method

We present a calculation of spectroscopic factors within coupled-cluster theory. Our derivation of algebraic equations for the one-body overlap functions are based on coupled-cluster equation-of-motion solutions for the ground and excited states of the doubly magic nucleus with mass number $A$ and the odd-mass neighbor with mass $A-1$. As a proof-of-principle calculation, we consider $^{16}$O and the odd neighbors $^{15}$O and $^{15}$N, and compute the spectroscopic factor for nucleon removal from $^{16}$O. We employ a renormalized low-momentum interaction of the $V_{\mathrm{low-}k}$ type derived from a chiral interaction at next-to-next-to-next-to-leading order. We study the sensitivity of our results by variation of the momentum cutoff, and then discuss the treatment of the center of mass.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Many-body interactions and nuclear structure

This article presents several challenges to nuclear many-body theory and our understanding of the stability of nuclear matte r. In order to achieve this, we present five different cases, starting with an idealized toy model. These cases expose problems that need to be understood in order to match recent advances in nuclear theory with current experimental programs in low-energy nuclear physics. In particular, we focus on our current understanding, or lack thereof, of many-body forces, and how they evolve as functions of the number of particles . We provide examples of discrepancies between theory and experiment and outline some selected perspectives for future research directions.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Gamow shell model and realistic nucleon-nucleon interactions

We present a new and efficient method to obtain a Gamow shell-model basis and matrix elements generated by realistic nucleon-nucleon interactions. We derive a self-consistent Hartree-Fock potential from the renormalized N3LO interaction model. The corresponding Gamow one-body eigenstates are generated in a plane wave basis in order to build a Gamow shell-model set of basis states for the closed shell nuclei 4He and 16O. We address also the problem of representing a realistic nucleon-nucleon interaction in a two-particle Berggren basis in the laboratory frame. To achieve this, an expansion of matrix elements of the residual nucleon-nucleon interaction in a finite set of harmonic oscillator wave functions is used. We show that all loosely bound and narrow resonant states converge surprisingly fast. Even broad resonances in these two-particle valence systems converge within a reasonable number of harmonic oscillator functions. Examples of 6He and 18O Gamow shell-model calculations using 4He and 16}O as closed shell cores are presented. This procedure allows Gamow shell-model calculations to be performed with all realistic nucleon-nucleon interactions and with either momentum or position space representations for the Gamow basis. The possibility to remove the center of mass spuriosity of Gamow shell-model nuclear states with this method is also discussed. Perspectives for nuclear structure calculations of dripline nuclei are outlined.

0 citations0 reviewsNo code linkedOpen access