
The correlation of the tidal polarizabilities $\Lambda_1$$\Lambda_2$ for
GW170817 is analyzed by combining dense matter equations of state (EOS) that
satisfy nuclear physics constraints and the posteriors of the LIGOVirgo (LV)
gravitational wave data analysis. Two EOS models are explored: a minimal model
(MM) relies on an extrapolation of a parametrized nuclear EOS, and a second
model, labelled CSM, is a general parametrization of the speed of sound that
also allows for phase transitions at high densities. Both models are
constrained by calculations of the neutron matter EOS using chiral effective
field theory (EFT) Hamiltonians up to a given density $n_{\text{tr}}$, and by
requiring that the neutron star maximum mass $M_{\text max} \ge 1.9~M_\odot$.
We find that the CSM model, which subsumes the predictions of MM, is
constrained by the gravitational wave observation from GW170817 if
$n_{\text{tr}}=n_{\text{sat}}$ where $n_{\text{sat}}=0.16$ fm$^{3}$. However,
if the EFT remains valid up to $n_{\text{tr}}=2 n_{\text{sat}}$ we predict the
combined tidal polarizability of the NSs in GW170817 to be $80\le
\tilde{\Lambda}\le 570$ ($280\le \tilde{\Lambda}\le 480$ for the minimal
model), which is compatible but smaller than the range suggested by LV data
analysis. Our analysis also shows that GW170817 requires a neutron star with
$M=1.4M_\odot$ to have a radius $R_{1.4}<13.6$~km.

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.

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.

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.

Regulator functions applied to two and threenucleon forces are a necessary
ingredient in manybody calculations based on chiral effective field theory
interactions. These interactions have been developed recently with a variety of
different cutoff forms, including regulating both the momentum transfer (local)
and the relative momentum (nonlocal). While in principle any regulator that
suppresses high momentum modes can be employed, in practice artifacts are
inevitable in current power counting schemes. Artifacts from particular
regulators may cause significant distortions of the physics or may affect
manybody convergence rates, so understanding their nature is important. Here
we characterize the differences between cutoff effects using uniform matter at
HartreeFock and secondorder in the interaction as a testbed. This provides a
clean laboratory to isolate phasespace effects of various regulators on both
two and threenucleon interactions. We test the normalordering approximation
for threenucleon forces in nuclear matter and find that the relative size of
the residual 3N contributions is sensitive to the employed regularization
scheme.

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 details of the derivation of local chiral effective field theory
interactions to nexttonexttoleading order, and show results for
nucleonnucleon phase shifts and deuteron properties for these potentials. We
then perform systematic auxiliaryfield diffusion Monte Carlo calculations for
neutron matter based on the developed local chiral potentials at different
orders. This includes studies of the effects of the spectralfunction
regularization and of the local regulators. For all orders, we compare the
quantum Monte Carlo results with perturbative manybody calculations and find
excellent agreement for low cutoffs.

Neutron matter presents a unique system in chiral effective field theory
(EFT), because all manybody forces among neutrons are predicted to
nexttonexttonexttoleading order (N3LO). We discuss perturbative and first
Quantum Monte Carlo (QMC) calculations of neutron matter with chiral EFT
interactions and their astrophysical impact for the equation of state and
neutron stars.

We calculate the chiral condensate in neutron matter at zero temperature
based on nuclear forces derived within chiral effective field theory. Two,
three and fournucleon interactions are included consistently to
nexttonexttonexttoleading order (N3LO) of the chiral expansion. We find
that the interaction contributions lead to a modest increase of the condensate,
thus impeding the restoration of chiral symmetry in dense matter and making a
chiral phase transition in neutronrich matter unlikely for densities that are
not significantly higher than nuclear saturation density.

The neutronmatter equation of state constrains the properties of many
physical systems over a wide density range and can be studied systematically
using chiral effective field theory (EFT). In chiral EFT, all manybody forces
among neutrons are predicted to nexttonexttonexttoleading order (N3LO).
We present details and additional results of the first complete N3LO
calculation of the neutronmatter energy, which includes the subleading
threenucleon as well as the leading fournucleon forces, and provides
theoretical uncertainties. In addition, we discuss the impact of our results
for astrophysics: for the supernova equation of state, the symmetry energy and
its density derivative, and for the structure of neutron stars. Finally, we
give a first estimate for the size of the N3LO manybody contributions to the
energy of symmetric nuclear matter, which shows that their inclusion will be
important in nuclear structure calculations.

We present the first quantum Monte Carlo (QMC) calculations with chiral
effective field theory (EFT) interactions. To achieve this, we remove all
sources of nonlocality, which hamper the inclusion in QMC calculations, in
nuclear forces to nexttonexttoleading order. We perform auxiliaryfield
diffusion Monte Carlo (AFDMC) calculations for the neutron matter energy up to
saturation density based on local leadingorder, nexttoleading order, and
nexttonexttoleading order nucleonnucleon interactions. Our results exhibit
a systematic orderbyorder convergence in chiral EFT and provide
nonperturbative benchmarks with theoretical uncertainties. For the softer
interactions, perturbative calculations are in excellent agreement with the
AFDMC results. This work paves the way for QMC calculations with systematic
chiral EFT interactions for nuclei and nuclear matter, for testing the
perturbativeness of different orders, and allows for matching to lattice QCD
results by varying the pion mass.

Neutron matter presents a unique system for chiral effective field theory
(EFT), because all manybody forces among neutrons are predicted to
nexttonexttonexttoleading order (N3LO). We present the first complete
N3LO calculation of the neutron matter energy. This includes the subleading
threenucleon (3N) forces for the first time and all leading fournucleon (4N)
forces. We find relatively large contributions from N3LO 3N forces. Our results
provide constraints for neutronrich matter in astrophysics with controlled
theoretical uncertainties.