2025ApJ...985..238Y

The recently developed nuclear effective interaction based on the so-called N3LO Skyrme pseudopotential is extended to include hyperon-nucleon and hyperon-hyperon interactions by assuming similar density, momentum, and isospin dependence as for the nucleon-nucleon interaction. The parameters in these interactions are determined from either experimental information, if any, or chiral effective field theory or lattice quantum chromodynamics calculations of the hyperon potentials in nuclear matter around nuclear saturation density ρ₀. We find that varying the high-density behavior of the symmetry energy E_sym(ρ) can significantly change the critical density for hyperon appearance in neutron stars and thus the maximum mass M_TOV of static hyperon stars. In particular, a symmetry energy that is soft around 2ρ₀-3ρ₀ but stiff above about 4ρ₀ can lead to M_TOV ≳ 2M_⊙ for hyperon stars and simultaneously be compatible with (1) the constraints on the equation of state of symmetric nuclear matter at suprasaturation densities obtained from flow data in heavy-ion collisions; (2) the microscopic calculations of the equation of state for pure neutron matter; (3) the tidal deformability of stars extracted from gravitational wave signal GW 170817; (4) the mass-radius relations of PSR J0030+0451, PSR J0740+6620, and PSR J0437-4715 measured from NICER; and (5) the observation of an unusually low mass and small radius in the central compact object of HESS J1731-347. Furthermore, the squared sound speed of the hyperon star matter naturally displays a strong peak structure around a baryon density of 3ρ₀-4ρ₀, consistent with a model-independent analysis of the multimessenger data. Our results suggest that the high-density symmetry energy could be a key to the solution of the hyperon puzzle in neutron star physics.