With the goal of developing predictive ab-initio capability for light and medium-mass nuclei, twonucleon 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 NNLOsat, 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. Introduction -Interactions from chiral effective field theory (EFT) [1][2][3][4] and modern applications of renormalization group techniques [5][6][7][8] have opened the door for a description of atomic nuclei consistent with the underlying symmetries of quantum chromodynamics, the theory of the strong interaction. Chiral nuclear forces can be constructed systematically from long-range pion physics augmented by contact interactions. Over the past decade, the renaissance of nuclear theory based on realistic nuclear forces and powerful computational methods has pushed the frontier of ab initio calculations from fewbody systems and light nuclei [6, 9, 10] to medium-mass nuclei [11][12][13][14][15][16][17][18][19].

SummaryWhat 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.

We optimize the nucleon-nucleon interaction from chiral effective field theory at next-to-next-to-leading order (NNLO). The resulting new chiral force NNLO(opt) yields χ(2)≈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 NNLO(opt) to study properties of key nuclei and neutron matter, and we demonstrate that many aspects of nuclear structure can be understood in terms of this nucleon-nucleon interaction, without explicitly invoking three-nucleon forces.

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 calcium isotopes have been successfully described by nuclear theory, it is still a challenge to predict the evolution of their charge radii. 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

Chiral effective field theory (χEFT) provides a systematic approach to describe low-energy nuclear forces. Moreover, χEFT 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 χEFT. Finally, we study the effect on other observables by demonstrating forwarderror-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. This is also vital for the regression analysis. We use power-counting arguments to estimate the systematic uncertainty that is inherent to χEFT, 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 χEFT. Furthermore, certain systematic uncertainties in the few-nucleon sector are shown to get substantially magnified in the many-body sector, in particular 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.

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

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.

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 coupledcluster calculations of nuclei up to 132 Sn converge fast for the ground-state energies and radii in feasible model spaces.