Strangeness in nuclei and neutron stars: a challenging puzzle
Résumé
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 threebody hyperonnucleon interactions based on the available experimental data for light and mediumheavy hypernuclei.
Strangeness in nuclei and neutron stars: a challenging puzzle
Diego Lonardoni, Alessandro Lovato, Stefano Gandolfi, et Francesco Pederiva
1Physics Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
\@textsuperscript2Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
\@textsuperscript3Department of Physics, University of Trento, Via Sommarive 14, I38123 Trento, Italy
\@textsuperscript4INFNTrento Institute for Fundamental Physics and Application, Trento, Italy

Résumé.
1 Introduction
Neutron stars are among the most compact and dense objects in the Universe, with typical masses up to and radii km. Their central densities can be several times larger than nuclear saturation density, . At densities larger than , 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 of the millisecond pulsars PSR J16142230 [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 manybody techniques. Its solution requires a more accurate theoretical analysis and a thorough experimental investigation of the hyperonnucleon 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 hyperonnucleon 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/isospindependent Hamiltonians. It has been successfully applied to solve the manybody problem in a non perturbative fashion for strongly correlated systems, both in the standard [4, 5, 6, 7] and strange nuclear sector [8, 9, 11, 10].
2 Hypernuclei: from few to manybody systems
As reported in Refs. [8, 9, 10], an accurate analysis of the separation energy of light and mediumheavy hypernuclei has been recently carried out using a phenomenological hypernuclear interaction analog to the ArgonneIllinois nucleonnucleon force [12, 13, 14, 15, 16, 17, 18, 19]. The twobody component of this interaction has been fit to the existing scattering data. The threebody sector instead has been refit 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 threebody hyperonnucleon interaction is the key to satisfactorily reproduce the ground state properties of hypernuclei from light to mediumheavy within a unique theoretical framework. Fig. 1 shows the AFDMC results for when using the complete 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 twobody charge symmetry breaking interaction has been also employed [13, 9]. As shown in Tab. 1, AFDMC results reproduce the hyperon separation energies and energy splitting of the mirror hypernuclei H and He, according to the old emulsion data [20]. H is weakly bound compared to deuterium and the neutral hypernucleus n is predicted to be unbound, in agreement with other theoretical results [21, 22, 23], contradicting the possible experimental evidence reported in Ref. [24].
System  

n  unbound  unknown  —  — 
H  0.23(9)  0.13(5)  
H  1.95(9)  2.04(4)  0.42(11)  0.35(6) 
He  2.37(9)  2.39(3) 
In order to validate the phenomenological approach adopted in quantum Monte Carlo calculations for hypernuclei, a benchmark among different fewbody methods has been recently proposed. In tab. 2 we report preliminary results for He for two different exact methods, AFDMC and Nonsymmetrized Hyperspherical Harmonics (HSHH) [25]. In both strange and nonstrange sectors the employed potentials are at twobody level only. Results agree among different methods whether using a soft (Minnesota [26]) or hard (Argonne V4 [27]) potential.
3 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 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 threebody hyperonnucleon force shows the difficulty in extracting the information on the Hamiltonian from currently available experimental information on hypernuclei [10].
Remerciements
This work was supported by the U.S. Department of Energy, Office of Science, Office of Nuclear Physics, under the NUCLEI SciDAC grant (D.L., A.L., S.G.), by the Department of Energy, Office of Science, Office of Nuclear Physics, under Contract No. DEAC0206CH11357 (A.L.), and by DOE under Contract No. DEAC5206NA25396 and Los Alamos LDRD grant (S.G.). This research used resources of the National Energy Research Scientific Computing Center (NERSC), which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DEAC0205CH11231, and computing time provided by Institutional Computing (IC) at LANL.
Références
 [1] P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. Roberts, and J. W. T. Hessels, Nature 467, 1081 (2010).
 [2] J. Antoniadis et al., Science 340, 1233232 (2013).
 [3] K. E. Schmidt and S. Fantoni, Phys. Lett. B 446 (1999) 99.
 [4] S. Gandolfi, J. Carlson, and S. C. Pieper, Phys. Rev. Lett. 106, 012501 (2011).
 [5] S. Gandolfi, J. Carlson, and S. Reddy, Phys. Rev. C 85, 032801 (2012).
 [6] S. Gandolfi, A. Lovato, J. Carlson, and K. E. Schmidt, Phys. Rev. C 90, 061306 (2014).
 [7] J. Carlson, S. Gandolfi, F. Pederiva, S. C. Pieper, R. Schiavilla, K. E. Schmidt, and R. B. Wiringa arXiv:1412.3081 (2015).
 [8] D. Lonardoni, S. Gandolfi, and F. Pederiva, Phys. Rev. C 87, 041303(R) (2013).
 [9] D. Lonardoni, F. Pederiva, and S. Gandolfi, Phys. Rev. C 89, 014314 (2014).
 [10] F. Pederiva, F. Catalano, D. Lonardoni, A. Lovato, and S. Gandolfi, arXiv:1506.04042 (2015).
 [11] D. Lonardoni, A. Lovato, S. Gandolfi, and F. Pederiva, Phys. Rev. Lett. 114, 092301 (2015).
 [12] A. R. Bodmer, Q. N. Usmani, and J. Carlson, Phys. Rev. C 29, 684 (1984).
 [13] A. R. Bodmer and Q. N. Usmani, Phys. Rev. C 31, 1400 (1985).
 [14] A. R. Bodmer and Q. N. Usmani, Nucl. Phys. A 477, 621 (1988).
 [15] A. A. Usmani, Phys. Rev. C 52, 1773 (1995).
 [16] A. A. Usmani, S. C. Pieper, and Q. N. Usmani, Phys. Rev. C 51, 2347 (1995).
 [17] Q. N. Usmani and A. R. Bodmer, Phys. Rev. C 60, 055215 (1999).
 [18] A. A. Usmani and F. C. Khanna, J. Phys. G 35, 025105 (2008).
 [19] M. Imran, A. A. Usmani, M. Ikram, Z. Hasan, and F. C. Khanna, J. Phys. G 41, 065101 (2014).
 [20] M. Jurič et al., Nucl. Phys. B 52, 130 (1973).
 [21] E. Hiyama, S. Ohnishi, B. F. Gibson, and Th. A. Rijken, Phys. Rev. C 89, 061302(R) (2014).
 [22] A. Gal and H. Garcilazo, Physics Letters B 736 (2014) 93–97.
 [23] J.M. Richard, Q. Wang, and Q. Zhao, Phys. Rev. C 91, 014003 (2015).
 [24] C. Rappold et al., Phys. Rev. C 88, 041001 (2013).
 [25] S. Deflorian, N. Barnea, W. Leidemann, and G. Orlandini, FewBody Syst. 54, 1879 (2013).
 [26] D. R. Thompson, M. Lemere, and Y. C. Tang, Nucl. Phys. A 286 (1977) 53 .
 [27] R. B. Wiringa and S. C. Pieper, Phys. Rev. Lett. 89,18 (2002).