Long-Term Multidimensional Models of Core-Collapse Supernovae: Progress and Challenges

While recent advancements in 3D multidimensional simulations have significantly validated the neutrino-driven explosion mechanism and enabled comprehensive modeling of core-collapse supernovae from collapse to remnant, critical uncertainties regarding explosion criteria, progenitor effects, magnetic fields, neutrino physics, and the equation of state continue to challenge definitive predictions of stellar outcomes.

The Gist

Imagine a star as a massive, cosmic pressure cooker. For millions of years, it balances the crushing weight of its own gravity with the outward push of nuclear fusion, like a tightrope walker balancing on a high wire. But eventually, the fuel runs out. The wire snaps. The star collapses inward in a fraction of a second, and then—BOOM—it explodes as a supernova, scattering the elements of life across the universe.

For decades, scientists have been trying to figure out exactly how that "BOOM" happens. This paper, written by expert Hans-Thomas Janka, is a progress report on the latest and most advanced computer simulations of these explosions. It's like moving from watching a blurry, black-and-white movie of a car crash to watching a high-definition, 3D slow-motion replay with every piece of debris tracked.

Here is the breakdown of what this paper tells us, using simple analogies.

1. The Problem: The "Stalled" Engine

When the star's core collapses, it bounces back like a super-dense rubber ball. This creates a shockwave that tries to blow the star apart. But in the old, simple 1D models (which only looked at the explosion from the center outward), this shockwave would hit a wall of falling debris, get tired, and stop. The explosion would fizzle out, and the star would just collapse into a black hole.

The New Insight:
The paper explains that the universe isn't a straight line; it's a messy, 3D room. When scientists started simulating these explosions in 3D (adding depth and width), they found that the shockwave doesn't just sit there. It gets a "boost" from neutrinos—tiny, ghost-like particles that flood out of the new, super-hot core (called a Proto-Neutron Star).

Think of the shockwave as a stalled car. In 1D, it just sits there. In 3D, the neutrinos act like a turbocharger. They heat up the gas behind the shockwave, creating massive bubbles of hot air (like a pot of boiling water) that push the shockwave forward, reviving the engine and blowing the star apart.

2. The "Long-Term" Challenge: The Marathon, Not the Sprint

The biggest breakthrough in this paper is the timeframe.

• Old Simulations: Ran for only a few milliseconds (a blink of an eye). They could see the engine start, but they couldn't see if the car would actually finish the race.
• New Simulations: These run for seconds, tens of seconds, and even minutes.

Why does this matter?
Imagine trying to judge a marathon runner by watching them take their first step. You don't know if they'll trip, sprint, or walk. By running the simulation for a long time, scientists can see:

• The Final Energy: How big is the explosion? (Is it a firecracker or a nuclear bomb?)
• The "Kicks": When the star explodes, the leftover core (a neutron star) often gets kicked backward, like a cannonball recoiling from a cannon. The paper shows that this kick happens over many seconds as the explosion pushes unevenly.
• The Black Hole Question: Sometimes, the explosion fails, and the star collapses into a black hole. The paper explores how this happens and even shows scenarios where a black hole forms after a successful explosion, as leftover debris falls back in.

3. The "Turbulent Kitchen" (Nucleosynthesis)

Inside the explosion, new elements are being cooked.

• Old View: Scientists thought the gas flowed out smoothly, like water from a hose.
• New View (3D): It's more like a violent kitchen blender. Hot gas shoots out in plumes, while cold gas falls back in. They crash into each other, creating shockwaves and mixing everything chaotically.

This chaos is actually good news for us. It explains how the universe creates heavy elements like Titanium-44 (which helps us date supernova remnants) and Nickel-56 (which makes the explosion shine bright). The 3D simulations show that this "turbulent kitchen" produces the right amount of these elements to match what we see in the sky.

4. The Mystery of the "Ghost" Neutrinos

In 1987, we detected neutrinos from a supernova (SN 1987A). It was a historic moment. But the signal stopped a bit earlier than our best models predicted.

• The Puzzle: The models say the core should cool down and stop emitting neutrinos quickly. But the detectors saw a few late "ghosts" (neutrinos) arriving much later.
• The Paper's Theory: Maybe the core didn't cool down as fast as we thought. Or, perhaps, a massive amount of falling debris (fallback) hit the core later on, reigniting the neutrino engine for a few more seconds. It's like a campfire that you thought was out, but then a gust of wind blows embers back to life.

5. The "Code Wars" (Why Scientists Disagree)

The paper is honest about a major headache: Different computer codes give different answers.

• One team's supercomputer (using the Fornax code) says a star of a certain size will explode.
• Another team's supercomputer (using the Prometheus-Vertex code) says that same star will collapse into a black hole.

It's like two chefs following the same recipe but getting different results because one uses a gas stove and the other uses an electric one. The paper admits we don't know which computer model is "right" yet. This is the biggest challenge for the future: we need to compare these codes side-by-side to find the truth.

Summary: What Have We Learned?

1. 3D is King: You cannot understand supernovae with simple, flat models. The 3D nature of the explosion is crucial.
2. Time Matters: We need to watch the explosion for a long time to see the final result. The "engine" keeps running for seconds, not milliseconds.
3. The Kick is Real: Neutron stars get kicked like a recoil, and black holes can get kicked too if the explosion is messy.
4. We Are Getting Closer: The simulations now match the light, the elements, and the neutrino signals we see in the real universe much better than before.
5. The Work Isn't Done: We still need to figure out why different computers disagree and how to include even more complex physics (like magnetic fields and quantum effects).

The Bottom Line:
This paper is a celebration of how far we've come in understanding the death of stars. We have moved from guessing to simulating the entire life cycle of an explosion in 3D. While we still have some puzzles to solve (like why the computers disagree), we are finally seeing the full picture of how the universe recycles itself, turning dead stars into the building blocks of new worlds.

Technical Summary

1. Problem Statement

Core-collapse supernovae (CCSNe) represent the explosive death of massive stars, a process governed by the complex interplay of neutrino physics, hydrodynamics, nuclear physics, and general relativity. Despite decades of research, several fundamental questions remain unresolved:

• Explodability: Which massive stars successfully explode as supernovae, and which collapse directly into black holes (BHs)? Current models show significant discrepancies depending on the simulation code and input physics.
• Explosion Mechanism: While the neutrino-driven mechanism is the leading theory, its quantitative viability in 3D and its ability to reproduce observed explosion energies and remnant properties (neutron star kicks, spins, and nucleosynthesis yields) are still being tested.
• Long-Term Evolution: Previous simulations often stopped shortly after shock revival (within $\sim$1 second). However, the final explosion energy, the formation of compact remnants, and the synthesis of heavy elements depend on evolution over many seconds to minutes.
• Observational Consistency: Models must explain the neutrino signal from SN 1987A (specifically the late-time events), the properties of SN remnants (e.g., Cas A), and the observed diversity of supernova types.

2. Methodology

The paper reviews progress in self-consistent, multidimensional (2D and 3D) hydrodynamic simulations that solve the coupled neutrino-hydrodynamics equations from the pre-collapse stage through core bounce, shock revival, and long-term post-bounce evolution.

• Codes and Transport: The review focuses on state-of-the-art codes, particularly Prometheus-Vertex (and its successor Prometheus-Nemesis) and Fornax.
  • Neutrino Transport: Uses multi-group, velocity-dependent transport schemes. The Ray-by-Ray-plus (RbR+) approximation is commonly employed, solving 1D transport in each angular direction but including non-radial advection and pressure terms.
  • Long-Term Strategy: To overcome the extreme computational cost of full neutrino transport over long timescales ($>10$ s), the Nemesis code is used. It replaces the full 3D transport with a 1D PNS cooling model (using Mixing-Length Theory for convection) to provide neutrino luminosities and spectra, which are then applied to the 3D hydrodynamics. This reduces computational cost by a factor of 10–30.
• Physics Inputs:
  • Equation of State (EoS): Various nuclear EoSs (e.g., LS220, SFHo, DD2, APR) are tested to understand their impact on PNS contraction and explosion dynamics.
  • Progenitors: Simulations use progenitors from stellar evolution models, including those with 3D perturbations from convective oxygen-shell burning prior to collapse.
  • Microphysics: Includes many-body corrections to neutrino-nucleon interactions, weak magnetism, recoil, and muon formation.

3. Key Contributions and Results

A. Explosion Dynamics and Energies

• 3D vs. 2D: 3D simulations reveal that the "explodability" of a star is highly sensitive to the initial conditions. Unlike 2D, where axisymmetry can artificially amplify shock expansion (especially with RbR+ transport), 3D models show that large-scale pre-collapse perturbations (from convective shells) are often necessary to trigger explosions in 3D for certain progenitors.
• Energy Saturation: Explosion energies do not saturate immediately. In massive progenitors, the explosion energy continues to grow for several seconds to over 10 seconds. This is driven by a cycle of accretion downflows (bringing fresh fuel close to the PNS) and neutrino-heated outflows. The energy gain is sustained as long as the mass infall rate is high.
• Code Discrepancies: Different codes (e.g., Prometheus-Vertex vs. Fornax) yield different explosion outcomes for the same progenitor. Fornax tends to produce more energetic and earlier explosions than Prometheus-Vertex, highlighting the need for systematic code comparisons.

B. Neutron Star (NS) and Black Hole (BH) Formation

• Kicks and Spins:
  • Kicks: Newborn NSs receive natal kicks from two main mechanisms: asymmetric mass ejection (hydrodynamic kick) and anisotropic neutrino emission (neutrino kick). Hydrodynamic kicks dominate in massive progenitors and can reach $>1000$ km/s, growing over 10+ seconds. Neutrino kicks (driven by the LESA phenomenon) dominate in low-mass explosions but are generally smaller ($\sim$50–150 km/s).
  • Spins: Asymmetric accretion and SASI spiral modes can transfer significant angular momentum, potentially spinning up NSs to periods of 10 ms to seconds.
• Black Hole Formation Channels: The paper identifies four distinct channels for BH formation, ranging from "direct collapse" (failed SN) to "fallback SNe" where a successful explosion is followed by BH formation due to fallback accretion. The outcome depends critically on the nuclear EoS and the accretion rate.

C. Nucleosynthesis

• Neutrino-Driven Winds (NDWs) vs. Turbulent Outflows: In low-mass explosions, quasi-spherical NDWs occur. However, in massive progenitors, the inner ejecta are dominated by turbulent, neutrino-heated outflows interacting with accretion downflows.
• 44Ti Production: Long-term 3D simulations resolve the underproduction of $^{44}$Ti found in 1D models. The turbulent environment allows for multiple reheating episodes and high-entropy $\alpha$-rich freeze-out, producing $^{44}$Ti yields consistent with SN 1987A and Cassiopeia A observations.
• 56Ni: There are discrepancies between codes regarding $^{56}$Ni yields. Some models (e.g., Fornax) may overproduce $^{56}$Ni compared to observed Type IIP SNe, suggesting uncertainties in the explosion mechanism or fallback accretion.

D. Multi-Messenger Signals

• Neutrinos:
  • 3D signals are smoother than 2D signals (which suffer from artificial equatorial accretion sheets).
  • SN 1987A Puzzle: Standard PNS cooling models with current EoSs predict cooling times that are too short to explain the last few neutrino events detected by Kamiokande II (9–12 s). Possible explanations include: (1) weaker PNS convection than modeled, (2) a late episode of massive fallback accretion, or (3) a phase transition (e.g., hadron-to-quark) reheating the PNS.
• Gravitational Waves (GWs):
  • GWs are generated by asymmetric mass motions (SASI, convection) and anisotropic neutrino emission (GW memory).
  • 3D matter signals are generally weaker (by a factor of 10–20) than 2D predictions but still detectable for galactic events. The signal contains distinct frequency components (SASI at $\sim$100–200 Hz, PNS oscillations at $>1000$ Hz).

4. Significance and Future Directions

• Validation of Mechanism: The paper confirms that the neutrino-driven mechanism is viable and capable of reproducing observed explosion energies and remnant properties, provided that 3D effects and long-term evolution are included.
• Computational Milestone: The ability to simulate CCSNe for tens of seconds with self-consistent physics (or efficient approximations like Nemesis) is a major breakthrough, allowing for direct comparison with late-time observables.
• Remaining Challenges:
  • Code Consensus: The lack of agreement between different simulation codes on "explodability" requires a community-wide, controlled code comparison project.
  • Microphysics: The nuclear EoS and neutrino flavor conversion (fast collective oscillations) remain major uncertainties.
  • Progenitor Physics: Better integration of rotation, magnetic fields, and 3D pre-collapse convection in stellar evolution models is needed.
  • Multi-Messenger Test: Future galactic supernovae will provide a definitive test of these models through simultaneous detection of neutrinos, gravitational waves, and electromagnetic radiation.

In summary, this review highlights that while 3D long-term simulations have matured into a powerful tool for understanding CCSNe, the field is currently in a phase of refining microphysics and resolving discrepancies between different numerical approaches to fully unlock the secrets of stellar death.