Pith. sign in

REVIEW 3 major objections 6 minor 35 references

Dust destruction and dust cooling are coupled: at 10^7 K and solar metallicity the cooling time halves, not fiftyfolds, and up to ~20% of dust survives supernova shocks when the initial grain sizes are large.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-01 01:03 UTC pith:JPKJOMOK

load-bearing objection A useful grid of dust cooling + sputtering calculations whose headline survival numbers are compromised by an inconsistent isochoric/isobaric density law in Eq. (1). the 3 major comments →

arxiv 2607.25934 v1 pith:JPKJOMOK submitted 2026-07-28 astro-ph.GA astro-ph.SR

The influence of dust destruction on gas cooling

S.A. Drozdov , M. P. Yudkevich This is my paper
classification astro-ph.GA astro-ph.SR
keywords dust destructionthermal sputteringgas coolingsupernova shocksinterstellar dustgrain size distributionearly Universe dustisochoric cooling
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper aims to resolve a puzzle: observations find more dust in the early Universe than models of supernova shocks predict, because those models assumed dust is efficiently destroyed. The authors show that when dust cooling and thermal sputtering are treated together, cooling is fast enough to quench the hot gas before dust is fully destroyed. At solar metallicity and 10^7 K, including dust shortens the cooling time by about a factor of two; neglecting destruction would overestimate this by a factor of ~50. The survival fraction depends strongly on the initial grain-size distribution: a flat distribution (α=1.5) preserves up to ~20% of dust mass, while steeper or evolved distributions preserve less than 10% once the initial temperature exceeds 3×10^6 K. These results suggest that dust destruction is self-limiting: the more dust cools the gas, the less dust is destroyed.

Core claim

The paper's central claim is that dust destruction and gas cooling form a coupled feedback loop: dust grains both cool the hot post-shock gas by absorbing collisions from electrons and protons, and are simultaneously destroyed by thermal sputtering. The authors solve the thermal evolution of a gas parcel with a non-equilibrium cooling function augmented by a dust cooling term, while tracking the time-dependent grain-size distribution as grains shrink and move to smaller size bins. They find that for initial temperatures above 3×10^6 K, almost all dust is destroyed in every model—but the fraction that survives is set by the initial size distribution. A power-law distribution with slope α=1.5,

What carries the argument

The engine of the calculation is the coupled system of the gas energy equation dT_g/dt = -(2 n_g / 3 k_B) Λ(Z, T_g), with n_g = 4 n0 T_g,0/T_g, and the thermal sputtering lifetime τ(a, n_g, T_g) ≈ (2.4×10^-19 µ_i η_i / ρ_i)^-1 (a / n_g) T_g^-1/2 e^(0.54/T_6), which moves grains to smaller size bins. The total cooling function Λ includes a gas term and a dust term Λ_d built from the Dwek (1987) heat-absorption formula, so as grains are sputtered away the cooling rate falls in concert with the temperature. The dust size distribution is discretized into 50 logarithmic bins from 10 Å to 3000 Å, and the time step is chosen from the shortest grain lifetime, so the feedback is captured explicitly.

Load-bearing premise

The quantitative survival fractions rest on the post-shock density law n_g = 4 n0 T_g,0/T_g, which is a constant-pressure relation even though the text calls the process isochoric; if the gas truly stays at constant volume, the density remains 4 n0 and the sputtering and dust-cooling rates evolve differently.

What would settle it

Rerun the same cooling calculation with strictly isochoric gas (n_g = 4 n0 constant) and compare the resulting dust survival fractions and cooling-time ratios to the reported values; if the factor-of-two cooling-time reduction and the survival percentages do not persist, the central claim fails.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • At metallicities below 0.1 Z_sun, dust can be dropped from cooling calculations without changing results, simplifying models of low-metallicity gas.
  • For high initial temperatures (>3×10^6 K) and solar metallicity, dust survival is bounded at ~10–20% in the models, so the 'few percent' destruction assumption underlying early-Universe dust budgets should be revisited.
  • In cloudy media, dust stripped from small clouds can locally increase the dust-to-gas ratio by up to a factor χ~10, accelerating cooling and pushing dust survival from ~1% to ~20% at 10^7 K.
  • The initial grain-size distribution is a controlling parameter: small-grain-rich models cool gas efficiently below 10^6 K but destroy dust faster at high temperature, while large-grain-dominated models preserve dust longer.

Where Pith is reading between the lines

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

  • If the gas expansion law were truly isochoric rather than the constant-pressure relation used in Eq. (1), the sputtering rate and dust cooling would both change, likely shifting the survival contours; the reported numbers are therefore tied to the assumed post-shock density behavior.
  • The feedback loop implies a natural thermostat: in regions where dust is abundant, cooling quenches the hot phase quickly, so dust destruction is self-limiting; this could be tested by comparing cooling times in SNRs of different ages and metallicities.
  • The model ignores grain-grain collisions, which would shatter large grains into small ones; adding shattering would likely increase cooling at low temperatures but also increase destruction, so the 20% survival might be an upper bound in dense environments.
  • Extending the 0-D parcel calculation to a 1-D shock structure with a cloud would let one check whether the dust-tail enhancement survives when the stripped dust is advected and mixed.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript studies post-shock gas cooling coupled to the thermal sputtering and destruction of dust grains. It uses a 0-D model with Eq. (1), a non-equilibrium gas cooling function, the Dwek (1987) dust cooling rate, and four initial grain-size distributions (MRN, α=1.5, WD, SC). Dust grains are evolved in 50 logarithmic size bins down to a 10 Å cutoff. The paper reports cooling-time ratios and surviving-dust-mass fractions on grids of initial gas temperature from 3×10^5 to 3×10^7 K and metallicity from 10^-2 to 2 Z_sun. The main claims are: dust cooling is negligible below 0.1 Z_sun; at 10^7 K and solar metallicity dust cooling with destruction shortens the cooling time by about a factor of two, but ignoring destruction would shorten it by a factor of ~50; less than 10% of dust survives for T_g,0 > 3×10^6 K in most models, with the α=1.5 model doing best; and dust stripped from clouds can accelerate cooling and increase dust survival. The central quantitative results, however, rest on Eq. (1), which is internally inconsistent: the text calls the process isochoric while the adopted density evolution is isobaric, and the prefactor is correct for neither interpretation.

Significance. The paper addresses a timely and important problem: the tension between the large dust masses observed at high redshift and the efficient destruction of dust by SN shocks. The qualitative proposal—that dust cooling shortens the cooling time and thereby increases the surviving dust fraction—is physically plausible and consistent with earlier work. The parameter grid over initial temperature, metallicity, and four dust models is a useful map of the regime, and coupling dust destruction with cooling is a step forward from fixed-abundance calculations. I also credit the authors for clearly separating runs with and without destruction. However, the quantitative claims as they stand are not reproducible from a physically consistent version of the model because of the thermodynamic inconsistency in Eq. (1). The paper is potentially valuable after a careful revision and re-running of the affected results; it is not acceptable in its current form.

major comments (3)
  1. [§2.1, Eq. (1)] The abstract and §2.1 call the model isochoric, but the density law n_g = 4n0 T_g,0/T_g is the constant-pressure post-shock relation; true isochoric cooling would keep n_g = 4n0 fixed. Conversely, if isobaric cooling is intended, the coefficient in Eq. (1) should be 2/(5k_B) for a monatomic gas, not 2/(3k_B). This is not cosmetic: the sputtering rate in Eq. (7) is ∝ n_g T_g^{1/2}, so with the adopted density law it scales as T_g^{-1/2} instead of T_g^{1/2}. Between 10^7 K and 10^6 K the adopted n_g grows by a factor of 10, making the low-temperature sputtering rate about three times larger than in a fixed-density run. The surviving-mass fractions (Fig. 5), cooling-time ratios (Fig. 4), and the factor-of-50 statement in Fig. 1 are all outputs of that inconsistent trajectory. The paper should be re-run under a single thermodynamic assumption, or a sensitivity test separating the two choice
  2. [§3.2 and Conclusion] The result that dust cooling is negligible for Z/Z_sun < 0.1 is built into the model through the assumption ζ = (1/120) Z/Z_sun in Eq. (3). Because the dust number density is scaled linearly with metallicity, the dust cooling term is forced to disappear at low Z. This conclusion is therefore an assumption-driven expectation rather than an independent finding. The authors should state this explicitly and, if the claim is intended as a physical conclusion, justify the linear scaling of the dust-to-gas ratio with Z in the environments considered.
  3. [§2.2 and Fig. 5] The numerical discretization of the size distribution is described (50 logarithmic bins, 10 Å cutoff, automatic time step from the minimum grain lifetime), but no convergence test is reported. The claimed survival fractions, especially the 10% versus 20% differences between models, are sensitive to the advection of grains between size bins and to the removal of grains below a_min. A convergence study varying bin number and time-step tolerance should be included to establish that the survival percentages are numerical, not just physical, results.
minor comments (6)
  1. [§2.1, Eq. (3)] The dust number density n_d(ζ) is not explicitly defined; please specify how it is computed from the size distribution and the dust-to-gas mass ratio.
  2. [Fig. 1] The axis label t[yr]/n0 is ambiguous. Clarify whether the cooling time is multiplied or divided by n0 and state the units of n0 explicitly.
  3. [Eq. (7)] The dimensions of the right-hand side are not transparent. Please make the units of µ_i, ρ_i, η_i, and the numerical prefactor explicit in one place.
  4. [§2.1] The assumption n_e = n_g is applied over the full integration down to 10^4 K, but the gas is not fully ionized at low temperatures. This may overestimate electron-driven dust cooling at late times; please justify or restrict the range.
  5. [§4] The sentence 'reduces the ratio ... by almost an order' should read 'by almost an order of magnitude'.
  6. [References] Several references are dated 2025–2026; please provide arXiv IDs or DOIs where available so the cooling functions and dust models can be checked.

Circularity Check

0 steps flagged

No significant circularity; central results are numerical outputs of an externally parameterized model. The only self-citation (Nath et al. 2023) is not load-bearing.

full rationale

The paper's derivation chain is a forward model: Eq. (1) integrates dT_g/dt from tabulated gas cooling functions (Vasiliev 2013) plus a dust cooling term (Dwek 1987), and dust destruction follows an adopted sputtering rate (Polikarpova & Shchekinov 2017) and initial size distributions (MRN, alpha=1.5, WD, SC). None of the central claims—cooling-time ratios, survival fractions, the factor-of-two/50 comparisons—are obtained by fitting those outputs back into the inputs; they are conditional numerical consequences of stated assumptions. The low-metallicity conclusion is not purely tautological: although zeta=(1/120)Z/Zsun makes the dust cooling amplitude proportional to Z, inclusion of destruction feedback means the low-Z behavior also depends on the longer cooling time, so it is a model result rather than a restatement of the scaling. The only self-citation is Nath et al. (2023), used to motivate the alpha=1.5 dust model; the same choice is independently supported by Maiolino et al. (2004), Nozawa et al. (2007), and Nishida et al. (2022), so it is not load-bearing. The isochoric/isobaric inconsistency in Eq. (1) is a physical/correctness concern, not a circularity: it does not reduce any predicted quantity to an assumed quantity by construction.

Axiom & Free-Parameter Ledger

4 free parameters · 6 axioms · 0 invented entities

The model relies on externally calibrated cooling functions, sputtering yields, and initial size distributions; the only parameters the paper itself sets are the normalization of dust abundance and a 10 Å destruction cutoff. No new physical entities are introduced.

free parameters (4)
  • Solar dust-to-gas mass ratio ζ_sun = 1/120 (scaled as Z/Zsun)
    Sets the normalization of dust cooling; because ζ is taken proportional to Z, the claim that dust cooling is negligible below 0.1 Z_sun is a direct consequence of this assumed scaling (Section 2.1).
  • Grain destruction cutoff a_min = 10 Å
    Grains reaching 10 Å are removed from the calculation; this cutoff impacts small-grain cooling in MRN/WD models and the reported survival fractions (Section 2.1).
  • Sputtering yield normalizations η_C, η_Si = 2.5, 7.5
    Taken from Polikarpova & Shchekinov (2017); survival times and hence survival fractions scale inversely with these numbers (Eq. 7).
  • Initial size distribution parameters = MRN α=3.5; flat α=1.5; WD b_C=6e-5, R_V=3.1; SC 15 Msun, 2000 days
    The central comparison is among these externally prescribed distributions; results are conditional on them.
axioms (6)
  • domain assumption The gas cools at constant volume (isochoric)
    Stated in abstract and Section 2; however Eq. (1) sets n_g=4n0Tg0/Tg, which is isobaric. If this axiom is taken literally it is violated by the paper's own equation.
  • domain assumption Fully ionized plasma with n_e=n_g at T>10^4 K
    Section 2.2; underpins both dust cooling and sputtering rates.
  • domain assumption Dust cooling follows Dwek (1987) collisional heat transfer with electron-dominated absorption
    Section 2.1, Eqs. (3)-(5).
  • domain assumption Thermal sputtering is the only destruction channel; kinetic sputtering (~15%) and shattering are neglected
    Section 2.2 and Discussion.
  • domain assumption Non-equilibrium gas cooling tables from Vasiliev (2013) accurately represent Λ_g over the grid
    Section 2.1.
  • domain assumption Dust-to-gas mass ratio scales as ζ=(1/120)Z/Z_sun
    Section 2.1; makes the Z<0.1 threshold essentially an assumption.

pith-pipeline@v1.3.0-alltime-deepseek · 10814 in / 16964 out tokens · 155199 ms · 2026-08-01T01:03:28.841103+00:00 · methodology

0 comments
read the original abstract

The observed dust abundance in the early Universe significantly exceeds the predictions of models assuming its efficient destruction at supernova shock wave fronts. We investigate the effect of dust on gas cooling behind the shock wave front, taking into account both dust cooling and thermal sputtering of dust grains for various interstellar dust models. Isochoric cooling of a gas element is considered for various initial temperatures (from $3\cdot 10^{5}$ K to $3\cdot 10^{7}$ K) and metallicities (from $10^{-2}$ Z$_{\odot}$ to 2 Z$_{\odot}$). It is shown that at metallicities below 0.1 Z$_{\odot}$, the effect of dust on gas cooling is negligible. The abundance of small grains in some models significantly accelerates cooling at temperatures below $10^{6}$ K, whereas large grains dominate cooling at high temperatures ($> 3\cdot 10^{6}$ K). For an initial temperature $T_{g,0} > 3\cdot 10^{6}$ K, less than 10% of the dust mass survives for some models of the initial size distribution. The maximum survival rate (up to 20% at $T_{g,0} = 3\cdot 10^{6}$ K and solar metallicity) is achieved for the model with the shallowest grain size distribution. The role of gas inhomogeneities is discussed: dust stripping from clouds by a shock wave can both decrease the cloud lifetime and contribute to the creation of dust tails, where the cooling of hot gas is accelerated, thereby increasing the dust survival rate.

Figures

Figures reproduced from arXiv: 2607.25934 by M. P. Yudkevich, S.A. Drozdov.

Figure 2
Figure 2. Figure 2: At such values of the initial gas temper [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

35 extracted references · 5 canonical work pages

  1. [1]

    , keywords =

    The Infrared Diagnostic of a Dusty Plasma with Applications to Supernova Remnants. , keywords =. doi:10.1086/165774 , adsurl =

    work page doi:10.1086/165774

  2. [2]

    , keywords =

    Dusty clump survival in supernova ejecta: Dust-mediated growth versus crushing by the reverse shock. , keywords =. doi:10.1051/0004-6361/202556389 , archivePrefix =. 2509.08887 , primaryClass =

    arXiv

  3. [3]

    , keywords =

    Dust Grain-Size Distributions and Extinction in the Milky Way, Large Magellanic Cloud, and Small Magellanic Cloud. , keywords =. doi:10.1086/318651 , archivePrefix =. astro-ph/0008146 , primaryClass =

    Pith/arXiv arXiv

  4. [4]

    , keywords =

    The size distribution of interstellar grains. , keywords =. doi:10.1086/155591 , adsurl =

    doi:10.1086/155591

  5. [5]

    Bulletin of the Lebedev Physics Institute , keywords =

    Dynamics of Gas and Dust during Interaction of Diffuse Clouds with a Shock Wave. Bulletin of the Lebedev Physics Institute , keywords =. doi:10.3103/S1068335624601560 , adsurl =

    work page doi:10.3103/s1068335624601560

  6. [6]

    New Astronomy , keywords =

    Inhibited destruction of dust by supernova in a clumpy medium. New Astronomy , keywords =. doi:10.1016/j.newast.2024.102293 , archivePrefix =. 2404.18317 , primaryClass =

    arXiv 2024

  7. [7]

    , keywords =

    Dust evolution in a supernova interacting with the ISM. , keywords =. doi:10.1093/mnras/stad3820 , archivePrefix =. 2308.03106 , primaryClass =

    Pith/arXiv arXiv

  8. [8]

    , keywords =

    Dust destruction by the supernova remnant forward shock in a turbulent interstellar medium. , keywords =. doi:10.1051/0004-6361/202556499 , archivePrefix =. 2512.05046 , primaryClass =

    arXiv

  9. [9]

    , keywords =

    Time-dependent Cooling and Grain Destruction in Hot Dusty Plasmas: A Simplified Model and Principal Results. , keywords =. doi:10.1086/178198 , adsurl =

    doi:10.1086/178198

  10. [10]

    , keywords =

    Non-equilibrium cooling rate for a collisionally cooled metal-enriched gas. , keywords =. doi:10.1093/mnras/stt189 , archivePrefix =. 1302.0159 , primaryClass =

    Pith/arXiv arXiv

  11. [11]

    Astronomy Reports , year = 2017, month = feb, volume =

    Dust in galaxy clusters. Astronomy Reports , year = 2017, month = feb, volume =. doi:10.1134/S1063772917020044 , adsurl =

    work page doi:10.1134/s1063772917020044 2017

  12. [12]

    Physics of SW and high-temp. hydr. phen

  13. [13]

    , keywords =

    Condensation of dust in the ejecta of Type II-P supernovae. , keywords =. doi:10.1051/0004-6361/201424969 , archivePrefix =. 1412.5522 , primaryClass =

    Pith/arXiv arXiv

  14. [14]

    Nature , keywords =

    A supernova origin for dust in a high-redshift quasar. Nature , keywords =. doi:10.1038/nature02930 , archivePrefix =. astro-ph/0409577 , primaryClass =

    Pith/arXiv arXiv

  15. [15]

    , keywords =

    Evolution of Dust in Primordial Supernova Remnants: Can Dust Grains Formed in the Ejecta Survive and Be Injected into the Early Interstellar Medium?. , keywords =. doi:10.1086/520621 , archivePrefix =. 0706.0383 , primaryClass =

    Pith/arXiv arXiv

  16. [16]

    Interstellar Dust , year = 1989, editor =

    Dust Destruction in the Interstellar Medium. Interstellar Dust , year = 1989, editor =

    1989

  17. [17]

    , keywords =

    Grain Destruction in Shocks in the Interstellar Medium. , keywords =. doi:10.1086/174689 , adsurl =

    doi:10.1086/174689

  18. [18]

    , keywords =

    Destruction of Interstellar Dust in Evolving Supernova Remnant Shock Waves. , keywords =. doi:10.1088/0004-637X/803/1/7 , archivePrefix =. 1502.00929 , primaryClass =

    Pith/arXiv arXiv

  19. [19]

    , keywords =

    Destruction mechanisms for interstellar dust. , keywords =. doi:10.1086/157206 , adsurl =

    doi:10.1086/157206

  20. [20]

    , keywords =

    Dust grains from the heart of supernovae. , keywords =. doi:10.1051/0004-6361/201527432 , archivePrefix =. 1601.06770 , primaryClass =

    Pith/arXiv arXiv

  21. [21]

    Dust sources

    The formation and cosmic evolution of dust in the early Universe: I. Dust sources. Astron. Astrophys. Rev. , keywords =. doi:10.1007/s00159-024-00151-2 , archivePrefix =. 2310.00053 , primaryClass =

    Pith/arXiv arXiv

  22. [22]

    , keywords =

    Two Massive, Compact, and Dust-obscured Candidate z ≃ 8 Galaxies Discovered by JWST. , keywords =. doi:10.3847/1538-4357/acef21 , archivePrefix =. 2304.12347 , primaryClass =

    Pith/arXiv arXiv

  23. [23]

    Nature , keywords =

    Spectroscopic confirmation of two luminous galaxies at a redshift of 14. Nature , keywords =. doi:10.1038/s41586-024-07860-9 , archivePrefix =. 2405.18485 , primaryClass =

    Pith/arXiv arXiv

  24. [24]

    Galaxies , keywords =

    Dust at the Cosmic Dawn. Galaxies , keywords =. doi:10.3390/galaxies13030064 , adsurl =

    work page doi:10.3390/galaxies13030064

  25. [25]

    Physics of the ISM and IGM

  26. [26]

    , keywords =

    Survival of dust in super-dusty galaxies at redshifts z ≍ 5-8. , keywords =. doi:10.1088/1475-7516/2025/11/030 , archivePrefix =. 2506.05591 , primaryClass =

    Pith/arXiv arXiv 2025

  27. [27]

    , keywords =

    Temperature Fluctuations and Infrared Emission from Dust Particles in a Hot Gas. , keywords =. doi:10.1086/163995 , adsurl =

    work page doi:10.1086/163995

  28. [28]

    , keywords =

    Infrared Observational Manifestations of Young Dusty Super Star Clusters. , keywords =. doi:10.3847/0004-637X/816/1/39 , archivePrefix =. 1511.03382 , primaryClass =

    Pith/arXiv arXiv

  29. [29]

    , keywords =

    The Physics of Grain-Grain Collisions and Gas-Grain Sputtering in Interstellar Shocks. , keywords =. doi:10.1086/174488 , adsurl =

    doi:10.1086/174488

  30. [30]

    , keywords =

    Dust-free starburst galaxies at redshifts z > 10. , keywords =. doi:10.1093/mnras/stad505 , archivePrefix =. 2211.12378 , primaryClass =

    Pith/arXiv arXiv

  31. [31]

    , keywords =

    A new galaxy spectral energy distribution model consistent with the evolution of dust. , keywords =. doi:10.1093/mnras/stac1355 , archivePrefix =. 2205.07591 , primaryClass =

    Pith/arXiv arXiv

  32. [32]

    , keywords =

    Soft X-ray spectrum of a hot plasma. , keywords =. doi:10.1086/190486 , adsurl =

    doi:10.1086/190486

  33. [33]

    , keywords =

    Cooling Functions for Low-Density Astrophysical Plasmas. , keywords =. doi:10.1086/191823 , adsurl =

    doi:10.1086/191823

  34. [34]

    , keywords =

    Ionization balance for optically thin plasmas: Rate coefficients for all atoms and ions of the elements H to NI. , keywords =. doi:10.1051/aas:1998330 , archivePrefix =. astro-ph/9806391 , primaryClass =

    Pith/arXiv arXiv

  35. [35]

    , keywords =

    The dust-scattering component of X-ray extinction: effects on continuum fitting and high-resolution absorption edge structure. , keywords =. doi:10.1093/mnras/stw376 , archivePrefix =. 1602.01100 , primaryClass =

    Pith/arXiv arXiv