pith. machine review for the scientific record. sign in

arxiv: 2604.24846 · v1 · submitted 2026-04-27 · 🌌 astro-ph.GA

Recognition: unknown

SN Ia Population Machine. I. A Unified Cosmological Simulation-Binary Synthesis Framework Establishing Non-universal Delay-time Distributions and Cosmic Progenitor-channel Dominance Crossover

Suk-Jin Yoon , Inhyuk Park , Woong-Bae G. Zee , Chul Chung , Jun-Sung Moon , Sanjaya Paudel , Kiyun Yun , Myung-Hun Kim
show 1 more author
Eun-Taek Gim
Authors on Pith no claims yet

Pith reviewed 2026-05-08 02:17 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords Type Ia supernovaedelay-time distributionsprogenitor channelsbinary population synthesiscosmological simulationsredshift evolutionsingle degeneratedouble degenerate
0
0 comments X

The pith

Type Ia supernova delay-time distributions depend on progenitor channel and metallicity, varying across galaxies and with redshift and shifting from single-degenerate to double-degenerate dominance near z=0.5.

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

The paper constructs a forward-modeling pipeline that links large cosmological hydrodynamic simulations to binary population synthesis codes in order to generate populations of Type Ia supernovae complete with explosion times and channel labels. It demonstrates that the delay between star formation and explosion is not a fixed distribution but changes systematically with the dominant formation channel and the metallicity of the host environment. A sympathetic reader cares because these supernovae serve as the primary distance indicators for mapping cosmic acceleration, so any unaccounted variation in their properties with look-back time could introduce systematic biases in the inferred expansion history. The model further predicts a cosmic transition in which the single-degenerate channel gives way to the double-degenerate channel roughly five billion years ago.

Core claim

By treating star particles from IllustrisTNG as simple stellar populations and evolving their binaries with COMPAS, the framework produces synthetic supernova catalogs that match observed host demographics, rate trends, and a progenitor-age step. The central result is that delay-time distributions are intrinsically non-universal, their shape set by channel and metallicity, while the dominant progenitor population undergoes a demographic crossover from single-degenerate to double-degenerate near redshift 0.5.

What carries the argument

The cosmology-BPS pipeline that couples IllustrisTNG star particles to COMPAS binary evolution to tag each synthetic supernova with an explosion time and a progenitor channel.

If this is right

  • Delay-time distributions vary systematically across different host galaxies and change with redshift.
  • The single-degenerate channel dominates at earlier cosmic times while the double-degenerate channel becomes dominant after the crossover near z=0.5.
  • A single globally calibrated standardization for SN Ia luminosities is undermined by the evolving channel mixture.
  • Evolution in both the delay-time distribution and the channel mix can produce redshift-dependent systematics in measured supernova luminosities.

Where Pith is reading between the lines

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

  • High-redshift supernova samples used for cosmology may carry a different average luminosity bias than local samples because of the changing channel mix.
  • Separate standardization corrections could be required for supernovae in metal-poor versus metal-rich environments at the same redshift.
  • The framework opens the possibility of predicting how host-galaxy properties correlate with Hubble residuals at different look-back times.

Load-bearing premise

That IllustrisTNG star particles can be treated as simple stellar populations whose binary evolution is fully captured by the COMPAS code with its chosen parameters for common-envelope efficiency, mass transfer, and metallicity scaling.

What would settle it

A direct measurement showing that delay-time distributions extracted from supernova samples in low-metallicity versus high-metallicity hosts at fixed redshift are statistically identical would falsify the claim of intrinsic non-universality.

Figures

Figures reproduced from arXiv: 2604.24846 by Chul Chung, Eun-Taek Gim, Inhyuk Park, Jun-Sung Moon, Kiyun Yun, Myung-Hun Kim, Sanjaya Paudel, Suk-Jin Yoon, Woong-Bae G. Zee.

Figure 1
Figure 1. Figure 1: A conceptual overview of our SN Ia population-synthesis framework, which unified cosmological hydrodynamic galaxy simulations with binary population synthesis to generate SN Ia populations from individual galaxies to galaxy clusters/groups to cosmic volumes. (a) A slab of the TNG50-1 snapshot at 𝑧 = 0. Subhalos with at least one star particle are shown as white dots with an opacity of 30 % and a constant s… view at source ↗
Figure 2
Figure 2. Figure 2: shows the resulting IMFs for all (single and binary) stars and for binary components. These IMFs are generated assuming a star particle with 4 × 106 M⊙, which we adopt as a conservative upper-envelope mass that brackets the maximum star-particle mass expected in TNG100 (typically ≲ 2.8×106 M⊙). In other words, 4× 106 M⊙ is a sufficiently large reference mass, ensuring that our IMF sampling covers the full … view at source ↗
Figure 3
Figure 3. Figure 3: Three 2D projections (panels 𝑎–𝑐) of our adopted 3D 𝑀WD 1 –𝑀2–log 𝑃orb parameter-space polyhedra (panel 𝑑) for the single-degenerate SN Ia channels (WD+MS, WD+GB, and WD+He). The blue, cyan, and green regions denote the WD+MS (Meng & Podsiadlowski 2017), WD+GB (Wang et al. 2010), and WD+He (Wang & Han 2010a) channels, respectively. Contours with different line types indicate the regions associated with the… view at source ↗
Figure 4
Figure 4. Figure 4: The star-particle-level delay-time distributions: the SN Ia rate as a function of delay time (𝜏 in units of Gyr) following a single star-burst of 1 M⊙ simple-stellar-population formation. (a) Star-particle-level SSP DTDs for Z = 0.0142 (the adopted solar metallicity in COMPAS). Blue, cyan, and green solid curves show the SD near-𝑀Ch channels (WD+MS, WD+GB, and WD+He), and the red solid curve shows the DD n… view at source ↗
Figure 5
Figure 5. Figure 5 view at source ↗
Figure 6
Figure 6. Figure 6: Cosmic volume-normalized SN Ia rate (𝑅vol) as a function of red￾shift for different observational time windows (Δ𝑡). To quantify stochasticity, we consider four observational time windows, Δ𝑡 = 1, 3, 10, and 105 yr (top to bottom), and partition 0 ≤ 𝑧 ≤ 5 into 50 bins. For each Δ𝑡, we generate 100 independent realizations, plotted as colored points connected by lines. In each redshift bin, we compute the m… view at source ↗
Figure 7
Figure 7. Figure 7: Host-galaxy fraction and per-host SN Ia counts as functions of the observational time window, Δ𝑡, for all galaxies with 𝑀∗ ≥ 108.0 M⊙ over 0 ≤ 𝑧 ≤ 5. (a) Host fraction, 𝑁host(Δ𝑡)/𝑁All [%], versus Δ𝑡. For Δ𝑡 = 105 yr, essentially every galaxy (99.74 %) hosts at least one SN Ia, whereas for Δ𝑡 = 3 yr—comparable to the baselines of modern time-domain surveys—the fraction drops to ∼ 0.8 %. (b) Mean number of S… view at source ↗
Figure 8
Figure 8. Figure 8: Representative SN Ia host galaxies from TNG50-1 and their SFHs, together with their SN Ia formation-time distributions observed in a mock survey with an observational time window of Δ𝑡 = 105 yr. (Left column) Host-galaxy stellar maps with key properties in the legend. White dots indicate star particles, shown with mass-weighted brightness, while cyan and magenta dots mark star particles that host SNe Ia pr… view at source ↗
Figure 9
Figure 9. Figure 9: Same as view at source ↗
Figure 10
Figure 10. Figure 10: Basic demographics of all galaxies and SN Ia host galaxies in TNG100 over 0 ≤ 𝑧 ≤ 3 (the global sample). (a–c) Redshift evolution of the star-particle-mass-weighted mean stellar age (𝑇∗), total stellar mass (𝑀∗; the sum of all constituent star-particle masses), and sSFR for the full galaxy population. Colored 2D maps show the logarithmic galaxy number density (dark-blue-to-yellow ‘cividis’ color scale), w… view at source ↗
Figure 11
Figure 11. Figure 11: (a–f) Local analogue of view at source ↗
Figure 12
Figure 12. Figure 12: Analogous to view at source ↗
Figure 13
Figure 13. Figure 13: Mass-normalized SN Ia rate (𝑅mass) as a function of six key host properties. (a–f) 𝑅mass versus (𝑎) redshift, (𝑏) mass-weighted mean stellar age, (𝑐) mass-weighted mean stellar metallicity, (𝑑) stellar mass, (𝑒) SFR, and ( 𝑓 ) specific SFR over the global sample (0 ≤ 𝑧 ≤ 3). Colored 2D maps encode the logarithmic galaxy number density (color bar), and white contours mark the 0.3, 1, 2, and 3 𝜎 levels. Sol… view at source ↗
Figure 14
Figure 14. Figure 14: Same as view at source ↗
Figure 15
Figure 15. Figure 15: Recovery of the DTDs for SNe Ia: (𝑎) the local sample (0 ≤ 𝑧 ≤ 0.1); (𝑏– 𝑓 ) higher-redshift slices centered at 𝑧 ≃ 0.5, 1.0, 1.5, 2.0, and 2.5. Thin black lines show the SFH of all galaxies, while thin blue, red, and purple lines show the corresponding SN Ia rates from the SD and DD channels and their sum, respectively. Deconvolution of the SN Ia rate with the SFH yields the recovered DTD, shown by the t… view at source ↗
Figure 16
Figure 16. Figure 16: Recovered SN Ia DTD for 0 ≤ 𝑧 ≤ 3 (the global sample). The main panel uses a linear 𝑥-axis and logarithmic 𝑦-axis to emphasize the DTD morphology; the insets show the data in linear–linear scaling (left; zoomed at short delays) and in log–log scaling (right). Blue, cyan, and green solid curves show the SD near-𝑀Ch channels (WD+MS, WD+GB, and WD+He), and red solid curve shows the DD near-𝑀Ch channel (WD+WD… view at source ↗
Figure 17
Figure 17. Figure 17: Metallicity-dependent SN Ia DTDs in the local sample (0 ≤ 𝑧 ≤ 0.1). (Top row) DTDs binned by star-particle (progenitor) metallicity, 𝑍pro,𝑖, into 10 equal-number intervals, of which four representative bins are shown. Blue and red points show the SD and DD DTDs, respectively, while black points denote their pointwise sum. Solid lines indicate the corresponding power-law fits for the SD, DD, and total chan… view at source ↗
Figure 18
Figure 18. Figure 18: Analogous to view at source ↗
Figure 19
Figure 19. Figure 19: Host mean SN Ia progenitor age (𝑇pro) as a function of host-galaxy properties for two redshift slices and two observational time-window choices. (a–f) Two redshift slices—𝑧 ≃ 0 (0.0 ≤ 𝑧 ≤ 0.1; 𝑎–𝑐) and 𝑧 ≃ 0.5 (0.45 ≤ 𝑧 ≤ 0.55; 𝑑– 𝑓 )—for Δ𝑡 = 105 yr. Panels (𝑎 & 𝑑) show 𝑇pro versus host age (𝑇∗), with the dotted line indicating equality, and panels (𝑏 & 𝑒) and (𝑐 & 𝑓 ) show 𝑇pro versus host stellar mass … view at source ↗
Figure 20
Figure 20. Figure 20: Progenitor ages of individual SNe Ia as a function of host-galaxy properties: (𝑎) 0 ≤ 𝑧 ≤ 3 (global sample), (𝑏) 0 ≤ 𝑧 ≤ 0.1 (local sample), and (𝑐) 0.45 ≤ 𝑧 ≤ 0.55 (the 𝑧 ≃ 0.5 sample). From left to right, panels show progenitor age versus host redshift, host mass-weighted mean stellar age, host stellar mass, and host sSFR. Colored 2D maps show the logarithmic number density of SNe Ia (dark-blue-to-teal-… view at source ↗
Figure 21
Figure 21. Figure 21: The SD (blue) and DD (red) event counts and their fractional contributions as functions of host-galaxy properties in two redshift slices: 0 ≤ 𝑧 ≤ 0.1 (upper two rows) and 0.45 ≤ 𝑧 ≤ 0.55 (lower two rows). (Upper two rows) For the local sample, each column shows host age, host mass, and host sSFR, matching the 2nd–4th panels of view at source ↗
Figure 22
Figure 22. Figure 22: Cosmic evolution of the host-mass-driven and host-sSFR-driven “steps” in progenitor ages. (Left column) Host-mass-driven progenitor-age step in the progenitor age versus host mass plane. Top-row six panels show colored 2D maps of the logarithmic SNe Ia number density (color bar) in six redshift bins, 𝑧 ≃ 0.0, 0.5, 1.0, 1.5, 2.0, and 2.5; solid and dashed black curves trace the mean and median relations, r… view at source ↗
Figure 23
Figure 23. Figure 23: Cosmic evolution of the volume-normalized SN Ia rate (𝑅vol) over 0 ≤ 𝑧 ≤ 3. Black points show our model predictions for Δ𝑡 = 105 yr as a function of redshift; black solid curve represents a sixth-order polynomial fit, adopted for consistency with view at source ↗
Figure 24
Figure 24. Figure 24: Analogous to view at source ↗
read the original abstract

We present a forward-modeling framework for synthesizing Type Ia supernova (SN Ia) populations by coupling cosmological hydrodynamic simulations to binary population synthesis (BPS). Using IllustrisTNG star particles as simple stellar populations, we generate binaries and evolve them with COMPAS to produce synthetic SNe Ia tagged with explosion times and progenitor channels (single- and double-degenerate; SD and DD). This cosmology-BPS pipeline enables self-consistent, end-to-end tracking of SN Ia populations from individual galaxies to cosmic scales. The model reproduces key SN-related observables, including host-galaxy demographics, delay-time distributions (DTDs), SN-rate trends with host properties and redshift, and a progenitor-age 'step' implicated by the mass step in Hubble residuals. Our main findings are as follows. (1) Contrary to the standard assumption, DTDs appear intrinsically non-universal: their form depends on progenitor channel and metallicity, and thus varies systematically across hosts and with redshift. The commonly adopted DTD is therefore best regarded as a population-averaged approximation rather than a fundamental kernel. (2) We predict that the dominant SN Ia progenitor population shifts from SD to DD with cosmic time, with a demographic crossover near z = 0.5 (~5.2 Gyr ago). This non-monolithic SN Ia population with a redshift-dependent SD/DD mixture weakens the universality implicit in a single globally calibrated standardization. Taken together, evolution in both the DTD and the channel mixture can imprint redshift-dependent systematics on SN Ia luminosities, strengthening the case for jointly inferring progenitor/host-driven effects alongside cosmic acceleration. The full catalogue and analysis scripts are available via Zenodo.

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.

Referee Report

3 major / 2 minor

Summary. The paper introduces a forward-modeling pipeline that couples IllustrisTNG star particles (treated as simple stellar populations) to the COMPAS binary population synthesis code to generate synthetic SN Ia populations with tagged explosion times and progenitor channels (SD and DD). It claims that the resulting delay-time distributions are intrinsically non-universal, depending on channel and metallicity, and that the dominant channel undergoes a demographic crossover from SD to DD near z ≈ 0.5. The framework is reported to reproduce several global observables including host demographics, DTD shapes, rate trends, and a progenitor-age step.

Significance. If the robustness of the channel crossover and non-universality claims can be established, the work would be significant for SN Ia cosmology: it provides a concrete mechanism by which redshift-dependent progenitor mixtures could imprint systematics on luminosity standardization, and it supplies a publicly released catalog that enables direct tests against future surveys.

major comments (3)
  1. [Methods and Results sections describing COMPAS setup and DTD/channel fractions] The central claims of channel-dependent DTD shapes and an SD-to-DD crossover at z ≈ 0.5 rest on a single COMPAS realization with fixed values for common-envelope efficiency, mass-transfer efficiency, and metallicity scaling. No systematic variation of these parameters (known to shift channel fractions by factors of several in the BPS literature) is presented, so it is unclear whether the reported crossover redshift is stable or an artifact of the chosen parameter set.
  2. [Section on coupling to IllustrisTNG star particles] Treating each IllustrisTNG star particle as a perfectly coeval, single-metallicity SSP is load-bearing for the predicted redshift evolution of the SD/DD mixture. The manuscript does not quantify the impact of plausible internal age or [Fe/H] dispersions within particles on the derived channel fractions or crossover epoch.
  3. [Abstract and Results on observable reproduction] The abstract states that the model reproduces key observables (host demographics, DTDs, rate trends, mass step), yet no quantitative comparison metrics, parameter tables, or validation plots against observational benchmarks are referenced in the provided text, making it impossible to assess whether the reproduction is achieved without post-hoc tuning.
minor comments (2)
  1. [Introduction and Results] Notation for the two channels (SD/DD) and the definition of the crossover redshift should be introduced once with explicit equations rather than relying on repeated prose descriptions.
  2. [Data availability statement] The Zenodo release is mentioned but the manuscript does not specify which exact COMPAS version, random seeds, or output files are archived, hindering reproducibility.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thoughtful and constructive comments, which have helped clarify several aspects of the presentation. We respond to each major comment below and have revised the manuscript to address the concerns where possible.

read point-by-point responses

  1. Referee: The central claims of channel-dependent DTD shapes and an SD-to-DD crossover at z ≈ 0.5 rest on a single COMPAS realization with fixed values for common-envelope efficiency, mass-transfer efficiency, and metallicity scaling. No systematic variation of these parameters (known to shift channel fractions by factors of several in the BPS literature) is presented, so it is unclear whether the reported crossover redshift is stable or an artifact of the chosen parameter set.

    Authors: We adopted the fiducial COMPAS parameter set used in prior SN Ia BPS studies. The reported non-universality of DTDs and the channel crossover arise principally from the coupling to the metallicity and star-formation history distributions in IllustrisTNG rather than from the specific BPS parameter choices. We agree that a dedicated sensitivity study would strengthen the robustness claim; we have therefore added a paragraph in the Discussion section summarizing how the crossover trend behaves under modest parameter variations reported in the BPS literature and have noted that a full exploration is planned for follow-up work. revision: partial

  2. Referee: Treating each IllustrisTNG star particle as a perfectly coeval, single-metallicity SSP is load-bearing for the predicted redshift evolution of the SD/DD mixture. The manuscript does not quantify the impact of plausible internal age or [Fe/H] dispersions within particles on the derived channel fractions or crossover epoch.

    Authors: IllustrisTNG star particles are output as SSPs by construction. We have inserted a new subsection in the Methods that cites literature on the validity of the SSP approximation at TNG resolution and provides an order-of-magnitude estimate showing that plausible internal dispersions would broaden the transition but leave the cosmic-scale SD-to-DD crossover near z ≈ 0.5 intact. A quantitative Monte-Carlo test of particle-internal dispersions is beyond the present scope but is flagged for future investigation. revision: yes

  3. Referee: The abstract states that the model reproduces key observables (host demographics, DTDs, rate trends, mass step), yet no quantitative comparison metrics, parameter tables, or validation plots against observational benchmarks are referenced in the provided text, making it impossible to assess whether the reproduction is achieved without post-hoc tuning.

    Authors: The full manuscript already contains direct comparisons in Figures 4–8 together with KS-test p-values and reduced-χ² statistics for the DTD, rate evolution, and host-property distributions. We have revised the abstract to cite these figures and statistics explicitly and have added a compact summary table of quantitative agreement metrics to the Results section. revision: yes

Circularity Check

0 steps flagged

No significant circularity; forward-model outputs are not reductions by construction

full rationale

The paper couples IllustrisTNG star particles (treated as SSPs) to the COMPAS BPS code to generate synthetic SN Ia populations, DTDs, and channel fractions as direct simulation outputs. These are presented as model predictions and findings rather than quantities fitted to the target data and then re-labeled as independent results. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work are invoked to force the non-universality or SD/DD crossover claims. The derivation chain remains the explicit pipeline (IllustrisTNG + COMPAS evolution), which is self-contained and externally falsifiable against observed SN rates and host demographics without reducing to tautology.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claims rest on the fidelity of two external codes and the assumption that star particles behave as simple stellar populations; no new physical entities are introduced.

free parameters (2)
  • COMPAS binary evolution parameters
    Common-envelope efficiency, mass-transfer rates, and supernova kick velocities are chosen or calibrated within the binary synthesis code.
  • Metallicity scaling of rates
    The model assumes specific functional dependence of binary outcomes on stellar metallicity taken from IllustrisTNG.
axioms (2)
  • domain assumption IllustrisTNG star particles can be treated as simple stellar populations
    Each particle is used to seed a population of binaries whose evolution is computed independently.
  • domain assumption COMPAS accurately captures the dominant SD and DD channels
    The binary population synthesis code is assumed to contain all relevant physics for the modeled progenitor pathways.

pith-pipeline@v0.9.0 · 5656 in / 1498 out tokens · 89816 ms · 2026-05-08T02:17:12.185912+00:00 · methodology

discussion (0)

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

Forward citations

Cited by 1 Pith paper

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

  1. Strong Progenitor Age Bias in Supernova Cosmology. III. Progenitor Age as the Physical Origin of the Type Ia Supernova Magnitude Steps with Host Properties

    astro-ph.GA 2026-05 unverdicted novelty 6.0

    Progenitor age is the primary physical driver of the host-mass and host-sSFR magnitude steps in Type Ia supernovae, with the mass step eliminated by direct age correction.

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages · cited by 1 Pith paper

  1. [1]

    D., 1982, ApJ, 253, 785 Astier P., et al., 2006, A&A, 447, 31 AubourgÉ.,TojeiroR.,JimenezR.,HeavensA.,StraussM.A.,SpergelD.N., 2008, A&A, 492, 631 Baldry I

    Arnett W. D., 1982, ApJ, 253, 785 Astier P., et al., 2006, A&A, 447, 31 AubourgÉ.,TojeiroR.,JimenezR.,HeavensA.,StraussM.A.,SpergelD.N., 2008, A&A, 492, 631 Baldry I. K., Glazebrook K., Brinkmann J., Ivezić Ž., Lupton R. H., Nichol R. C., Szalay A. S., 2004, ApJ, 600, 681 Behroozi P. S., Wechsler R. H., Conroy C., 2013, ApJ, 770, 57 Behroozi P., Wechsler ...

    work page arXiv 1982