Dust enrichment from core-collapse supernovae and extinction curves in the high-redshift universe

Joris Witstok; Koki Otaki; Luca Graziani; Marco Limongi; Raffaella Schneider; Roberto Maiolino; Simone Bianchi

arxiv: 2604.26750 · v1 · submitted 2026-04-29 · 🌌 astro-ph.GA

Dust enrichment from core-collapse supernovae and extinction curves in the high-redshift universe

Koki Otaki , Raffaella Schneider , Simone Bianchi , Joris Witstok , Luca Graziani , Marco Limongi , Roberto Maiolino This is my paper

Pith reviewed 2026-05-07 12:12 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords dustcurvesextinctionmathrmshockreverseeffectenrichment

0 comments

The pith

Reverse shocks destroy small grains from supernovae, leaving mostly amorphous carbon and yielding flatter UV extinction curves consistent with high-redshift galaxy observations.

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

Core-collapse supernovae create dust grains of many sizes in the first galaxies. The model follows how these grains grow or shrink over the first 30 million years in a star-forming burst, accounting for the reverse shock that sweeps back through the exploded material. This shock destroys nearly all grains smaller than 10 nanometers, leaving larger amorphous carbon grains as the main survivor. The resulting dust mixture absorbs and scatters ultraviolet light differently, producing a flatter attenuation curve with a wide bump near 2500 angstroms instead of the sharper features seen without shocks. The total dust mass relative to stars falls to about one ten-thousandth due to this destruction. These modeled curves and dust amounts line up with what JWST measures in distant galaxies.

Core claim

We find that our models are consistent with the observed attenuation curve and emissivity of high-redshift galaxies and show that the reverse shock processing significantly affects dust enrichment and grain properties such as extinction curves and emissivity in supernova yields for high-redshift galaxies.

Load-bearing premise

The assumption that a single starburst system with a standard initial mass function and varied but fixed progenitor parameters (rotation, metallicity, ISM density) accurately captures the conditions and dust processing in real high-redshift galaxies.

read the original abstract

Recent JWST observations have revealed that some galaxies at $z \gtrsim 7$ generally exhibit relatively flat ultraviolet (UV) attenuation curves and a weak UV bump. These features suggest that the first dust grains formed rapidly, possibly originating from core-collapse supernovae (SNe). We investigate the time evolution of grain size distributions and extinction curves in the early phase of dust enrichment for different parameters of progenitor stars, rotation velocities, metallicity, and interstellar medium densities, including the effect of the reverse shock. We model a single starburst system assuming an initial mass function. Extinction curves are calculated from the grain size distribution for each dust species. The total dust-to-stellar mass ratio at $30 \,\mathrm{Myr}$ is $M_\mathrm{dust}/M_\star \sim 10^{-3}$ before the passage of the reverse shock, but we find it to be at most $M_\mathrm{dust}/M_\star \sim 10^{-5}$ due to the destruction effect of the reverse shock. This effect destroys grains smaller than $\sim 10\,\mathrm{nm}$ and makes amorphous carbon the dominant species, resulting in a flatter extinction curve with a wide bump at $2500\,\mathrm{\mathring{A}}$ compared to the no-reverse shock models. We find that our models are consistent with the observed attenuation curve and emissivity of high-redshift galaxies and show that the reverse shock processing significantly affects dust enrichment and grain properties such as extinction curves and emissivity in supernova yields for high-redshift galaxies.

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

Reverse shocks in core-collapse SN models can flatten high-z extinction curves to match JWST data, but the single-burst assumption limits how far the result generalizes.

read the letter

The main point is that including reverse-shock destruction of small grains in supernova dust yields produces flatter UV attenuation curves with a broad 2500 Å bump, plus a dust-to-stellar mass ratio down to 10^{-5}, that line up with JWST observations at z greater than or equal to 7. They reach this by evolving grain size distributions from core-collapse SNe in a single starburst, varying progenitor rotation, metallicity, and ISM density, then computing extinction curves for each species. Amorphous carbon ends up dominant after small grains are destroyed, which changes both the curve shape and emissivity in the direction of the data. This is a straightforward extension of existing SN dust models rather than a new framework, but the explicit reverse-shock step on grain sizes for high-redshift conditions is the concrete addition. The parameter sweeps show the effect is robust within the explored range. The soft spot is the fixed single instantaneous starburst at the 30 Myr snapshot. Real z greater than or equal to 7 galaxies have longer or bursty star formation histories with overlapping supernova generations, which would shift the timing of dust injection versus destruction and change the final grain-size distribution. They vary progenitor and ISM parameters but keep the star formation history fixed, so the reported consistency is not yet shown to survive more extended histories. The abstract claims a match to observed attenuation and emissivity, and the physical setup looks internally consistent, though quantitative fit details and exact evolution equations would need checking in the full text. This work is aimed at people modeling early dust enrichment or interpreting JWST attenuation curves. A reader focused on supernova contributions to high-redshift dust would get direct value from the grain-size results and parameter study. It deserves a serious referee because it tackles a timely observational puzzle with explicit calculations rather than hand-waving, even if revisions for broader star formation histories are likely needed. I would send it to peer review.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard astrophysical assumptions about supernova dust formation and shock processing, with several free parameters varied to explore outcomes.

free parameters (2)

progenitor star parameters
Rotation velocities, metallicities, and masses of progenitor stars are varied across models.
ISM density
Interstellar medium densities are included as a variable parameter affecting grain evolution.

axioms (2)

domain assumption Standard initial mass function governs the distribution of progenitor stars in the starburst system.
Invoked to model the collective supernova contribution in a single starburst.
domain assumption Reverse shock destroys grains smaller than approximately 10 nm.
Used to explain the shift to amorphous carbon dominance and flatter curves.

pith-pipeline@v0.9.0 · 5606 in / 1397 out tokens · 51895 ms · 2026-05-07T12:12:01.067111+00:00 · methodology

discussion (0)

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

The Lifecycle and Emission Properties of PAHs in Cosmological Hydrodynamic Galaxy Formation Simulations
astro-ph.GA 2026-06 unverdicted novelty 7.0

Cosmological zoom-in simulations find that grain-grain shattering in diffuse ISM gas drives rising PAH mass fraction with time, naturally producing the observed PAH-metallicity relation and inverse qPAH-molecular gas trends.
Grain-size evolution and rapid dust growth in high-redshift galaxies
astro-ph.GA 2026-06 conditional novelty 5.0

A multiphase ISM grain-size model with low supernova dust yield reproduces observed dust-to-stellar mass ratios and UV luminosity functions at z=7-12 by letting small grains seed rapid metal accretion.