Logo
ScienceStack
Table of Contents
Fetching ...
Sign In

Beyond $ρ^{2/3}$ Scaling: Microscopic Origins and Multimessengers of High-Density Nuclear Symmetry Energy

Bao-An Li

TL;DR

This work reviews the density dependence of the nuclear symmetry energy $E_{ m sym}(\rho)$, highlighting Siemens' $\rho^{2/3}$ scaling as a useful benchmark near nuclear saturation but not universal at suprasaturation densities. It derives the microscopic links between $E_{ m sym}$ and single-particle potentials via the HVH theorem, and discusses the origin of the scaling from the momentum dependence of the isoscalar potential and weak density dependence of the isovector sector. The authors summarize substantial empirical and theoretical support for the scaling up to about $2\rho_0$, while detailing multiple microscopic mechanisms (tensor forces, SRC, three-body forces, relativistic effects) that can break it at higher densities; they advocate a multimessenger program combining neutron-star observations and heavy-ion experiments to constrain high-density $E_{ m sym}(\rho)$ and the supradense EOS. The review further surveys neutron-star inversion analyses, Bayesian inferences, and heavy-ion observables (flow, pion, and strange particle ratios) as complementary probes, acknowledging model dependencies and the need for coordinated, high-precision measurements in the coming era of advanced detectors and radioactive-beam facilities.

Abstract

Nuclear symmetry energy $E_{\mathrm{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 affects significantly 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 structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\mathrm{sym}}(ρ)$ scales as $(ρ/ρ_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}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2ρ_0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\mathrm{sym}}(ρ)$ up to $2ρ_0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\mathrm{sym}}$ and higher-order parameters, (3) Above $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 multimessengers from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining precisely the high-density $E_{\mathrm{sym}}(ρ)$ and the EOS of supradense neutron-rich matter.

Beyond $ρ^{2/3}$ Scaling: Microscopic Origins and Multimessengers of High-Density Nuclear Symmetry Energy

TL;DR

This work reviews the density dependence of the nuclear symmetry energy , highlighting Siemens' scaling as a useful benchmark near nuclear saturation but not universal at suprasaturation densities. It derives the microscopic links between and single-particle potentials via the HVH theorem, and discusses the origin of the scaling from the momentum dependence of the isoscalar potential and weak density dependence of the isovector sector. The authors summarize substantial empirical and theoretical support for the scaling up to about , while detailing multiple microscopic mechanisms (tensor forces, SRC, three-body forces, relativistic effects) that can break it at higher densities; they advocate a multimessenger program combining neutron-star observations and heavy-ion experiments to constrain high-density and the supradense EOS. The review further surveys neutron-star inversion analyses, Bayesian inferences, and heavy-ion observables (flow, pion, and strange particle ratios) as complementary probes, acknowledging model dependencies and the need for coordinated, high-precision measurements in the coming era of advanced detectors and radioactive-beam facilities.

Abstract

Nuclear symmetry energy 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 affects significantly 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 structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that scales as near the saturation density 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' scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to but breaks down at higher densities, (2) Experimental and theoretical findings about up to are broadly consistent, but uncertainties remain large for its curvature and higher-order parameters, (3) Above , 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 multimessengers from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining precisely the high-density and the EOS of supradense neutron-rich matter.

Paper Structure

This paper contains 30 sections, 36 equations, 32 figures, 1 table.

Table of Contents

  1. Introduction
  2. Symmetry energy in terms of the Lane potential according to the Hugenholtz-Van Hove (HVH) theorem
  3. Siemens' $\rho^{2/3}$ scaling for nuclear symmetry energy
  4. Empirical evidence for the $\rho^{2/3}$ scaling of nuclear symmetry energy around $\rho_0$
  5. Neutron-proton effective mass splitting from theories, experiments and the $\rho^{2/3}$ scaling for $E_{\rm sym}(\rho)$
  6. Neutron-proton effective mass splitting based on many-body perturbation theory using $\chi$EFT forces
  7. Neutron-proton effective mass splitting from experiments
  8. Neutron-proton effective mass splitting from $\rho^{2/3}$ scaling of $E_{\rm sym}(\rho)$ within HVH theorem
  9. Neutron-proton effective mass splitting from SHF
  10. Neutron-proton effective mass splitting as a surrogate of isospin-quartic symmetry energy
  11. Expected breakdown of the $\rho^{2/3}$ scaling at high densities and its microscopic causes
  12. Example-1: Density and momentum dependences of nucleon isoscalar and isovector potentials as well as the corresponding $E_{\rm{sym}}(\rho)$ from Gogny-Hartree-Fock (GHF)
  13. Example-2: Effects of short-range correlation (SRC) on nuclear symmetry energy
  14. High-density nuclear symmetry energy from neutron star observations
  15. Example-1: Bayesian statistical inference of high-density symmetry energy from neutron star observables
  16. ...and 15 more sections

Figures (32)

  • Figure 1: (color online) The magnitude $E_{\rm{sym}}(\rho_0)$ (upper window) and slope parameter $L$ (lower window) of nuclear symmetry energy extracted by the community from analyzing various terrestrial nuclear experiments and neutron star observations up to 2013 LiBA13
  • Figure 2: (Color Online). Updated constraints on the slope parameter $L$ of symmetry energy up to year 2023 including analyses of several recent terrestrial experiments and NS observables since GW170817 in comparison with earlier systematics and the chiral EFT ($\chi$EFT) prediction. Starting from the left are the $L$ values from (1) 24 independent analyses of neutron star observables carried out by various groups between 2017 and 2021, they give an average $L\approx58\pm 19$ MeV (thick horizontal black line) LCXZ2021; (2) the original analysis of the PREX-II data Adh21 by Reed et al. Reed:2021nqk and 3 independent analyses of PREX-II data together with different combinations of terrestrial and/or astrophysical data by Reinhard et al. Reinhard:2021utvReinhard2, Essick et al. Essick21 and Yue et al. Yue:2021yfx, respectively; (3) liquid drop model analyses using separately the PREX-I and II data together, CREX data only crex, and the combination of all data assuming they are equally reliable by Lattimer Lattimer:2023rpe; (4) charged pion ration in heavy-ion reactions at RIKEN by Estee et al. SpiRIT:2021gtq; (5) the ratio of average transverse momentum (sky blue star) and the ratio of charged particle multiplicities (black star) in isobar collisions $^{96}$Zr+$^{96}$Zr and $^{96}$Ru+$^{96}$Ru) from STAR/RHIC experiments analyzed by Xu et al. Xu:2022ikx; (6) using neutron-skin thickness of $^{208}$Pb inferred by Giacalone et al. from $^{208}$Pb + $^{208}$Pb collisions measured by the ALICE/LHC Collaboration Giacalone:2023cet, (7) the difference of charge radii of the mirror pair $^{54}$Ni-$^{54}$Fe by Pineda et al. mirror; (8) the $\chi$EFT prediction by Dirschler et al. Ohio20. The horizontal band covering $L\approx59\pm 28$ MeV is the 2016 average of 53 earlier analyses of various data mostly from terrestrial nuclear experiments Oer17LiBA13. Figure modified from those in Refs. Zhang:2022sepCai:2025nxn
  • Figure 8: Nucleon effective masses computed using Skyrme interaction SLy4, the modified Skyrme interaction BSk20, and Brueckner-Hartree-Fock results with V18+TBF, V18+UIX, and CDB+UIX interactions of Ref. Baldo2014, as functions of density in asymmetric nuclear matter with proton fraction $Y_p=0.1$. Upper panel: neutron effective mass; lower panel: proton effective mass. Taken from Ref. Duan.
  • Figure 9: A prediction for the density dependence of neutron-proton effective mass splitting Whitehead with a Many-Body Perturbative Theory (MBPT) using $\chi$EFT forces in comparison with the empirical value (red) extracted from nucleon-nucleus scattering data XHLi. The figure is modified from a plot provided to the present author by Dr. T.R. Whitehead in 2021.
  • Figure 10: A survey of neutron-proton effective mass splitting at $\rho_0$ from analyzing nuclear reaction and structure experiments in comparison with the latest chiral effective field theory prediction (its 68% confidence range is indicated by the horizontal violet lines). Taken from Ref. Li-Italy. Similar plots possibly with additional data and/or calculations can also be found in Refs. SWangYang25.
  • ...and 27 more figures