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.
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.
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.
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.
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.
The results will impact neutron star physics, nuclear physics, and gravitational-wave astrophysics by providing a new tool for studying superdense matter.
Apply the method to real data from future observations and refine the correlation with more realistic simulations.
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.
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
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.
Nearby works in our archive
The archive is sparse around this work: the closest review is only 0.60, against a usual 0.64. It looks like few have drilled nearby.
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.