🔭 astrophysics

Causality alone bounds the maximum radius difference between different-mass neutron stars

This paper demonstrates that the assumption of a common causal equation of state, anchored to chiral effective field theory, establishes a strict upper bound on the radius difference between neutron stars of different masses, thereby significantly constraining observational uncertainties and removing large-radius tails from independent NICER data without relying on specific equation-of-state priors.

Original authors: Aleksi Kurkela, Tuhin Malik

Published 2026-06-25

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Aleksi Kurkela, Tuhin Malik

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Idea: One Rulebook for All Neutron Stars

Imagine the universe has a giant library of "rulebooks" that explain how matter behaves under extreme pressure. These rulebooks are called Equations of State (EoS).

Scientists have been measuring neutron stars (the super-dense corpses of exploded stars) and trying to figure out which rulebook is the right one. They have measured two main things:

How big a "small" neutron star is (about 1.4 times the mass of our Sun).
How big a "heavy" neutron star is (about 2.0 times the mass of our Sun).

The problem is that these measurements come from different stars and are statistically independent. It's like trying to guess the height of a basketball player and a sumo wrestler by asking two different people who have never met them. Usually, you'd get a wide range of guesses.

The Paper's Insight:
The authors realized that all neutron stars must follow the same rulebook. Furthermore, there is one universal law that cannot be broken: Causality. In physics, this means nothing can travel faster than the speed of light. If a rulebook suggests that sound or pressure waves could move faster than light, that rulebook is fake and must be thrown away.

The authors asked: "If we only assume that all stars follow the same rulebook and that nothing breaks the speed of light, how much can we narrow down the possible sizes of these stars?"

The Analogy: The Elevator and the Roof

Think of a neutron star like an elevator shaft.

The Bottom (Low Density): We know exactly how the elevator works near the ground floor because we have good data (from a theory called Chiral Effective Field Theory).
The Top (High Density): As the elevator goes up into the heavy star, the pressure gets so intense that our current theories break down. We don't know exactly how the elevator behaves at the very top.
The Roof (Causality): There is a hard roof at the very top. No matter how you build the elevator, it cannot punch through the roof. That roof is the speed of light.

The authors built a mathematical model to find the highest possible ceiling (the "causal ceiling") for the elevator. They found that if you fix the size of the elevator at the 1.4-sun mark, there is a strict mathematical limit on how big the 2.0-sun elevator can be.

They discovered a simple formula for this limit:

The size of the heavy star cannot be more than 1.16 times the size of the light star, minus a small constant.

If a measurement suggests the heavy star is bigger than this limit, that measurement is physically impossible under the rules of causality.

The "Magic" Construction

To prove this limit exists, the authors didn't just guess. They built a "perfect" set of rulebooks (Equations of State) that push right up against the speed-of-light limit.

Imagine you are trying to stretch a rubber band as far as possible without snapping it.

They started with the known rules at the bottom (the light star).
They then stretched the rubber band as hard as physics allowed (making the matter as "stiff" as possible) until they hit the speed-of-light limit.
They found that the "stiffest" possible material behaves like a sudden jump (a phase transition) followed by a straight line where sound travels at the speed of light.

This specific construction creates the absolute maximum size for the heavy star for any given size of the light star.

What This Means for Real Data (The "Filter")

The authors took real data from NASA's NICER telescope, which has measured three specific neutron stars:

Two "light" ones (PSR J0437 and PSR J0614).
One "heavy" one (PSR J0740).

Previously, when scientists looked at the heavy star, their data had a "long tail" of possibilities. It suggested the star might be huge (up to 16 km wide), though it was unlikely.

The Result:
When the authors applied their "Causality Filter" to this data:

The "Impossible" Tail Got Cut Off: The huge, unlikely sizes for the heavy star were instantly removed because they would require physics to break the speed of light.
The Uncertainty Shrank: The range of possible sizes for the heavy star became much smaller.
The Survival Rate: Only 7.5% of the original possible combinations of sizes survived the filter. The other 92.5% were mathematically impossible if the stars share a common, causal rulebook.

The Hidden Clue: The Speed of Sound

By looking at the remaining valid data, the authors found a surprising clue about the inside of these stars.

They compared the real data against a "boring" scenario where the speed of sound inside the star stays constant at a low value (called the "conformal limit"). They found that the real data does not fit this boring scenario.

The Conclusion:
For the heavy star to exist at the sizes we observe, the speed of sound deep inside its core must be much faster than the "boring" limit. In fact, it must be very close to the speed of light. This proves that the matter inside these stars is incredibly stiff and resistant to being squished.

Summary

This paper didn't discover a new star or a new telescope. Instead, it acted like a mathematical sieve.

It took messy, independent measurements of neutron stars.
It applied the strict rule of "nothing travels faster than light."
It filtered out the impossible answers, leaving a much clearer, smaller range of what these stars can actually be.
It proved that the inside of these stars is made of material so tough that sound waves race through them at nearly the speed of light.

Technical Summary: Causality Bounds on Neutron Star Radius Differences

Problem Statement
Neutron stars are governed by a universal equation of state (EoS), $p(\varepsilon)$ , linking their masses, radii, and tidal deformabilities. While independent observations (e.g., from NICER and gravitational wave events) probe different mass regimes, they constrain the same underlying EoS. Consequently, the radii of neutron stars at different masses are correlated. However, existing inferences regarding the radius difference $\Delta R = R(2.0 M_\odot) - R(1.4 M_\odot)$ —a key probe of high-density stiffening and phase transitions—are typically tied to specific EoS models or priors. The central question addressed is the extent to which $\Delta R$ is constrained by causality alone, independent of any specific EoS parametrization, assuming only a common causal EoS anchored near saturation density.

Methodology
The authors employ a two-pronged approach to determine the "causal ceiling" for $R(2.0 M_\odot)$ given a fixed $R(1.4 M_\odot)$ :

Agnostic Numerical Construction: They utilize a "fractal-bridge" prior based on thermodynamic bounds. This method hierarchically samples the allowed EoS space between two anchors:
- Low-density anchor: Chiral effective field theory ( $\chi$ EFT) constraints between $0.5 n_s$ and $1.1 n_s$ (where $n_s \approx 0.16 \text{ fm}^{-3}$ ).
- High-density anchor: Perturbative QCD (pQCD) results at $\mu_H = 2.6 \text{ GeV}$ .
  The construction uses ten refinement layers to generate $\sim 8.6 \times 10^5$ EoSs, ensuring thermodynamic consistency and causality ( $c_s^2 \leq 1$ ) without diffusive filtering.
Analytic Characterization: They analytically identify the specific family of EoSs that saturate the causal bound.
Astrophysical Filtering: All generated EoSs are integrated via the Tolman–Oppenheimer–Volkoff (TOV) equations. EoSs failing to support a maximum mass $M_{\text{TOV}} > 1.97 M_\odot$ are discarded.

Key Contributions and Results

Closed-Form Causal Bound: The study derives a model-independent upper bound on the radius of a $2.0 M_\odot$ star as a function of the $1.4 M_\odot$ radius. The bound is saturated by a one-parameter family of "direct-causal saturators," characterized by a phase transition at the $\chi$ EFT endpoint followed by a maximally stiff segment ( $c_s^2 = 1$ ). The resulting linear relation is:
$R_{\text{max}}(2.0 M_\odot) \leq 1.16 R(1.4 M_\odot) - 1.1 \text{ km}$
(Coefficients vary slightly based on soft/stiff $\chi$ EFT endpoints, as detailed in Table I).
Limitations of Parametric Models: The authors demonstrate that widely used sound-speed parametrizations (e.g., piecewise linear $c_s^2(\varepsilon)$ ) do not fully explore the causal EoS space. The $c_s^2$ model falls short of the true causal ceiling by approximately $0.2\text{--}0.5$ km, depending on $R(1.4 M_\odot)$ , due to its limited flexibility compared to the fractal-bridge prior.
Observational Refinement: Applying this causal ceiling as a prior-independent filter to independent NICER posteriors for PSR J0437−4715, PSR J0614−3329, and PSR J0740+6620 yields significant constraints:
- Only 7.5% of the joint product distribution of the three pulsars satisfies the causal bound.
- The large-radius tail of the PSR J0740+6620 posterior (which previously extended to $\sim 16$ km) is truncated to $R \lesssim 12$ km.
- The low-mass pulsars shift slightly toward larger radii, as the constraint preferentially excludes combinations where low-mass stars have small radii.
Lower Bound on Sound Speed: By comparing the NICER data against a "hybrid" construction where the EoS is forced to have a conformal sound speed ( $c_s^2 = 1/3$ ) above the $\chi$ EFT endpoint, the authors show that such models cannot support the observed radii. This implies that the squared sound speed in the densest stellar cores must exceed the conformal value ( $c_s^2 > 1/3$ ) at densities above the $\chi$ EFT endpoint.

Significance
The paper establishes that the assumption of a common causal EoS, anchored only to $\chi$ EFT near saturation density, provides a robust, prior-independent framework for interpreting neutron star observations. Unlike full EoS inference, which inherits assumptions from specific model choices, this approach cleanly isolates the consequences of generic physical principles (causality and universality).

The derived causal ceiling serves as a "transparent benchmark" or "robust null hypothesis." It allows for the rapid assessment of observational data without requiring re-analysis of raw likelihoods or access to original data. Deviations from this bound could indicate non-standard physics (e.g., dark matter components or modified gravity), while consistency with the bound validates the standard causal EoS framework. The authors note that while the bounds depend mildly on current $\chi$ EFT uncertainties, they are insensitive to high-density pQCD constraints, making the results robust against current theoretical limitations in those regimes.