On the Possibility of a Strong First-Order Phase Transition in Neutron Stars

By performing Bayesian inference on neutron star data from gravitational waves and X-ray observations alongside theoretical constraints from chiral effective field theory and perturbative QCD, this study finds evidence favoring a strong first-order phase transition in dense matter that likely occurs above the central density of the most massive neutron stars, thereby reconciling the need for a stiff equation of state with asymptotic softening.

The Gist

Imagine the universe is filled with a mysterious, super-dense "cosmic dough" found only inside neutron stars—the collapsed cores of dead stars. For decades, physicists have been trying to figure out exactly how this dough behaves when you squeeze it tighter and tighter.

This paper is like a high-stakes detective story where the authors try to solve a specific mystery: Does this cosmic dough suddenly change its texture in a violent, abrupt way (a "strong first-order phase transition"), or does it just slowly get denser and smoother?

Here is the breakdown of their investigation, using simple analogies:

1. The Mystery: The "Texture" of the Universe

Think of the matter inside a neutron star like a block of Jell-O.

• The "Smooth" Theory (NPT): Some scientists think that as you squeeze the Jell-O, it just gets harder and harder to compress, but it stays Jell-O the whole time. It's a smooth transition.
• The "Abrupt" Theory (FOPT): Others think that at a certain pressure, the Jell-O might suddenly snap into a completely different state—like turning instantly into a rock or a gas. In physics terms, this is a "strong first-order phase transition." The paper defines this as a moment where the material loses its ability to "bounce back" (sound speed drops to zero) for a short stretch of density.

2. The Clues: Listening to the Stars

The authors couldn't just go inside a neutron star to check. Instead, they acted like detectives gathering clues from two main sources:

• The "Squish" Test (Gravitational Waves): When two neutron stars crashed into each other (an event called GW170817), they sent ripples through space. How much the stars "squished" before crashing tells us how stiff or soft their internal dough is.
• The "Flashlight" Measurements (NICER): A space telescope named NICER took pictures of several pulsars (spinning neutron stars). By measuring their size and weight, the team got a better idea of how the dough behaves under pressure.
• The "Lab" Rules: They also used two sets of theoretical rules:
  • Low-Density Rules: Based on experiments with atomic nuclei (Chiral Effective Field Theory).
  • High-Density Rules: Based on math that describes how particles behave when squeezed to extreme limits (Perturbative QCD).

3. The Investigation: A Digital Simulation

The authors built a massive computer simulation using a method called "Bayesian inference." Think of this as running millions of different scenarios to see which ones fit the clues best.

• They created two groups of scenarios: one where the dough changes smoothly (No Phase Transition) and one where it snaps abruptly (Phase Transition).
• They fed all the real-world data (the crash waves and the star measurements) into the simulation to see which group of scenarios was more likely to be true.

4. The Verdict: The "Snap" is Likely, But Hidden

The results were surprising and specific:

• The "Snap" is Real: The data slightly favors the idea that the abrupt "snap" (the phase transition) does happen. It's not a smooth Jell-O all the way through.
• The "Snap" is Deep: Here is the twist. The transition doesn't happen in the outer layers of the star where we can easily see it. The data suggests the "snap" happens deep inside, in the very center of the heaviest neutron stars.
  • Analogy: Imagine a heavy metal ball. The outside is smooth and hard. The "snap" only happens if you crush the ball so hard that the very core turns into something else. Since our current observations only see the outside of the ball, we don't see the change directly.
• Why This Matters: This finding solves a puzzle. The "smooth" theory struggles to explain how neutron stars can be so heavy without collapsing, while the "abrupt" theory usually makes them too soft to hold up that weight. By placing the "snap" deep in the center (where it doesn't affect the star's outer shape much), the authors found a way to have a heavy star that still obeys the high-density rules of physics.

5. What This Means for the Future

The paper concludes that while we can't see this "snap" in the stars we currently observe, it likely exists just beyond our reach.

• The Twin Star Myth: The study found that this "snap" probably doesn't create "twin stars" (two stars with the same weight but different sizes), which some people thought might happen.
• The Next Clue: To actually see this transition, we need to look at the aftermath of neutron star crashes. When two stars merge, they briefly create a super-dense remnant that goes deeper than any stable star ever could. Future detectors listening to the "ringing" of these crashes might finally catch the sound of this cosmic dough snapping.

In short: The authors used star data to guess the recipe of the universe's densest matter. They found that the matter likely undergoes a sudden, dramatic change deep inside the heaviest stars, a secret that keeps the stars from collapsing while satisfying the laws of physics.

Technical Summary

Technical Summary: On the Possibility of a Strong First-Order Phase Transition in Neutron Stars

Problem Statement
The phase structure of strong interaction matter at densities several times the nuclear saturation density ($n_0 \approx 0.16 \text{ fm}^{-3}$) remains a central open problem in nuclear physics and astrophysics. While lattice QCD and heavy-ion experiments have established that the transition from hadronic matter to quark-gluon plasma is an analytic crossover at high temperature and low baryon chemical potential, the nature of the transition at low temperature and high baryon density is unknown due to the sign problem in first-principles lattice calculations. Specifically, it is unclear whether cold dense QCD matter undergoes a strong first-order phase transition (FOPT), characterized by a vanishing sound speed ($c_s^2 = 0$) over a finite density interval. Previous studies have generally not directly compared the FOPT hypothesis against the no-phase-transition (NPT) hypothesis using Bayesian model selection, and those that have often restricted the FOPT onset to densities already probed within stable neutron stars.

Methodology
The authors perform a Bayesian inference using a non-parametric Gaussian-process (GP) equation of state (EOS) for cold, $\beta$-equilibrated neutron-star matter. The methodology involves:

1. EOS Construction: The squared sound speed $c_s^2(n)$ is modeled via GP regression on an auxiliary variable $\phi(n) \equiv -\ln[1/c_s^2(n) - 1]$ to enforce causality ($0 \le c_s^2 \le 1$).
   • NPT Hypothesis: The GP yields a smooth $c_s^2(n)$ throughout the density range.
   • FOPT Hypothesis: The model samples two densities, $n_S$ and $n_E$, from a uniform distribution. The sound speed is set to zero ($c_s^2=0$) for $n_S \le n \le n_E$, with independent GP realizations on either side of the transition.
2. Constraints: The inference combines multiple data sources:
   • Observational: Tidal deformability from the GW170817 binary neutron-star merger and mass-radius measurements from NICER for PSR J0740+6620, PSR J0030+0451, PSR J0437−4715, and PSR J0614−3329.
   • Theoretical: Chiral effective field theory (ChEFT) constraints below $1.5 n_0$ and perturbative QCD (pQCD) constraints at high density.
3. Density Range: The GP is terminated at a density $n_L$. The primary analysis uses $n_L = 25 n_0$ to span the full window between ChEFT and pQCD regimes, with a secondary analysis at $n_L = 12 n_0$ for comparison.
4. Statistical Framework: A hierarchical Bayesian framework is used to compute the posterior distributions and the Bayes factor ($B_{NPT}^{FOPT}$) to compare the two hypotheses. The analysis distinguishes between "FOPT-in" (transition occurs inside the most massive stable neutron star, $n_S < n_c$) and "FOPT-out" ($n_S \ge n_c$).

Key Contributions and Results

• Evidence for FOPT: For $n_L = 25 n_0$, the data provide moderate evidence in favor of the FOPT hypothesis over the NPT hypothesis, with a Bayes factor $B_{NPT}^{FOPT} = 3.2$.
• Location of Transition: Within the FOPT hypothesis, the data strongly favor the "FOPT-out" scenario ($B_{in}^{out} = 5.1$), indicating that the onset of the phase transition ($n_S$) most likely lies above the central density ($n_c$) of the most massive neutron star.
  • The median onset density is inferred to be $n_S \approx 7.35 n_0$ (with a 68% credible interval of $5.30$–$11.85 n_0$), while the central density of the maximum-mass star is $n_c \approx 6.13 n_0$.
  • When the analysis is restricted to $n_L = 12 n_0$, the evidence for FOPT weakens to anecdotal levels ($B_{NPT}^{FOPT} = 2.8$), but the preference for the transition occurring above $n_c$ persists ($B_{in}^{out} = 2.6$).
• Stellar Observables:
  • Canonical Masses: Because the transition is inferred to occur above the central density of stable stars, the canonical-mass observables (radius $R_{1.4}$ and tidal deformability $\Lambda_{1.4}$) are statistically indistinguishable between the FOPT and NPT hypotheses at the 68% credible level.
  • Maximum Mass: The FOPT hypothesis allows for a modest increase in the maximum non-rotating neutron star mass ($M_{TOV}$). For $n_L = 25 n_0$, $M_{TOV}$ shifts from $2.07^{+0.10}_{-0.08} M_\odot$ (NPT) to $2.15^{+0.13}_{-0.11} M_\odot$ (FOPT).
  • Twin Stars: The posterior probability of a "twin-star" solution (disconnected branches) is found to be $\le 0.1\%$, strongly disfavoring this signature.
• Sound Speed and Conformality: Under the FOPT hypothesis, the sound speed $c_s^2$ exhibits a peak near $4 n_0$, drops to zero across the transition plateau, and rebounds to a smaller peak before approaching the pQCD limit. The trace anomaly $\Delta$ indicates that matter inside neutron stars remains non-conformal and strongly coupled, even under the FOPT scenario.

Significance and Claims
The paper claims that the inferred FOPT, occurring predominantly above the central density of the most massive stable neutron star, offers a natural resolution to a physical tension: the need for a stiff EOS to support $\sim 2 M_\odot$ neutron stars versus the softening favored by pQCD at asymptotically high densities. By placing the softening (the $c_s^2=0$ plateau) outside the stable stellar interior, the model maintains the necessary stiffness for massive stars while satisfying high-density pQCD constraints.

The authors conclude that current stable neutron star observations cannot directly detect this specific type of FOPT. Instead, they suggest that direct tests should target systems that briefly access densities beyond $n_c$, such as the transient remnants of binary neutron-star mergers. Future gravitational-wave observations of post-merger signals by next-generation detectors are identified as the natural candidate for testing this conclusion. The paper does not claim to definitively prove the existence of an FOPT but presents it as a hypothesis moderately favored by current data when tested over a broad density range.