🔭 astrophysics

Constraining hydrodynamic model of nearby type IIP SN 2023ixf

This paper presents a refined hydrodynamic model of the nearby Type IIP supernova SN 2023ixf that utilizes a broader set of observables, particularly early ejecta velocities, to constrain its explosion energy, ejected mass, and circumstellar environment while simulating its evolution from explosion to hard X-ray emergence.

Original authors: V. P. Utrobin (NRC "Kurchatov Institute", Moscow, Russia, Institute of Astronomy, Moscow, Russia), N. N. Chugai (Institute of Astronomy, Moscow, Russia)

Published 2026-06-16

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: V. P. Utrobin (NRC "Kurchatov Institute", Moscow, Russia, Institute of Astronomy, Moscow, Russia), N. N. Chugai (Institute of Astronomy, Moscow, Russia)

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 supernova as a massive, cosmic firework exploding in a neighboring galaxy. In May 2023, a specific firework called SN 2023ixf went off in the galaxy M101, just 6.7 million light-years away. Because it was so close and so bright, astronomers had a treasure trove of data to study it.

However, there was a problem. When scientists tried to build computer models to explain how this explosion happened, they couldn't agree on the basics. Some models said the explosion was weak; others said it was incredibly powerful. It was like trying to guess the size of a car by looking at its shadow, but everyone was using a different light source, leading to wildly different guesses about the car's actual size.

This paper by Utrobin and Chugai is like bringing in a master detective with a better flashlight. They didn't just look at the usual clues; they focused on a specific, fleeting detail that previous models had ignored.

The Missing Clue: The "Speed Limit"

The authors realized that to solve the mystery, they needed to know exactly how fast the outer layers of the star were flying away right at the beginning of the explosion.

Think of the explosion like a cannonball hitting a thick fog. Usually, the fog (the gas surrounding the star) hides the cannonball, making it impossible to see how fast it's going. But in this case, the "fog" was slightly transparent in a specific way. The team looked at very early, high-quality snapshots of the explosion and found faint "shadows" (absorption lines) in the light. These shadows acted like a speed trap, revealing that the outer edge of the explosion was moving at a blistering 12,000 kilometers per second (about 27 million miles per hour).

This single number was the key. Once they locked this speed in, the "degeneracy" (the confusion where many different answers seemed possible) vanished. Suddenly, the computer model could only fit one specific set of numbers.

The New Portrait of the Explosion

With this new speed limit as a guide, the authors rebuilt the model and found a much clearer picture of what happened:

The Star's Size: The star that exploded was a Red Supergiant, a cosmic giant with a radius about 1,540 times larger than our Sun. If you put this star in our solar system, its surface would swallow the orbit of Jupiter.
The Explosion Power: The energy released was massive, about 2.8 times 10 to the 51st ergs. To put that in perspective, this is a staggering amount of energy, perhaps even more than standard "neutrino-driven" explosions usually produce.
The Debris: The explosion threw out about 13.2 times the mass of our Sun in debris.
The "Fuel": Inside the debris, there was about 0.07 solar masses of a radioactive element called Nickel-56. This is the "battery" that keeps the supernova glowing for months after the initial flash.

The Cosmic "Fog" and the Hard X-Ray Mystery

The paper also solved a puzzle about the environment surrounding the star.

Imagine the star was coughing up a dense cloud of gas (a "circumstellar shell") just before it died. When the explosion hit this cloud, it created a shockwave. Usually, when a shockwave hits gas, it gets hot and emits X-rays. But for SN 2023ixf, the X-rays appeared earlier and behaved differently than expected.

The authors realized the "fog" wasn't a smooth, uniform cloud. Instead, it was clumpy, like a cloud made of dense cotton balls floating in a thin mist.

The Cotton Balls: These dense clumps absorbed the X-rays, explaining why the surrounding gas looked thick when measured by X-ray telescopes.
The Mist: The space between the clumps was thin. The shockwave raced through this thin mist, heating it up so fast that it emitted the hard X-rays we saw.

This "clumpy" structure allowed the model to perfectly match the observations: the gas was dense enough to block light in some ways, but thin enough in other spots to let the X-rays escape early.

The "Zig-Zag" in the Light

The model also explained a weird blip in the supernova's brightness. About 5 days after the explosion, the light dipped slightly before rising again. The authors explain this as a "traffic jam" of radiation. The outer shell of the explosion got so cold that it became very good at trapping light (like a thick blanket). The light built up pressure behind this blanket, then suddenly burst through, causing the dip and the subsequent rise.

The Bottom Line

By focusing on the early speed of the debris and the specific structure of the gas surrounding the star, the authors finally pinned down the true nature of SN 2023ixf. They moved from a range of confusing possibilities to a single, detailed story: a massive, 1500-sun-radius star exploded with immense energy, throwing its debris into a clumpy, pre-existing cloud of gas, creating a complex dance of light and X-rays that we can now finally understand.

This is the first time such a detailed simulation, from the initial explosion all the way to the emergence of hard X-rays, has been successfully modeled for this type of supernova.

Technical Summary: Constraining Hydrodynamic Model of Nearby Type IIP SN 2023ixf

Problem Statement
Despite the proximity of the Type IIP supernova SN 2023ixf and the availability of extensive observational data, previously released hydrodynamic models have failed to converge on a reliable set of explosion parameters. Existing studies (e.g., Bersten et al., 2024; Moriya and Singh, 2024) yield a broad range of inferred explosion energies ( $0.5–3 \times 10^{51}$ erg) and ejected masses, undermining the reliability of current parameter derivations. This "parameter degeneracy" arises because standard observables, such as the light curve plateau luminosity and photospheric velocity, are insufficient to uniquely constrain the model. Furthermore, previous models often invoked massive circumstellar (CS) shells ( $\sim 0.4–0.5 M_\odot$ ) to explain early luminosity peaks, a feature not necessarily supported by all data.

Methodology
The authors revisit the hydrodynamic modeling of SN 2023ixf using the Lagrangian radiation hydrodynamic code CRAB. To resolve the parameter degeneracy, the study incorporates a broader and more specific set of observational constraints than previous works:

Early Maximum Ejecta Velocity: Utilizing shallow Balmer line absorptions detected in early spectra (days 8–14) to determine the maximum expansion velocity ( $v_{sn}$ ). The authors model the formation of these absorptions in the layer between the opaque ejecta core and the opaque circumstellar shell (CDS), accounting for Thomson scattering and CDS optical depth.
$^{56}$ Ni Mixing Extent: Incorporating the velocity extent of $^{56}$ Ni ejecta derived from the blue width of the [Co II] 10.521 $\mu$ m emission line in the JWST spectrum (day 252).
Circumstellar Structure: Constraining the CS shell density and extent using the early H $\alpha$ profile (Thomson scattering wings) and the column density ( $N_H$ ) inferred from X-ray data.
Radiation Hydrodynamics: Simulating the interaction between the forward shock and the CS medium to model the transition from a radiation-dominated shock to a matter-dominated shock, specifically addressing the early emergence of hard X-rays.

The pre-supernova (pre-SN) model is based on a hydrostatic Red Supergiant (RSG) but is modified to include steeper density gradients in the hydrogen envelope and a low-mass CS shell, consistent with early ionization flash spectra.

Key Results
The optimized hydrodynamic model yields the following specific parameters for SN 2023ixf:

Explosion Energy: $2.8 \times 10^{51}$ erg ( $\pm 0.4 \times 10^{51}$ erg).
Ejected Mass: $13.2 M_\odot$ ( $\pm 1.3 M_\odot$ ).
Pre-SN Radius: $1540 R_\odot$ ( $\pm 150 R_\odot$ ).
$^{56}$ Ni Mass: $0.07 M_\odot$ ( $\pm 0.004 M_\odot$ ).
$^{56}$ Ni Mixing Extent: Up to $4000$ km s $^{-1}$ .
Circumstellar Matter: Composed of a dense internal CS shell (mass $\approx 0.01 M_\odot$ , radius $\approx 5 \times 10^{14}$ cm) and a rarefied external wind. The internal shell has a density parameter $w_{int} = 3.7 \times 10^{16}$ g cm $^{-1}$ , while the external wind has $w_{ext} = 1.7 \times 10^{15}$ g cm $^{-1}$ .

The model successfully reproduces the bolometric light curve, including a brief luminosity minimum around day 5 caused by the interplay of CDS optical depth and Rosseland opacity, as well as the evolution of photospheric velocities. The derived pre-SN mass of $14.8 M_\odot$ (ejecta + neutron star) aligns with progenitor estimates of $16–17 M_\odot$ .

Regarding the early hard X-rays, the radiation hydrodynamic simulation indicates that a uniform CS shell density derived from H $\alpha$ and X-ray column density overpredicts the X-ray luminosity on day 5. A better fit is achieved by assuming a clumpy structure for the CS shell, where a rarefied intercloud medium ( $\eta \approx 0.1$ ) allows for hard X-ray emission, while dense clumps provide the necessary occultation to match X-ray absorption data.

Significance and Claims
The authors claim that this study provides the first hydrodynamic simulation of SN 2023ixf that spans from the explosion to the emergence of hard X-ray radiation, explicitly modeling the transition from a radiation-dominated forward shock to an adiabatic, matter-dominated regime.

The primary significance of the work lies in:

Resolving Parameter Degeneracy: By prioritizing the early maximum ejecta velocity (derived from shallow Balmer absorptions) and the $^{56}$ Ni velocity extent, the authors significantly narrow the uncertainty in explosion energy and ejected mass compared to previous models.
Revising CS Shell Models: The study demonstrates that a massive CS shell is unnecessary; instead, a low-mass ( $\sim 0.01 M_\odot$ ) shell is sufficient to explain early spectral features and X-ray data, provided the shell has a clumpy structure.
Radiation Hydrodynamics of Shock Interaction: The paper presents a detailed description of the shock transition mechanism in a dense CS environment, a phenomenon previously un-described for supernovae within the framework of radiation hydrodynamics.

The authors note that the derived explosion energy ( $2.8 \times 10^{51}$ erg) is high and may exceed the energetic limits of standard neutrino-driven explosion mechanisms, suggesting potential complexities in the explosion physics or progenitor structure.