Astrophysical constraints on the cold equation of state of the strongly interacting matter

This study utilizes astrophysical observations, including massive pulsar measurements, NICER data, and GW170817 tidal deformability constraints alongside perturbative QCD calculations, to significantly restrict the admissible parameter space for the equation of state of cold, dense strongly interacting matter.

The Gist

Imagine the universe as a giant kitchen with a very specific, impossible recipe: cold, dense matter.

On Earth, we can't cook this dish. Our most powerful particle colliders (like the LHC) are like high-heat ovens; they can smash atoms together, but they create a hot, chaotic soup that doesn't tell us much about what happens when matter is squeezed cold and tight.

The only place in the universe where this "cold, dense" recipe actually exists is inside Neutron Stars. These are the cosmic leftovers of massive stars that have collapsed. They are so heavy that a teaspoon of their material would weigh a billion tons on Earth. Because they are so dense, they act as the universe's only natural laboratory for studying how matter behaves under extreme pressure.

The Mystery: The "Equation of State"

Physicists want to know the "Equation of State" (EOS). Think of the EOS as the instruction manual for this dense matter. It tells us: If you squeeze this matter harder, how much does it resist? Does it get squishy, or does it turn into something harder than diamond?

The problem is, we don't have the manual. We have to guess the rules by looking at the Neutron Stars and seeing how they behave.

The Detective Work: Using Clues to Narrow the Guess

The authors of this paper acted like detectives trying to solve a mystery. They started with a huge library of 10,000 possible instruction manuals (theories about how matter behaves). Most of these manuals were just guesses based on math and physics principles.

Then, they used real-world clues from space to cross out the manuals that didn't fit. Here are the clues they used:

1. The "Heavyweight" Clue (Mass):
   We know there is a Neutron Star called a "Black Widow" pulsar that is incredibly heavy (about 2.22 times the mass of our Sun).

   • The Analogy: Imagine you have a stack of 10,000 different bridges. You know for a fact that a truck weighing 2.22 tons drove across one of them without it collapsing. Any bridge design that would have crumbled under that weight is immediately thrown in the trash.
   • Result: This single clue knocked out about 80% of the possible manuals.

2. The "Speed Limit" Clue (pQCD):
   At the very center of a Neutron Star, matter is so dense that the rules of physics change, and we can use a specific type of math (perturbative QCD) to predict what happens.

   • The Analogy: It's like knowing that no matter how you design a car, it cannot legally drive faster than the speed of light. If a bridge design implies the car would break the speed of light, it's invalid.
   • Result: This ruled out a few more manuals that were physically impossible.

3. The "Squishiness" Clue (Tidal Deformability):
   When two Neutron Stars crash into each other (like the event GW170817), they stretch each other like taffy before merging. This "stretchiness" is called tidal deformability.

   • The Analogy: Imagine two marshmallows colliding. If they are very stiff, they barely change shape. If they are soft, they squish out a lot. The gravitational waves from the crash tell us exactly how much they squished.
   • Result: This was the biggest filter. It turned out that most of the remaining manuals predicted Neutron Stars that were either too stiff or too soft compared to what we saw in the crash. This clue alone reduced the list of valid manuals to less than 2%.

4. The "Size" Clue (NICER):
   The NICER telescope on the International Space Station takes X-ray pictures of Neutron Stars to measure their size (radius).

   • The Analogy: This is like measuring the circumference of the marshmallow.
   • Result: While helpful, the measurements from NICER still have a bit of "fuzziness" (uncertainty). They helped narrow the list, but they weren't as strict as the "squishiness" clue.

What Did They Find?

After applying all these filters, the authors found that the "instruction manual" for dense matter is much more specific than we thought.

• The "Sweet Spot": The matter inside these stars seems to undergo a transition. It starts as normal atomic matter (hadrons) and then morphs into a soup of quarks (the building blocks of protons and neutrons).
• The Transition: This change doesn't happen instantly like a light switch (a sharp jump); it happens gradually, like a smooth fade. The authors found this transition likely happens at a density about 4.8 times the density of a normal atomic nucleus.
• The Size: The valid manuals suggest Neutron Stars are generally quite large (around 12–13 km in radius) and not as small as some other theories suggested.

The "What If" Scenarios

The authors also tested two wild cards:

1. The "Tiny" Star: There is a candidate object that might be a very light Neutron Star. If this is real, it would force the rules to change even more. However, the authors note this object is controversial and might not even be a Neutron Star.
2. The "Gap" Star: There was a mysterious object detected in a crash (GW190814) that is heavier than any known Neutron Star but lighter than a Black Hole. If this object is a Neutron Star, it would be a massive constraint, forcing the "instruction manual" to be very stiff to support that weight.

The Bottom Line

The paper concludes that Neutron Star observations are the ultimate filter. While we have many theories about how matter works, the universe is very picky. The combination of the heaviest known stars and the "squishiness" observed in crashing stars has narrowed down the possibilities significantly.

Currently, the most restrictive clues are the mass of the heaviest stars and the tidal deformability from collisions. The "size" measurements from telescopes are useful but still a bit too fuzzy to be the deciding factor. The authors are left with a specific set of rules that matter must follow, but they admit there is still work to do to understand exactly why the matter behaves this way.

Technical Summary

1. Problem Statement

The fundamental properties of cold, dense strongly interacting matter remain one of the central unsolved problems in nuclear physics. While Quantum Chromodynamics (QCD) is well-understood at low densities (via nuclear physics) and at asymptotically high densities (via perturbative QCD), the intermediate regime—where the critical end point is expected to reside—is inaccessible to terrestrial experiments. Heavy-ion collisions probe this region but at high temperatures ($\sim 100$ MeV), making it difficult to isolate the cold Equation of State (EOS).

Neutron stars (NSs) serve as the only natural laboratories for studying matter at densities several times higher than nuclear saturation density ($\rho_0$) at near-zero temperatures. However, the EOS of dense matter is not uniquely determined by current observations. The paper aims to quantify how various astrophysical observables constrain the EOS, specifically investigating the transition from hadronic to quark matter and the resulting constraints on NS mass, radius, and tidal deformability.

2. Methodology

2.1 EOS Parametrization

The authors construct a unified EOS by connecting three distinct density regimes:

1. Low Density (Crust): Modeled using the unified crust EOS of Douchin and Haensel.
2. Hadronic Phase ($\rho \lesssim \rho_{BL}$): Modeled using the SFHo relativistic mean-field EOS (based on nucleons and $\omega, \rho, \sigma$ mesons). The authors also test robustness using the DD2 EOS.
3. Quark Phase ($\rho \gtrsim \rho_{BU}$): Modeled using a quark-meson model based on $U(3) \times U(3)$ chiral symmetry (eLSM), which includes constituent quarks and meson nonets.
4. Asymptotic High Density: Constrained by perturbative QCD (pQCD) calculations at $\mu \approx 2.6$ GeV ($\rho \approx 40\rho_0$).

Interpolation: The transition between the hadronic and quark phases is handled via a fifth-order polynomial interpolation of the energy density $\varepsilon(\rho_B)$ in an intermediate density region defined by a central density $\bar{\rho}$ and width $\Gamma$. This ensures continuity of pressure and the speed of sound ($c_s$).

Free Parameters: The analysis varies three key parameters:

• $g_V$: Vector coupling constant between vector mesons and constituent quarks ($0 \le g_V \le 10$).
• $\bar{\rho}$: Central density of the transition region ($2\rho_0 \le \bar{\rho} \le 5\rho_0$).
• $\Gamma$: Width of the transition region ($\rho_0 \le \Gamma \le 4\rho_0$).

2.2 Constraints and Data

The study employs Bayesian inference to constrain the parameter space using a dataset of approximately 10,000 randomly sampled EOSs. Constraints are applied as either strict conditions (hard cuts) or probabilistic likelihoods:

• Strict Conditions:

  • Causality & Stability: $c_s < 1$ and thermodynamic stability.
  • pQCD Connection: The EOS must connect causally to the pQCD reference point at $\mu = 2.6$ GeV.
  • Maximum Mass: The EOS must support a non-rotating (TOV) star with $M_{TOV} \ge 2.22 M_\odot$ (based on population analysis including PSR J0952–0607).

• Probabilistic Constraints (Likelihoods):

  • NICER: Mass-radius measurements for five pulsars (PSR J0030+0451, J0740+6620, J0614-3329, J1231-1411, J0437-4715).
  • GW170817: Tidal deformability constraints. The study uses a conservative bound $\tilde{\Lambda} < 720$ and a stricter range $70 < \Lambda_{1.4} < 580$.
  • Hypothetical Scenarios:
    • HESS J1731–347: A very light compact object ($M \approx 0.77 M_\odot$).
    • GW190814: A "mass-gap" object ($M \approx 2.59 M_\odot$) interpreted as a neutron star.

3. Key Contributions

1. Systematic Bayesian Analysis: The paper provides a rigorous statistical framework to map the allowed EOS manifold by varying the transition density and width between hadronic and quark matter, rather than assuming a fixed transition model.
2. Quantification of Constraint Strength: It explicitly ranks the constraining power of different astrophysical observations, demonstrating that tidal deformability and maximum mass are the most restrictive factors, while current NICER radius measurements (due to large uncertainties) have a lesser impact on narrowing the EOS manifold.
3. Model Robustness: The authors validate their results by comparing the SFHo-eLSM model against alternative constructions (DD2-eLSM, SFHo-MHT, and a first-order phase transition model BBKF), showing that the preference for specific transition parameters is robust across different hadronic models.
4. Transition Characteristics: The study identifies that the data favors a broad crossover transition rather than a sharp first-order phase transition.

4. Key Results

• EOS Manifold Reduction: Starting from ~9,261 initial EOSs, the application of strict causality and stability reduces this to ~1,993. The minimum mass constraint ($M_{TOV} \ge 2.22 M_\odot$) is the most powerful single filter, reducing the viable set to ~600.
• Tidal Deformability: The GW170817 tidal deformability constraint is the second most restrictive. Applying the conservative bound $\tilde{\Lambda} < 720$ reduces the allowed EOSs to <2% of the original set. The stricter bound ($70 < \Lambda_{1.4} < 580$) further shrinks the viable space.
• Preferred Parameters: The most probable parameter set found is:
  • Transition density: $\bar{\rho} \approx 4.6 - 5.0 \rho_0$
  • Transition width: $\Gamma \approx 2.8 - 3.0 \rho_0$
  • Vector coupling: $g_V \approx 4.0 - 6.0$
• Radius Predictions: The model predicts radii for a $1.4 M_\odot$ neutron star in the range of 12.5–12.8 km. This creates a tension with the NICER measurement for PSR J0614–3329, which suggests a smaller radius ($\sim 10.3$ km), though the large uncertainties in NICER data prevent this from being a definitive exclusion.
• Mass Gap Implications: If the GW190814 object ($2.59 M_\odot$) is a neutron star, the EOS constraints become extremely severe, effectively ruling out the current model's preferred parameter space.
• HESS Object: The existence of a very light NS ($0.77 M_\odot$) in HESS J1731–347 has negligible impact on the EOS constraints within this framework.

5. Significance

This work demonstrates that astrophysical observations are now the primary driver for constraining the cold EOS of strongly interacting matter. The study highlights a critical tension: current models that satisfy the high-mass constraint ($>2.2 M_\odot$) and pQCD limits tend to predict larger radii and higher tidal deformabilities than some specific NICER measurements suggest.

The finding that a broad crossover transition ($\Gamma \approx 3\rho_0$) is favored over a sharp first-order transition or a narrow crossover suggests that the phase transition from hadronic to quark matter in neutron star cores is likely a smooth crossover occurring at high densities ($\sim 5\rho_0$). This provides a crucial benchmark for future theoretical developments in QCD and guides the interpretation of upcoming multi-messenger data from gravitational wave detectors and X-ray telescopes.