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.

We report benchmark calculations of the energy per particle of pure neutron matter as a function of the baryon density using three independent many-body methods: Brueckner-Bethe-Goldstone, Fermi hypernetted chain/single-operator chain, and auxiliary-field diffusion Monte Carlo. Significant technical improvements are implemented in the latter two methods. The calculations are made for two distinct families of realistic coordinate-space nucleon-nucleon potentials fit to scattering data, including the standard Argonne $v_{18}$ interaction and two of its simplified versions, and four of the new Norfolk $\Delta$-full chiral effective field theory potentials. The results up to twice nuclear matter saturation density show some divergence among the methods, but improved agreement compared to earlier work. We find that the potentials fit to higher-energy nucleon-nucleon scattering data exhibit a much smaller spread of energies.

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.

In a core-collapse supernova, a huge amount of energy is released in the Kelvin-Helmholtz phase subsequent to the explosion, when the proto-neutron star cools and deleptonizes as it loses neutrinos. Most of this energy is emitted through neutrinos, but a fraction of it can be released through gravitational waves. We model the evolution of a proto-neutron star in the Kelvin-Helmholtz phase using a general relativistic numerical code, and a recently proposed finite temperature, many-body equation of state; from this we consistently compute the diffusion coefficients driving the evolution. To include the many-body equation of state, we develop a new fitting formula for the high density baryon free energy at finite temperature and intermediate proton fraction. We estimate the emitted neutrino signal, assessing its detectability by present terrestrial detectors, and we determine the frequencies and damping times of the quasi-normal modes which would characterize the gravitational wave signal emitted in this stage.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

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 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.

We studied the Fermi and Gamow-Teller responses of cold symmetric nuclear matter within a unified dynamical model, suitable to account for both short- and long-range correlation effects. The formalism of correlated basis functions has been used to construct two-body effective interactions and one-body effective weak operators. The inclusion of the three-body cluster term allowed for incorporating in the effective interaction a realistic model of three- nucleon forces, namely the UIX potential. Moreover, the sizable unphysical dependence of the effective weak operator is removed once the three-body cluster term is taken into account.

In this thesis, the ground state properties of nuclear matter, namely the energy per particle and the response to weak probes, are computed, studying the effects of three nucleon interactions. Both the variational approach, based on the formalism of correlated basis function, and the auxiliary field diffusion Monte Carlo method have been used. A scheme suitable to construct a density-dependent two-nucleon potential in correlated basis approach is discussed. The density dependent potential resulting from UIX three-nucleon force has been employed in auxiliary field diffusion Monte Carlo calculations that turned out to be in very good agreement with correlated basis variational results. Hence, the underbinding of symmetric nuclear matter has to be ascribed to deficiencies of the UXI potential. A comparative analysis of the equations of state of both pure neutron matter and symmetric nuclear matter obtained using a new generation of "chiral inspired" local three-body potentials has been performed. These potentials provide an excellent description of the properties of light nuclei, as well as of the neutron-deuteron doublet scattering length. The weak response of symmetric nuclear matter has been computed at three-body cluster level. Two-body effective interactions and one-body effective operators have been derived within the formalism of correlated basis functions. The inclusion of the three-body cluster term in the effective interaction allowed for a direct inclusion of the UIX three-nucleon potential. Moreover, the sizable unphysical dependence of the effective weak operator is removed once the three-body cluster term is taken into account.

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.

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.