In the last years auxiliary field diffusion Monte Carlo has been used to assess the properties of hypernuclear systems, from light- to medium-heavy hypernuclei and hyper-neutron matter. One of the main findings is the key role played by the three-bod
y hyperon-nucleon-nucleon interaction in the determination of the hyperon separation energy of hypernuclei and as a possible solution to the hyperon puzzle. However, there are still aspects of the employed hypernuclear potential that remain to be carefully investigated. In this paper we show that the isospin dependence of the Lambda-NN force, which is crucial in determining the NS structure, is poorly constrained by the available experimental data.
Auxiliary Field Diffusion Monte Carlo (AFDMC) calculations have been employed to revise the interaction between $Lambda$-hyperons and nucleons in hypernuclei. The scheme used to describe the interaction, inspired by the phenomenological Argonne-Urban
a forces, is the $Lambda N+Lambda NN$ potential firstly introduced by Bodmer, Usmani et al.. Within this framework, we performed calculations on light and medium mass hypernuclei in order to assess the extent of the repulsive contribution of the three-body part. By tuning this contribution in order to reproduce the $Lambda$ separation energy in $^5_Lambda$He and $^{17}_{~Lambda}$O, experimental findings are reproduced over a wide range of masses. Calculations have then been extended to $Lambda$-neutron matter in order to derive an analogous of the symmetry energy to be used in determining the equation of state of matter in the typical conditions found in the inner core of neutron stars.
We study the problem of an impurity in fully polarized (spin-up) low density neutron matter with the help of an accurate quantum Monte Carlo method in conjunction with a realistic nucleon-nucleon interaction derived from chiral effective field theory
at next-to-next-to-leading-order. Our calculations show that the behavior of the proton spin-down impurity is very similar to that of a polaron in a fully polarized unitary Fermi gas. We show that our results can be used to put tight constraints on the time-odd parts of the energy density functional, independent of the time-even parts, in the density regime relevant to neutron-rich nuclei and compact astrophysical objects such as neutron stars and supernovae.
The dynamic response of asymmetric nuclear matter is studied by using a Time-Dependent Local Isospin Density (TDLIDA) approximation approach. Calculations are based on a local density energy functional derived by an Auxiliary Field Diffusion Monte Ca
rlo (AFDMC) calculation of bulk nuclear matter. Three types of excited states emerge: collective states, a continuum of quasi-particle-quasi-hole excitations and unstable solutions. These states are analyzed and discussed for different values of the nuclear density $rho$ and isospin asymmetry $xi=(N-Z)/A$. An analytical expression of the compressibility as a function of $rho$ and $xi$ is derived which show explicitly an instability of the neutron matter around $rhosimeq 0.09 fm^{-3}$ when a small fraction of protons is added to the system.
An ab-initio method for determining the dynamical structure function of an interacting many--body quantum system has been devised by combining a generalized integral transform method with Quantum Monte Carlo methods. As a first application, the coher
ent and, separately, the incoherent excitation spectrum of bulk atomic 4He has been computed, both in the low and intermediate momentum range. The peculiar form of the kernel in the integral transform of the dynamical structure function allows to predict, without using any model, both position and width of the collective excitations in the maxon--roton region, as well as the second collective peak. A prediction of the dispersion of the single--particle modes described by the incoherent part is also presented.
