
The determination of nucleusnucleus potentials is important not only to
describe the properties of the colliding system, but also to extract
nuclearstructure information and for modelling nuclear reactions for
astrophysics. We present the first determination of doublefolding potentials
based on chiral effective field theory at leading, nexttoleading, and
nexttonexttoleading 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 fewbody resonances can be extracted from
studying the volume dependence of the calculated energy spectra. Using a
planewavebased 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
multibody resonances. To benchmark our method we use twobody calculations,
where resonance properties can be determined with other methods, as well as a
threeboson model interaction known to generate a threeboson resonance state.
Finding good agreement for these cases, we then predict threebody and
fourbody resonances for models using a shifted Gaussian potential. Our results
establish fewbody finitevolume calculations as a new tool to study fewbody
resonances. In particular, the approach can be used to study fewneutron
systems, where such states have been conjectured to exist.

We discuss the uncertainties in constraining lowenergy constants of chiral
effective field theory from $^3$H $\beta$ decay. The halflife is very
precisely known, so that the GamowTeller matrix element has been used to fit
the coupling $c_D$ of the axialvector current to a shortrange twonucleon
pair. Because the same coupling also describes the leading onepionexchange
threenucleon force, this in principle provides a very constraining fit,
uncorrelated with the $^3$H binding energy fit used to constrain another
lowenergy coupling in threenucleon forces. However, so far such $^3$H
halflife fits have only been performed at a fixed cutoff value. We show that
the cutoff dependence due to the regulators in the axialvector twobody
current can significantly affect the GamowTeller 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 threenucleon 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 halflife and estimate the resulting uncertainties for
manybody 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 threenucleon interactions up to
nexttonexttoleading 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 threenucleon interactions. While in
light nuclei the cutoff variation and systematic uncertainties are rather
small, in $^{16}$O these can be significant for large coordinatespace cutoffs.
Overall, we show that chiral interactions constructed to reproduce properties
of very light systems and nucleonnucleon scattering give an excellent
description of binding energies, charge radii, and form factors for all these
nuclei, including openshell 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 threenucleon interactions up to nexttonexttoleading 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 threebody force for
different cutoffs. The outcomes suggest that local chiral interactions fit to
fewbody observables give a very good description of the groundstate
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 neutronrich titanium isotopes
$^{5155}$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 $^{5255}$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 stateoftheart \textit{abinitio} 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 MultipleReflection TimeofFlight Mass Spectrometer
(MRTOFMS), 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 lowlying
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 nucleonnucleon
($NN$) and threenucleon ($NNN$) forces constrained by data of fewnucleon
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 threeparticletwohole 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 longrange pion physics explicitly, while shorterrange physics is
expanded in a general operator basis. The number of lowenergy 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 shortrange
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 nexttoleading order (NLO) potentials
for all possible LOoperator pairs and study their reproduction of phase
shifts, the ${}^4$He groundstate energy, and the neutronmatter 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 threenucleon interactions.

Background: Oddodd nuclei, around doubly closed shells, have been
extensively used to study protonneutron 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 protonneutron 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 oneproton 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 $A1$ 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}$ protonhole 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 independentparticle shell model, a
phenomenological shellmodel, and the ab initio valencespace inmedium
similarity renormalization group method.Conclusions: The deduced effective
protonneutron interaction is weakened by about 3040\% 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 nucleonnucleon 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 longrange 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 manybody methods for
nuclear physics. In this work, we go over detailed features of local chiral
nucleonnucleon interactions and examine their effect on properties of the
deuteron, paying special attention to the perturbativeness of the expansion. We
then turn to threenucleon 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 wavefunction correlations and
approximations for the two and threebody 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 fewneutron systems confined
in external potentials based on local chiral interactions at
nexttonexttoleading order in chiral effective field theory. The energy and
radial densities for these systems are calculated in different external
WoodsSaxon potentials. We assume that their extrapolation to zero
externalpotential depth provides a quantitative estimate of three and
fourneutron resonances. The validity of this assumption is demonstrated by
benchmarking with an exact diagonalization in the twobody 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
threeneutron resonance exists below a fourneutron resonance in nature and is
potentially measurable. To confirm that the relative ordering of three and
fourneutron resonances is not an artifact of the external confinement, we test
that the oddeven 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
inmedium similarity renormalization group (IMSRG), we study groundstate
energies and charge radii of closedshell nuclei from $^4$He to $^{78}$Ni,
based on a set of lowresolution two and threenucleon interactions that
predict realistic saturation properties. We first investigate in detail the
convergence properties of these Hamiltonians with respect to modelspace
truncations for both two and threebody interactions. We find one particular
interaction that reproduces well the groundstate energies of all closedshell
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 groundstate energies. Extending our calculations to
complete isotopic chains in the $sd$ and $pf$ shells with the valencespace
IMSRG, the same interaction reproduces not only experimental ground states but
twoneutronseparation energies and first excited $2^+$ states. We also
calculate radii with the valencespace IMSRG 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
nextgeneration chiral interactions.

Neutron matter is an ideal laboratory for nuclear interactions derived from
chiral effective field theory since all contributions are predicted up to
nexttonexttonexttoleading order (N$^3$LO) in the chiral expansion. By
making use of recent advances in the partialwave decomposition of three
nucleon (3N) forces, we include for the first time N$^3$LO 3N interactions in
manybody perturbation theory (MBPT) up to third order and in selfconsistent
Green's function theory (SCGF). Using these two complementary manybody
frameworks we provide improved predictions for the equation of state of neutron
matter at zero temperature and also analyze systematically the manybody
convergence for different chiral EFT interactions. Furthermore, we present an
extension of the normalordering framework to finite temperatures. These
developments open the way to improved calculations of neutronrich 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 auxiliaryfield 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 finitevolume effects, whose understanding is crucial for
determining observables from the calculated energies. Using the L\"uscher
formula, we extract the lowenergy Swave scattering parameters from ground
and excitedstate 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 longwavelength limit, the
virial expansion makes modelindependent predictions for neutrinonucleon
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 corecollapse 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 nucleonnucleon (NN) plus threenucleon (3N)
interactions up to nexttonexttonexttoleading 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
selfconsistent solutions of the gap equation. Besides the widelyused
quasilinear 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 shellmodel 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 oneneutron 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
shellmodel calculations based on two and threenucleon forces in the neutron
$pf$ model space, implying a fragmentation of the $l$=3 strength to
higherlying 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 nucleusdependent valencespace approach for calculating ground
and excited states of nuclei, which generalizes the shellmodel inmedium
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 valencespace Hamiltonian
specifically targeted to each nucleus of interest. Predicted groundstate
energies from carbon through nickel agree with results of other largespace 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^+$ groundstate 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
mediummass nuclei.

We present ab initio predictions for ground and excited states of doubly
openshell fluorine and neon isotopes based on chiral two and threenucleon
interactions. We use the inmedium similarity renormalization group, in both
flowequation and Magnus formulations, to derive massdependent sd
valencespace Hamiltonians. The experimental groundstate energies are
reproduced through neutron number N=14, beyond which a new targeted
normalordering procedure improves agreement with data and largescale
multireference calculations. For spectroscopy, we focus on neutronrich 2326F
and 2426Ne 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 mediummass region.

Despite being a complex manybody 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 doublymagic
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 stateoftheart
theoretical calculations. The large and unexpected increase of the size of the
neutronrich calcium isotopes beyond N = 28 challenges the doublymagic 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
neutronrich atomic nuclei.

We present quantum Monte Carlo calculations of light nuclei, neutron$\alpha$
scattering, and neutron matter using local two and threenucleon ($3N$)
interactions derived from chiral effective field theory up to
nexttonexttoleading order (N$^2$LO). The two undetermined $3N$ lowenergy
couplings are fit to the $^4$He binding energy and, for the first time, to the
spinorbit 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
lowenergy 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 nucleonnucleon (NN) interactions up to
nexttonexttoleading order (N$^2$LO). Chiral EFT naturally predicts
consistent manybody forces. In this paper, we consider the leading chiral
threenucleon (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 HartreeFock
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 InMedium Similarity Renormalization
Group (IMSRG), a novel ab inito method for nuclei. The IMSRG employs a
continuous unitary transformation of the manybody Hamiltonian to decouple the
ground state from all excitations, thereby solving the manybody problem.
Starting from a pedagogical introduction of the underlying concepts, the IMSRG
flow equations are developed for systems with and without explicit spherical
symmetry. We study different IMSRG generators that achieve the desired
decoupling, and how they affect the details of the IMSRG flow. Based on
calculations of closedshell 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 centerofmass factorization and
demonstrate that the IMSRG groundstate wave function exhibits an approximate
decoupling of intrinsic and centerofmass degrees of freedom, similar to
Coupled Cluster (CC) wave functions. To put the IMSRG in context with other
manybody methods, in particular manybody perturbation theory and
nonperturbative approaches like CC, a detailed perturbative analysis of the
IMSRG flow equations is carried out. We conclude with a discussion of ongoing
developments, including IMSRG calculations with threenucleon forces, the
multireference IMSRG for openshell nuclei, first nonperturbative
derivations of shell model interactions, and the consistent evolution of
operators in the IMSRG. We dedicate this review to the memory of Gerry Brown,
one of the pioneers of manybody calculations of nuclei.