🔭 astrophysics

Resolution Dependence in Magnetohydrodynamic Simulations of Neutrino-Driven Core-Collapse Supernovae

This study demonstrates that while shock revival times in neutrino-driven core-collapse supernovae are largely independent of resolution and initial magnetic field strength, higher resolutions and stronger initial magnetic fields lead to more efficient small-scale dynamo amplification, resulting in higher explosion energies and modified proto-neutron star properties that compensate for reduced neutrino heating.

Original authors: Vishnu Varma, Bernhard Müller

Published 2026-04-02

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Vishnu Varma, Bernhard Müller

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 13 times heavier than our Sun, reaching the end of its life. It runs out of fuel, its core collapses under its own weight, and it bounces back in a spectacular explosion called a supernova. For decades, scientists have tried to understand exactly how this happens. The leading theory is that invisible ghost particles called neutrinos act like a blowtorch, heating the material behind the shockwave and pushing the star apart.

However, there's a catch: stars also have magnetic fields, like giant, invisible magnets. For a long time, scientists mostly ignored these fields in their computer simulations, or assumed the stars weren't spinning fast enough for the fields to matter.

This paper is a "stress test" for our computer models. The authors asked: "Does it matter how detailed our computer simulation is, and how strong the star's magnetic field is?"

To find out, they ran the same explosion scenario six different times, changing two things:

The Resolution: How many "pixels" they used to draw the star. Some models were blurry (low resolution), and some were ultra-sharp (high resolution).
The Magnetism: Some stars started with a weak magnetic field (like a fridge magnet), and others started with a super-strong one (like a magnetar).

Here is what they discovered, explained through simple analogies:

1. The "Pixel" Problem (Resolution)

Think of the simulation like a video game. If you play on "Low Graphics," the water looks like a flat blue sheet. If you play on "Ultra," you see individual droplets and splashes.

The Finding: When the scientists used "Low Graphics" (low resolution), the computer smoothed out the chaos. It missed the tiny, swirling eddies of gas.
The Result: In the high-resolution models, the computer could see the turbulence clearly. This turbulence didn't just sit there; it acted like a generator, turning the chaotic motion of the gas into magnetic energy.
The Analogy: Imagine rubbing your hands together. If you do it slowly (low resolution), you don't get much heat. If you rub them furiously and quickly (high resolution), you generate a lot of heat. In the star, the "rubbing" (turbulence) generated strong magnetic fields that helped push the explosion harder.

2. The Magnetic "Sidekick"

The authors expected the strong magnetic fields to be the main hero, shooting out jets of energy to blow the star apart (like a magnetar-driven explosion).

The Finding: Surprisingly, the magnetic fields did not start the explosion. The neutrinos (the ghost particles) were still the ones doing the heavy lifting to revive the shockwave.
The Twist: However, once the explosion started, the strong magnetic fields acted like a turbocharger. They didn't start the car, but they made it go faster and with more force.
The Analogy: Think of the explosion as a rocket. The neutrinos are the main engine. The magnetic fields are like a secondary booster pack. If the rocket is already flying, the booster pack makes it go much faster and carry more energy.

3. The "Deformed Balloon" (The Neutron Star)

When the star explodes, the core collapses into a tiny, super-dense ball called a Neutron Star.

The Finding: In the models with strong magnetic fields, this new Neutron Star wasn't a perfect sphere. The magnetic pressure squeezed it, making it bulge at the poles, like a balloon being squeezed by a strong hand.
The Consequence: Because the "balloon" was squished, the ghost particles (neutrinos) escaping from it had to travel through different temperatures and densities. This changed the "color" and intensity of the neutrino light we would see from Earth.
The Analogy: If you shout through a round megaphone, the sound goes straight. If you squeeze the megaphone into a weird shape, the sound gets distorted. The strong magnetic fields distorted the Neutron Star, distorting the neutrino signal.

4. The Spin and The Kick

When the star explodes, the leftover Neutron Star often spins incredibly fast and gets kicked away from the center of the explosion (like a spinning top that wobbles and shoots off).

The Finding:
- Spin: In the models with strong magnetic fields, the magnetic forces acted like a giant hand grabbing the spinning star and pulling it faster. The resulting Neutron Stars spun incredibly fast (some in just 0.02 seconds per rotation!).
- Kick: The "kick" velocity depended on how lopsided the explosion was. Interestingly, the model with the strongest magnetic field actually had a more symmetrical (balanced) explosion, resulting in a smaller kick. The model with a weaker field had a more chaotic, one-sided explosion, which gave the Neutron Star a harder kick.

The Big Takeaway

The most important lesson from this paper is about trust.

In the past, scientists might have said, "Our simulation is good enough." But this paper shows that if you don't have enough "pixels" (resolution) in your simulation, you miss the tiny details that turn chaos into powerful magnetic fields.

Low Resolution: Misses the magnetic boost. The explosion looks weaker.
High Resolution: Sees the magnetic boost. The explosion is stronger, the Neutron Star spins faster, and the neutrino signal is different.

In short: To understand how stars die and how the universe gets its heavy elements, we need to look at the stars with the sharpest possible "eyes" (computer simulations), because the invisible magnetic forces are quietly running the show in the background.

1. Problem Statement

Core-collapse supernovae (CCSNe) are complex events where the explosion mechanism is not yet fully understood. While the neutrino-driven mechanism is the leading theory for non-rotating stars, the role of magnetic fields in these scenarios remains debated. Previous studies have focused heavily on magnetorotational explosions (rapidly rotating stars) or purely hydrodynamic models.

The Gap: There is a lack of systematic resolution studies for 3D magnetohydrodynamic (MHD) simulations of neutrino-driven supernovae.
The Question: How do numerical grid resolution and initial magnetic field strength ( $10^{10}$ G vs. $10^{12}$ G) affect shock revival, explosion dynamics, neutrino emission, and the properties of the resulting compact remnant (proto-neutron star, PNS)?
Specific Challenge: In turbulent environments like the CCSN gain region, small-scale dynamos can amplify magnetic fields. The growth rate of these dynamos is resolution-dependent, potentially altering the explosion outcome if the grid is too coarse.

2. Methodology

The authors performed a suite of 3D MHD simulations using the CoCoNuT-FMT code.

Progenitor Model: A non-rotating, magnetized $13 M_\odot$ star (model 13SH21 from Whitehead et al. 2026), evolved using MESA with metallicity $Z=0.001$ .
Initial Conditions:
- Magnetic Fields: Two initial central dipole field strengths were tested:
  - Strong field: $B_{pol,tor} \approx 10^{12}$ G.
  - Weak field: $B_{pol,tor} \approx 10^{10}$ G.
- Geometry: A simple dipole geometry defined by a vector potential (to avoid complications of complex pre-supernova geometries).
Resolution Study: Six models were run with varying grid resolutions ( $N_r \times N_\theta \times N_\phi$ $N_{r} \times N_{θ} \times N_{ϕ}$ ):
- vLRes: $275 \times 64 \times 128$ (Very Low Res)
- LRes: $550 \times 64 \times 128$ (Low Res, radial refinement)
- MRes: $550 \times 128 \times 256$ (Medium Res, standard)
- HRes: $550 \times 256 \times 512$ (High Res)
- Note: MRes and vLRes were run with both strong ( $10^{12}$ G) and weak ( $10^{10}$ G) initial fields.
Physics:
- Hydro/MHD: Newtonian MHD with HLLC Riemann solver, piecewise parabolic reconstruction, and hyperbolic divergence cleaning.
- Neutrinos: Fast Multi-group Transport (FMT) scheme solving the zeroth moment equation for three species ( $\nu_e, \bar{\nu}_e, \nu_x$ ).
- Equation of State: SFHo at high densities; ideal gas + NSE solver at lower densities.

3. Key Contributions

This paper provides the first systematic resolution study of 3D MHD neutrino-driven supernovae. It isolates the effects of numerical resolution on magnetic field amplification via the small-scale dynamo (SSD) and quantifies how these changes influence the explosion mechanism and remnant properties.

4. Key Results

A. Explosion Dynamics and Morphology

Shock Revival: The time of shock revival is largely independent of resolution and initial magnetic field strength for this specific progenitor (occurring around 0.15–0.20 s post-bounce).
Explosion Energy:
- Higher resolution leads to higher diagnostic explosion energies.
- Stronger initial magnetic fields also lead to higher explosion energies.
- The effect of resolution on energy is comparable to the effect of changing the initial magnetic field strength by two orders of magnitude.
Mechanism: Unlike pure hydrodynamic studies where higher resolution increases turbulent pressure, in these MHD models, the increased turbulent kinetic energy is efficiently converted into magnetic energy via the small-scale dynamo. The magnetic flux aids the shock, but the magnetic field does not directly drive the initial revival.

B. Neutrino Emission

Inverse Correlation: Counter-intuitively, higher resolution models exhibit lower neutrino luminosities and mean energies (by $\approx 5-8\%$ ) compared to lower resolution models.
Cause: Higher resolution allows for more efficient magnetic field amplification. Stronger magnetic fields deform the PNS (making it oblate/extended at the poles). This deformation causes neutrinos to be emitted from larger radii (cooler regions), reducing the overall luminosity and mean energy.
Convergence: Pure hydrodynamic simulations often show convergence in neutrino properties at lower resolutions; this study shows that in MHD, properties are not converged even at the highest resolution tested ( $0.7^\circ$ angular resolution).

C. Gain Region Dynamics

Heating Efficiency: While neutrino heating rates ( $\dot{Q}_\nu$ ) are slightly lower in high-resolution models (due to lower neutrino luminosity), the convective efficiency ( $\eta_{conv}$ ) is significantly higher (up to 50% higher).
Energy Conversion: High-resolution grids resolve turbulence better, allowing the small-scale dynamo to convert turbulent kinetic energy into magnetic energy more effectively.
Magnetic Energy Ratios: The ratio of magnetic to kinetic energy ( $E_B/E_{kin}$ ) increases with resolution, reaching up to 35% in the HRes model (compared to ~12% in vLRes).

D. Proto-Neutron Star (PNS) Properties

Convection Zones: Strong magnetic fields and high resolution lead to more extended and vigorous PNS convection zones.
Stress Tensors: In strongly magnetized, high-resolution models, magnetic stresses become comparable to kinetic stresses in the inner regions, altering the energy redistribution.
Spin and Kick:
- Spin: Magnetic torques dominate the spin-up of the PNS. Strongly magnetized models result in rapid spin periods ($0.024$ s to $0.178$ s). Even the weakly magnetized MRes10 model spins up significantly ($8.12$ s) due to field amplification, whereas the low-resolution vLRes10 does not ($280$ s).
- Kick: Kick velocities are stochastic and driven by global asymmetry rather than total explosion energy. MRes10 (weak field, medium res) had the highest kick ($71$ km/s) due to asymmetric flow, while the symmetric MRes had a lower kick despite a stronger explosion.
Magnetar Formation: All strongly magnetized models reached magnetar-level field strengths ( $>10^{14}$ G) at the PNS surface by the end of the simulation.

5. Significance and Conclusions

Resolution Sensitivity: MHD supernova simulations are highly sensitive to spatial resolution. Under-resolving the grid leads to excessive numerical diffusion, which suppresses the small-scale dynamo, preventing the development of strong magnetic fields and altering the explosion energetics.
Non-Convergence: Unlike pure hydrodynamic CCSN simulations which may converge at $\sim 1^\circ$ resolution, these MHD simulations have not converged even at $0.7^\circ$ . Bulk properties (explosion energy, PNS spin, magnetic energy) continue to change with resolution.
Physical Implications:
- Strong magnetic fields in non-rotating stars act as a "subsidiary" driver, compensating for reduced neutrino heating by enhancing magnetic flux and energy redistribution.
- The initial magnetic field geometry and strength significantly impact the birth properties of neutron stars (spin, kick, and potential magnetar formation).
Future Work: The authors emphasize the need for long-duration simulations and code comparisons to fully understand the saturation of these properties and the late-time evolution of the neutrino-driven wind and magnetic fields.

In summary, this study demonstrates that numerical resolution is a critical physical parameter in MHD supernova modeling, not just a numerical artifact, as it dictates the efficiency of magnetic field amplification, which in turn fundamentally alters the explosion dynamics and the nature of the compact remnant.