---
title: 'SISSI Project: Supernovae in a Stratified ISM'
url: https://www.emergentmind.com/topics/sissi-project
type: topic
---

# SISSI Project: Supernovae in a Stratified ISM

The **SISSI Project**—**“Supernovae In a Stratified, Shearing Interstellar medium”**—is a simulation program devoted to the evolution of supernova remnants (SNRs) in a galactic interstellar medium shaped by **galactic rotation, gravity, turbulence, and vertical stratification**, rather than in an idealized homogeneous background. Its central purpose is to quantify how realistic Galactic environments alter remnant morphology, dynamics, and interpretation, with particular emphasis on geometric diagnostics and on comparison with nearby structures such as the **Local Bubble** [2503.12977][2509.04221].

## 1. Scientific rationale and project scope

SISSI was introduced in response to a mismatch between classical SNR theory and the observational reality of evolved Galactic cavities. Standard analytic and numerical treatments are highly successful for spherical blast waves in uniform media, but they are poorly matched to remnants and superbubbles embedded in a live, anisotropic ISM. The first SISSI paper states this explicitly by targeting SNRs evolving in a **complex interstellar medium structured under the influence of galactic rotation, gravity and turbulence**, and by framing geometry as a primary observable rather than a secondary by-product [2503.12977].

This framing matters because nearby cavities are now observed with genuine three-dimensional information. In the SISSI interpretation, structures such as the Local Bubble are not merely energy-driven shells whose size can be read off from one-zone models; they are objects whose **axis ratios, orientations, and deviations from spherical symmetry** encode the effects of stratification, shear, low-density channels, and ambient substructure. The second SISSI paper extends that program from general remnant morphology to a specific local application, using the Local Bubble as a constraint on recent supernova activity and on the recent evolution of the solar-neighborhood ISM [2509.04221].

A plausible implication is that SISSI is best understood not only as a simulation suite, but also as a methodological program: it treats remnant geometry as a way to infer the physical processes shaping the Galactic ISM. That emphasis distinguishes it from work focused only on global energy or momentum budgets.

## 2. Numerical framework and simulation suite

SISSI uses a deliberately multiscale design. The parent environment is an **isolated Milky-Way-like disk galaxy** from the **AVALON** project, evolved with **RAMSES** in a **48 kpc** cubic box with **outflow boundaries** and a **\(256^3\)** coarse grid. The galaxy-scale ISM reaches a finest resolution of approximately **\(11.7\,\mathrm{pc}\)**, while the zoomed SNR calculations reach approximately **\(0.18\,\mathrm{pc}\)** near the explosion site. The galaxy includes gas self-gravity, a static axisymmetric background potential, radiative cooling and heating, and star formation under the conditions **\(n_{\rm H}>100\,\mathrm{cm^{-3}}\)** and **\(T<150\,\mathrm{K}\)** with **\(\epsilon_{\rm ff}=1\%\)** [2503.12977].

The ISM is first relaxed for less than **\(500\,\mathrm{Myr}\)**, during which star particles deposit **\(2\times 10^{52}\,\mathrm{erg}\)** and **\(200\,M_\odot\)** after **\(8\,\mathrm{Myr}\)**, with cooling disabled for approximately **\(500\,\mathrm{kyr}\)** in affected cells. SISSI then selects **30** future explosion sites at galactocentric radii **\(R\in\{2,\;4.5,\;8\}\,\mathrm{kpc}\)**, close to the midplane and roughly equally spaced in azimuth. Around each site, cells within
$$
r_{\rm zoom},l}=N_{\rm zoom}\,\Delta x_l,\qquad N_{\rm zoom}=15,
$$
are refined up to **\(l_{\max}=18\)**, corresponding to **\(\Delta x_{\min}\approx 0.18\,\mathrm{pc}\)**. The local region is then relaxed for less than **\(50\,\mathrm{kyr}\)** before the supernova is injected [2503.12977].

The suite contains control runs without new supernovae and three explosion prescriptions that define the main SISSI comparison set:

| Model | \((N_{\rm SN},\Delta t_{\rm SN})\) | Interpretation |
|---|---:|---|
| N1 | \((1,\infty)\) | single-SN remnant |
| N10 | \((10,\infty)\) | 10 SNe exploding all at once |
| N1x10 | \((1,1\,\mathrm{Myr})\) | 10 SNe exploding sequentially |

Each SN injects **\(10^{51}\,\mathrm{erg}\)** of thermal energy and **\(5\,M_\odot\)** of ejecta inside a sphere of radius **\(R_{\rm inj}=5\Delta x_{\min}\approx 0.92\,\mathrm{pc}\)**. Cooling remains active during the zoom-in SNR runs. Refinement is also tied to the shell-formation scale through
$$
\Delta x > 0.1\,R_{\rm sf}^{\rm KO15}(n_{\rm H}) = 2.3\left(\frac{n_{\rm H}}{\mathrm{cm^{-3}}}\right)^{-0.42}\,\mathrm{pc},
$$
and the Sedov-Taylor phase is resolved up to
$$
n_{\rm H,max}\sim 430 \left(\frac{\Delta x}{\Delta x_{\min,{\rm zoom}}}\right)^{-2.4}\,\mathrm{cm^{-3}}.
$$
The first paper is also explicit about what is absent: it does **not** include magnetic fields, cosmic rays, thermal conduction, or non-equilibrium chemistry [2503.12977].

## 3. Geometric formalism and analytic baselines

The first SISSI paper is fundamentally a geometry paper, and it defines a formal shape analysis for each remnant. SNR material is identified with passive ejecta scalars, and the geometry is summarized by the volume-weighted shape tensor
$$
S_{ij} = V_{\rm SNR}^{-1}\int_{\rm SNR}\left(\lVert \mathbf{x}\rVert^2\delta_{ij}-x_i x_j\right)\,d^3\mathbf{x}.
$$
If the remnant is approximated as an ellipsoid, the three equivalent ellipsoidal radii follow from the eigenvalues \(S_i\) as
$$
r_i=\sqrt{2.5\left({\rm tr}(S)-2S_i\right)},
$$
with minor, intermediate, and major axes denoted by \(a\), \(b\), and \(c\), and an effective radius
$$
r_{\rm eff}=(abc)^{1/3}.
$$
Orientation is tracked by a pitch angle \(\alpha\) relative to Galactic rotation and by \(|\cos\theta|\), where \(0\) is in-plane and \(1\) is perpendicular to the disk [2503.12977].

SISSI compares these geometric diagnostics to analytic expectations for isolated remnants in uniform media. The paper explicitly quotes the standard shell-formation and snowplow scalings, including
$$
t_{\rm sf}\sim 44\,E_{51}^{0.22}n_0^{-0.55}\,\mathrm{kyr},
$$
$$
R_{\rm sf}=22.6\,E_{51}^{0.29}n_0^{-0.42}\,\mathrm{pc},
$$
and
$$
p_{\rm sf}=2.17\times10^5\,E_{51}^{0.93}n_0^{-0.13}\,M_\odot\,\mathrm{km\,s^{-1}}.
$$
For the sequential-supernova case, it also uses the momentum-driven superbubble scaling
$$
R = 40\,t_6^{1/2}L_{45}^{0.23}n_0^{-0.28}\,\mathrm{pc}.
$$
To isolate the role of differential rotation, SISSI introduces a simple shearing-sphere model, whose characteristic deformation time is
$$
t_{\rm deform}\approx 0.065\,t_{\rm orb}(\bar R).
$$
This model becomes the reference against which the simulated loss of sphericity is measured [2503.12977].

The methodological point is important: SISSI does not use geometry descriptively. It treats geometric quantities as dynamical observables that can be compared quantitatively with both analytic theory and observational reconstructions.

## 4. Principal findings on remnant evolution and morphology

One of the clearest SISSI results is that **young** remnants remain well described by standard theory even in a complex ISM. The first paper reports that shell-formation time, pressure-driven snowplow duration, characteristic radii, and momentum per SN all agree well with isolated-uniform-medium expectations at early stages, especially for moderate and high ambient densities. This means the realistic galactic environment does not erase the validity of classical SNR theory; rather, it limits where that theory remains applicable [2503.12977].

The main departure occurs at later times. SISSI finds that remnants begin to lose spherical symmetry after roughly **a percent of the local orbital timescale**, and often even earlier. The paper defines the deformation time through the condition
$$
a/c<2/3,
$$
and reports that most remnants satisfy it within **\(\lesssim 1\%\)** of an orbit, with some as early as **\(10^{-3}\)** of an orbit. By **\(10\,\mathrm{Myr}\)**, only **4 of the 30** simulated remnants remain spherical. Most have become distinctly **prolate** or **oblate**, with dense-environment cases reaching approximately
$$
a/c\sim 0.2.
$$
In the N1 sample, single-SN remnants tend toward lower \(a/b\), whereas the sequential model N1x10 is periodically “re-sphericalized” by later explosions and remains closer to the simple shear trajectory [2503.12977].

This deformation is too rapid to be explained by direct galactic shear acting on an initially spherical cavity. SISSI therefore concludes that **galactic shear alone cannot explain** the observed timing and degree of asymmetry. Yet the preferred orientations remain consistent with a shear-related origin: oblate remnants tend to have minor axes pointing out of the disk, while prolate remnants tend to have major axes lying in or near the disk plane with pitch angles clustered around values expected from the shearing-sphere model. The paper’s interpretation is that remnants are expanding into ambient density structures that were themselves pre-shaped by galactic shear. That reading is reinforced by the reported correlation
$$
t_{\rm deform}\propto (\delta\rho/\rho)^{-3},
$$
which ties earlier deformation to larger ambient density anisotropy [2503.12977].

The result is conceptually significant. In SISSI, late-time remnant geometry is not dominated by the remnant’s own internal symmetry; it is dominated by the anisotropic structure of the surrounding ISM.

## 5. The Local Bubble program

The second SISSI paper applies the suite to the **Local Bubble** by reconstructing its geometry from the public **3D dust maps** of **Edenhofer et al. (2024)**. It uses the **12 sample dust maps** from that release, follows the shell-finding strategy of **O’Neill et al. (2024)** with modifications, and smooths the density field with **\(\sigma_{\rm smth}=9\,\mathrm{pc}\)** using a mass-conserving kernel. The shell geometry is then summarized with the same shape-tensor formalism used in the simulations [2509.04221].

The reconstructed Local Bubble properties are highly specific. SISSI reports an effective radius
$$
212.3\pm1.0\,\mathrm{pc},
$$
a mass
$$
(6.24\pm0.07)\times10^5\,M_\odot,
$$
and a mean hydrogen number density
$$
0.494\pm0.006\,\mathrm{cm^{-3}}.
$$
Its axis ratios are **minor-to-major \(=0.469\pm0.007\)** and **semi-major-to-major \(=0.562\pm0.011\)**, with a **major-axis pitch angle \(-19.0^\circ\pm1.0^\circ\)** and a **minor-axis pitch angle \(-66.1^\circ\pm2.7^\circ\)**. Relative to earlier reconstructions, the paper states that this Local Bubble is **larger** and embedded in a **much lower mean ambient density** [2509.04221].

SISSI then combines this geometry with a momentum estimator. Assuming homologous expansion with \(R\propto t^{1/2}\),
$$
p_{\rm LB}(t_{\rm age}) \approx \frac{1}{2t_{\rm age}} \int\!\!\int_0^{R_{\rm out}(\Omega)} \mu\,n_{\rm H}(r,\Omega)\,r^3\,\mathrm{d}r\,\mathrm{d}\Omega.
$$
After calibrating this estimator against the SISSI sample, the corrected Local Bubble momentum becomes
$$
p_{\rm LB} \simeq (2.96\pm0.05)\times10^7\left(\frac{\mathrm{Myr}}{t_{\rm age}}\right)\,M_\odot\,\mathrm{km\,s^{-1}}.
$$
Using the SISSI momentum-per-supernova relation
$$
\hat{p}_{\rm SN}\sim 2.6\times10^5\,n_0^{-0.13}\,M_\odot\,\mathrm{km\,s^{-1}},
$$
the paper obtains, for the Local Bubble density,
$$
\hat{p}_{\rm SN,LB}\sim (2.854\pm0.005)\times10^5\,M_\odot\,\mathrm{km\,s^{-1}},
$$
and therefore
$$
N_{\rm SN,LB}(t_{\rm age})\sim (104\pm2)\left(\frac{t_{\rm age}}{\mathrm{Myr}}\right)^{-1}.
$$
The central conclusion is that, **in contrast to previous estimates**, the Local Bubble requires **\(\gtrsim 20\)** SNe over approximately **\(4\,\mathrm{Myr}\)** to explain both its size and momentum [2509.04221].

In the paper’s extrapolated simultaneous solution, the Local Bubble age is
$$
t_{\rm age,LB}\sim 3.5\text{--}5.5\,\mathrm{Myr},
$$
with
$$
N_{\rm SN,LB}\sim 19\text{--}30,
$$
a mean SN interval
$$
\Delta t_{\rm SN}\sim 0.1\text{--}0.3\,\mathrm{Myr},
$$
and an expansion speed
$$
v_{\rm LB}\sim 38\text{--}59\,\mathrm{km\,s^{-1}}.
$$
This substantially revises the older picture of a **\(\gtrsim 14\,\mathrm{Myr}\)**-old bubble powered by roughly **15–20** supernovae [2509.04221].

## 6. Interpretation, limitations, and significance

The revised Local Bubble chronology has direct astrophysical consequences. The second SISSI paper argues that the new age and SN count are in tension with the assumption that the Local Bubble was powered **solely** by supernovae from the nearby **Scorpius-Centaurus** OB association. It instead points to a broader and more recent episode of star formation and supernova activity in the **\(\alpha\)-Persei family**. On that basis, the paper casts serious doubt on the claim that star formation in the solar neighborhood was **triggered** by Local Bubble expansion and argues that the expansion may instead have **quenched** nearby star formation [2509.04221].

At the same time, SISSI is careful about its own limits. The Local Bubble reconstruction depends on identifying the shell as the **first prominent density peak** along each line of sight, which cannot perfectly separate a true shell from unrelated or entrained structures. The central **\(\sim 69\,\mathrm{pc}\)** of the dust map lacks differential extinction information and must be filled by assumption. More importantly, the preferred Local Bubble solution lies partly **outside the parameter space covered directly by the simulations**, so the age estimate relies on extrapolation from the simulated tracks. The first paper also notes that the underlying SNR calculations omit magnetic fields, cosmic rays, thermal conduction, and non-equilibrium chemistry, all of which could affect late-time structure [2503.12977][2509.04221].

Even with those caveats, the project establishes a coherent result: early SNR evolution remains close to classical theory, but late-time morphology is strongly controlled by the anisotropic ISM assembled by galactic dynamics. That conclusion makes geometry an astrophysical diagnostic in its own right. The first paper explicitly states that **future studies targeting the geometry of Galactic SNRs may use this insight to obtain a clearer picture of the processes shaping the Galactic ISM** [2503.12977].

The project also emphasizes reproducibility. The second paper makes the Local Bubble analysis code public in **Julia** at **https://doi.org/10.5281/zenodo.17054923**, which suggests that SISSI is intended as an open comparison framework as well as a simulation suite [2509.04221]. In that sense, SISSI’s significance is twofold: it provides a specific high-resolution sample of supernova remnants in a **stratified, shearing** ISM, and it advances a broader claim that realistic remnant geometry can be used to infer the recent dynamical history of the Galactic environment.

Source: https://www.emergentmind.com/topics/sissi-project