⚛️ general relativity

Universal relations applied to proto-neutron star generated gravitational waves from three-dimensional core collapse supernova simulations

This study evaluates the reliability of universal relations linking proto-neutron star oscillation frequencies to gravitational wave signals using three-dimensional core-collapse supernova simulations, finding that while these relations show promise, significant discrepancies necessitate caution in interpreting gravitational wave detections to infer proto-neutron star properties and highlight the need for further refinement of these models.

Original authors: R. Daniel Murphy, Anthony Mezzacappa, Colter J. Richardson, Pedro Marronetti, Eric J. Lentz, Ryan E. Landfield

Published 2026-07-28

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: R. Daniel Murphy, Anthony Mezzacappa, Colter J. Richardson, Pedro Marronetti, Eric J. Lentz, Ryan E. Landfield

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a cosmic symphony, where the loudest, most dramatic instruments are exploding stars. When a massive star runs out of fuel, it doesn't just fade away; it collapses in on itself with a violent crunch, creating a core collapse supernova. This event is so energetic that it ripples through the very fabric of space and time, sending out gravitational waves—invisible waves that stretch and squeeze the universe itself. For decades, scientists have been trying to "listen" to these ripples. But here's the tricky part: the sound isn't just a simple boom. It's a complex, shifting melody that changes as the newborn star, called a proto-neutron star, settles down.

Think of this newborn star like a giant, glowing bell that has just been struck. As it rings, it vibrates in specific patterns, much like a guitar string vibrating at a certain pitch. Scientists have discovered a fascinating idea: if you know how fast the bell is vibrating (the frequency of the gravitational waves), you might be able to figure out how heavy or dense the bell is, without ever seeing it. These are called "universal relations." They are like a secret codebook that claims to translate the music of the star directly into its physical properties, regardless of what kind of star it was before it exploded. The big question is: does this codebook actually work in the real, messy universe, or is it just a pretty theory that falls apart when things get complicated?

This paper is a reality check for that secret codebook. The authors, a team of astrophysicists, decided to test these universal relations using the most detailed computer simulations they could find. They took data from complex, three-dimensional models of exploding stars—simulations that mimic the chaotic, swirling physics of a real supernova—and asked: "If we pretend we are detecting these gravitational waves from Earth, will the codebook tell us the right story about the star?"

The results are a bit of a "yes, but..." story. The researchers found that some of the universal relations work surprisingly well for a short while, like a translator who gets the first few sentences of a conversation perfectly right. For example, one specific relation (based on the work of Rodriguez et al.) managed to track the star's density quite accurately for the first several hundred milliseconds after the explosion. However, as time went on, the codebook started to get confused. Many of the relations began to predict the wrong properties, or they only worked if you defined the "surface" of the star in a very specific way that didn't match the simulation.

The authors discovered that the "universal" nature of these relations is more fragile than hoped. They found that differences in how the simulations were built—like whether they used a simplified version of gravity or a more complex one, or whether the simulation was 2D or 3D—could throw the predictions off. In fact, when they tried to use the gravitational wave frequencies to predict the star's surface gravity, the codebook often gave answers that were way off the mark, especially in the later stages of the explosion.

The main takeaway is a warning to future astronomers: while these universal relations are a powerful tool, they shouldn't be used blindly. If we detect a gravitational wave from a supernova in the future, we can't just plug the numbers into these equations and expect a perfect answer. The paper suggests that we need to be very careful, understanding that the "music" of the star is influenced by many complex factors that these simple formulas might miss. It's not that the codebook is useless, but it needs a lot more editing and testing before we can trust it to tell us the true story of a dying star. The authors conclude that while we see some agreement between the predictions and the simulations, there is still a long way to go to make these relations truly reliable for the real universe.

Technical Summary: Universal Relations Applied to Proto-Neutron Star Generated Gravitational Waves from Three-Dimensional Core Collapse Supernova Simulations

Problem Statement
The detection of gravitational waves (GWs) from a core-collapse supernova (CCSN) offers a unique opportunity to probe the birth and early evolution of the proto-neutron star (PNS). The high-frequency component of the GW signal (starting $\sim$ 100 ms after core bounce) is dominated by PNS oscillations. Previous research has established "universal relations" that link the frequencies of these quasi-normal oscillation modes to global PNS properties, such as mean density or surface gravity. These relations are termed "universal" because they are parameterized solely by PNS properties, theoretically independent of progenitor characteristics (e.g., Zero Age Main Sequence mass, metallicity) or specific simulation physics (e.g., nuclear equation of state). However, the extent to which these relations, derived from various simulation codes and dimensionalities, hold true for high-fidelity three-dimensional (3D) CCSN simulations remains an open question. This work investigates the validity of applying these externally developed universal relations to GW data generated by the Chimera CCSN simulation code.

Methodology
The authors utilized GW emission data from four specific models generated by the Chimera code:

D-series (3D): Two 3D simulations of a 15 $M_\odot$ solar-metallicity progenitor (D15) and a 25 $M_\odot$ zero-metallicity progenitor (D25).
E-series (2D): Two 2D simulations of the 15 $M_\odot$ progenitor using two different equations of state (LSBCK and SFHo).

The study applied seven distinct sets of universal relations derived from seven different studies (Torres-Forné et al. 2021; Sotani and Sumiyoshi 2019; Sotani et al. 2021, 2025; Mori et al. 2023; Rodriguez et al. 2023; Bizouard et al. 2021). These relations generally take the form of a polynomial in the logarithm of a PNS property $x$ (e.g., $f = c_1 + c_2\ln(x) + c_3x + c_4x^2 + c_5x^3$ ), where $f$ is the oscillation frequency.

The analysis proceeded via two complementary methods:

Forward Prediction: Using time-evolving PNS mass ( $M$ ) and radius ( $R$ ) extracted directly from the Chimera simulations to calculate predicted oscillation frequencies via the universal relations. These predictions were overlaid onto GW spectrograms constructed from the simulation data to compare against the observed peak frequencies (the $g/f$ -mode feature, or $gfF$).
Inverse Prediction: Using the peak $gfF$ frequencies identified in the spectrograms as inputs to the universal relations to solve for the predicted PNS properties (mean density or surface gravity). These predicted properties were then compared against the true properties derived from the simulation data.

The spectrograms were constructed using Kaiser windowing with specific overlap parameters to resolve the high-frequency signal, and the $gfF$ was defined as the collection of frequencies with maximum power spectral density above 250 Hz.

Key Results
The application of universal relations to the Chimera 3D and 2D models yielded mixed results, indicating that agreement is not guaranteed across all models or timeframes:

Dimensionality and Gravity Mismatch: The universal relations derived from 1D or 2D simulations, or those utilizing different gravitational treatments (e.g., full General Relativity in the modal analysis vs. pseudo-Newtonian in the CCSN simulation), frequently failed to track the $gfF$ frequencies accurately, particularly at late times ( $>300$ ms).
Rodriguez et al. (2023) Performance: The relations from Rodriguez et al. (2023), which were derived using 3D CCSN data (albeit with pseudo-Newtonian gravity in the simulation and relativistic gravity in the modal analysis), showed the best agreement. Specifically, the $2p_1$ mode relation tracked the $gfF$ frequencies in the D15 model well for nearly the entire signal, deviating only slightly after $\sim$ 500 ms.
Surface Definition Sensitivity: A significant discrepancy was observed regarding the definition of the PNS surface. Relations derived using a $10^{10}$ g cm $^{-3}$ density contour often failed to track the mean density defined by a $10^{11}$ g cm $^{-3}$ contour, and vice versa. Notably, the $2f$ mode relation from Rodriguez et al. (2023) tracked the mean density defined by the $10^{11}$ g cm $^{-3}$ contour in the D15 model, despite being derived using the $10^{10}$ g cm $^{-3}$ contour.
Stochasticity in D25: In the D25 model, characterized by stronger accretion and more stochastic GW emission, the universal relations showed less initial agreement and diverged more rapidly than in the D15 model.
Inverse Analysis: When using observed frequencies to predict PNS properties, the relations often overestimated or underestimated the true mean density or surface gravity, particularly at late times. The "best-case" scenario (noise-free detection) still revealed significant deviations from the true simulation evolution.

Significance and Claims
The paper concludes that while universal relations offer a promising framework for asteroseismology of PNSs, caution must be exercised when applying them to real or simulated GW detections. The authors emphasize that:

No Unique Definition: There is no unique definition for the PNS surface boundary used in linear perturbation analyses, leading to inconsistencies between relations.
Gravity Treatment Discrepancies: Differences between the gravity treatment in CCSN simulations (often pseudo-Newtonian) and the linear modal analysis (often relativistic) introduce systematic frequency shifts (up to 15%) that degrade the "universality."
Dimensionality Limits: Most relations were derived from 1D or 2D models; their applicability to 3D simulations is not guaranteed.
Risk of Misinterpretation: A naive application of these relations could lead to incorrect conclusions about the specific oscillation mode generating the signal or the physical properties of the PNS.

The authors assert that while there is some agreement between asteroseismological predictions and simulation outcomes, the observed discrepancies highlight the need for further development of universal relations that account for these systematic differences, particularly regarding dimensionality and gravitational treatments, before they can be reliably used to interpret future CCSN GW detections.

More like this