pith. sign in

arxiv: 2605.23512 · v1 · pith:JBD4QH7Mnew · submitted 2026-05-22 · 🌌 astro-ph.GA · astro-ph.CO

Examining extinction distributions for type Ia supernovae in simulated 3D galaxies

Jo\~ao Duarte , Santiago Gonz\'alez-Gait\'an , Ana M. Mour\~ao , Rita P. Santos , Radoslaw Wojtak This is my paper

Pith reviewed 2026-05-25 03:50 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO
keywords type Ia supernovaedust extinctionextinction PDFradiative transfergalaxy morphologyintrinsic colorWeibull distributionhost galaxy properties
0
0 comments X

The pith

Two-parameter Weibull and exponentiated exponential distributions fit simulated type Ia supernova extinction more accurately than the standard exponential PDF across galaxy environments.

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

The paper generates synthetic extinction distributions for type Ia supernovae by running radiative transfer calculations through three-dimensional galaxy models with varying morphologies and dust content. It tests the common single-parameter exponential model against two-parameter generalizations and finds that the exponential form systematically underpredicts the number of low-extinction events while overpredicting high-extinction ones. The two-parameter forms capture the simulated distributions well enough that differences in host morphology and dust mass produce distinguishable shifts in separate PDF parameters. When the same two-parameter forms are applied to published SN color data under a Gaussian assumption for intrinsic color, the inferred mean intrinsic color moves redward by roughly two standard deviations relative to the exponential case.

Core claim

Simulations of SN Ia lines of sight through realistic 3D galaxies show that the extinction probability density function is not adequately described by the single-parameter exponential distribution, which underestimates low-extinction events and overestimates high-extinction ones. Two-parameter generalizations, particularly the Weibull and exponentiated exponential PDFs, parameterize the distributions more accurately in all tested environments. Variations in host morphology or dust mass affect different parameters of these two-parameter forms, allowing the two effects to be separated. Applying the two-parameter forms to observed SN colors from the literature, under the assumption of a normal

What carries the argument

Two-parameter generalizations of the exponential PDF (Weibull and exponentiated exponential) fitted to extinction values generated by SKIRT radiative transfer along lines of sight in simulated galaxies.

If this is right

  • The single-parameter exponential PDF cannot distinguish extinction distributions arising from different host morphologies or dust masses.
  • Two-parameter PDFs allow host morphology and dust mass to be disentangled through their separate effects on distinct parameters.
  • Use of two-parameter extinction PDFs in cosmological analyses shifts the inferred mean intrinsic color of SNe Ia redward by approximately 2 sigma.
  • SN simulation pipelines and distance-ladder calibrations that adopt the exponential form will misrepresent the frequency of low- and high-extinction events.

Where Pith is reading between the lines

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

  • Cosmological analyses that marginalize over extinction may obtain tighter constraints on dark energy parameters once the two-parameter forms replace the exponential.
  • Host-galaxy demographic studies could use the separate PDF parameters to map dust-mass and morphology trends across large SN samples without additional radiative-transfer runs.
  • Future high-resolution galaxy simulations with different dust grain models could test whether the preference for Weibull or exponentiated exponential persists or depends on microphysical assumptions.

Load-bearing premise

The SKIRT radiative transfer simulations with the chosen 3D galaxy models, dust properties, and morphologies accurately reproduce the true extinction distributions along lines of sight in real galaxies.

What would settle it

A large, homogeneous sample of observed SN Ia colors whose host-galaxy properties are also measured; if the best-fit extinction PDF parameters recovered from that sample match those from the simulations only when the exponential form is used, the central claim is falsified.

Figures

Figures reproduced from arXiv: 2605.23512 by Ana M. Mour\~ao, Jo\~ao Duarte, Radoslaw Wojtak, Rita P. Santos, Santiago Gonz\'alez-Gait\'an.

Figure 1
Figure 1. Figure 1: Distributions of AV for three samples of simulated SNe Ia, occurring in the spiral disk (left), spiral bulge (center), and elliptical galaxy environments (right). The distributions were obtained from radiative transfer simulated observations along random lines of sight, for host galaxies with dust mass MD = 107M⊙. Best-fit curves for the exponentiated exponential (green), exponential-logarithmic (orange), … view at source ↗
Figure 2
Figure 2. Figure 2: Best-fit parameters for the EE (left), EL (center) and W (right) PDFs as a function of host dust mass MD, with error bars defined by a 68% credible region. Results for two samples of SNe in spiral disks are shown, corresponding to host galaxies with dust disks of thickness h D z = 250pc (lighter squares) and h D z = 350pc (darker diamonds). The dashed lines represent the median α, θ or γ values for each fi… view at source ↗
Figure 3
Figure 3. Figure 3: Same as [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Same as [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Color distribution for Pantheon+SH0ES SNe Ia. Different best-fit curves to the color distribution are shown in black, assuming a Gaussian distribution for the intrinsic color cint (blue) and one of the following distributions for the color excess E(B − V): exponential PDF (red, top lef), exponentiated PDF (green, top right), exponential-logarithmic PDF (orange, bottom left), or Weibull PDF (purple, bottom … view at source ↗
Figure 6
Figure 6. Figure 6: For these reasons, we conclude that they o [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
read the original abstract

Dust extinction and reddening greatly contribute to type Ia supernovae (SNe Ia) observed color and magnitude variations. The models used to describe the extinction probability density function (PDF) are often simplistic, which can negatively impact SN simulations and cosmology. We present an analysis of simulated SN Ia extinction in galaxies along realistic lines of sight and investigate the parameterization of its PDF, as well as its dependence on host properties. We employed SKIRT, a radiative transfer code, to simulate observations of SNe Ia in different environments and generate synthetic extinction distributions. To parameterize and fit these distributions, we used both the commonly assumed single-parameter exponential PDF and some of its two-parameter generalizations. We find that the standard exponential PDF does not adequately describe simulated SN extinction: It underestimates low-extinction events and overestimates high-extinction ones. 2D KS tests show significant differences between the simulated extinction distributions for SNe in different environments, which the exponential parameterization cannot properly distinguish. In contrast, the two-parameter PDFs parameterize SN extinction distributions more accurately across all simulated environments. Variations in host morphology or dust mass relate to variations in different PDF parameters, meaning that the two effects can effectively be disentangled. We conclude that the two-parameter Weibull or exponentiated exponential PDFs offer the best parameterizations of SN Ia extinction for a wide range of simulated environments. Analyzing observed SN colors from the literature and assuming a Gaussian distribution for the intrinsic component, we conclude that a two-parameter extinction PDF results in intrinsically redder SNe, with their mean intrinsic color shifted ~2$\sigma$ in relation to the standard exponential extinction PDF.

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 employs SKIRT radiative-transfer simulations of Type Ia supernovae in 3D galaxy models to generate synthetic line-of-sight extinction (A_V) distributions. It shows via 2D KS tests that the standard single-parameter exponential PDF underfits low-extinction events and overfits high-extinction ones, while two-parameter generalizations (Weibull, exponentiated exponential) provide superior fits across varied host morphologies and dust masses. Reanalysis of literature SN colors under a Gaussian intrinsic-color assumption then yields a ~2σ redder mean intrinsic color when the two-parameter extinction PDF is adopted instead of the exponential.

Significance. If the simulated distributions faithfully represent real SN Ia hosts, the work supplies improved, environment-aware PDF parameterizations that could reduce systematic errors in SN Ia standardization and cosmology. The demonstration that morphology and dust-mass effects can be disentangled via distinct PDF parameters is a concrete advance. The manuscript also supplies machine-generated synthetic data and explicit statistical comparisons, which are strengths.

major comments (3)
  1. [Methods (galaxy models and dust setup)] The load-bearing step is the untested transfer from SKIRT runs to observed SN colors. No section compares the simulated A_V distributions against independent empirical constraints (e.g., extinction statistics from real SN Ia host samples or other radiative-transfer codes). Without such validation, both the PDF ranking and the inferred ~2σ intrinsic-color shift remain conditional on the adopted dust geometries, clumpiness, and grain properties.
  2. [Results (observed-color reanalysis)] The reanalysis of observed SN colors assumes a Gaussian intrinsic-color component whose mean and width are not varied; the abstract states the ~2σ shift but does not report a sensitivity test to this choice or to the precise functional form of the intrinsic distribution.
  3. [Results (KS tests)] Table or figure reporting the 2D KS p-values across environments: the claim that the exponential PDF “cannot properly distinguish” environments is central, yet the quantitative thresholds and sample sizes used for the tests are not stated, making it impossible to judge whether the reported significance is robust to binning or sample variance.
minor comments (2)
  1. [Abstract] Notation for the two-parameter PDFs (shape/scale parameters) should be defined once in the methods and used consistently; the abstract introduces “Weibull or exponentiated exponential” without explicit functional forms.
  2. [Figures] Figure captions for the synthetic extinction histograms should state the number of lines of sight per environment and the adopted dust-mass range so that readers can assess statistical power.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their constructive and detailed report. We address each major comment below with honest responses and indicate where revisions will be made to the manuscript.

read point-by-point responses

  1. Referee: [Methods (galaxy models and dust setup)] The load-bearing step is the untested transfer from SKIRT runs to observed SN colors. No section compares the simulated A_V distributions against independent empirical constraints (e.g., extinction statistics from real SN Ia host samples or other radiative-transfer codes). Without such validation, both the PDF ranking and the inferred ~2σ intrinsic-color shift remain conditional on the adopted dust geometries, clumpiness, and grain properties.

    Authors: We agree this is a substantive limitation. The manuscript is a simulation study whose primary goal is to derive and compare extinction PDF forms from realistic 3D SKIRT models; it does not claim direct empirical validation of the resulting A_V distributions against real SN Ia host samples. The ~2σ shift is presented explicitly as an illustration under the simulated PDFs. We will add a new subsection in the discussion that (i) states the conditional nature of the results, (ii) cites relevant observational extinction studies for context, and (iii) outlines how future work could perform such validation. No new simulations are added at this stage. revision: partial

  2. Referee: [Results (observed-color reanalysis)] The reanalysis of observed SN colors assumes a Gaussian intrinsic-color component whose mean and width are not varied; the abstract states the ~2σ shift but does not report a sensitivity test to this choice or to the precise functional form of the intrinsic distribution.

    Authors: The referee correctly identifies that the reanalysis fixes the Gaussian parameters. While the Gaussian form is the standard assumption in the SN Ia literature, we acknowledge that a sensitivity test would strengthen the claim. In the revised manuscript we will add a short sensitivity subsection (or appendix) that varies the mean and dispersion of the intrinsic-color Gaussian over a plausible range and reports the resulting change in the inferred mean intrinsic color and its significance. This will be done using the same observed color sample. revision: yes

  3. Referee: [Results (KS tests)] Table or figure reporting the 2D KS p-values across environments: the claim that the exponential PDF “cannot properly distinguish” environments is central, yet the quantitative thresholds and sample sizes used for the tests are not stated, making it impossible to judge whether the reported significance is robust to binning or sample variance.

    Authors: We agree that the quantitative details of the 2D KS tests must be provided. The current text states that the tests show significant differences but omits the actual p-values, effective sample sizes, and any checks on binning or variance. We will add a new table (or expanded figure caption) that reports the 2D KS p-values for each pair of environments, the number of lines of sight used per environment, and a brief statement on robustness to binning choices. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation relies on independent simulations and statistical fits

full rationale

The paper's chain proceeds from SKIRT radiative-transfer runs on chosen 3D galaxy models to generate line-of-sight extinction distributions, followed by direct fitting of single- and two-parameter PDFs and 2D KS tests on those simulated distributions, then a separate application to literature SN colors under an explicit Gaussian intrinsic-color assumption. No step equates a claimed result to its own inputs by construction, renames a fit as a prediction, or loads the central claim on a self-citation whose content is unverified. The simulations and tests are external to the final interpretive claim about observed colors, satisfying the criteria for a self-contained, non-circular analysis.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The analysis depends on the fidelity of radiative transfer simulations and the Gaussian assumption for intrinsic colors when interpreting observations; PDF parameters are fitted to simulation outputs.

free parameters (1)
  • shape and scale parameters of Weibull and exponentiated exponential PDFs
    Fitted directly to the simulated extinction distributions for different host environments.
axioms (1)
  • domain assumption Intrinsic color distribution of SNe Ia is Gaussian
    Invoked explicitly when analyzing observed SN colors to isolate the extinction component.

pith-pipeline@v0.9.0 · 5850 in / 1294 out tokens · 56639 ms · 2026-05-25T03:50:55.988580+00:00 · methodology

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

52 extracted references · 52 canonical work pages

  1. [1]

    P., James, P

    Anderson, J. P., James, P. A., Förster, F., et al. 2015, Monthly Notices of the Royal Astronomical Society, 448, 732

    work page 2015

  2. [2]

    2000, Monthly Notices of the Royal Astronomical Society, 318, 798

    Baes, M., Dejonghe, H., & de Rijcke, S. 2000, Monthly Notices of the Royal Astronomical Society, 318, 798

    work page 2000

  3. [3]

    M., & Zhu, Y

    Beifiori, A., Courteau, S., Corsini, E. M., & Zhu, Y . 2012, MNRAS, 419, 2497

    work page 2012

  4. [4]

    & Scolnic, D

    Brout, D. & Scolnic, D. 2021, ApJ, 909, 26

    work page 2021

  5. [5]

    & Charlot, S

    Bruzual, G. & Charlot, S. 2003, MNRAS, 344, 1000

    work page 2003

  6. [6]

    & Baes, M

    Camps, P. & Baes, M. 2015, Astronomy and Computing, 9, 20

    work page 2015

  7. [7]

    & Baes, M

    Camps, P. & Baes, M. 2020, Astronomy and Computing, 31, 100381

    work page 2020

  8. [8]

    P., Bianchi, S., et al

    Casasola, V ., Cassarà, L. P., Bianchi, S., et al. 2017, A&A, 605, A18

    work page 2017

  9. [9]

    2003, PASP, 115, 763

    Chabrier, G. 2003, PASP, 115, 763

    work page 2003

  10. [10]

    Commins, E. D. 2004, New Astronomy Reviews, 48, 567, proceedings of the Workshop on Supernovae and Dust De Geyter, G., Baes, M., Camps, P., et al. 2014, MNRAS, 441, 869

    work page 2004

  11. [11]

    2023, A&A, 680, A56

    Duarte, J., González-Gaitán, S., Mourão, A., et al. 2023, A&A, 680, A56

    work page 2023

  12. [12]

    2025, A&A, 700, A169

    Duarte, J., González-Gaitán, S., Mourão, A., et al. 2025, A&A, 700, A169

    work page 2025

  13. [13]

    W., Lang, D., & Goodman, J

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 González-Gaitán, S., de Jaeger, T., Galbany, L., et al. 2021, MNRAS, 508, 4656

    work page 2013

  14. [14]

    2007, A&A, 466, 11

    Guy, J., Astier, P., Baumont, S., et al. 2007, A&A, 466, 11

    work page 2007

  15. [15]

    2010, A&A, 523, A7

    Guy, J., Sullivan, M., Conley, A., et al. 2010, A&A, 523, A7

    work page 2010

  16. [16]

    Hallgren, L., Wojtak, R., Hjorth, J., & Steinhardt, C. L. 2026, A&A, 706, A129

    work page 2026

  17. [17]

    1998, The Astrophysical Journal, 502, 177

    Hatano, K., Branch, D., & Deaton, J. 1998, The Astrophysical Journal, 502, 177

    work page 1998

  18. [18]

    W., Reynolds, A., Smith, M., & Kraan-Korteweg, R

    Holwerda, B. W., Reynolds, A., Smith, M., & Kraan-Korteweg, R. C. 2014, Monthly Notices of the Royal Astronomical Society, 446, 3768

    work page 2014

  19. [19]

    G., & Kirshner, R

    Jha, S., Riess, A. G., & Kirshner, R. P. 2007, ApJ, 659, 122

    work page 2007

  20. [20]

    L., Hicken, M., Burke, D

    Kelly, P. L., Hicken, M., Burke, D. L., Mandel, K. S., & Kirshner, R. P. 2010, ApJ, 715, 743–756

    work page 2010

  21. [21]

    & Scolnic, D

    Kessler, R. & Scolnic, D. 2017, The Astrophysical Journal, 836, 56

    work page 2017

  22. [22]

    C., et al

    Lampeitl, H., Smith, M., Nichol, R. C., et al. 2010, ApJ, 722, 566–576 Le´sniewska, A., Michałowski, M. J., Gall, C., et al. 2023, ApJ, 953, 27

    work page 2010

  23. [23]

    S., Scolnic, D

    Mandel, K. S., Scolnic, D. M., Shariff, H., Foley, R. J., & Kirshner, R. P. 2017, The Astrophysical Journal, 842, 93

    work page 2017

  24. [24]

    S., Thorp, S., Narayan, G., Friedman, A

    Mandel, K. S., Thorp, S., Narayan, G., Friedman, A. S., & Avelino, A. 2022, MNRAS, 510, 3939

    work page 2022

  25. [25]

    P., Kessler, R., et al

    Marriner, J., Bernstein, J. P., Kessler, R., et al. 2011, ApJ, 740, 72

    work page 2011

  26. [26]

    Martins, G., González-Gaitán, S., Duarte, J., & Mourão, A. M. 2025, arXiv e- prints, arXiv:2511.14332 Michałowski, M. J., Hjorth, J., Gall, C., et al. 2019, A&A, 632, A43

    work page arXiv 2025

  27. [27]

    D., Sullivan, M., Balam, D., et al

    Neill, J. D., Sullivan, M., Balam, D., et al. 2006, The Astronomical Journal, 132, 1126

    work page 2006

  28. [28]

    1999, ApJS, 517, 565

    Perlmutter, S., Aldering, G., Goldhaber, G., et al. 1999, ApJS, 517, 565

    work page 1999

  29. [29]

    2012, The Astronomical Journal, 144, 59

    Perrett, K., Sullivan, M., Conley, A., et al. 2012, The Astronomical Journal, 144, 59

    work page 2012

  30. [30]

    Phillips, M. M. 1993, ApJ, 413, L105

    work page 1993

  31. [31]

    S., Zinchenko, I

    Pilyugin, L. S., Zinchenko, I. A., Lara-López, M. A., Nefedyev, Y . A., & Vílchez, J. M. 2021, A&A, 646, A54

    work page 2021

  32. [32]

    Plummer, H. C. 1911, Monthly Notices of the Royal Astronomical Society, 71, 460

    work page 1911

  33. [33]

    2023, The Astrophysical Jour- nal, 945, 84

    Popovic, B., Brout, D., Kessler, R., & Scolnic, D. 2023, The Astrophysical Jour- nal, 945, 84

    work page 2023

  34. [34]

    2021, The Astrophysi- cal Journal, 913, 49

    Popovic, B., Brout, D., Kessler, R., Scolnic, D., & Lu, L. 2021, The Astrophysi- cal Journal, 913, 49

    work page 2021

  35. [35]

    2024, MNRAS, 534, 2263

    Popovic, B., Wiseman, P., Sullivan, M., et al. 2024, MNRAS, 534, 2263

    work page 2024

  36. [36]

    2024, The Astronomical Jour- nal, 167, 131

    Pritchet, C., Thanjavur, K., Bottrell, C., & Gao, Y . 2024, The Astronomical Jour- nal, 167, 131

    work page 2024

  37. [37]

    & Patat, F

    Riello, M. & Patat, F. 2005, Monthly Notices of the Royal Astronomical Society, 362, 671

    work page 2005

  38. [38]

    G., Filippenko, A

    Riess, A. G., Filippenko, A. V ., Challis, P., et al. 1998, AJ, 116, 1009

    work page 1998

  39. [39]

    G., Yuan, W., Macri, L

    Riess, A. G., Yuan, W., Macri, L. M., et al. 2022, ApJ, 934, L7

    work page 2022

  40. [40]

    A., Riess, A

    Rodney, S. A., Riess, A. G., Strolger, L.-G., et al. 2014, AJ, 148, 13

    work page 2014

  41. [41]

    2012, MNRAS, 419, 2545

    Rowlands, K., Dunne, L., Maddox, S., et al. 2012, MNRAS, 419, 2545

    work page 2012

  42. [42]

    2022, The Astrophysical Journal, 938, 113 Sérsic, J

    Scolnic, D., Brout, D., Carr, A., et al. 2022, The Astrophysical Journal, 938, 113 Sérsic, J. L. 1963, Boletin de la Asociacion Argentina de Astronomia La Plata Argentina, 6, 41

    work page 2022

  43. [43]

    Smith, M. W. L., Gomez, H. L., Eales, S. A., et al. 2012, ApJ, 748, 123

    work page 2012

  44. [44]

    A., et al

    Sullivan, M., Conley, A., Howell, D. A., et al. 2010, MNRAS, 406, 782

    work page 2010

  45. [45]

    2021, Monthly Notices of the Royal Astronomical Society, 505, 2819

    Vincenzi, M., Sullivan, M., Graur, O., et al. 2021, Monthly Notices of the Royal Astronomical Society, 505, 2819

    work page 2021

  46. [46]

    E., et al

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261

    work page 2020

  47. [47]

    M., Dhawan, S., Mandel, K

    Ward, S. M., Dhawan, S., Mandel, K. S., Grayling, M., & Thorp, S. 2023, Monthly Notices of the Royal Astronomical Society, 526, 5715

    work page 2023

  48. [48]

    2021, MNRAS, 506, 3330

    Wiseman, P., Sullivan, M., Smith, M., et al. 2021, MNRAS, 506, 3330

    work page 2021

  49. [49]

    2022, MNRAS, 515, 4587

    Wiseman, P., Vincenzi, M., Sullivan, M., et al. 2022, MNRAS, 515, 4587

    work page 2022

  50. [50]

    & Hjorth, J

    Wojtak, R. & Hjorth, J. 2025, A&A, 702, A176

    work page 2025

  51. [51]

    Wojtak, R., Hjorth, J., & Hjortlund, J. O. 2023, Monthly Notices of the Royal Astronomical Society, 525, 5187

    work page 2023

  52. [52]

    Zubko, V ., Dwek, E., & Arendt, R. G. 2004, ApJS, 152, 211 Article number, page 12 of 13 J. Duarte: Examining extinction distributions for type Ia supernovae in simulated 3D galaxies Appendix A: Mass distributions for Dustpedia galaxies The dust mass distributions for both spiral and elliptical galaxies from the Dustpedia dataset are plotted in Fig. A. Fo...

    work page 2004