• The determination of nucleus-nucleus potentials is important not only to describe the properties of the colliding system, but also to extract nuclear-structure information and for modelling nuclear reactions for astrophysics. We present the first determination of double-folding potentials based on chiral effective field theory at leading, next-to-leading, and next-to-next-to-leading order. To this end, we construct new soft local chiral effective field theory interactions. We benchmark this approach in the $^{16}$O-$^{16}$O system, and present results for cross sections computed for elastic scattering up to 700 MeV in energy, as well as for the astrophysical $S$-factor of the fusion reaction.
  • We study systems of bosons and fermions in finite periodic boxes and show how the existence and properties of few-body resonances can be extracted from studying the volume dependence of the calculated energy spectra. Using a plane-wave-based discrete variable representation to conveniently implement periodic boundary conditions, we establish that avoided level crossings occur in the spectra of up to four particles and can be linked to the existence of multi-body resonances. To benchmark our method we use two-body calculations, where resonance properties can be determined with other methods, as well as a three-boson model interaction known to generate a three-boson resonance state. Finding good agreement for these cases, we then predict three-body and four-body resonances for models using a shifted Gaussian potential. Our results establish few-body finite-volume calculations as a new tool to study few-body resonances. In particular, the approach can be used to study few-neutron systems, where such states have been conjectured to exist.
  • We discuss the uncertainties in constraining low-energy constants of chiral effective field theory from $^3$H $\beta$ decay. The half-life is very precisely known, so that the Gamow-Teller matrix element has been used to fit the coupling $c_D$ of the axial-vector current to a short-range two-nucleon pair. Because the same coupling also describes the leading one-pion-exchange three-nucleon force, this in principle provides a very constraining fit, uncorrelated with the $^3$H binding energy fit used to constrain another low-energy coupling in three-nucleon forces. However, so far such $^3$H half-life fits have only been performed at a fixed cutoff value. We show that the cutoff dependence due to the regulators in the axial-vector two-body current can significantly affect the Gamow-Teller matrix elements and consequently also the extracted values for the $c_D$ coupling constant. The degree of the cutoff dependence is correlated with the softness of the employed NN interaction. As a result, present three-nucleon forces based on a fit to $^3$H $\beta$ decay underestimate the uncertainty in $c_D$. We explore a range of $c_D$ values that is compatible within cutoff variation with the experimental $^3$H half-life and estimate the resulting uncertainties for many-body systems by performing calculations of symmetric nuclear matter.
  • We report accurate quantum Monte Carlo calculations of nuclei up to $A=16$ based on local chiral two- and three-nucleon interactions up to next-to-next-to-leading order. We examine the theoretical uncertainties associated with the chiral expansion and the cutoff in the theory, as well as the associated operator choices in the three-nucleon interactions. While in light nuclei the cutoff variation and systematic uncertainties are rather small, in $^{16}$O these can be significant for large coordinate-space cutoffs. Overall, we show that chiral interactions constructed to reproduce properties of very light systems and nucleon-nucleon scattering give an excellent description of binding energies, charge radii, and form factors for all these nuclei, including open-shell systems in $A=6$ and 12.
  • Quantum Monte Carlo methods have recently been employed to study properties of nuclei and infinite matter using local chiral effective field theory interactions. In this work, we present a detailed description of the auxiliary field diffusion Monte Carlo algorithm for nuclei in combination with local chiral two- and three-nucleon interactions up to next-to-next-to-leading order. We show results for the binding energy, charge radius, charge form factor, and Coulomb sum rule in nuclei with $3\le A\le16$. Particular attention is devoted to the effect of different operator structures in the three-body force for different cutoffs. The outcomes suggest that local chiral interactions fit to few-body observables give a very good description of the ground-state properties of nuclei up to $^{16}$O, with the exception of one fit for the softer cutoff which predicts overbinding in larger nuclei.
  • A precision mass investigation of the neutron-rich titanium isotopes $^{51-55}$Ti was performed at TRIUMF's Ion Trap for Atomic and Nuclear science (TITAN). The range of the measurements covers the $N=32$ shell closure and the overall uncertainties of the $^{52-55}$Ti mass values were significantly reduced. Our results confirm the existence of a weak shell effect at $N=32$, establishing the abrupt onset of this shell closure. Our data were compared with state-of-the-art \textit{ab-initio} shell model calculations which, despite very successfully describing where the $N=32$ shell gap is strong, overpredict its strength and extent in titanium and heavier isotones. These measurements also represent the first scientific results of TITAN using the newly commissioned Multiple-Reflection Time-of-Flight Mass Spectrometer (MR-TOF-MS), substantiated by independent measurements from TITAN's Penning trap mass spectrometer.
  • We calculate the frequencies of collective modes of neutrons, protons and electrons in the outer core of neutron stars. The neutrons and protons are treated in a hydrodynamic approximation and the electrons are regarded as collisionless. The coupling of the nucleons to the electrons leads to Landau damping of the collective modes and to significant dispersion of the low-lying modes. We investigate the sensitivity of the mode frequencies to the strength of entrainment between neutrons and protons, which is not well characterized. The contribution of collective modes to the thermal conductivity is evaluated.
  • We link the structure of nuclei around $^{100}$Sn, the heaviest doubly magic nucleus with equal neutron and proton numbers ($N=Z=50$), to nucleon-nucleon ($NN$) and three-nucleon ($NNN$) forces constrained by data of few-nucleon systems. Our results indicate that $^{100}$Sn is doubly magic, and we predict its quadrupole collectivity. We present precise computations of $^{101}$Sn based on three-particle--two-hole excitations of $^{100}$Sn, and reproduce the small splitting between the lowest $J^\pi=7/2^+$ and $5/2^+$ states. Our results are consistent with the sparse available data.
  • Chiral effective field theory is a framework to derive systematic nuclear interactions. It is based on the symmetries of quantum chromodynamics and includes long-range pion physics explicitly, while shorter-range physics is expanded in a general operator basis. The number of low-energy couplings at a particular order in the expansion can be reduced by exploiting the fact that nucleons are fermions and therefore obey the Pauli exclusion principle. The antisymmetry permits the selection of a subset of the allowed contact operators at a given order. When local regulators are used for these short-range interactions, however, this "Fierz rearrangement freedom" is violated. In this paper, we investigate the impact of this violation at leading order (LO) in the chiral expansion. We construct LO and next-to-leading order (NLO) potentials for all possible LO-operator pairs and study their reproduction of phase shifts, the ${}^4$He ground-state energy, and the neutron-matter energy at different densities. We demonstrate that the Fierz rearrangement freedom is partially restored at NLO where subleading contact interactions enter. We also discuss implications for local chiral three-nucleon interactions.
  • Background: Odd-odd nuclei, around doubly closed shells, have been extensively used to study proton-neutron interactions. However, the evolution of these interactions as a function of the binding energy, ultimately when nuclei become unbound, is poorly known. The $^{26}$F nucleus, composed of a deeply bound $\pi0d\_{5/2}$ proton and an unbound $\nu0d\_{3/2}$ neutron on top of an $^{24}$O core, is particularly adapted for this purpose. The coupling of this proton and neutron results in a $J^{\pi} = 1^{+}\_1 - 4^{+}\_1$ multiplet, whose energies must be determined to study the influence of the proximity of the continuum on the corresponding proton-neutron interaction. The $J^{\pi} = 1^{+}\_1, 2^{+}\_1,4^{+}\_1$ bound states have been determined, and only a clear identification of the $J^{\pi} =3^{+}\_1$ is missing.Purpose: We wish to complete the study of the $J^{\pi} = 1^{+}\_1 - 4^{+}\_1$ multiplet in $^{26}$F, by studying the energy and width of the $J^{\pi} =3^{+}\_1$ unbound state. The method was firstly validated by the study of unbound states in $^{25}$F, for which resonances were already observed in a previous experiment.Method: Radioactive beams of $^{26}$Ne and $^{27}$Ne, produced at about $440A$\,MeV by the FRagment Separator at the GSI facility, were used to populate unbound states in $^{25}$F and $^{26}$F via one-proton knockout reactions on a CH$\_2$ target, located at the object focal point of the R$^3$B/LAND setup. The detection of emitted $\gamma$-rays and neutrons, added to the reconstruction of the momentum vector of the $A-1$ nuclei, allowed the determination of the energy of three unbound states in $^{25}$F and two in $^{26}$F. Results: Based on its width and decay properties, the first unbound state in $^{25}$F is proposed to be a $J^{\pi} = 1/2^-$ arising from a $p\_{1/2}$ proton-hole state. In $^{26}$F, the first resonance at 323(33)~keV is proposed to be the $J^{\pi} =3^{+}\_1$ member of the $J^{\pi} = 1^{+}\_1 - 4^{+}\_1$ multiplet. Energies of observed states in $^{25,26}$F have been compared to calculations using the independent-particle shell model, a phenomenological shell-model, and the ab initio valence-space in-medium similarity renormalization group method.Conclusions: The deduced effective proton-neutron interaction is weakened by about 30-40\% in comparison to the models, pointing to the need of implementing the role of the continuum in theoretical descriptions, or to a wrong determination of the atomic mass of $^{26}$F.
  • We perform a comprehensive Weinberg eigenvalue analysis of a representative set of modern nucleon-nucleon interactions derived within chiral effective field theory. Our set contains local, semilocal, and nonlocal potentials, developed by Gezerlis, Tews et al. (2013); Epelbaum, Krebs, and Mei{\ss}ner (2015); and Entem, Machleidt, and Nosyk (2017) as well as Carlsson, Ekstr\"om et al. (2016), respectively. The attractive eigenvalues show a very similar behavior for all investigated interactions, whereas the magnitudes of the repulsive eigenvalues sensitively depend on the details of the regularization scheme of the short- and long-range parts of the interactions. We demonstrate that a direct comparison of numerical cutoff values of different interactions is in general misleading due to the different analytic form of regulators; for example, a cutoff value of $R=0.8$ fm for the semilocal interactions corresponds to about $R=1.2$ fm for the local interactions. Our detailed comparison of Weinberg eigenvalues provides various insights into idiosyncrasies of chiral potentials for different orders and partial waves. This shows that Weinberg eigenvalues could be used as a helpful monitoring scheme when constructing new interactions.
  • Local chiral effective field theory interactions have recently been developed and used in the context of quantum Monte Carlo few- and many-body methods for nuclear physics. In this work, we go over detailed features of local chiral nucleon-nucleon interactions and examine their effect on properties of the deuteron, paying special attention to the perturbativeness of the expansion. We then turn to three-nucleon interactions, focusing on operator ambiguities and their interplay with regulator effects. We then discuss the nuclear Green's function Monte Carlo method, going over both wave-function correlations and approximations for the two- and three-body propagators. Following this, we present a range of results on light nuclei: Binding energies and distribution functions are contrasted and compared, starting from several different microscopic interactions.
  • We present quantum Monte Carlo calculations of few-neutron systems confined in external potentials based on local chiral interactions at next-to-next-to-leading order in chiral effective field theory. The energy and radial densities for these systems are calculated in different external Woods-Saxon potentials. We assume that their extrapolation to zero external-potential depth provides a quantitative estimate of three- and four-neutron resonances. The validity of this assumption is demonstrated by benchmarking with an exact diagonalization in the two-body case. We find that the extrapolated trineutron resonance, as well as the energy for shallow well depths, is lower than the tetraneutron resonance energy. This suggests that a three-neutron resonance exists below a four-neutron resonance in nature and is potentially measurable. To confirm that the relative ordering of three- and four-neutron resonances is not an artifact of the external confinement, we test that the odd-even staggering in the helium isotopic chain is reproduced within this approach. Finally, we discuss similarities between our results and ultracold Fermi gases.
  • We explore the impact of nuclear matter saturation on the properties and systematics of finite nuclei across the nuclear chart. Using the ab initio in-medium similarity renormalization group (IM-SRG), we study ground-state energies and charge radii of closed-shell nuclei from $^4$He to $^{78}$Ni, based on a set of low-resolution two- and three-nucleon interactions that predict realistic saturation properties. We first investigate in detail the convergence properties of these Hamiltonians with respect to model-space truncations for both two- and three-body interactions. We find one particular interaction that reproduces well the ground-state energies of all closed-shell nuclei studied. As expected from their saturation points relative to this interaction, the other Hamiltonians underbind nuclei, but lead to a remarkably similar systematics of ground-state energies. Extending our calculations to complete isotopic chains in the $sd$ and $pf$ shells with the valence-space IM-SRG, the same interaction reproduces not only experimental ground states but two-neutron-separation energies and first excited $2^+$ states. We also calculate radii with the valence-space IM-SRG for the first time. Since this particular interaction saturates at too high density, charge radii are still too small compared with experiment. Except for this underprediction, the radii systematics is, however, well reproduced. Our results highlight the importance of nuclear matter as a theoretical benchmark for the development of next-generation chiral interactions.
  • Neutron matter is an ideal laboratory for nuclear interactions derived from chiral effective field theory since all contributions are predicted up to next-to-next-to-next-to-leading order (N$^3$LO) in the chiral expansion. By making use of recent advances in the partial-wave decomposition of three- nucleon (3N) forces, we include for the first time N$^3$LO 3N interactions in many-body perturbation theory (MBPT) up to third order and in self-consistent Green's function theory (SCGF). Using these two complementary many-body frameworks we provide improved predictions for the equation of state of neutron matter at zero temperature and also analyze systematically the many-body convergence for different chiral EFT interactions. Furthermore, we present an extension of the normal-ordering framework to finite temperatures. These developments open the way to improved calculations of neutron-rich matter including estimates of theoretical uncertainties for astrophysical applications.
  • Ab initio calculations provide direct access to the properties of pure neutron systems that are challenging to study experimentally. In addition to their importance for fundamental physics, their properties are required as input for effective field theories of the strong interaction. In this work, we perform auxiliary-field diffusion Monte Carlo calculations of the ground and first excited state of two neutrons in a finite box, considering a simple contact potential as well as chiral effective field theory interactions. We compare the results against exact diagonalizations and present a detailed analysis of the finite-volume effects, whose understanding is crucial for determining observables from the calculated energies. Using the L\"uscher formula, we extract the low-energy S-wave scattering parameters from ground- and excited-state energies for different box sizes.
  • We generalize our virial approach to study the neutral current neutrino response of nuclear matter at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for neutrino-nucleon scattering rates and the density S_V and spin S_A responses. We find S_A is significantly reduced from one even at low densities. We provide a simple fit S_A^f(n,T,Y_p) of the axial response as a function of density n, temperature T and proton fraction Y_p. This fit reproduces our model independent virial results at low densities and reproduces the Burrows and Sawyer random phase approximation (RPA) results at high densities. Our fit can be incorporated into supernova simulations in a straight forward manner. Preliminary one dimensional supernova simulations suggest that the reduction in the axial response may enhance neutrino heating rates in the gain region during the accretion phase of a core-collapse supernovae.
  • We present solutions of the BCS gap equation in the channels ${}^1S_0$ and ${}^3P_2-{}^3F_2$ in neutron matter based on nuclear interactions derived within chiral effective field theory (EFT). Our studies are based on a representative set of nonlocal nucleon-nucleon (NN) plus three-nucleon (3N) interactions up to next-to-next-to-next-to-leading order (N$^3$LO) as well as local and semilocal chiral NN interactions up to N$^2$LO and N$^4$LO, respectively. In particular, we investigate for the first time the impact of subleading 3N forces at N$^3$LO on pairing gaps and also derive uncertainty estimates by taking into account results for pairing gaps at different orders in the chiral expansion. Finally, we discuss different methods for obtaining self-consistent solutions of the gap equation. Besides the widely-used quasi-linear method by Khodel et al. we demonstrate that the modified Broyden method is well applicable and exhibits a robust convergence behavior. In contrast to Khodel's method it is based on a direct iteration of the gap equation without imposing an auxiliary potential and is straightforward to implement.
  • The calcium isotopes have emerged as an important testing ground for new microscopically derived shell-model interactions, and a great deal of focus has been directed toward this region. We investigate the relative spectroscopic strengths associated with $1f_{7/2}$ neutron hole states in $^{47, 49}$Ca following one-neutron knockout reactions from $^{48,50}$Ca. The observed reduction of strength populating the lowest 7/2$^{-}_{1}$ state in $^{49}$Ca, as compared to $^{47}$Ca, is consistent with the description given by shell-model calculations based on two- and three-nucleon forces in the neutron $pf$ model space, implying a fragmentation of the $l$=3 strength to higher-lying states. The experimental result is inconsistent with both the GXPF1 interaction routinely used in this region of the nuclear chart and with microscopic calculations in an extended model space including the $\nu1g_{9/2}$ orbital.
  • We present a nucleus-dependent valence-space approach for calculating ground and excited states of nuclei, which generalizes the shell-model in-medium similarity renormalization group to an ensemble reference with fractionally filled orbitals. Because the ensemble is used only as a reference, and not to represent physical states, no symmetry restoration is required. This allows us to capture 3N forces among valence nucleons with a valence-space Hamiltonian specifically targeted to each nucleus of interest. Predicted ground-state energies from carbon through nickel agree with results of other large-space ab initio methods, generally to the 1\% level. In addition, we show that this new approach is required in order to obtain convergence for nuclei in the upper $p$ and $sd$ shells. Finally, we address the $1^+$/$3^+$ ground-state inversion problem in $^{22}\text{Na}$ and $^{46}\text{V}$. This approach extends the reach of ab initio nuclear structure calculations to essentially all light- and medium-mass nuclei.
  • We present ab initio predictions for ground and excited states of doubly open-shell fluorine and neon isotopes based on chiral two- and three-nucleon interactions. We use the in-medium similarity renormalization group, in both flow-equation and Magnus formulations, to derive mass-dependent sd valence-space Hamiltonians. The experimental ground-state energies are reproduced through neutron number N=14, beyond which a new targeted normal-ordering procedure improves agreement with data and large-scale multi-reference calculations. For spectroscopy, we focus on neutron-rich 23-26F and 24-26Ne isotopes near N=14,16 magic numbers. In all cases we find an agreement with experiment competitive with established phenomenology. Moreover, yrast states are well described in 20Ne and 24Mg, providing an ab initio description of deformation in the medium-mass region.
  • Despite being a complex many-body system, the atomic nucleus exhibits simple structures for certain "magic" numbers of protons and neutrons. The calcium chain in particular is both unique and puzzling: evidence of doubly-magic features are known in 40,48Ca, and recently suggested in two radioactive isotopes, 52,54Ca. Although many properties of experimentally known Ca isotopes have been successfully described by nuclear theory, it is still a challenge to predict their charge radii evolution. Here we present the first measurements of the charge radii of 49,51,52Ca, obtained from laser spectroscopy experiments at ISOLDE, CERN. The experimental results are complemented by state-of-the-art theoretical calculations. The large and unexpected increase of the size of the neutron-rich calcium isotopes beyond N = 28 challenges the doubly-magic nature of 52Ca and opens new intriguing questions on the evolution of nuclear sizes away from stability, which are of importance for our understanding of neutron-rich atomic nuclei.
  • We present quantum Monte Carlo calculations of light nuclei, neutron-$\alpha$ scattering, and neutron matter using local two- and three-nucleon ($3N$) interactions derived from chiral effective field theory up to next-to-next-to-leading order (N$^2$LO). The two undetermined $3N$ low-energy couplings are fit to the $^4$He binding energy and, for the first time, to the spin-orbit splitting in the neutron-$\alpha$ $P$-wave phase shifts. Furthermore, we investigate different choices of local $3N$-operator structures and find that chiral interactions at N$^2$LO are able to simultaneously reproduce the properties of $A=3,4,5$ systems and of neutron matter, in contrast to commonly used phenomenological $3N$ interactions.
  • Chiral effective field theory (EFT) enables a systematic description of low-energy hadronic interactions with controlled theoretical uncertainties. For strongly interacting systems, quantum Monte Carlo (QMC) methods provide some of the most accurate solutions, but they require as input local potentials. We have recently constructed local chiral nucleon-nucleon (NN) interactions up to next-to-next-to-leading order (N$^2$LO). Chiral EFT naturally predicts consistent many-body forces. In this paper, we consider the leading chiral three-nucleon (3N) interactions in local form. These are included in auxiliary field diffusion Monte Carlo (AFDMC) simulations. We present results for the equation of state of neutron matter and for the energies and radii of neutron drops. In particular, we study the regulator dependence at the Hartree-Fock level and in AFDMC and find that present local regulators lead to less repulsion from 3N forces compared to the usual nonlocal regulators.
  • We present a comprehensive review of the In-Medium Similarity Renormalization Group (IM-SRG), a novel ab inito method for nuclei. The IM-SRG employs a continuous unitary transformation of the many-body Hamiltonian to decouple the ground state from all excitations, thereby solving the many-body problem. Starting from a pedagogical introduction of the underlying concepts, the IM-SRG flow equations are developed for systems with and without explicit spherical symmetry. We study different IM-SRG generators that achieve the desired decoupling, and how they affect the details of the IM-SRG flow. Based on calculations of closed-shell nuclei, we assess possible truncations for closing the system of flow equations in practical applications, as well as choices of the reference state. We discuss the issue of center-of-mass factorization and demonstrate that the IM-SRG ground-state wave function exhibits an approximate decoupling of intrinsic and center-of-mass degrees of freedom, similar to Coupled Cluster (CC) wave functions. To put the IM-SRG in context with other many-body methods, in particular many-body perturbation theory and non-perturbative approaches like CC, a detailed perturbative analysis of the IM-SRG flow equations is carried out. We conclude with a discussion of ongoing developments, including IM-SRG calculations with three-nucleon forces, the multi-reference IM-SRG for open-shell nuclei, first non-perturbative derivations of shell- model interactions, and the consistent evolution of operators in the IM-SRG. We dedicate this review to the memory of Gerry Brown, one of the pioneers of many-body calculations of nuclei.