Precision constraints on the neutron star equation of state with third-generation gravitational-wave observatories

This paper demonstrates that a network of third-generation gravitational-wave observatories, specifically the Cosmic Explorer and Einstein Telescope, will be able to constrain the neutron star radius to within approximately 75–200 meters (90% credibility) across a relevant mass range by analyzing the loudest events from an expected 300,000 binary neutron star mergers in a single year, representing a tenfold improvement over current constraints.

The Gist

The Big Picture: Squeezing the Universe's Toughest Stuff

Imagine you have a block of matter so dense that a single teaspoon of it would weigh as much as a mountain. This is what exists inside neutron stars—the collapsed cores of dead stars.

Scientists have a huge problem: we can't recreate this kind of density in a lab on Earth. It's like trying to understand how a diamond feels by only looking at a pile of sand. We don't know the "rules" (the Equation of State) that govern how this matter behaves. Does it act like a super-hard rock? Or does it squish like a soft sponge?

This paper is about how we are going to finally figure out those rules, not by building a bigger lab, but by listening to the universe with super-powerful new ears.

The New Ears: Cosmic Explorer and Einstein Telescope

Currently, we have gravitational wave detectors (like LIGO) that can hear the "chirp" of two neutron stars crashing into each other. But these current detectors are a bit like listening to a whisper from across a noisy room. They can tell us something, but the details are fuzzy.

The authors of this paper are looking ahead to Third-Generation Observatories: the Cosmic Explorer (in the US) and the Einstein Telescope (in Europe).

• The Analogy: If current detectors are a pair of standard hearing aids, these new machines are like super-hearing implants that can hear a pin drop from a mile away. They will be 10 times more sensitive.

The Experiment: A Year of Listening

The researchers ran a massive computer simulation. They asked: "If we turn on these new super-detectors for one year, what will we learn?"

1. The Crowd: In just one year, this network is expected to hear 300,000 neutron star collisions. That's a lot of noise!
2. The VIPs: Instead of trying to analyze all 300,000 events (which would take forever), the team decided to focus only on the 75 loudest, closest events.
   • The Analogy: Imagine trying to guess the average height of people in a stadium. You could measure everyone, but it's easier and almost just as accurate if you just measure the 75 people standing right in front of the camera. The loudest events give us the clearest signal.

The Discovery: Measuring the "Squishiness"

When two neutron stars spiral toward each other, they stretch and squeeze each other like taffy before they crash. This stretching is called tidal deformability.

• If the star is made of "stiff" matter, it resists stretching (like a rock).
• If it's made of "soft" matter, it stretches easily (like jelly).

By listening to the gravitational waves, we can measure exactly how much they squish. This tells us the size (radius) of the star.

The Results:

• Current Status: With our current detectors, we can guess the radius of a neutron star within about 2.8 kilometers. That's like trying to measure a person's height and being off by the length of a school bus.
• Future Status: With the new detectors, the authors predict we can pin the radius down to within 200 meters (and even 75 meters for the most common types of stars).
• The Analogy: We are going from guessing someone's height within the length of a bus, to guessing it within the length of a single car. That is a 10x improvement.

The "Sweet Spot"

The study found that you don't need to listen to all 300,000 events to get this amazing result.

• The Analogy: It's like trying to figure out the average temperature of a city. You don't need to check every single thermometer; you just need the first 20 or so readings from the most obvious spots.
• The data showed that after about 20 of the loudest events, the improvement starts to level off. The "loudest" events do the heavy lifting.

Why This Matters

Knowing the exact size and "squishiness" of neutron stars is like finding the missing piece of a puzzle for nuclear physics.

• It tells us how matter behaves at the most extreme densities in the universe.
• It helps us understand if the core of a neutron star is made of normal neutrons, or if it turns into a soup of exotic particles (like "strange" matter or free-floating quarks).

The Bottom Line

This paper is a "roadmap" for the future. It tells us that in just a few years, when these new telescopes come online, we will move from guessing what neutron stars are made of to knowing it with incredible precision. We will finally be able to read the "instruction manual" for the densest matter in the universe.

Technical Summary

1. Problem Statement

The internal composition of neutron stars (NS) involves matter at densities far exceeding those achievable in terrestrial laboratories. The behavior of this matter is described by the Equation of State (EoS), which relates pressure to energy density. While current gravitational wave (GW) observations (e.g., GW170817) and X-ray timing (e.g., NICER) have provided initial constraints on the NS EoS, they lack the precision to definitively distinguish between competing nuclear physics models.

The primary challenge is quantifying how much third-generation (3G) GW observatories (specifically Cosmic Explorer and the Einstein Telescope) will improve these constraints. Previous studies have struggled to model the full range of binary neutron star (BNS) masses due to the computational complexity of analyzing long, high signal-to-noise ratio (SNR) signals. This paper aims to determine the precision with which a 3G network can constrain the NS radius and pressure across a plausible mass range ($1 - 1.97 M_\odot$).

2. Methodology

The authors employed a hierarchical Bayesian inference framework using simulated data from a network of one Cosmic Explorer (CE) and one Einstein Telescope (ET).

• Data Simulation:

  • Event Selection: They simulated the 75 loudest BNS merger events expected from one year of coincident operation.
  • Population Model: Merger distances were sampled from a rate density model. Masses were drawn from a Gaussian approximation of the observed Galactic BNS mass distribution.
  • Physics Assumptions: Neutron stars were assumed to have zero spin (simplification based on low observed spins). The "true" underlying EoS used to generate the mock data was the SLy model.
  • Waveforms: Gravitational waveforms were generated using the IMRPhenomPv2 NRTidal approximant within the LALSuite software suite.

• Inference Pipeline:

  1. Parameter Estimation: For each of the 75 events, Bayesian inference was performed using the Bilby package and dynesty (dynamic nested sampling). The ROQ (Reduced Order Quadrature) likelihood was used to accelerate computation.
  2. Marginalization: The analysis focused on component masses ($m_1, m_2$) and tidal deformabilities ($\Lambda_1, \Lambda_2$). Other parameters were marginalized out to obtain marginal posteriors.
  3. Continuous Likelihood Construction: Since individual events yield discrete posterior samples, the authors used Kernel Density Estimation (KDE) with an Epanechnikov kernel to create a continuous representation of the marginal likelihood for each event.
  4. Hierarchical Inference: The total likelihood for the dataset was calculated by integrating the single-event likelihoods along the predicted $\Lambda(m)$ curve defined by a specific EoS.
  5. EoS Parametrization: The EoS was modeled using a four-parameter spectral decomposition of the adiabatic index $\Gamma(p)$ as a function of pressure. This allows the EoS to be uniquely determined by the coefficients $\Upsilon = \{\gamma_0, \gamma_1, \gamma_2, \gamma_3\}$.

3. Key Contributions

• Quantitative Precision Forecast: The study provides the first rigorous, hierarchical Bayesian forecast of EoS constraints using a realistic 3G network (CE + ET) over a full year of operation, specifically focusing on the "loudest" events.
• Demonstration of Diminishing Returns: The paper identifies that the constraint precision on the NS radius plateaus after approximately the 20 loudest events, suggesting that the highest SNR events dominate the information gain regarding the EoS.
• Methodological Advancement: It successfully implements a scalable hierarchical inference pipeline that integrates discrete parameter estimation posteriors (via KDE) with a spectral EoS parametrization, overcoming previous computational bottlenecks.

4. Key Results

• Radius Constraints:
  • Using the 75 loudest events, the network can constrain the NS radius ($R$) to within $\lesssim 200$ meters (90% credibility) across the mass range $1 - 1.97 M_\odot$.
  • In the mass range where the BNS distribution peaks ($1.4 - 1.6 M_\odot$), the constraint tightens to **$\lesssim 75$ meters**.
  • This represents an improvement of roughly 10 times over current constraints from LIGO-Virgo-KAGRA and NICER.
• Pressure Constraints:
  • The pressure ($p$) is constrained to within an average of $\sim 18\%$ in the energy density range $2 \times 10^{34} - 2 \times 10^{35} \text{ J m}^{-3}$.
  • At $\epsilon \approx 5.2 \times 10^{34} \text{ J m}^{-3}$, the constraint reaches $\sim 4\%$.
• Event Count Sensitivity:
  • The width of the credible interval for the radius shows little improvement beyond the first 20 loudest events.
  • While adding hundreds of thousands of weaker events (expected total $\approx 3 \times 10^5$ per year) might further improve constraints, the marginal gain per event is small. The authors note that the improvement curve fits both a decaying exponential (negligible gain) and a power law (factor of 2 gain), making their current estimates conservative.
• Model Recovery: The analysis successfully recovered the true SLy EoS used to generate the data, with the true curve falling well within the 1-sigma credible intervals.

5. Significance and Limitations

• Significance:

  • The results demonstrate that 3G observatories will transform NS physics from qualitative constraints to precision measurements, potentially ruling out or confirming specific nuclear physics models (e.g., phase transitions to quark matter).
  • The predicted precision ($\lesssim 2\%$ on radius at $1.4 M_\odot$) surpasses current combined constraints from GWs and X-ray observations ($\sim 16%$).
  • This precision will serve as a critical cross-check for future X-ray missions like STROBE-X.

• Limitations & Future Work:

  • Waveform Accuracy: The results rely on the accuracy of the IMRPhenomPv2 NRTidal approximant. Systematic biases from missing higher-order post-Newtonian terms or resonant modes (e.g., r-modes) could affect the inferred EoS.
  • Population Models: The analysis assumes a Galactic mass distribution and a specific merger rate model. The true extragalactic distribution (potentially broader, as suggested by GW190425) could alter the constraints.
  • EoS Complexity: The spectral decomposition assumes a smooth EoS. If the true EoS contains discontinuities (e.g., due to phase transitions), smooth parametric models might underestimate uncertainties.

In conclusion, this paper establishes that a network of Cosmic Explorer and the Einstein Telescope will provide unprecedented constraints on the neutron star equation of state, primarily driven by the highest signal-to-noise ratio events, effectively opening a new era of precision nuclear astrophysics.