Nuclear Symmetry Energy, Neutron Skin Thickness, and the Equation of State of Dense Matter: Constraints from PREX-II, CREX, NICER, and Gravitational-Wave Observations

The nuclear symmetry energy and its density dependence are among the most consequential yet least constrained quantities in nuclear and astrophysical physics. Governed by the isospin asymmetry of nuclear matter, the symmetry energy connects laboratory measurements of finite nuclei to the macroscopic structure and composition of neutron stars. This review examines the current state of constraints on the symmetry energy parameters — the saturation value J, the slope parameter L, and the curvature Ksym — arising from four independent and complementary sources: (i) parity-violating electron scattering experiments PREX-II and CREX, which furnish model-independent measurements of neutron skin thicknesses in ²⁰⁸Pb and ⁴⁸Ca respectively; (ii) X-ray pulse-profile modeling by NASA’s NICER mission, which has now delivered mass-radius measurements of four millisecond pulsars; (iii) gravitational-wave observations of binary neutron star mergers, most notably GW170817, constraining tidal deformability; and (iv) chiral effective field theory (χEFT) calculations of neutron matter. The PREX-II and CREX results present an enduring theoretical puzzle — favoring stiff and soft equations of state respectively — that no existing energy density functional can simultaneously reconcile within 1σ. We analyze proposed resolutions, including isovector spin-orbit interactions and density-dependent meson couplings, and discuss their implications for neutron star radii, maximum mass, and tidal deformability. The emerging multimessenger picture constrains the radius of a 1.4 M☉ neutron star to R₁.₄ ≈ 11.4–2.5 km and the tidal deformability to Λ₁.₄ ≈ 70–580, pointing toward a moderate-stiffness equation of state with possible phase softening near twice nuclear saturation density.