Listening to the Long Ring: How Gravitational Waves Reveal the Equation of State of Neutron Stars ✛

Original: "Listening to the long ringdown: A novel way to pinpoint the EOS in neutron-star cores"
arXiv:2509.18665 · 2025-09-23 · CC BY 4.0 · 3 min
The long ring after a neutron star merger lets us directly measure the pressure in their cores.

Context

In neutron star cores, matter is compressed to densities several times that of an atomic nucleus. The equation of state of such matter remains one of the major unsolved problems in strong-interaction physics. Existing constraints come from observations of neutron star mergers (e.g., GW170817), mass and radius measurements, and theoretical calculations based on quantum chromodynamics. However, all of these pertain to densities up to about three times nuclear saturation, whereas in the centers of the most massive neutron stars, densities can reach five to six times saturation. This work joins a line of research using post-merger gravitational waves to probe these extreme conditions, but it proposes a new observable parameter that had not been considered before.

Methods

The authors built an ensemble of equations of state using Gaussian processes constrained by data from chiral effective field theory and perturbative quantum chromodynamics. Then, using principal component analysis, they selected six 'golden' equations of state that cover the 68% confidence region. For each, they performed numerical simulations of neutron star mergers with parameters consistent with GW170817: a fixed chirp mass of 1.18 solar masses and three mass ratios of the components (0.7, 0.85, 1.0). From the resulting gravitational-wave signals, they extracted the energy and angular momentum carried away by the waves, as well as the instantaneous frequency. In the long ring phase (from 1 to 15 ms after merger), they computed the slope of the energy versus angular momentum relation, which turned out to be constant.

Results

Key result: the slope dE_GW/dJ_GW in the long ring phase correlates with the pressure and density at the center of the most massive neutron star (the TOV point). For all six 'golden' equations of state and all three mass ratios, this correlation is described by a bilinear model. For example, as the pressure at the TOV point increases from 0.2 to 0.6 GeV/fm³, the slope changes from about 3.0 to 3.5 (in dimensionless units). Simulations of measurements with a ±4% uncertainty showed that combining the slope with the peak frequency f2 can reduce the uncertainty in the equation of state at high densities by a factor of several compared to using f2 alone.

Implications

This work offers a new observable parameter of the gravitational-wave signal after a neutron star merger that is directly linked to the equation of state at maximum densities. Unlike previous studies that focused on spectral peaks (e.g., f2), here they use the energy-to-angular-momentum ratio, which is easier to extract from the signal and less sensitive to simulation details. This paves the way for direct constraints on the equation of state in the density range of 4–6 times saturation, where currently no data exist.

Future development

With the advent of third-generation detectors such as the Einstein Telescope and Cosmic Explorer, we expect to detect many neutron star merger events with high signal-to-noise ratios. This will allow the proposed method to be applied to real data, yielding increasingly precise constraints on the equation of state. Furthermore, the correlation could be refined with more detailed simulations that include magnetic fields, neutrino emission, and other effects.

Impact

The results will impact neutron star physics, nuclear physics, and gravitational-wave astrophysics by providing a new tool for studying superdense matter.

Next steps

Apply the method to real data from future observations and refine the correlation with more realistic simulations.

Key open problems

This work directly addresses the unsolved problem of the equation of state of superdense matter: so far, there are no direct constraints at densities above 3 times saturation, and the proposed method could provide them.

The knowledge machine's view

This section is addressed to the author of the paper and their colleagues. Neighbouring works were found by meaning-based search across our archive, not by word overlap. We do not review: we show the field around the work.

What the work does

The authors studied the long ringdown phase following a neutron star merger, extracting from numerical simulations the slope of the gravitational-wave energy dependence on angular momentum. They showed that this slope correlates with the central pressure of the most massive neutron star, and that future third-generation detectors could use it to narrow down the uncertainty in the equation of state at supranuclear densities.

Where its strength is

Introducing a new observable, the slope dE_GW/dJ_GW, which is directly linked to the equation of state at the highest densities, and demonstrating its potential for future observations.

Where it could go next

Worth trying

  • Test whether the slope-pressure correlation persists when using a broader set of equations of state, including those generated by machine learning methods.
  • Assess how systematic errors in numerical simulations (e.g., grid resolution) affect the extracted slope.
  • Explore the possibility of combining the slope with other post-merger parameters to improve constraints on the equation of state.

Why this may matter

If the correlation holds, it would provide a new tool for probing the equation of state in a density regime inaccessible to other methods, which is important for understanding strong interactions.

The long ring lasts only a few milliseconds, but during that time the merger remnant loses energy equivalent to the mass of several Jupiters.

In Carl Sagan's novel 'Contact', gravitational waves from a neutron star merger serve as the key to deciphering an alien message; here, too, they carry information about the most extreme states of matter.

\frac{dE_{\rm GW}}{dJ_{\rm GW}} = \frac{\dot{E}_{\rm GW}}{\dot{J}_{\rm GW}} = \frac{\dot{h}_+^2 + \dot{h}_\times^2}{h_+ \dot{h}_\times - \dot{h}_+ h_\times}
Ratio of the energy loss rate to the angular momentum loss rate, expressed in terms of the gravitational wave amplitudes h+ and h×
f_{\rm GW} = \frac{1}{2\pi} \frac{h_+ \dot{h}_\times - \dot{h}_+ h_\times}{h_+^2 + h_\times^2}
Frequency at which the gravitational wave oscillates at a given moment
P({\rm EOS}, {\rm NSs} | {\rm data}) = \frac{P({\rm data} | {\rm EOS}, {\rm NSs}) \, P({\rm EOS}, {\rm NSs})}{P({\rm data})}
Posterior probability of the equation of state and neutron star parameters given the observational data

Key numbers

  • nuclear saturation density: 0.16 baryons/fm³
  • maximum density at the center of a neutron star: 5–6 times saturation
  • duration of the long ring: 1–15 ms
  • chirp mass in simulations: 1.18 solar masses
  • uncertainty of the modeled slope measurement: ±4%
Scientists
Albert EinsteinBarry BarishKip ThorneRichard FeynmanRussell HulseHenri Poincaré
Concepts
neutron star merger neutron star binary merger gw170817 gravitational wave echoes direct waves ringdown mass ratio gw190814 gravitational wave resonance kilonova overtones nuclear fusion neural network regression nanograv stellar remnant binary stellar evolution nuclear reaction mass gap core collapse event binning