Central Speed of Sound, Trace Anomaly and Observables of Neutron Stars from Perturbative Analyses of Scaled TOV Equations

TL;DR

This paper uses perturbative solutions of scaled TOV equations to analyze neutron star properties in a model-independent way, revealing bounds on the speed of sound, pressure, and radius that align with observations and suggest a continuous crossover in NS cores.

Contribution

It introduces a perturbative approach to study neutron star properties independently of specific nuclear EOS models, providing new bounds and insights into NS structure and phase transitions.

Findings

Speed of sound increases with central pressure and exceeds the conformal bound in massive NSs.

Upper bound on P/ε ratio is approximately 0.374 at NS centers, indicating intrinsic strong-field gravity effects.

A new causality boundary relates maximum NS radius and mass, consistent with observational data.

Abstract

The central speed of sound (SS) measures the stiffness of the Equation of State (EOS) of superdense neutron star (NS) matter. Its variations with density and radial coordinate in NSs in conventional analyses often suffer from uncertainties of the specific nuclear EOSs used. Using the central SS and NS mass/radius scaling obtained from solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations, we study the variations of SS, trace anomaly and several closely related properties of NSs in an EOS-model independent manner. We find that the SS increases with the reduced central pressure $\widehat{P}_{\rm{c}}\equiv P_{\rm{c}}/\varepsilon_{\rm{c}}$ (scaled by the central energy density $\varepsilon_{\rm{c}}$), and the conformal bound for SS tends to break down for NSs with masses higher than about 1.9$M_{\odot}$. The ratio $P/\varepsilon$ is upper bounded as $P/\varepsilon\lesssim0.374$ around the centers of stable NSs. We demonstrate that it is an intrinsic property of strong-field gravity and is more relevant than the perturbative QCD bound on it. While a sharp phase transition at high densities characterized by a sudden vanishing of SS in cores of massive NSs are basically excluded, the probability for a continuous crossover signaled by a peaked radial profile of SS is found to be enhanced as $\widehat{P}_{\rm{c}}$ decreases, implying it likely happens near the centers of massive NSs. Moreover, a new and more stringent causality boundary as $R_{\max}/\rm{km}\gtrsim 4.73M_{\rm{NS}}^{\max}/M_{\odot}+1.14$ for NS M-R curve is found to be excellently consistent with observational data on NS masses and radii. Furthermore, new constraints on the ultimate energy density and pressure allowed in NSs before collapsing into black holes are obtained and compared with earlier predictions in the literature.