Bolometric Light Curve Research Articles

Optical and Radio Analysis of Systematically Classified Broad-lined Type Ic Supernovae from the Zwicky Transient Facility

We study a magnitude-limited sample of 36 broad-lined type Ic supernovae (SNe Ic-BL) from the Zwicky Transient Facility Bright Transient Survey (detected between 2018 March and 2021 August), which is the largest systematic study of SNe Ic-BL done in literature thus far. We present the light curves (LCs) for each of the SNe and analyze the shape of the LCs to derive empirical parameters, along with the explosion epochs for every event. The sample has an average absolute peak magnitude in the r band of M¯r,max=−18.51±0.15 mag. Using spectra obtained around peak light, we compute expansion velocities from the Fe ii 5169 Å line for each event with high enough signal-to-noise ratio spectra, and find an average value of vph¯=16,100±1100 km s−1. We also compute bolometric LCs, study the blackbody temperature and radii evolution over time, and derive the explosion properties of the SNe. The explosion properties of the sample have average values of M¯Ni=0.37−0.06+0.08M⊙ , M¯ej=2.45−0.41+0.47M⊙ , and E¯K=(4.02−1.00+1.37)×1051 erg. Thirteen events have radio observations from the Very Large Array, with eight detections and five non-detections. We find that the populations that have radio detections and radio non-detections are indistinct from one another with respect to their optically inferred explosion properties, and there are no statistically significant correlations present between the events’ radio luminosities and optically inferred explosion properties. This provides evidence that the explosion properties derived from optical data alone cannot give inferences about the radio properties of SNe Ic-BL and likely their relativistic jet formation mechanisms.

Evolution of Tidal Disruption Event Disks with Magnetically Driven Winds

We present a time-dependent, one-dimensional, magnetically driven disk wind model based on magnetohydrodynamic (MHD) equations, in the context of tidal disruption events (TDEs). We assume that the disk is geometrically thin, is gas pressure dominated, and explicitly accounts for magnetic braking and turbulent viscosity through an extended α-viscosity prescription. We find a particular wind solution for a set of basic equations that satisfies the necessary and sufficient conditions for vertically unbound MHD flows. The solution shows that the disk evolves with mass loss due to wind and accretion from the initial Gaussian density distribution. We confirm that the mass accretion rate follows the power law of time t −19/16 at late times in the absence of wind, which matches the classical solution of J. K. Cannizzo et al. We find that the mass accretion rate is steeper than the t −19/16 curve when the wind is present. Mass accretion is also induced by magnetic braking, known as the wind-driven accretion mechanism, which results in a faster decay with time of both the mass accretion and mass-loss rates. In the disk emission, the ultraviolet (UV) luminosity is the highest among the optical, UV, and X-ray luminosities. While the optical and X-ray emission is observationally insignificant without magnetic braking, the X-ray emission is brighter at late times, especially in the presence of magnetic braking. This provides a possible explanation for observed delayed X-ray flares. Our model predicts that late-time bolometric light curves steeper than t −19/16 in UV-bright TDEs are potentially compelling indicators of magnetically driven winds.

The Primary Flare Following a Stellar Collision in a Galactic Nucleus

High-velocity stellar collisions near supermassive black holes may result in a complete disruption of the stars. The initial disruption can have energies on par with supernovae and power a very fast transient. In this work, we examine the primary flare that follows the initial transient, which arises when streams of gas from the disrupted stars travel around the central black hole and collide with each other on the antipodal side with respect to the original collision. We present a simple analytic estimate for the properties of the flare, which depends on the distance of the collision from the central black hole and on the center of mass velocity of the colliding stars. We also present the first-of-their-kind radiation-hydrodynamics simulations of a few examples of stellar collisions and postcollision flow of the ejected gas and calculate the expected bolometric light curves. We find that such postcollision flares are expected to be similar to flares that arise in tidal disruptions events of single stars.

Accretion onto a Supermassive Black Hole Binary before Merger

While supermassive binary black holes (SMBBHs) inspiral toward merger they may also accrete matter from a surrounding disk. To study the dynamics of this system requires simultaneously describing the evolving spacetime and the magnetized plasma. We present the first relativistic calculation simulating two equal-mass, nonspinning black holes as they inspiral from a 20 M (G = c = 1) initial separation almost to merger. Our results imply important observational consequences: for instance, the accretion rate Ṁ onto the black holes first decreases and then plateaus, dropping by only a factor of ∼3 despite the rapid inspiral. An estimated bolometric light curve follows the same profile, suggesting some merging SMBBHs may be significantly luminous past the predicted circumbinary disk decoupling. The minidisks are nonstandard: Reynolds, not Maxwell, stresses dominate, and they oscillate between two states. In one part of the cycle, “sloshing” streams transfer mass between minidisks, carrying kinetic energy at a rate sometimes as high as the peak minidisk bolometric luminosity. We also discover that episodic accretion drives time-varying minidisk tilts. These complex dynamics all contribute to unique cyclical behavior in the light curves of late-time inspiraling SMBBHs. The poloidal magnetic flux on the black holes is roughly constant at a dimensionless level ϕ ∼ 2–3, but doubles just before merger; for significant black hole spin, this flux predicts powerful jets with variability driven by binary dynamics, another potentially unique electromagnetic signature. This simulation is the first to employ our multipatch infrastructure PatchworkMHD, decreasing the computational expense to ∼3% of conventional single-grid methods’ cost.

SN 2021dbg: A Luminous Type IIP–IIL Supernova Exploding from a Massive Star with a Layered Shell

We present extensive observations and analysis of supernova (SN) SN 2021dbg, utilizing optical photometry and spectroscopy. For approximately 385 days following the explosion, SN 2021dbg exhibited remarkable luminosity, surpassing most Type II SNe (SNe II). This initial high luminosity is potentially attributed to interaction between the ejected material and the surrounding circumstellar material (CSM), as evidenced by the pronounced interaction signatures observed in its spectra. The subsequent high luminosity is primarily due to the significant 56Ni mass (0.17 ± 0.05 M ⊙) produced in the explosion. Based on the flux of flash emission lines detected in the initial spectra, we estimate that the CSM mass near the progenitor amounted to ∼(1.0–2.0) × 10−3 M ⊙, likely resulting from intense stellar wind activity 2–3 yr preceding the explosion. Considering the bolometric light curve, nebular spectrum modeling, and mass-loss rate, we suggest that the progenitor of SN 2021dbg was a red supergiant (RSG) with a mass of ∼20 M ⊙ and a radius of 1200 R ⊙. This RSG featured a thick hydrogen shell, which may have contained a region with a sharp decrease in material density, electron density, and temperature, contributing to its layered structure. This object demonstrates mixed features of Type IIP and IIL SNe, making it a transitional event linking the above two subclasses of SNe II.

The story of SN 2021aatd: A peculiar 1987A-like supernova with an early-phase luminosity excess

Context. There is a growing number of peculiar events that cannot be assigned to any of the main classes. SN 1987A and a handful of similar objects, thought to be explosive outcomes of blue supergiant stars, is one of them: while their spectra closely resemble those of H-rich (IIP) SNe, their light curve (LC) evolution is very different. Aims. Here we present the detailed photometric and spectroscopic analysis of SN 2021aatd, a peculiar Type II explosion. While its early-time evolution resembles that of the slowly evolving double-peaked SN 2020faa (although at a lower luminosity scale), after ∼40 days its LC shape becomes similar to that of SN 1987A-like explosions. Methods. In addition to comparing LCs, color curves, and spectra of SN 2021aatd to those of SNe 2020faa, 1987A, and other objects, we compared the observed spectra with our own SYN++ models and with the outputs of published radiative transfer models. We also carried out a detailed modeling of the pseudo-bolometric LCs of SNe 2021aatd and 1987A with a self-developed semi-analytical code, assuming a two-component ejecta (core + shell), and involving the rotational energy of a newborn magnetar in addition to radioactive decay. Results. We find that the photometric and the spectroscopic evolution of SN 2021aatd can be well described with the explosion of a ∼15 M⊙ blue supergiant star. Nevertheless, SN 2021aatd shows higher temperatures and weaker Na I D and Ba II 6142 Å lines than SN 1987A, which is instead reminiscent of IIP-like atmospheres. With the applied two-component ejecta model (accounting for decay and magnetar energy), we can successfully describe the bolometric LC of SN 2021aatd, including the first ∼40-day phase showing an excess compared to 87A-like SNe, but being strikingly similar to that of the long-lived SN 2020faa. Nevertheless, finding a unified model that also explains the LCs of more luminous events (e.g., SN 2020faa) is still a matter of debate.

STELLA Lightcurves of Energetic Pair-instability Supernovae in the Context of SN2018ibb

SN2018ibb is a recently observed hydrogen-poor superluminous supernova that appears to be powered by the decay of 30 M ⊙ of radioactive nickel. This supernova has been suggested to show hybrid signatures of a pair-instability supernova and an interacting supernova. In a previous paper, we found that rotating, metal-enriched pair-instability supernova progenitors appeared to check both of these boxes. In this paper, we model the lightcurves of the pair-instability supernovae using STELLA. We find that the STELLA models can explain the overall shape of the bolometric lightcurve of SN2018ibb, though not specific morphological features such as the luminosity peak or the bump at roughly 300 days after the peak. We also estimate the contribution from interaction and find that with relatively low wind velocities, the circumstellar medium originating from the stellar winds is consistent with the evidence for interaction in the spectra. The observed values of the photosphere velocity in the 100 days after peak luminosity are similar to the STELLA models, but the deceleration is lower. This leads to the biggest inconsistency, which is the blackbody temperature of SN2018ibb being much hotter than any of the STELLA models. We note that this high temperature (and the flat velocity) may be difficult to reconcile with the long rise time of SN2018ibb, but nevertheless conclude that if it is accurate, this discrepancy represents a challenge for SN2018ibb being a robust PISN candidate. This result is noteworthy given the lack of other scenarios for this supernova.

SN 2019tua: A Type IIb Supernova with Multiple Bumps in the Light Curves

We present photometric and spectroscopic observations and analysis of the Type IIb supernova (SN) SN 2019tua, which exhibits multiple bumps in its declining light curves between 40 and 65 days after discovery. SN 2019tua shows a time to peak of about 25 days, similar to other Type IIb SNe. Our observations indicate a decrease in its brightness of about 1 mag in the 60 days after the peak. At about days 50 and 60, its multiband light curves exhibit bumpy behavior. The complex luminosity evolution of SN 2019tua could not be well modeled with a single currently popular energy source model, e.g., radioactive decay of 56Ni, magnetar, interaction between the ejecta and a circumstellar shell. Even though the magnetar model has a smaller χ 2/dof value, the complex changes in SN 2019tua’s brightness suggest that more than one physical process might be involved. We propose a hybrid circumstellar medium (CSM) interaction plus 56Ni model to explain the bolometric light curve of SN 2019tua. The fitting results show that the ejecta mass M ej ≈ 2.4 M ⊙, the total CSM mass M CSM ≈ 1.0 M ⊙, and the 56Ni mass M Ni ≈ 0.4 M ⊙. The total kinetic energy of the ejecta is E k ≈ 0.5 × 1051 erg. Preexisting multiple shells suggest that the progenitor of SN 2019tua experienced mass ejections within approximately ∼6–44 yr prior to the explosion.

Progenitor mass and ejecta asymmetry of supernova 2023ixf from nebular spectroscopy

Context. Supernova (SN) 2023ixf was discovered in the galaxy M 101 in May 2023. Its proximity provided the scientific community an extremely valuable opportunity to study the characteristics of the SN and its progenitor. A point source detected on archival images and hydrodynamical modeling of the bolometric light curve have been used to constrain the former star’s properties. There is a significant variation in the published results regarding the initial mass of the progenitor. Nebular spectroscopy can be used to enhance our understanding of the SN and its progenitor. Aims. We determined the SN progenitor mass by studying the first published nebular spectrum, taken 259 days after the explosion. Methods. We analyzed the nebular spectrum taken with GMOS at the Gemini North Telescope. We identified typical emission lines, such as [O I], Hα, and [Ca II], among others. Some species’ line profiles show broad and narrow components, indicating two ejecta velocities and an asymmetric ejecta. We inferred the progenitor mass of SN 2023ixf by comparing its spectra with synthetic spectra and by measuring the forbidden oxygen doublet flux. Results. Based on the flux ratio and the direct comparison with spectra models, the progenitor star of SN 2023ixf had a MZAMS between 12 and 15 M⊙. We find that using the [O I] doublet flux provides a less tight constraint on the progenitor mass. Our results agree with those from hydrodynamical modeling of the early light curve and pre-explosion image estimates that point to a relatively low-mass progenitor.

Light-curve Structure and Hα Line Formation in the Tidal Disruption Event AT 2019azh

AT 2019azh is a H+He tidal disruption event (TDE) with one of the most extensive ultraviolet and optical data sets available to date. We present our photometric and spectroscopic observations of this event starting several weeks before and out to approximately 2 yr after the g-band's peak brightness and combine them with public photometric data. This extensive data set robustly reveals a change in the light-curve slope and a possible bump in the rising light curve of a TDE for the first time, which may indicate more than one dominant emission mechanism contributing to the pre-peak light curve. Indeed, we find that the MOSFiT-derived parameters of AT 2019azh, which assume reprocessed accretion as the sole source of emission, are not entirely self-consistent. We further confirm the relation seen in previous TDEs whereby the redder emission peaks later than the bluer emission. The post-peak bolometric light curve of AT 2019azh is better described by an exponential decline than by the canonical t −5/3 (and in fact any) power-law decline. We find a possible mid-infrared excess around the peak optical luminosity, but cannot determine its origin. In addition, we provide the earliest measurements of the Hα emission-line evolution and find no significant time delay between the peak of the V-band light curve and that of the Hα luminosity. These results can be used to constrain future models of TDE line formation and emission mechanisms in general. More pre-peak 1–2 days cadence observations of TDEs are required to determine whether the characteristics observed here are common among TDEs. More importantly, detailed emission models are needed to fully exploit such observations for understanding the emission physics of TDEs.

Magnetars as powering sources of gamma-ray burst associated supernovae, and unsupervized clustering of cosmic explosions

ABSTRACT We present the semi-analytical light curve modelling of 13 supernovae associated with gamma-ray bursts (GRB-SNe) along with two relativistic broad-lined (Ic-BL) SNe without GRB association (SNe 2009bb and 2012ap), considering millisecond magnetars as central-engine-based power sources for these events. The bolometric light curves of all 15 SNe in our sample are well-regenerated utilizing a χ2-minimization code, MINIM, and numerous parameters are constrained. The median values of ejecta mass (Mej), magnetar’s initial spin period (Pi), and magnetic field (B) for GRB-SNe are determined to be ≈5.2 M⊙, 20.5 ms, and 20.1 × 1014 G, respectively. We leverage machine learning (ML) algorithms to comprehensively compare the three-dimensional parameter space encompassing Mej, Pi, and B for GRB-SNe determined herein to those of H-deficient superluminous SNe (SLSNe-I), fast blue optical transients (FBOTs), long GRBs (LGRBs), and short GRBs (SGRBs) obtained from the literature. The application of unsupervized ML clustering algorithms on the parameters Mej, Pi, and B for GRB-SNe, SLSNe-I, and FBOTs yields a classification accuracy of ∼95 per cent. Extending these methods to classify GRB-SNe, SLSNe-I, LGRBs, and SGRBs based on Pi and B values results in an accuracy of ∼84 per cent. Our investigations show that GRB-SNe and relativistic Ic-BL SNe presented in this study occupy different parameter spaces for Mej, Pi, and B than those of SLSNe-I, FBOTs, LGRBs, and SGRBs. This indicates that magnetars with different Pi and B can give birth to distinct types of transients.

Diversity of Early Kilonova with the Realistic Opacities of Highly Ionized Heavy Elements

We investigate early (t < 1 day) kilonova from a neutron star merger by deriving atomic opacities for all the elements from La to Ra (Z = 57–88) ionized to the states V–XI. The opacities at high temperatures for the elements with open f-shells (e.g., lanthanides) are exceptionally high, reaching κexp∼3×103cm2g−1 at λ ≤ 1000 Å at T ∼ 70,000 K, whereas the opacities at the same temperature and wavelengths for the elements with open d-, p-, and s-shells reach κexp∼1 , 0.1, and 0.01 cm2 g−1, respectively. Using the new opacity data set, we derive early kilonovae for various compositions and density structures expected for neutron star merger ejecta. The bolometric luminosity of the lanthanide-rich ejecta shows distinct signatures and is fainter than that of the lanthanide-free ejecta. Early luminosity is suppressed by the presence of a thin outer layer, agreeing with the results of Kasen et al. and Banerjee et al. The early brightness in the Swift UVOT filters and in the optical g, r, i, and z filters for a source at 100 Mpc are about ∼22–19.5 and ∼21–20 mag, respectively, at t ∼ 0.1 day. Such kilonovae are ideal targets for the upcoming UV satellites, such as ULTRASAT, UVEX, and DORADO, and the upcoming surveys, e.g., the Vera Rubin Observatory. We suggest that the gray opacities that reproduce the bolometric light curves with and without lanthanides are ∼1–10 and ∼0.8 cm2 g−1.

An analytic description of electron thermalization in kilonovae ejecta

ABSTRACT A simple analytic description is provided of the rate of energy deposition by β-decay electrons in the homologously expanding radioactive plasma ejected in neutron star mergers, valid for a wide range of ejecta parameters – initial entropy, electron fraction {s0, Ye}, and scaled (time-independent) density ρt3. The formulae are derived using detailed numerical calculations following the time-dependent composition and β-decay emission spectra (including the effect of delayed deposition). The deposition efficiency depends mainly on ρt3 and only weakly on {s0, Ye}. The time te at which the ratio between the rates of electron energy deposition and energy production drops to 1 − e−1, is given by $t_e=t_{0e}\Big (\frac{\rho t^3}{0.5(\rho t^3)_0}\Big)^a$, where $(\rho t^3)_0=\frac{0.05\, {\rm M}_{\odot }}{4\pi (0.2c)^3}$, t0e(s0, Ye) ≈ 17 d, and 0.4 ≤ a(s0, Ye) ≤ 0.5. The fractional uncertainty in te due to nuclear physics uncertainties is $\approx 10~{{\ \rm per\ cent}}$. The result a ≤ 0.5 reflects the fact that the characteristic β-decay electron energies do not decrease with time (largely due to ‘inverted decay chains’ in which a slowly decaying isotope decays to a rapidly decaying isotope with higher end-point energy). We provide an analytic approximation for the time-dependent electron energy deposition rate, reproducing the numerical results to better than 50 per cent (typically $\lt 30~{{\ \rm per\ cent}}$, well within the energy production rate uncertainty due to nuclear physics uncertainties) over a 3–4 orders-of-magnitude deposition rate decrease with time. Our results may be easily incorporated in calculations of kilonovae light curves (with general density and composition structures), eliminating the need to numerically follow the time-dependent electron spectra. Identifying te, e.g. in the bolometric light curve, will constrain the ejecta density distribution.

The complex circumstellar environment of supernova 2023ixf.

The early evolution of a supernova (SN) can reveal information about the environment and the progenitor star. When a star explodes in vacuum, the first photons to escape from its surface appear as a brief, hours-long shock-breakout flare1,2, followed by a cooling phase of emission. However, for stars exploding within a distribution of dense, optically thick circumstellar material (CSM), the first photons escape from the material beyond the stellar edge and the duration of the initial flare can extend to several days, during which the escaping emission indicates photospheric heating3. Early serendipitous observations2,4 that lacked ultraviolet (UV) data were unable to determine whether the early emission is heating or cooling and hence the nature of the early explosion event. Here we report UV spectra of the nearby SN 2023ixf in the galaxy Messier 101 (M101). Using the UV data as well as a comprehensive set of further multiwavelength observations, we temporally resolve the emergence of the explosion shock from a thick medium heated by the SN emission. We derive a reliable bolometric light curve that indicates that the shock breaks out from a dense layer with a radius substantially larger than typical supergiants.

Abundance stratification in type Ia supernovae – VII. The peculiar, C-rich iPTF16abc: highlighting diversity among luminous events

ABSTRACT Observations of Type Ia supernovae (SNe Ia) reveal diversity, even within assumed subcategories. Here, the composition of the peculiar iPTF16abc (SN 2016bln) is derived by modelling a time series of optical spectra. iPTF16abc’s early spectra combine traits of SNe 1999aa and 1991T known for weak Si ii λ 6355 and prominent Fe iii features. However, it differs with weak early Fe iii lines, and persistent C ii lines post-peak. It also exhibits a weak Ca ii H&amp;K feature aligning it with SN 1991T, an observation supported by their bolometric light curves. The early attenuation of Fe iii results from abundance effect. The weakening of the Si ii λ 6355 line, stems from silicon depletion in the outer shells, a characteristic shared by both SNe 1999aa and 1991T, indicating a common explosion mechanism that terminates nuclear burning at around 12 000 km s−1 unseen in normal events. Beneath a thin layer of intermediate mass elements (IMEs) with a total mass of 0.18 M⊙, extends a 56Ni rich shell totaling 0.76 M⊙ and generating a bolometric luminosity as high as Lpeak = 1.60 ± 0.1 × 1043 ergs s−1. Inner layers, typical of SNe Ia, hold neutron-rich elements, (54Fe and 58Ni), totaling 0.20 M⊙. Stable iron, exceeding solar abundance, and carbon, coexist in the outermost layers, challenging existing explosion models. The presence of carbon down to v ≈ 9000 km s−1, totalling ∼0.01 M⊙ unprecedented in this class, links iPTF16abc to SN 2003fg-like events. The retention of 91T-like traits in iPTF16abc underscores its importance in understanding the diversity of SNe Ia.

The carbon-rich type Ic supernova 2016adj in the iconic dust lane of Centaurus A: Potential signatures of an interaction with circumstellar hydrogen

We present a comprehensive data set of supernova (SN) 2016adj located within the central dust lane of Centaurus A. SN 2016adj is significantly reddened and after correcting the peak apparent $B$-band magnitude ($m_B = 17.48 for Milky Way reddening and our inferred host-galaxy reddening parameters (i.e. V host and V host we estimated it reached a peak absolute magnitude of $M_B -18$. A detailed inspection of the optical and near-infrared (NIR) spectroscopic time series reveals a carbon-rich SN Ic and not a SN Ib/IIb as previously suggested in the literature. The NIR spectra show prevalent carbon-monoxide formation occurring already by $+$41 days past $B$-band maximum, which is $ 11$ days earlier than previously reported in the literature for this object. Interestingly, around two months past maximum, the NIR spectrum of SN 2016adj begins to exhibit H features, with a $+$97 days medium resolution spectrum revealing both Paschen and Bracket lines with absorption minima of $ 2,000$ km $, full-width-half-maximum emission velocities of $ 1,000$ km $, and emission line ratios consistent with a dense emission region. We speculate that these attributes are due to a circumstellar interaction (CSI) between the rapidly expanding SN ejecta and a H-rich shell of material that formed during the pre-SN phase. A bolometric light curve was constructed and a semi-analytical model fit suggests the SN synthesized 0.5 $M_ sun $ of 56Ni and ejected 4.7 $M_ sun $ of material, though these values should be approached with caution given the large uncertainties associated with the adopted reddening parameters and known light echo emission. Finally, inspection of the Hubble Space Telescope archival data yielded no progenitor detection.

Circumstellar interaction models for the early bolometric light curve of SN 2023ixf

Type II supernovae (SNe II) show growing evidence of an interaction with circumstellar material (CSM) surrounding their progenitor stars as a consequence of enhanced mass loss during the last years of the progenitor’s life, although the exact mechanism is still unknown. We present an analysis of the progenitor mass-loss history of SN 2023ixf, a nearby SN II showing signs of an interaction. First, we calculated the early-time (&lt; 19 days) bolometric light curve for SN 2023ixf based on the integration of the observed flux covering ultraviolet, optical and near-infrared bands, and black-body extrapolations for the unobserved flux. Our calculations detected the sudden increase to maximum luminosity and temperature, in addition to the subsequent fall, displaying an evident peak. This is the first time that this phase can be precisely estimated for a SN II. We used the early-time bolometric light curve of SN 2023ixf to test the calibrations of bolometric corrections against colours from the literature. In addition, we included the observations of SN 2023ixf into some of the available calibrations to extend their use to earlier epochs. A comparison of the observed bolometric light curve to SN II explosion models with CSM interaction suggests a progenitor mass-loss rate of Ṁ = 3 × 10−3 M⊙ yr−1 confined to 12 000 R⊙ (∼8 × 1014 cm) and a wind acceleration parameter of β = 5. This model reproduces the early bolometric light curve, expansion velocities, and the epoch of disappearance of interacting lines in the spectra. This model indicates that the wind was launched ∼80 yr before the explosion. If the effect of the wind acceleration is not taken into account, the enhanced wind must have developed over the final months to years prior to the SN, which may not be consistent with the lack of outburst detection in pre-explosion images over the last ∼20 yr before explosion.

RAINBOW: A colorful approach to multipassband light-curve estimation

Context. Time series generated by repeatedly observing astronomical transients are generally sparse, irregularly sampled, noisy, and multidimensional (obtained through a set of broad-band filters). In order to fully exploit their scientific potential, it is necessary to use this incomplete information to estimate a continuous light-curve behavior. Traditional approaches use ad hoc functional forms to approximate the light curve in each filter independently (hereafter, the MONOCHROMATIC method). Aims. We present RAINBOW, a physically motivated framework that enables simultaneous multiband light-curve fitting. It allows the user to construct a 2D continuous surface across wavelength and time, even when the number of observations in each filter is significantly limited. Methods. Assuming the electromagnetic radiation emission from the transient can be approximated by a blackbody, we combined an expected temperature evolution and a parametric function describing its bolometric light curve. These three ingredients allow the information available in one passband to guide the reconstruction in the others, thus enabling a proper use of multisurvey data. We demonstrate the effectiveness of our method by applying it to simulated data from the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC) as well as to real data from the Young Supernova Experiment (YSE DR1). Results. We evaluate the quality of the estimated light curves according to three different tests: goodness of fit, peak-time prediction, and ability to transfer information to machine-learning (ML) based classifiers. The results confirm that RAINBOW leads to an equivalent goodness of fit (supernovae II) or to a goodness of fit that is better by up to 75% (supernovae Ibc) than the MONOCHROMATIC approach. Similarly, the accuracy improves for all classes in our sample when the RAINBOW best-fit values are used as a parameter space in a multiclass ML classification. Conclusions. Our approach enables a straightforward light-curve estimation for objects with observations in multiple filters and from multiple experiments. It is particularly well suited when the light-curve sampling is sparse. We demonstrate its potential for characterizing supernova-like events here, but the same approach can be used for other classes by changing the function describing the light-curve behavior and temperature representation. In the context of the upcoming large-scale sky surveys and their potential for multisurvey analysis, this represents an important milestone in the path to enable population studies of photometric transients.

1100 days in the life of the supernova 2018ibb

Stars with zero-age main sequence masses between 140 and 260 M⊙ are thought to explode as pair-instability supernovae (PISNe). During their thermonuclear runaway, PISNe can produce up to several tens of solar masses of radioactive nickel, resulting in luminous transients similar to some superluminous supernovae (SLSNe). Yet, no unambiguous PISN has been discovered so far. SN 2018ibb is a hydrogen-poor SLSN at z = 0.166 that evolves extremely slowly compared to the hundreds of known SLSNe. Between mid 2018 and early 2022, we monitored its photometric and spectroscopic evolution from the UV to the near-infrared (NIR) with 2–10 m class telescopes. SN 2018ibb radiated &gt; 3 × 1051 erg during its evolution, and its bolometric light curve reached &gt; 2 × 1044 erg s−1 at its peak. The long-lasting rise of &gt; 93 rest-frame days implies a long diffusion time, which requires a very high total ejected mass. The PISN mechanism naturally provides both the energy source (56Ni) and the long diffusion time. Theoretical models of PISNe make clear predictions as to their photometric and spectroscopic properties. SN 2018ibb complies with most tests on the light curves, nebular spectra and host galaxy, and potentially all tests with the interpretation we propose. Both the light curve and the spectra require 25–44 M⊙ of freshly nucleosynthesised 56Ni, pointing to the explosion of a metal-poor star with a helium core mass of 120–130 M⊙ at the time of death. This interpretation is also supported by the tentative detection of [Co II] λ 1.025 μm, which has never been observed in any other PISN candidate or SLSN before. We observe a significant excess in the blue part of the optical spectrum during the nebular phase, which is in tension with predictions of existing PISN models. However, we have compelling observational evidence for an eruptive mass-loss episode of the progenitor of SN 2018ibb shortly before the explosion, and our dataset reveals that the interaction of the SN ejecta with this oxygen-rich circumstellar material contributed to the observed emission. That may explain this specific discrepancy with PISN models. Powering by a central engine, such as a magnetar or a black hole, can be excluded with high confidence. This makes SN 2018ibb by far the best candidate for being a PISN, to date.

Extrabol: A Python Package for Estimating Bolometric Light Curves of Thermal Transients

We introduce a new, open-source, Python-based package, extrabol, for inferring the bolometric light curve evolution of extragalactic thermal transients. extrabol uses non-parametric Gaussian Process regression for light curve estimation that requires minimal user interaction. extrabol is available via GitHub.