A Novel Parameterization for Rapid Cooling in Supernova Remnants, with applications to the Pa 30 nebula

Abigail Polin; Miranda Pikus; Paul Duffell; Soham Mandal

arxiv: 2511.08519 · v2 · submitted 2025-11-11 · 🌌 astro-ph.HE

A Novel Parameterization for Rapid Cooling in Supernova Remnants, with applications to the Pa 30 nebula

Miranda Pikus , Paul Duffell , Soham Mandal , Abigail Polin This is my paper

Pith reviewed 2026-05-17 23:15 UTC · model grok-4.3

classification 🌌 astro-ph.HE

keywords supernova remnantsrapid coolingfilamentary morphologyPa 30 nebulaRayleigh-Taylor instabilityKelvin-Helmholtz instabilitybeta parameterizationSedov time

0 comments

The pith

A single cooling parameter β produces filamentary structures in supernova remnants like Pa 30 when cooling is faster than 1/400 of the Sedov time.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper introduces a parameterization for rapid cooling in supernova remnants that uses one adjustable parameter called β to set how quickly the thermal energy is radiated away relative to the Sedov expansion time. By running simulations with different values of β the authors map out a sequence of remnant shapes. When β reaches or exceeds 400 the ejecta develops a network of filaments that closely resembles the observed Pa 30 nebula. The filaments arise because rapid cooling lets Rayleigh-Taylor fingers grow while suppressing the Kelvin-Helmholtz mixing that would otherwise smooth them out. In this regime the material in the filaments continues to move at nearly its original free-expansion speed, which in turn implies a modest total explosion energy and a specific cooling luminosity.

Core claim

The paper claims that for β ≳ 400, or when the cooling timescale is shorter than approximately 1/400 of the Sedov time, the ejecta is shaped into a filamentary structure similar to Pa 30. Filament creation is explained by the formation of Rayleigh-Taylor Instability fingers where cooling has prevented the Kelvin-Helmholtz Instability from overturning and mixing out the tips. The ejecta in these filaments have not decelerated and are moving almost completely ballistically at approximately 95-100% their free expansion speed. In this rapid cooling regime an explosion energy of approximately 3.5 × 10^47 erg is inferred, and the cooling mechanism must remove energy at a rate of 2% of E_ej/t, for

What carries the argument

The β parameterization of cooling timescale relative to the Sedov time, which sets the relative strength of Rayleigh-Taylor finger growth versus Kelvin-Helmholtz mixing in the expanding ejecta.

If this is right

Filamentary structures form in supernova remnants once cooling becomes faster than about 1/400 of the Sedov time.
Material inside the filaments continues moving at 95-100 percent of free-expansion speed.
The total explosion energy consistent with such filaments is approximately 3.5 × 10^47 erg.
The implied cooling luminosity required to sustain the structures is roughly 10^36 erg/s.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

This single-parameter approach may allow simpler modeling of other strongly cooling flows such as certain planetary nebulae or galactic cooling flows.
High-precision proper-motion measurements of Pa 30 filaments could directly test whether the material is still moving ballistically.
Comparing the β threshold against full microphysical cooling simulations would show whether one parameter is adequate to capture the instability balance.

Load-bearing premise

That a single tunable parameter β can capture the essential cooling physics sufficiently to control the competition between Rayleigh-Taylor and Kelvin-Helmholtz instabilities without needing detailed radiative transfer or microphysical cooling functions.

What would settle it

A simulation that includes explicit radiative cooling functions and either produces or fails to produce filamentary structures when the cooling time reaches roughly 1/400 of the Sedov time.

Figures

Figures reproduced from arXiv: 2511.08519 by Abigail Polin, Miranda Pikus, Paul Duffell, Soham Mandal.

**Figure 1.** Figure 1: Three-dimensional visualization of the remnant structure for no cooling implemented and for our most rapid cooling with β = 800. Both models are at t = 0.1 tSedov. The yellow-colored iso-surface is a choice of density intended to target the shape of the ejecta, in units of ρCSM. Furthermore, the transparent blue-colored iso-surface represents the pressure near the forward shock displayed in units of P0. Se… view at source ↗

**Figure 2.** Figure 2: Integrated shocked ejecta density along the z axis for our model suite at t = 0.1 tSedov. The forward shock (whitedashed line) and the reverse shock (black-dashed line) are overlaid for clarity and are determined using the conditions in subsection 2.3. The integrated shocked column density is computed by an integration along the line of sight only in the ejecta at radii greater than the reverse shock radi… view at source ↗

**Figure 3.** Figure 3: Cross-section slices of pressure (erg/cm3 ) for four different β values. The side lengths are scaled arbitrary to easily compare between different β values at t = 0.1 tSedov. The forward shock (outer white-dashed line) and the reverse shock (inner white-dashed line) are overlaid for clarity and are determined using the conditions in subsection 2.3. We suggest that these pressure slices demonstrate the cont… view at source ↗

**Figure 4.** Figure 4: Cross-section slices of the homologous expansion fraction k = v/(r/t) for the velocity in the domain at t = 0.1 tSedov. The distribution of k demonstrates that as β increases, or when a more rapid cooling timescale is implemented, the velocity of the material in the filaments gets more ballistic (k ≈ 0.9 − 1.0). Stronger cooling prevents the shocked ejecta from decelerating and the Rayleigh-Taylor instabil… view at source ↗

**Figure 5.** Figure 5: Radially averaged ejecta velocities vej(r) and homologous expansion fraction kej(r) = vej(r)/t for our model suite. Each column denotes a different time. The average is constrained to regions where the ejecta is at least 90% abundant and weighted by the ejecta density in that region. Each column denotes a different time and the color scheme denotes a different β. We remark that as β is increased, the eject… view at source ↗

**Figure 6.** Figure 6: The reverse to forward shock position ratio rRS/rFS with time. The color scheme differentiates between the β runs. The reverse to forward shock position ratio provides a way to dynamically match our model suite with other systems. As expected, SNRs that are cooled off more efficiently evolve faster (i.e. higher β runs show their shock position ratio reaching lower values before the others) since they los… view at source ↗

**Figure 7.** Figure 7: In the left panel, we plot cooling luminosities as a function of time for each β model. The cooling luminosity at each time step is computed by summing the total thermal energy in the system multiplied by a factor of β/tSedov. The cooling luminosity curves rise and fall with a peak value that occurs earlier in time with larger βs. In the right panel, we plot the ratio of the cooling luminosity Lcool to Eej… view at source ↗

**Figure 8.** Figure 8: Scatter plot of selected zones for β = 800 at t = 0.1 tSedov ≈ 844 years. The x-axis is the radius of the selected point and the y-axis is the total velocity for that point. The points are colored by their ballistic fraction the ejecta fraction, or passive scalar value. The forward and reverse shock positions are highlighted by a light grey line, which are found by our estimated tracking algorithm discuss… view at source ↗

**Figure 9.** Figure 9: First row: normalized ejecta density binned by angular coordinates ϕ and θ at t = 0.1 tSedov ≈ 844 yr for β = 800. The darker purple color denotes a higher density region, thus mapping the locations of the filaments in angular space. The black crosses denote the location of the filaments and are plotted again as purple crosses in the left panel of the second row. Second row: locations of filaments generat… view at source ↗

read the original abstract

We systematically study how cooling creates structural changes in supernova remnants as they evolve. Inspired by the peculiar morphology of the Pa 30 nebula, we adopt a framework in which to characterize supernova remnants under different degrees of cooling. Our cooling framework characterizes remnants with a singular parameter called $\beta$ that sets how rapidly the system's thermal energy is radiated or emitted away. A continuum of morphologies is created by the implementation of different cooling timescales. For $\beta \gtrsim 400$, or when the cooling timescale is shorter than $\approx \frac{1}{400}$ of the Sedov time, the ejecta is shaped into a filamentary structure similar to Pa 30. We explain the filament creation by the formation of Rayleigh-Taylor Instability fingers where cooling has prevented the Kelvin-Helmholtz Instability from overturning and mixing out the tips. The ejecta in these filaments have not decelerated and are moving almost completely ballistically at $\approx 95-100\%$ their free expansion speed. In this rapid cooling regime, an explosion energy $\approx 3.5 \times 10^{47}$ erg is inferred. We also propose the cooling mechanism required to create these structures necessitates removing energy at a rate of $2\%$ of $E_{\rm ej}/t$, which implies a cooling luminosity of $\approx 10^{36}$ erg/s.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper's main contribution is a single-parameter β cooling model that produces filamentary SNR morphologies by letting RT fingers grow without KH mixing when cooling is rapid.

read the letter

The key point is that this work introduces a simple β parameter to control cooling timescale relative to the Sedov time, and shows through simulations that β ≳ 400 creates filamentary structures like those in Pa 30. The filaments arise because cooling keeps RT fingers from being mixed away by KH instabilities, leaving material moving ballistically at nearly free-expansion speeds. That link between cooling rate and observable shape is the clearest new piece here, and the simulations do a decent job mapping out a continuum of morphologies as β changes. They also give a straightforward physical story for why the tips survive and stay fast-moving. The paper does well at laying out how different cooling strengths reshape the remnant without needing full radiative transfer in every run. That makes the framework easy to apply to other objects. The soft spots sit in the quantitative parts. The explosion energy of 3.5 × 10^47 erg and the 2% energy-removal rate are obtained by tuning β to match Pa 30's observed speeds and shape, so the numbers are more like a fit than an independent prediction. No error bars or baseline runs against other cooling prescriptions are mentioned, which leaves the strength of those claims unclear. The stress-test concern is fair: a single tunable energy-loss term may not reproduce the exact pressure and density gradients that real microphysical cooling would produce at post-shock conditions, so the instability competition could shift once more detailed functions are used. This paper is aimed at researchers modeling supernova remnant evolution and morphology, especially anyone looking at Pa 30 or similar filamentary nebulae. A reader who wants a practical handle on how cooling affects structure will find the β scans useful. It deserves serious referee time because the parameterization is new and the qualitative morphology results look reproducible from the described setup, even if the numbers need tightening.

Referee Report

2 major / 2 minor

Summary. The paper introduces a single-parameter β to characterize rapid cooling in supernova remnant evolution, where β is the ratio of the cooling timescale to the Sedov time. Hydrodynamic simulations show that for β ≳ 400 the ejecta develops a filamentary morphology resembling the Pa 30 nebula. This is attributed to Rayleigh-Taylor instability fingers whose tips survive because cooling suppresses Kelvin-Helmholtz overturning and mixing; the filaments retain near-ballistic velocities (95–100 % of free-expansion speed). From matching to Pa 30 observations the authors infer an explosion energy of ≈ 3.5 × 10^47 erg and a required cooling rate of 2 % of E_ej/t, corresponding to a luminosity ≈ 10^36 erg s^−1.

Significance. If the β parameterization can be shown to reproduce the essential effects of realistic radiative cooling on instability competition, the approach would supply a computationally inexpensive way to explore cooling-driven morphologies in SNR simulations. The hydrodynamic results that demonstrate the morphological transition at high β constitute a clear strength and provide qualitative support for the filament-formation mechanism. However, the quantitative inferences for explosion energy and cooling luminosity rest on observational tuning without reported uncertainties or baseline comparisons, which reduces the immediate impact of the work.

major comments (2)

[Abstract and results section describing the β ≳ 400 regime] The threshold β ≳ 400 is chosen to reproduce the target filamentary morphology of Pa 30; the subsequent inference of the 2 % energy-removal rate and the 3.5 × 10^47 erg explosion energy then follows from matching simulated ballistic speeds to the observed nebula. This makes the quantitative claims dependent on the same observational tuning that defined the regime, rather than an independent test of the model.
[Section explaining filament creation via RT/KH competition] The central physical claim—that cooling selectively prevents KH overturning of RT fingers while leaving ballistic ejecta—requires that the single-parameter energy-loss term modulate local pressure and density gradients in a manner equivalent to microphysical radiative losses. No growth-rate measurements or side-by-side comparisons with explicit cooling functions at post-shock temperatures and metallicities relevant to Pa 30 are reported.

minor comments (2)

[Methods section introducing the cooling framework] The definition of β as a ratio of timescales is introduced without an explicit equation or comparison table to existing cooling prescriptions (e.g., those based on Sutherland & Dopita or other tabulated cooling curves).
[Figure captions] Figure captions and axis labels should explicitly state the value of β used for each panel and whether the plotted velocities are in the lab frame or normalized to free-expansion speed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments. We address each major comment below and have revised the manuscript to clarify our approach and strengthen the presentation of results.

read point-by-point responses

Referee: [Abstract and results section describing the β ≳ 400 regime] The threshold β ≳ 400 is chosen to reproduce the target filamentary morphology of Pa 30; the subsequent inference of the 2 % energy-removal rate and the 3.5 × 10^47 erg explosion energy then follows from matching simulated ballistic speeds to the observed nebula. This makes the quantitative claims dependent on the same observational tuning that defined the regime, rather than an independent test of the model.

Authors: We agree that the β ≳ 400 regime was identified from the simulations as the point at which filamentary structures resembling Pa 30 appear. The explosion energy is subsequently inferred from the near-ballistic velocities (95–100 % of free-expansion speed) once that regime is selected. This is not fully independent, as the morphology match guides the choice of regime. In the revised manuscript we have added text in the results section explaining how the β threshold emerges from the hydrodynamic behavior across a range of values, independent of the specific Pa 30 observations, and we include a brief sensitivity study showing how the inferred energy varies with β near the threshold. revision: partial
Referee: [Section explaining filament creation via RT/KH competition] The central physical claim—that cooling selectively prevents KH overturning of RT fingers while leaving ballistic ejecta—requires that the single-parameter energy-loss term modulate local pressure and density gradients in a manner equivalent to microphysical radiative losses. No growth-rate measurements or side-by-side comparisons with explicit cooling functions at post-shock temperatures and metallicities relevant to Pa 30 are reported.

Authors: The β parameterization is an effective description intended to isolate the dominant dynamical consequence of rapid cooling on the RT–KH competition. Our simulations demonstrate that sufficiently high β suppresses KH overturning at the RT finger tips while preserving ballistic motion. We acknowledge that quantitative growth-rate measurements and direct comparisons against explicit cooling functions at the relevant temperatures and metallicities are not provided. In the revised manuscript we have inserted a dedicated paragraph in the methods and discussion sections that justifies the parameterization on physical grounds, states its limitations explicitly, and notes that future work will include such comparisons. revision: yes

Circularity Check

1 steps flagged

Tunable β parameter tuned to Pa 30 morphology; inferred energy and cooling rate reduce to data matching

specific steps

fitted input called prediction [Abstract]
"For β ≳ 400, or when the cooling timescale is shorter than ≈1/400 of the Sedov time, the ejecta is shaped into a filamentary structure similar to Pa 30. ... In this rapid cooling regime, an explosion energy ≈ 3.5 × 10^{47} erg is inferred. We also propose the cooling mechanism required to create these structures necessitates removing energy at a rate of 2% of E_ej/t, which implies a cooling luminosity of ≈ 10^{36} erg/s."

β is introduced as a tunable ratio of cooling to Sedov time whose threshold is chosen to reproduce the observed Pa 30 filamentary morphology; the explosion energy and 2% energy-removal rate are subsequently inferred by matching simulated ballistic speeds to the same observed nebula, so the quantitative results are forced by tuning against the target data.

full rationale

The paper introduces β as a single free parameter controlling the cooling timescale relative to the Sedov time. The central claim that β ≳ 400 produces filamentary structures like Pa 30 is obtained by selecting the threshold that reproduces the target morphology in simulations. The quantitative results—an explosion energy of ≈3.5×10^47 erg and a required energy removal rate of 2% of E_ej/t—are then inferred by calibrating simulated ballistic ejecta speeds against the observed Pa 30 nebula. This reduces the 'predictions' to a fit against the same observational data used to define the regime of interest, constituting fitted-input-called-prediction circularity. No independent microphysical validation or external benchmark is reported for the instability competition.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claim rests on introducing β as a free parameter to set cooling timescale and on the domain assumption that standard hydrodynamic instabilities govern structure formation once cooling is added; no new physical entities are postulated.

free parameters (1)

β
Single parameter controlling cooling timescale relative to Sedov time; threshold of 400 selected to produce filamentary morphology matching Pa 30.

axioms (1)

domain assumption Hydrodynamic evolution of supernova remnants with parameterized cooling is governed by the competition between Rayleigh-Taylor and Kelvin-Helmholtz instabilities.
Invoked to explain why rapid cooling produces intact filaments rather than mixed structures.

pith-pipeline@v0.9.0 · 5558 in / 1509 out tokens · 61333 ms · 2026-05-17T23:15:34.021496+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel contradicts

?

contradicts
CONTRADICTS: the theorem conflicts with this paper passage, or marks a claim that would need revision before publication.

Our cooling framework characterizes remnants with a singular parameter called β that sets how rapidly the system’s thermal energy is radiated or emitted away... SE = −β ε / tSedov
IndisputableMonolith/Foundation/DimensionForcing.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

For β ≳ 400... the ejecta is shaped into a filamentary structure

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

5 extracted references · 5 canonical work pages

[1]

2006, The Astrophysical Journal, 645, 1373, doi: 10.1086/504399

Badenes, C., Borkowski, K. J., Hughes, J. P., Hwang, U., & Bravo, E. 2006, The Astrophysical Journal, 645, 1373, doi: 10.1086/504399 Beloborodov, A. M. 2014, Monthly Notices of the Royal Astronomical Society, 438, 169, doi: 10.1093/mnras/stt2140 Bietenholz, M. F. 2006, The Astrophysical Journal, 645, 1180, doi: 10.1086/504584 Chevalier, R. A., Blondin, J....

work page doi:10.1086/504399 2006
[2]

M., Wright, E

https://baas.aas.org/pub/2025n2i322p04/release/1 Cutri, R. M., Wright, E. L., Conrow, T., et al. 2012, Explanatory Supplement to the WISE All-Sky Data Release Products. https://ui.adsabs.harvard.edu/abs/2012wise.rept....1C Dubner, G., Castelletti, G., Kargaltsev, O., et al. 2017, The Astrophysical Journal, 840, 82, doi: 10.3847/1538-4357/aa6983 Duffell, P...

work page doi:10.3847/1538-4357/aa6983 2012
[3]

2018, The Astrophysical Journal, 869, 140, doi: 10.3847/1538-4357/aaedb7 Ko, T., Suzuki, H., Kashiyama, K., et al

https://ui.adsabs.harvard.edu/abs/ 1957SCoA....2..109H/abstract Kashyap, R., Haque, T., Lor´ en-Aguilar, P., Garc´ ıa-Berro, E., & Fisher, R. 2018, The Astrophysical Journal, 869, 140, doi: 10.3847/1538-4357/aaedb7 Ko, T., Suzuki, H., Kashiyama, K., et al. 2024, The Astrophysical Journal, 969, 116, doi: 10.3847/1538-4357/ad4d99 Kothes, R. 2013, Astronomy ...

work page doi:10.3847/1538-4357/aaedb7 2018
[4]

edu/abs/1971QJRAS..12...10S/abstract Stephenson, F

https://ui.adsabs.harvard. edu/abs/1971QJRAS..12...10S/abstract Stephenson, F. R., & Green, D. A. 2002, International Series in Astronomy and Astrophysics,

work page 2002
[5]

1990, The Astrophysical Journal, 357, 138, doi: 10.1086/168899 Vaytet, N

https: //ui.adsabs.harvard.edu/abs/2002ISAA....5.....S/abstract Van Den Bergh, S. 1990, The Astrophysical Journal, 357, 138, doi: 10.1086/168899 Vaytet, N. M. H., O’Brien, T. J., & Rushton, A. P. 2007, Monthly Notices of the Royal Astronomical Society, 380, 175, doi: 10.1111/j.1365-2966.2007.12057.x Warren, J. S., Hughes, J. P., Badenes, C., et al. 2005, ...

work page doi:10.1086/168899 1990

[1] [1]

2006, The Astrophysical Journal, 645, 1373, doi: 10.1086/504399

Badenes, C., Borkowski, K. J., Hughes, J. P., Hwang, U., & Bravo, E. 2006, The Astrophysical Journal, 645, 1373, doi: 10.1086/504399 Beloborodov, A. M. 2014, Monthly Notices of the Royal Astronomical Society, 438, 169, doi: 10.1093/mnras/stt2140 Bietenholz, M. F. 2006, The Astrophysical Journal, 645, 1180, doi: 10.1086/504584 Chevalier, R. A., Blondin, J....

work page doi:10.1086/504399 2006

[2] [2]

M., Wright, E

https://baas.aas.org/pub/2025n2i322p04/release/1 Cutri, R. M., Wright, E. L., Conrow, T., et al. 2012, Explanatory Supplement to the WISE All-Sky Data Release Products. https://ui.adsabs.harvard.edu/abs/2012wise.rept....1C Dubner, G., Castelletti, G., Kargaltsev, O., et al. 2017, The Astrophysical Journal, 840, 82, doi: 10.3847/1538-4357/aa6983 Duffell, P...

work page doi:10.3847/1538-4357/aa6983 2012

[3] [3]

2018, The Astrophysical Journal, 869, 140, doi: 10.3847/1538-4357/aaedb7 Ko, T., Suzuki, H., Kashiyama, K., et al

https://ui.adsabs.harvard.edu/abs/ 1957SCoA....2..109H/abstract Kashyap, R., Haque, T., Lor´ en-Aguilar, P., Garc´ ıa-Berro, E., & Fisher, R. 2018, The Astrophysical Journal, 869, 140, doi: 10.3847/1538-4357/aaedb7 Ko, T., Suzuki, H., Kashiyama, K., et al. 2024, The Astrophysical Journal, 969, 116, doi: 10.3847/1538-4357/ad4d99 Kothes, R. 2013, Astronomy ...

work page doi:10.3847/1538-4357/aaedb7 2018

[4] [4]

edu/abs/1971QJRAS..12...10S/abstract Stephenson, F

https://ui.adsabs.harvard. edu/abs/1971QJRAS..12...10S/abstract Stephenson, F. R., & Green, D. A. 2002, International Series in Astronomy and Astrophysics,

work page 2002

[5] [5]

1990, The Astrophysical Journal, 357, 138, doi: 10.1086/168899 Vaytet, N

https: //ui.adsabs.harvard.edu/abs/2002ISAA....5.....S/abstract Van Den Bergh, S. 1990, The Astrophysical Journal, 357, 138, doi: 10.1086/168899 Vaytet, N. M. H., O’Brien, T. J., & Rushton, A. P. 2007, Monthly Notices of the Royal Astronomical Society, 380, 175, doi: 10.1111/j.1365-2966.2007.12057.x Warren, J. S., Hughes, J. P., Badenes, C., et al. 2005, ...

work page doi:10.1086/168899 1990