🔭 astrophysics

Supernova interactions with aspherical circumstellar material I: calculations of light curves, AB magnitudes, spectra, and polarisation

This paper presents upgraded 2D radiation-hydrodynamic simulations of supernova interactions with aspherical circumstellar material, specifically discs and bipolar lobes, to calculate and compare light curves, spectra, and polarization profiles over two years, revealing distinct observational signatures that help identify the underlying CSM geometry in real supernova events.

Original authors: Petr Kurfürst, Georgi Bless, Jakub Fišák, Filip Holoubek, Jiří Krtička, Brankica Kubátová, Jiří Kubát, Michal Zajaček

Published 2026-04-02

📖 5 min read🧠 Deep dive

CC0 1.0

Original authors: Petr Kurfürst, Georgi Bless, Jakub Fišák, Filip Holoubek, Jiří Krtička, Brankica Kubátová, Jiří Kubát, Michal Zajaček

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.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 not just as a star exploding, but as a massive firework going off inside a complex, invisible maze. This paper is about simulating what happens when that firework explodes inside two very different types of mazes: a flat, spinning disc (like a pizza dough being tossed) and bipolar lobes (like a dumbbell or the hourglass shape of a nebula).

Here is a breakdown of the research using simple analogies:

1. The Setup: The Exploding Star and the "Maze"

When a massive star dies, it explodes. Usually, scientists assume the debris (the "ejecta") flies out in a perfect sphere, like a balloon inflating evenly in all directions.

But in reality, stars often spin or have companions, which can fling material out before they explode. This creates a "circumstellar material" (CSM) cloud around the star.

The Disc Model: Imagine the star was spinning so fast it spun off a ring of gas, like a figure skater spinning and flinging their arms out.
The Lobes Model: Imagine the star had a violent cough before dying, shooting out two giant puffs of gas in opposite directions, like a dumbbell.

The researchers wanted to see what happens when the explosion hits these specific shapes.

2. The Simulation: A High-Tech Video Game

The team used powerful supercomputers to run a "video game" of physics.

The Engine (CASTRO): This part simulates the explosion itself. It tracks how the gas moves, how hot it gets, and how the shockwave crashes into the surrounding gas. They added a new feature: radiation. In previous versions, they ignored the heat and light pressure; now, they included it, which acts like a "brake" on the explosion, slowing it down more realistically.
The Camera (SEDONA & SIROCCO): Once the explosion happens in the simulation, they use these tools to act as "virtual cameras." They calculate what an observer on Earth would actually see from different angles. They don't just look at the explosion; they look at the light curves (how bright it gets over time), the colors (spectra), and even how the light is polarized (the direction the light waves vibrate).

3. The Results: What Happens When the Shockwave Hits?

Scenario A: The Disc (The Pizza Dough)

The Interaction: When the explosion hits the disc, it's like a cannonball hitting a flat wall. The gas flies easily out the top and bottom (the poles), but it gets stuck and slowed down when it tries to go through the thick disc (the equator).
The Light: The brightness changes smoothly. It's a steady glow. If you look from the side (through the disc), it might look slightly different than if you look from the top, but the overall story is a smooth, steady fade.
The Sound (Spectra): The light from the explosion creates a "double-horn" shape in the spectrum (like a W or M shape). This is a signature that tells astronomers, "Hey, there's a disc here!"

Scenario B: The Lobes (The Dumbbell)

The Interaction: This is much more chaotic. The explosion hits the inner part of the lobes, then the outer part. It's like running through a hallway with doors that open and close at different times.
The Light: The brightness is "bumpy." It goes up, drops down, and then spikes up again. These bumps happen because the shockwave hits different dense layers of gas at different times.
The Sound (Spectra): The light patterns are more complex and "spiky," reflecting the uneven, lumpy nature of the lobes.

4. The "Polarization" Clue

The researchers also looked at polarization. Think of light as a rope being shaken. If you shake it up and down, it's polarized vertically. If the explosion is perfectly round, the light is a mix of all directions (unpolarized). But if the explosion is squashed or stretched (like hitting a disc or lobes), the light waves get lined up in a specific direction.

The Finding: The amount of polarization acts like a fingerprint. It tells astronomers exactly how "squashed" the explosion is and which way the debris is pointing, even if they can't see the shape directly.

5. Why Does This Matter?

Imagine you are a detective trying to solve a crime, but you only see the smoke and the sound, not the weapon.

The Detective Work: Astronomers see a supernova light curve that goes up and down, or a spectrum that looks weird. They don't know why.
The Solution: This paper provides a "cheat sheet" or a library of simulations. If an astronomer sees a "bumpy" light curve, they can say, "Ah, that matches our 'Lobes' simulation! The star must have ejected gas in a dumbbell shape before it died."
The Goal: By comparing real observations with these simulations, we can figure out how massive stars die, how they lose mass, and what kind of environment they lived in.

Summary

This paper is like building a library of "what-if" scenarios for exploding stars. By using supercomputers to simulate explosions hitting different shapes of gas clouds, the team created a guidebook. This guidebook helps real-world astronomers decode the messy, colorful, and bright signals they see from dying stars, turning a blurry explosion into a clear story about the star's final moments.

1. Problem Statement

The interaction between supernovae (SNe) and dense, aspherical circumstellar matter (CSM) significantly alters observable characteristics such as light curves (LCs), spectra, and polarization. While previous studies (e.g., Kurfürst et al. 2020, hereafter PK20) explored these interactions using adiabatic hydrodynamic models, they lacked a self-consistent treatment of radiation hydrodynamics (RHD) and detailed non-local thermodynamic equilibrium (NLTE) radiative transfer.

The primary objectives of this work are to:

Perform upgraded, self-consistent 2D radiation-hydrodynamic (RHD) simulations of spherical SNe interacting with two distinct aspherical CSM morphologies: circumstellar discs and bipolar lobes.
Calculate robust observables (LCs, AB magnitudes, broadband spectra, spectral line evolution, and polarization) using multidimensional Monte Carlo radiation transfer (MC-RT) codes.
Investigate how varying CSM densities and viewing angles affect the evolution of these observables over a period of up to two years post-shock breakout.
Provide a theoretical framework to help identify the physical configuration (morphology and mass loss rate) of real SN events, particularly Type IIn and super-luminous SNe.

2. Methodology

Numerical Setup and Hydrodynamics

Code: The authors use the CASTRO code, a 2D cylindrical (axisymmetric) RHD solver.
Physics: The simulations solve the conservation equations for mass, momentum, and total energy, coupled with a grey radiation solver (flux-limited diffusion). Radiative cooling and heating are included, unlike previous adiabatic models.
Initial Conditions:
- Progenitor: A red supergiant (RSG) with $M_\star \approx 12.3 M_\odot$ and $R_\star \approx 1000 R_\odot$ . The explosion is initiated with $10^{51}$ erg energy and $0.05 M_\odot$ of radioactive $^{56}$ Ni (though radioactive heating is not included in the RHD evolution phase, only in the comparison LCs).
- CSM Configurations: Two morphologies are modeled with varying density levels:
  1. Circumstellar Discs (DMs): Models DMIN (low density) and DMAX (high density, $100\times$ DMIN).
  2. Bipolar Lobes (LMs): Models LMIN, LMOD, and LMAX, scaled to represent increasing mass loss rates ( $10^{-5}$ to $10^{-1} M_\odot \text{yr}^{-1}$ ).
Grid: A $4200 \times 8400$ grid in radial ( $\varpi$ ) and vertical ( $z$ ) directions, extending to $R_{max} \approx 500 R_\star$ .
Gravity: The Poisson solver is used to include gravitational forces for denser models.

Radiative Transfer and Observables

Codes: Post-processed MC-RT calculations are performed using SEDONA (for LCs and broadband spectra) and SIROCCO (for specific spectral line profiles).
SEDONA:
- Uses 250 time steps from CASTRO as input.
- Calculates monochromatic luminosity $L_\nu(t, \mu)$ and AB magnitudes for 9 filters (UV to IR).
- Frequency range: $10^{14} - 2 \times 10^{16}$ Hz.
- Includes NLTE effects for electron scattering, bound-bound, free-free, and bound-free opacities.
SIROCCO:
- Used specifically for the temporal evolution of the Ly- $\alpha$ line (and H $\alpha$ in analysis).
- Simulates 2.5D axisymmetric outflows with a central compact object (neutron star remnant approximation).
Polarization: Calculated semi-analytically using the Brown & McLean (1977) formalism for Thomson scattering in optically thin regions ( $\tau \le 1$ ).

3. Key Contributions

Transition to RHD: The study moves beyond adiabatic approximations to include full RHD, revealing that radiative cooling leads to significantly slower expansion velocities (approx. 60% of adiabatic models) and different density structures.
Multidimensional Observables: Unlike previous semi-analytical approaches, this work provides direction-dependent (viewing angle) LCs, magnitudes, and spectra derived directly from hydrodynamic snapshots.
Systematic Morphology Comparison: A direct comparison between disc-like and lobe-like CSM interactions across a wide range of densities, establishing distinct signatures for each morphology.
Spectral Line Evolution: Detailed modeling of the Ly- $\alpha$ and H $\alpha$ line profiles, specifically tracking the transition from blue-shifted to red-shifted asymmetry due to resonant scattering.

4. Key Results

Hydrodynamic Evolution

Disc Interaction (DMs): The SN ejecta propagate relatively unimpeded in the polar direction but are decelerated in the equatorial plane. In the high-density DMAX model, a spiral-arm-like clump of accumulated material forms within the disc opening angle.
Lobe Interaction (LMs): The ejecta are strongly confined by the lobes. In the LMAX model, expansion is almost exclusively channeled along the polar axis. The LMIN model shows partial deformation, while LMOD exhibits intermediate behavior.
Velocity: RHD models show a significantly slower shock propagation compared to PK20's adiabatic models (e.g., reaching the domain edge in ~700 days vs. ~400 days).

Light Curves (LCs) and Magnitudes

Disc Models (DMs): Exhibit smooth evolution without significant drops or re-brightening. The DMAX model is brighter than DMIN. Polar and equatorial views diverge only in the early phase (<100 days); later, they converge as UV radiation is reprocessed into longer wavelengths.
Lobe Models (LMs): Show a "bumpy" evolution with distinct peaks and drops.
- First Peak: Caused by the shock hitting the inner lobes/polar regions.
- Drop: Occurs around 100 days as the shock traverses lower-density regions.
- Second Peak: Caused by the shock hitting the outer lobes in the equatorial direction.
Viewing Angle: Shorter wavelengths (UV) show much more pronounced variability and angle dependence than bolometric or visual bands, making UV observations crucial for diagnosing CSM geometry.

Spectral Patterns

Broadband Spectra: Show clear Balmer jumps and continuum features. SEDONA currently displays absorption curves for lines but lacks full emission physics for nebular phases.
Ly- $\alpha$ Evolution:
- Early Phase (<100 days): Broad, slightly blue-shifted emission plateau.
- Late Phase (>100 days): Transitions to a double-peaked (horn-like) profile.
- Asymmetry: A distinct transition from blue-shifted to red-shifted dominance occurs over time. This is attributed to resonant scattering in the expanding neutral hydrogen envelope (photons diffusing to the red wing where opacity is lower).
H $\alpha$ Evolution: Initially a narrow emission peak (Type IIn signature), later broadening into a double-peaked profile. The models predict H $\alpha$ to be weaker in equivalent width than Ly- $\alpha$ .

Polarization

All models exhibit a prolate ejecta distribution (positive shape factor $\gamma > 1/3$ ).
Maximum relative polarization ( $P_R$ ) occurs in the equatorial view, ranging from 0.5% to 2.5% depending on the model density and time.
Polarization values are consistent with observed Type IIn SNe.

5. Significance and Future Work

Diagnostic Tool: The distinct LC shapes (smooth vs. bumpy) and spectral line profiles (double-peaked vs. flat-topped) provide a diagnostic toolkit to distinguish between disc and lobe CSM morphologies in real SN observations.
Observational Guidance: The strong angle dependence in UV bands suggests that future UV missions (e.g., ULTRASAT, QUVIK, UVEX) are critical for constraining CSM geometry.
Limitations: Current models are 2D axisymmetric. 3D simulations would likely show faster instability growth (Rayleigh-Taylor) and more complex filamentary structures. Dust is currently excluded but is planned for future iterations.
Next Steps: The authors plan to refine progenitor and CSM parameters to match specific observed events (e.g., iPTF14hls, SN 1988Z) and extend simulations to parsec scales to include interactions with galactic nuclei and jets.

In summary, this paper establishes a rigorous, multi-messenger framework for understanding SN-CSM interactions, demonstrating that radiation hydrodynamics and viewing angle are critical for interpreting the complex light curves and spectra of interacting supernovae.