Strangeness in nuclei and neutron stars: a challenging puzzle

The prediction of neutron stars properties is strictly connected to the employed nuclear interactions. The appearance of hyperons in the inner core of the star is strongly dependent on the details of the underlying hypernuclear force. We summarize our recent quantum Monte Carlo results on the development of realistic two- and three-body hyperon-nucleon interactions based on the available experimental data for light- and medium-heavy hypernuclei.


Introduction
Neutron stars are among the most compact and dense objects in the Universe, with typical masses up to ∼ 2 M and radii R ∼ 10 km. Their central densities can be several times larger than nuclear saturation density, ρ 0 = 0.16 fm −3 . At densities larger than ∼ 2ρ 0 , when the nucleon chemical potential is large enough, the conversion of nucleons into hyperons might become energetically favorable, and Pauli blocking would prevent hyperons from decaying by limiting the phase space available to nucleons. This would lead to a reduction of the Fermi pressure exerted by the baryons and to a softening of the equation of state (EOS). As a consequence, the maximum mass determined by the equilibrium condition between gravitational and nuclear forces would be reduced.
Currently there is no general agreement (even qualitative) among the predicted results for the EOS of strange baryonic matter and the maximum mass of a neutron star (NS), that violates the constraints given by the recently measured value of ∼ 2 M of the millisecond pulsars PSR J1614-2230 [1] and PSR J0348+0432 [2], see Ref. [11]. This apparent inconsistency between NS mass observations and theoretical calculations is a long standing problem known as the hyperon puzzle. It has to be ascribed to a combination of an incomplete knowledge of the forces governing the system and to the concurrent use of approximate theoretical many-body techniques. Its solution requires a more accurate theoretical analysis and a thorough experimental investigation of the hyperon-nucleon force in a variety of systems, ranging from light to medium and heavy hypernuclei.
We employ a quantum Monte Carlo computational technique to tackle the general problem of the hyperon-nucleon interaction. The algorithm, auxiliary field diffusion Monte Carlo (AFDMC), has been introduced by Schmidt and Fantoni [3] as an extension of the usual diffusion Monte Carlo method to deal in an efficient way with spin/isospin-dependent Hamiltonians. It has been successfully applied to solve the many-body problem in a non perturbative fashion for strongly correlated systems, both in the standard [4][5][6][7] and strange nuclear sector [8][9][10][11]. arXiv:1508.04722v1 [nucl-th] 19 Aug 2015 EPJ Web of Conferences 2 Hypernuclei: from few-to many-body systems As reported in Refs. [8][9][10], an accurate analysis of the Λ separation energy B Λ of light-and medium-heavy hypernuclei has been recently carried out using a phenomenological hypernuclear interaction analog to the Argonne-Illinois nucleon-nucleon force [12][13][14][15][16][17][18][19]. The ΛN two-body component of this interaction has been fit to the existing Λp scattering data. The ΛNN three-body sector instead has been re-fit to the available Λ separation energies of single-Λ hypernuclei over a wide mass range by means of AFDMC calculations. It has been shown that the repulsive nature of the three-body hyperon-nucleon interaction is the key to satisfactorily reproduce the ground state properties of hypernuclei from light-to medium-heavy within a unique theoretical framework. Fig. 1 shows the AFDMC results for B Λ when using the complete ΛN + ΛNN force. The agreement with experimental data is good over all the mass range investigated and for the Λ particle in different single particle states. For more details see Refs. [8][9][10].
In the light sector, a two-body ΛN charge symmetry breaking interaction has been also employed [9,13]. As shown in Tab. 1, AFDMC results reproduce the hyperon separation energies and energy splitting of the mirror hypernuclei 4 Λ H and 4 Λ He, according to the old emulsion data [20]. 3 Λ H is weakly bound compared to deuterium and the neutral hypernucleus 3 Λ n is predicted to be unbound, in agreement with other theoretical results [21][22][23], contradicting the possible experimental evidence reported in Ref. [24].
In order to validate the phenomenological approach adopted in quantum Monte Carlo calculations for hypernuclei, a benchmark among different few-body methods has been recently proposed. In tab. 2 we report preliminary results for 5 Λ He for two different exact methods, AFDMC and Non- Measured (see Ref. [10] for the complete list of references) and computed Λ separation energies as a function of A −2/3 . Results for the Λ particle in different single particle states are shown.
21st International Conference on Few-Body Problems in Physics Table 1. AFDMC and experimental [20] Λ separation energies for the lightest hypernuclei. In the case of A = 4 the separation energy difference ∆B Λ is also shown. Evidence for the neutral hypernucleus 3 Λ n has been reported in Ref. [24], but no B Λ value has been reported. All energies are in MeV.  (3) symmetrized Hyperspherical Harmonics (HSHH) [25]. In both strange and non-strange sectors the employed potentials are at two-body level only. Results agree among different methods whether using a soft (Minnesota [26]) or hard (Argonne V4 [27]) NN potential.

Neutron stars: from hypernuclei to the infinite medium
Precise AFDMC calculations have been performed to first assess the parameters of the interaction from measured hyperon separation energies, and then to extrapolate the results to the case of neutron matter with strangeness. However, parametrizations of the potential predicting relatively small differences in the Λ separation energies of hypernuclei give very different results for the properties of the infinite medium [11]. The resulting EOS spans the whole regime extending from the appearance of a substantial fraction of hyperons at ∼ 2ρ 0 to the absence of Λ particles in the entire density range of the star. This yields a sizable effect on the predicted NS structure. Therefore, the derivation of realistic hypernuclear potential models is of paramount importance to properly assess the role of hyperons in NSs and reconcile theoretical predictions with astrophysical observations. This demands a precise and systematic experimental investigation of hypernuclear properties over a wide range of masses. For instance, a recent study of the isospin dependence of the present three-body hyperon-nucleon force shows the difficulty in extracting the information on the Hamiltonian from currently available experimental information on hypernuclei [10].