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 practical problem faced by the IceCube neutrino observatory: how to establish a precise temporal coincidence between a neutrino burst and the optical signature of a core‑collapse supernova (CCSN). While IceCube already runs an online neutrino‑burst alert system that forwards candidate directions to a network of optical telescopes, the identification of a supernova in the follow‑up images still leaves the explosion time ambiguous. Traditionally, the explosion time is approximated by the time of the first optical detection, which can be off by days to weeks, severely limiting the statistical significance of any neutrino‑photon coincidence.

To overcome this limitation, the authors propose a simple, physically motivated model that fits the early optical light curve of a CCSN and directly yields the explosion epoch (t₀). The model consists of two phases. The first is the shock‑breakout (or “shock‑cooling”) phase, during which the emergent luminosity rises roughly as a power law L ∝ (t‑t₀)ⁿ with n≈2, reflecting the diffusion of the shock through the stellar envelope. The second phase is the cooling envelope phase, where the luminosity decays as a power law L ∝ (t‑t₀)⁻α, with α typically between 0.5 and 0.8 depending on the progenitor’s radius, opacity, and ejecta mass. The two regimes are joined smoothly at a transition time t₁, and the full functional form is fitted to multi‑band photometry (UBVRI and, when available, X‑ray) after correcting for distance modulus and Galactic extinction.

The fitting procedure employs a non‑linear least‑squares optimizer (Levenberg‑Marquardt) to solve for the parameters A, B, n, α, t₁, and most importantly t₀. To assess uncertainties, the authors run Markov Chain Monte Carlo (MCMC) simulations around the best‑fit solution, extracting credible intervals for t₀.

Three well‑studied supernovae are used as test cases: SN 1987A (the prototypical nearby event in the Large Magellanic Cloud), SN 2006aj (associated with GRB 060218), and SN 2008D (discovered via an X‑ray flash). For each, extensive photometric data spanning the first few days after explosion are compiled. The model reproduces the observed light curves with reduced χ² values close to unity, and the inferred explosion times are:

• SN 1987A: t₀ = 1987‑02‑23 07:35 UT ± 0.3 day (≈ 7 h)
• SN 2006aj: t₀ = 2006‑02‑18 03:21 UT ± 0.12 day (≈ 3 h)
• SN 2008D: t₀ = 2008‑01‑09 13:07 UT ± 0.09 day (≈ 2 h)

These uncertainties are an order of magnitude smaller than those obtained by simply adopting the first detection time. Monte‑Carlo experiments further demonstrate that a cadence of ≤ 0.5 day (i.e., at least one observation per 12 hours) is sufficient to keep the t₀ error below 5 hours, provided the photometric errors are ≤ 0.05 mag.

The authors discuss systematic limitations. The model assumes spherical symmetry and a constant opacity, which may not hold for highly aspherical explosions or for progenitors surrounded by dense circum‑stellar material (CSM). In such cases, the rise index n and decay index α can deviate from the canonical values, potentially biasing t₀. Nevertheless, the authors argue that for the majority of CCSNe—especially those without strong CSM interaction—the model remains robust.

From the perspective of IceCube, the ability to constrain the explosion time to a few hours dramatically narrows the temporal search window for associated neutrinos. Since the background neutrino rate scales linearly with the size of the time window, reducing the window from days to hours can improve the signal‑to‑noise ratio by a factor of ~10–20. Consequently, the probability of a statistically significant detection of a CCSN neutrino burst—if such a burst occurs within IceCube’s field of view—is substantially increased.

In conclusion, the paper presents a practical, data‑driven method to extract CCSN explosion times from early optical (and optionally X‑ray) light curves. Validation on three historic events shows that the method can achieve sub‑day, often sub‑hour, precision. The approach is directly applicable to real‑time multi‑messenger campaigns, providing IceCube and other neutrino observatories with a powerful tool to tighten photon‑neutrino temporal coincidences and enhance the discovery potential for astrophysical neutrino sources.