Showing 1–50 of 50 results for author: Ekström, A

Uncertainty Quantification of Collective Nuclear Observables From the Chiral Potential Parametrization

Authors: Kevin S. Becker, Kristina D. Launey, Andreas Ekström, Tomáš Dytrych, Daniel Langr, Grigor H. Sargsyan, Jerry P. Draayer

Abstract: We perform an uncertainty estimate of quadrupole moments and B(E2) transition rates that inform nuclear collectivity. In particular, we study the low-lying states of 6Li and 12C using the ab initio symmetry-adapted no-core shell model. For a narrow standard deviation of approximately 1% on the low-energy constants which parametrize high-precision chiral potentials, we find output standard deviatio… ▽ More We perform an uncertainty estimate of quadrupole moments and B(E2) transition rates that inform nuclear collectivity. In particular, we study the low-lying states of 6Li and 12C using the ab initio symmetry-adapted no-core shell model. For a narrow standard deviation of approximately 1% on the low-energy constants which parametrize high-precision chiral potentials, we find output standard deviations in the collective observables ranging from approximately 3-6%. The results mark the first step towards a rigorous uncertainty quantification of collectivity in nuclei that aims to account for all sources of uncertainty in ab initio descriptions of challenging collective and clustering observables. △ Less

Submitted 27 March, 2024; originally announced April 2024.

Comments: 12 pages, 10 figures, proceedings article for the International Workshop "Shapes and Dynamics of Atomic Nuclei: Contemporary Aspects" (SDANCA-23)

Multiscale physics of atomic nuclei from first principles

Authors: Z. H. Sun, A. Ekström, C. Forssén, G. Hagen, G. R. Jansen, T. Papenbrock

Abstract: Atomic nuclei exhibit multiple energy scales ranging from hundreds of MeV in binding energies to fractions of an MeV for low-lying collective excitations. As the limits of nuclear binding is approached near the neutron- and proton driplines, traditional shell-structure starts to melt with an onset of deformation and an emergence of coexisting shapes. It is a long-standing challenge to describe thi… ▽ More Atomic nuclei exhibit multiple energy scales ranging from hundreds of MeV in binding energies to fractions of an MeV for low-lying collective excitations. As the limits of nuclear binding is approached near the neutron- and proton driplines, traditional shell-structure starts to melt with an onset of deformation and an emergence of coexisting shapes. It is a long-standing challenge to describe this multiscale physics starting from nuclear forces with roots in quantum chromodynamics. Here we achieve this within a unified and non-perturbative framework that captures both short- and long-range correlations starting from modern nucleon-nucleon and three-nucleon forces from chiral effective field theory. The short-range correlations which accounts for the bulk of the binding energy is included within a symmetry-breaking framework, while long-range correlations (and fine details about the collective structure) are included via symmetry projection. Our calculations accurately reproduce available experimental data for low-lying collective states and the electromagnetic quadrupole transitions in $^{20-30}$Ne. We also reveal coexisting spherical and deformed shapes in $^{30}$Ne, which indicates the breakdown of the magic neutron number $N=20$ as the key nucleus $^{28}$O is approached, and we predict that the dripline nuclei $^{32,34}$Ne are strongly deformed. By developing reduced-order-models for symmetry-projected states, we perform a global sensitivity analysis and find that the subleading singlet S-wave contact and a pion-nucleon coupling strongly impact nuclear deformation in chiral effective-field-theory. The techniques developed in this work clarify how microscopic nuclear forces generate the multiscale physics of nuclei spanning collective phenomena as well as short-range correlations and allow to capture emergent and dynamical phenomena in finite fermion systems. △ Less

Submitted 26 March, 2024; originally announced April 2024.

Comments: This work significantly expands upon ideas from arXiv:2305.06955. While it includes some overlap in text to maintain context, substantial new content, methodology, and analysis are presented. 26 pages, 16 figures

Perturbative computations of neutron-proton scattering observables using renormalization-group invariant $χ$EFT up to N$^3$LO

Authors: Oliver Thim, Andreas Ekström, Christian Forssén

Abstract: We predict neutron-proton scattering cross-sections and polarization observables up to next-to-next-to-next-to leading order in a renormalization-group invariant description of the strong nucleon-nucleon interaction. Low-energy constants are calibrated to phase shifts, sub-leading corrections are computed in distorted-wave perturbation theory, and we employ momentum-cutoff values 500 and 2500 MeV.… ▽ More We predict neutron-proton scattering cross-sections and polarization observables up to next-to-next-to-next-to leading order in a renormalization-group invariant description of the strong nucleon-nucleon interaction. Low-energy constants are calibrated to phase shifts, sub-leading corrections are computed in distorted-wave perturbation theory, and we employ momentum-cutoff values 500 and 2500 MeV. We find a steady order-by-order convergence and realistic descriptions of scattering observables up to a laboratory scattering energy of approximately 100 MeV. We also compare perturbative and non-perturbative calculations for phase shifts and cross sections and quantify how unitarity is gradually restored at higher orders. The perturbative approach offers an important diagnostic tool for any power counting and our results suggest that the breakdown scale in chiral effective field theory might be significantly lower than estimates obtained in non-perturbative calculations. △ Less

Submitted 11 June, 2024; v1 submitted 23 February, 2024; originally announced February 2024.

Comments: 18 pages, 6 figures

Journal ref: Phys. Rev. C 109, 064001 (2024)

Understanding the Effect of Chiral NN Parametrization on Nuclear Shapes From an Ab Initio Perspective

Authors: Kevin S. Becker, Kristina D. Launey, Andreas Ekström, Grigor H. Sargsyan, Darin C. Mumma, Tomáš Dytrych, Daniel Langr, Jerry P. Draayer

Abstract: The ab initio symmetry-adapted no-core shell model naturally describes nuclear deformation and collectivity, and is therefore well-suited to studying the dynamics and coexistence of shapes in atomic nuclei. For the first time, we analyze how these features in low-lying states of 6Li and 12C are impacted by the underlying realistic nucleon-nucleon interaction. We find that the interaction parametri… ▽ More The ab initio symmetry-adapted no-core shell model naturally describes nuclear deformation and collectivity, and is therefore well-suited to studying the dynamics and coexistence of shapes in atomic nuclei. For the first time, we analyze how these features in low-lying states of 6Li and 12C are impacted by the underlying realistic nucleon-nucleon interaction. We find that the interaction parametrization has a notable but limited effect on collective shapes in the lowest 6Li and 12C states, while collective structures in the excited 2+ state of 12C are significantly more sensitive to the interaction parameters and exhibits emergent shape coexistence. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Comments: 4 pages + references, 3 figures, CGS17 (2023) conference proceedings contribution

Eigenvector Continuation and Projection-Based Emulators

Authors: Thomas Duguet, Andreas Ekström, Richard J. Furnstahl, Sebastian König, Dean Lee

Abstract: Eigenvector continuation is a computational method for parametric eigenvalue problems that uses subspace projection with a basis derived from eigenvector snapshots from different parameter sets. It is part of a broader class of subspace-projection techniques called reduced-basis methods. In this colloquium article, we present the development, theory, and applications of eigenvector continuation an… ▽ More Eigenvector continuation is a computational method for parametric eigenvalue problems that uses subspace projection with a basis derived from eigenvector snapshots from different parameter sets. It is part of a broader class of subspace-projection techniques called reduced-basis methods. In this colloquium article, we present the development, theory, and applications of eigenvector continuation and projection-based emulators. We introduce the basic concepts, discuss the underlying theory and convergence properties, and present recent applications for quantum systems and future prospects. △ Less

Submitted 17 June, 2024; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: Final version to appear as colloquium article in Rev. Mod. Phys., 22 pages, 17 figures

How chiral forces shape neutron-rich Ne and Mg nuclei

Authors: Andreas Ekström, Christian Forssén, G. Hagen, G. R. Jansen, T. Papenbrock, Z. H. Sun

Abstract: We compute the structure of the exotic even nuclei $^{20-34}$Ne and $^{34-40}$Mg using interactions from chiral effective field theory (EFT). Our results for the ground-state rotational bands in $^{20-32}$Ne and $^{36-40}$Mg agree with data. We predict a well-deformed $^{34}$Ne and find that $^{40}$Mg exhibits an oblate deformed band close to the prolate ground-state, indicating the emergence of s… ▽ More We compute the structure of the exotic even nuclei $^{20-34}$Ne and $^{34-40}$Mg using interactions from chiral effective field theory (EFT). Our results for the ground-state rotational bands in $^{20-32}$Ne and $^{36-40}$Mg agree with data. We predict a well-deformed $^{34}$Ne and find that $^{40}$Mg exhibits an oblate deformed band close to the prolate ground-state, indicating the emergence of shape co-existence at the neutron dripline. A global sensitivity analysis shows that the subleading singlet $S$-wave contact and a pion-nucleon coupling strongly impact deformation in chiral EFT. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Comments: 13 pages, 8 figures, supplemental material

Inference of the low-energy constants in delta-full chiral effective field theory including a correlated truncation error

Authors: Isak Svensson, Andreas Ekström, Christian Forssén

Abstract: We sample the posterior probability distributions of the low-energy constants (LECs) in delta-full chiral effective field theory ($χ$EFT) up to third order. We use eigenvector continuation for fast and accurate emulation of the likelihood and Hamiltonian Monte Carlo to draw effectively independent samples from the posteriors. Our Bayesian inference is conditioned on the Granada database of neutron… ▽ More We sample the posterior probability distributions of the low-energy constants (LECs) in delta-full chiral effective field theory ($χ$EFT) up to third order. We use eigenvector continuation for fast and accurate emulation of the likelihood and Hamiltonian Monte Carlo to draw effectively independent samples from the posteriors. Our Bayesian inference is conditioned on the Granada database of neutron-proton ($np$) cross sections and polarizations. We use priors grounded in $χ$EFT assumptions and a Roy-Steiner analysis of pion-nucleon scattering data. We model correlated EFT truncation errors using a two-feature Gaussian process, and find correlation lengths for $np$ scattering energies and angles in the ranges 40--120 MeV and 25--45 degrees, respectively. These correlations yield a non-diagonal covariance matrix and reduce the number of independent scattering data with a factor of eight and four at the second and third chiral orders, respectively. The relatively small difference between the second and third order predictions in delta-full $χ$EFT suppresses the marginal variance of the truncation error and the effects of its correlation structure. Our results are particularly important for analyzing the predictive capabilities in \textit{ab initio} nuclear theory. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 17 pages, 11 figures

Ab Initio Symmetry-Adapted Emulator for Studying Emergent Collectivity and Clustering in Nuclei

Authors: Kevin S. Becker, Kristina D. Launey, Andreas Ekström, Tomas Dytrych

Abstract: We discuss emulators from the ab initio symmetry-adapted no-core shell-model framework for studying the formation of alpha clustering and collective properties without effective charges. We present a new type of an emulator, one that utilizes the eigenvector continuation technique but is based on the use of symplectic symmetry considerations. This is achieved by using physically relevant degrees o… ▽ More We discuss emulators from the ab initio symmetry-adapted no-core shell-model framework for studying the formation of alpha clustering and collective properties without effective charges. We present a new type of an emulator, one that utilizes the eigenvector continuation technique but is based on the use of symplectic symmetry considerations. This is achieved by using physically relevant degrees of freedom, namely, the symmetry-adapted basis, which exploits the almost perfect symplectic symmetry in nuclei. Specifically, we study excitation energies, point-proton root-mean-square radii, along with electric quadrupole moments and transitions for 6Li and 12C. We show that the set of parameterizations of the chiral potential used to train the emulators has no significant effect on predictions of dominant nuclear features, such as shape and the associated symplectic symmetry, along with cluster formation, but slightly varies details that affect collective quadrupole moments, asymptotic normalization coefficients, and alpha partial widths up to a factor of two. This makes these types of emulators important for further constraining the nuclear force for high-precision nuclear structure and reaction observables. △ Less

Submitted 1 March, 2023; originally announced March 2023.

Comments: 11 pages including references, 47 cited references, 1 table, 6 figures

Journal ref: JOURNAL={Frontiers in Physics}, VOLUME={11}, YEAR={2023}, ISSN={2296-424X}

Bayesian Analysis of $χ$EFT at Leading Order in a Modified Weinberg Power Counting Approach

Authors: Oliver Thim, Eleanor May, Andreas Ekström, Christian Forssén

Abstract: We present a Bayesian analysis of renormalization-group invariant nucleon-nucleon interactions at leading order in chiral effective field theory ($χ$EFT) with momentum cutoffs in the range 400--4000 MeV. We use history matching to identify relevant regions in the parameter space of low-energy constants (LECs) and subsequently infer the posterior probability density of their values using Markov cha… ▽ More We present a Bayesian analysis of renormalization-group invariant nucleon-nucleon interactions at leading order in chiral effective field theory ($χ$EFT) with momentum cutoffs in the range 400--4000 MeV. We use history matching to identify relevant regions in the parameter space of low-energy constants (LECs) and subsequently infer the posterior probability density of their values using Markov chain Monte Carlo. All posteriors are conditioned on experimental data for neutron-proton scattering observables and we estimate the $χ$EFT truncation error in an uncorrelated limit. We do not detect any significant cutoff dependence in the posterior predictive distributions for two-nucleon observables. For all cutoff values we find a multimodal LEC posterior with an insignificant mode harboring a bound $^1{S}_0$ state. The $^3P_0$ and $^3P_2$ phase shifts emerging from the Bayesian analysis are less constrained and typically more repulsive compared to the results of a phase shift optimization. We expect that our inference will impact predictions for nuclei. This work demonstrates how to perform inference in the presence of limit-cycle-like behavior and spurious bound states, and lays the foundation for a Bayesian analysis of renormalization-group invariant $χ$EFT interactions beyond leading order. △ Less

Submitted 19 December, 2023; v1 submitted 24 February, 2023; originally announced February 2023.

Comments: 20 pages, 16 figures

Journal ref: Phys. Rev. C 108, 054002 (2023)

What is ab initio in nuclear theory?

Authors: A. Ekström, C. Forssén, G. Hagen, G. R. Jansen, W. Jiang, T. Papenbrock

Abstract: Ab initio has been used as a label in nuclear theory for over two decades. Its meaning has evolved and broadened over the years. We present our interpretation, briefly review its historical use, and discuss its present-day relation to theoretical uncertainty quantification. Ab initio has been used as a label in nuclear theory for over two decades. Its meaning has evolved and broadened over the years. We present our interpretation, briefly review its historical use, and discuss its present-day relation to theoretical uncertainty quantification. △ Less

Submitted 28 February, 2023; v1 submitted 21 December, 2022; originally announced December 2022.

Comments: Contribution to Research Topic "Uncertainty Quantification in Nuclear Physics" of Frontiers in Physics

Journal ref: Front. Phys. 11:1129094 (2023)

Posterior predictive distributions of neutron-deuteron cross sections

Authors: Sean B. S. Miller, Andreas Ekström, Christian Forssén

Abstract: We quantify the posterior predictive distributions (PPDs) of elastic neutron-deuteron ($nd$) scattering cross sections using nucleon-nucleon ($NN$) interactions from chiral effective field theory ($χ$EFT) up to and including next-to-next-to-next-to-leading order (N$^3$LO). These PPDs quantify the spread in $nd$ predictions due to the variability of the low-energy constants (LECs) inferred from… ▽ More We quantify the posterior predictive distributions (PPDs) of elastic neutron-deuteron ($nd$) scattering cross sections using nucleon-nucleon ($NN$) interactions from chiral effective field theory ($χ$EFT) up to and including next-to-next-to-next-to-leading order (N$^3$LO). These PPDs quantify the spread in $nd$ predictions due to the variability of the low-energy constants (LECs) inferred from $NN$ scattering data. We use the wave-packet continuum discretization method to solve the Alt-Grassberger-Sandhas form of the Faddeev equations for elastic scattering. We draw 100 samples from the PPDs of $nd$ cross sections up to 67 MeV in scattering energy, i.e., in the energy region where the effects of three-nucleon forces are expected to be small. We find that the uncertainty about $NN$ LECs inferred from $NN$ scattering data, when assuming uncorrelated errors, does not translate to significant uncertainty in the low-energy $nd$ continuum. Based on our estimates, the uncertainty of $nd$ predictions are dominated by the $χ$EFT truncation error, at least below N$^3$LO. At this order, the 90% credible interval of the PPD and the truncation error are comparable, although both are very small on an absolute scale. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 10 pages, 6 figures

Bayesian estimation of the low-energy constants up to fourth order in the nucleon-nucleon sector of chiral effective field theory

Authors: Isak Svensson, Andreas Ekström, Christian Forssén

Abstract: We use Bayesian methods and Hamiltonian Monte Carlo (HMC) sampling to infer the posterior probability density (PDF) for the low-energy constants (LECs) up to next-to-next-to-next-to-leading order (N3LO) in a chiral effective field theory ($χ$EFT) description of the nucleon-nucleon interaction. In a first step, we condition the inference on neutron-proton and proton-proton scattering data and accou… ▽ More We use Bayesian methods and Hamiltonian Monte Carlo (HMC) sampling to infer the posterior probability density (PDF) for the low-energy constants (LECs) up to next-to-next-to-next-to-leading order (N3LO) in a chiral effective field theory ($χ$EFT) description of the nucleon-nucleon interaction. In a first step, we condition the inference on neutron-proton and proton-proton scattering data and account for uncorrelated $χ$EFT truncation errors. We demonstrate how to successfully sample the 31-dimensional space of LECs at N3LO using a revised HMC inference protocol. In a second step we extend the analysis by means of importance sampling and an empirical determination of the neutron-neutron scattering length to infer the posterior PDF for the leading charge-dependent contact LEC in the $^{1}S_0$ neutron-neutron interaction channel. While doing so we account for the $χ$EFT truncation error via a conjugate prior. We use the resulting posterior PDF to sample the posterior predictive distributions for the effective range parameters in the $^{1}S_0$ wave as well as the strengths of charge-symmetry breaking and charge-independence breaking. We conclude that empirical point-estimate results of isospin breaking in the $^{1}S_0$ channel are consistent with the PDFs obtained in our Bayesian analysis and that, when accounting for $χ$EFT truncation errors, one must go to next-to-next-to-leading order to confidently detect isospin breaking effects. △ Less

Submitted 21 June, 2022; v1 submitted 16 June, 2022; originally announced June 2022.

Comments: 16 pages, 8 figures. One additional figure in the ancillary material. [Minor typo fixed.]

Nuclear Forces for Precision Nuclear Physics -- a collection of perspectives

Authors: Ingo Tews, Zohreh Davoudi, Andreas Ekström, Jason D. Holt, Kevin Becker, Raúl Briceño, David J. Dean, William Detmold, Christian Drischler, Thomas Duguet, Evgeny Epelbaum, Ashot Gasparyan, Jambul Gegelia, Jeremy R. Green, Harald W. Grießhammer, Andrew D. Hanlon, Matthias Heinz, Heiko Hergert, Martin Hoferichter, Marc Illa, David Kekejian, Alejandro Kievsky, Sebastian König, Hermann Krebs, Kristina D. Launey , et al. (20 additional authors not shown)

Abstract: This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges… ▽ More This is a collection of perspective pieces contributed by the participants of the Institute of Nuclear Theory's Program on Nuclear Physics for Precision Nuclear Physics which was held virtually from April 19 to May 7, 2021. The collection represents the reflections of a vibrant and engaged community of researchers on the status of theoretical research in low-energy nuclear physics, the challenges ahead, and new ideas and strategies to make progress in nuclear structure and reaction physics, effective field theory, lattice QCD, quantum information, and quantum computing. The contributed pieces solely reflect the perspectives of the respective authors and do not represent the viewpoints of the Institute for Nuclear theory or the organizers of the program. △ Less

Submitted 2 February, 2022; originally announced February 2022.

Comments: Perspective pieces of the virtual INT program 21-1b "Nuclear Forces for Precision Nuclear Physics", 107 pages

Report number: INT-PUB-22-002, LA-UR-22-20419, UMD-PP-022-02, FERMILAB-PUB-22-090-T

Journal ref: Few-Body Systems 63, 67 (2022)

Neutron-deuteron scattering cross-sections with chiral $NN$ interactions using wave-packet continuum discretization

Authors: Sean B. S. Miller, Andreas Ekström, Kai Hebeler

Abstract: In this work we present a framework that allows to solve the Faddeev equations for three-nucleon scattering using the wave-packet continuum-discretization method. We perform systematic benchmarks using results in the literature and study in detail the convergence of this method with respect to the number of wave packets. We compute several different elastic neutron-deuteron scattering cross-sectio… ▽ More In this work we present a framework that allows to solve the Faddeev equations for three-nucleon scattering using the wave-packet continuum-discretization method. We perform systematic benchmarks using results in the literature and study in detail the convergence of this method with respect to the number of wave packets. We compute several different elastic neutron-deuteron scattering cross-section observables for a variety of energies using chiral nucleon-nucleon interactions. For the interaction $\text{N2LO}_{\text{opt}}$ we find good agreement with data for nucleon scattering-energies $E_\text{Lab} \leq 70$ MeV and a slightly larger maximum of the neutron analyzing power $A_y(n)$ at $E_\text{Lab} = 10$ MeV and 21 MeV compared with other interactions. This work represents a first step towards a systematic inclusion of three-nucleon scattering observables in the construction of next-generation nuclear interactions. △ Less

Submitted 24 January, 2022; originally announced January 2022.

Comments: 17 pages, 7 figures

Ab initio predictions link the neutron skin of ${}^{208}$Pb to nuclear forces

Authors: Baishan Hu, Weiguang Jiang, Takayuki Miyagi, Zhonghao Sun, Andreas Ekström, Christian Forssén, Gaute Hagen, Jason D. Holt, Thomas Papenbrock, S. Ragnar Stroberg, Ian Vernon

Abstract: 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 experim… ▽ More 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 parameterisations 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. △ Less

Submitted 22 August, 2022; v1 submitted 2 December, 2021; originally announced December 2021.

Comments: Published in Nature Physics. 21 pages including Methods and Extended data sections, 3 figures in the main text plus 10 extended data display items. Source data is included as ancillary files

Journal ref: Nature Phys. 18, 1196 (2022)

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

Authors: Chinmay Mishra, A. Ekström, G. Hagen, T. Papenbrock, L. Platter

Abstract: 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. Cu… ▽ More 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. △ Less

Submitted 30 November, 2021; originally announced November 2021.

Comments: 8 pages, 8 figures

Journal ref: Phys. Rev. C 106, 024004 (2022)

Bayesian parameter estimation in $χ$EFT using Hamiltonian Monte Carlo

Authors: Isak Svensson, Andreas Ekström, Christian Forssén

Abstract: The number of low-energy constants (LECs) in chiral effective field theory ($χ$EFT) grows rapidly with increasing chiral order, necessitating the use of Markov chain Monte Carlo techniques for sampling their posterior probability density function. For this we introduce a Hamiltonian Monte Carlo (HMC) algorithm and sample the LEC posterior up to next-to-next-to-leading order (NNLO) in the two-nucle… ▽ More The number of low-energy constants (LECs) in chiral effective field theory ($χ$EFT) grows rapidly with increasing chiral order, necessitating the use of Markov chain Monte Carlo techniques for sampling their posterior probability density function. For this we introduce a Hamiltonian Monte Carlo (HMC) algorithm and sample the LEC posterior up to next-to-next-to-leading order (NNLO) in the two-nucleon sector of $χ$EFT. We find that the sampling efficiency of HMC is three to six times higher compared to an affine-invariant sampling algorithm. We analyze the empirical coverage probability and validate that the NNLO model yields predictions for two-nucleon scattering data with largely reliable credible intervals, provided that one ignores the leading order EFT expansion parameter when inferring the variance of the truncation error. We also find that the NNLO truncation error dominates the error budget. △ Less

Submitted 8 October, 2021; originally announced October 2021.

Comments: 22 pages, 12 figures

The importance of few-nucleon forces in chiral effective field theory

Authors: C. -J. Yang, A. Ekström, C. Forssén, G. Hagen, G. Rupak, U. van Kolck

Abstract: We study the importance of few-nucleon forces in chiral effective field theory for describing many-nucleon systems. A combinatorial argument suggests that three-nucleon forces -- which are conventionally regarded as next-to-next-to-leading order -- should accompany the two-nucleon force already at leading order (LO) starting with mass number $A\approx 10-20$. We find that this promotion enables th… ▽ More We study the importance of few-nucleon forces in chiral effective field theory for describing many-nucleon systems. A combinatorial argument suggests that three-nucleon forces -- which are conventionally regarded as next-to-next-to-leading order -- should accompany the two-nucleon force already at leading order (LO) starting with mass number $A\approx 10-20$. We find that this promotion enables the first realistic description of the $^{16}$O ground state based on a renormalization-group-invariant LO interaction. We also performed coupled-cluster calculations of the equation of state for symmetric nuclear matter and our results indicate that LO four-nucleon forces could play a crucial role for describing heavy-mass nuclei. The enhancement mechanism we found is very general and could be important also in other many-body problems. △ Less

Submitted 5 October, 2023; v1 submitted 27 September, 2021; originally announced September 2021.

Comments: 6 pages, 4 figures

Fast & rigorous predictions for $A=6$ nuclei with Bayesian posterior sampling

Authors: T. Djärv, A. Ekström, C. Forssén, H. T. Johansson

Abstract: We make ab initio predictions for the A = 6 nuclear level scheme based on two- and three-nucleon interactions up to next-to-next-to-leading order in chiral effective field theory ($χ$EFT). We utilize eigenvector continuation and Bayesian methods to quantify uncertainties stemming from the many-body method, the $χ$EFT truncation, and the low-energy constants of the nuclear interaction. The construc… ▽ More We make ab initio predictions for the A = 6 nuclear level scheme based on two- and three-nucleon interactions up to next-to-next-to-leading order in chiral effective field theory ($χ$EFT). We utilize eigenvector continuation and Bayesian methods to quantify uncertainties stemming from the many-body method, the $χ$EFT truncation, and the low-energy constants of the nuclear interaction. The construction and validation of emulators is made possible via the development of JupiterNCSM -- a new M-scheme no-core shell model code that uses on-the-fly Hamiltonian matrix construction for efficient, single-node computations up to $N_\mathrm{max} = 10$ for ${}^{6}\mathrm{Li}$. We find a slight underbinding of ${}^{6}\mathrm{He}$ and ${}^{6}\mathrm{Li}$, although consistent with experimental data given our theoretical error bars. As a result of incorporating a correlated $χ$EFT-truncation errors we find more precise predictions (smaller error bars) for separation energies: $S_d({}^{6}\mathrm{Li}) = 0.89 \pm 0.44$ MeV, $S_{2n}({}^{6}\mathrm{He}) = 0.20 \pm 0.60$ MeV, and for the beta decay Q-value: $Q_{β^-}({}^{6}\mathrm{He}) = 3.71 \pm 0.65$ MeV. We conclude that our error bars can potentially be reduced further by extending the model space used by JupiterNCSM. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 13 pages, 9 figures

Journal ref: Phys. Rev. C 105, 014005 (2022)

Wave-packet continuum discretisation for nucleon-nucleon scattering predictions

Authors: Sean B. S. Miller, Andreas Ekström, Christian Forssén

Abstract: In this paper we analyse the efficiency, precision, and accuracy of computing elastic nucleon-nucleon (NN) scattering amplitudes with the wave-packet continuum discretisation method (WPCD). This method provides approximate scattering solutions at multiple scattering energies simultaneously. We therefore utilise a graphics processing unit (GPU) to explore the benefits of this inherent parallelism.… ▽ More In this paper we analyse the efficiency, precision, and accuracy of computing elastic nucleon-nucleon (NN) scattering amplitudes with the wave-packet continuum discretisation method (WPCD). This method provides approximate scattering solutions at multiple scattering energies simultaneously. We therefore utilise a graphics processing unit (GPU) to explore the benefits of this inherent parallelism. From a theoretical perspective, the WPCD method promises a speedup compared to a standard matrix-inversion method. We use the chiral NNLO$_{\rm opt}$ interaction to demonstrate that WPCD enables efficient computation of NN scattering amplitudes provided one can tolerate an averaged method error of $~1-5$ mb in the total cross section at scattering energies $0-350$ MeV in the laboratory frame of reference. Considering only scattering energies $\sim40-350$ MeV, we find a smaller method error of $\lesssim 1-2$ mb. By increasing the number of wave-packets we can further reduce the overall method error. However, the parallel leverage of the WPCD method will be offset by the increased size of the resulting discretisation mesh. In practice, a GPU-implementation is mainly advantageous for matrices that fit in the fast on-chip shared memory. We find that WPCD is a promising method for computationally efficient, statistical analyses of nuclear interactions from effective field theory, where we can utilise Bayesian inference methods to incorporate relevant uncertainties. △ Less

Submitted 17 June, 2022; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: 33 pages, 9 figures, v2 corresponds to published version

Journal ref: J.Phys.G 49 (2022) 2, 024001

Normal-ordering approximations and translational (non) invariance

Authors: T. Djärv, A. Ekström, C. Forssén, G. R. Jansen

Abstract: Normal-ordering provides an approach to approximate three-body forces as effective two-body operators and it is therefore an important tool in many-body calculations with realistic nuclear interactions. The corresponding neglect of certain three-body terms in the normal-ordered Hamiltonian is known to influence translational invariance, although the magnitude of this effect has not yet been system… ▽ More Normal-ordering provides an approach to approximate three-body forces as effective two-body operators and it is therefore an important tool in many-body calculations with realistic nuclear interactions. The corresponding neglect of certain three-body terms in the normal-ordered Hamiltonian is known to influence translational invariance, although the magnitude of this effect has not yet been systematically quantified. In this work we study in particular the normal-ordering two-body approximation applied to a single harmonic-oscillator reference state. We explicate the breaking of translational invariance and demonstrate the magnitude of the approximation error as a function of model space parameters for $^4\rm{He}$ and $^{16}\rm{O}$ by performing full no-core shell-model calculations with and without three-nucleon forces. We combine two different diagnostics to better monitor the breaking of translational invariance. While the center-of-mass effect is shown to become potentially very large for $^4\rm{He}$, it is also shown to be much smaller for $^{16}\rm{O}$ although full convergence is not reached. These tools can be easily implemented in studies using other many-body frameworks and bases. △ Less

Submitted 30 April, 2021; originally announced April 2021.

Comments: 13 pages, 10 figures

Journal ref: Phys. Rev. C 104, 024324 (2021)

Rigorous constraints on three-nucleon forces in chiral effective field theory from fast and accurate calculations of few-body observables

Authors: S. Wesolowski, I. Svensson, A. Ekström, C. Forssén, R. J. Furnstahl, J. A. Melendez, D. R. Phillips

Abstract: We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory ($χ$EFT) that are provided by bound-state observables in the $A=3$ and $A=4$ sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the $χ$EFT expansion at next-to-next-to-leading order. A consistent solution… ▽ More We explore the constraints on the three-nucleon force (3NF) of chiral effective field theory ($χ$EFT) that are provided by bound-state observables in the $A=3$ and $A=4$ sectors. Our statistically rigorous analysis incorporates experimental error, computational method uncertainty, and the uncertainty due to truncation of the $χ$EFT expansion at next-to-next-to-leading order. A consistent solution for the ${}^3$H binding energy, the ${}^4$He binding energy and radius, and the ${}^3$H $β$-decay rate can only be obtained if $χ$EFT truncation errors are included in the analysis. All of these except the $β$-decay rate give essentially degenerate constraints on the 3NF low-energy constants, so it is crucial for estimating these parameters. We use eigenvector continuation for fast and accurate emulation of No-Core Shell Model calculations of the considered few-nucleon observables. This facilitates sampling of the posterior probability distribution, allowing us to also determine the distributions of the hyperparameters that quantify the truncation error. We find a $χ$EFT expansion parameter of $Q=0.33 \pm 0.06$ for these observables. △ Less

Submitted 11 January, 2022; v1 submitted 9 April, 2021; originally announced April 2021.

Comments: 15 pages, 7 figures; same as published version

Charge radii of exotic potassium isotopes challenge nuclear theory and the magic character of $N = 32$

Authors: Á. Koszorús, X. F. Yang, W. G. Jiang, S. J. Novario, S. W. Bai, J. Billowes, C. L. Binnersley, M. L. Bissell, T. E. Cocolios, B. S. Cooper, R. P. de Groote, A. Ekström, K. T. Flanagan, C. Forssén, S. Franchoo, R. F. Garcia Ruiz, F. P. Gustafsson, G. Hagen, G. R. Jansen, A. Kanellakopoulos, M. Kortelainen, W. Nazarewicz, G. Neyens, T. Papenbrock, P. -G. Reinhard , et al. (4 additional authors not shown)

Abstract: Nuclear charge radii are sensitive probes of different aspects of the nucleon-nucleon interaction and the bulk properties of nuclear matter; thus, they provide a stringent test and challenge for nuclear theory. The calcium region has been of particular interest, as experimental evidence has suggested a new magic number at $N = 32$ [1-3], while the unexpectedly large increases in the charge radii [… ▽ More Nuclear charge radii are sensitive probes of different aspects of the nucleon-nucleon interaction and the bulk properties of nuclear matter; thus, they provide a stringent test and challenge for nuclear theory. The calcium region has been of particular interest, as experimental evidence has suggested a new magic number at $N = 32$ [1-3], while the unexpectedly large increases in the charge radii [4,5] open new questions about the evolution of nuclear size in neutron-rich systems. By combining the collinear resonance ionization spectroscopy method with $β$-decay detection, we were able to extend the charge radii measurement of potassium ($Z =19$) isotopes up to the exotic $^{52}$K ($t_{1/2}$ = 110 ms), produced in minute quantities. Our work provides the first charge radii measurement beyond $N = 32$ in the region, revealing no signature of the magic character at this neutron number. The results are interpreted with two state-of-the-art nuclear theories. For the first time, a long sequence of isotopes could be calculated with coupled-cluster calculations based on newly developed nuclear interactions. The strong increase in the charge radii beyond $N = 28$ is not well captured by these calculations, but is well reproduced by Fayans nuclear density functional theory, which, however, overestimates the odd-even staggering effect. These findings highlight our limited understanding on the nuclear size of neutron-rich systems, and expose pressing problems that are present in some of the best current models of nuclear theory. △ Less

Submitted 3 December, 2020; originally announced December 2020.

Comments: submitted version; revision accepted in Nature Physics

Journal ref: Nature Physics 17, 439 (2021)

Power counting in chiral effective field theory and nuclear binding

Authors: C. -J. Yang, A. Ekström, C. Forssén, G. Hagen

Abstract: Chiral effective field theory ($χ$EFT), as originally proposed by Weinberg, promises a theoretical connection between low-energy nuclear interactions and quantum chromodynamics (QCD). However, the important property of renormalization-group (RG) invariance is not fulfilled in current implementations and its consequences for predicting atomic nuclei beyond two- and three-nucleon systems has remaine… ▽ More Chiral effective field theory ($χ$EFT), as originally proposed by Weinberg, promises a theoretical connection between low-energy nuclear interactions and quantum chromodynamics (QCD). However, the important property of renormalization-group (RG) invariance is not fulfilled in current implementations and its consequences for predicting atomic nuclei beyond two- and three-nucleon systems has remained unknown. In this work we present a first and systematic study of recent RG-invariant formulations of $χ$EFT and their predictions for the binding energies and other observables of selected nuclear systems with mass-numbers up to $A =16$. Specifically, we have carried out ab initio no-core shell-model and coupled cluster calculations of the ground-state energy of $^3$H, $^{3,4}$He, $^{6}$Li, and $^{16}$O using several recent power-counting (PC) schemes at leading order (LO) and next-to-leading order (NLO), where the subleading interactions are treated in perturbation theory. Our calculations indicate that RG-invariant and realistic predictions can be obtained for nuclei with mass number $A \leq 4$. We find, however, that $^{16}$O is either unbound with respect to the four $α$-particle threshold, or deformed, or both. Similarly, we find that the $^{6}$Li ground-state resides above the $α$-deuteron separation threshold. These results are in stark contrast with experimental data and point to either necessary fine-tuning of all relevant counterterms, or that current state-of-the-art RG-invariant PC schemes at LO in $χ$EFT lack necessary diagrams -- such as three-nucleon forces -- to realistically describe nuclei with mass number $A>4$. △ Less

Submitted 13 May, 2021; v1 submitted 23 November, 2020; originally announced November 2020.

Comments: 18 pages, 12 figures, published version

Journal ref: Phys. Rev. C 103, 054304 (2021)

Accurate bulk properties of nuclei from $A = 2$ to $\infty$ from potentials with $Δ$ isobars

Authors: W. G. Jiang, A. Ekström, C. Forssén, G. Hagen, G. R. Jansen, T. Papenbrock

Abstract: We optimize $Δ$-full nuclear interactions from chiral effective field theory. The low-energy constants of the contact potentials are constrained by two-body scattering phase shifts, and by properties of bound-state of $A=2$ to $4$ nucleon systems and nuclear matter. The pion-nucleon couplings are taken from a Roy-Steiner analysis. The resulting interactions yield accurate binding energies and radi… ▽ More We optimize $Δ$-full nuclear interactions from chiral effective field theory. The low-energy constants of the contact potentials are constrained by two-body scattering phase shifts, and by properties of bound-state of $A=2$ to $4$ nucleon systems and nuclear matter. The pion-nucleon couplings are taken from a Roy-Steiner analysis. The resulting interactions yield accurate binding energies and radii for a range of nuclei from $A=16$ to $A=132$, and provide accurate equations of state for nuclear matter and realistic symmetry energies. Selected excited states are also in agreement with data. △ Less

Submitted 9 November, 2020; v1 submitted 30 June, 2020; originally announced June 2020.

Comments: 7 pages, 6 figures

Journal ref: Phys. Rev. C 102, 054301 (2020)

New Ideas in Constraining Nuclear Forces

Authors: I. Tews, Z. Davoudi, A. Ekström, J. D. Holt, J. E. Lynn

Abstract: In recent years, nuclear physics has benefited greatly from the development of powerful ab initio many-body methods and their combination with interactions from chiral effective field theory. With increasing computational power and continuous development of these methods, we are entering an era of precision nuclear physics. Indeed, uncertainties from nuclear Hamiltonians now dominate over uncertai… ▽ More In recent years, nuclear physics has benefited greatly from the development of powerful ab initio many-body methods and their combination with interactions from chiral effective field theory. With increasing computational power and continuous development of these methods, we are entering an era of precision nuclear physics. Indeed, uncertainties from nuclear Hamiltonians now dominate over uncertainties from many-body methods. This review summarizes the current status of, and future directions in, deriving and constraining nuclear Hamiltonians. △ Less

Submitted 16 August, 2020; v1 submitted 10 January, 2020; originally announced January 2020.

Comments: 44 pages, 16 figures, Topical Review of ECT* workshop "New Ideas in Constraining Nuclear Forces", published version

Report number: LA-UR-19-32540, UMD-PP-019-09

Journal ref: J. Phys. G: Nucl. Part. Phys. 47, 103001 (2020)

Ab initio models of atomic nuclei: challenges and new ideas

Authors: Andreas Ekström

Abstract: This review presents some of the challenges in constructing models of atomic nuclei starting from theoretical descriptions of the strong interaction between nucleons. The focus is on statistical computing and methods for analyzing the link between bulk properties of atomic nuclei, such as radii and binding energies, and the underlying microscopic description of the nuclear interaction. The importa… ▽ More This review presents some of the challenges in constructing models of atomic nuclei starting from theoretical descriptions of the strong interaction between nucleons. The focus is on statistical computing and methods for analyzing the link between bulk properties of atomic nuclei, such as radii and binding energies, and the underlying microscopic description of the nuclear interaction. The importance of careful model calibration and uncertainty quantification of theoretical predictions is highlighted. △ Less

Submitted 4 December, 2019; originally announced December 2019.

Comments: 25 pages, 5 figures, Submitted to Frontiers in Physics (Research topic) The Long-Lasting Quest for Nuclear Interactions: The Past, the Present and the Future

Improved many-body expansions from eigenvector continuation

Authors: Pepijn Demol, Thomas Duguet, Andreas Ekström, Mikael Frosini, Kai Hebeler, Sebastian König, Dean Lee, Achim Schwenk, Vittorio Somà, Alexander Tichai

Abstract: Quantum many-body theory has witnessed tremendous progress in various fields, ranging from atomic and solid-state physics to quantum chemistry and nuclear structure. Due to the inherent computational burden linked to the ab initio treatment of microscopic fermionic systems, it is desirable to obtain accurate results through low-order perturbation theory. In atomic nuclei however, effects such as s… ▽ More Quantum many-body theory has witnessed tremendous progress in various fields, ranging from atomic and solid-state physics to quantum chemistry and nuclear structure. Due to the inherent computational burden linked to the ab initio treatment of microscopic fermionic systems, it is desirable to obtain accurate results through low-order perturbation theory. In atomic nuclei however, effects such as strong short-range repulsion between nucleons can spoil the convergence of the expansion and make the reliability of perturbation theory unclear. Mathematicians have devised an extensive machinery to overcome the problem of divergent expansions by making use of so-called resummation methods. In large-scale many-body applications such schemes are often of limited use since no a priori analytical knowledge of the expansion is available. We present here eigenvector continuation as an alternative resummation tool that is both efficient and reliable because it is based on robust and simple mathematical principles. △ Less

Submitted 28 November, 2019; originally announced November 2019.

Comments: 6 pages, 2 figures

Journal ref: Phys. Rev. C 101, 041302 (2020)

Global sensitivity analysis of bulk properties of an atomic nucleus

Authors: Andreas Ekström, Gaute Hagen

Abstract: We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus $^{16}$O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose we develop a subspace-projected coupled-cluster method using eigenvector continuation [Frame D. et al., Phys. Rev. Lett. 121, 03… ▽ More We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus $^{16}$O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose we develop a subspace-projected coupled-cluster method using eigenvector continuation [Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)]. With this method we compute the binding energy and charge radius of $^{16}$O at more than one million different values of the 16 low-energy constants in one hour on a standard laptop. For relatively small subspace projections, the root-mean-square error is about 1% compared to full space coupled-cluster results. We find that 58(1)% of the variance in the energy can be apportioned to a single contact-term in the $^3S_1$-wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation, and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up for an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segré chart. △ Less

Submitted 7 October, 2019; originally announced October 2019.

Comments: 6 pages, 3 figures, Supplemental Material provided as ancillary material

Journal ref: Phys. Rev. Lett. 123, 252501 (2019)

Eigenvector Continuation as an Efficient and Accurate Emulator for Uncertainty Quantification

Authors: S. König, A. Ekström, K. Hebeler, D. Lee, A. Schwenk

Abstract: First principles calculations of atomic nuclei based on microscopic nuclear forces derived from chiral effective field theory (EFT) have blossomed in the past years. A key element of such ab initio studies is the understanding and quantification of systematic and statistical errors arising from the omission of higher-order terms in the chiral expansion as well as the model calibration. While there… ▽ More First principles calculations of atomic nuclei based on microscopic nuclear forces derived from chiral effective field theory (EFT) have blossomed in the past years. A key element of such ab initio studies is the understanding and quantification of systematic and statistical errors arising from the omission of higher-order terms in the chiral expansion as well as the model calibration. While there has been significant progress in analyzing theoretical uncertainties for nucleon-nucleon scattering observables, the generalization to multi-nucleon systems has not been feasible yet due to the high computational cost of evaluating observables for a large set of low-energy couplings. In this Letter we show that a new method called eigenvector continuation (EC) can be used for constructing an efficient and accurate emulator for nuclear many-body observables, thereby enabling uncertainty quantification in multi-nucleon systems. We demonstrate the power of EC emulation with a proof-of-principle calculation that lays out all correlations between bulk ground-state observables in the few-nucleon sector. On the basis of ab initio calculations for the ground-state energy and radius in 4He, we demonstrate that EC is more accurate and efficient compared to established methods like Gaussian processes. △ Less

Submitted 15 October, 2020; v1 submitted 18 September, 2019; originally announced September 2019.

Comments: 8 pages, 6 figures, Python code and input files provided as ancillary material, published version

Journal ref: Phys. Lett. B 810 (2020) 135814

A statistical analysis of the nuclear structure uncertainties in $μ$D

Authors: Oscar J. Hernandez, Sonia Bacca, Nir Barnea, Nir Nevo-Dinur, Andreas Ekström, Chen Ji

Abstract: The charge radius of the deuteron (D), was recently determined to three times the precision compared with previous measurements using the measured Lamb shift in muonic deuterium (muD). However, the muD value is 5.6 $σ$ smaller than the world averaged CODATA-2014 value [1]. To shed light on this discrepancy we analyze the uncertainties of the nuclear structure calculations of the Lamb shift in muD… ▽ More The charge radius of the deuteron (D), was recently determined to three times the precision compared with previous measurements using the measured Lamb shift in muonic deuterium (muD). However, the muD value is 5.6 $σ$ smaller than the world averaged CODATA-2014 value [1]. To shed light on this discrepancy we analyze the uncertainties of the nuclear structure calculations of the Lamb shift in muD and conclude that nuclear theory uncertainty is not likely to be the source of the discrepancy. △ Less

Submitted 6 March, 2019; originally announced March 2019.

Comments: 4 pages, 1 figure, accepted proceedings of the "XXII International Conference on Few-Body Problems in Physics", 9-13 July 2018

Report number: MITP/19-017

Bayesian optimization in ab initio nuclear physics

Authors: A. Ekström, C. Forssén, C. Dimitrakakis, D. Dubhashi, H. T. Johansson, A. S. Muhammad, H. Salomonsson, A. Schliep

Abstract: Theoretical models of the strong nuclear interaction contain unknown coupling constants (parameters) that must be determined using a pool of calibration data. In cases where the models are complex, leading to time consuming calculations, it is particularly challenging to systematically search the corresponding parameter domain for the best fit to the data. In this paper, we explore the prospect of… ▽ More Theoretical models of the strong nuclear interaction contain unknown coupling constants (parameters) that must be determined using a pool of calibration data. In cases where the models are complex, leading to time consuming calculations, it is particularly challenging to systematically search the corresponding parameter domain for the best fit to the data. In this paper, we explore the prospect of applying Bayesian optimization to constrain the coupling constants in chiral effective field theory descriptions of the nuclear interaction. We find that Bayesian optimization performs rather well with low-dimensional parameter domains and foresee that it can be particularly useful for optimization of a smaller set of coupling constants. A specific example could be the determination of leading three-nucleon forces using data from finite nuclei or three-nucleon scattering experiments. △ Less

Submitted 3 February, 2019; originally announced February 2019.

Comments: 33 pages, 14 figures

Simulations of Subatomic Many-Body Physics on a Quantum Frequency Processor

Authors: Hsuan-Hao Lu, Natalie Klco, Joseph M. Lukens, Titus D. Morris, Aaina Bansal, Andreas Ekström, Gaute Hagen, Thomas Papenbrock, Andrew M. Weiner, Martin J. Savage, Pavel Lougovski

Abstract: Simulating complex many-body quantum phenomena is a major scientific impetus behind the development of quantum computing, and a range of technologies are being explored to address such systems. We present the results of the largest photonics-based simulation to date, applied in the context of subatomic physics. Using an all-optical quantum frequency processor, the ground-state energies of light nu… ▽ More Simulating complex many-body quantum phenomena is a major scientific impetus behind the development of quantum computing, and a range of technologies are being explored to address such systems. We present the results of the largest photonics-based simulation to date, applied in the context of subatomic physics. Using an all-optical quantum frequency processor, the ground-state energies of light nuclei including the triton ($^3$H), $^{3}$He, and the alpha particle ($^{4}$He) are computed. Complementing these calculations and utilizing a 68-dimensional Hilbert space, our photonic simulator is used to perform sub-nucleon calculations of the two-body and three-body forces between heavy mesons in the Schwinger model. This work is a first step in simulating subatomic many-body physics on quantum frequency processors---augmenting classical computations that bridge scales from quarks to nuclei. △ Less

Submitted 9 October, 2018; originally announced October 2018.

Comments: 23 pages, 5 figures, 20 pages supplemental material

Journal ref: Phys. Rev. A 100, 012320 (2019)

Effective-field-theory predictions of the muon-deuteron capture rate

Authors: Bijaya Acharya, Andreas Ekström, Lucas Platter

Abstract: We quantify the theoretical uncertainties of chiral effective-field-theory predictions of the muon-deuteron capture rate. Theoretical error estimates of this low-energy process is important for a reliable interpretation of forthcoming experimental results by the MuSun collaboration. Specifically, we estimate the three dominant sources of uncertainties that impact theoretical calculations of this r… ▽ More We quantify the theoretical uncertainties of chiral effective-field-theory predictions of the muon-deuteron capture rate. Theoretical error estimates of this low-energy process is important for a reliable interpretation of forthcoming experimental results by the MuSun collaboration. Specifically, we estimate the three dominant sources of uncertainties that impact theoretical calculations of this rate: those resulting from uncertainties in the pool of fit data used to constrain the coupling constants in the nuclear interaction, those due to the truncation of the effective field theory, and those due to uncertainties in the axial radius of the nucleon. For the capture rate into the ${}^1S_0$ channel, we find an uncertainty of approximately $4.6~s^{-1}$ due to the truncation in the effective field theory and an uncertainty of $3.9~s^{-1}$ due to the uncertainty in the axial radius of the nucleon, both of which are similar in size to the targeted experimental precision. △ Less

Submitted 3 January, 2019; v1 submitted 25 June, 2018; originally announced June 2018.

Comments: 4 pages, 2 figures

Journal ref: Phys. Rev. C 98, 065506 (2018)

Pion-less effective field theory for atomic nuclei and lattice nuclei

Authors: A. Bansal, S. Binder, A. Ekström, G. Hagen, G. R. Jansen, T. Papenbrock

Abstract: We compute the medium-mass nuclei $^{16}$O and $^{40}$Ca using pionless effective field theory (EFT) at next-to-leading order (NLO). The low-energy coefficients of the EFT Hamiltonian are adjusted to experimantal data for nuclei with mass numbers $A=2$ and $3$, or alternatively to results from lattice quantum chromodynamics (QCD) at an unphysical pion mass of 806 MeV. The EFT is implemented throug… ▽ More We compute the medium-mass nuclei $^{16}$O and $^{40}$Ca using pionless effective field theory (EFT) at next-to-leading order (NLO). The low-energy coefficients of the EFT Hamiltonian are adjusted to experimantal data for nuclei with mass numbers $A=2$ and $3$, or alternatively to results from lattice quantum chromodynamics (QCD) at an unphysical pion mass of 806 MeV. The EFT is implemented through a discrete variable representation in the harmonic oscillator basis. This approach ensures rapid convergence with respect to the size of the model space and facilitates the computation of medium-mass nuclei. At NLO the nuclei $^{16}$O and $^{40}$Ca are bound with respect to decay into alpha particles. Binding energies per nucleon are 9-10 MeV and 30-40 MeV at pion masses of 140 MeV and 806 MeV, respectively. △ Less

Submitted 29 December, 2017; originally announced December 2017.

Comments: 26 pages

Journal ref: Phys. Rev. C 98, 054301 (2018)

The deuteron-radius puzzle is alive: A new analysis of nuclear structure uncertainties

Authors: Oscar Javier Hernandez, Andreas Ekström, Nir Nevo Dinur, Chen Ji, Sonia Bacca, Nir Barnea

Abstract: To shed light on the deuteron radius puzzle we analyze the theoretical uncertainties of the nuclear structure corrections to the Lamb shift in muonic deuterium. We find that the discrepancy between the calculated two-photon exchange correction and the corresponding experimentally inferred value by Pohl et al. [1] remain. The present result is consistent with our previous estimate, although the dis… ▽ More To shed light on the deuteron radius puzzle we analyze the theoretical uncertainties of the nuclear structure corrections to the Lamb shift in muonic deuterium. We find that the discrepancy between the calculated two-photon exchange correction and the corresponding experimentally inferred value by Pohl et al. [1] remain. The present result is consistent with our previous estimate, although the discrepancy is reduced from 2.6 $σ$ to 2 $σ$. The error analysis includes statistic as well as systematic uncertainties stemming from the use of nucleon-nucleon interactions derived from chiral effective field theory at various orders. We therefore conclude that nuclear theory uncertainty is more likely not the source of the discrepancy. △ Less

Submitted 11 June, 2018; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: updated according to referee's comments

Corrections to nucleon capture cross sections computed in truncated Hilbert spaces

Authors: B. Acharya, A. Ekström, D. Odell, T. Papenbrock, L. Platter

Abstract: Nucleon capture cross sections enter various astrophysical processes. The measurement of proton capture on nuclei at astrophysically relevant low energies is a challenge, and theoretical computations in this long-wavelength regime are sensitive to the long-distance asymptotics of the wave functions. A theoretical foundation for estimating and correcting errors introduced in capture cross sections… ▽ More Nucleon capture cross sections enter various astrophysical processes. The measurement of proton capture on nuclei at astrophysically relevant low energies is a challenge, and theoretical computations in this long-wavelength regime are sensitive to the long-distance asymptotics of the wave functions. A theoretical foundation for estimating and correcting errors introduced in capture cross sections due to Hilbert space truncation has so far been lacking. We derive extrapolation formulas that relate the infrared regularized capture amplitudes to the infinite basis limit and demonstrate their efficacy for proton-proton fusion. Our results are thus relevant to current calculations of few-body capture reactions such as proton-proton fusion or proton capture on the deuteron, and they also open the way for the use of {\it ab initio} many-body wave functions represented in finite Hilbert spaces in precision calculations of nucleon capture on heavier nuclei. △ Less

Submitted 1 March, 2017; v1 submitted 16 August, 2016; originally announced August 2016.

Comments: 6 pages, 4 figures, published version

Journal ref: Phys. Rev. C 95, 031301(R) (2017)

Uncertainty quantification for proton-proton fusion in chiral effective field theory

Authors: B. Acharya, B. D. Carlsson, A. Ekström, C. Forssén, L. Platter

Abstract: We compute the $S$-factor of the proton-proton ($pp$) fusion reaction using chiral effective field theory ($χ$EFT) up to next-to-next-to-leading order (NNLO) and perform a rigorous uncertainty analysis of the results. We quantify the uncertainties due to (i) the computational method used to compute the $pp$ cross section in momentum space, (ii) the statistical uncertainties in the low-energy coupl… ▽ More We compute the $S$-factor of the proton-proton ($pp$) fusion reaction using chiral effective field theory ($χ$EFT) up to next-to-next-to-leading order (NNLO) and perform a rigorous uncertainty analysis of the results. We quantify the uncertainties due to (i) the computational method used to compute the $pp$ cross section in momentum space, (ii) the statistical uncertainties in the low-energy coupling constants of $χ$EFT, (iii) the systematic uncertainty due to the $χ$EFT cutoff, and (iv) systematic variations in the database used to calibrate the nucleon-nucleon interaction. We also examine the robustness of the polynomial extrapolation procedure, which is commonly used to extract the threshold $S$-factor and its energy-derivatives. By performing a statistical analysis of the polynomial fit of the energy-dependent $S$-factor at several different energy intervals, we eliminate a systematic uncertainty that can arise from the choice of the fit interval in our calculations. In addition, we explore the statistical correlations between the $S$-factor and few-nucleon observables such as the binding energies and point-proton radii of $^{2,3}$H and $^3$He as well as the $D$-state probability and quadrupole moment of $^2$H, and the $β$-decay of $^{3}$H. We find that, with the state-of-the-art optimization of the nuclear Hamiltonian, the statistical uncertainty in the threshold $S$-factor cannot be reduced beyond 0.7%. △ Less

Submitted 26 July, 2016; v1 submitted 4 March, 2016; originally announced March 2016.

Comments: peer-reviewed version: some passages modified, typographic errors fixed

Journal ref: Physics Letters B 760 (2016) 584-589

Effective field theory in the harmonic oscillator basis

Authors: S. Binder, A. Ekström, G. Hagen, T. Papenbrock, K. A. Wendt

Abstract: 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 eva… ▽ More 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. △ Less

Submitted 29 April, 2016; v1 submitted 11 December, 2015; originally announced December 2015.

Comments: 16 pages, 6 figures; published version

Journal ref: Phys. Rev. C 93, 044332 (2016)

Charge, neutron, and weak size of the atomic nucleus

Authors: G. Hagen, A. Ekström, C. Forssén, G. R. Jansen, W. Nazarewicz, T. Papenbrock, K. A. Wendt, S. Bacca, N. Barnea, B. Carlsson, C. Drischler, K. Hebeler, M. Hjorth-Jensen, M. Miorelli, G. Orlandini, A. Schwenk, J. Simonis

Abstract: 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… ▽ More 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. △ Less

Submitted 23 September, 2015; originally announced September 2015.

Journal ref: Nature Physics 12, 186-190 (2016)

Uncertainty analysis and order-by-order optimization of chiral nuclear interactions

Authors: B. D. Carlsson, A. Ekström, C. Forssén, D. Fahlin Strömberg, G. R. Jansen, O. Lilja, M. Lindby, B. A. Mattsson, K. A. Wendt

Abstract: Chiral effective field theory (chiEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, chiEFT is able to provide well-founded estimates of statistical and systematic uncertainties -- although this unique advantage has not yet been fully exploited. We fill this gap by performing an optimization and statistical analysis of all the low-energy constants (LECs) up to nex… ▽ More Chiral effective field theory (chiEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, chiEFT is able to provide well-founded estimates of statistical and systematic uncertainties -- although this unique advantage has not yet been fully exploited. We fill this gap by performing an optimization and statistical analysis of all the low-energy constants (LECs) up to next-to-next-to-leading order. Our optimization protocol corresponds to a simultaneous fit to scattering and bound-state observables in the pion-nucleon, nucleon-nucleon, and few-nucleon sectors, thereby utilizing the full model capabilities of chiEFT. We study the effect on other observables by demonstrating error-propagation methods that can easily be adopted by future works. We employ mathematical optimization and implement automatic differentiation to attain efficient and machine-precise first- and second-order derivatives of the objective function with respect to the LECs. We use power-counting arguments to estimate the systematic uncertainty that is inherent to chiEFT and we construct chiral interactions at different orders with quantified uncertainties. Statistical error propagation is compared with Monte Carlo sampling showing that statistical errors are in general small compared to systematic ones. In conclusion, we find that a simultaneous fit to different sets of data is critical to (i) identify the optimal set of LECs, (ii) capture all relevant correlations, (iii) reduce the statistical uncertainty, and (iv) attain order-by-order convergence in chiEFT. Furthermore, certain systematic uncertainties in the few-nucleon sector are shown to get substantially magnified in the many-body sector; in particlar when varying the cutoff in the chiral potentials. The methodology and results presented in this Paper open a new frontier for uncertainty quantification in ab initio nuclear theory. △ Less

Submitted 27 January, 2016; v1 submitted 8 June, 2015; originally announced June 2015.

Comments: Accepted for publication in Physical Review X. This version includes Supplemental Material

Journal ref: Phys. Rev. X 6, 011019 (2016)

Uncertainty Quantification of the Pion-Nucleon Low-Energy Coupling Constants up to Fourth Order in Chiral Perturbation Theory

Authors: K. A. Wendt, B. D. Carlsson, A. Ekström

Abstract: We extract the statistical uncertainties for the pion-nucleon ($πN$) low energy constants (LECs) up to fourth order $\mathcal{O}(Q^4)$ in the chiral expansion of the nuclear effective Lagrangian. The LECs are optimized with respect to experimental scattering data. For comparison, we also present an uncertainty quantification that is based solely on \pin{} scattering phase shifts. Statistical error… ▽ More We extract the statistical uncertainties for the pion-nucleon ($πN$) low energy constants (LECs) up to fourth order $\mathcal{O}(Q^4)$ in the chiral expansion of the nuclear effective Lagrangian. The LECs are optimized with respect to experimental scattering data. For comparison, we also present an uncertainty quantification that is based solely on \pin{} scattering phase shifts. Statistical errors on the LECs are critical in order to estimate the subsequent uncertainties in \textit{ab initio} modeling of light and medium mass nuclei which exploit chiral effective field theory. As an example of the this, we present the first complete predictions with uncertainty quantification of peripheral phase shifts of elastic proton-neutron scattering. △ Less

Submitted 2 October, 2014; originally announced October 2014.

Statistical uncertainties of a chiral interaction at next-to-next-to leading order

Authors: A. Ekström, B. D. Carlsson, K. A. Wendt, C. Forssén, M. Hjorth-Jensen, R. Machleidt, S. M. Wild

Abstract: We have quantified the statistical uncertainties of the low-energy coupling-constants (LECs) of an optimized nucleon-nucleon (NN) interaction from chiral effective field theory ($χ$EFT) at next-to-next-to-leading order (NNLO). In addition, we have propagated the impact of the uncertainties of the LECs to two-nucleon scattering phase shifts, effective range parameters, and deuteron observables. We have quantified the statistical uncertainties of the low-energy coupling-constants (LECs) of an optimized nucleon-nucleon (NN) interaction from chiral effective field theory ($χ$EFT) at next-to-next-to-leading order (NNLO). In addition, we have propagated the impact of the uncertainties of the LECs to two-nucleon scattering phase shifts, effective range parameters, and deuteron observables. △ Less

Submitted 30 June, 2014; v1 submitted 26 June, 2014; originally announced June 2014.

Comments: Submitted to JPhysG Focus issue "Enhancing the interaction between nuclear experiment and theory through information and statistics" (ISNET) (updated author list)

Journal ref: J. Phys. G: Nucl. Part. Phys. 42 034003 (2015)

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

Authors: A. Ekström, G. R. Jansen, K. A. Wendt, G. Hagen, T. Papenbrock, S. Bacca, B. Carlsson, D. Gazit

Abstract: 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 c… ▽ More 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. △ Less

Submitted 18 June, 2014; originally announced June 2014.

Journal ref: Phys. Rev. Lett. 113, 262504 (2014)

Coupled-cluster calculations of nucleonic matter

Authors: G. Hagen, T. Papenbrock, A. Ekström, K. A. Wendt, G. Baardsen, S. Gandolfi, M. Hjorth-Jensen, C. J. Horowitz

Abstract: 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 be… ▽ More 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. △ Less

Submitted 12 November, 2013; originally announced November 2013.

Journal ref: Phys. Rev. C 89, 014319 (2014)

Quadrupole collectivity in neutron-deficient Sn nuclei: \nuc{104}{Sn} and the role of proton excitations

Authors: V. M. Bader, A. Gade, D. Weisshaar, B. A. Brown, T. Baugher, D. Bazin, J. S. Berryman, A. Ekstrom, M. Hjorth-Jensen, S. R. Stroberg, W. B. Walters, K. Wimmer, R. Winkler

Abstract: We report on the experimental study of quadrupole collectivity in the neutron-deficient nucleus \nuc{104}{Sn} using intermediate-energy Coulomb excitation. The $B(E2; 0^+_1 \rightarrow 2^+_1)$ value for the excitation of the first $2^+$ state in \nuc{104}{Sn} has been measured to be $0.180(37)~e^2$b$^2$ relative to the well-known $B(E2)$ value of \nuc{102}{Cd}. This result disagrees by more than o… ▽ More We report on the experimental study of quadrupole collectivity in the neutron-deficient nucleus \nuc{104}{Sn} using intermediate-energy Coulomb excitation. The $B(E2; 0^+_1 \rightarrow 2^+_1)$ value for the excitation of the first $2^+$ state in \nuc{104}{Sn} has been measured to be $0.180(37)~e^2$b$^2$ relative to the well-known $B(E2)$ value of \nuc{102}{Cd}. This result disagrees by more than one sigma with a recently published measurement \cite{Gua13}. Our result indicates that the most modern many-body calculations remain unable to describe the enhanced collectivity below mid-shell in Sn approaching $N=Z=50$. We attribute the enhanced collectivity to proton particle-hole configurations beyond the necessarily limited shell-model spaces and suggest the asymmetry of the $B(E2)$-value trend around mid-shell to originate from enhanced proton excitations across $Z=50$ as $N=Z$ is approached. △ Less

Submitted 22 October, 2013; originally announced October 2013.

Comments: Accepted for publication as rapid communication in Physical Review C

Journal ref: Phys. Rev. C 88, 051301(R) (2013)

Coupled-cluster studies of infinite nuclear matter

Authors: G. Baardsen, A. Ekström, G. Hagen, M. Hjorth-Jensen

Abstract: 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 expan… ▽ More 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. △ Less

Submitted 15 November, 2013; v1 submitted 24 June, 2013; originally announced June 2013.

Comments: 18 pages, 14 figures

Journal ref: Phys. Rev. C 88, 054312 (2013)

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

Authors: A. Ekström, G. Baardsen, C. Forssén, G. Hagen, M. Hjorth-Jensen, G. R. Jansen, R. Machleidt, W. Nazarewicz, T. Papenbrock, J. Sarich, S. M. Wild

Abstract: 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 us… ▽ More 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. △ Less

Submitted 19 March, 2013; originally announced March 2013.

Comments: 6 pages, 4 figures

Journal ref: Phys. Rev. Lett. 110, 192502 (2013)

Universal properties of infrared oscillator basis extrapolations

Authors: S. N. More, A. Ekstrom, R. J. Furnstahl, G. Hagen, T. Papenbrock

Abstract: 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 spa… ▽ More 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. △ Less

Submitted 24 April, 2013; v1 submitted 15 February, 2013; originally announced February 2013.

Comments: 16 pages, 21 figures, published version

Journal ref: Phys. Rev. C 87, 044326 (2013)

Coulomb excitation of 73Ga

Authors: J. Diriken, I. Stefanescu, D. Balabanski, N. Blasi, A. Blazhev, N. Bree, J. Cederkäll, T. E. Cocolios, T. Davinson, J. Eberth, A. Ekström, D. V. Fedorov, V. N. Fedosseev, L. M. Fraile, S. Franchoo, G. Georgiev, K. Gladnishki, M. Huyse, O. V. Ivanov, V. S. Ivanov, J. Iwanicki, J. Jolie, T. Konstantinopoulos, Th. Kröll, R. Krücken , et al. (18 additional authors not shown)

Abstract: The B(E2; Ii -> If) values for transitions in 71Ga and 73Ga were deduced from a Coulomb excitation experiment at the safe energy of 2.95 MeV/nucleon using post-accelerated beams of 71,73Ga at the REX-ISOLDE on-line isotope mass separator facility. The emitted gamma rays were detected by the MINIBALL-detector array and B(E2; Ii->If) values were obtained from the yields normalized to the known stren… ▽ More The B(E2; Ii -> If) values for transitions in 71Ga and 73Ga were deduced from a Coulomb excitation experiment at the safe energy of 2.95 MeV/nucleon using post-accelerated beams of 71,73Ga at the REX-ISOLDE on-line isotope mass separator facility. The emitted gamma rays were detected by the MINIBALL-detector array and B(E2; Ii->If) values were obtained from the yields normalized to the known strength of the 2+ -> 0+ transition in the 120Sn target. The comparison of these new results with the data of less neutron-rich gallium isotopes shows a shift of the E2 collectivity towards lower excitation energy when adding neutrons beyond N = 40. This supports conclusions from previous studies of the gallium isotopes which indicated a structural change in this isotopical chain between N = 40 and N = 42. Combined with recent measurements from collinear laser spectroscopy showing a 1/2- spin and parity for the ground state, the extracted results revealed evidence for a 1/2-; 3/2- doublet near the ground state in 73 31Ga42 differing by at most 0.8 keV in energy. △ Less

Submitted 25 November, 2010; v1 submitted 13 October, 2010; originally announced October 2010.

Journal ref: Phys.Rev.C82:064309,2010