pith. machine review for the scientific record. sign in

arxiv: 2605.14081 · v1 · submitted 2026-05-13 · 🌌 astro-ph.SR · astro-ph.GA· astro-ph.HE

Recognition: no theorem link

SN2023ixf: ultraviolet-to-infrared radiative-transfer modeling of the nebular-phase evolution until 1000 days

Luc Dessart , Wynn V. Jacobson-Galan , K. Azalee Bostroem , Alexei V. Filippenko , Weikang Zheng , Thomas G. Brink , Lluis Galbany , Claudia Gutierrez
show 1 more author
Stefano Valenti
Authors on Pith no claims yet

Pith reviewed 2026-05-15 02:25 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GAastro-ph.HE
keywords dustejectamsunemissionmassfluxlinessn2023ixf
0
0 comments X

The pith

SN2023ixf's nebular light curve to 1000 days requires dust formation in the cold dense shell plus enhanced gamma-ray escape after 200 days.

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

The paper models the ultraviolet-to-infrared emission of supernova SN2023ixf from the nebular phase through 1000 days with non-local thermodynamic equilibrium radiative transfer. It uses the same ejecta parameters that fit the photospheric phase: 7-8 solar masses ejected at 1.2 times 10 to the 51 ergs with 0.05 solar masses of nickel-56, plus a 0.2 solar mass cold dense shell at 8000 km/s from circumstellar interaction. Matching the observed V-band decline needs both greater gamma-ray leakage and dust condensation starting after 200 days, first in the cold dense shell and later in the inner ejecta. Circumstellar interaction powers the ultraviolet flux at all times and the optical light after 600 days, while the cold dense shell slows to 6500 km/s by 998 days, implying ongoing mass growth. An external cold dust component is also needed to explain the mid-infrared, and the models produce cooler, dust-attenuated spectra compared with SN1993J.

Core claim

Non-local thermodynamic equilibrium radiative transfer models of SN2023ixf during the nebular phase to 1000 days, based on ejecta with 7-8 solar masses, 1.2 x 10^51 erg, 0.05 solar masses of nickel-56, and a 0.2 solar mass cold dense shell at 8000 km/s, match the data when enhanced gamma-ray escape and dust formation after 200 days are included. Dust forms first in the cold dense shell and subsequently in the inner ejecta, reaching up to 10^{-4} solar masses by 700 days, attenuating emission lines with possible blue-red asymmetries depending on location. The ultraviolet flux is powered by circumstellar material interaction, influenced by iron lines and strong blueshifted Ly alpha and Mg II

What carries the argument

The cold dense shell of 0.2 solar masses at 8000 km/s produced by ejecta-circumstellar material interaction, together with subsequent dust formation inside the shell and inner ejecta.

If this is right

  • The cold dense shell velocity decreases continuously from 8000 km/s at 112 days to 6500 km/s at 998 days, indicating mass growth of several 0.1 solar masses.
  • Dust mass in the cold dense shell and inner ejecta rises to possibly 10^{-4} solar masses at 700 days as a carbon-rich and silicon-rich mixture.
  • Emission lines are uniformly attenuated or show blue-red asymmetry depending on whether the dust lies interior or exterior to the line-forming gas.
  • Ultraviolet radiation remains largely unaffected by dust and is powered by circumstellar interaction at all epochs, with strengthening blueshifted Ly alpha and Mg II emission after 200 days.
  • The supernova is fainter and cooler than SN1993J at 1-3 years due to greater cold dense shell and ejecta masses plus dust extinction.

Where Pith is reading between the lines

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

  • Persistent circumstellar interaction through the nebular phase may be typical for Type IIb events and could bias nickel-mass estimates derived from late-time luminosity.
  • If the same ejecta parameters work from early through late phases, early-time data become more reliable for inferring progenitor structure without needing major late-time revisions.
  • The required external cold dust component implies either pre-existing circumstellar dust or dust formation outside the main ejecta, which could be tested with spatially resolved mid-infrared observations.
  • Greater gamma-ray escape at late times suggests the ejecta become progressively more transparent, potentially verifiable if future gamma-ray telescopes can detect such events.
  • keywords:[
  • supernova
  • SN2023ixf
  • nebular phase

Load-bearing premise

The ejecta mass, kinetic energy, and nickel-56 mass fixed from the photospheric phase continue to apply without modification through the entire nebular phase to 1000 days.

What would settle it

A direct measurement of cold dense shell velocity or dust mass at 300 days showing no deceleration or no dust increase would contradict the adjustments needed to match the V-band light curve.

Figures

Figures reproduced from arXiv: 2605.14081 by Alexei V. Filippenko, Claudia Gutierrez, K. Azalee Bostroem, Lluis Galbany, Luc Dessart, Stefano Valenti, Thomas G. Brink, Weikang Zheng, Wynn V. Jacobson-Galan.

Figure 1
Figure 1. Figure 1: Profiles of the density (main panel), composition (left inset), and volume filling factor (right inset) versus velocity at 112 d for the ejecta model x6p0 used in this work. SN 2023ixf, but here the focus is on the later evolution out to 1000 d. Because at nebular times the influence of interac￾tion power is significant (e.g., at 120 d) and eventually domi￾nates (at ∼ > 600 d) over decay power, only model … view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of gas and radiation properties for the interacting s15p2 model from Dessart et al. (2023) at 350 d (violet) and the interacting model x6p0 used in this work for SN 2023ixf at 320 d (red). Both models are subject to the same interaction power of 1040 erg s−1 , but differ in Ekin/Mej. From left to right and top to bottom, we show the mass density (uncorrected for clumping), the total energy depos… view at source ↗
Figure 4
Figure 4. Figure 4: Comparison of gas and radiation properties for the tailored model x6p0 with an interaction power of 1040 erg s−1 but injected as high-energy particles within the CDS at 8000 km s−1 (red) or in the form of X-rays throughout the outer ejecta and CDS regions (green). From left to right and top to bottom, we show the gas temperature, the H, C, and Mg ionization (zero corresponds to a neutral state), and finall… view at source ↗
Figure 5
Figure 5. Figure 5: shows a comparison between SN 2023ixf at 112, 147, and 175 d, and our model x6p0 with a decreasing interaction power from 4, 3, to 2×1040 erg s−1 . Although not obvious, inter￾action power alters the overall brightness of the model at those times, increasing the V-band brightness by about 0.3 mag (see Dessart et al. 2026), meaning that at the onset of the nebular phase, the SN optical brightness is not jus… view at source ↗
Figure 6
Figure 6. Figure 6: Comparison between UV, optical, and NIR observations of SN 2023ixf at 208 d (black) and model x6p0 with interaction power of 2 × 1040 erg s−1 , injected within the CDS (blue and red curves) or in the form of X-rays (green). Dust is included is one of the models (red curve). The top (bottom) panel used a log-linear (log-log) scale to better emphasize local and global agreements and mismatches. See Section 5… view at source ↗
Figure 7
Figure 7. Figure 7: Comparison between optical and IR observations of SN 2023ixf at 265 d (black) and model x6p0 with interaction power of 2 × 1040 erg s−1 , including dust (red) or not (blue). For the dusty model, the radiative-transfer calculations account for both cool dust in the inner ejecta and warm dust in the CDS, but a match to the IR emission requires an additional, external 600 K 3 × 10−4 M⊙ dust emission component… view at source ↗
Figure 8
Figure 8. Figure 8: Comparison between UV, optical, and NIR observations of SN 2023ixf at 300 d (black) and model x6p0 with interaction power of 1040 erg s−1 , injected within the CDS (blue and red curves) or in the form of X-rays (green). For the dusty model (red), the radiative-transfer calculations account for both cool dust in the inner ejecta and warm dust in the CDS. The top panel cuts the very strong Ly α for better vi… view at source ↗
Figure 9
Figure 9. Figure 9: Comparison between optical observations of SN 2023ixf at 329 d (black) and model x6p0 with interaction power of 1040 erg s−1 , including dust (red) or not (blue). We used the same dust parameters in the dusty model as at 300 d apart from an increase in the CDS dust mass to 2 × 10−5 M⊙. See Section 5.5 for discussion. 5.7. SN 2023ixf at 442 d [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Comparison between UV, optical, and IR observations of SN 2023ixf at 379 d (black) and model x6p0 with interaction power of 1040 erg s−1 , injected within the CDS (blue and red curves) or in the form of X-rays (green). The dusty model has a 10−5 M⊙ 500 K dust uni￾formly distributed in the inner 3000 km s−1 and 7 × 10−6 M⊙ 1000 K dust confined within the CDS. An external, optically thin dust component (6 ×… view at source ↗
Figure 11
Figure 11. Figure 11: Comparison between optical observations of SN 2023ixf at 442 d (black) and model x6p0 with an interaction power of 1040 erg s−1 , including dust (red) or not (blue). The dusty model has a 2 × 10−5 M⊙ 500 K dust uniformly distributed in the inner 3000 km s−1 and 3 × 10−5 M⊙ 500 K dust confined within the CDS. See Section 5.7 for dis￾cussion. the SN radiation arises preferentially from the CDS (i.e., locate… view at source ↗
Figure 12
Figure 12. Figure 12: Comparison between UV, optical, and IR observations of SN 2023ixf at 620 d (black) and model x6p0 with interaction power of 1040 erg s−1 , injected within the CDS (blue and red curves) or in the form of X-rays (green). The dusty model (red) has a 10−4 M⊙ 350 K dust uniformly distributed in the inner 3000 km s−1 and within the CDS. A full description of the bound-bound transitions contributing to the (dust… view at source ↗
Figure 13
Figure 13. Figure 13: Comparison between UV, optical, and IR observations of SN 2023ixf at 708 d (black) and model x6p0 with interaction power of 1040 erg s−1 , injected within the CDS (blue and red curves) or in the form of X-rays (green). The dusty model has a 10−4 M⊙ 350 K dust uni￾formly distributed in the inner 3000 km s−1 and 2 × 10−4 M⊙ 350 K dust confined within the CDS. For the IR emission, an additional contribution … view at source ↗
Figure 14
Figure 14. Figure 14: Comparison between optical observations of SN 2023ixf at 998 d (black) and model x6p0 at 1000 d with an interaction power of 5 × 1039 erg s−1 without dust. See Section 5.10 for discussion. (for details, see Sec. A). In [PITH_FULL_IMAGE:figures/full_fig_p013_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Comparison of the V-band light curve of SN 2023ixf with the x6p0 model having a time-variable interaction power (blue), as well as with allowance for dust at ∼ > 200 d (solid red). All models were pre￾sented by Dessart et al. (2026) as well as in Section 5. The model pho￾tometry was corrected for the distance and extinction to SN 2023ixf. The inset illustrates the evolution of quantities related to the ra… view at source ↗
Figure 16
Figure 16. Figure 16: Top: Illustration of the spatial origin of the emergent flux in the x6p0 model with interaction power but no dust at 388 d (see also [PITH_FULL_IMAGE:figures/full_fig_p015_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Comparison between the observations of SN 2023ixf (black) and the model with interaction power and allowance for dust (red) discussed in Section 5. Each row corresponds to one epoch, and each column corresponds to a wavelength range, with the NIR at left, the 7 µm region in the center (i.e., lines of [Ni ii] 6.634 µm, [Ar ii] 6.983 µm, and a blend of H i 7.458 µm and [Ni i] 7.505 µm), and the 10 µm region… view at source ↗
Figure 18
Figure 18. Figure 18: Top: Spectral comparison between the noninteracting model s15p2 at 375 d from Dessart et al. (2021b), used for comparison with SN 2024ggi by Dessart et al. (2025b), with the interacting model x6p0 at 388 d used in this work for SN 2023ixf. Here, the emergent luminosity has been scaled by the total decay power absorbed in each ejecta model to illustrate how interaction power influences primarily the UV. Bo… view at source ↗
Figure 19
Figure 19. Figure 19: Variation of spectral properties in the optical (left), the 7 µm region (center), and [Ne ii] 12.810 µm (right) resulting from changes in the clumping of the O-rich material (top row) or progenitor mass (lower three rows). All models are at 388 d and employ the same interaction power of 1040 erg s−1 and the same ejecta structure (e.g., density versus velocity) as model x6p0. et al. (2025) with more simula… view at source ↗
Figure 20
Figure 20. Figure 20: Comparison of the UV observations of SN 2023ixf at 200, 308, 619, and 723 d from Bostroem et al. (2025) and the model x6p0 with power injected within the CDS as high-energy electrons (red) or in the form of X-rays (green). The latter model ignores dust, although its effect on the UV spectrum is small. Observations and spectra have been rebinned to a fixed resolution of 100 km s−1 . no continuum flux to bl… view at source ↗
Figure 21
Figure 21. Figure 21: Wavelength dependence of the ejecta velocity V(τλ = 2/3) at which the radially (i.e., impact parameter p = 0) inward integrated total optical depth is 2/3 for model x6p0 with interaction power and shown as a function of time from about 200 to 760 d. A rebinning to a resolution of 600 (i.e., 500 km s−1 ) was applied to reduce the jaggedness of the curves. The abscissa covers from the UV, which is overall o… view at source ↗
Figure 22
Figure 22. Figure 22: Illustration of the cumulative flux integrated from 0.1 to 15 µm for model x6p0 with interaction power. The dashed gray line shows the result at 24 d (from Dessart et al. 2026). The thick colored lines correspond to the case wherein shock power is injected as high-energy electrons within the CDS, covering the nebular-phase evolution from 112 to 756 d. The thick dashed lines correspond to model counterpart… view at source ↗
Figure 23
Figure 23. Figure 23: Exploration of the impact on the emergent spectrum of various amounts of dust in the inner ejecta and in the CDS for the interacting model x6p0 at 265 d (see also Sec. 5.3 and [PITH_FULL_IMAGE:figures/full_fig_p020_23.png] view at source ↗
Figure 25
Figure 25. Figure 25: Same as [PITH_FULL_IMAGE:figures/full_fig_p021_25.png] view at source ↗
Figure 24
Figure 24. Figure 24: Evolution of the region encompassing the [O i] λλ 6300, 6364 and Hα profiles from 112 until 998 d for SN 2023ixf (top) as well as the best-matching models from Section 5 including the dust-free (middle) and dusty (bottom) counterpart (not all epochs coincide). The filled dots indicate the location where the flux is half that at the red edge of the Hα emission. The rest wavelength of Hα corresponds to 0 km… view at source ↗
Figure 26
Figure 26. Figure 26: Influence of changes in N, Mg, or Fe composition (top row and bottom-left panels) and CDS mass (bottom right) on the spectral properties of the dust-free x6p0 model with an interaction power of 1040 erg s−1 at 625 d. All spectra were rebinned to a resolution of 1000 (i.e., 300 km s−1 ). See Section 13 for discussion. and at ∼ 6500 km s−1 at 998 d (a filled dot indicates this veloc￾ity as the midpoint in f… view at source ↗
Figure 27
Figure 27. Figure 27: Comparison of SN 2023ixf with SN 1993J at about 380, 660, and 980 d after explosion. Each spectrum, which agrees within a few 0.01 mag of the corresponding observed V-band magnitude, is shown as luminosity after correcting for the distance and extinction, as well as the redshift. Additional scaling is applied to SN 1993J at the last two epochs for better visibility. The inset provides a zoom-in view of th… view at source ↗
read the original abstract

We present non-local thermodynamic equilibrium radiative-transfer modeling of SN2023ixf during the nebular phase out to 1000d, using the same ejecta that matched its photospheric evolution, namely a partially stripped red-supergiant star of initially 15Msun whose terminal explosion yielded ejecta with 7-8Msun, kinetic energy of 1.2e51erg, and 56Ni mass of 0.05Msun, augmented with a cold dense shell (CDS) of 0.2Msun at 8000km/s. Interaction with circumstellar material persists at all epochs, powering the ultraviolet (UV) flux at all times, but dominating the optical only after ~600d. Matching the V-band light curve requires invoking both enhanced gamma-ray escape and dust formation after ~200d, first in the CDS and eventually in the inner ejecta as well. Depending on where they form relative to the dust, emission lines are uniformly attenuated or skewed with a blue-red asymmetry. Our models suggest a rising dust mass (chosen as an C-rich and Si-rich mixture) in the CDS and inner ejecta, possibly reaching 1e-4Msun at 700d, while an external cold dust component is required to match the mid-infrared emission. The UV radiation, largely unaffected by dust, is influenced by the emission and absorption from Fe lines, together with strong, blueshifted emission from Lyalpha and MgII2800, both present at >~200d and with a strengthening fractional flux thereafter. Optical-depth effects play a critical role for the UV flux, and most notably on Lyalpha whose strength depends strongly on the CDS structure (mass and extent) and the treatment of power injection. The CDS is continuously slowing down from 8000km/s at 112d to ~6500km/s at 998d, suggesting a growth in mass of several 0.1Msun. SN2023ixf shares many similarities with SN1993J at 1-3yr, but it is eventually fainter due to dust extinction and cooler (i.e., weak [NII] and no [OIII] lines) likely as a result of greater CDS and ejecta masses.

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

2 major / 2 minor

Summary. The manuscript presents non-LTE radiative-transfer modeling of SN 2023ixf in the nebular phase out to 1000 days, employing the same ejecta parameters (7–8 M⊙, 1.2 × 10^51 erg, 0.05 M⊙ of 56Ni) that matched the photospheric phase, augmented by a 0.2 M⊙ cold dense shell (CDS) at 8000 km s⁻¹. Persistent CSM interaction is claimed to power the UV flux at all epochs and the optical after ~600 d. Matching the V-band decline is stated to require enhanced γ-ray escape plus dust formation (C-rich/Si-rich mixture) after ~200 d, first in the CDS and later in the inner ejecta, with an external cold dust component needed for the mid-IR; the CDS is reported to decelerate, implying mass growth of several 0.1 M⊙. The models also address UV line formation (Lyα, Mg II) and compare the event to SN 1993J.

Significance. If the necessity of the additional components can be shown to be robust rather than compensatory, the work would provide a valuable multi-epoch radiative-transfer framework for interacting supernovae, quantifying the roles of CDS dust formation, γ-ray leakage, and ongoing mass accretion in shaping nebular light curves and spectra. It would also strengthen the observational link between SN 2023ixf and events like SN 1993J while highlighting how dust extinction and CDS structure affect UV and optical diagnostics at late times.

major comments (2)
  1. [Abstract] Abstract: The claim that matching the V-band light curve 'requires' both enhanced gamma-ray escape and dust formation after ~200 d is load-bearing for the central conclusion, yet rests on holding the photospheric ejecta parameters (7–8 M⊙, 1.2 × 10^51 erg, 0.05 M⊙ 56Ni) exactly fixed. No sensitivity study or quantitative error analysis is presented to show that modest re-adjustments (e.g., ±0.01 M⊙ in 56Ni or small changes in energy redistribution) fail to reproduce the decline without these additions.
  2. [Abstract] Abstract: The reported CDS deceleration from 8000 km s⁻¹ at 112 d to ~6500 km s⁻¹ at 998 d is interpreted as evidence for mass growth of several 0.1 M⊙. Because CDS mass, velocity, and dust mass are themselves free parameters adjusted to fit the light curve and line profiles, this physical interpretation is circular; the fitted parameters define the conclusion rather than providing an independent test.
minor comments (2)
  1. [Abstract] Abstract: Numerical notation such as '1.2e51erg' and '1e-4Msun' should be rendered in standard scientific format (e.g., 1.2 × 10^{51} erg) for consistency with journal style.
  2. The manuscript would benefit from explicit reporting of goodness-of-fit metrics (e.g., reduced χ² or residual statistics) for the model–data comparisons rather than qualitative statements of agreement.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment below, providing clarifications on our modeling choices and noting revisions made to strengthen the presentation.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The claim that matching the V-band light curve 'requires' both enhanced gamma-ray escape and dust formation after ~200 d is load-bearing for the central conclusion, yet rests on holding the photospheric ejecta parameters (7–8 M⊙, 1.2 × 10^51 erg, 0.05 M⊙ 56Ni) exactly fixed. No sensitivity study or quantitative error analysis is presented to show that modest re-adjustments (e.g., ±0.01 M⊙ in 56Ni or small changes in energy redistribution) fail to reproduce the decline without these additions.

    Authors: The ejecta parameters were deliberately held fixed to maintain consistency with our earlier photospheric-phase modeling of SN 2023ixf, which already constrained the 7–8 M⊙ mass, 1.2 × 10^51 erg energy, and 0.05 M⊙ 56Ni mass. The nebular models demonstrate that these fixed values cannot reproduce the observed V-band decline rate or the UV-optical spectral evolution without the additional components. We agree that a quantitative sensitivity study would be valuable; we have added a short paragraph in the revised discussion section exploring the effects of modest (±0.01 M⊙) variations in 56Ni and noting that they do not remove the requirement for enhanced γ-ray escape and dust formation. revision: partial

  2. Referee: [Abstract] Abstract: The reported CDS deceleration from 8000 km s⁻¹ at 112 d to ~6500 km s⁻¹ at 998 d is interpreted as evidence for mass growth of several 0.1 M⊙. Because CDS mass, velocity, and dust mass are themselves free parameters adjusted to fit the light curve and line profiles, this physical interpretation is circular; the fitted parameters define the conclusion rather than providing an independent test.

    Authors: The CDS velocity at each epoch is not a free parameter but is directly constrained by matching the observed Doppler shifts and line-profile shapes in the multi-epoch spectra. The model then evolves the CDS under ongoing CSM interaction, producing the deceleration that matches the data. This velocity evolution is therefore an independent spectral diagnostic that supports the inferred mass growth, rather than being defined solely by the light-curve fit. We have clarified this distinction in the revised methods and results sections. revision: yes

Circularity Check

2 steps flagged

CDS mass/velocity and dust mass fitted to match V-band and line profiles; deceleration then interpreted as mass growth, with 'requires' language applied to the fitted adjustments

specific steps
  1. fitted input called prediction [Abstract]
    "augmented with a cold dense shell (CDS) of 0.2Msun at 8000km/s. ... Matching the V-band light curve requires invoking both enhanced gamma-ray escape and dust formation after ~200d, first in the CDS and eventually in the inner ejecta as well. ... Our models suggest a rising dust mass ... possibly reaching 1e-4Msun at 700d"

    CDS mass, initial velocity, and dust mass (C-rich/Si-rich mixture) are chosen/adjusted to reproduce the light curve and profiles; the statement that matching 'requires' enhanced escape and dust formation therefore follows by construction from those adjustments rather than emerging as a prediction.

  2. self definitional [Abstract]
    "The CDS is continuously slowing down from 8000km/s at 112d to ~6500km/s at 998d, suggesting a growth in mass of several 0.1Msun."

    The model is constructed with a fixed initial CDS mass of 0.2 Msun whose velocity evolution is computed under interaction; the observed deceleration is then interpreted as implying additional mass growth of several 0.1 Msun, making the physical conclusion a direct re-description of the fitted velocity trajectory.

full rationale

The derivation fixes photospheric ejecta (7-8 Msun, 1.2e51 erg, 0.05 Msun 56Ni) and augments with a CDS (0.2 Msun at 8000 km/s), then adjusts gamma-ray escape and dust mass to reproduce the observed V-band decline and line profiles. The resulting model velocity slowdown is presented as evidence for additional mass growth. This reduces the central claim (that enhanced escape and dust 'require' invocation) to a direct consequence of the fitting choices rather than an independent prediction. No sensitivity study on parameter variations is reported to test uniqueness.

Axiom & Free-Parameter Ledger

5 free parameters · 2 axioms · 3 invented entities

The model depends on numerous fitted quantities (ejecta mass, kinetic energy, nickel mass, CDS mass and velocity, dust composition and growth rate) chosen to match both photospheric and nebular data; standard non-LTE radiative transfer and the existence of a CDS are taken as given without new justification.

free parameters (5)
  • ejecta mass
    Set to 7-8 solar masses to match earlier photospheric modeling
  • kinetic energy
    Fixed at 1.2e51 erg from photospheric phase
  • 56Ni mass
    0.05 solar masses chosen to power the light curve
  • CDS mass
    0.2 solar masses at 8000 km/s introduced to power late UV and optical
  • dust mass
    Rising to 1e-4 solar masses at 700 days, C-rich and Si-rich mixture
axioms (2)
  • domain assumption Non-local thermodynamic equilibrium radiative transfer accurately describes the nebular-phase emission
    Invoked throughout the modeling without new validation
  • domain assumption The same ejecta structure from the photospheric phase remains valid at 1000 days
    Central modeling choice stated in the abstract
invented entities (3)
  • Cold dense shell (CDS) no independent evidence
    purpose: To sustain UV flux and explain optical light curve after 600 days
    Postulated 0.2 solar-mass shell at 8000 km/s whose deceleration is observed
  • Internal dust component no independent evidence
    purpose: To attenuate optical lines and explain V-band decline after 200 days
    Formed first in CDS then in inner ejecta; mass chosen to fit data
  • External cold dust component no independent evidence
    purpose: To account for mid-infrared emission
    Required separately from internal dust

pith-pipeline@v0.9.0 · 5780 in / 1932 out tokens · 115353 ms · 2026-05-15T02:25:52.085188+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

104 extracted references · 104 canonical work pages · 1 internal anchor

  1. [1]

    M., et al

    Baron, E., Ashall, C., DerKacy, J. M., et al. 2025, ApJ, 994, 249

    work page 2025

  2. [2]

    C., Orellana, M., Folatelli, G., et al

    Bersten, M. C., Orellana, M., Folatelli, G., et al. 2024, A&A, 681, L18

    work page 2024

  3. [3]

    J., & Milisavljevic, D

    Bevan, A., Barlow, M. J., & Milisavljevic, D. 2017, MNRAS, 465, 4044

    work page 2017

  4. [4]

    I., Eastman, R., Bartunov, O

    Blinnikov, S. I., Eastman, R., Bartunov, O. S., Popolitov, V . A., & Woosley, S. E. 1998, ApJ, 496, 454

    work page 1998

  5. [5]

    A., Pearson, J., Shrestha, M., et al

    Bostroem, K. A., Pearson, J., Shrestha, M., et al. 2023, ApJ, 956, L5

    work page 2023

  6. [6]

    A., Sand, D

    Bostroem, K. A., Sand, D. J., Dessart, L., et al. 2024, ApJ, 973, L47

    work page 2024

  7. [7]

    A., Valenti, S., Sand, D

    Bostroem, K. A., Valenti, S., Sand, D. J., et al. 2025, submitted to ApJ, arXiv:2508.11756

    work page arXiv 2025

  8. [8]

    A., Laplace, E., Schneider, F

    Bronner, V . A., Laplace, E., Schneider, F. R. N., & Podsiadlowski, P. 2025, A&A, 703, A61

    work page 2025

  9. [9]

    A., Maeda, K., Ray, A

    Chandra, P., Chevalier, R. A., Maeda, K., Ray, A. K., & Nayana, A. J. 2024, ApJ, 963, L4

    work page 2024

  10. [10]

    & Blondin, J

    Chevalier, R. & Blondin, J. M. 1995, ApJ, 444, 312

    work page 1995

  11. [11]

    A., Blondin, J

    Chevalier, R. A., Blondin, J. M., & Emmering, R. T. 1992, ApJ, 392, 118

    work page 1992

  12. [12]

    M., Ashall, C., Baron, E., et al

    DerKacy, J. M., Ashall, C., Baron, E., et al. 2026, ApJ, 997, 179

    work page 2026

  13. [13]

    2024, A&A, 692, A204

    Dessart, L. 2024, A&A, 692, A204

    work page 2024

  14. [14]

    2025, A&A, 703, A287

    Dessart, L. 2025, A&A, 703, A287

    work page 2025

  15. [15]

    & Audit, E

    Dessart, L. & Audit, E. 2019, A&A, 629, A17

    work page 2019

  16. [16]

    P., Kuncarayakti, H., Fox, O

    Dessart, L., Gutiérrez, C. P., Kuncarayakti, H., Fox, O. D., & Filippenko, A. V . 2023, A&A, 675, A33

    work page 2023

  17. [17]

    & Hillier, D

    Dessart, L. & Hillier, D. J. 2019, A&A, 625, A9

    work page 2019

  18. [18]

    & Hillier, D

    Dessart, L. & Hillier, D. J. 2020, A&A, 643, L13

    work page 2020

  19. [19]

    & Hillier, D

    Dessart, L. & Hillier, D. J. 2022, A&A, 660, L9

    work page 2022

  20. [20]

    J., & Audit, E

    Dessart, L., Hillier, D. J., & Audit, E. 2017, A&A, 605, A83

    work page 2017

  21. [21]

    J., & Wilk, K

    Dessart, L., Hillier, D. J., & Wilk, K. D. 2018, A&A, 619, A30

    work page 2018

  22. [22]

    & Jacobson-Galán, W

    Dessart, L. & Jacobson-Galán, W. V . 2023, A&A, 677, A105

    work page 2023

  23. [23]

    SN2023ixf: Radiative-transfer modeling of the photospheric phase evolution from the ultraviolet to the infrared

    Dessart, L., Jacobson-Galan, W. V ., Bostroem, K. A., et al. 2026, arXiv:2604.05145

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  24. [24]

    2025, ApJ, 984, 71

    Dickinson, D., Milisavljevic, D., Garretson, B., et al. 2025, ApJ, 984, 71

    work page 2025

  25. [25]

    J., Valenti, S., et al

    Dong, Y ., Sand, D. J., Valenti, S., et al. 2023, ApJ, 957, 28

    work page 2023

  26. [26]

    Draine, B. T. & Li, A. 2007, ApJ, 657, 810

    work page 2007

  27. [27]

    J., et al

    Fabbri, J., Otsuka, M., Barlow, M. J., et al. 2011, MNRAS, 418, 1285

    work page 2011

  28. [28]

    J., Ferrari, L., et al

    Fang, Q., Moriya, T. J., Ferrari, L., et al. 2025, ApJ, 978, 36

    work page 2025

  29. [29]

    Ferrari, L., Folatelli, G., Ertini, K., Kuncarayakti, H., & Andrews, J. E. 2024, A&A, 687, L20

    work page 2024

  30. [30]

    Filippenko, A. V . 1997, ARA&A, 35, 309

    work page 1997

  31. [31]

    V ., Matheson, T., & Ho, L

    Filippenko, A. V ., Matheson, T., & Ho, L. C. 1993, ApJ, 415, L103

    work page 1993

  32. [32]

    L., Hillier, D

    Flores, B. L., Hillier, D. J., & Dessart, L. 2023, MNRAS, 518, 5001

    work page 2023

  33. [33]

    2025, A&A, 698, A213

    Folatelli, G., Ferrari, L., Ertini, K., Kuncarayakti, H., & Maeda, K. 2025, A&A, 698, A213

    work page 2025

  34. [34]

    & Goldberg, J

    Forde, S. & Goldberg, J. A. 2025, Research Notes of the American Astronomical Society, 9, 135

    work page 2025

  35. [35]

    M., Chevalier, R

    Fransson, C., Challis, P. M., Chevalier, R. A., et al. 2005, ApJ, 622, 991

    work page 2005

  36. [36]

    A., Bildsten, L., & Paxton, B

    Goldberg, J. A., Bildsten, L., & Paxton, B. 2019, ApJ, 879, 3

    work page 2019

  37. [37]

    W., Brightman, M., Earnshaw, H

    Grefenstette, B. W., Brightman, M., Earnshaw, H. P., Harrison, F. A., & Margutti, R. 2023, ApJ, 952, L3

    work page 2023

  38. [38]

    Hillier, D. J. & Dessart, L. 2012, MNRAS, 424, 252

    work page 2012

  39. [39]

    Hillier, D. J. & Dessart, L. 2019, A&A, 631, A8

    work page 2019

  40. [40]

    Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407

    work page 1998

  41. [41]

    A., et al

    Hsu, B., Smith, N., Goldberg, J. A., et al. 2025, ApJ, 990, 148

    work page 2025

  42. [42]

    2023, Transient Name Server Discovery Report, 2023-1158, 1

    Itagaki, K. 2023, Transient Name Server Discovery Report, 2023-1158, 1

    work page 2023

  43. [43]

    2025, ApJ, 978, 138 Jacobson-Galán, W

    Iwata, Y ., Akimoto, M., Matsuoka, T., et al. 2025, ApJ, 978, 138 Jacobson-Galán, W. V ., Dessart, L., Davis, K. W., et al. 2025a, ApJ, 992, 100 Jacobson-Galán, W. V ., Dessart, L., Kilpatrick, C. D., et al. 2025b, ApJ, 994, L14 Jacobson-Galán, W. V ., Dessart, L., Margutti, R., et al. 2023, ApJ, 954, L42

    work page 2025

  44. [44]

    E., Pearson, J., Beasor, E

    Jencson, J. E., Pearson, J., Beasor, E. R., et al. 2023, ApJ, 952, L30

    work page 2023

  45. [45]

    2012, A&A, 546, A28

    Jerkstrand, A., Fransson, C., Maguire, K., et al. 2012, A&A, 546, A28

    work page 2012

  46. [46]

    D., Foley, R

    Kilpatrick, C. D., Foley, R. J., Jacobson-Galán, W. V ., et al. 2023, ApJ, 952, L23

    work page 2023

  47. [47]

    2006, ApJ, 651, L117

    Kotak, R., Meikle, P., Pozzo, M., et al. 2006, ApJ, 651, L117

    work page 2006

  48. [48]

    D., Höflich, P

    Kotak, R., Meikle, P., van Dyk, S. D., Höflich, P. A., & Mattila, S. 2005, ApJ, 628, L123

    work page 2005

  49. [49]

    Kotak, R., Meikle, W. P. S., Farrah, D., et al. 2009, ApJ, 704, 306

    work page 2009

  50. [50]

    2025, A&A, 694, A319

    Kozyreva, A., Caputo, A., Baklanov, P., Mironov, A., & Janka, H.-T. 2025, A&A, 694, A319

    work page 2025

  51. [51]

    R., Singleton, A

    Kumar, A., Dastidar, R., Maund, J. R., Singleton, A. J., & Sun, N.-C. 2025, MNRAS, 538, 659

    work page 2025

  52. [52]

    A., Schneider, F

    Laplace, E., Bronner, V . A., Schneider, F. R. N., & Podsiadlowski, P. 2026, ApJ, 998, L40

    work page 2026

  53. [53]

    C., Filippenko, A

    Leonard, D. C., Filippenko, A. V ., Gates, E. L., et al. 2002, PASP, 114, 35 Article number, page 25 of 29 A&A proofs:manuscript no. ms

    work page 2002

  54. [54]

    2024, Nature, 627, 754

    Li, G., Hu, M., Li, W., et al. 2024, Nature, 627, 754

    work page 2024

  55. [55]

    2020, A&A, 642, A135

    Liljegren, S., Jerkstrand, A., & Grumer, J. 2020, A&A, 642, A135

    work page 2020

  56. [56]

    & Dalgarno, A

    Liu, W. & Dalgarno, A. 1995, ApJ, 454, 472

    work page 1995

  57. [57]

    B., Danziger, I

    Lucy, L. B., Danziger, I. J., Gouiffes, C., & Bouchet, P. 1989, in IAU Colloq. 120: Structure and Dynamics of the Interstellar Medium, ed. G. Tenorio-Tagle, M. Moles, & J. Melnick, V ol. 350, 164

    work page 1989

  58. [58]

    2025, submitted to ApJL, arXiv:2510.14875

    Ma, J.-Z., Justham, S., Pakmor, R., et al. 2025, submitted to ApJL, arXiv:2510.14875

    work page arXiv 2025

  59. [59]

    V ., Barth, A

    Matheson, T., Filippenko, A. V ., Barth, A. J., et al. 2000, AJ, 120, 1487

    work page 2000

  60. [60]

    J., & Dessart, L

    McLeod, C., Hillier, D. J., & Dessart, L. 2024, MNRAS, 532, 549

    work page 2024

  61. [61]

    2025, ApJ, 993, 191

    Medler, K., Ashall, C., Hoeflich, P., et al. 2025, ApJ, 993, 191

    work page 2025

  62. [62]

    Meikle, W. P. S., Kotak, R., Farrah, D., et al. 2011, ApJ, 732, 109

    work page 2011

  63. [63]

    2026, ApJ, 997, 330

    Mera, T., Ashall, C., Hoeflich, P., et al. 2026, ApJ, 997, 330

    work page 2026

  64. [64]

    2015, ApJ, 815, 120

    Milisavljevic, D., Margutti, R., Kamble, A., et al. 2015, ApJ, 815, 120

    work page 2015

  65. [65]

    J., Chevalier, R

    Milisavljevic, D., Patnaude, D. J., Chevalier, R. A., et al. 2018, ApJ, 864, L36

    work page 2018

  66. [66]

    Moriya, T. J. & Singh, A. 2024, PASJ, 76, 1050

    work page 2024

  67. [67]

    J., Margutti, R., Wiston, E., et al

    Nayana, A. J., Margutti, R., Wiston, E., et al. 2025, ApJ, 985, 51

    work page 2025

  68. [68]

    Neustadt, J. M. M., Kochanek, C. S., & Smith, M. R. 2024, MNRAS, 527, 5366

    work page 2024

  69. [69]

    J., Bevan, A., et al

    Niculescu-Duvaz, M., Barlow, M. J., Bevan, A., et al. 2022, MNRAS, 515, 4302

    work page 2022

  70. [70]

    R., et al

    Niu, Z., Sun, N.-C., Maund, J. R., et al. 2023, ApJ, 955, L15

    work page 2023

  71. [71]

    1993, Nature, 364, 507

    Nomoto, K., Suzuki, T., Shigeyama, T., et al. 1993, Nature, 364, 507

    work page 1993

  72. [72]

    H., Rho, J., Yoon, S.-C., et al

    Park, S. H., Rho, J., Yoon, S.-C., et al. 2025, A&A, 703, A227

    work page 2025

  73. [73]

    J., et al

    Pearson, J., Subrayan, B., Sand, D. J., et al. 2025, ApJ, 993, 213

    work page 2025

  74. [74]

    A., Gal-Yam, A., Irani, I., & Zimmerman, E

    Perley, D. A., Gal-Yam, A., Irani, I., & Zimmerman, E. 2023, Transient Name Server AstroNote, 119, 1

    work page 2023

  75. [75]

    Podsiadlowski, P., Hsu, J. J. L., Joss, P. C., & Ross, R. R. 1993, Nature, 364, 509

    work page 1993

  76. [76]

    X., Hennawi, J., Cooke, R., et al

    Prochaska, J. X., Hennawi, J., Cooke, R., et al. 2020, pypeit/PypeIt: Release 1.0.0

    work page 2020

  77. [77]

    X., Hennawi, J

    Prochaska, J. X., Hennawi, J. F., Westfall, K. B., et al. 2020, Journal of Open Source Software, 5, 2308

    work page 2020

  78. [78]

    2024, MNRAS, 534, 271

    Qin, Y .-J., Zhang, K., Bloom, J., et al. 2024, MNRAS, 534, 271

    work page 2024

  79. [79]

    L., Villar, V

    Ransome, C. L., Villar, V . A., Tartaglia, A., et al. 2024, ApJ, 965, 93

    work page 2024

  80. [80]

    W., Treffers, R

    Richmond, M. W., Treffers, R. R., Filippenko, A. V ., & Paik, Y . 1996, AJ, 112, 732

    work page 1996

Showing first 80 references.