⚛️ nuclear theory

Strongly interacting matter with criticality induced by modified excluded volume in core-collapse supernova simulations

This paper develops a novel modified excluded volume equation of state that incorporates a first-order phase transition with a critical point to simulate core-collapse supernovae, revealing a distinct, burst-like neutrino signature and gravitational wave modes that differ from those predicted by traditional hybrid hadron-quark models.

Original authors: Anil Kumar, Noshad Khosravi Largani, Stefan Typel, Pablo Cerdá-Durán, Alejandro Torres-Forné, Tobias Fischer

Published 2026-07-14

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Anil Kumar, Noshad Khosravi Largani, Stefan Typel, Pablo Cerdá-Durán, Alejandro Torres-Forné, Tobias Fischer

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, a cosmic giant weighing 40 times our Sun, reaching the end of its life. It's like a giant, overfilled balloon that suddenly loses its air. The core collapses inward at supersonic speeds, crushing everything until it hits a point where it can't get any tighter. Usually, this creates a "bounce," sending a shockwave out that might blow the star apart in a spectacular supernova. But sometimes, that shockwave stalls, and the star just collapses into a black hole.

For a long time, scientists thought that if the star's core got hot and dense enough, the normal matter inside (made of protons and neutrons) would suddenly melt into a new, exotic soup called "quark matter." They imagined this change happening like a sudden switch flipping: one moment it's solid, the next it's liquid, with a sharp, jarring jump in density. This "two-phase" idea was the standard recipe for explaining how some stars might explode.

But in this study, the team led by Anil Kumar tried a different recipe. Instead of a sharp switch, they used a "Modified Excluded Volume" (MEV) approach. Think of it like a crowded dance floor. In the old model, dancers (particles) suddenly stop being people and turn into something else entirely, creating a sudden gap. In the new MEV model, the dancers just get more crowded and start bumping into each other more, changing how they move and interact gradually, like a gas behaving strangely under pressure. This model mimics the behavior of real gases (like the ones in your soda can) that can have "critical points" where the rules of physics get a bit wobbly.

The Simulation Results: A Tale of Two Stars

The researchers ran super-computer simulations of these 40-solar-mass stars using their new "MEV" rules to see what would happen.

Scenario 1: The "Sharp Switch" (Gibbs Construction)
First, they tested what happens if they force the new MEV model to act like the old "sharp switch" idea. They built a barrier where the matter jumps from one state to another.

The Outcome: The star failed to explode.
What happened: The core started to turn into quark matter, but the jump in density was too small to trigger a massive explosion. The star's core just reorganized itself a little bit, then continued to collapse. It didn't bounce back hard enough. Instead of a supernova, the simulation ended with the star turning into a black hole about 2.9 seconds after the bounce.
The Signal: Because there was no big explosion, there was no "burst" of neutrinos (ghostly particles). The signal was just a quiet, steady hum that eventually faded as the star died.

Scenario 2: The "Wobbly Gas" (Van der Waals Behavior)
Next, they let the MEV model do what it naturally does: behave like a real gas with a "wobbly" region where pressure can actually drop as density increases (a bit like a spring that gets loose before it snaps tight).

The Outcome: The star exploded!
What happened: As the core collapsed, it hit this "wobbly" region. Instead of a sharp jump, the core collapsed supersonically, creating a second, powerful shockwave. This shockwave grew stronger, pushing the outer layers of the star away.
The Signal: This explosion sent out a massive burst of neutrinos. But here's the twist: this burst was much longer than the ones predicted by the old "sharp switch" models. Instead of a quick 1–2 millisecond flash, this one lasted about 5–10 milliseconds. It was a broad, slow-rising wave of energy.
The Catch: Even though it exploded, the explosion was a bit messy. A lot of the material that was thrown out fell back onto the star, and the final energy of the explosion dropped significantly over time. The simulation showed it was still evolving, so the final result isn't fully settled yet.

Listening to the Star's Heartbeat

The team also listened to the "gravitational waves"—ripples in space-time caused by the star's wiggling. They looked at the star's "heartbeat" (its oscillation modes).

In the failed explosion (Scenario 1), the heartbeat changed in a specific way that signaled a collapse.
In the successful explosion (Scenario 2), the heartbeat went wild and fluctuated wildly when the shockwave formed.
Interestingly, the standard formulas scientists use to predict these heartbeats based on the star's size and density stopped working once the quark matter appeared. The star's internal physics had changed so much that the old math couldn't describe the new rhythm.

What This Means

This study suggests that the way we model the transition from normal matter to quark matter matters a lot. If the transition is a sharp, two-step jump, the star might just fail to explode. But if the transition is a smoother, more complex change (like the "wobbly gas" in their model), it could trigger an explosion.

However, the authors are careful to say this is all based on simulations in a simplified, spherical world. They haven't proven this happens in real life yet. They also point out that their "wobbly" model has some weird physics (like imaginary sound speeds) that they had to ignore to keep the computer running. They need to test more versions of this model and run even more complex, 3D simulations to see if this "wobbly" explosion really works in the messy, real universe.

So, while the old "sharp switch" idea might be too simple to explain some explosions, this new "wobbly gas" idea offers a fresh, exciting possibility for how a dying star might go out with a bang.

Technical Summary: Strongly Interacting Matter with Criticality Induced by Modified Excluded Volume in Core-Collapse Supernova Simulations

Problem Statement
The explosion mechanism of core-collapse supernovae (CCSN) remains one of the most significant open questions in astrophysics, heavily dependent on the high-density Equation of State (EOS). A specific scenario posits that a sufficiently strong first-order phase transition from hadronic matter to deconfined quark matter can trigger a successful explosion. However, previous studies of this mechanism have relied on the "two-phase approach," where distinct hadronic and quark matter models are joined via Gibbs conditions. This approach introduces uncertainties regarding the density jump, the onset density, and the pressure slope between phases, often lacking a critical point at finite density and temperature. Furthermore, recent constraints from neutron star observations (masses $\approx 2 M_\odot$ , radii $\approx 9\text{--}13$ km) and gravitational wave events (GW170817) necessitate EOS models that are consistent with both nuclear physics and astrophysical data.

Methodology
The authors develop a novel class of multi-purpose EOS, termed DD2-MEV, based on the Relativistic Mean Field (RMF) framework with density-dependent meson-nucleon couplings (specifically the DD2 model). Instead of a traditional two-phase construction, the authors employ a Modified Excluded Volume (MEV) approach.

MEV Mechanism: The model replaces the geometric interpretation of excluded volume with a change in the effective number of degrees of freedom. A medium-dependent excluded-volume functional, $\Phi_N(T, n_B)$ , modulates the degeneracy factors of nucleons. This functional is designed to mimic Van der Waals behavior, featuring over- and under-critical regimes and a critical endpoint at high temperatures.
Thermodynamic Consistency: The implementation includes rearrangement terms arising from the density dependence of the couplings and the MEV functional to ensure thermodynamic consistency.
Comparative Models: The study compares two versions of the DD2-MEV EOS:
1. DD2-MEV (Gibbs): A version where the Van der Waals instability region is removed by enforcing Gibbs phase equilibrium conditions (constant pressure and chemical potentials).
2. DD2-MEV (Van der Waals): The continuous model retaining the negative pressure slope and instability regions characteristic of the Van der Waals equation.
Simulations: The authors perform general relativistic neutrino radiation hydrodynamics simulations in spherical symmetry using the AGILE-BOLTZTRAN code. They utilize a six-species Boltzmann transport scheme and a $40 M_\odot$ progenitor (s40a28).
Analysis: The study includes a detailed analysis of the neutrino signal and a gravitational wave (GW) mode analysis using the GREAT code (general relativistic astroseismology) to track fundamental ( $f$ ) and gravity ( $g$ ) modes.

Key Contributions and Results

Failure of the Gibbs Construction:
Simulations using the DD2-MEV (Gibbs) EOS, which enforces a standard phase transition construction, fail to produce a supernova explosion. The density jump between the hadronic and quark phases in this model is significantly smaller than in previous RDF hybrid models. Consequently, the protoneutron star (PNS) undergoes only a mild hydrodynamical reconfiguration rather than a supersonic collapse. The simulation proceeds toward black hole formation without shock revival.
Explosion via Continuous Van der Waals Behavior:
In contrast, simulations using the continuous DD2-MEV (Van der Waals) EOS result in a successful explosion. The phase transition induces a sharp rise in central density and temperature, leading to a supersonic PNS collapse. However, unlike previous models where a strong second shock forms immediately at the phase boundary, this model exhibits a delayed shock formation. A sub-sonic perturbation wave propagates outward, reaching the standing bounce shock after several milliseconds before turning into a strong shock wave.
- Ejecta: The explosion ejects a larger mass initially ( $\sim 0.2 M_\odot$ ) compared to previous reports, with a diagnostic energy exceeding 10 Bethe ( $10^{52}$ erg), though significant fallback reduces the final ejecta mass to $\sim 1 M_\odot$ .
Distinct Neutrino Signatures:
- Gibbs Case: No millisecond burst is observed. The neutrino signal shows only short-time fluctuations corresponding to the structural readjustment of the PNS, followed by canonical accretion behavior leading to black hole formation.
- Van der Waals Case: A burst-like neutrino signature is released, but it is substantially broader (lasting 5–10 ms) compared to the 1–2 ms bursts reported in previous studies using bag models or RDF EOS. This broadening is attributed to the delayed formation and propagation of the second shock. The burst is characterized by a two-stage rise in luminosity: a slow initial rise due to PNS contraction, followed by a sharp spike as the shock crosses the neutrinospheres.
Gravitational Wave Mode Analysis:
The study performs a mode analysis of the PNS oscillations ( $f$ and $2g_1$ modes).
- Breakdown of Standard Fits: The authors find that the standard fitting relations linking mode frequencies to PNS mean density ( $\bar{\rho}$ ) and compactness ( $\kappa$ ) break down after the appearance of quark matter degrees of freedom.
- Distinct Behaviors: In the failed (Gibbs) case, the $f$ -mode slope changes abruptly due to structural changes, while the $2g_1$ mode decreases rapidly. In the successful (Van der Waals) case, both modes exhibit large fluctuations during the phase transition and shock formation, indicating that their excitation is governed by physical phenomena distinct from simple density/compactness scaling.

Significance and Claims
The paper claims to overcome the caveats of the traditional two-phase approach by utilizing a continuous EOS that mimics a first-order phase transition through a modified excluded volume mechanism. The primary significance lies in demonstrating that the details of the phase transition construction are critical for the CCSN explosion mechanism.

The study suggests that the "density jump" is a crucial parameter; a small jump (as in the Gibbs construction of this specific MEV model) prevents explosion, whereas the continuous Van der Waals behavior allows for it.
It highlights that the neutrino burst duration and the GW mode evolution are sensitive probes capable of distinguishing between different phase transition scenarios (failed vs. successful, continuous vs. Gibbs construction).
The authors modestly conclude that while these results highlight large differences compared to previous QCD-driven explosion models, they rely on a specific parameterization. They emphasize the need for multi-dimensional simulations to confirm these findings and to compute direct GW signals, noting that the current analysis avoids spurious solutions by nulling the sound speed in unstable regions.

Technical Summary: Strongly Interacting Matter with Criticality Induced by Modified Excluded Volume in Core-Collapse Supernova Simulations

More like this