Destruction of the interstellar dust by a supernova
Pith reviewed 2026-05-16 19:13 UTC · model grok-4.3
The pith
Supernovae destroy several times less interstellar dust in dense high-star-formation regions than in the Milky Way diffuse medium.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper argues that interstellar dust destruction in supernova remnants proceeds in two distinct regimes determined by the ambient medium's density: rapid and almost complete destruction occurs in low-mass compact remnants expanding into dense gas, while gradual and weak destruction takes place in massive remnants evolving in low-density environments. The mass of destroyed dust is maximized when the time available for thermal sputtering equals the dynamical age of the remnant. Increasing the ambient gas density reduces the destroyed dust mass only logarithmically, dust cooling suppresses the destruction by a factor of up to 1.6, and the destroyed dust mass depends linearly on gas metadata.
What carries the argument
The two-regime model of dust grain destruction in supernova remnants, where the comparison between thermal sputtering time and remnant dynamical time determines the efficiency, modulated by ambient density and metallicity.
Load-bearing premise
The modeling assumes thermal sputtering rates and dust cooling efficiencies can be applied uniformly without detailed grain-size distribution evolution or magnetic field effects, and that the two-regime classification applies across all observed supernova environments.
What would settle it
A measurement showing that the dust mass destroyed per supernova in a high-density environment like an ultraluminous infrared galaxy is comparable to or exceeds that in the Milky Way diffuse medium would falsify the claim of several times smaller destruction.
Figures
read the original abstract
Destruction of the interstellar dust proceeds primary behind supernova shocks. The previous estimates of the mass of the interstellar dust destroyed in the SN remnant do not take into account the physical properties of the ambient medium. Here we consider how some parameters, i.e. gas density and metallicity, can influence the destruction of the interstellar dust. We show that there are two regimes of the interstellar dust grains destruction in SN remnants: rapid and almost complete in compact low-mass SN remnants expanding in dense medium, and gradual and weak destruction in massive remnants evolving in the low-dense environment. When time for thermal sputtering is close to the dynamical one, i.e. to the SN remnant age, the mass of the interstellar dust destroyed in the SN remnant reaches its maximum value. We find that change of the ambient gas density results in the reduction of the dust mass logarithmically. We argue that dust cooling suppresses the interstellar dust destruction up to a factor of 1.6 by mass. This factor decreases for higher density of the ambient medium. We found that the dust mass depends linearly on gas metallicity as ${\rm log}~M_d \sim {\rm [Z/H]}$ or, in other words, on the dust-to-gas ratio as $M_d \sim \zeta_d$. In turn, the destruction efficiency is higher in low-metallicity environments due to relatively longer adiabatic phase. We point out that the mass of the interstellar dust destroyed per one SN in a high density environment of the high star formation regions like in local ultraluminous infrared and high-redshift massive galaxies is about several times smaller than that in the Milky Way diffuse medium.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript examines the destruction of interstellar dust in supernova remnants, accounting for ambient gas density and metallicity. It delineates two destruction regimes: rapid and nearly complete destruction in compact remnants expanding in dense media, versus gradual and weak destruction in massive remnants in low-density environments. The maximum destroyed dust mass occurs when the thermal sputtering timescale matches the dynamical age of the remnant. The destroyed mass is found to decrease logarithmically with increasing ambient density, with dust cooling providing an additional suppression factor of up to 1.6 (decreasing at higher densities). The dust mass scales linearly with gas metallicity. The authors conclude that the dust mass destroyed per supernova in high-density environments typical of ultraluminous infrared galaxies and high-redshift massive galaxies is several times smaller than in the Milky Way's diffuse medium.
Significance. If the scalings hold after explicit derivation, the result would refine dust survival estimates in dense star-forming regions, implying higher net dust retention in ULIRGs and high-redshift galaxies than Milky Way-calibrated models predict. The two-regime distinction offers a useful organizing principle for incorporating environmental effects into galactic chemical evolution codes.
major comments (3)
- [Abstract and regime discussion] The central claim that destroyed dust mass is 'several times smaller' in high-density environments rests on the stated logarithmic density dependence and the ~1.6 cooling suppression, yet no explicit equations for t_sput, remnant radius evolution, or the numerical evaluation of the multiplier appear; without these the factor cannot be verified independently.
- [Modeling assumptions] The modeling applies fixed thermal sputtering rates uniformly across remnant phases without integrating over an evolving grain-size distribution; in dense media the remnant enters the sputtering regime at smaller radius, so preferential removal of small grains could flatten the claimed logarithmic density scaling or change the 'several times' factor.
- [Metallicity dependence] The linear metallicity scaling log M_d ~ [Z/H] is tied to a longer adiabatic phase at low Z, but the quantitative link between metallicity, cooling, and the duration of the adiabatic phase is not shown with specific remnant models or parameter values.
minor comments (2)
- [Abstract] The abstract states the 1.6 suppression factor decreases with density but provides no table or plot showing the density dependence of this factor.
- [Introduction] Notation for dust-to-gas ratio (ζ_d) and metallicity ([Z/H]) should be defined at first use with a brief reminder of their relation.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments. We address each major point below and have revised the manuscript to improve clarity and verifiability where possible.
read point-by-point responses
-
Referee: [Abstract and regime discussion] The central claim that destroyed dust mass is 'several times smaller' in high-density environments rests on the stated logarithmic density dependence and the ~1.6 cooling suppression, yet no explicit equations for t_sput, remnant radius evolution, or the numerical evaluation of the multiplier appear; without these the factor cannot be verified independently.
Authors: We agree that the absence of explicit equations limits independent verification. In the revised manuscript we have inserted the thermal sputtering timescale formula t_sput = (a / 0.1 μm) * (n / 1 cm^{-3})^{-1} * (T / 10^6 K)^{-1/2} yr (with appropriate constants), the Sedov-Taylor radius evolution R(t) = (E t^2 / ρ_0)^{1/5}, and a short appendix section that walks through the numerical integration of destroyed mass over the remnant lifetime, showing how the logarithmic density dependence and the factor of ~1.6 arise for the fiducial parameters. revision: yes
-
Referee: [Modeling assumptions] The modeling applies fixed thermal sputtering rates uniformly across remnant phases without integrating over an evolving grain-size distribution; in dense media the remnant enters the sputtering regime at smaller radius, so preferential removal of small grains could flatten the claimed logarithmic density scaling or change the 'several times' factor.
Authors: The referee correctly identifies a modeling simplification. Our calculation adopts effective, size-averaged sputtering rates calibrated to standard MRN distributions. A full time-dependent grain-size integration would be more complete and could modestly alter the numerical prefactor in dense media. We have added a paragraph in the methods section acknowledging this limitation and noting that the dominant logarithmic density scaling is set by the dynamical time available for sputtering rather than by the detailed size evolution; we retain the original scaling but flag the possible quantitative adjustment. revision: partial
-
Referee: [Metallicity dependence] The linear metallicity scaling log M_d ~ [Z/H] is tied to a longer adiabatic phase at low Z, but the quantitative link between metallicity, cooling, and the duration of the adiabatic phase is not shown with specific remnant models or parameter values.
Authors: We accept that the quantitative connection was not demonstrated explicitly. The revised text now includes the metallicity-dependent cooling function Λ(Z,T) used, the resulting adiabatic-phase duration t_ad ∝ Z^{-0.4} for the relevant temperature range, and two example remnant models (Z = 0.1 Z_⊙ and Z = Z_⊙) that illustrate how the longer adiabatic phase at low metallicity increases the integrated sputtering time and produces the reported linear scaling. revision: yes
Circularity Check
No significant circularity; scalings derived from time-equality and cooling efficiency without self-referential reduction
full rationale
The paper's central results follow from equating thermal sputtering time to remnant dynamical age to locate the maximum destroyed dust mass, then showing logarithmic density dependence and a ~1.6 suppression from dust cooling. These steps are presented as direct consequences of the adopted sputtering and cooling rates applied to the two-regime classification, without the maximum mass or the logarithmic factor being defined in terms of the output itself. No self-citation is invoked to justify uniqueness or to import an ansatz; the metallicity linearity is stated as a direct proportionality to dust-to-gas ratio. The derivation chain remains self-contained against external benchmarks and does not reduce any prediction to a fitted input by construction.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Thermal sputtering is the dominant destruction mechanism behind supernova shocks and its timescale can be directly compared to the remnant dynamical age.
- domain assumption Dust cooling suppresses destruction by a fixed factor of up to 1.6 that decreases with ambient density.
Reference graph
Works this paper leans on
-
[1]
" write newline "" before.all 'output.state := FUNCTION fin.entry write newline FUNCTION new.block output.state before.all = 'skip after.block 'output.state := if FUNCTION new.sentence output.state after.block = 'skip output.state before.all = 'skip after.sentence 'output.state := if if FUNCTION not #0 #1 if FUNCTION and 'skip pop #0 if FUNCTION or pop #1...
- [2]
- [3]
-
[4]
Blinnikov S. I., Imshennik V. S., Utrobin V. P., 1982, Soviet Astronomy Letters, 8, 361
work page 1982
- [5]
- [6]
- [7]
- [8]
- [9]
-
[10]
Decarli R., Pensabene A., Venemans B., Walter F., Ba \ n ados E., Bertoldi F., Carilli C. L., Cox P., Fan X., Farina E. P., Ferkinhoff C., Groves B. A., Li J., Mazzucchelli C., Neri R., Riechers D. A., Uzgil B., Wang F., Wang R., Weiss A., Winters J. M., Yang J., 2022, , 662, A60
work page 2022
-
[11]
Decarli R., Walter F., Venemans B. P., Ba \ n ados E. u., Bertoldi F., Carilli C., Fan X., Farina E. P., Mazzucchelli C., Riechers D., Rix H.-W., Strauss M. A., Wa ng R., Yang Y., 2018, , 854, 97
work page 2018
- [12]
- [13]
-
[14]
T., 2009, in Henning T., Gr \"u n E., Steinacker J., eds, Cosmic Dust - Near and Far Vol
Draine B. T., 2009, in Henning T., Gr \"u n E., Steinacker J., eds, Cosmic Dust - Near and Far Vol. 414 of Astronomical Society of the Pacific Conference Series, Interstellar Dust Models and Evolutionary Implications . p. 453
work page 2009
- [15]
- [16]
- [17]
- [18]
-
[19]
Dwek E., 1986, , 302, 363
work page 1986
- [20]
- [21]
- [22]
- [23]
-
[24]
Ferrara A., Peroux C., 2021, , 503, 4537
work page 2021
- [25]
-
[26]
Ginolfi M., Graziani L., Schneider R., Marassi S., Valiante R., Dell'Agli F., Ventura P., Hunt L. K., 2018, , 473, 4538
work page 2018
-
[27]
Gnat O., Sternberg A., 2007, , 168, 213
work page 2007
- [28]
-
[29]
Hirashita H., Yan H., 2009, , 394, 1061
work page 2009
- [30]
-
[31]
Huang Z. P., Thuan T. X., Chevalier R. A., Condon J. J., Yin Q. F., 1994, , 424, 114
work page 1994
-
[32]
Jones A. P., Tielens A. G. G. M., 1994, in Tielens A. G. G. M., ed., The Diffuse Interstellar Bands The formation of small grains in shocks in the ISM . pp 79--83
work page 1994
- [33]
-
[34]
Jones A. P., Tielens A. G. G. M., Hollenbach D. J., McKee C. F., 1994, , 433, 797
work page 1994
- [35]
- [36]
-
[37]
Klingenberg C., Schmidt W., Waagan K., 2007, Journal of Computational Physics, 227, 12
work page 2007
- [38]
- [39]
-
[40]
Mart \' nez-Gonz \'a lez S., W \"u nsch R., Silich S., Tenorio-Tagle G., Palou s J., Ferrara A., 2019, , 887, 198
work page 2019
- [41]
-
[42]
Mattsson L., 2021, Research Notes of the American Astronomical Society, 5, 288
work page 2021
-
[43]
A., Villarroel B., 2019, , 483, 5623
Mattsson L., Bhatnagar A., Gent F. A., Villarroel B., 2019, , 483, 5623
work page 2019
-
[44]
Mattsson L., Hedvall R., 2022, , 509, 3660
work page 2022
-
[45]
McKee C., 1989, in Allamandola L. J., Tielens A. G. G. M., eds, Interstellar Dust Vol. 135, Dust Destruction in the Interstellar Medium . p. 431
work page 1989
-
[46]
R., Matsuura M., Sarangi A., 2018, , 214, 53
Micelotta E. R., Matsuura M., Sarangi A., 2018, , 214, 53
work page 2018
-
[47]
Micha owski M. J., Murphy E. J., Hjorth J., Watson D., Gall C., Dunlop J. S., 2010, , 522, A15
work page 2010
-
[48]
J., Watson D., Hjorth J., 2010, , 712, 942
Micha owski M. J., Watson D., Hjorth J., 2010, , 712, 942
work page 2010
-
[49]
Mignone A., Flock M., Vaidya B., 2019, , 244, 38
work page 2019
- [50]
-
[51]
Moseley E. R., Teyssier R., Draine B. T., 2023, , 518, 2825
work page 2023
- [52]
-
[53]
Nozawa T., Kozasa T., Habe A., 2006, , 648, 435
work page 2006
-
[54]
Oort J. H., 1951, in Problems of Cosmical Aerodynamics Interaction of Nova and Supernova Shells With the Interstellar Medium . p. 118
work page 1951
-
[55]
Ostriker J., Silk J., 1973, , 184, L113
work page 1973
- [56]
-
[57]
Rowlands K., Gomez H. L., Dunne L., Arag \'o n-Salamanca A., Dye S., Maddox S., da Cunha E., van der Werf P., 2014, , 441, 1040
work page 2014
-
[58]
Sakamoto K., Wang J., Wiedner M. C., Wang Z., Peck A. B., Zhang Q., Petitpas G. R., Ho P. T. P., Wilner D. J., 2008, , 684, 957
work page 2008
-
[59]
Santini P., Maiolino R., Magnelli B., Silva L., Grazian A., Altieri B., Andreani P., Aussel H., Berta S., Bongiovanni A., etc. 2010, , 518, L154
work page 2010
-
[60]
S., Kirchschlager F., De Looze I., Barlow M
Scheffler T., Sartorio N. S., Kirchschlager F., De Looze I., Barlow M. J., Schmidt F. D., 2025, arXiv e-prints, p. arXiv:2512.05046
-
[61]
Schneider R., Maiolino R., 2024, , 32, 2
work page 2024
-
[62]
C., Barnes J., Kartaltepe J., Davies R., van der Werf P., Fomalont E., 2015, , 800, 70
Scoville N., Sheth K., Walter F., Manohar S., Zschaechner L., Yun M., Koda J., Sanders D., Murchikova L., Thompson T., Robertson B., Genzel R., Hernquist L., Tacconi L., Brown R., Narayanan D., Hayward C. C., Barnes J., Kartaltepe J., Davies R., van der Werf P., Fomalont E., 2015, , 800, 70
work page 2015
-
[63]
I., 1959, Similarity and Dimensional Methods in Mechanics
Sedov L. I., 1959, Similarity and Dimensional Methods in Mechanics
work page 1959
-
[64]
Y., Koo B.-C., Hirashita H., 2015, , 807, 100
Seok J. Y., Koo B.-C., Hirashita H., 2015, , 807, 100
work page 2015
-
[65]
Y., Koo B.-C., Onaka T., 2013, , 779, 134
Seok J. Y., Koo B.-C., Onaka T., 2013, , 779, 134
work page 2013
-
[66]
B., Shchekinov Y., 2014, , 443, 3463
Sharma P., Roy A., Nath B. B., Shchekinov Y., 2014, , 443, 3463
work page 2014
- [67]
- [68]
- [69]
- [70]
- [71]
-
[72]
Smith H. E., Lonsdale C. J., Lonsdale C. J., Diamond P. J., 1998, , 493, L17
work page 1998
-
[73]
Smith R. K., Krzewina L. G., Cox D. P., Edgar R. J., Miller W. Warren I., 1996, , 473, 864
work page 1996
-
[74]
Toro E., 2009, Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction. Springer Berlin Heidelberg
work page 2009
- [75]
-
[76]
Valiante R., Schneider R., Salvadori S., Bianchi S., 2011, , 416, 1916
work page 2011
- [77]
- [78]
-
[79]
Vasiliev E. O., Drozdov S. A., Nath B. B., Dettmar R.-J., Shchekinov Y. A., 2023, , 520, 2655
work page 2023
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.