Estimating the Explosion Time of Core-Collapse Supernovae from Their Optical Light Curves

Core-collapse supernovae are among the prime candidate sources of high energy neutrinos. Accordingly, the IceCube collaboration has started a program to search for such a signal. IceCube operates an online search for neutrino bursts, forwarding the directions of candidate events to a network of optical telescopes for immediate follow-up observations. If a supernova is identified from the optical observations, in addition to a directional coincidence a temporal photon-neutrino coincidence also needs to be established. To achieve this, we present a method for estimating the supernova explosion time from its light curve using a simple model. We test the model with supernova light curve data from SN1987A, SN2006aj and SN2008D and show that the explosion times can be determined with an accuracy of better than a few hours.

💡 Research Summary

The paper addresses a critical need in multimessenger astronomy: determining the explosion time (t₀) of core‑collapse supernovae (SNe) from their optical light curves with sufficient precision to enable a robust temporal coincidence with high‑energy neutrino detections by IceCube. IceCube’s real‑time neutrino‑burst alert system forwards the direction of candidate events to optical telescopes (e.g., ROTSE, PTF). While a directional coincidence reduces background, a temporal coincidence is essential because the background rate of accidental neutrino doublets is non‑negligible. Previously, analyses have assumed that the explosion time can be known to about one day, but this assumption lacked observational validation.

The authors propose a simple analytic model for the early optical light curve, consisting of two physically motivated phases: (1) the shock‑breakout phase, modeled after Waxman et al., where the emission is approximated as a blackbody with radius r ∝ Δt^0.8 and temperature T ∝ Δt^‑0.5 (Δt = t − t₀). This yields a flux expression Φ_BB = a₁ exp(a₂ Δt^0.5) − 1 · Δt^1.6, with a₁, a₂, and t₀ as free parameters. (2) The subsequent expansion phase, for which two alternatives are tested: a simple t² dependence (Φ ∝ a₃ Δt²) corresponding to a constant‑temperature photosphere expanding at constant velocity, and the more sophisticated Arnett model that solves a time‑dependent diffusion equation and includes ⁵⁶Ni heating.

The method is applied to three well‑studied SNe with independently known explosion times: SN 2006aj (Swift UVOT U, B, V data), SN 2008D (FLWO B, V, R, I data), and SN 1987A (Hamuy et al. B, V, R, I data). For each SN, the authors fit the early light‑curve points (typically within the first few days after explosion) using the combined shock‑breakout plus expansion model, allowing each photometric band its own set of parameters.

Key results:

• SN 2006aj: The fitted t₀ is on average –0.04 days relative to the X‑ray flash, with a statistical uncertainty of ±0.005 days. Both the t² and Arnett expansion models give essentially identical t₀, indicating that the early shock‑breakout dominates the timing constraint.
• SN 2008D: Three of the four bands (B, R, I) yield t₀ consistent with zero offset, while the V‑band gives a later value of 0.24 ± 0.08 days. The simple t² expansion model provides the best overall fit (χ²/NDF = 15.9/16), whereas the Arnett model yields a significantly larger χ². The average statistical error across bands is ≈0.06 days (≈1.5 hours).
• SN 1987A: Using data starting 1.14 days after the neutrino burst, the fit returns a t₀ with a statistical error of ≈0.2 days and a systematic shift of about 0.3 days, reflecting the crudeness of the model for a Type II‑P light curve. The redder bands (R, I) show larger deviations because the shock‑breakout signature is weaker at longer wavelengths.

The authors also explore the impact of early data acquisition. By artificially removing the earliest points, they demonstrate that the uncertainty on t₀ grows rapidly, confirming that observations within the first few hours are crucial for achieving sub‑day precision.

Overall, the study shows that a minimal analytic description—blackbody shock breakout plus a t² expansion—can recover the explosion time of core‑collapse SNe to better than a few hours for Type Ib/c events and to ≈0.2 days for a Type II‑P event, even when the first photometric point arrives >1 day after the true explosion. This level of precision dramatically narrows the temporal window for neutrino–optical coincidences, reducing the probability of accidental matches and enhancing the statistical significance of any true multimessenger detection.

Limitations include the neglect of detailed radiative‑transfer effects, possible contributions from ⁵⁶Ni decay at later times, and the reliance on a single‑instrument photometric set to avoid cross‑calibration issues. Future work should incorporate multi‑wavelength (UV/X‑ray) data, more sophisticated radiation‑hydrodynamics models, and a broader sample of SNe to quantify systematic uncertainties across different progenitor types.

In conclusion, the paper provides a practical, validated tool for the IceCube community and other multimessenger facilities to estimate SN explosion times from early optical light curves, thereby enabling tighter temporal correlations with high‑energy neutrino alerts and improving the prospects for discovering neutrino emission from core‑collapse supernovae.