An analysis of the Type Ia SN 2024gy and a comparison of different host extinction estimation techniques

Technical Summary: An analysis of the Type Ia SN 2024gy and a comparison of different host extinction estimation techniques

Problem Statement
Type Ia supernovae (SNe Ia) serve as standardisable candles for measuring cosmic distances, yet their utility is compromised by interstellar extinction caused by dust, which dims and reddens the observed light. Accurately correcting for this extinction requires estimating the column density of intervening material and the wavelength dependence of the extinction (parameterized by $R_V$). However, direct measurement of dust properties is often impossible, necessitating reliance on indirect tracers. Different estimation techniques—ranging from photometric light curve fitting to spectroscopic absorption line analysis—often yield disparate results, highlighting the challenge of accurately quantifying host galaxy extinction. This paper addresses this challenge by presenting a comprehensive, multi-wavelength analysis of the nearby, well-observed SN 2024gy to compare the efficacy of various extinction estimation methods.

Methodology
The authors present a heterogeneous dataset for SN 2024gy, located in the edge-on spiral galaxy NGC 4216 ($z \approx 0.001183$), collected from over a dozen ground- and space-based facilities. The dataset includes:

• Photometry: Multi-band optical light curves from ASAS-SN, ATLAS, GOTO, Pan-STARRS, and ZTF, covering the pre-peak, peak, and post-peak phases (up to ~400 days).
• Spectroscopy: 28 optical spectra and 5 near-infrared (NIR) spectra obtained from -15.2 to +316.8 days relative to peak brightness, utilizing instruments such as ALFOSC, FIES, X-Shooter, and NIR spectrographs.
• Polarimetry: Linear polarimetry in $gri$ bands at four epochs.

The analysis proceeds through several distinct technical pathways:

1. Light Curve Analysis: The authors fit the light curves using the SALT2 model to derive stretch ($x_1$) and colour ($c$) parameters. They also apply the "Lira law" (comparing the $g-r$ colour evolution to a dust-free sample) and the Bayesian framework BayeSN to estimate extinction.
2. Spectroscopic Analysis:
   • Line Measurements: Absorption features (Si II, Ca II NIR triplet, O I) were fitted to track velocity evolution. Late-time nebular spectra were used to measure the [Ni II]/[Fe II] flux ratio to infer the Ni/Fe abundance ratio.
   • ISM Tracers: High-resolution spectra (FIES, X-Shooter) were analyzed to detect narrow interstellar absorption lines (Na I D, Ca II H&K, K I, and Diffuse Interstellar Bands). Equivalent widths (EW) were converted to $E(B-V)$ using empirical relations (e.g., Poznanski et al. 2012; Munari & Zwitter 1997).
   • Modeling: Spectral synthesis was performed using the TARDIS radiative transfer code to compare the observed spectra against double-detonation and delayed-detonation explosion models.
3. Polarimetry: Stokes parameters were measured and corrected for Milky Way foreground polarization to isolate the host galaxy contribution, utilizing the Serkowski relation ($P_{max} \approx 9 \times E(B-V)$) to derive a lower limit on extinction.

Key Results

• Extinction Estimates: The study yields a range of host extinction estimates depending on the method used:
  • Photometric Methods: The Lira law, BayeSN, and spectral matching techniques generally converge on a mean $E(B-V)_{host} \approx 0.22 \pm 0.04$ mag (assuming $R_V = 3.1$).
  • Spectroscopic Methods: Estimates based on narrow ISM absorption lines show significant scatter. The mean of non-saturated lines (DIBs, K I) suggests a lower value of $E(B-V)_{host} = 0.12 \pm 0.02$ mag. The Na I D doublet yielded a higher, likely overestimated value ($0.26 \pm 0.09$ mag) due to line saturation in high-resolution data, which was not apparent in lower-resolution spectra.
  • Polarimetry: The observed polarization implies a lower limit of $E(B-V)_{host} \ge 0.08 \pm 0.01$ mag, consistent with the ISM origin of the polarization.
• Explosion Mechanism:
  • High Velocity (HV) Features: SN 2024gy exhibits a persistent, strong HV component in the Ca II NIR triplet (velocities $>30,000$ km s$^{-1}$) and Si II $\lambda 6355$ (phase $< -10$ days). TARDIS modeling suggests these features are better reproduced by double-detonation scenarios (specifically a sub-Chandrasekhar mass WD with a He-shell detonation) than by the standard delayed-detonation N100 model, which lacks the HV Ca II component.
  • Nucleosynthesis: Conversely, the measured stable Ni/Fe ratio ($n_{Ni}/n_{Fe} = 0.11 \pm 0.05$) in the nebular phase favors a delayed-detonation scenario involving a near-Chandrasekhar mass progenitor.
  • Conclusion on Mechanism: The authors conclude that no single explosion model is definitively favored; the data presents a mix of features supporting both double-detonation (HV Ca II) and delayed-detonation (Ni/Fe ratio) scenarios.
• Precursor Activity: Deep binning of pre-explosion data revealed no precursor activity brighter than $M \approx -8.75$ mag, and late-time light curves showed no flux excess indicative of circumstellar material (CSM) interaction.

Significance and Claims
The paper emphasizes that the spread in extinction estimates derived from different methods underscores the difficulty in accurately quantifying dust extinction. Key technical claims include:

1. Limitations of Na I D: The study demonstrates that Na I D can be an inadequate tracer for extinction in specific cases if the line is saturated, a condition that may be missed in lower-resolution spectra, leading to overestimation of $E(B-V)$.
2. Methodological Comparison: While photometric methods (Lira law, BayeSN) and spectral matching provide consistent results for this object, ISM line tracers show larger dispersion. The authors advocate for using weighted means of multiple weak lines (like DIBs and K I) when high-resolution data is available, rather than relying solely on strong lines like Na I D.
3. Polarimetry as a Lower Limit: The authors assert that while polarization data is difficult to obtain and reduce, it provides a robust, albeit non-constraining, lower limit on extinction, supporting the hypothesis that the observed polarization is induced by the host ISM rather than intrinsic to the SN.
4. Explosion Model Ambiguity: The conflicting evidence regarding the explosion mechanism (HV features vs. Ni/Fe ratio) suggests that SN 2024gy may not be easily categorized into a single theoretical framework, or that current models require further refinement to explain the coexistence of these features.

The authors do not claim to have resolved the explosion mechanism debate but rather highlight the complexity of interpreting multi-messenger data for individual SNe Ia and the necessity of cross-validating extinction estimates using multiple independent techniques.