A density-dependent two-nucleon potential has been derived in the formalism of correlated basis function. The effects of 3-particle interactions has been included by integrating out the degrees of freedom of the third nucleon. The potential can be easily employed in nuclear matter calculations. It yields results in agreement with those obtained from the underlying three-body potential. The use of the density dependent potential allowed us to study the effects of three-nucleon interactions in symmetric nuclear matter within the Auxiliary Field Diffusion Monte Carlo (AFDMC) computational scheme.

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.

A new generation of local three-body potentials providing an excellent description of the properties of light nuclei, as well as of the neutron-deuteron doublet scattering length, has been recently derived. We have performed a comparative analysis of the equations of state of both pure neutron matter and symmetric nuclear matter obtained using these models of three-nucleon forces. None of the considered potentials simultaneously explains the empirical equilibrium density and binding energy of symmetric nuclear matter. However, two of them provide reasonable values of the saturation density. The ambiguity concerning the treatment of the contact term of the chiral inspired potentials is discussed.

Neutron star crustal matter, whose properties are relevant in many models aimed at explaining observed astrophysical phenomena, has so far always been studied using a mean field approach. In order to check the results obtained in this way, a sensible next step is to make use of a realistic nuclear potential. The present paper extends the periodic-box Fermi HyperNetted Chain method to include longitudinal-isospin dependence of the correlations, making feasible a study of asymmetric crustal matter. Results are presented for the symmetry energy, the low-density neutron star equation of state and the single particle neutron and proton energies.

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.

The $^1S_0$ pairing in neutron matter has been investigated in presence of realistic two-- and three--nucleon interactions. We have adopted the Argonne $v_{8^\prime}$ NN and the Urbana IX 3N potentials. Quantum Monte Carlo theory, specifically the Auxiliary Field Diffusion Monte Carlo method, and Correlated Basis Function theory are employed in order to get quantitative and reliable estimates of the gap. They both fully take into account the medium modifications due to the interaction induced correlations. The two methods are in good agreement up to the maximum gap density and both point to a slight reduction with respect to the standard BCS value. In fact, the maximum gap is about $2.5 \text{MeV}$ at $k_F \sim 0.8 \text{fm}^{-1}$ in BCS and 2.3--$2.4 \text{MeV}$ at $k_F \sim 0.6 \text{fm}^{-1}$ in correlated matter. At higher densities the Quantum Monte Carlo gap becomes close to BCS. In general, the computed medium polarization effects are much smaller than those previously estimated within \emph{all theories}. Truncations of Argonne $v_{8^\prime}$ to simpler forms give the same gaps in BCS, provided the truncated potentials have been refitted to the same NN data set. Differences among the models appear in the correlated theories, most of the reduction being attributable to the tensor force. The three--nucleon interaction provides an additional increase of the gap of about 0.35 MeV.