Supernovae Ic-BL Research Articles

Gravitational-wave Radiation from the Magnetar-driven Supernovae

Rapidly spinning magnetars are potential candidates for the energy source of supernovae (SNe) and gamma-ray bursts and the most promising sources for continuous gravitational waves (GWs) detected by ground-based GW detectors. Continuous GWs can be radiated from magnetars due to magnetic-induced deformation or fluid oscillations, compatible with magnetic dipole (MD) radiation for spin-down energy. In this paper, we investigate the diverse light curves of magnetar-driven SNe in the scenario that the spin-down is dominated by GW radiation and/or MD radiation. By simulating the light curves of SNe and employing the Markov Chain Monte Carlo method, we constrain the parameters of the magnetars and SN explosions and show that the signature of GW radiation may be indicated by the bolometric luminosity curves of SNe Ic-BL 2007ru and 2009bb. We find that the ellipticity of magnetars in the order of 10−3 can be induced by the magnetic field of ∼1016 G. If such continuous GWs associated with SNe can be detected in the future by the Advanced LIGO and Virgo detectors, this would be a smoking gun for a magnetar engine powering SNe.

Collapsars as Sites of r-process Nucleosynthesis: Systematic Photometric Near-infrared Follow-up of Type Ic-BL Supernovae

One of the open questions following the discovery of GW170817 is whether neutron star (NS) mergers are the only astrophysical sites capable of producing r-process elements. Simulations have shown that 0.01–0.1 M ⊙ of r-process material could be generated in the outflows originating from the accretion disk surrounding the rapidly rotating black hole that forms as a remnant to both NS mergers and collapsing massive stars associated with long-duration gamma-ray bursts (collapsars). The hallmark signature of r-process nucleosynthesis in the binary NS merger GW170817 was its long-lasting near-infrared (NIR) emission, thus motivating a systematic photometric study of the light curves of broad-lined stripped-envelope (Ic-BL) supernovae (SNe) associated with collapsars. We present the first systematic study of 25 SNe Ic-BL—including 18 observed with the Zwicky Transient Facility and 7 from the literature—in the optical/NIR bands to determine what quantity of r-process material, if any, is synthesized in these explosions. Using semi-analytic models designed to account for r-process production in SNe Ic-BL, we perform light curve fitting to derive constraints on the r-process mass for these SNe. We also perform independent light curve fits to models without the r-process. We find that the r-process-free models are a better fit to the light curves of the objects in our sample. Thus, we find no compelling evidence of r-process enrichment in any of our objects. Further high-cadence infrared photometric studies and nebular spectroscopic analysis would be sensitive to smaller quantities of r-process ejecta mass or indicate whether all collapsars are completely devoid of r-process nucleosynthesis.

Diagnosing the ejecta properties of engine-driven supernovae from observables in their initial phase

ABSTRACT Engine-driven explosions with continuous energy input from the central system have been suggested for supernovae (SNe) associated with a Gamma-Ray Burst (GRB), superluminous SNe (SLSNe), and at least a fraction of broad-lined SNe Ic (SNe Ic-BL) even without an associated GRB. In the present work, we investigate observational consequences in this scenario, focusing on the case where the energy injection is sufficiently brief, which has been suggested for GRB-SNe. We construct a simplified, spherical ejecta model sequence taking into account the major effects of the central engine; composition mixing, density structure, and the outermost ejecta velocity. Unlike most of the previous works for GRB-SNe, we solve the formation of the photosphere self-consistently, with which we can predict the photometric and spectroscopic observables. We find that these ejecta properties strongly affect their observational appearance in the initial phase (≲ a week since the explosion), highlighted by blended lines suffering from higher-velocity absorptions for the flatter density distribution and/or higher outermost ejeca velocity. This behaviour also affects the multiband light curves in a non-monotonic way. Prompt follow-up observations starting immediately after the explosion thus provides key diagnostics to unveil the nature of the central engine behind GRB-SNe and SNe Ic-BL. For SN 2017iuk associated with GRB 171205A these diagnosing observational data are available, and we show that the expected structure from the engine-driven explosion, i.e. a flat power-law density structure extending up to ≳100 000 km s−1, can explain the observed spectral evolution reasonably well.

Supernova Luminosity Powered by Magnetar–Disk System

Abstract Magnetars are one of the potential power sources for some energetic supernova explosions such as type I superluminous supernovae (SLSNe I) and broad-lined type Ic supernovae (SNe Ic-BL). In order to explore the possible link between these two subclasses of supernovae (SNe), we study the effect of fallback accretion disk on magnetar evolution and magnetar-powered SNe. In this scenario, the interaction between a magnetar and a fallback accretion disk would accelerate the spin of the magnetar in the accretion regime but could result in substantial spin-down of the magnetars in the propeller regime. Thus, the initial rotation of the magnetar plays a less significant role in the spin evolution. Such a magnetar–disk interaction scenario can explain well the light curves of both SNe Ic-BL and SLSNe I, for which the observed differences are sensitive to the initial magnetic field of the magnetar and the fallback mass and timescale for the disk. Compared to the magnetars powering the SNe Ic-BL, those accounting for more luminous SNe usually maintain faster rotation and have relatively lower effective magnetic fields around peak time. In addition, the association between SLSNe I and long gamma-ray bursts, if observed in the future, could be explained in the context of a magnetar–disk system.

ASASSN-18am/SN 2018gk: an overluminous Type IIb supernova from a massive progenitor

ABSTRACT ASASSN-18am/SN 2018gk is a newly discovered member of the rare group of luminous, hydrogen-rich supernovae (SNe) with a peak absolute magnitude of MV ≈ −20 mag that is in between normal core-collapse SNe and superluminous SNe. These SNe show no prominent spectroscopic signatures of ejecta interacting with circumstellar material (CSM), and their powering mechanism is debated. ASASSN-18am declines extremely rapidly for a Type II SN, with a photospheric-phase decline rate of ∼6.0 mag (100 d)−1. Owing to the weakening of H i and the appearance of He i in its later phases, ASASSN-18am is spectroscopically a Type IIb SN with a partially stripped envelope. However, its photometric and spectroscopic evolution shows significant differences from typical SNe IIb. Using a radiative diffusion model, we find that the light curve requires a high synthesized 56Ni mass $M_{\rm Ni} \sim 0.4\, \rm {M_{\odot }}$ and ejecta with high kinetic energy Ekin = (7–10) × 1051 erg. Introducing a magnetar central engine still requires $M_{\rm Ni} \sim 0.3\, \rm {M_{\odot }}$ and Ekin = 3 × 1051 erg. The high 56Ni mass is consistent with strong iron-group nebular lines in its spectra, which are also similar to several SNe Ic-BL with high 56Ni yields. The earliest spectrum shows ‘flash ionization’ features, from which we estimate a mass-loss rate of $\dot{M}\approx 2\times 10^{-4} \, \rm \rm {M_{\odot }}\,yr^{-1}$. This wind density is too low to power the luminous light curve by ejecta–CSM interaction. We measure expansion velocities as high as 17 000 $\rm {\, km\, s^{-1}}$ for Hα, which is remarkably high compared to other SNe II. We estimate an oxygen core mass of 1.8–3.4 M⊙ using the [O i] luminosity measured from a nebular-phase spectrum, implying a progenitor with a zero-age main-sequence mass of 19–26 M⊙.

Connection of supernovae 2002ap, 2003gd, 2013ej, and 2019krl in M 74 with atomic gas accretion and spiral structure

Studying the nature of various types of supernovae (SNe) is important for our understanding of stellar evolution. Observations of atomic and molecular gas in the host galaxies of gamma-ray bursts (GRBs) and SNe have recently been used to learn about the nature of the explosions themselves and the star formation events during which their progenitors were born. Based on archival data for M 74, which previously has not been investigated in the context of SN positions, we report the gas properties in the environment of the broad-lined type Ic (Ic-BL) SN 2002ap and the type II SNe 2003gd, 2013ej, and 2019krl. The SN 2002ap is located at the end of an off-centre, asymmetric, 55 kpc-long HIextension containing 7.5% of the total atomic gas in M 74, interpreted as a signature of external gas accretion. It is the fourth known case of an explosion of a presumably massive star located close to a concentration of atomic gas (after GRBs 980425, 060505, and SN 2009bb). It is unlikely that all these associations are random (at a 3σsignificance), so the case of SN 2002ap adds to the evidence that the birth of the progenitors of type Ic-BL SNe and GRBs is connected with the accretion of atomic gas from the intergalactic medium. The HIextension could come from tidally disrupted companions of M 74, or be a remnant of a galaxy or a gas cloud that accreted entirely from the intragroup medium. The other (type II) SNe in M 74 are located at the outside edge of a spiral arm. This suggests that either their progenitors were born when gas was piling up there or that the SN progenitors moved away from the arm due to their orbital motions. These type II SNe do not seem to be related to gas accretion.

Host Galaxies of Type Ic and Broad-lined Type Ic Supernovae from the Palomar Transient Factory: Implications for Jet Production

Abstract Unlike ordinary supernovae (SNe), some of which are hydrogen and helium deficient (called Type Ic SNe), broad-lined Type Ic SNe (SNe Ic-bl) are very energetic events, and only SNe Ic-bl are coincident with long-duration gamma-ray bursts (GRBs). Understanding the progenitors of SN Ic-bl explosions versus those of their SN Ic cousins is key to understanding the SN–GRB relationship and jet production in massive stars. Here we present the largest existing set of host galaxy spectra of 28 SNe Ic and 14 SNe Ic-bl, all discovered by the same galaxy-untargeted survey, namely, the Palomar Transient Factory (PTF). We carefully measure their gas-phase metallicities, stellar masses (M *), and star formation rates (SFRs). We further reanalyze the hosts of 10 literature SN–GRBs using the same methods and compare them to our PTF SN hosts with the goal of constraining their progenitors from their local environments. We find that the metallicities, SFRs, and M * values of our PTF SN Ic-bl hosts are statistically comparable to those of SN–GRBs but significantly lower than those of the PTF SNe Ic. The mass–metallicity relations as defined by the SNe Ic-bl and SN–GRBs are not significantly different from the same relations as defined by Sloan Digital Sky Survey galaxies, contradicting claims by earlier works. Our findings point toward low metallicity as a crucial ingredient for SN Ic-bl and SN–GRB production since we are able to break the degeneracy between high SFR and low metallicity. We suggest that the PTF SNe Ic-bl may have produced jets that were choked inside the star or were able to break out of the star as unseen low-luminosity or off-axis GRBs.

Broad-lined type Ic supernova iPTF16asu: A challenge to all popular models

ABSTRACT It is well known that ordinary supernovae (SNe) are powered by 56Ni cascade decay. Broad-lined type Ic SNe (SNe Ic-BL) are a subclass of SNe that are not all exclusively powered by 56Ni decay. It was suggested that some SNe Ic-BL are powered by magnetar spin-down. iPTF16asu is a peculiar broad-lined type Ic supernova discovered by the intermediate Palomar Transient Factory. With a rest-frame rise time of only 4 d, iPTF16asu challenges the existing popular models, for example, the radioactive heating (56Ni-only) and the magnetar +56Ni models. Here we show that this rapid rise could be attributed to interaction between the SN ejecta and a pre-existing circumstellar medium ejected by the progenitor during its final stages of evolution, while the late-time light curve can be better explained by energy input from a rapidly spinning magnetar. This model is a natural extension to the previous magnetar model. The mass-loss rate of the progenitor and ejecta mass are consistent with a progenitor that experienced a common envelope evolution in a binary. An alternative model for the early rapid rise of the light curve is the cooling of a shock propagating into an extended envelope of the progenitor. It is difficult at this stage to tell which model (interaction+magnetar + 56Ni or cooling+magnetar + 56Ni) is better for iPTF16asu. However, it is worth noting that the inferred envelope mass in the cooling+magnetar + 56Ni is very high.

The Carnegie Supernova Project I

Stripped-envelope (SE) supernovae (SNe) include H-poor (Type IIb), H-free (Type Ib), and He-free (Type Ic) events thought to be associated with the deaths of massive stars. The exact nature of their progenitors is a matter of debate with several lines of evidence pointing towards intermediate mass (Minit< 20 M⊙) stars in binary systems, while in other cases they may be linked to single massive Wolf-Rayet stars. Here we present the analysis of the light curves of 34 SE SNe published by the Carnegie Supernova Project (CSP-I) that are unparalleled in terms of photometric accuracy and wavelength range. Light-curve parameters are estimated through the fits of an analytical function and trends are searched for among the resulting fit parameters. Detailed inspection of the dataset suggests a tentative correlation between the peak absolute B-band magnitude and Δm15(B), while the post maximum light curves reveals a correlation between the late-time linear slope and Δm15. Making use of the full set of optical and near-IR photometry, combined with robust host-galaxy extinction corrections, comprehensive bolometric light curves are constructed and compared to both analytic and hydrodynamical models. This analysis finds consistent results among the two different modeling techniques and from the hydrodynamical models we obtained ejecta masses of 1.1–6.2M⊙, 56Ni masses of 0.03–0.35M⊙, and explosion energies (excluding two SNe Ic-BL) of 0.25–3.0 × 1051 erg. Our analysis indicates that adopting κ = 0.07 cm2 g-1 as the mean opacity serves to be a suitable assumption when comparing Arnett-model results to those obtained from hydrodynamical calculations. We also find that adopting He i and O i line velocities to infer the expansion velocity in He-rich and He-poor SNe, respectively, provides ejecta masses relatively similar to those obtained by using the Fe ii line velocities, although the use of Fe ii as a diagnostic does imply higher explosion energies. The inferred range of ejecta masses are compatible with intermediate mass (MZAMS ≤ 20M⊙) progenitor stars in binary systems for the majority of SE SNe. Furthermore, our hydrodynamical modeling of the bolometric light curves suggests a significant fraction of the sample may have experienced significant mixing of 56Ni, particularly in the case of SNe Ic.

A UNIFIED ENERGY-RESERVOIR MODEL CONTAINING CONTRIBUTIONS FROM ${}^{56}\mathrm{Ni}$ AND NEUTRON STARS AND ITS IMPLICATION FOR LUMINOUS TYPE Ic SUPERNOVAE

Most type-Ic core-collapse supernovae (CCSNe) produce $^{56}$Ni and neutron stars (NSs) or black holes (BHs). The dipole radiation of nascent NSs has usually been neglected in explaining supernovae (SNe) with peak absolute magnitude $M_{\rm peak}$ in any band are $\gtrsim -19.5$~mag, while the $^{56}$Ni can be neglected in fitting most type-Ic superluminous supernovae (SLSNe Ic) whose $M_{\rm peak}$ in any band are $\lesssim -21$~mag, since the luminosity from a magnetar (highly magnetized NS) can outshine that from a moderate amount of $^{56}$Ni. For luminous SNe Ic with $-21 \lesssim M_{\rm peak}\lesssim -19.5$~mag, however, both contributions from $^{56}$Ni and NSs cannot be neglected without serious modeling, since they are not SLSNe and the $^{56}$Ni mass could be up to $\sim 0.5 M_{\odot}$. In this paper we propose a unified model that contain contributions from both $^{56}$Ni and a nascent NS. We select three luminous SNe Ic-BL, SN~2010ay, SN~2006nx, and SN~14475, and show that, if these SNe are powered by $^{56}$Ni, the ratio of $M_{\rm Ni}$ to $M_{\rm ej}$ are unrealistic. Alternatively, we invoke the magnetar model and the hybrid ($^{56}$Ni + NS) model and find that they can fit the observations, indicating that our models are valid and necessary for luminous SNe Ic. Owing to the lack of late-time photometric data, we cannot break the parameter degeneracy and thus distinguish among the model parameters, but we can expect that future multi-epoch observations of luminous SNe can provide stringent constraints on $^{56}$Ni yields and the parameters of putative magnetars.

THE METAL AVERSION OF LONG-DURATION GAMMA-RAY BURSTS

Recently, it has been suggested that the metallicity aversion of long-duration gamma-ray bursts (LGRBs) is not intrinsic to their formation, but rather a consequence of the anti-correlation between star-formation and metallicity seen in the general galaxy population. To investigate this proposal, we compare the metallicity of the hosts of LGRBs, broad-lined Type Ic (Ic-bl) supernovae (SNe), and Type II SNe to each other and to the metallicity distribution of star-forming galaxies using the SDSS to represent galaxies in the local universe and the TKRS for galaxies at intermediate redshifts. The differing metallicity distributions of the LGRB hosts and the star formation in local galaxies forces us to conclude that the low-metallicity preference of LGRBs is not primarily driven by the anti-correlation between star-formation and metallicity, but rather must be overwhelmingly due to the astrophysics of the LGRBs themselves. Three quarters of our LGRB sample are found at metallicities below 12+log(O/H) < 8.6, while less than a tenth of local star-formation is at similarly low metallicities. However, our SN samples are statistically consistent with the metallicity distribution of the general galaxy population. Using the TKRS population of galaxies, we are able to exclude the possibility that the LGRB host metallicity aversion is caused by the decrease in galaxy metallicity with redshift. The presence of the strong metallicity difference between LGRBs and Ic-bl SNe largely eliminates the possibility that the observed LGRB metallicity bias is a byproduct of a difference in the initial mass functions of the galaxy populations. Rather, metallicity below half-solar must be a fundamental component of the evolutionary process that separates LGRBs from the vast majority of Ic-bl SNe and from the bulk of local star-formation.

A new method for estimating the bolometric properties of Ibc supernovae

The bolometric properties (nickel mass, ejecta mass and kinetic energies) of 61 Ibc supernovae (SNe), including 20 gamma-ray burst and X-ray flash (GRB/XRF), 19 Ib, 13 Ic and nine Ic-BL (broad-lined) SNe are presented. All of the available BVRI photometry in the literature have been collected and used in a new method that utilizes a template supernova (SN 1998bw) and an analytical model to accurately estimate the bolometric properties of each SN. A statistical analysis of the bolometric properties is then performed, where it is found that GRB/XRF SNe are the most energetic, and eject more mass (including nickel content) than Ib, Ic and Ic-BL SNe. The results are then compared to the existing progenitor models of Ibc SNe, where it is concluded that it is highly likely that at least two progenitor channels exist for producing a Ibc SN: most Ibc SNe arise via binary interactions, where the mass of the stellar progenitor is less than what is attributed to a Wolf–Rayet star. Conversely, the progenitors of Ic-BL and GRB/XRF are more massive than those of Ib and Ic SNe, and a key difference between GRB/XRF SNe and Ic-BL SNe is progenitor metallicity, where it is observed that the latter arise from more metal-rich progenitors. As mass loss in massive stars is influenced by metal content, the progenitors of Ic-BL SNe lose more mass, and therefore more angular momentum, before exploding. It is expected that the explosion mechanism in Ic-BL and GRB/XRF SNe is ‘engine-driven’ (i.e. an accreting black hole, or a millisecond magnetar), but the increased mass loss of Ic-BL SNe means the central engine is less powerful than in GRB/XRF SNe. Finally, it is found that the SNe that accompany GRBs and XRFs are statistically indistinguishable, and some mechanism other than metallicity is needed to explain the differences in the high-energy components in these events.

On The Metal Aversion of LGRBs

Recently, it has been suggested that the metallicity aversion of long-duration gamma-ray bursts (LGRBs) is not intrinsic to their formation, but rather a consequence of the anti-correlation between star-formation and metallicity seen in the general galaxy population. To investigate this proposal, we compare the metallicity of the hosts of LGRBs, broad-lined Type Ic (Ic-bl) supernovae (SNe), and Type II SNe to each other and to the metallicity distribution of star-forming galaxies using the SDSS to represent galaxies in the local universe and the TKRS for galaxies at intermediate redshifts. The differing metallicity distributions of the LGRB hosts and the star-formation in local galaxies forces us to conclude that the low-metallicity preference of LGRBs is not primarily driven by the anti-correlation between star-formation and metallicity, but rather must be overwhelmingly due to the astrophysics of the LGRBs themselves. Three quarters of our LGRB sample are found at metallicities below 12+log(O/H) < 8.6, while less than a tenth of local star-formation is at similarly low metallicities. However, our SN samples are statistically consistent with the metallicity distribution of the general galaxy population. Using the TKRS population of galaxies, we are able to exclude the possibility that the LGRB host metallicity aversion is caused by the decrease in galaxy metallicity with redshift. The presence of the strong metallicity difference between LGRBs and Ic-bl SNe largely eliminates the possibility that the observed LGRB metallicity bias is a byproduct of a difference in the initial mass functions of the galaxy populations. Rather, metallicity below half-solar must be a fundamental component of the evolutionary process that separates LGRBs from the vast majority of Ic-bl SNe and from the bulk of local star-formation.

CORE-COLLAPSE SUPERNOVAE FROM THE PALOMAR TRANSIENT FACTORY: INDICATIONS FOR A DIFFERENT POPULATION IN DWARF GALAXIES

We use the first compilation of 72 core-collapse supernovae (SNe) from the Palomar Transient Factory (PTF) to study their observed subtype distribution in dwarf galaxies compared to giant galaxies. Our sample is the largest single-survey, untargeted, spectroscopically classified, homogeneous collection of core-collapse events ever assembled, spanning a wide host-galaxy luminosity range (down to M_r ~ -14 mag) and including a substantial fraction (>20%) of dwarf (M_r >= -18 mag) hosts. We find more core-collapse SNe in dwarf galaxies than expected and several interesting trends emerge. We use detailed subclassifications of stripped-envelope core-collapse SNe and find that all Type I core-collapse events occurring in dwarf galaxies are either SNe Ib or broad-lined SNe Ic (SNe Ic-BL), while "normal" SNe Ic dominate in giant galaxies. We also see a significant excess of SNe IIb in dwarf hosts. We hypothesize that in lower metallicity hosts, metallicity-driven mass loss is reduced, allowing massive stars that would have appeared as "normal" SNe Ic in metal-rich galaxies to retain some He and H, exploding as Ib/IIb events. At the same time, another mechanism allows some stars to undergo extensive stripping and explode as SNe Ic-BL (and presumably also as long-duration gamma-ray bursts). Our results are still limited by small-number statistics, and our measurements of the observed N(Ib/c)/N(II) number ratio in dwarf and giant hosts (0.25_{-0.15}^{+0.3} and 0.23_{-0.08}^{+0.11}, respectively; 1 sigma uncertainties) are consistent with previous studies and theoretical predictions. As additional PTF data accumulate, more robust statistical analyses will be possible, allowing the evolution of massive stars to be probed via the dwarf-galaxy SN population.