Long-rising Type II supernovae from PTF and CCCP
Key Words.:supernovae: general – supernovae: individual: SN 2004ek, SN 2004em, SN 2005ci, SN 2005dp, PTF09gpn, PTF12gcx, PTF12kso, SN 2012gg, LSQ12fjm, SN 1987A, SN 1998A, SN 2009E, SN 2006V, SN 2006au, SN 2000cb, SN 1909A
Context:Supernova (SN) 1987A was a peculiar hydrogen-rich event with a long-rising ( d) light curve, stemming from the explosion of a compact blue supergiant star. Only a few similar events have been presented in the literature in recent decades.
Aims:We present new data for a sample of six long-rising Type II SNe (SNe II), three of which were discovered and observed by the Palomar Transient Factory (PTF) and three observed by the Caltech Core-Collapse Project (CCCP). Our aim is to enlarge this small family of long-rising SNe II, characterizing their differences in terms of progenitor and explosion parameters. We also study the metallicity of their environments.
Methods:Optical light curves, spectra, and host-galaxy properties of these SNe are presented and analyzed. Detailed comparisons with known SN 1987A-like events in the literature are shown, with particular emphasis on the absolute magnitudes, colors, expansion velocities, and host-galaxy metallicities. Bolometric properties are derived from the multiband light curves. By modeling the early-time emission with scaling relations derived from the SuperNova Explosion Code (SNEC) models of MESA progenitor stars, we estimate the progenitor radii of these transients. The modeling of the bolometric light curves also allows us to estimate other progenitor and explosion parameters, such as the ejected Ni mass, the explosion energy, and the ejecta mass.
Results:We present PTF12kso, a long-rising SN II that is estimated to have the largest amount of ejected Ni mass measured for this class. PTF09gpn and PTF12kso are found at the lowest host metallicities observed for this SN group. The variety of early light-curve luminosities depends on the wide range of progenitor radii of these SNe, from a few tens of R (SN 2005ci) up to thousands (SN 2004ek) with some intermediate cases between 100 R (PTF09gpn) and 300 R (SN 2004em).
Conclusions:We confirm that long-rising SNe II with light-curve shapes closely resembling that of SN 1987A generally arise from blue supergiant (BSG) stars. However, some of them, such as SN 2004em, likely have progenitors with larger radii ( R, typical of yellow supergiants) and can thus be regarded as intermediate cases between normal SNe IIP and SN 1987A-like SNe. Some extended red supergiant (RSG) stars such as the progenitor of SN 2004ek can also produce long-rising SNe II if they synthesized a large amount of Ni in the explosion. Low host metallicity is confirmed as a characteristic of the SNe arising from compact BSG stars.
Supernova (SN) 1987A marked the explosion of a blue supergiant (BSG) star (Sanduleak 202, Gilmozzi et al. 1987), and it showed a peculiar long-rising light curve (84 days from explosion to peak) and a hydrogen-rich (Type II SN) spectrum. Given its extraordinary proximity [merely 50 kpc, in the Large Magellanic Cloud (LMC)], SN 1987A is to date the best-studied SN. It changed our understanding of SN explosions (see the reviews by Arnett et al. 1989 and McCray 1993), and it continues to trigger scientific activity; it is the only SN that can be spatially resolved in the optical and was the first one whose transition to a young SN remnant could be observed (Larsson et al. 2011).
In the last decade, a handful of SNe resembling SN 1987A have been presented in the literature. Most of these SN 1987A-like SNe were collected by Pastorello et al. (2012). The best-observed objects are SNe 1998A (Pastorello et al. 2005), 2000cb (Kleiser et al. 2011), 2006V, 2006au (Taddia et al. 2012), and 2009E (Pastorello et al. 2012). Well-sampled optical light curves of the long-rising ( d) SNe II 2004ek, 2004em, and 2005ci by the Caltech Core-Collapse Project (CCCP; Gal-Yam et al. 2007) were shown by Arcavi et al. (2012), and here we present their full datasets, including unpublished spectroscopy. Gonzalez-Gaitan et al. (2015) also discussed slow risers in their sample of SNe II. Some of their SNe lack spectral classifications and might be of Type IIn; however, at least one spectroscopically confirmed SN 1987A-like event is presented (SNLS-07D2an), with a rise time of d. The first SN 1987A-like event was likely SN 1909A (Young & Branch 1988), but for this object there is no spectral classification.
The analysis of the light curves and spectra of the SNe discussed by Pastorello et al. (2012) led to the conclusion that they arise from the explosions of compact ( R) BSG stars, with initial masses of about M, explosion energies of a few erg, and Ni masses of M.
As BSG stars are not expected to explode in the standard stellar evolution models, following the explosion of Sanduleak 202 there were many attempts to explain this surprising fact (see Podsiadlowski 1992 for a review). Models showed that a low metallicity (like that measured in the LMC) and/or fast rotation could make a single M star explode in the blue part of the Hertzsprung-Russell (HR) diagram. We (Taddia et al. 2013) presented a detailed study of the metallicity at the locations of the known SN 1987A-like events, finding that these SNe indeed tend to occur at lower (LMC-like) metallicity than normal SNe IIP. Alternatively, a binary scenario may also explain BSG explosions, and in the specific case of SN 1987A (Podsiadlowski 1992) it may also explain the complex shape of the circumstellar medium (CSM).
Given the rarity of long-rising SNe II, it is important to extend the number of well-observed objects, in order to study the variety of their properties. The observations from the Palomar Transient Factory (PTF; Rahmer et al. 2008; Rau et al. 2009) and from the CCCP that we present here are therefore significant additions to this SN subgroup.
The PTF discoveries are particularly important as they are not biased toward bright (i.e., metal-rich) host galaxies (Arcavi et al. 2010). Our small PTF sample of SN 1987A-like events allowed us to find these SNe at unprecedented low metallicity, with important implications for the progenitor scenario. We notice that other rare and peculiar core-collapse (CC) SNe have been observed in low or very low metallicity environments (e.g., the Type Ic-BL PTF12gzk; Ben-Ami et al. 2012). Other recent examples are the class of hydrogen-poor superluminous SNe (SLSNe; Quimby et al. 2011; Gal-Yam 2012; Lunnan et al. 2014; Leloudas et al. 2015), and that of the SN impostors (Taddia et al. 2015a).
The PTF survey and its continuation (the intermediate PTF, iPTF; Kulkarni 2013) not only extends the domain of SN discovery to low-luminosity hosts, but thanks to the high cadence allows regular discoveries of very young transients (e.g., iPTF13ast; Gal-Yam et al. 2014). Indeed, some of our PTF (but also CCCP) objects were likely discovered and observed soon after explosion. Therefore, some of them show possible signatures of the shock-breakout cooling tail in their optical light curves (Chevalier 1992; Chevalier & Fransson 2008; Rabinak & Waxman 2011; Piro & Nakar 2013). This feature is rarely observed and contains important information to understand the nature of the SN progenitor. First, it allows an estimate of the explosion epoch with small uncertainty, which is crucial for deriving the properties of the progenitor from light-curve and spectral modeling. In addition, especially for objects that are radioactively powered at peak (like our long-rising SNe II), the modeling of the early-time light curve can constrain the radius of the SN progenitor.
Here we present new observations of six long-rising SNe II, which appear rather similar to SN 1987A. The paper is structured as follows. In Sec. 2 we describe the SN sample, and in Sec. 3 we discuss data acquisition and reduction. Section 4 includes the analysis of the host galaxies. We describe and analyze SN photometry and spectroscopy in Sec. 5 and Sec. 6, respectively. Bolometric properties are derived in Sec. 7. Data are modeled in Sec. 8, in order to derive the progenitor parameters. We discuss the results in Sec. 9 and give our conclusions in Sec. 10.
2 The supernova sample
We define a SN II as “long-rising” based on the epoch of the (second) peak of its -band light curve. In Fig. 1 (top panel) we show the distribution of -band rise times of SNe II, with data for normal SNe II (IIP and IIL) and SNe IIb from Gall et al. (2015), Rubin et al. (2015), and Taddia et al. (2015a), as well as data for long-rising SNe II from the light curves discussed in Sec. 5 and displayed in the bottom panel of Fig. 1. SNe showing rise times d are considered “long risers” in this work. Interestingly, among SNe II (excluding SNe IIn, which are powered by CSM interaction), there is a lack of SNe rising in 30–50 d. This is likely due to the fact that the explosion energy and the ejecta mass (up to M) of SNe II appear to be linearly correlated (Utrobin et al. 2010), and their ratios imply long rise times for compact progenitor stars (see Eq. 2).
SNe II with rise times of –25 d are presented by Rubin et al. (2015) (e.g., iPTF14bas), and in that work these are named “extremely slow risers.” However, we include in our definition of slow risers only those SNe with rise times d, to be sure to include only objects that are mainly powered by radioactive decay at peak, like SN 1987A. In SNe IIP and IIL, a rise of only a few weeks in the optical could instead be due to a temperature effect, with most of the flux in the ultraviolet (UV) part of the spectrum until –25 d. This is why we exclude from our sample another CCCP SN II with relatively slow rise (24 d), SN 2005dp (Arcavi et al. 2012). The rising part of the -band light curve of this transient is shown in the bottom panel of Fig. 1, and it is clearly distinct from the family of slow risers resembling SN 1987A. For SNe IIb, the peak is powered by radioactive decay as in SN 1987A, but their shorter rise times are due to the lower ejecta mass as compared to SN 1987A-like SNe. SNe IIb have lost their H envelopes almost entirely prior to exploding.
As shown in the bottom panel of Fig. 1, the slow risers can exhibit a variety of light-curve shapes at early epochs. Their rise can be only a few tenths of a magnitude from the early minimum to the second peak (e.g., PTF09gpn), or it can be several magnitudes (e.g., SN 1987A). We include both of these types of objects in the category of long-rising SNe II or SN 1987A-like events. Their common feature is the late-time peak, which we will show also appears in their bolometric light curves; it is therefore likely to be powered mainly by the decay chain of Ni. We will propose that the variety of the early-time light curves is a result of the explosion of progenitors with different radii (see Sec. 8.3 and the results by Young 2004).
Here we present data for six SNe II with long-rising light curves.
Three of them were discovered and followed by PTF (PTF09gpn, PTF12gcx, PTF12kso), and three were observed by CCCP (SNe 2004ek, 2004em, 2005ci). For the PTF objects we present unpublished light curves and spectra. For the CCCP objects we present unpublished spectra and make use of the light curves already included in Arcavi et al. (2012).
In Table 1 the basic information on the targets is given, including SN coordinates, host-galaxy name and type, redshift, distance, and extinction. The redshifts are obtained from NED
3 Data acquisition and reduction
3.1 Discovery and explosion epochs
In Table 2 we report the last nondetection date and limiting magnitude, the discovery date and magnitude (along with their references), and the assumed explosion epoch for each SN. It is important to constrain the explosion date for each SN, in order to compare the different SNe at similar phases and to derive the progenitor parameters. The last nondetection and the discovery epochs can be used to constrain the explosion epoch, if they occur close in time and the nondetection magnitude limit is sufficiently deep. Alternatively, a power-law (PL) fit to the early-time light curves can be used (e.g., Cao et al. 2013).
For SNe 2004em, 2005ci, PTF09gpn, and PTF12gcx, we used a PL fit to the early light curves to estimate the explosion epoch. For SN 2004em we fit the first six epochs of the bolometric light curve (see Sec. 7). For SN 2005ci, PTF09gpn, and PTF12gcx, we fit the first 11, 3, and 7 epochs of the -band light curves (see Sec. 5), respectively. The uncertainty in the explosion date is determined by the fit uncertainty and/or by the difference between the discovery and the last nondetection epochs.
For SN 2004ek we derive the explosion date from the average between discovery and the last nondetection, and the uncertainty from the interval between these two epochs. We adopt the same method for PTF12kso as we do not have observations of this SN at early times that would allow us to fit a PL.
In the rest of the paper the phase of each object will be defined in days since explosion.
3.2 CCCP supernova data
The CCCP SN light curves were obtained from Arcavi et al. (2012) via Wiserep (Yaron & Gal-Yam 2012). Additional data for SNe 2004ek and 2005ci were obtained from Kleiser et al. (2011, unfiltered) and Tsvetkov (2008, ), respectively. The photometry from Tsvetkov (2008) was shifted by 0.1, 0.2, 0.2, and 0.1 mag (respectively), in order to match the CCCP photometry of SN 2004ek. We report all the magnitudes in Table 3.
Five spectra of each of SNe 2004ek and 2004em were obtained with the Double Spectrograph (DBSP; Oke & Gunn 1982) mounted on the 200-inch Hale telescope at Palomar Observatory (P200), the Low Resolution Imaging Spectrometer (LRIS; Oke et al. 1995) mounted on the Keck 1 10-m telescope, and the Kast spectrograph (Miller & Stone 1993) on the 3-m Shane reflector at Lick Observatory. Three spectra of SN 2005ci were obtained with P200+DBSP. For a detailed log of the CCCP spectral observations, see Table 5. The spectra were reduced in a standard manner, including (i) bias and flatfield corrections, (ii) one-dimensional (1D) spectral extraction from the 2D frames after tracing and choosing proper apertures and background regions, (iii) wavelength calibration with spectra of calibration lamps, and (iv) flux calibration with spectrophotometric standard stars observed during the same night.
3.3 PTF supernova data
In Fig. 2 we show the discovery images, along with reference and subtracted images, of the three PTF SNe. For these SNe -band photometry was obtained with the 48-inch Samuel Oschin Telescope at Palomar Observatory (P48), equipped with the 96 Mpixel mosaic camera CFH12K (Rahmer et al. 2008) and a Mould filter (Ofek et al. 2012). Additional -band photometry of PTF12kso was obtained with the same telescope and instrument. PTF09gpn, PTF12gcx, and PTF12kso were also imaged in the bands with the automated Palomar 60-inch telescope (P60; Cenko et al. 2006). PTF12gcx was not detected in the band, whereas for PTF09gpn -band photometry was also obtained. For PTF12kso, a large number of photometric epochs with the 1-m telescopes of the Las Cumbres Observatory Global Telescope Network (LCOGT; Brown et al. 2013) in the bands were also obtained. For the same SN additional photometry was obtained at Mount Laguna Observatory (MLO; Smith & Nelson 1969) using the telescope, camera, and reduction technique described by Smith et al. (2015).
Point-spread-function (PSF) photometry was obtained on template-subtracted P48 and P60 images using the Palomar Transient Factory Image Differencing and Extraction (PTFIDE) pipeline
A single optical spectrum of PTF09gpn was obtained using the P200+DBSP. Keck-1+LRIS was used to acquire an optical spectrum of the host galaxy.
For PTF12gcx, three spectra were obtained with three different facilities: Keck 2DEep Imaging Multi-Object Spectrograph (DEIMOS; Faber et al. 2003), P200+DBSP, and Lick 3-m+Kast. Four spectra of PTF12kso were obtained with P200+DBSP and Keck 2+DEIMOS, and one with the Copernico 1.82m+AFOSC (obtained from the Asiago Supernova Classification program
4 Host galaxies
We have already examined the host galaxies of the CCCP SN 1987A-like SNe in Taddia et al. (2013) along with those of the entire literature sample of long-rising SNe II. There we found that SN 1987A-like SN locations are more metal-poor than those of normal SNe II. These peculiar SNe are located either in the metal-poor outskirts of bright galaxies or in dwarf hosts. Here we present and analyze spectra and images of the PTF objects, which are shown in Fig. 3.
PTF09gpn was located in a faint host galaxy characterized by an apparent magnitude of 19.13 in the band (Arcavi et al. 2010). A color image of the galaxy is shown in Fig. 3 (top-left panel). The host-galaxy spectrum (top-right panel of Fig. 3) reveals several emission lines (Balmer lines, [O iii] 4363, 4959, 5007, [O ii] 3727, [N ii] 6584, [S ii] 6717, 6731; see Fig. 3).
The absolute magnitude of the host galaxy in the band is mag, which is that of a dwarf galaxy. By using the value of the absolute magnitude in the band and Eq. 1 in Arcavi et al. (2010), we estimated the metallicity to be 0.13 Z, or 12log(O/H) = 7.80 (where 12log(O/H) = 8.69 dex; Asplund et al. 2009). This would place PTF09gpn at the lowest end of the known metallicity distribution of SN 1987A-like events (Taddia et al. 2013).
To provide a more direct measurement of the host metallicity, we made use of the emission-line fluxes in the host-galaxy spectrum. We first modeled the spectrum to account for the Balmer absorption lines due to the underlying stellar population, using the STARLIGHT code (Cid Fernandes et al. 2005)
The bright [O iii] emission makes the spectrum of this host similar to those of the so-called “green pea” (G.P.) galaxies (Cardamone et al. 2009) and H-poor SLSN hosts (Leloudas et al. 2015), which often show a preference for extreme emission-line galaxies (Amorín et al. 2015). The [O iii] 5007/[O ii] 3727 ratio is 4.2, consistent with those of G.P. galaxies (Jaskot & Oey 2013). The diameter of this galaxy is ″, which corresponds to kpc, also consistent with the size of G.P. galaxies.
The H luminosity is (H) = erg s. This gives an estimate of the star-formation rate (SFR) by applying the formula of Kennicutt & Evans (2012):log(SFR) = log(H)41.27. We obtain SFR = 0.07 M yr. The spectrum was calibrated against the -band magnitude of the galaxy.
From the STARLIGHT fit we also obtained a value for the stellar mass ( M), which implies a specific star-formation rate (sSFR) of yr.
PTF12gcx was hosted by a spiral galaxy, SDSS J154417.02095743.8, characterized by SDSS apparent magnitudes of 18.08(0.04), 16.92(0.01), 16.51(0.01), 16.25(0.01), 16.06(0.02) mag. A color image of the galaxy from SDSS is shown in Fig. 3 (center-left panel). The spectrum of its nucleus (see Fig. 3, center-right panel) reveals the same emission lines as shown by the host of PTF09gpn, with the exception of [O iii] 4363 and H, which are too faint to be detected in this case.
The absolute magnitude of the host galaxy in the band is mag, which is that of a bright spiral galaxy. Again, by using the value of the absolute magnitude in the band and Arcavi et al. (2010), we estimated a metallicity of 0.75 Z, or 12log(O/H) = 8.56. This places PTF12gcx at the highest end of the known metallicity distribution of SN 1987A-like events (Taddia et al. 2013), together with SNe 1998A, 2004em, and 2006au.
We made use of the emission-line fluxes (see Table 6) in the host-galaxy spectrum and strong-line diagnostics to further constrain the metallicity. The SN exploded relatively close to the host center (at 4 separation), so the SDSS spectrum (obtained with a 3-wide fiber) can provide a good estimate of the metallicity at the position of PTF12gcx. N2 and O3N2 indicate 12log(O/H) = and 12log(O/H) = (respectively), confirming that PTF12gcx is at relatively high metallicity. R23 indicates a slightly lower oxygen abundance of 12log(O/H) = (see Table 6).
The H luminosity is (H) = erg s, which corresponds to a SFR of 0.13 M yr (Kennicutt & Evans 2012). The spectrum was calibrated against the SDSS -band magnitude of the galaxy. With M, we obtain sSFR yr.
PTF12kso was hosted by an extremely faint galaxy, barely detected in the SDSS images (see Fig. 3, bottom-left panel). Aperture photometry on stacked P48 -band images reveals a host apparent magnitude of 20.64 mag. The absolute magnitude of mag implies a metallicity of Z, or 12log(O/H) = 7.76, which makes the host of PTF12kso one of the most metal-poor SN 1987A-like SN hosts observed so far.
The spectrum of this faint host galaxy exhibits strong H emission and a faint continuum. We use the flux of H and the limit on the flux of [N ii] to estimate an oxygen abundance of 12log(O/H) with the N2 method. The SFR obtained from the H luminosity is SFR = 0.015 M yr. The spectral continuum is too faint and noisy to obtain a reliable estimate from the STARLIGHT fit.
In summary, two out of our three new metallicity measurements push the oxygen-abundance distribution of the host galaxies of SN 1987A-like SNe shown in Taddia et al. (2013) to even lower values, making the SN 1987A-like SNe one of the most metal-poor populations of CC SNe, together with SN impostors (Taddia et al. 2015b) and H-poor SLSN hosts (Leloudas et al. 2015).
5 Supernova light curves and colors
5.1 Light-curve shape
PTF09gpn was imaged in the bands, and the light curves (with the exception of the late -band epochs) are reported in Fig. 4 (left-hand panel). The -band light curve of PTF09gpn covers the first d. The -band light curve sampling is rather frequent (2.4 d average cadence) until d. Several epochs were taken at late times (300 to 350 d). In the other bands the follow-up observations started weeks after discovery, with an average cadence of 4.1, 3.6, 5.7, and 5.7 d for , , , and (respectively). In these filters the last photometric epoch was taken at d () and d (). The -band light curve shows a rise of mag in the first 10 d. Then it settles onto an early plateau between days 10 and 30, and it finally starts a long rise to peak ( mag brighter than the plateau) that occurs at d. All of the other bands reach maximum around the same time, and between days 10 and 30 they also exhibit an early-time plateau, preceded by a decline in the case of and . Since the follow-up observations started 2 weeks later for these filters than in the band, no early-time rise could be detected. The overplotted and scaled light curves of SN 2006au (from Taddia et al. 2012) clearly show the close resemblance of PTF09gpn to this SN 1987A-like SN.
PTF12gcx was imaged in the band with P48. The light curve (central panel of Fig. 4) is characterized by frequent sampling, with the SN being imaged every 1.3 nights on average during the rise ( mag from discovery to peak). Three epochs (the last at d) were taken after peak, which occurred at day 57. After peak the decline is of mag in d. The light-curve shape is very similar to that of SN 1987A, which is also shown in the figure. A difference is that PTF12gcx rises faster in the first days after discovery and thus peaks earlier than SN 1987A. A few epochs with P60 in were also taken after the peak.
The light curves of PTF12kso are shown in the right-hand panel of Fig. 4. This SN was imaged in the bands. The first observation in with P48 occurred at d. Then, after a gap of d in the light-curve coverage, at d intense follow-up observations started in all the bands with P60 and in with LCOGT. The SN was observed and detected by LCOGT almost every night from a few days before peak (which occurred at 61 d in ) until d ( d in and ). P60 took additional images up to 200 d. After peak the light curves decline by mag in d, but then after 120 d the decline rate becomes smaller by about 50%, with the light curve setting onto a tail very similar to that observed in SN 1987A.
Rise times and peak magnitudes for the PTF SNe are summarized in Table 7.
SN 2004ek was observed in by Arcavi et al. (2012) and Tsvetkov (2008). The light curves are shown in Fig. 5 (top-left panel). The monitoring started at d and, before maximum (which occurred at d in ), the light-curve sampling was rather frequent, with the SN observed every 2.7 nights on average. After peak, a few epochs were obtained between and d. The last observations occurred at almost 290 d in all bands except , for which the last epoch was obtained at d. The light curves of SN 2004ek show an early decline up to –32 d in all the bands except , which shows an early plateau of comparable length. The decline is faster in the bluer bands, with and exhibiting a decline of and 0.2 mag from the early peak, respectively. Following this phase, the bands show a slow rise to peak ( mag in days), which occurred at days. At the same epochs, the band shows a plateau. The early, long-lasting declining phase and the slow rise to maximum shown by SN 2004ek make this SN a very peculiar object — sort of an intermediate case between SN 1987A and a normal SN IIP such as SN 1999em, as we show in Fig. 5. The post-maximum decline up to 160 d is also quite slow in all of the redder bands, as compared with SNe 1987A and 1999em. However, the last epoch obtained at 290 d shows a difference in magnitude from peak comparable to that of SN 1999em and larger than that of SN 1987A.
SN 2004em has light-curve coverage and shape that are rather similar to those of SN 2004ek (see Fig. 5, top-right panel). For this SN the monitoring began at weeks, and most of the epochs were obtained before maximum (which occurred at d). Before maximum the SN was observed once per week on average. A few additional epochs were obtained at d, and the last one at 275 d. This SN shows a 0.6 mag decline in the -band light curve in the first days. At the same epochs, the and bands appear constant, whereas the band first rises by 0.5 mag for weeks, and then also shows a short plateau up to days. After the early phase, all of the light curves rise to peak. The rise to peak is steeper than in the case of SN 2004ek, reaching about 1 mag over 50 d in . After peak, the decline is slower than for SN 1987A, as shown in the figure. The last epoch shows that the difference in magnitude from peak to tail is similar to that of SN 1987A.
The light curves of SN 2005ci are rather similar to those of SN 1987A, as shown in Fig. 5 (bottom-left panel). This SN was observed by Arcavi et al. (2012) in and unfiltered () by Kleiser et al. (2011). The CCCP follow-up observations started at 11 d, with an average cadence of 4 nights and the last epoch taken at 130 d. The unfiltered light curve has the first epoch only a few hours after explosion, and covers the first 44 d. The rise to peak, which occurs at d in , is faster than for SN 1987A in the first –40 days in all bands. Thereafter, the slope becomes less steep ( mag in d in ). The band shows hints of an early plateau before day 20. Its shape is flatter than that of SN 1987A. The post-peak decline ( mag in d in the band) is similar to that of SN 1987A.
Rise times and peak magnitudes for the CCCP SNe are summarized in Table 7.
5.2 Absolute magnitudes
In Fig. 6 we plot the absolute magnitudes in , , , , and filters of all our SNe (colored symbols) and of other well-observed long-rising SNe II resembling SN 1987A (black symbols and lines). The literature data were taken from Kleiser et al. (2011), Pastorello et al. (2005), Pastorello et al. (2012), Taddia et al. (2012), and reference therein. For our SNe we adopted the distances and extinctions listed in Table 1. For the SNe in the literature we obtained redshift-independent distances from NED when available; otherwise we adopted WMAP5 cosmological parameters and a redshift from NED, to be consistent with the distances computed for our objects. The Galactic extinctions were obtained from NED. The explosion epochs and the host-galaxy extinctions were obtained from the literature.
The only band in which all our objects and those from the literature have a frequently sampled and extended light curve is . The comparison in this band reveals that the magnitude of the (second) peak ranges from mag (as in the case of SNe 2005ci, 2009E, and PTF09gpn) to mag (SN 2004ek). One PTF SN (PTF12kso) and two CCCP objects (SNe 2004ek and 2004em) are more luminous than all of the other SN 1987A-like events previously published (among those SNe, SN 2006V was hitherto the most luminous), peaking at mag. PTF09gpn and SN 2005ci are instead among the least luminous SN 1987A-like SNe, being only marginally brighter than SN 2009E. On the tail (i.e., after d), the SNe which have brighter peaks also have brighter tails. However, at very late epochs ( d), the absolute magnitudes of SNe 2004ek and 2004em are slightly brighter (up to 0.4 mag) than those of SN 1998A, which at peak was up to 1.5 mag fainter. Compared to the band, the and bands exhibit wider ranges of peak magnitudes (3.5 and 2.5 mag, respectively).
At early phases, when most of the flux is in the bluer bands, ranges from mag (SN 2004ek) down to mag (SN 2005ci). In between, we find SNe 1987A, 2004em, 2006au, and PTF09gpn, showing early declines.
The large range of absolute magnitudes has important implications for the progenitor parameters, in particular for the Ni mass and the progenitor radius (see Sec. 8).
Adopting the light curves and the extinctions used to compute the absolute magnitudes, we also compared the colors of our SNe to those of other long-rising SNe. In Fig. 7 we plot three different optical colors: , , and . The and color evolutions are similar for all of the SNe, with bluer colors at early epochs which progressively become redder. Among the reddest SNe we find SN 2005ci, similar to SN 1987A. SNe 2004ek, 2004em, PTF09gpn, and PTF12kso are bluer, similar to SNe 2006V and 2006au.
The color looks flatter for all of the SNe, as its value is more influenced by the different strengths of the spectral lines (e.g., H and Ca II) than by the continuum.
6 Supernova spectra
6.1 PTF SN spectra
In Fig. 8 we plot the spectral sequences of PTF09gpn, PTF12gcx, and PTF12kso. The spectra are made publicly available via WISeREP (Yaron & Gal-Yam 2012). All of the spectra were calibrated against -band photometry and corrected for reddening.
In Fig. 8 (top panel) we show the only spectrum obtained of PTF09gpn, which was taken at d, during the early peak of the -band light curve, as indicated by the vertical dashed line in Fig. 4 (left-hand panel). The Balmer lines (H, H, H, H, H) dominate the spectrum of PTF09gpn and are labeled in the figure. All of them exhibit a broad P-Cygni profile. The continuum appears quite blue, as the spectrum was obtained at early times.
In the central panel of Fig. 8 three spectra of PTF12gcx are reported. The first two spectra, obtained on the rise of the light curve at 39 and 49 d, show clear P-Cygni profiles for the Balmer lines (from H to H), for Na i D, and for the Ca ii triplet. Redward of Na i D there may be traces of high-velocity H absorption or Si ii 6355. There are also signatures of Fe ii 5169 in absorption. The continuum spectral energy distribution (SED) peaks in the band. The last spectrum, obtained after peak, has low signal-to-noise ratio and only reveals a faint H feature.
The spectral sequence of PTF12kso is shown in the bottom panel of Fig. 8. Five spectra were obtained, three of them before peak (at 61 to 68 d) and two after peak (at 94 and 124 d). In the early-time spectra, besides the Balmer lines down to H, we detect P-Cygni features associated with Na i D (possibly blended with He i 5875), Ca ii triplet, Fe ii triplet, and probably Ca ii H&K. There are also traces of Ba ii at 6142 Å and 4554 Å. Ba ii 6142 may be blended with high-velocity H absorption or Si II 6355. The continuum SED peaks in the band. In the last spectrum, the absorption part of the P-Cygni features decreases its strength in all the lines except for Na i D. At these late epochs, Ca ii 7291, 7323 emerges in emission, the continuum is faint and almost constant, and the lines dominate the SED.
6.2 CCCP SN spectra
In Fig. 9 we plot the spectral sequences of SNe 2004ek, 2004em, and 2005ci. The spectra are made publicly available via WISeREP (Yaron & Gal-Yam 2012). All the spectra were calibrated against -band photometry and corrected for reddening.
The spectra of SN 2004ek were obtained at 21, 22, 78 (pre-peak), 103 and 134 (post-peak) days. The first two spectra were taken at the epoch of the early cooling phase, and they show a blue continuum where diluted P-Cygni profiles of Balmer lines are barely visible. Fe ii 5169 and possibly He i 5875 in absorption are also visible. The later spectra peak in the continuum around 5000–5500 Å, and show strong broad H emission with a marginal and narrow P-Cygni absorption. The other Balmer lines and in particular H show prominent P-Cygni features. Na i D emerges and Fe ii 5169 (along with Fe ii 5018) becomes more evident. The Ca ii triplet emerges and traces of Ba ii 4554, 6142 are observable. Possibly Ca ii H&K is also detected.
The first spectrum of SN 2004em was obtained at 8 d, at the cooling phase. Because of the early phase, it appears much bluer than the other spectra, and shows broad Balmer P-Cygni profiles. The next three spectra were also obtained before maximum, at 33, 65, and 81 d. Although the second spectrum was taken on the decline of the -band light curve, it does not appear as blue as the first two spectra of SN 2004ek or as the previous spectrum of SN 2004em. These spectra have a continuum that peaks in the band, and are dominated by the typical SN II P-Cygni Balmer profiles. As in the spectra of the SNe that were previously discussed, they exhibit broad P-Cygni profiles for Na i D and for the Ca ii and Fe ii triplets. In the last spectrum, taken at 202 d, Ca ii 7291, 7323 emerges in emission, the Ca ii triplet does not show absorption, the continuum becomes weak, and H dominates the SED with a strong broad emission and residual absorption.
SN 2005ci has three spectra obtained at 44, 72, and 98 d (only one after the peak and before the tail). These spectra are particularly rich in features and are dominated by P-Cygni Balmer lines. All of the other lines mentioned for the other SNe are visible in emission and absorption, including Na i D, Ca ii triplet, Fe ii 4924, 5018, 5169, Ca ii 7291, 7323 (in emission), and Ba ii 4554, 6142. The spectral continua are quite red, peaking around 6000 Å.
6.3 Spectral comparison and expansion velocities
For a more direct spectral comparison among our SNe and other SN 1987A-like events, we plot a sequence of spectra taken at early epochs in the top panel of Fig. 10, one taken around peak in the central panel, and three spectra obtained at late epochs in the bottom panel.
In the top panel it is evident that at early epochs these long-rising SNe II show Balmer and (later) Fe features diluted by a strong black-body continuum. The case of SN 2004ek is particularly extreme, with the spectral lines almost completely diluted even after 3 weeks from the explosion. As the ejecta expand, the continuum becomes less prominent and more features emerge. These lines characterize the spectra taken around peak, shown in the central panel. Here we can see the SNe exhibit different line broadening, corresponding to different expansion velocities.
The late-time spectra shown in the bottom panel highlight the transition to the nebular phase, with a faint continuum and strong emission lines, in particular H and Ca ii.
We measured the expansion velocities of our SNe from the P-Cygni minima of different lines at all available epochs, and display the results in Fig. 11. Balmer lines, Fe ii 5169, Na i D, Ba ii 6142, and Ca ii 8662 are analyzed.
Starting with the Fe ii 5169 line velocities, which better represent the photospheric velocities (Dessart & Hillier 2005), we can see how our SNe are faster than SN 1987A and SN 2009E, and similar to the more energetic SNe 1998A, 2006V and 2006au. In particular, SNe 2004ek and 2004em exhibit the fastest photospheres. SN 2005ci is the slowest in our sample. These results are confirmed when we look at the Ba ii lines, another good indicator of the photospheric velocity.
Looking at the bright H line, we observe that PTF12kso and PTF12gcx have larger velocities compared to the other SNe, which show similar expansion velocities as did SN 1987A. This is also observed in H, but the difference is less significant.
Na i D and especially the Ca ii lines show similar velocities in all of the SNe.
7 Bolometric properties
For the CCCP and PTF SNe (except for PTF12gcx), we constructed the SEDs after having interpolated the light curves to the same epochs and converted the extinction-corrected magnitudes into fluxes at the effective wavelength of each filter. The bolometric light curves were computed following the prescriptions by Lyman et al. (2014), to account for the missing near-infrared (NIR) and UV flux: (a) the optical flux was obtained as the trapezoidal integral of the fluxes. (b) The UV flux was computed as the integral of the black-body (BB) fit to the entire SED blueward of the effective -band wavelength during the early cooling phase, and as the trapezoidal integral blueward of the effective -band wavelength with Å) = 0 at later epochs. The cooling phase coincides with the epochs when the early -band light curve is observed declining. (c) the NIR flux was obtained as the integral of the best BB fit in the region redward of the effective -band wavelength. The BB was fit to the entire SED during the cooling phase, and only to at later epochs, as the -band flux is potentially affected by line blanketing. The sum of optical, UV, and NIR fluxes was then multiplied by , with being the distance of the SN, to obtain the bolometric luminosity (). This is shown in Fig. 12 (top panel). In the case of SN 2005ci, we verified that obtaining the bolometric light curve by simply integrating the best BB fit to the entire () SED at all epochs would not have changed the result by much, with an average difference in luminosity of %. For this reason, when the band was missing at late epochs (e.g., SN 2004em), we adopted the integral of the BB fit to as the bolometric flux. For PTF12gcx, where we only have a well-sampled -band light curve, the bolometric light curve was obtained using the bolometric corrections from Lyman et al. (2014) assuming equal to that of SN 1987A (see Fig. 7).
With the BB fit to the SED we also determined the BB radius () and the temperature (), which are shown in the central and bottom panels of Fig. 12. Again, fitting only or the entire SED does not have a large impact on the final results, with of SN 2005ci differing by 3% on average and by 12% at most. The radius differs by 6% on average and by 27% at most.
With the same method outlined for CCCP and PTF SNe, we also derived , , and for other well-observed SN 1987A-like SNe from the literature (data from Pastorello et al. 2005; Kleiser et al. 2011; Taddia et al. 2012; Pastorello et al. 2012). The results are shown in Fig. 12, where the top panel confirms that among the brightest long-rising SNe II we find SN 2004ek, PTF12kso, and SN 2004em. The faintest are SN 2009E, PTF09gpn, and SN 2005ci. Long rising SNe II from the literature (black lines and symbols) are on average fainter than the SNe in our sample (colored symbols), with SNe 2006V and 2006au being the most luminous. At late epochs, all of the light curves show a decline rate identical to that of SN 1987A, which is consistent with the Co decay rate, with the exception of SNe 2004ek and 2004em (which decline faster). SNe 2004ek, 2004em, 2005ci, and PTF09gpn have a decreasing luminosity and a decreasing temperature (central panel of Fig. 12) at early epochs. The temperature profiles are rather similar for all the SNe, with an early cooling (up to d) followed by an almost constant temperature between and K depending on the object. The BB radii (bottom panel of Fig. 12) are on the order of cm, with PTF12kso and SN 2004ek exhibiting the largest values. After a fast expansion, they peak around the epoch of maximum luminosity and then slowly decrease. In summary, our SN sample of long-rising SNe II increases the spread of the observed bolometric properties and exhibits large variety.
8.1 Ni mass from the linear decaying tail of the bolometric light curve
We can use the late part of the bolometric light curves to estimate the amount of Ni ejected in the SN. The light-curve decline at this phase is linear in magnitude and mag every 100 days (the decay rate of Co). As the decay times of Ni and Co ( and ) are known, as well as their specific energy generation rates ( and ), it is possible to constrain the Ni mass () by fitting the following expression to the bolometric light curve:
For PTF12kso, SN 2004em, and SN 2005ci, we can fit at least two epochs on the linear decline, while SN 2004ek was observed only once at late epochs. For PTF09gpn and PTF12gcx, we can set only an upper limit on the Ni mass, as our last bolometric epochs coincide with the beginning of the light-curve tails.
For PTF09gpn we also have -band measurements at very late epochs. If we compare the absolute magnitude at these epochs with that of SN 2009E, the two SNe turn out to have similar Ni masses. Pastorello et al. (2012) estimated a Ni mass of 0.04 M from the quasi-bolometric light curve of SN 2009E, which we adopt as the best value for PTF09gpn.
The objects with the largest and smallest Ni masses are PTF12kso (0.23 M) and PTF09gpn (0.04 M), respectively. All of the results are reported in Table 8.
The uncertainties in these values are mainly affected by the uncertainty in the distance and in the explosion date. There could also be uncertainty related to the host-galaxy extinction. Therefore, we estimate the total error in the Ni mass to be at least 10%.
Our sample of long-rising SNe II extends the range of Ni masses inferred for SN 1987A-like SNe (see Table 5 of Pastorello et al. 2012), with PTF12kso being the most Ni rich, about four times more than for SN 1987A. When compared to the Ni mass distribution of normal SNe II from Hamuy (2003), Inserra et al. (2013), Anderson et al. (2014), and Rubin et al. (2015), it appears that long-rising SNe II produce more Ni on average (see Fig. 13). However, finding SN 1987A-like SNe with low Ni masses is more difficult than discovering SNe IIP with the same Ni masses, as their lower cooling-envelope luminosity at early epochs makes them fainter on the rise. Even accounting for this bias, long-rising SNe II seem to have a different distribution of Ni masses as compared to SNe IIP. A likely explanation for this result is that the progenitors of SN 1987A-like events are on average more massive than those of SNe IIP, and therefore they also produce more Ni, as these two quantities show a correlation (see, e.g., Utrobin & Chugai 2011). We also notice that for seven SN 1987A-like events the Ni mass is derived from at most two points on the radioactive tail. This might introduce a bias in the Ni mass estimates if gamma-rays are escaping.
Based on the Ni mass estimates from the tail of the bolometric light curves of our SNe and other SN 1987A-like SNe, we can conclude that the amount of synthesized Ni in long-rising SNe II can vary widely, between 0.04 M and 0.23 M. The lower value is consistent with those of normal SNe IIP such as SN 1999em, the higher one with those of luminous SNe IIP (e.g., SN 1992af; Nadyozhin 2003) or moderately faint SNe Ibc (e.g., Taddia et al. 2015a). In Fig. 14 we can see how our sample (colored symbols) extends the range of Ni masses up to 0.23 M compared to the previously known long-rising SNe II.
8.2 Explosion energy and ejecta mass from scaling relations
In order to estimate the explosion energy () and the ejecta mass (), we need to combine the information about the bolometric rise time and the photospheric velocity, as done by Taddia et al. (2012). For radioactively powered SNe, the diffusion time given by Arnett (1979) is , and . We assume the same mean opacity for all our SNe, and adopt the Fe ii 5169 velocity (see Fig. 11) around peak as (for PTF09gpn we adopted the velocity of SN 1987A) and the bolometric rise time as . Using values from SN 1987A, (87A) erg, (87A) = 14 M (Blinnikov et al. 2000), (87A) = 84 d, and (87A) = 2200 km s, we obtain estimates of and for the other SNe by scaling. These are listed in Table 8.
The long rise times and fast expansion velocities of SNe 2004ek and 2004em imply large ejecta masses (–40 M) and energies ( erg). Lower values characterize the other SNe, which however show larger ejecta masses and energies than SN 1987A, with the exception of PTF09gpn. The SNe observed in our sample are found to be more energetic and have larger masses (on average) than the old sample, as shown in Fig. 14. It appears from the same figure that the ejecta mass and the explosion energy are correlated, as well as the Ni mass and the explosion energy. The estimates of and are based on simple scaling relations and therefore must be taken with caution. We will show in Sec. 8.3 that the rise times of hydrodynamical models show a more complicated dependence on the progenitor parameters as compared to what we used here, including a dependence on the degree of Ni mixing and on the Ni mass.
8.3 Hydrodynamical modeling
In Figs. 1 and 12 we showed how the early-time shape of our long-rising SNe II is different for the different events. It is known, also from analytic modeling (see, e.g., Chevalier & Fransson 2008), that SN progenitor radii influence the early SN emission, with the more luminous early light curves corresponding to the SNe with larger progenitors.
Here we aim to demonstrate that we need progenitor radii between that of a RSG ( R, typical of SN IIP progenitors) and that of a compact BSG ( R, progenitor of a canonical SN 1987A-like SN) to explain the early light-curve dispersion of these objects. To do this, we made use of open-source codes to evolve massive and hydrogen-rich progenitor stars and explode them, producing bolometric light curves that were compared to our observations.
We used the Modules for Experiments in Stellar Astrophysics (MESA; Paxton et al. 2011) to evolve , 20, 25, and 30 M stars with different metallicities (from to ) and rotation velocities (–600 km s), until they end their lives at different positions in the HR diagram. We thus constructed a sequence of models with stellar radii from to R for stars with M. We decided to produce stars with compact radii for the other initial masses (15, 25, and 30 M). We prohibited mass loss to preserve a large H envelope and did not include overshooting in our computations. The HR diagram including our models is shown in Fig. 15. Low metallicity and fast rotation favor the production of BSG stars in the HR diagram, whereas high metallicity and low rotation velocity produce RSG stars as SN progenitors (e.g., Podsiadlowski 1992).
The aim of this exercise was to construct a sequence of similar stars of different final radii to use as input for the SuperNova Explosion Code (SNEC; Morozova et al. 2015), in order to produce bolometric light curves and analyze the effect of the different parameters on the SN light curves. We produced a set of hydrodynamical models with a range of progenitor radii (), explosion energies (), final masses ( M, as we inhibited the mass loss and consider the formation of a central compact object of 1.4 M), Ni masses (), and degrees of Ni mixing (). The degree of Ni mixing is expressed as a fraction of the final mass, and the Ni is uniformly distributed in the region where it is present. For example, a of 0.9 implies that the Ni mass is uniformly distributed up to the radius that includes 90% of the final progenitor mass.
If we keep the same mass ( M), the same Ni mass (0.075 M), and the same explosion energy ( erg = 1 B), we can see in Fig. 16 (top-left panel) that for compact radii, the light curve is very similar to that of SN 1987A. For larger radii, the luminosity of the early part (before the light-curve rise) becomes larger, until the SN light curve reaches that of a normal SN IIP. This occurs for radii R for our given parameter set. The same result was found by Young (2004). For large radii, the decline in luminosity is also slower than for small radii.
In Fig. 16 (top-central panel) we show that, for a range of final masses from 15 to 30 M, if we keep the same ratio between and , as well as approximately the same , then the rise time and the early luminosity of the light curves are generically very similar.
For larger Ni mass (keeping the same M, R, and erg = 1 B), the rise time is longer; for larger Ni mixing it is shorter (bottom panels of Fig. 16). The effects of Ni mixing on the peak was also investigated by Young (2004).
As we could not produce all the models needed to scan the entire space of the progenitor parameters, we made use of our set of models to derive scaling relations for the time of the peak, luminosity at 10 d, and time of the rise. This allowed us to constrain the degree of Ni mixing, , and the progenitor radius for our SNe.
We measured the peak epochs of these models and we investigated how they change depending on each progenitor parameter. We found that
The luminosity of the early part (measured at 10 d after explosion) scales approximately like
for R. The time at which the light curve starts to rise () after the early cooling depends mainly on the degree of Ni mixing as
We also obtained estimates of the progenitor radius from these relations by scaling with the parameters of the hydrodynamical model of SN 2005ci, shown in the bottom-left panel of Fig. 16. The bolometric light curve and the photospheric velocity of SN 2005ci are well reproduced by a model with M, R, B, M, and . The light curve of SN 2005ci is rather similar to that of SN 2000cb (Kleiser et al. 2011), which was modeled by Utrobin & Chugai (2011), who found a similar mass ( M) and progenitor radius ( R), but higher energy ( B).
For each of our SNe we make use of the measured in Sec. 8.1. We measure the epoch of the light-curve rise () to assess the degree of via Eq. 4 and scaling to SN 2005ci. By measuring the rise time we make use of Eq. 2 to estimate by scaling to SN 2005ci. Then we can use Eq. 3 and the obtained values of to estimate the radius by scaling to SN 2005ci.
In Table 8 we report the parameter estimates obtained from these scaling relations. We obtain radii of about 320–350 R for SN 2004em and PTF12gcx, of R for PTF09gpn, and of a few thousands R for SN 2004ek. Therefore, we suggest that we have observed H-rich SNe with progenitor radii that are intermediate between those of RSGs and BSGs, implying that H-rich stars in the A to K spectral classes can also explode and become SNe II. Moreover, we found that there is a large variety in the degree of Ni mixing, ranging from 25% to 100% of the final mass.
In this work we have added three new and well-observed long-rising SNe II from PTF and analyzed three more SNe from CCCP. These objects are an important addition to this SN class as long-rising SNe II are extremely rare events. Pastorello et al. (2012) estimated that SN 1987A-like SNe make up 1–3% of CC SNe based on the data from the Asiago Supernova Catalogue. Also, Kleiser et al. (2011) derived a similar value based on the data from the Lick Observatory SN Search (Filippenko et al. 2001; Leaman et al. 2011; Li et al. 2011a, b). If we consider the CCCP SN program, which observed nearby CC SNe, three SN 1987A-like SNe (from a total of 62 CC SNe) were followed, suggesting a fraction of 5%. The PTF/iPTF survey is better suited to estimate the rate as it is not a targeted search. Within a distance modulus of mag, we found only PTF09gpn out of a total of 108 CC SNe observed within the same volume. This suggests a fraction of 1%, consistent with the value suggested by Pastorello et al. (2012).
The progenitor parameters suggest stars with a radius smaller than that of a typical RSG, which is R (with the exception of SN 2004ek).
We estimated the progenitor radii by using hydrodynamical models. In particular, we obtained a good fit to the bolometric luminosity and photospheric velocity of SN 2005ci. Thereafter, using the hydrodynamical models, we derived scaling relations for the peak epoch, the early luminosity, and the epoch of the light-curve rise, which depends on , , , and . We confirmed that the early luminosity depends on the ratio and on the with similar power-law indexes. However, compared to the scaling relation based on the diffusion time from the model of Arnett, we found important differences. For instance, we found that depends on instead of . This was also found by Utrobin (2005) for his set of hydrodynamical models in the envelope-mass range 15–21 M. We also considered the dependence of on and , which were previously ignored.
With the relations estimated from this series of hydrodynamical models, we obtained progenitor radii of R for SN 2005ci, R for PTF09gpn, and –350 R for SN 2004em and PTF12gcx. Considering that the ejecta masses range from to M, the radii of these stars imply a range of progenitor stars from BSGs (SN 2005ci and PTF09gpn) to yellow supergiants (YSGs; e.g., SN 2004em and PTF12gcx), as can be seen in the HR diagram of Taddia et al. (2012) and in Fig. 15.
For SN 2004ek, we obtained a large progenitor radius of a few thousand R, typical of a RSG. Its bolometric light curve shows a late-time peak, but the rise occurring after the initial luminosity drop is very shallow compared to that of SN 1987A, and the early-time light curve is very bright. We can interpret the light curve of SN 2004ek as the combination of radioactive emission from its large Ni mass (important only at peak and even more at later epochs) and cooling emission from its large envelope (dominating the early epochs).
For two host galaxies of our PTF SNe (PTF09gpn and PTF12kso), we have measured very low metallicity ( Z), lower than the metallicity measured in the Small Magellanic Cloud (SMC). A low metallicity ( Z) can favor a blue solution in the HR diagram for a single star, with the star exploding as a BSG (see Podsiadlowski 1992 and Weiss 1989). In this case, the SN progenitors could have spent their entire lives in the BSG stage, and did not have a RSG phase such as that of SN 1987A. These two SNe confirm the trend of long-rising SNe II at low metallicity found by Taddia et al. (2013).
We have added and analyzed six new long-rising SNe II to this rare family, three from PTF and three from CCCP. These objects present similarities in the light curves, but also interesting differences, especially at early epochs.
Our SN sample expands the range of Ni masses, explosion energies, and ejecta masses for the long-rising SNe II.
The radii and the large ejecta masses found for our SNe suggest a BSG to YSG origin for these events. Canonical SN 1987A-like SNe arise from compact BSGs (e.g., SN 2005ci), while SNe with more luminous early phases arise from stars with a radius of a few hundred R (e.g., SN 2004em). We also found that SN 2004ek comes from an even more extended RSG progenitor star, but shows a long-rising light curve owing to the large amount of ejected Ni.
We found PTF09gpn and PTF12kso at the lowest host metallicities ever estimated for long-rising SNe II, confirming that low host metallicity is an important ingredient for producing peculiar SNe II, whose radii are typically shorter than those of normal SN IIP progenitors.
Acknowledgements.We thank the staffs of the various observatories (Palomar, Lick, Keck, etc.) where data for this study were obtained. The Oskar Klein Centre is funded by the Swedish Research Council. We gratefully acknowledge the support from the Knut and Alice Wallenberg Foundation. A.G.-Y. is supported by the EU/FP7 via ERC grant No. 307260, the Quantum Universe I-Core program by the Israeli Committee for Planning and Budgeting and the ISF; by Minerva and ISF grants; by the Weizmann-UK “making connections” program; and by Kimmel and ARCHES awards. A.V.F.’s research is supported by the Christopher R. Redlich Fund, the TABASGO Foundation, and NSF grant AST-1211916. D.C.L. and S.F.A. acknowledge support from NSF grants AST-1009571 and AST-1210311, under which part of this research (photometry data collected at MLO) was carried out. We thank Joseph Fedrow, Alyssa Del Rosario, Chuck Horst, and David Jaimes for assistance with the MLO observations. J.M.S. is supported by an NSF Astronomy and Astrophysics Postdoctoral Fellowship under award AST-1302771. This research has made use of the APASS database, located at the AAVSO web site; funding for APASS has been provided by the Robert Martin Ayers Sciences Fund. Research at Lick Observatory is partially supported by a generous gift from Google. Some of the data presented herein were obtained at the W. M. Keck Observatory, which is operated as a scientific partnership among the California Institute of Technology, the University of California, and NASA; the observatory was made possible by the generous financial support of the W. M. Keck Foundation. We acknowledge contributions to CCCP by S. B. Cenko, D. Fox, D. Sand, and A. Soderberg. We acknowledge M. Sullivan and K. Sharon for helping with the CCCP spectral observations. LANL participation in iPTF is supported by the US Department of Energy as part of the Laboratory Directed Research and Development program.
Note. – MW extinctions are obtained via NED (Fitzpatrick 1999; Schlafly & Finkbeiner 2011). Redshifts are from NED except for PTF09gpn and PTF12kso, where we measured it from the narrow lines of their host-galaxy spectra. Distances are redshift-independent for the three CCCP SNe, and obtained from the redshift for the three PTF SNe (assuming WMAP5 cosmological parameters; Komatsu et al. 2009). PTF12kso is also known as SN 2012gg, CSS J043611.05-050452.1, and LSQ12fjm (Sternberg et al. 2012).
||Last nondetection||Limit||Discovery||Discovery||Explosion date|
|2004ek||IAUC 8405||CBET 86||3242.5||19.5||3258.4730||17.1||3250.58.0|
|2004em||CBET 82||IAUC 8411||3214.5||19.0||3263.4680||17.5||3263.3|
|2005ci||CBET 168||CBET 169/Kleiser et al. (2011)||3523.81||19.0||3523.86||20.0||3523.69|
|PTF09gpn||PTF (Arcavi et al. 2010)||PTF (Arcavi et al. 2010)||5127.91 (P48)||21.00 (P48)||5137.77 (P48)||20.30 (P48)||5137.23|
|PTF12gcx||PTF (ATel 4293)||PTF (ATel 4293)||6085.39 (P48)||20.93 (P48)||6085.84 (P48)||21.14 (P48)||6085.72|
|PTF12kso||PTF / CBET 3296,1||PTF / CBET 3296.2||6153.00 (P48)||20.85 (P48,)||6198.98 (P48)||18.41 (P48)||6176.023.0|
Note. – Magnitudes are unfiltered or band except for the limiting magnitude of PTF12kso. For PTF12kso there is a nondetection at JD = 2,456,169.5 by LSQ (Aug. 30), which is later than the last PTF nondetection, but the limiting magnitude is unknown. LSQ discovered it at JD = 2,456,202.7324, later than PTF. For SN 2005ci the explosion date (3522.3) derived from the expanding photospheric method (EPM; see Jones et al. 2009) using the bands and the first two spectra matches that derived from the PL fit which we adopted. We notice that the explosion date of SN 2005ci occurs before the last nondetection, but the magnitude limit at this phase was not deep enough (19 mag) to detect SN 2005ci.
Note. – Besides CCCP photometry from Arcavi et al. (2012), the table includes photometry from Tsvetkov (2008) (marked by “*” and shifted by +0.1 (), +0.2(), +0.2(), and +0.1() mag with respect to the original values of Tsvetkov 2008 to match the CCCP photometry) of SN 2004ek and unfiltered photometry of SN 2005ci from Kleiser et al. (2011) (marked by “**”).
Note. – Photometric epochs obtained with P48 are marked with “*”, those obtained by LCOGT are marked with “”**”. For PTF12kso the epochs marked with “” correspond to photometry from MLO, which was standardized to the Johnson-Cousins filter system (no correction has been made to convert the obtained magnitudes to the system). The rest of the data were obtained by P60.
|08 Nov. 2009||5143.9||6.7||PTF09gpn||P200||DBSP||335710207|
|26 Dec. 2011||5921.5||…||PTF09gpn host galaxy||Keck-1||LRIS||358010198|
|06 Nov. 2012||6237.6||61.6||PTF12kso||Copernico 1.82m||AFOSC||46737386|
|07 Nov. 2012||6238.9||62.9||PTF12kso||P200||DBSP||32029946|
|13 Nov. 2012||6242.5||68.5||PTF12kso||P200||DBSP||36309782|
|09 Dec. 2012||6270.5||94.5||PTF12kso||Keck-2||DEIMOS||49027467|
|08 Jan. 2013||6301.8||124.8||PTF12kso||Keck-2||DEIMOS||45009632|
|27 Sep. 2015||7292.5||…||PTF12kso host galaxy||Keck-2||DEIMOS||46509500|
|16 Jul. 2012||6124.5||38.8||PTF12gcx||Keck-2||DEIMOS||44029593|
|26 Jul. 2012||6134.5||48.8||PTF12gcx||P200||DBSP||35029992|
|23 Aug. 2012||6162.7||77.0||PTF12gcx||Lick 3-m||Kast||345810332|
|15 Jun. 2007||4266.5||…||PTF12gcx host galaxy||SDSS 2.5-m||SDSS spectrograph||38019221|
|23 Sep. 2004||3271.5||21.0||SN 2004ek||P200||DBSP||41909600|
|24 Sep. 2004||3272.5||22.0||SN 2004ek||Lick 3-m||Kast||3320.0010500|
|19 Nov. 2004||3328.8||78.3||SN 2004ek||P200||DBSP||40009776|
|14 Dec. 2004||3353.5||103.0||SN 2004ek||Keck-1||LRIS||34509200|
|14 Jan. 2005||3385.1||133.6||SN 2004ek||P200||DBSP||33009900|
|24 Sep. 2004||3272.5||8.2||SN 2004em||Lick 3-m||Kast||3320.0010500|
|18 Oct. 2004||3296.5||32.2||SN 2004em||Keck-1||LRIS||34509200|
|19 Nov. 2004||3328.6||65.3||SN 2004em||P200||DBSP||35009780|
|05 Dec. 2004||3344.6||81.3||SN 2004em||P200||DBSP||35009796|
|05 Apr. 2005||3466.0||201.7||SN 2004em||P200||DBSP||36009991|
|16 Jul. 2005||3567.5||43.8||SN 2005ci||P200||DBSP||34009700|
|14 Aug. 2005||3596.7||73.0||SN 2005ci||P200||DBSP||35009801|
|08 Sep. 2005||3621.6||97.9||SN 2005ci||P200||DBSP||33009654|
|Band/SN||SN 2004ek||SN 2004em||SN 2005ci|
Note. – Rise times and peak magnitudes are obtained from low-order polynomial fits to the SN light curves.
Note. – Errors in the Ni mass are mainly due to the uncertainty in the distance and in the explosion epoch. We estimate the total uncertainty to be % for each SN. The Ni mass for PTF09gpn is determined by the -band light-curve comparison to SN 2009E at late epochs.
- The NASA/IPAC Extragalactic Database, http://ned.ipac.caltech.edu
IAUC 8405: Monard et al. (2004);
CBET 86: Modjaz et al. (2004);
CBET 82: Armstrong (2004);
IAUC 8411: Filippenko et al. (2004);
CBET 168: Madison & Li (2005);
CBET 169: Modjaz et al. (2005);
ATel 4293: Gal-Yam et al. (2012);
CBET 3296,1: Mahabal et al. (2012);
CBET 3296,2: Tomasella et al. (2012);
IAUC 8405: Monard et al. (2004);
CBET 86: Modjaz et al. (2004);
CBET 82: Armstrong (2004);
IAUC 8411: Filippenko et al. (2004);
CBET 168: Madison & Li (2005);
CBET 169: Modjaz et al. (2005);
ATel 4293: Gal-Yam et al. (2012);
CBET 3296,1: Mahabal et al. (2012);
CBET 3296,2: Tomasella et al. (2012);
- From explosion date.
- H flux is 2.4210 erg s cm.
- H flux is 6.7210 erg s cm.
- H flux is 1.3610 erg s cm.
- footnotetext: From Arnett et al. (1989).
- footnotetext: From the modeling of the bolometric light curve tail.
- footnotetext: From the hydrodynamical models.
- Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al. 2014, ApJS, 211, 17
- Amorín, R., Pérez-Montero, E., Contini, T., et al. 2015, A&A, 578, A105
- Anderson, J. P., González-Gaitán, S., Hamuy, M., et al. 2014, ApJ, 786, 67
- Arcavi, I., Gal-Yam, A., Cenko, S. B., et al. 2012, ApJ, 756, LL30
- Arcavi, I., Gal-Yam, A., Kasliwal, M. M., et al. 2010, ApJ, 721, 777
- Armstrong, M. 2004, Central Bureau Electronic Telegrams, 82, 1
- Arnett, W. D. 1979, ApJ, 230, L37
- Arnett, W. D., Bahcall, J. N., Kirshner, R. P., & Woosley, S. E. 1989, ARA&A, 27, 629
- Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
- Ben-Ami, S., Gal-Yam, A., Filippenko, A. V., et al. 2012, ApJ, 760, L33
- Blinnikov, S., Lundqvist, P., Bartunov, O., Nomoto, K., & Iwamoto, K. 2000, ApJ, 532, 1132
- Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, PASP, 125, 1031
- Cao, Y., Kasliwal, M. M., Arcavi, I., et al. 2013, ApJ, 775, LL7
- Cardamone, C., Schawinski, K., Sarzi, M., et al. 2009, MNRAS, 399, 1191
- Cenko, S. B., Fox, D. B., Moon, D.-S., et al. 2006, PASP, 118, 1396
- Chevalier, R. A. 1992, ApJ, 394, 599
- Chevalier, R. A., & Fransson, C. 2008, ApJ, 683, L135
- Cid Fernandes, R., Mateus, A., Sodré, L., Stasińska, G., & Gomes, J. M. 2005, MNRAS, 358, 363
- Dessart, L., & Hillier, D. J. 2005, A&A, 439, 671
- Elmhamdi, A., Danziger, I. J., Chugai, N., et al. 2003, MNRAS, 338, 939
- Faber, S. M., Phillips, A. C., Kibrick, R. I., et al. 2003, Proc. SPIE, 4841, 1657
- Filippenko, A. V. 1997, ARAA, 35, 307
- Filippenko, A. V., Ganeshalingam, M., & Swift, B. J. 2004, IAU Circ., 8411, 1
- Filippenko, A. V., Li, W. D., Treffers, R. R., & Modjaz, M. 2001, in Small-Telescope Astronomy on Global Scales, ed. W. P. Chen, C. Lemme, & B. Paczyński (San Francisco: ASP Vol. 246), 121
- Fitzpatrick, E. L. 1999, PASP, 111, 63
- Fremling, C., et al. 2016, in prep.
- Gal-Yam, A. 2012, Science, 337, 927
- Gal-Yam, A., Arcavi, I., Ofek, E. O., et al. 2014, Nature, 509, 471
- Gal-Yam, A., Ben-Ami, S., Arcavi, I., et al. 2012, The Astronomer’s Telegram, 4293, 1
- Gal-Yam, A., Cenko, S. B., Fox, D. B., et al. 2007, The Multicolored Landscape of Compact Objects and Their Explosive Origins, 924, 297
- Gall, E. E. E., Polshaw, J., Kotak, R., et al. 2015, arXiv:1502.06034
- Gilmozzi, R., Cassatella, A., Clavel, J., et al. 1987, Nature, 328, 318
- Gonzalez-Gaitan, S., Tominaga, N., Molina, J., et al. 2015, arXiv:1505.02988
- Hamuy, M. 2003, ApJ, 582, 905
- Inserra, C., Pastorello, A., Turatto, M., et al. 2013, A&A, 555, A142
- Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531
- Kewley, L. J., & Ellison, S. L. 2008, ApJ, 681, 1183
- Kleiser, I. K. W., Poznanski, D., Kasen, D., et al. 2011, MNRAS, 415, 372
- Komatsu, E., Dunkley, J., Nolta, M. R., et al. 2009, ApJS, 180, 330
- Kulkarni, S. R. 2013, The Astronomer’s Telegram, 4807, 1
- Jaskot, A. E., & Oey, M. S. 2013, ApJ, 766, 91
- Jones, M. I., Hamuy, M., Lira, P., et al. 2009, ApJ, 696, 1176
- Larsson, J., Fransson, C., Östlin, G., et al. 2011, Nature, 474, 484
- Leaman, J., Li, W., Chornock, R., & Filippenko, A. V. 2011, MNRAS, 412, 1419
- Leloudas, G., Schulze, S., Krühler, T., et al. 2015, MNRAS, 449, 917
- Li, W., et al. 2011a, MNRAS, 412, 1441
- Li, W., et al. 2011b, MNRAS, 412, 1473
- Lunnan, R., Chornock, R., Berger, E., et al. 2014, ApJ, 787, 138
- Lyman, J. D., Bersier, D., & James, P. A. 2014, MNRAS, 437, 3848
- Madison, D. R., & Li, W. 2005, Central Bureau Electronic Telegrams, 168, 1
- Mahabal, A. A., Drake, A. J., Djorgovski, S. G., et al. 2012, Central Bureau Electronic Telegrams, 3296, 1
- McCray, R. 1993, ARA&A, 31, 175
- McGaugh, S. S. 1991, ApJ, 380, 140
- Miller, J. S., & Stone, R. P. S. 1993, Lick Observatory Technical Report 66 (Santa Cruz, CA)
- Modjaz, M., Kirshner, R., Challis, P., & Calkins, M. 2005, Central Bureau Electronic Telegrams, 169, 1
- Modjaz, M., Kirshner, R., Challis, P., Matheson, T., & Calkins, M. 2004, Central Bureau Electronic Telegrams, 86, 1
- Monard, L. A. G., Boles, T., Puckett, T., & Cox, L. 2004, IAU Circ., 8405, 1
- Morozova, V., Piro, A. L., Renzo, M., et al. 2015, arXiv:1505.06746
- Nadyozhin, D. K. 2003, MNRAS, 346, 97
- Ofek, E. O., Laher, R., Law, N., et al. 2012, PASP, 124, 62
- Oke, J. B., Cohen, J. G., Carr, M., et al. 1995, PASP, 107, 375
- Oke, J. B., & Gunn, J. E. 1982, PASP, 94, 586
- Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei, 2nd ed. (Mill Valley, CA: University Science Books)
- Pastorello, A., Baron, E., Branch, D., et al. 2005, MNRAS, 360, 950
- Pastorello, A., Pumo, M. L., Navasardyan, H., et al. 2012, A&A, 537, A141
- Paxton, B., Bildsten, L., Dotter, A., et al. 2011, ApJS, 192, 3
- Pettini, M., & Pagel, B. E. J. 2004, MNRAS, 348, L59
- Phillips, M. M., Heathcote, S. R., Hamuy, M., & Navarrete, M. 1988, AJ, 95, 1087
- Piro, A. L., & Nakar, E. 2013, ApJ, 769, 67
- Podsiadlowski, P. 1992, PASP, 104, 717
- Quimby, R. M., Kulkarni, S. R., Kasliwal, M. M., et al. 2011, Nature, 474, 487
- Rabinak, I., & Waxman, E. 2011, ApJ, 728, 63
- Rahmer, G., Smith, R., Velur, V., et al. 2008, Proc. SPIE, 7014,
- Rau, A., Kulkarni, S. R., Law, N. M., et al. 2009, PASP, 121, 1334
- Rubin, A., et al. 2015, submitted
- Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103
- Smith, N., Mauerhan, J. C., Cenko, S. B., et al. 2015, MNRAS, 449, 1876
- Smith, C. E., & Nelson, B. 1969, PASP, 81, 74
- Sternberg, A., El-Hage, P., Elias-Rosa, N., et al. 2012, The Astronomer’s Telegram, 4493, 1
- Taddia, F., Sollerman, J., Fremling, C., et al. 2015, A&A, 580, A131
- Taddia, F., Sollerman, J., Leloudas, G., et al. 2015, A&A, 574, A60
- Taddia, F., Sollerman, J., Razza, A., et al. 2013, A&A, 558, A143
- Taddia, F., Stritzinger, M. D., Sollerman, J., et al. 2012, A&A, 537, A140
- Taubenberger, S., Pastorello, A., Mazzali, P. A., et al. 2006, MNRAS, 371, 1459
- Tomasella, L., Benetti, S., Cappellaro, E., et al. 2012, Central Bureau Electronic Telegrams, 3296, 2
- Tsvetkov, D. Y. 2008, Peremennye Zvezdy, 28, 3
- Utrobin, V. P. 2005, Astronomy Letters, 31, 806
- Utrobin, V. P., & Chugai, N. N. 2011, A&A, 532, A100
- Utrobin, V. P., Chugai, N. N., & Botticella, M. T. 2010, ApJ, 723, L89
- Weiss, A. 1989, ApJ, 339, 365
- Yaron, O., & Gal-Yam, A. 2012, PASP, 124, 668
- Young, T. R. 2004, ApJ, 617, 1233
- Young, T. R., & Branch, D. 1988, Nature, 333, 305