Neutrino-matter interactions play a key role in binary neutron star mergers. Thermodynamics conditions at the surfaces where neutrinos decouple from matter influence neutrino spectra, ultimately affecting the evolution of the remnant and the properties of the ejecta. In this work, we post-process results of general relativistic merger simulations employing microphysical equations of state and approximate neutrino transport to investigate the thermodynamics conditions at which weak and thermal equilibrium freezes out (equilibrium surfaces), as well as conditions at which the transition between diffusion and free-streaming regime occurs (diffusion surfaces). We find that the rest mass density and the neutrino energy are the most relevant quantities in determining the location of the decoupling surfaces. For mean energy neutrinos ($\langle E_{{\nu}_e} \rangle \approx 9~{\rm MeV}$, $\langle E_{{\bar{\nu}}_e} \rangle \approx 15~{\rm MeV}$, $\langle E_{{\nu}_{\mu,\tau}} \rangle \approx 25~{\rm MeV}$), diffusion surfaces are located around $10^{11}{\rm g~cm^{-3}}$ for all neutrino species, while equilibrium surfaces for heavy flavor neutrinos are significantly deeper (several $10^{12}{\rm g~cm^{-3}}$) than the ones of $\bar{\nu}_e$ and $\nu_e$ ($\gtrsim 10^{11}{\rm g~cm^{-3}}$). The resulting decoupling temperatures are in good agreement with the average neutrino energies ($\langle E_{\nu} \rangle \sim 3.15~T$), with the softer equation of state characterized by systematically larger decoupling temperatures ($\Delta T \lesssim 1~{\rm MeV}$). Neutrinos streaming at infinity with different energies come from very different regions of the remnant. The presence of a massive NS or of a BH in the remnant influences the neutrino thermalization process.

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.

The presence of exotic states of matter in neutron stars (NSs) is currently an open issue in physics. The appearance of muons, kaons, hyperons, and other exotic particles in the inner regions of the NS, favored by energetic considerations, is considered to be an effective mechanism to soften the equation of state (EoS). In the so-called two-families scenario, the softening of the EoS allows for NSs characterized by very small radii, which become unstable and convert into a quark stars (QSs). In the process of conversion of a NS into a QS material can be ablated by neutrinos from the surface of the star. Not only neutron-rich nuclei, but also more exotic material, such as hypernuclei or deconfined quarks, could be ejected into the atmosphere. In the NS atmosphere, atoms like H, He, and C should exist, and attempts to model the NS thermal emission taking into account their presence, with spectra modified by the extreme magnetic fields, have been done. However, exotic atoms, like muonic hydrogen $(p\,\mu^-)$ or the so-called Sigmium $(\Sigma^+\,e^-)$, could also be present during the conversion process or in its immediate aftermath. At present, analytical expressions of the wave functions and eigenvalues for these atoms have been calculated only for H. In this work, we extend the existing solutions and parametrizations to the exotic atoms $(p\,\mu^-)$ and $(\Sigma^+\,e^-)$, making some predictions on possible transitions. Their detection in the spectra of NS would provide experimental evidence for the existence of hyperons in the interior of these stars.

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 present fully non-perturbative quantum Monte Carlo calculations with non-local chiral effective field theory (EFT) interactions for the ground state properties of neutron matter. The equation of state, the nucleon chemical potentials and the momentum distribution in pure neutron matter up to one and a half times the nuclear saturation density are computed with a newly optimized chiral EFT interaction at next-to-next-to-leading order. This work opens the way to systematic order by order benchmarking of chiral EFT interactions, and \emph{ab initio} prediction of nuclear properties while respecting the symmetries of quantum chromodynamics.

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 Carlo (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 $\rho\simeq 0.09 fm^{-3}$ when a small fraction of protons is added to the system.

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.