Causality as a Universal Limiter on Neutron Star Radii ✛ ✦
Context
Neutron stars are the densest objects in the Universe, where matter is compressed to nuclear and supranuclear densities. Their structure is described by a unified equation of state (EoS) linking pressure and energy density. However, until now, observation-based radius estimates have strongly depended on the choice of a specific EoS model. The breakthrough of the new work is that it for the first time derives radius constraints based solely on the fundamental principle of the constancy of the speed of light (sound cannot outrun light) and the assumption that all neutron stars obey the same EoS. This allows separation of universal physical consequences from artifacts of model assumptions. Predicted by Fritz Zwicky and discovered as radio pulsars by Jocelyn Bell Burnell, these objects are now becoming even more precise cosmic laboratories.
Methods
The authors constructed two types of equation-of-state ensembles: a maximally agnostic one using a fractal bridge covering the entire space of thermodynamically stable and causal EoS, and a commonly used parameterization of the sound speed $c_s^2$(ε). As the low-density baseline, chiral effective field theory (χEFT) calculations were used, and at high densities, perturbative QCD. For each EoS, the hydrostatic equilibrium equation Tolman–Oppenheimer–Volkoff (TOV) was solved numerically within general relativity, and mass and radius were computed. Only EoS that allow a neutron star with a maximum mass of at least 1.97 M⊙, corresponding to observations of the most massive pulsars, were selected. This approach allowed isolating the pure effect of causality without additional model assumptions.
Results
The central result is the upper bound on the radius of a neutron star of mass 2.0 M⊙, expressed via the radius of a canonical 1.4 M⊙ star: $R_2$.0 ≤ 1.16 $R_1$.4 − 1.1 km. This limit is saturated by a special family of EoS in which, immediately after the end of the χEFT region, a phase transition occurs, and then the sound speed is maximal ($c_s^2$ = 1). The bound does not depend on the choice of prior distribution and is a direct consequence of the causality principle. When this bound is imposed on independent posterior radius distributions of three pulsars obtained by the [instrument:NICER]NICER[/instrument] observatory (PSR J0437−4715, PSR J0614−3329, and PSR J0740+6620), it turns out that the joint probability shrinks to 7.5%. This mainly happens due to the removal of an unphysical large-radius tail for the most massive pulsar J0740+6620 (initially its median radius was about 12.5 km, and after the constraint it dropped to ~11.4 km). Moreover, a comparison with a hypothetical EoS family with a conformal sound speed ($c_s^2$ = 1/3) shows that such a soft equation of state is incompatible with the data: it lies below the observational contours. Hence, in the interiors of neutron stars, the sound speed must necessarily exceed the conformal limit, and in the most likely scenarios, it approaches the speed of light.
Implications
The obtained bound for the first time provides a model-independent tool to test the consistency of neutron star radius measurements. It allows, without resorting to full Bayesian inference using EoS priors, to sharply narrow the allowed parameter region and cut off unphysical values. This is especially important for identifying possible deviations caused by exotic physics, such as the presence of dark matter in interiors. The method also demonstrates that even minimal and indisputable physical principles can yield strong constraints.
Future development
As data accumulate from [instrument:NICER]NICER[/instrument] and future missions such as STROBE-X or eXTP, the constraint will become ever tighter. Moreover, the approach can be directly generalized to tidal deformabilities measured in gravitational-wave events involving neutron stars. In the future, it will become a standard test for any new EoS model claiming to describe superdense matter.
Impact
The results will have a direct impact on the interpretation of data from X-ray observatories like [instrument:NICER]NICER[/instrument] and future telescopes, on the analysis of gravitational waves from neutron star mergers (LIGO/Virgo/KAGRA), and on the theory of nuclear matter at high densities, including searches for phase transitions.
Next steps
Immediate next steps include applying the method to other known pulsars (e.g., PSR J0030+0451) and refining the low-density input using improved χEFT calculations, which will slightly tighten the constraint.
Key open problems
The work is directly connected to the unsolved problem of the equation of state of cold superdense matter and to the question of the existence of phase transitions in neutron star cores, which remains one of the key issues in astrophysics and nuclear physics.
The knowledge machine's view
This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.
What the work does
The authors derive an upper bound on the radius of a neutron star with mass 2.0 M⊙ in terms of the radius of a canonical 1.4 M⊙ star: R2.0 ≤ 1.16 R1.4 − 1.1 km, relying only on causality and a general equation of state. They construct ensembles of EoS, solve the TOV equations, select configurations with maximum mass ≥1.97 M⊙, and show that the bound is saturated by EoS with a phase transition and maximum sound speed. When applied to NICER data, the joint probability drops to 7.5% because the large-radius tail of PSR J0740+6620 is removed.
Where its strength is
The universality of the bound, independent of model assumptions and prior distributions, is a strong result that can serve as a filter for any future radius estimates.
Where it could go next
- It might be interesting to compare this bound with the EoS stiffness analysis in a neighboring work on the stiffest limit of neutron stars, to see how causality interacts with other extreme properties.based on: The Ultimate Stiffness Limit of Neutron Stars arXiv:2601.02980v1
- Nearby is a study linking neutron skin thickness to the internal structure of neutron stars; perhaps your universal radius bound could help refine skin predictions and vice versa.based on: How the Thickness of a Neutron Skin Reveals the Heart of Dead Stars arXiv:2604.26876
- There is an archive paper analyzing how nebular lines reflect different densities; comparing it with your radius bound could provide additional checks for atmosphere and surface composition models.based on: Nebula Orchestra: Why Lines Hear Different Densities arXiv:2607.07973
Worth trying
- One could test how the bound R2.0 ≤ 1.16 R1.4 − 1.1 km affects posterior radius distributions from future missions such as eXTP or Athena.
- It would be interesting to apply the same method to rotating neutron stars, adding spin corrections, and see if the universality of the bound holds.
- The saturating EoS family could be used as a test case to compare with microwave pulsar observations, indirectly confirming the presence of a phase transition.
Why this may matter
This bound provides a model-independent filter that reduces uncertainties in neutron star radius estimates, which is important for understanding matter behavior at supra-nuclear densities and for interpreting gravitational wave data.
Nearby works in our archive
- The Ultimate Stiffness Limit of Neutron Starssimilarity 0.67
- How the Thickness of a Neutron Skin Reveals the Heart of Dead Starssimilarity 0.66
- The Unheard Hum of Neutron Starssimilarity 0.66
- A Stellar Recipe on Papersimilarity 0.65
- Stellar Interiors: A Clash of Equationssimilarity 0.64
- A Star with a Double Heartsimilarity 0.63
- Nebula Orchestra: Why Lines Hear Different Densitiessimilarity 0.63
The closest work in the archive sits at 0.67 by meaning — the usual level for our corpus.
Original abstract
as written by the authors arXiv:2606.25186
We investigate how the assumption of a common causal equation of state (EoS) correlates the radii of neutron stars at different masses and thereby reduces the uncertainties inferred from independent observations. We show that causality, anchored only to the chiral effective field theory ($$\chi$$EFT) EoS near saturation density, places a closed-form upper bound on the radius difference, $$R(2.0\,M_\odot)\le 1.16\,R(1.4\,M_\odot)-1.1\,$$km. The bound is saturated exactly by a one-parameter family of EoSs that we construct analytically. Imposing this prior-independent causal ceiling on the independent NICER posteriors of PSR J0437-4715, PSR J0614-3329, and PSR J0740+6620 retains only 7.5% of their joint product distribution and removes the large-radius tail of the PSR J0740+6620 posterior. Unlike full EoS-informed inferences, our construction cleanly isolates the consequences of the generic physical assumptions of a common causal EoS from those associated with a particular choice of EoS prior, providing a transparent benchmark for interpreting neutron-star observations.
Interestingly, if a neutron star had a radius forbidden by causality, then inside it, sound would have to travel faster than light, violating special relativity. Nature wisely avoids such paradoxes, and now we can measure that.
In the science fiction novel 'Dragon's Egg' by Robert Forward, intelligent life on the surface of a neutron star is described, where due to monstrous gravity, time runs a million times faster. Our causal constraints on the radii of such objects add physical rigor to such fantasies, showing what sizes are actually possible.
Key numbers
- sound speed in core: > 1/3 c^2
- fraction of models passing with causality imposed: 7.5%
- coefficient of linear radius relation: 1.16
- intercept in linear constraint: -1.1 km
- median radius of J0740+6620 after constraint: ~11.4 km