Impact of Symmetry Energy on Sound Speed and Spinodal Decomposition in Dense Neutron-Rich Matter

Using a meta model for nuclear Equation of State (EOS) with its parameters constrained by astrophysical observations and terrestrial nuclear experiments, we examine effects of nuclear EOS especially its symmetry energy \esym term on the speed of sound squared $C^2_s(\rho)$ and the critical density $\rho_t$ where $C^2_s(\rho_t)$ vanishes (indicating the onset of spinodal decomposition) in both dense neutron-rich nucleonic matter relevant for relativistic heavy-ion collisions and the cold $n+p+e+\mu$ matter in neutron stars at $\beta$-equilibrium. Unlike in nucleonic matter with fixed values of the isospin asymmetry $\delta$, in neutron stars with a density dependent isospin profile $\delta(\rho)$ determined consistently by the $\beta$ equilibrium and charge neutrality conditions, the $C^2_s(\rho)$ almost always show a peak and then vanishes at $\rho_t$. The latter strongly depends on the high-density behavior of \esym if the skewness parameter $J_0$ characterizing the stiffness of high-density symmetric nuclear matter (SNM) EOS is not too far above its currently known most probable value of about $-190$ MeV inferred from recent Bayesian analyses of neutron star observables. Moreover, in the case of having a super-soft \esym that is decreasing with increasing density above about twice the saturation density of nuclear matter, the $\rho_t$ is significantly lower than the density where the \esym vanishes (indicating the onset of isospin-separation instability and pure neutron matter formation) in neutron star cores.