New Simulation Reveals Supernova Explosions Are Three Times Weaker Than Thought

When massive stars collapse and explode as core-collapse supernovae, they release monumental torrents of thermal and kinetic energy that shape the evolution of entire galaxies. For decades, standard galaxy simulations have relied on a simplifying assumption that every dying massive star packs a punch of exactly one billion times ten to the power of 51 ergs. A team of astrophysicists has now challenged this standard yardstick by simulating a low-mass dwarf galaxy using a more realistic, variable energy model. Their findings suggest that the average explosion energy of these stellar blasts is actually only about one-third of the traditional value, a revelation that alters how scientists understand the cosmic matter cycle.

The study, available on the arXiv preprint server, was conducted by Natalia Lahén and her colleagues at the Max Planck Institute for Astrophysics, the Zentrum für Astronomie at Heidelberg University, and the Max Planck Institute for Astronomy. To capture the complex physics of dying stars, the researchers implemented a novel prescription that ties a star’s explosion energy and explodability directly to the mass of its carbon-oxygen core and its chemical composition, or metallicity. Unlike older models where every star between eight and forty solar masses detonates with identical force, this advanced approach accounts for failed supernovae that collapse directly into black holes, as well as partial fallback events where stellar material tumbles back onto the newly formed compact remnant.

Why Fading Stellar Detonations Change Galaxy Growth

When the team integrated this new energy framework into their high-resolution hydrodynamical simulations of a dwarf galaxy modeled after the Wolf-Lundmark-Melotte system, the consequences rippled through the entire simulated 1 billion-year evolution. Under the new model, the average explosion energy dropped to approximately 0.33 times the traditional unit. Because these stellar blasts were fundamentally weaker on average than previously assumed, they drove less aggressive galactic outflows of gas. Consequently, the mass-loading factor, which measures the ratio of gas pushed out of the galaxy compared to the rate of new star formation, fell below unity.

This quieter, less violent feedback loop meant that gas lingered longer within the dwarf galaxy rather than being forcefully ejected into intergalactic space. With more fuel readily available, the simulated galaxy maintained a steadier and higher rate of star formation over time. The increased star formation naturally led to enhanced chemical enrichment, allowing heavy elements forged within stars to recycle more efficiently into newly forming stellar generations. The authors noted that this steady production of metals resulted in an average increase in stellar metallicity by a factor of over two compared to runs using the traditional, highly energetic supernova feedback prescription.

Surprising Consequences for Hidden Black Holes and Star Clusters

The revised explosion model also reshaped the population of dark remnants left behind by dead stars. By incorporating specific core-mass thresholds for stellar collapse, the simulation produced nearly three times as many black holes as the canonical model. Interestingly, the researchers observed a pronounced deficit in the mass distribution of these black holes between roughly ten and fifteen solar masses, a feature that aligns well with recent gravitational wave observations of binary black hole mergers reported by LIGO and Virgo collaborations.

Furthermore, the team analyzed young star clusters containing stars younger than 10 million years old to see how the altered feedback affected cluster formation efficiency. They discovered that the reduced supernova energies fostered a higher total number of star clusters, roughly doubling the count seen in standard 1-billion-erg simulations. These clusters formed in denser ambient gas environments because weaker supernova shocks were less effective at carving out massive, galaxy-wide bubbles early in the star formation cycle.

While the study successfully bridges the gap between detailed three-dimensional supernova physics and galactic-scale evolution, the authors emphasize that certain limitations remain. For instance, the current simulations exclude the complex binary star interactions that frequently occur in nature and alter stellar envelopes before collapse. Future research will aim to incorporate stellar binary evolution models to further refine how these cosmic powerhouses govern the lifecycle of galaxies across the universe.

The study was posted to the arXiv preprint server on November 10, 2026.