Cold Dense Shell (CDS) in Supernovae

• Cold Dense Shell (CDS) is the radiatively cooled, high-density layer formed between the reverse and forward shocks in supernova interactions.
• It is characterized by rapid cooling, fragmentation, and conditions favorable for dust condensation, as demonstrated in SN 2005ip studies.
• The CDS acts as a kinematic intermediary that affects spectral diagnostics and explosion energy estimates in interacting supernovae.

The Cold Dense Shell (CDS) is the radiatively cooled, compressed interaction layer that forms in a supernova when freely expanding ejecta collide with circumstellar material (CSM). In the standard double-shock structure, the reverse shock propagates back into the ejecta in the Lagrangian sense, the forward shock propagates outward into the CSM, and the CDS lies between them around the contact discontinuity. In the literature summarized here, the CDS is not merely a geometric byproduct of interaction: it is treated as a dynamically important shell whose density, cooling state, fragmentation, and optical depth can control line formation, dust condensation, and even the inferred explosion energy of an interacting event (Shahbandeh et al., 2024, Chugai, 2020).

1. Physical definition and shock geometry

In interacting supernovae, the CDS occupies the region between the reverse shock and the forward shock. The forward shock processes the pre-existing circumstellar gas, the reverse shock processes ejecta, and the contact discontinuity is the interface between those two shocked zones. When post-shock cooling is efficient, the shocked gas loses thermal support and collapses into a thin, compressed layer rather than remaining hot and extended; that cooled compressed layer is the CDS (Shahbandeh et al., 2024, Chugai, 2020).

The terminology is specific. The shell is “dense” because shock compression combined with radiative cooling raises the post-shock density strongly. It is “cold” not in an everyday thermodynamic sense, but because it is much cooler than the immediate hot post-shock gas. In this usage, “cold” means radiatively cooled. The shell can contain shocked ejecta and shocked CSM accumulated in the interaction zone, and in analytic treatments it is often represented with the thin-shell approximation (Chugai, 2020).

The geometric placement of the CDS is essential for interpretation. In the SN 2005ip analysis, the newly detected cold dust is inferred to lie “behind the forward shock,” that is, in the compressed post-shock shell. In the SN 2013fs analysis, the CDS is the shell whose velocity is inferred from the broad He II $4686$ Å profile. These two cases illustrate complementary uses of the same structure: one as a condensation site for dust and the other as a kinematic tracer of ejecta–CSM interaction (Shahbandeh et al., 2024, Chugai, 2020).

2. Formation conditions, cooling, and instability

A CDS forms when the ejecta–CSM interaction is sufficiently strong and the shocks are radiative. In that regime, the post-shock gas cools efficiently enough to collapse into a geometrically thin, high-density layer. This is why both the dust analysis of SN 2005ip and the dynamical analysis of SN 2013fs treat the shell as a natural consequence of dense circumstellar interaction rather than as an optional phenomenological component (Shahbandeh et al., 2024, Chugai, 2020).

The SN 2013fs model provides an explicit example of the surrounding environment that can generate such a shell. The explosion is assumed to occur inside a confined dense circumstellar shell of radial extent

$R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$

with approximately uniform density inside that radius, then

$$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$

and beyond

$r>5\times10^{15}\ \mathrm{cm},$

a normal red-supergiant wind with

$$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$

The ejecta are taken to expand homologously with density profile

$\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$

where $q=7.6$ is adopted from the hydrodynamic model of SN 2008in (Chugai, 2020).

The CDS is not necessarily smooth. In SN 2013fs, the shell is argued to be decelerating and therefore subject to Rayleigh–Taylor instability. The paper states that “the decelerating CDS is subject to the Rayleigh-Taylor instability,” with fragmentation into filaments and folded sheets. That instability is not a peripheral detail: it becomes central to the line-profile interpretation, because the observed broad He II profile cannot be reproduced by either a smooth optically thin shell or a smooth optically thick shell (Chugai, 2020).

A recurring implication is that the CDS is favored whenever high density and efficient cooling coexist. In the SN 2005ip discussion, this is taken further: the cooled compressed shell is presented as a favorable dust-formation site because high density raises collision and nucleation rates, radiative cooling lowers temperatures into the condensation regime, and the shell can provide shielding from destructive radiation fields. The paper explicitly notes that in the CDS “all dust compositions can then recondense,” in contrast with the pre-shock CSM, where silicates may be more easily destroyed near the shock because their vaporization temperature is lower than that of graphite (Shahbandeh et al., 2024).

3. Spectroscopic and radiative diagnostics

The CDS is often inferred indirectly from observables whose morphology depends on where emission forms and how the shell is structured. In SN 2013fs, the decisive diagnostic is the broad, low-contrast emission near $4500$ Å on day $2.42$, identified as He II $4686$ Å. An origin in unshocked ejecta ionized by reverse-shock X-rays is rejected because the emitting layer would be too thin: $R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$0 with

$R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$1

and homologous expansion would then imply

$R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$2

Such a layer should produce a narrow-in-velocity, boxy profile with blue skew from occultation, whereas the observed line is dome-like (Chugai, 2020).

A smooth spherical thin shell also fails. In the SN 2013fs treatment, an optically thick smooth shell gives an M-shaped profile, while an optically thin smooth shell gives a boxy profile. The observed line instead requires a fragmented CDS, parameterized by the area ratio

$R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$3

where $R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$4 is the cumulative surface area of fragments and $R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$5 is the shell radius. For $R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$6, corresponding to strong mixing and a highly folded fragmented shell, the line profile becomes parabolic; for $R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$7, it remains boxy. The dome-like He II feature therefore implies a fragmented, well-mixed CDS with effectively $R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$8. With a continuum optical depth

$R_{\rm ds}\sim 5\times 10^{14}\ \mathrm{cm},$9

in the shell and circumstellar Thomson optical depth

$$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$0

the best fit yields

$$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$1

at day $$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$2 (Chugai, 2020).

In SN 2005ip, the radiative diagnostics are dominated by dust emission rather than a single optical line. The dust-mass analysis uses

$$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$3

where the escape probability for a spherical geometry is

$$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$4

and the optical depth is

$$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$5

These expressions matter because a dense region such as the CDS can be optically thick at some wavelengths, so one cannot assume that all dust emission escapes freely. The authors note that scenarios “such as the CDS” are better represented by a shell rather than a filled sphere, but they retain the spherical expressions because they include photon trapping and therefore provide robust estimates relative to a purely optically thin approximation (Shahbandeh et al., 2024).

The same paper also emphasizes prior line-profile evidence. Earlier optical work on SN 2005ip had found blue-shifted asymmetries in both intermediate-width and broad components. The intermediate-width asymmetry was associated with the post-shock interaction region and therefore favored CDS dust, while broad-component asymmetry suggested some ejecta dust. The mid-IR data do not replace that line-profile context; they refine it by locating the dominant late-time dust component more directly (Shahbandeh et al., 2024).

4. The CDS as a dust-formation site in SN 2005ip

The SN 2005ip study treats the CDS as the most plausible site of a large, newly formed dust reservoir detected nearly $$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$6 years after explosion. The data combine a Spitzer IRS spectrum obtained in 2008, about $$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$7 days after discovery, covering $$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$8, with a JWST/MIRI MRS spectrum obtained on 2023 April 20, $$\rho \propto r^{-3}\qquad (r>5\times10^{14}\ \mathrm{cm}),$$9 days after discovery, with $r>5\times10^{15}\ \mathrm{cm},$0 spectroscopy covering $r>5\times10^{15}\ \mathrm{cm},$1. The long baseline is central to the interpretation because a single epoch can accommodate multiple geometries, whereas two mid-IR spectra separated by roughly $r>5\times10^{15}\ \mathrm{cm},$2 years reveal genuine evolution in the dust population (Shahbandeh et al., 2024).

The most important new result is the appearance of a dominant cool silicate component in the JWST spectrum. The earlier Spitzer spectrum was featureless and had “no obvious silicate features,” while the JWST data display dominant Mg-silicate features plus a shorter-wavelength, featureless warm continuum attributed to amorphous carbon. The paper considers two possibilities for the absence of silicate features at the earlier epoch: either the silicate component was not yet present, or it was already present but suppressed by optical thickness. The favored interpretation is that the cold component was earlier more compact and optically thick, then expanded, became optically thinner, and continued forming dust (Shahbandeh et al., 2024).

The fitted properties are explicit. For the Spitzer epoch, the cool component is Mg-silicate with

$r>5\times10^{15}\ \mathrm{cm},$3

$r>5\times10^{15}\ \mathrm{cm},$4

$r>5\times10^{15}\ \mathrm{cm},$5

The warm component is amorphous carbon with

$r>5\times10^{15}\ \mathrm{cm},$6

$r>5\times10^{15}\ \mathrm{cm},$7

and a radius of order a few $r>5\times10^{15}\ \mathrm{cm},$8. For the JWST epoch, the cool Mg-silicate component is

$r>5\times10^{15}\ \mathrm{cm},$9

$$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$0

with fitted

$$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$1

but because the component is optically thin the authors instead emphasize a blackbody lower-limit radius

$$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$2

The warm amorphous-carbon component in the JWST spectrum has

$$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$3

$$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$4

$$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$5

The authors’ interpretation is correspondingly specific: the dominant cold Mg-silicate component is attributed to the CDS, while the smaller warm carbonaceous component is attributed either to ejecta dust or possibly to dust near the shock in the pre-shocked CSM (Shahbandeh et al., 2024).

The physical argument for the CDS as a dust factory rests on the conditions in the shell. The paper quotes approximate sublimation or vaporization temperatures of $$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$6 K for silicates and $$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$7 K for graphites. This helps explain why a warm carbonaceous component can persist near an active shock while the large silicate component is inferred to condense farther into the cooled post-shock shell. The shell’s high density and cooling are presented as conditions that ordinary expanding ejecta do not always maintain long enough or at high enough density (Shahbandeh et al., 2024).

5. Spatial localization, temporal evolution, and heating

The CDS interpretation in SN 2005ip is not based on composition alone; it depends on a spatial argument. The modeled cool dust radii are described as “too high to be consistent with the ejecta radii,” so the bulk of the cold component cannot be in the freely expanding ejecta. At the same time, the updated models place the cool dust “within (not external to) the forward shock radius,” which argues against a distant, static pre-existing CSM shell. Outside the ejecta scale but inside the forward shock scale is precisely the spatial domain expected for the CDS (Shahbandeh et al., 2024).

Temporal evolution strengthens that inference. The cold component changes from a compact, featureless or feature-suppressed component at about $$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$8 K in the 2008 Spitzer spectrum to a much cooler, clearly silicate-rich component at about $$\rho \propto r^{-2}, \qquad w=\dot M/u_w \sim 10^{15}\ \mathrm{g\,cm^{-1}}.$$9 K in the 2023 JWST spectrum. By contrast, the warm carbon component cools from about $\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$0 K to $\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$1 K while keeping roughly the same low mass, $\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$2. The dominant late-time mass is therefore carried by the cold component, with

$\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$3

The paper summarizes the result as a “new high-mass dust component ($\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$4) that is not present in the earlier Spitzer spectrum” (Shahbandeh et al., 2024).

The optical-depth evolution is a key technical point. The cold component is inferred to evolve from optically thick at the Spitzer epoch, with average

$\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$5

to optically thin by the JWST epoch. Since Equation (3) shows that $\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$6 decreases as the dust distribution expands, this implies a substantial growth in characteristic radius. The lower limit on the cold-dust radius increases by about a factor of $\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$7, consistent with the transition from optically thick to optically thin. The preferred interpretation is therefore not simply cooling, but continued CDS dust formation combined with expansion: “as the CDS expanded, the dust shell likely evolved from optically thick to optically thin while also continuing to form new dust” (Shahbandeh et al., 2024).

The heating mechanism is treated as another discriminator between newly formed CDS dust and a distant IR echo. SN 2005ip remains an actively interacting supernova nearly $\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$8 years after explosion. The late-time optical spectra correspond to luminosities of about $\rho(v)=\frac{\rho_0}{1+(v/v_0)^q},$9, and strong H$q=7.6$0 emission indicates sustained shock interaction. The forward shock can generate shock power of order

$q=7.6$1

with much of the thermalized emission emerging in the UV and only a few percent in the optical; Dessart & Hillier models are cited as showing that an observed optical luminosity of $q=7.6$2 can be compatible with a UV luminosity of $q=7.6$3. The dust is therefore interpreted as being radiatively heated by ongoing interaction rather than by radioactive decay, which would be inadequate by day $q=7.6$4, or by a one-time light echo from the original explosion (Shahbandeh et al., 2024).

6. The CDS as a kinematic intermediary in SN 2013fs

In SN 2013fs, the CDS serves a different but equally important role. Chugai’s reinterpretation argues that the broad early He II line traces the interaction-driven shell rather than freely expanding outer ejecta, and that this distinction lowers the inferred explosion energy. The preferred model has

$q=7.6$5

with confined-shell mass

$q=7.6$6

The alternative higher-energy model has

$q=7.6$7

with

$q=7.6$8

and the energy–mass scaling used for fixed CSM density is

$q=7.6$9

The uncertainty range favored by the light-curve comparison is

$4500$0

This revised energy estimate contrasts with the earlier value

$4500$1

for $4500$2, which had appeared to exceed the usual neutrino-mechanism upper limit of about

$4500$3

By reassigning the early fast material to the CDS rather than directly to unshocked ejecta, the paper places SN 2013fs back within the range of neutrino-driven core-collapse explosions (Chugai, 2020).

The logic combines several observables. The first is the CDS velocity from the broad He II line at day $4500$4: $4500$5 The second is the maximum velocity of unshocked ejecta from H$4500$6 absorption on day $4500$7: $4500$8 The third is the early photospheric radius, assumed to match the CDS radius because at early times the photosphere is taken to coincide essentially with the inner boundary of the shell. These constraints are evolved within the thin-shell approximation (Chugai, 2020).

Kinematics alone are insufficient because the paper explicitly states that a change in energy can be compensated by a corresponding change in CSM density so that the CDS velocity at a given age is preserved. The degeneracy is broken with the observed H$4500$9 luminosity from the photoionized CSM at day $2.42$0: $2.42$1 For case C recombination,

$2.42$2

using $2.42$3 K. Model A gives

$2.42$4

while model B gives

$2.42$5

Model A matches the observed H$2.42$6 luminosity, whereas model B overpredicts it by a factor of $2.42$7. The CDS is thus the observable intermediary that links the early line profile to the revised explosion-energy estimate (Chugai, 2020).

7. Conceptual distinctions, caveats, and broader significance

Several distinctions are necessary to avoid conflating the CDS with adjacent regions. Ejecta dust refers to dust condensed in the freely expanding supernova ejecta interior to the reverse shock. Circumstellar dust or pre-existing dust refers to dust formed before explosion in the progenitor’s wind or eruptive mass loss, exterior to the forward shock. CDS dust lies between the shocks in the cooled compressed interaction shell. In SN 2005ip, the dominant cool Mg-silicate reservoir is assigned to the CDS, while the smaller warm carbonaceous component remains only tentatively assigned to ejecta or possibly to dust near the forward shock in the pre-shocked CSM (Shahbandeh et al., 2024).

A common misconception is to interpret “cold” literally. In both papers, “cold” means radiatively cooled relative to the immediate hot shocked gas, not absolutely cold material in isolation. Another misconception is that any broad interaction-powered line must measure freely expanding ejecta. The SN 2013fs case shows that a broad profile can instead trace a fragmented CDS, and that this reinterpretation can materially change the inferred explosion energy (Chugai, 2020).

The limitations are substantial and explicitly acknowledged. For SN 2005ip, the dust modeling assumes a homogeneous, spherically symmetric dust distribution even though the true CDS is a shell and may have a temperature gradient, nonuniform density, and asymmetric geometry; the exact heating source remains uncertain; the location of the warm carbonaceous component is not definitively established; and the longest-wavelength JWST channel is downweighted above $2.42$8. For SN 2013fs, the line identification and line-formation region are model-dependent; the fragmented-shell requirement with $2.42$9 is physically motivated but not directly observed; the interaction is modeled spherically; the photosphere is assumed to coincide with the CDS inner boundary; the thin-shell approximation replaces full radiation hydrodynamics; and the dominant structural uncertainty is the degeneracy between explosion energy and CSM density (Shahbandeh et al., 2024, Chugai, 2020).

Within those limits, the broader significance is twofold. First, the SN 2005ip result suggests that Type IIn supernovae can form large dust masses in the CDS, offering what the paper calls an “alternative to the SN ejecta” as a dust-production channel. The inferred dust mass of about $4686$0 approaches the $4686$1 scale relevant to discussions of rapid dust enrichment, although the paper does not claim that Type IIn events alone solve the high-redshift dust problem. Second, the SN 2013fs analysis shows that the CDS can function as a dynamical observable that reduces apparent tension between data and the neutrino-driven explosion paradigm. Taken together, these studies establish the CDS as both a physically consequential structure and an interpretive pivot in the analysis of interacting supernovae (Shahbandeh et al., 2024, Chugai, 2020).