We report a quantum Monte Carlo calculation of the equation of state of symmetric nuclear matter using local interactions derived from chiral effective field theory up to next-to-next-to-leading order fit to few-body observables only. The empirical saturation density and energy are well reproduced within statistical and systematic uncertainties. We have also derived the symmetry energy as a function of the density, finding good agreement with available experimentally derived constraints at saturation and twice saturation density. We find that the corresponding pressure is also in excellent agreement with recent constraints extracted from gravitational waves of the neutron-star merger GW170817 by the LIGO-Virgo detection.

Recent observations of neutron-star properties, in particular the recent detection of gravitational waves emitted from binary neutron stars, GW 170817, open the way to put strong constraints on nuclear interactions. In this paper, we review the state of the art in calculating the equation of state of strongly interacting matter from first principle calculations starting from microscopic interactions among nucleons. We then review selected properties of neutron stars that can be directly compared with present and future observations.

The longitudinal and transverse density- and spin-density response functions in Pure Neutron Matter (PNM) are derived over a wide range of densities within the Time Dependent Local Spin Density Approximation (TDLSDA). The underlying density functional was derived from an Auxiliary Field Diffusion Monte Carlo (AFDMC) computation of the equation of state of unpolarized and fully spin polarized pure neutron matter. In order to assess the dependence of the results on the specific underlying nucleon-nucleon Hamiltonian, we used both the phenomenological Argonne AV8'+UIX force, and local chiral forces up to N2LO. The resulting response function has then been applied to the study of the neutrino mean free path in PNM.

Quantum Monte Carlo methods have proved very valuable to study the structure and reactions of light nuclei and nucleonic matter starting from realistic nuclear interactions and currents. These ab-initio calculations reproduce many low-lying states, moments and transitions in light nuclei, and simultaneously predict many properties of light nuclei and neutron matter over a rather wide range of energy and momenta. We review the nuclear interactions and currents, and describe the continuum Quantum Monte Carlo methods used in nuclear physics. These methods are similar to those used in condensed matter and electronic structure but naturally include spin-isospin, tensor, spin-orbit, and three-body interactions. We present a variety of results including the low-lying spectra of light nuclei, nuclear form factors, and transition matrix elements. We also describe low-energy scattering techniques, studies of the electroweak response of nuclei relevant in electron and neutrino scattering, and the properties of dense nucleonic matter as found in neutron stars. A coherent picture of nuclear structure and dynamics emerges based upon rather simple but realistic interactions and currents.

The onset of hyperons in the core of neutron stars and the consequent softening of the equation of state have been questioned for a long time. Controversial theoretical predictions and recent astrophysical observations of neutron stars are the grounds for the so-called hyperon puzzle. We calculate the equation of state and the neutron star mass-radius relation of an infinite systems of neutrons and $\Lambda$ particles by using the auxiliary field diffusion Monte Carlo algorithm. We find that the three-body hyperon-nucleon interaction plays a fundamental role in the softening of the equation of state and for the consequent reduction of the predicted maximum mass. We have considered two different models of three-body force that successfully describe the binding energy of medium mass hypernuclei. Our results indicate that they give dramatically different results on the maximum mass of neutron stars, not necessarily incompatible with the recent observation of very massive neutron stars. We conclude that stronger constraints on the hyperon-neutron force are necessary in order to properly assess the role of hyperons in neutron stars.

We perform a systematic assessment of models for the equation of state (EOS) of dense matter in the context of recent neutron star mass and radius measurements to obtain a broad picture of the structure of neutron stars. We demonstrate that currently available neutron star mass and radius measurements provide strong constraints on moments of inertia, tidal deformabilities, and crust thicknesses. A measurement of the moment of inertia of PSR J0737-3039A with 10% error, without any other information from observations, will constrain the EOS over a range of densities to within 50%$-$60%. We find tidal deformabilities between 0.6 and $6\times 10^{36}$ g cm$^{2}$ s$^{2}$ (to 95% confidence) for $M=1.4~\mathrm{M}_{\odot}$, and any measurement which constrains this range will provide an important constraint on dense matter. The crustal fraction of the moment of inertia can be as large as 10% for $M=1.4~\mathrm{M}_{\odot}$ permitting crusts to have a large enough moment of inertia reservoir to explain glitches in the Vela pulsar even with a large amount of superfluid entrainment. Finally, due to the uncertainty in the equation of state, there is at least a 40% variation in the thickness of the crust for a fixed mass and radius, which implies that future simulations of the cooling of a neutron star crust which has been heated by accretion will need to take this variation into account.

Neutron matter is an intriguing nuclear system with multiple connections to other areas of physics. Considerable progress has been made over the last two decades in exploring the properties of pure neutron fluids. Here we begin by reviewing work done to explore the behavior of very low density neutron matter, which forms a strongly paired superfluid and is thus similar to cold Fermi atoms, though at energy scales differing by many orders of magnitude. We then increase the density, discussing work that ties the study of neutron matter with the determination of the properties of neutron-rich nuclei and neutron-star crusts. After this, we review the impact neutron matter at even higher densities has on the mass-radius relation of neutron stars, thereby making contact with astrophysical observations.

We present microscopic calculations of light and medium mass nuclei and the equation of state of symmetric and asymmetric nuclear matter using different nucleon-nucleon forces, including a new Argonne version that has the same spin/isospin structure as local chiral forces at next-to-next-to-leading order (N2LO). The calculations are performed using Auxiliary Field Diffusion Monte Carlo (AFDMC) combined with an improved variational wave function. We show that the AFDMC method can now be used to successfully calculate the energies of very light to medium mass nuclei as well as the energy of isospin-asymmetric nuclear matter, demonstrating microscopically the quadratic dependence of the energy on the symmetry energy.

We report the results of a calculation of the response of cold neutron matter to neutral-current interactions with low energy neutrinos, carried out using an effective interaction and effective operators consistently derived within the formalism of Correlated Basis Functions. The neutrino mean free path obtained from the calculated responses turns out to be strongly affected by both short and long range correlations, leading to a sizable increase with respect to the prediction of the Fermi gas model. The consistency between the proposed approach and Landau theory of normal Fermi liquids also has been investigated, using a set of Landau parameters obtained from the matrix elements of the effective interaction.

Recent progress in quantum Monte Carlo with modern nucleon-nucleon interactions have enabled the successful description of properties of light nuclei and neutron-rich matter. As a demonstration, we show that the agreement between theoretical calculations of the charge form factor of 12C and the experimental data is excellent. Applying similar methods to isospin-asymmetric systems allows one to describe neutrons confined in an external potential and homogeneous neutron-rich matter. Of particular interest is the nuclear symmetry energy, the energy cost of creating an isospin asymmetry. Combining these advances with recent observations of neutron star masses and radii gives insight into the equation of state of neutron-rich matter near and above the saturation density. In particular, neutron star radius measurements constrain the derivative of the symmetry energy.

We review the calculation of the equation of state of pure neutron matter using quantum Monte Carlo (QMC) methods. QMC algorithms permit the study of many-body nuclear systems using realistic two- and three-body forces in a nonperturbative framework. We present the results for the equation of state of neutron matter, and focus on the role of three-neutron forces at supranuclear density. We discuss the correlation between the symmetry energy, the neutron star radius and the symmetry energy. We also combine QMC and theoretical models of the three-nucleon interactions, and recent neutron star observations to constrain the value of the symmetry energy and its density dependence.

We study the spin response of cold dense neutron matter in the limit of zero momentum transfer, and show that the frequency dependence of the long-wavelength spin response is well constrained by sum-rules and the asymptotic behavior of the two-particle response at high frequency. The sum-rules are calculated using Auxiliary Field Diffusion Monte Carlo technique and the high frequency two-particle response is calculated for several nucleon-nucleon potentials. At nuclear saturation density, the sum-rules suggest that the strength of the spin response peaks at $\omega \simeq$ 40--60 MeV, decays rapidly for $\omega \geq $100 MeV, and has a sizable strength below 40 MeV. This strength at relatively low energy may lead to enhanced neutrino production rates in dense neutron-rich matter at temperatures of relevance to core-collapse supernova.

We present an overview of microscopical calculations of the Equation of State (EOS) of neutron matter performed using Quantum Monte Carlo techniques. We focus to the role of the model of the three-neutron force in the high-density part of the EOS up to a few times the saturation density. We also discuss the interplay between the symmetry energy and the neutron star mass-radius relation. The combination of theoretical models of the EOS with recent neutron stars observations permits us to constrain the value of the symmetry energy and its slope. We show that astrophysical observations are starting to provide important insights into the properties of neutron star matter.

The symmetry energy contribution to the nuclear Equation of State (EoS) impacts various phenomena in nuclear astrophysics, nuclear structure, and nuclear reactions. Its determination is a key objective of contemporary nuclear physics with consequences for the understanding of dense matter within neutron stars. We examine the results of laboratory experiments that have provided initial constraints on the nuclear symmetry energy and its density dependence at and somewhat below normal nuclear matter density. Some of these constraints have been derived from properties of nuclei. Others have been derived from the nuclear response to electroweak and hadronic probes. We also examine the most frequently used theoretical models that predict the symmetry energy and its slope. By comparing existing constraints on the symmetry pressure to theories, we demonstrate how the contribution of the three-body force, an essential ingredient in neutron matter models, can be determined.

We calculate the equation of state of neutron matter with realistic two- and three-nucleon interactions using quantum Monte Carlo techniques, and illustrate that the short-range three-neutron interaction determines the correlation between neutron matter energy at nuclear saturation density and higher densities relevant to neutron stars. Our model also makes an experimentally testable prediction for the correlation between the nuclear symmetry energy and its density dependence -- determined solely by the strength of the short-range terms in the three neutron force. The same force provides a significant constraint on the maximum mass and radius of neutron stars.

Using a phenomenological form of the equation of state of neutron matter near the saturation density which has been previously demonstrated to be a good characterization of quantum Monte Carlo simulations, we show that currently available neutron star mass and radius measurements provide a significant constraint on the equation of state of neutron matter. At higher densities we model the equation of state using polytropes and a quark matter model, and we show that our results do not change strongly upon variation of the lower boundary density where these polytropes begin. Neutron star observations offer an important constraint on a coefficient which is directly connected to the strength of the three-body force in neutron matter, and thus some theoretical models of the three-body may be ruled out by currently available astrophysical data. In addition, we obtain an estimate of the symmetry energy of nuclear matter and its slope that can be directly compared to the experiment and other theoretical calculations.

The problem of accurately determining the equation of state of nuclear and neutron matter at density near and beyond saturation is still an open challenge. In this paper we will review the most recent progress made by means of Quantum Monte Carlo calculations, which are at present the only ab-inito method capable to treat a sufficiently large number of particles to give meaningful estimates depending only on the choice of the nucleon-nucleon interaction. In particular, we will discuss the introduction of density-dependent interactions, the study of the temperature dependence of the equation of state, and the possibility of accurately studying the effect of the onset of hyperons by developing an accurate hyperon-nucleon and hyperon-nucleon-nucleon interaction.

We present results for neutron star models constructed with a new equation of state for nuclear matter at zero temperature. The ground state is computed using the Auxiliary Field Diffusion Monte Carlo (AFDMC) technique, with nucleons interacting via a semi-phenomenological Hamiltonian including a realistic two-body interaction. The effect of many-body forces is included by means of additional density-dependent terms in the Hamiltonian. In this letter we compare the properties of the resulting neutron-star models with those obtained using other nuclear Hamiltonians, focusing on the relations between mass and radius, and between the gravitational mass and the baryon number.

We calculate the equation of state of neutron matter at zero temperature by means of the auxiliary field diffusion Monte Carlo method (AFDMC) combined with a fixed-phase approximation. The calculation of the energy is carried out by simulating up to 114 neutrons in a periodic box. Special attention was made to reduce finite size effects at the energy evaluation by adding to the interaction the effect due to the truncation of the simulation box, and by performing several simulations using different number of neutrons. The finite size effects due to the kinetic energy were also checked by employing the twist--averaged boundary conditions. We considered a realistic nuclear Hamiltonian containing modern two-- and three--body interactions of the Argonne and Urbana family. The equation of state can be used to compare and to calibrate other many-body calculations and to predict properties of neutron stars.

We present a Quantum Monte Carlo study of the zero temperature equation of state of neutron matter and the computation of the $^1S_0$ pairing gap in the low-density regime with $\rho<0.04$ fm$^{-3}$. The system is described by a non-relativistic nuclear Hamiltonian including both two-- and three--nucleon interactions of the Argonne and Urbana type. This model interaction provides very accurate results in the calculation of the binding energy of light nuclei. A suppression of the gap with respect to the pure BCS theory is found, but sensibly weaker than in other works that attempt to include polarization effects in an approximate way.

In this thesis, I discuss the use of the Auxiliary Field Diffusion Monte Carlo method to compute the ground state of nuclear Hamiltonians, and I show several applications to interesting problems both in nuclear physics and in nuclear astrophysics. In particular, the AFDMC algorithm is applied to the study of several nuclear systems, finite, and infinite matter. Results about the ground state of nuclei ($^4$He, $^8$He, $^{16}$O and $^{40}$Ca), neutron drops (with 8 and 20 neutrons) and neutron rich-nuclei (isotopes of oxygen and calcium) are discussed, and the equation of state of nuclear and neutron matter are calculated and compared with other many-body calculations. The $^1S_0$ superfluid phase of neutron matter in the low-density regime was also studied.