Which Neutron Stars Reach the Stiffening Regime? Multimessenger Constraints on Core Sound Speed and Stellar-Mass Thresholds

By combining gravitational wave and X-ray observations, this study infers that neutron star cores likely exhibit sound-speed stiffening above $1/3$ of the speed of light, a regime typically reached by stars with masses around $1.6\,M_\odot$ and fully accessed only by the most massive pulsars near $2.1\,M_\odot$.

The Gist

Imagine the universe is filled with the most extreme "stuff" you can think of: neutron stars. These are the collapsed cores of dead stars, so dense that a single teaspoon of their material would weigh a billion tons on Earth.

For decades, physicists have been trying to figure out the "recipe" for this super-dense stuff. Specifically, they want to know: How does it react when you squeeze it?

In this paper, the authors act like cosmic detectives. They combine two types of clues to solve the mystery:

1. The "Squeeze" Clue: Data from a cosmic crash (GW170817) where two neutron stars smashed together, sending ripples through space-time.
2. The "Size" Clue: Precise measurements from the NICER telescope, which weighs and measures the size of three specific, heavy neutron stars.

Here is the simple breakdown of what they found, using some everyday analogies.

1. The Mystery of the "Stiffening"

Think of neutron star matter like a sponge.

• At first, when you squeeze a sponge, it's easy to compress.
• But if you keep squeezing, it eventually gets hard to push down. It "stiffens."

In physics, this "stiffening" is measured by the speed of sound inside the star. If the matter is very stiff, sound travels through it incredibly fast.

• The Big Question: Does the matter get so stiff that the speed of sound exceeds a universal speed limit (called the "conformal limit")?
• The Answer: The data suggests YES. Inside these stars, the matter gets incredibly stiff, much stiffer than we expected, right in the middle layers of the star.

2. The "Goldilocks" Zone of Mass

This is the paper's most important discovery. It's not just about what the stuff is made of, but which stars are actually showing us this stuff.

Imagine a ladder of neutron stars, sorted by weight:

• The Light Ones (1.4 Solar Masses): These are like a child trying to lift a heavy box. They are too light to reach the "stiff" part of the recipe. They only feel the easy, squishy bottom layers.
• The Heavy Ones (2.0+ Solar Masses): These are like strong weightlifters. They are heavy enough to crush the star down deep enough to hit the "stiff" zone.

The Finding:
The authors found that the "stiffening" starts happening in stars that are about 1.6 times the mass of our Sun.

• The famous heavy pulsar PSR J0740 (which is about 2.1 times the Sun's mass) is almost certainly inside this stiff zone. It has "cracked the code."
• However, even this heavy star hasn't quite reached the peak of the stiffness yet. It's like the weightlifter has lifted the heavy box, but hasn't quite reached the top of the shelf yet.

3. The "Smooth" vs. "Bumpy" Road

The authors tested different theories about how the matter behaves.

• The Smooth Road: Imagine the matter gets stiffer and stiffer in a smooth, predictable curve. The data likes this idea.
• The Bumpy Road: Imagine the matter gets stiff, then suddenly gets soft again (like hitting a pothole), before getting stiff again.
• The Result: The data strongly prefers the "Smooth Road." While we can't rule out the "Bumpy Road" entirely, the evidence points to a steady, strong stiffening as you go deeper into the star.

4. Why This Matters (The "Observational Roadmap")

Before this paper, scientists were guessing about the "recipe" of neutron stars in a very abstract way, talking about densities and numbers that no one can see.

This paper changes the game by giving us a target list.

• Old Way: "We think the stiff stuff happens at density X." (Hard to test).
• New Way: "We think the stiff stuff starts in stars that weigh between 1.6 and 2.1 Suns." (Easy to test!).

The Takeaway:
The universe has already handed us the perfect test subjects: the heavy neutron stars we can see right now (like PSR J0740). They are the "bridge" between the easy-to-understand physics and the extreme, super-stiff physics.

What's Next?
The authors are telling astronomers: "Stop looking at the light stars. Go measure the heavy ones (between 1.9 and 2.2 solar masses) even more precisely."

• If we find that these heavy stars are indeed behaving like the "stiff" model predicts, we will have solved the mystery of how matter behaves at its most extreme.
• If they don't, we'll know the "recipe" is different, perhaps involving a "bumpy" road or a sudden change in the material.

In a nutshell: We used cosmic crash data and telescope measurements to realize that heavy neutron stars are the key. They are the only ones heavy enough to crush matter into a super-stiff state, and by studying them, we are finally learning the true nature of the universe's densest material.

Technical Summary

1. Problem Statement

The paper addresses a fundamental question in dense matter physics: Do the cores of observed neutron stars (NS) reach densities where the speed of sound ($c_s$) exceeds the conformal limit ($c_s^2 > 1/3$), and if so, which specific observed stars probe this regime?

While perturbative QCD predicts $c_s^2 \to 1/3$ at asymptotically high densities, multimessenger analyses suggest that intermediate-density stiffening (where $c_s^2 > 1/3$) may occur within the cores of observed neutron stars. However, previous studies often focused on the probability of this phenomenon occurring at a specific density without clearly translating that inference into the stellar masses required to access those densities. The authors aim to bridge the gap between abstract density-space constraints and observable stellar-mass thresholds.

2. Methodology

The authors employ a multimessenger inference framework combining gravitational wave (GW) and X-ray timing data to constrain the Equation of State (EOS) of dense matter.

• Data Sources:
  • GW170817: Tidal deformability information used to constrain the EOS at lower-to-intermediate densities.
  • NICER Mass-Radius Posteriors: Data from three pulsars: PSR J0030+0451, PSR J0740+6620 (a massive ~2.07 $M_\odot$ pulsar), and PSR J0437−4715.
  • Mass Constraint: The EOS must support radio pulsars with masses $> 2 M_\odot$.
• EOS Parameterization:
  • The baseline model uses a smooth, five-parameter sound-speed profile ($c_s^2(n)$) defined by:
    • $n_{rise}$: Onset density where stiffening begins.
    • $\delta n$: Width of the stiffening transition.
    • $c_{s,peak}^2$: Peak sound speed.
    • $c_{s,high}^2$: High-density relaxation value.
  • The profile is constructed using hyperbolic tangent functions to ensure smoothness and thermodynamic consistency.
• Computational Approach:
  • Machine Learning Emulator: A neural network emulator trained on 7,256 exact Tolman–Oppenheimer–Volkoff (TOV) solutions accelerates the likelihood evaluation.
  • Validation: The main posterior results are validated using exact, posterior-resampled TOV solutions to ensure emulator errors do not bias the conclusions.
  • Likelihood Construction: Uses kernel-density estimates (KDE) directly from public posterior samples (avoiding Gaussian approximations) for both GW170817 (in the $M_c, \tilde{\Lambda}$ plane) and NICER sources (in the $M, R$ plane).
• Model Comparison: The baseline "smooth peaked" family is compared against three alternative families:
  1. Monotonic: Smooth stiffening without a post-peak drop.
  2. Piecewise: Six-knot interpolation allowing more flexibility.
  3. Transition-capable: Allows for a strong intermediate-density softening (dip) before restiffening.

3. Key Contributions

• Translation to Stellar Mass: The primary novelty is converting density-space inferences (e.g., "stiffening occurs at $3.5 n_{sat}$") into stellar-mass thresholds (e.g., "stars above $1.6 M_\odot$ enter the stiffening regime").
• Model-Averaged Probabilities: Instead of relying on a single EOS family, the authors provide an equal-prior averaged probability across four distinct EOS families to offer a more conservative and robust estimate of supra-conformal behavior.
• Differentiation of Onset vs. Peak: The study distinguishes between the probability of a star entering the stiffening regime (onset) versus traversing the peak of the sound speed profile.

4. Key Results

• Supra-Conformal Probability:
  • In the baseline smooth peaked family, the probability that $c_s^2 > 1/3$ at $3.5 n_{sat}$ is 85.4%.
  • When averaging over all four families (including transition-capable models), the conservative probability drops to 79.0%.
  • The data favor intermediate-density stiffening, but the strength of this statement is model-dependent; allowing for strong softening dips (transition-capable models) reduces the probability to 55%.
• Mass Thresholds:
  • Onset Density ($n_{rise}$): Corresponds to a threshold mass $M_{rise} \approx 1.59 M_\odot$ (90% CI: 0.72–2.14 $M_\odot$).
  • Peak Density: Corresponds to a threshold mass $M_{peak} \approx 2.08 M_\odot$ (90% CI: 1.47–2.46 $M_\odot$).
• Access Probabilities for Specific Stars:
  • PSR J0740+6620 ($2.07 M_\odot$): Reaches the onset of stiffening in 91% of posterior draws but reaches the peak region in only 46% of draws.
  • Canonical $1.4 M_\odot$ Stars: Reach the onset in only 34% of draws and the peak in just 2%.
• Implication: Current multimessenger data primarily constrain whether the most massive observed stars have entered the stiffening regime, rather than whether they have fully traversed its peak.

5. Significance and Future Outlook

• Reframing the EOS Problem: The paper argues that the most valuable future constraints are not generic additional observations, but targeted mass-radius measurements for stars in the $1.8 - 2.2 M_\odot$ range.
• Observational Roadmap:
  • If future NICER-like observations show that stars near $2 M_\odot$ consistently exceed the $M_{rise}$ threshold and begin approaching $M_{peak}$, the case for broad intermediate-density stiffening will be solidified.
  • If high-mass stars remain consistent with EOS families that delay onset or feature strong softening, the current smooth-family interpretation will weaken.
• Conclusion: The central result is not just a probability value, but an observational program. Current data suggest that the "entry" into the supra-conformal regime is accessible to the heaviest known pulsars (like J0740), but determining the full morphology (rise-peak-fall) of the sound speed requires resolving the mass-radius relations of stars near the upper end of the observed mass distribution.

In summary, this work shifts the focus of neutron star EOS research from abstract density constraints to stellar-mass thresholds, providing a clear, falsifiable target for next-generation multimessenger astronomy.