Beyond ρ 2/3 scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

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 Esym(ρ ) 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 Esym(ρ ) E sym ( ρ ) scales as (ρ /ρ ₀)2/3 ( ρ / ρ 0 ) 2 / 3 near the saturation density ρ ₀ ρ 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’ ρ 2/3 ρ 2 / 3 scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to 2ρ ₀ 2 ρ 0 but breaks down at higher densities, (2) Experimental and theoretical findings about Esym(ρ ) E sym ( ρ ) up to 2ρ ₀ 2 ρ 0 are broadly consistent, but uncertainties remain large for its curvature Ksym K sym and higher-order parameters, (3) Above 2ρ ₀ 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 Esym(ρ ) E sym</mml