🔭 astrophysics

An Exploration of the Equation of State Dependence of Core-Collapse Supernova Explosion Outcomes and Signatures

This study utilizes 3D simulations to demonstrate that the choice of nuclear equation of state (comparing SFHo and DD2) significantly influences core-collapse supernova outcomes, including explosion energies, neutrino and gravitational-wave signatures, recoil kicks, and nucleosynthetic yields, thereby highlighting the need for broader multi-progenitor investigations to fully quantify these dependencies.

Original authors: Aleksandr Rusakov, Adam S. Burrows, Tianshu Wang, David Vartanyan

Published 2026-04-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Aleksandr Rusakov, Adam S. Burrows, Tianshu Wang, David Vartanyan

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 a massive star, about nine times the weight of our Sun, reaching the end of its life. It runs out of fuel, its core collapses under its own weight, and then—BOOM—it explodes as a supernova. This explosion is one of the most violent events in the universe, creating the heavy elements that make up planets and people.

But here's the mystery: What exactly happens inside that collapsing core?

Scientists have built super-computer simulations to watch these explosions in 3D, but there's a huge unknown variable: the Equation of State (EOS). Think of the EOS as the "rulebook" for how dense, hot nuclear matter behaves when you squeeze it. It's like asking: "If I squeeze a sponge, does it squish easily, or does it turn into a rock-hard diamond?"

In this paper, the authors (Rusakov, Burrows, and team) ran two different simulations using two different "rulebooks" (called SFHo and DD2) to see how the outcome of the explosion changes.

Here is the story of what they found, explained simply:

1. The Two Rulebooks: The Sponge vs. The Spring

The SFHo Rulebook: This one suggests that nuclear matter is a bit "softer." When you squeeze it, it compresses more easily.
The DD2 Rulebook: This one says nuclear matter is "stiffer." It resists being squeezed and acts more like a stiff spring.

2. The Explosion: A Faster, Stronger Kick vs. A Sluggish Push

When the star collapses, it forms a tiny, super-dense ball called a Protoneutron Star (PNS).

With the "Softer" Rulebook (SFHo): The core compresses tightly. This creates a very hot, high-pressure environment that acts like a powerful rocket booster. The explosion happens sooner and is stronger. The star kicks out its outer layers with more energy.
With the "Stiffer" Rulebook (DD2): The core stays a bit larger and less dense. It's like trying to launch a rocket with a weaker engine. The explosion is slower to start and weaker. The star doesn't push as hard.

Analogy: Imagine two people trying to pop a balloon.

SFHo is a person with strong, quick hands. They squeeze the balloon, and pop! It explodes instantly with a loud bang.
DD2 is a person with stiff, rubbery hands. They squeeze, but the balloon stretches out more before finally giving way. The pop is quieter and slower.

3. The Neutrino "Flashlight"

Supernovas shine with a special kind of light called neutrinos (ghostly particles that barely interact with anything).

Because the SFHo core is smaller and hotter, it shines a brighter, hotter flashlight. The neutrinos come out with more energy.
Because the DD2 core is larger and cooler, its flashlight is dimmer and cooler.
Why it matters: The explosion is driven by these neutrinos heating up the gas behind the shockwave. Since the DD2 flashlight is dimmer, it doesn't heat the gas as well, which is why the explosion is weaker.

4. The "Kick" and the "Shake"

The Kick: When a star explodes, the leftover neutron star often gets kicked away at high speed (like a cannonball firing). The SFHo model gave the star a faster kick (about 132 km/s), while the DD2 model gave it a slower kick (about 68 km/s). A weaker explosion just doesn't push the star as hard.
The Shake (Gravitational Waves): The explosion creates ripples in space-time called gravitational waves. The authors found that the "sound" of the shake is different for each rulebook.
- The SFHo shake has a higher pitch (frequency).
- The DD2 shake has a lower pitch.
- They also found a weird "gap" in the sound waves that happens at different times for each model. It's like listening to two different musical instruments play the same song; the notes are slightly off, telling you which instrument is being played.

5. The Cosmic Kitchen: Making New Elements

Supernovas are cosmic kitchens that cook up heavy elements (like gold, silver, and iron).

The SFHo explosion happened so fast and violently that it created a slightly different mix of ingredients. It made a bit more of the "heavy" elements (atomic weights between 60 and 90) and a slightly more "neutron-rich" soup.
The DD2 explosion was slower, resulting in a slightly different recipe.
The Catch: The authors realized that the speed of the explosion matters more for the recipe than the rulebook itself. If the explosion is fast, you get one set of elements; if it's slow, you get another.

The Big Takeaway

This paper is a preliminary study. It's like testing two different recipes for a cake with just one type of flour. The authors found that the "stiffness" of the nuclear matter changes everything:

How big the leftover neutron star is.
How hard the star explodes.
How fast the leftover star flies away.
What kind of "sound" (gravitational waves) it makes.
What elements are created.

The Future: The authors say, "We did this for one star. Now we need to do it for many different stars with many different rulebooks." They want to build a complete map of how the universe's "rulebook" determines whether a star becomes a neutron star, a black hole, or explodes with a bang.

In short: The universe is like a giant physics lab. By changing the "rules" of how matter behaves, the scientists showed that the entire drama of a star's death—from the size of the bang to the music of the ripples—changes completely.

1. Problem Statement

Core-collapse supernova (CCSN) theory has advanced significantly, with the neutrino-driven mechanism (aided by turbulence) now widely accepted as the primary explosion engine. However, significant theoretical uncertainties remain, particularly regarding the nuclear Equation of State (EOS). The EOS determines the stiffness of nuclear matter at and above nuclear densities, which dictates the structure, radius, and compressibility of the protoneutron star (PNS).

While previous studies have examined how the EOS affects "explodability" (whether a star explodes), there is a lack of comprehensive 3D simulations running to late times that compare final observables (explosion energy, nucleosynthetic yields, recoil kicks, and multi-messenger signatures) across different viable EOS models. This paper addresses that gap by comparing two distinct, empirically consistent EOS models: SFHo (Steiner et al. 2013) and DD2 (Typel et al. 2010).

2. Methodology

The authors utilized the state-of-the-art 3D simulation code Fornax to simulate the collapse of a 9 $M_\odot$ solar-metallicity progenitor (from Sukhbold et al. 2016).

EOS Models:
- SFHo: Softer at high densities, higher effective nucleon mass ( $m^* \approx 0.76$ ).
- DD2: Stiffer at high densities, lower effective nucleon mass ( $m^* \approx 0.56$ ).
- Key Difference: The DD2 EOS produces a more extended PNS with lower central densities compared to SFHo.
Simulation Setup:
- Geometry: 3D spherical coordinates ( $r \times \theta \times \phi$ ) with a "dendritic" grid ( $1024 \times 128 \times 256$ ) extending to 100,000 km.
- Physics: Multi-group neutrino transport (12 energy groups for $\nu_e, \bar{\nu}_e, \nu_\mu$ ) including energy redistribution via neutrino-electron and neutrino-nucleon scattering (Kompaneets method). General Relativity is treated approximately using the TOV equation.
- Duration: Simulations ran to late times (several seconds post-bounce) to allow observables to asymptote.
Comparison: The study contrasts the hydrodynamic evolution, neutrino emission, explosion energetics, recoil kicks, nucleosynthesis, and gravitational wave (GW) signatures between the two EOS models.

3. Key Contributions

Late-Time 3D Comparison: This is one of the first studies to compare two distinct EOS models in full 3D to the asymptotic state, moving beyond simple "explodability" metrics.
Diagnostic Signatures: The paper identifies specific, potentially observable differences in neutrino light curves, gravitational wave spectrograms, and recoil kicks that could serve as diagnostics for the nuclear EOS.
Convection Dynamics: It highlights how the EOS alters the timing and nature of PNS convection, which in turn modulates neutrino emission and GW signals.
Nucleosynthesis Nuance: It clarifies the interplay between explosion speed, neutron richness ( $Y_e$ ), and the production of weak r-process elements.

4. Key Results

A. Hydrodynamics and Explosion Dynamics

Explodability: Both models explode, but the SFHo model explodes slightly earlier and more vigorously.
Shock Velocity: The SFHo shock reaches 15,000–18,000 km/s, whereas the DD2 shock is slower (12,000–15,000 km/s).
PNS Structure: The DD2 EOS (stiffer) results in a larger-radius PNS with lower central densities ( $4.485 \times 10^{14}$ g cm $^{-3}$ ) compared to SFHo ( $5.474 \times 10^{14}$ g cm $^{-3}$ ).
Convection: The DD2 model exhibits an earlier onset of global PNS convection (penetrating the core within $\sim$ 1 second vs. $\sim$ 1.7 seconds for SFHo) due to lower optical depths and faster lepton loss.

B. Neutrino Signatures

Luminosity & Energy: The DD2 model produces lower mean neutrino energies and lower luminosities for $\nu_e$ and $\bar{\nu}_e$ during the early explosion phase. This is a direct consequence of the lower PNS temperatures and densities.
Heating Rates: Lower neutrino energies and luminosities in the DD2 model lead to reduced neutrino heating in the gain region, explaining the weaker explosion.
Diagnostic Feature: A temporary $\sim$ 10-15% enhancement in $\nu_e$ luminosity occurs in the DD2 model at $\sim$ 1 second post-bounce, coinciding with the onset of deep PNS convection.

C. Explosion Energy and Recoil Kicks

Explosion Energy: The SFHo model yields a significantly higher asymptotic explosion energy ( $\sim$ 0.123 B, where $1 \text{ B} = 10^{51}$ erg) compared to DD2 ( $\sim$ 0.08 B). This is roughly a 50% difference.
Recoil Kicks: The final neutron star kick speed is higher for SFHo ( $\sim$ 132 km/s) than for DD2 ( $\sim$ 68 km/s), correlating with the explosion energy and asymmetry.

D. Nucleosynthesis

Yields: Both models produce similar amounts of $^{56}\text{Ni}$ and $^{44}\text{Ti}$ .
Neutron Richness: The faster SFHo explosion leads to slightly more neutron-rich ejecta.
Weak r-process: The SFHo model produces more species with atomic weights $A=60-90$ and a stronger weak r-process signature. This suggests that faster explosions (often associated with lower compactness progenitors) are more favorable sites for the weak r-process, though the effect is modulated by initial seed perturbations.

E. Gravitational Waves (GW)

F-mode Frequency: The SFHo model (more compact PNS) exhibits a slightly higher f-mode frequency.
Power Gaps & Avoided Crossings: The spectrograms show distinct differences in "power gaps" (regions of low power) and avoided crossings between g-modes and f-modes.
- DD2: Power gap at $\sim$ 550 Hz.
- SFHo: Power gap at $\sim$ 1000 Hz.
Neutrino Memory: The DD2 model produces a weaker "neutrino memory" GW signal due to lower early-time neutrino luminosity and reduced quadrupolar asymmetry.

5. Significance and Future Directions

EOS Diagnostics: The study demonstrates that the nuclear EOS leaves distinct imprints on multi-messenger observables (GW frequencies, neutrino light curves, kick speeds). Future observations of galactic supernovae could potentially constrain the EOS.
Progenitor Dependence: The authors note that for a more massive progenitor ( $18.5 M_\odot$ ), the DD2 EOS fails to explode entirely, while SFHo succeeds. This implies the EOS critically influences the black hole vs. neutron star birth mass function.
Limitations: This is a preliminary study limited to a single progenitor mass. The authors emphasize the need for a comprehensive survey across a wide range of progenitor masses (9–80 $M_\odot$ ), metallicities, and EOS models to build a fully predictive theory of CCSN observables.

In conclusion, the paper establishes that the stiffness of the nuclear equation of state is a fundamental driver of supernova outcomes, affecting not just whether a star explodes, but the energy, speed, chemical yield, and multi-messenger signatures of the event.