Beyond $$\rho ^{2/3}$$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

Abstract The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $$E_{\textrm{sym}}(\rho )$$ E sym ( ρ ) encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $$E_{\textrm{sym}}(\rho )$$ E sym ( ρ ) scales as $$(\rho /\rho _0)^{2/3}$$ ( ρ / ρ 0 ) 2 / 3 near the saturation density $$\rho _0$$ ρ 0 of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $$\rho ^{2/3}$$ ρ 2 / 3 scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $$2\rho _0$$ 2 ρ 0 but breaks down at higher densities, (2) Experimental and theoretical findings about $$E_{\textrm{sym}}(\rho )$$ E sym ( ρ ) up to $$2\rho _0$$ 2 ρ 0 are broadly consistent, but uncertainties remain large for its curvature $$K_{\textrm{sym}}$$ K sym and higher-order parameters, (3) Above $$2\rho _0$$ 2 ρ 0 , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $$E_{\textrm{sym}}(\rho )$$ E sym</mml