No Arabic abstract
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-Urbana 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 consider the problem of including $Lambda$ hyperons into the ab initio framework of nuclear lattice effective field theory. In order to avoid large sign oscillations in Monte Carlo simulations, we make use of the fact that the number of hyperons is typically small compared to the number of nucleons in the hypernuclei of interest. This allows us to use the impurity lattice Monte Carlo method, where the minority species of fermions in the full nuclear Hamiltonian is integrated out and treated as a worldline in Euclidean projection time. The majority fermions (nucleons) are treated as explicit degrees of freedom, with their mutual interactions described by auxiliary fields. This is the first application of the impurity lattice Monte Carlo method to systems where the majority particles are interacting. Here, we show how the impurity Monte Carlo method can be applied to compute the binding energy of the light hypernuclei. In this exploratory work we use spin-independent nucleon-nucleon and hyperon-nucleon interactions to test the computational power of the method. We find that the computational effort scales approximately linearly in the number of nucleons. The results are very promising for future studies of larger hypernuclear systems using chiral effective field theory and realistic hyperon-nucleon interactions, as well as applications to other quantum many-body systems.
The spin susceptibility in pure neutron matter is computed from auxiliary field diffusion Monte Carlo calculations over a wide range of densities. The calculations are performed for different spin asymmetries, while using twist-averaged boundary conditions to reduce finite-size effects. The employed nuclear interactions include both the phenomenological Argonne AV8$^prime$+UIX potential and local interactions that are derived from chiral effective field theory up to next-to-next-to-leading order.
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.
In this contribution I will review some of the researches that are currently being pursued in Padova (mainly within the In:Theory and Strength projects), focusing on the interdisciplinary applications of nuclear theory to several other branches of physics, with the aim of contributing to show the centrality of nuclear theory in the Italian scientific scenario and the prominence of this fertile field in fostering new physics. In particular, I will talk about: i) the recent solution of the long-standing electron screening puzzle that settles a fundamental controversy in nuclear astrophysics between the outcome of lab experiments on earth and nuclear reactions happening in stars; the application of algebraic methods to very diverse systems such as: ii) the supramolecular complex H2@C60, i.e. a diatomic hydrogen molecule caged in a fullerene and iii) to the spectrum of hypernuclei, i.e. systems made of a Lambda particles trapped in (heavy) nuclei.
The slope of the nuclear symmetry energy at saturation density $L$ is pointed out as a crucial quantity to determine the mass and width of neutron-star crusts. This letter clarifies the relation between $L$ and the core-crust transition. We confirm that the transition density is soundly correlated with $L$ despite differences between models, and we propose a clear understanding of this correlation based on a generalised liquid drop model. Using a large number of nuclear models, we evaluate the dispersion affecting the correlation between the transition pressure $P_t$ and $L$. From a detailed analysis it is shown that this correlation is weak due to a cancellation between different terms. The correlation between the isovector coefficients $K_{rm sym}$ and $L$ plays a crucial role in this discussion.