🔭 astrophysics

Probing the 3D Structures of Supernovae through IR Signatures of CO and SiO

This paper introduces MOFAT, a new data-driven tool utilizing multidimensional radiative transfer simulations to analyze near- and mid-infrared CO and SiO signatures, which was applied to SN2024ggi to reveal ongoing SiO formation triggered by CO, significant mass evaporation, and clumped structures while finding no evidence for dust formation.

Original authors: T. Mera, P. Hoeflich, C. R. Burns, C. Ashall, K. Medler, E. Fereidouni, W. B. Hoogendam, M. Shahbandeh, S. Shiber, C. M. Pfeffer, E. Baron, J. Lu, N. Morrell, E. Y. Hsiao, M. M. Phillips

Published 2026-04-21

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: T. Mera, P. Hoeflich, C. R. Burns, C. Ashall, K. Medler, E. Fereidouni, W. B. Hoogendam, M. Shahbandeh, S. Shiber, C. M. Pfeffer, E. Baron, J. Lu, N. Morrell, E. Y. Hsiao, M. M. Phillips

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

The Big Picture: Peeking Inside a Cosmic Explosion

Imagine a massive star dying. It doesn't just fade away; it explodes in a spectacular supernova. For a long time, astronomers have been trying to understand what happens inside that explosion after the initial flash. Specifically, they want to know how dust forms in the universe.

Think of a supernova like a giant, chaotic kitchen. When the star explodes, it throws out all its ingredients (elements like carbon, oxygen, and silicon). As the debris cools down, these ingredients start to stick together to form molecules (like Carbon Monoxide and Silicon Monoxide), which eventually clump together to make dust.

But here's the problem: The "kitchen" is messy. The ingredients aren't mixed evenly; they are swirling in huge, turbulent clouds with pockets of high density and low density. Traditional computer models tried to simulate this by pretending the explosion was a perfect, smooth ball. But when astronomers looked at real data from the James Webb Space Telescope (JWST), those smooth ball models didn't match what they saw. The light coming from the explosion looked different than the models predicted.

The New Tool: MOFAT

To fix this, the authors created a new tool called MOFAT (MOlecular Fitting Analysis Tool).

The Analogy:
Imagine you are trying to figure out what a hidden object looks like by listening to the sound it makes when you tap it.

The Old Way: You assume the object is a perfect sphere. You tap it, hear the sound, and try to guess the material. But the sound doesn't match your guess because the object is actually lumpy and irregular.
The MOFAT Way: Instead of guessing the shape first, MOFAT listens to the sound (the light from the telescope) and works backward. It asks, "What kind of lumpy, irregular shape would make this specific sound?"

MOFAT is a "data-driven" detective. It takes the actual light we see from the supernova and runs complex simulations to find the specific 3D structure (the "lumps" or "clumps") that would create that exact pattern of light.

How It Works: The "Picket Fence" Effect

The paper focuses on two specific molecules: CO (Carbon Monoxide) and SiO (Silicon Monoxide). These molecules glow in infrared light (heat), which JWST can see very well.

The authors discovered that to explain the light, the gas in the supernova can't be a smooth fog. It has to be made of clumps, like a picket fence.

The Overtone Bands (The "Thin" Light): Some of the light passes through easily, like looking through the gaps in a picket fence. This tells us how much total gas is there.
The Fundamental Bands (The "Thick" Light): Other light gets blocked by the slats of the fence. This tells us about the shape and density of the clumps.

By analyzing how much light gets blocked versus how much gets through, MOFAT can reconstruct the 3D shape of the explosion's interior.

The Case Study: SN 2024ggi

The team tested MOFAT on a real supernova called SN 2024ggi, looking at it 285 days and 385 days after the explosion. Here is what they found:

The Domino Effect: It turns out that CO is the "big brother" that has to do the heavy lifting first. The CO molecules form, cool down the surrounding gas, and then allow the SiO molecules to form. Without the CO cooling things down, the SiO couldn't survive.
The Inner Edge is Moving: The region where SiO exists is shrinking. The inner edge of the SiO cloud is moving inward (slowing down), which suggests that SiO is still being formed in the deep layers, but the outer layers are being destroyed by harsh radiation.
No Dust Yet: Even though we see these molecules, they aren't cooling the gas enough to start building dust grains yet. The gas is still too hot.
Clumps are Real: The data strongly suggests that the gas is clumpy. However, because the SiO light is a bit "fuzzy" (not thick enough to show the clumps clearly), the team couldn't say exactly what shape those clumps are (like whether they are flat pancakes or long cigars).

Why This Matters

This paper is a breakthrough because it moves us from "guessing" what the inside of a supernova looks like to measuring it.

Before: We had to assume the explosion was a smooth ball, which led to wrong answers about how much dust is made.
Now: We have a tool (MOFAT) that can map the 3D structure of the explosion, revealing the "instabilities" (the turbulence and mixing) that happened during the blast.

The Bottom Line:
Just like a doctor uses an X-ray to see inside a body without cutting it open, MOFAT uses infrared light to see the 3D structure of a dying star. This helps us understand how the universe recycles old stars into new dust, which is the building block for new stars, planets, and eventually, us.

The authors conclude that while MOFAT is powerful, it needs more time to work its magic. To get the clearest picture, we need to watch these explosions from the very beginning of their cooling phase all the way to when they turn into dust.

1. Problem Statement

Core-collapse supernovae (CC SNe) are primary sources of interstellar dust and heavy elements. However, there is a longstanding discrepancy between theoretical predictions and observations regarding dust masses in the early universe. A critical step in dust formation is the synthesis of molecular coolants, specifically Carbon Monoxide (CO) and Silicon Monoxide (SiO), which radiate in the near- and mid-infrared (NIR/MIR).

Key Challenges:

Limitations of 1D Models: Traditional time-dependent, spherical forward-modeling fails to reproduce the observed flux ratios between the fundamental and overtone vibrational bands of CO and SiO (often off by factors of 2–5). These models cannot account for the complex, clumpy 3D structures resulting from hydrodynamic instabilities (e.g., Rayleigh-Taylor, Richtmyer-Meshkov).
Optical Depth Effects: The fundamental bands are optically thick, while overtone bands are optically thin. Spherical models struggle to simultaneously fit both regimes because the temperature profile is dictated by optical depth, which is directly tied to the degree of clumping.
Computational Complexity: Fully 3D time-dependent radiation-hydrodynamic simulations are computationally prohibitive for systematic parameter exploration, while 1D/2D models lack the necessary morphological realism.

2. Methodology: MOFAT

The authors introduce MOFAT (MOlecular Fitting Analysis Tool), a novel inverse, data-driven approach designed to probe the 3D structure of molecular regions in SNe.

Core Concept: Instead of testing random explosion instabilities against observations, MOFAT uses high-precision JWST spectra to drive stationary (time-independent) radiative transfer simulations. It iteratively converges on a set of physical parameters that best reproduce the observed spectra.
Radiative Transfer Framework:
- Assumes a spherically symmetric large-scale geometry with a homologously expanding envelope ( $\rho \sim r^{-n}$ ).
- Incorporates multidimensional clump-like structures using a probabilistic "picket-fence" paradigm. Clumps are modeled as spheroids with specific shape ( $\epsilon$ ), relative size ( $r_R$ ), and density enhancement ( $f_c$ ).
- Solves the radiation transport equations in the observer's frame using the impact parameter ray system and Sobolev approximation for rest-frame opacities.
Fitting Procedure:
- Two-Step Sequential Fitting: First, CO features are fitted to determine the global temperature structure ( $T(v)$ ) and density slope. This temperature profile is then used as a fixed input for fitting SiO features.
- Optimization: Uses a nested-loop iteration system minimizing a $\chi^2$ statistic (and a shape-based $\chi^2_s$ ) to converge on seven key parameters: inner/outer velocity boundaries ( $v_1, \Delta v$ ), density slope ( $n$ ), inner temperature ( $T_1$ ), and three clump parameters ( $f_c, r_R, \epsilon$ ).
- Constraints: The method relies on the absolute flux ratio between fundamental and overtone bands to uniquely constrain the solution space.

3. Key Contributions

Development of MOFAT: A robust tool that bridges the gap between complex 3D hydrodynamics and observational data, allowing for the reconstruction of small-scale structures without full 3D hydro simulations.
Resolution of Flux Ratio Discrepancy: Demonstrates that including small-scale 3D clumping is essential to reproduce the observed flux ratios between CO/SiO fundamental and overtone bands, overcoming the failures of traditional 1-zone models.
Diagnostic Framework: Establishes a method to distinguish between optically thin (mass/temperature diagnostics) and optically thick (clump morphology diagnostics) regimes.

4. Results: Application to SN 2024ggi

The tool was applied to JWST NIRSpec and MIRI-LRS observations of the core-collapse supernova SN 2024ggi at +285 and +385 days post-explosion.

CO Analysis:

Clumping Confirmation: Clumped models (prolate and oblate) were required to fit the data. The best-fit models indicate significant small-scale structure.
Parameter Convergence: The models converged on a density slope $n \approx 6$ , an inner temperature $T_1 \approx 1950$ K, and specific clump geometries ( $f_c \approx 1.8$ , $r_R \approx 0.14$ , $\epsilon \approx 14$ ).
Evolution: CO mass decreased from $\sim 8.7 \times 10^{-3} M_\odot$ to $\sim 1.3 \times 10^{-3} M_\odot$ , indicating active destruction due to exposure to hard radiation as the ejecta expanded.

SiO Analysis:

Formation Trigger: SiO formation is triggered by CO cooling. The inner edge of the SiO-emitting region receded from $v_1 \approx 1500$ km/s to $1000$ km/s, indicating continued formation in inner layers.
Mass and Cooling: SiO mass decreased by an order of magnitude (from $\sim 2-6 \times 10^{-3} M_\odot$ to $\sim 10^{-4} M_\odot$ ). Crucially, SiO does not significantly cool the ejecta; CO dominates the cooling process.
Dust Formation: No evidence for dust formation was found at these epochs; temperatures remained too high ( $>1900$ K) for dust nucleation.
Clump Ambiguity: While SiO features suggest clumping, the fundamental band is semi-optically thick, leading to shallow convergence on clump parameters. The results are inconclusive regarding the specific 3D morphology of SiO clumps.

5. Significance and Future Outlook

Tomography of SNe: MOFAT enables "spectral tomography," allowing astronomers to reconstruct the chemical stratification (C/O/Si layers) and instability scales of the progenitor star by tracking the evolution of molecular regions over time.
Dust Budget Implications: By clarifying the cooling mechanisms (CO vs. SiO) and the destruction rates of molecules, the study refines our understanding of how much dust can ultimately form in SNe, addressing the "missing dust" problem in the early universe.
Observational Strategy: The paper emphasizes the need for low-latency, extended time-series observations. Early observations are required to capture the onset of molecule formation and mixing, while late-time observations are necessary to track the evolution of clumps and the transition to dust formation.
Limitations: The current implementation assumes a uniform distribution of clumps with identical shapes/sizes. Future work requires high-spectral coverage across the entire formation period and potentially full 3D time-dependent simulations to distinguish between different instability types (e.g., Rayleigh-Taylor vs. shear flows) and viewing angle effects.

In summary, this paper presents a paradigm shift in supernova modeling, moving from forward-modeling limitations to an inverse, data-driven approach that successfully utilizes JWST data to probe the 3D physical structure of supernova ejecta.