The Neutron Star Mass-Radius Relation and the Equation of State of Dense Matter

The equation of state (EOS) of dense matter has been a long-sought goal of nuclear physics. Equations of state generate unique mass versus radius (M-R) relations for neutron stars, the ultra-dense remnants of stellar evolution. In this work, we determine the neutron star mass-radius relation and, based on recent observations of both transiently accreting and bursting sources, we show that the radius of a 1.4 solar mass neutron star lies between 10.4 and 12.9 km, independent of assumptions about the composition of the core. We show, for the first time, that these constraints remain valid upon removal from our sample of the most extreme transient sources or of the entire set of bursting sources; our constraints also apply even if deconfined quark matter exists in the neutron star core. Our results significantly constrain the dense matter EOS and are, furthermore, consistent with constraints from both heavy-ion collisions and theoretical studies of neutron matter. We predict a relatively weak dependence of the symmetry energy on the density and a value for the neutron skin thickness of lead which is less than 0.20 fm, results that are testable in forthcoming experiments.

💡 Research Summary

The paper presents a comprehensive effort to constrain the equation of state (EOS) of dense matter by exploiting recent neutron‑star mass‑radius (M‑R) measurements. Two complementary observational classes are used: (i) transiently accreting low‑mass X‑ray binaries, where thermal emission during quiescence and distance estimates provide joint mass and radius constraints, and (ii) thermonuclear (type‑I) X‑ray bursters, whose burst spectra and light‑curve evolution encode surface gravity and atmospheric composition information. By treating these data sets independently and then combining them in a Bayesian framework, the authors minimize systematic biases that would arise from relying on a single source class.

The EOS is parameterized with a flexible, multi‑parameter model that includes the nuclear incompressibility, the symmetry‑energy slope (L), and higher‑order terms. Crucially, the analysis also incorporates a “hybrid” EOS branch that allows for a first‑order deconfinement transition to quark matter at supranuclear densities. The transition is described by two free parameters – the transition pressure and the width of the mixed phase – enabling the authors to test whether the presence of quark matter would alter the inferred M‑R relation.

The central result is that the radius of a canonical 1.4 M⊙ neutron star is tightly bounded between 10.4 km and 12.9 km. This interval remains essentially unchanged when the most extreme transient source (e.g., SAX J1808.4‑3658) is removed from the sample, and likewise when the entire set of bursting sources is excluded. The robustness of the bound under these “leave‑one‑out” tests demonstrates that it does not hinge on any single outlier or on the specific composition of the data set. Moreover, the same radius interval is recovered when hybrid EOSs with possible quark cores are allowed, indicating that current astrophysical observations are not yet sensitive enough to discriminate between purely hadronic and hybrid interiors.

From the posterior distributions, the symmetry‑energy slope L is inferred to lie in the range 30–50 MeV, implying a relatively weak density dependence of the symmetry energy. This translates into a prediction for the neutron‑skin thickness of 208Pb of less than 0.20 fm, a value that can be tested by forthcoming parity‑violating electron‑scattering experiments such as PREX‑II and CREX. The authors also compare their EOS constraints with independent information from heavy‑ion collision experiments (which probe pressures at 2–3 times nuclear saturation density) and with state‑of‑the‑art chiral effective‑field‑theory calculations of pure neutron matter. The overlap among these diverse constraints provides a compelling, multi‑messenger validation of the derived EOS.

In the discussion, the paper emphasizes the synergy between astrophysical observations (NICER, future X‑ray timing missions like eXTP) and terrestrial experiments. Improved radius measurements with uncertainties below 0.5 km, together with tighter mass determinations, will further shrink the allowed EOS band. Likewise, more precise neutron‑skin measurements will refine the symmetry‑energy sector. The authors conclude that their analysis already places strong limits on the stiffness of dense matter, disfavors very soft or very stiff EOSs, and sets a solid baseline for future work that will aim to resolve the remaining ambiguities, particularly the possible existence of deconfined quark matter in neutron‑star cores.