Letters

DISCOVERY AND EARLY MULTI-WAVELENGTH MEASUREMENTS OF THE ENERGETIC TYPE IC SUPERNOVA PTF12GZK: A MASSIVE-STAR EXPLOSION IN A DWARF HOST GALAXY

, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2012 November 9 © 2012. The American Astronomical Society. All rights reserved.
, , Citation Sagi Ben-Ami et al 2012 ApJL 760 L33 DOI 10.1088/2041-8205/760/2/L33

2041-8205/760/2/L33

ABSTRACT

We present the discovery and extensive early-time observations of the Type Ic supernova (SN) PTF12gzk. Our light curves show a rise of 0.8 mag within 2.5 hr. Power-law fits (f(t)∝(tt0)n) to these data constrain the explosion date to within one day. We cannot rule out a quadratic fireball model, but higher values of n are possible as well for larger areas in the fit parameter space. Our bolometric light curve and a dense spectral sequence are used to estimate the physical parameters of the exploding star and of the explosion. We show that the photometric evolution of PTF12gzk is slower than that of most SNe Ic. The high ejecta expansion velocities we measure (∼30, 000 km s−1 derived from line minima four days after explosion) are similar to the observed velocities of broad-lined SNe Ic associated with gamma-ray bursts (GRBs) rather than to normal SN Ic velocities. Yet, this SN does not show the persistent broad lines that are typical of broad-lined SNe Ic. The host-galaxy characteristics are also consistent with GRB-SN hosts, and not with normal SN Ic hosts. By comparison with the spectroscopically similar SN 2004aw, we suggest that the observed properties of PTF12gzk indicate an initial progenitor mass of 25–35 M and a large ((5–10) × 1051 erg) kinetic energy, the later being close to the regime of GRB-SN properties.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Type Ic CCSNe are heterogeneous (for a review on CCSNe, see Filippenko 1997). Their luminosity, ejected mass, and kinetic energy span over an order of magnitude, from the subluminous SN 2004aw to the overluminous SN 1998bw (Mazzali et al. 2009; Drout et al. 2011). The light-curve shapes of different events are also quite diverse. A subclass of SNe Ic whose spectra are characterized by broad lines (Type Ic-BL; prototype SN 1998bw) is the only one for which clear evidence of an association with gamma-ray bursts (GRBs) exists (GRB-SNe; see Woosley & Bloom 2006 for a review). Superluminous SNe (SLSNe) of Type Ic are even more powerful (Gal-Yam 2012 and references therein), but these probably result from a different physical mechanism. While SNe Ic are common in the center of high-metallicity galaxies (Anderson et al. 2012), SLSNe-I and broad-lined GRB-SNe tend to be found in dwarf hosts (e.g., Stanek et al. 2006; Modjaz et al. 2008; Arcavi et al. 2010), giving untargeted sky surveys an advantage over targeted surveys in detecting these types of cosmic explosions.

The Palomar Transient Factory (PTF; Law et al. 2009; Rau et al. 2009) is a wide-field untargeted sky survey which explores the transient optical sky. It uses the PTF CFH12k camera mounted on the Palomar 48 inch telescope (P48). PTF's short observing cadence and real-time capability (e.g., Gal-Yam et al. 2011) enables the discovery and study of SNe at early stages of the explosion. In this Letter we report the discovery and study of PTF12gzk, a peculiar SN Ic in a dwarf star-forming galaxy located at redshift z = 0.0137 (distance 57.8 Mpc, distance modulus 33.8 mag, assuming H0 = 71 km s−1 Mpc−1).

2. DISCOVERY

PTF12gzk was discovered on 2012 July 24.3 (UT dates are used herein) at α(J2000) = 22h12m41fs53 and δ(J2000) = +00°30'43farcs1, in the dwarf galaxy SDSS J221241.53+003042.7 (within the SDSS Stripe 82 footprint).32 The discovery magnitude was 20.66 in the r band,33 and it was not detected down to mag 21.6 (3σ) in previous PTF images obtained on July 19 (Ben-Ami et al. 2012).

Shortly after discovery, we initiated an extensive follow-up campaign, including our Hubble Space Telescope (HST) Target-of-Opportunity (ToO) program for STIS ultraviolet (UV) spectroscopy of a stripped-envelope SN (Cycle 19, GO-12530; PI Filippenko) and Swift X-ray and UV photometry (Cycle 8, PID 8110099; PI Kasliwal), the results of which are presented herein (see Figures 1 and 2 for photometry and spectroscopy, respectively). We also triggered radio and millimeter observations using the JVLA (program 12A-363; PI Horesh) and the Combined Array for Research in Millimeter-wave Astronomy (program 12A-c0945; PI Horesh); see A. Horesh et al. (in preparation).

Figure 1.

Figure 1. PTF12gzk photometry. Top: the long rise time and the large delay in peak magnitudes between the B and r bands is evident in the light curve. Inset: a bolometric light curve derived from a fourth-degree polynomial fit to our UV-optical-IR data. Bottom left: a χ2 probability-density-function contour plot comparing power-law models with the first 2.5 hr of observations (see the text). The value of the index n is not well constrained. Discontinuities in the contour curves are caused by finite sampling of the parameter space, and have no statistical significance. Bottom right: early photometry. The flux is scaled with respect to the first data point. Fits of the scaled flux to f∝(tt0)2 (red) and f∝(tt0)5 (blue) are also shown. Considering the flux normalization as an additional parameter, there is no significant difference between the two fits.

Standard image High-resolution image
Figure 2.

Figure 2. Top: PTF12gzk optical spectra; telluric features are marked. Bottom: PTF12gzk host galaxy image/spectrum taken from SDSS. Note the blue color and small spatial size (1'' = 280 pc), as well as the strong emission lines and blue continuum. The August 1 spectrum was taken with the FLOYDS spectrograph mounted on the FTN, a low-dispersion spectrograph (R ≈ 400) with a single, fixed grating and a cross-dispersing prism, placing the first- and second-order spectra onto the CCD with a single-exposure wavelength coverage of 3200–10,000 Å. While the identical twin FLOYDS spectrograph on FTN and FTS will ultimately be robotically operated, the present spectrum was taken during hardware commissioning under human control.

Standard image High-resolution image

3. OBSERVATIONS

3.1. Photometry

Optical photometry of PTF12gzk was obtained using multiple telescopes (Table 1). All data were calibrated with respect to the Sloan Digital Sky Survey (SDSS) catalog. Light curves of PTF12gzk are shown in Figure 1.

Table 1. PTF12gzk Photometry

Date (JD) Instrument Filter Apparent Magnitude 1σ Uncertainty
2456132.763 P48+CFH12k R 20.66 0.12
2456132.779 P48+CFH12k R 20.24 0.15
2456132.808 P48+CFH12k R 20.13 0.08
2456132.809 P48+CFH12k R 20.11 0.12
2456132.848 P48+CFH12k R 20.02 0.07

Note. The full set of photometric data is available in the online journal, as well as from WISeREP (http://www.weizmann.ac.il/astrophysics/wiserep/; Yaron & Gal-Yam 2012).

Only a portion of this table is shown here to demonstrate its form and content. A machine-readable version of the full table is available.

Download table as:  DataTypeset image

The optical data were reduced using standard IRAF procedures for aperture photometry (FTS data were reduced using point-spread function photometry via DoPHOT). We subtract reference templates from the P48 and P60 data to remove contamination from the host.34 Pre-explosion templates were not used for other data, but the contribution from the underlying galaxy (g, r, z = 19.05, 19.03, 18.75 mag from SDSS) is negligible, as was also demonstrated by reducing P60 data both with and without template subtraction. The data were calibrated to SDSS stars in the field, using the transformation equations given by Jordi et al. (2006) to place the local standards on the Johnson–Cousins system.

Infrared (IR) photometry of PTF12gzk (Table 1) was obtained using the Wide Field Camera mounted on the United Kingdom Infrared Telescope (UKIRT-WFCAM) using SExtractor, and calibrated with respect to the 2MASS catalog (magnitude errors <0.07 mag) using the relation of Hodgkin et al. (2009).

Based on the absence of strong Na i D lines (see Section 4), as well as the blue early-time spectrum (see below), we take a conservative approach assuming negligible extinction by the host galaxy at the SN location.

During the first night of observation, PTF12gzk brightened by ∼0.8 mag in less than 2.5 hr to 19.85 mag in r; we thus obtained remarkably early coverage of an SN Ic. We find that the SN peaked at r = 15.2 mag on August 14, g = 15.55 on August 8, and B = 16 on August 4. On August 15, PTF12gzk peaked in the i band.

PTF12gzk was observed with the X-Ray Telescope (XRT) and the Ultraviolet/Optical Telescope (UVOT) on board the Swift satellite. XRT measurements, beginning at 13:39 on July 31, detected no source at the location of PTF12gzk; we estimate a dead-time-corrected limit on the XRT count rate of <2 × 10−3 cps. Assuming a power-law spectrum with a photon index of 2, this corresponds to a limit on the X-ray flux of <7 × 10−14 erg cm−2 s−1. UVOT data were reduced using a 5'' aperture to measure the counts for the coincidence-loss correction, a 3'' or 5'' source aperture (based on the signal-to-noise ratio) for the aperture photometry, and applying an aperture correction as appropriate. We used zero points from Breeveld et al. (2011), including a time-dependent sensitivity loss to put the magnitudes on the UVOT photometric system. Results are given in Figure 1.

We have calculated a bolometric light curve by integrating the flux in the UVW2, UVM2, UVW1, B, g, V, r, R, i, I, J, H, K filters and fitting a fourth-degree polynomial to the obtained data points. When lacking IR photometry, we assume a constant fraction of IR flux, found to be ∼18% from synthetic photometry using TSPEC IR spectra (see Section 4.2). We estimate that the lack of IR (UV) coverage prior to August 7 introduces an uncertainty of ∼10% from the small variation seen in the IR contribution (<5% with respect to the overall flux) between August 4 (the first TSPEC IR spectrum) and August 12. Uncertainties introduced by the lack of UV photometry are ∼5% from the even smaller variation in the UV contribution (<2% relative to the overall flux). The bolometric light curve is given in Figure 1 and shows a rise time of 18 ± 1 day, similar to that in the r band.

3.2. Spectroscopy

Extensive spectroscopy of PTF12gzk was performed, and detailed analysis will be presented in a future paper. A selection of optical spectra is shown in Figure 2.35

The first spectrum was obtained on July 27, revealing the object to be a peculiar SN I with a blue continuum, reminiscent of the blue continuum observed in the early-time spectra of SN 2006aj (Mazzali et al. 2006b). Initially it was difficult to classify the spectrum securely, but later optical spectra resemble those of SNe Ic, with an absence of prominent He i lines (Ben-Ami et al. 2012).

Classification with the help of SNID (Blondin & Tonry 2007) of the August 9 Lick 3 m/Kast spectrum suggests that the best fit is to the peculiar SN Ic SN 2004aw near peak brightness (Taubenberger et al. 2006). A good match is also obtained for PTF10qts, a Type Ic-BL SN, at ∼30 days after explosion. All spectra show prominent, broad absorption lines of Ca ii, Si ii, Fe ii, and O i, with SYNOW fits yielding maximum velocities of ∼35, 000 km s−1 for the July 27 spectrum, and ∼20, 000 km s−1 for the August 9 spectrum. See Figure 3 and Section 5 for further discussion.

Figure 3.

Figure 3. Left: spectroscopic analysis; top to bottom. The DBSP spectrum taken on July 27. An HST STIS/MAMA and STIS/CCD (G430 and G750) spectrum taken on August 7. A Lick/Kast spectrum taken on August 9. The continuum is consistent with blackbody temperatures of 11,000 K and 8000 K in the July 27 and August 9 spectra, respectively. A SYNOW fit shows that the July 27 spectrum is dominated by Ca ii, Ti II, and Fe ii lines at velocities of 35,000 km s−1, and by O i, Mg ii, and Si ii at 28,000 km s−1; the corresponding velocities for the August 9 spectrum are 21,000 and 18,000 km s−1. The best match suggested by SNID (Blondin & Tonry 2007) is to a spectrum of SN 2004aw near peak brightness (blue curve). A good match is also obtained for PTF10qts, a Type Ic-BL SN, at ∼30 days after explosion. Top right: the Palomar 5 m TSPEC IR spectrum obtained on August 7. Bottom right: the HST STIS/MAMA UV spectrum taken on August 7. The flux deficit with respect to a blackbody spectrum indicates strong line blanketing, evidence for a highly mixed outer envelope.

Standard image High-resolution image

The absence of the 2.1 μm He i feature in the IR TSPEC spectrum taken on August 7 (Figure 3, bottom-right panel) suggests that He is not abundant in the envelope (Hachinger et al. 2012). While a strong absorption line is present at 1 μm, which could be identified with He i 1.0830 μm, in this case it is likely contaminated by other elements such as C, Mg, S, and Ca (e.g., Mazzali& Lucy 1998; Sauer et al. 2006; Hachinger et al. 2012).

The first HST UV spectrum, obtained by the Space Telescope Imaging Spectrometer with the near-UV Multi-Anode Microchannel Array (STIS/MAMA), was taken between 20:06 August 6 and 03:35 August 7, with a total exposure time of 11,278 s. We also obtained optical spectra with the STIS/CCD, using grisms centered at 4300 and 7500 Å. A mosaic of all three HST spectra is given in Figure 3, with a focus on the UV spectrum in the bottom-right panel. The flux deficit with respect to blackbody spectra seen in the STIS/MAMA UV spectrum is characteristic of all SNe I (Panagia 2007) and indicates strong line blanketing, evidence for a highly mixed outer envelope devoid of hydrogen. Full analysis of the HST UV spectra will be presented by S. Ben-Ami et al. (in preparation).

4. HOST-GALAXY ANALYSIS

An image and a spectrum of the host galaxy obtained by the SDSS prior to explosion are shown in Figure 2. The SN exploded very close to the center of the host galaxy (offset 0farcs0 N and 0farcs4 E). We downloaded the host-galaxy spectrum from the SDSS 9th Data Release (DR9; Ahn et al. 2012). After correcting for the host-galaxy redshift, we measure integrated emission-line fluxes using standard procedures via splot in IRAF, and follow Pérez-Montero & Diaz (2003) to compute statistical errors.36 Using the SDSS Petrosian magnitudes, correcting for Galactic and host-galaxy extinction (see below), and applying K-corrections via kcorrect (v4.2, Blanton & Roweis 2007) at the host-galaxy redshift, we derive an absolute magnitude of Mr = −14.8 mag. This is an extremely underluminous dwarf galaxy, well below the luminosities of the Large and Small Magellanic Clouds. It is one of the least luminous host galaxies of CCSNe discovered by PTF, and it is not a common host galaxy for an SN Ic, even considering those found via untargeted surveys (Arcavi et al. 2010), except for SNe Ic-BL. In addition, the clear detections of numerous bright emission lines allow us to apply standard diagnostics of the star formation rate (SFR), extinction, and metallicity to this galaxy.

4.1. Star Formation Properties

Using the observed intensity ratio of Hα/Hβ, and assuming the Case B recombination value of 2.86 and the Cardelli et al. (1989) extinction law with RV = 3.1, we estimate a reddening of E(BV) = 0.14 ± 0.05 mag. The values of the intensity ratios [N ii]/Hα and [O iii]/Hβ indicate that the prominent emission lines are due to recent vigorous star formation rather than to an active galactic nucleus (Baldwin et al. 1981). Furthermore, we do not detect any absorption features or a Balmer break that may indicate the presence of an older stellar population.

From the measured rest-frame equivalent width of Hβ, EW = 23.3 ± 1 Å, we derive an age of 6.4 ± 0.1 Myr for the young stellar population using the relations of Levesque et al. (2010, Section 3.4 and Table 4) and the measured metallicity of the host (for Z = 0.004; see below). From the measured integrated Hα emission-line flux, corrected for extinction, we obtain L(Hα) = (16 ± 2) × 1038 erg s−1, which translates to SFR(Hα) =0.013 ± 0.002 M yr−1 using the conversion from Kennicutt (1998). Since the SDSS spectrum was obtained through a 3'' diameter fiber encompassing an area of 0.55 kpc2 at the host-galaxy redshift, we measure a normalized SFR of 0.023 ± 0.003 M yr−1 kpc−2. Since the Petrosian radius (from DR9) is 3farcs4, the fraction of global host-galaxy light entering the fiber, which is also at the SN position, is 20%.

4.2. Oxygen Abundance

For computing the oxygen abundance (Modjaz et al. 2011 and references therein), we correct the detected emission-line fluxes of [O iii], [N ii], Hβ, and Hα for reddening, and employ the scales of Pettini & Pagel (2004, PP04-O3N2) and of Kewley & Dopita (2002, KD02) to obtain values of 12 + log(O/H)PP04-O3N2 = 8.12+0.04−0.02 and 12 + log(O/H)KD02 = 8.13+0.05−0.04, respectively. We conclude that the metallicity of the host galaxy is 0.2–0.3 Z, having used the solar oxygen abundance of 12 + log(O/H) = 8.69 (Asplund et al. 2009).

The oxygen abundance of the PTF12gzk host is well below that of the hosts of normal SNe Ic found via untargeted surveys as presented by Modjaz et al. (2011; mean 12 + log(O/H)PP04 = 8.7 ± 0.1) and Sanders et al. (2012, 12 + log(O/H)PP04 = 8.61 ± 0.2). Indeed, it is much closer to that of the hosts of SNe Ic-BL and GRB-SNe (Modjaz et al. 2008; Sanders et al. 2012), and of SLSNe-I (Young et al. 2010; Stoll et al. 2011).

5. DISCUSSION

PTF12gzk is a luminous SN Ic, at the high end of the SN Ic luminosity distribution (Drout et al. 2011). It exhibits a slow rise of 18 days to its peak r-band magnitude, with B peaking ∼10 days earlier. This is a large gap relative to other SN I, though similar to SN 2004aw (Taubenberger et al. 2006); it is caused by metal-line absorption from heavy elements in the outer layers of the ejecta, as is evident from spectra taken after August 1.

A least-squares fit to a f(t)∝(tt0)n behavior of our well-sampled early photometry places the explosion date between 10 and 40 hr prior to our discovery at the 95% confidence level.37 We cannot rule out the popular quadratic fireball model, but higher values of n are possible as well for larger areas in the fit parameter space (Figure 1, bottom panels).

Spectroscopically, PTF12gzk exhibits high expansion velocities, ∼30, 000 km s−1 (Si ii absorption velocity). Other SNe Ic with similar velocities are broad-lined SNe Ic (Figure 4), some of which are associated with GRBs (Woosley & Bloom 2006 and references within), while no such association was determined for PTF12gzk (see also SN 2009bb; Soderberg et al. 2010). Most similar is SN 2003lw, an SN associated with a GRB (Mazzali et al. 2006a). A possible explanation is a burst misaligned with our line of sight, or a failed GRB. Such a scenario is further supported by the host-galaxy characteristics, resembling those of a broad-lined SN Ic host galaxy. We know of no typical SN Ic exploding in a host with similar luminosity and oxygen abundance. The observed relatively narrow lines give a dispersion of Δv/v ≈ 0.25, compared to ∼1 in the case of broad-lined SNe Ic, and may suggest a nonspherical explosion geometry (Leonard et al. 2006), or that the ejecta mass is high or has a very steep density gradient (a discussion on the effects of a steep density gradient on the light curve (LC) can be found in Piro & Nakar 2012). Late-time, nebular spectra will probe the geometry of the explosion in more detail.

Figure 4.

Figure 4. PTF12gzk has characteristic velocities of a broad-lined SN Ic. All SNe above the dashed line besides PTF12gzk are GRB-SNe (Xs) or Type Ic-BL with no GRB association (circles), while those below it are normal SNe Ic (squares). Velocities are obtained through modeling of the spectrum or through direct measurements of the Si ii 6355  Å line (SNe 2010bh, 2002ap, 2003jd, 2004aw, and 2007gr).

Standard image High-resolution image

From the Si ii line velocity at peak brightness of PTF12gzk (15,300 km s−1 from the August 12 spectrum) and SN 2004aw (12,400 km s−1; J. Deng et al., in preparation), and the rise time of these two SNe, we use the following scaling relations (Arnett 1982; Mazzali et al. 2009; see also P. A. Mazzali et al., in preparation) to estimate the physical properties of PTF12gzk: τ ≈ κ1/2M3/4E−1/4 and v = (2E/M)1/2, where τ is the light-curve rise time, E is the kinetic energy, and κ is the opacity. The derived ejecta mass is 7.5 M (6–12 M), pointing to a large initial progenitor mass of 25–35 M, though the latter values are highly uncertain (Mazzali et al. 2000). We derive a kinetic energy of 7.5 × 1051 erg [(5–10) × 1051  erg]. Using Lmax and tmax we get an estimated 56Ni mass of 0.37 M (scaling the PTF12gzk light curve to that of SN 2003dh; Mazzali et al. 2003). Using the V-band peak magnitude versus nickel mass relation presented by Perets et al. (2010), we derive a 56Ni mass of 0.35 M, in agreement with the results derived from the scaling relations. These physical properties, as well as the high expansion velocities and the host galaxy, are unlike those of normal SNe Ic, which typically occur in large hosts and have low ejecta masses, and kinetic energies, (1.7 M, and 1051 erg Drout et al. 2011. For nickel mass in Type Ic SNe see for example Taubenberger et al. 2006, and Sauer et al. 2006). Instead, they are reminiscent of GRB-SNe (Mazzali et al. 2009).

PTF12gzk is a remarkable example of an SN Ic in terms of expansion velocities, evolution timescale, the ejected mass, and the kinetic energy released in the explosion. We conclude that these properties point to the explosion of a massive star deficient in H and He, at the higher-mass end of SN Ic progenitors. This further illustrates the peculiar population of SNe Ic exploding in dwarf hosts (Arcavi et al. 2010), as seen also in the case of GRB-SNe and most SLSNe-I.

PTF12gzk demonstrates the advantages of using an untargeted sky survey such as PTF with an extensive network of instruments and telescopes in various wavebands to detect and rapidly characterize unusual cases of cosmic explosions.

S.B. is supported by a Ramon Fellowship from ISA. A.G. acknowledges support by grants from the ISF, BSF, GIF, Minerva, and the EU FP7/ERC. A.V.F. and his group benefit from financial assistance from Gary & Cynthia Bengier, the Richard & Rhoda Goldman Fund, the Sylvia & Jim Katzman Foundation, the Christopher R. Redlich Fund, the TABASGO Foundation, NSF grants AST-0908886 and AST-1211916, and NASA/HST grant GO-12530 from STScI (which is operated by the AURA, Inc., under NASA contract NAS 05-26555). P.A.M. and E.P. acknowledge financial support from grants INAF PRIN 2011 and ASI/INAF I/088/06/0. M.I. and Y.J. were supported by the Creative Initiative program of the NRFK. M.M.K. acknowledges Hubble and Carnegie-Princeton Fellowships. D.C.L. is supported by NSF grant AST-1009571. E.O.O. acknowledges the Arye Dissentshik career development chair and a grant from the Israeli MOST.

PTF is a collaboration of Caltech, LCOGT, the Weizmann Institute, LBNL/NERSC, Oxford, Columbia, IPAC, and UC Berkeley. We are grateful for the assistance of the staff at the various observatories used to obtain data. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by JPL, Caltech, under contract with NASA.

Footnotes

  • 32 

    It was also independently discovered by the La-Silla Quest (LSQ 2012dwl) and Pan-Starrs1 (PS1-12baa) surveys.

  • 33 

    PTF magnitudes are given in the PTF natural-magnitude system (Ofek et al. 2012), with respect to the SDSS r-band magnitudes.

  • 34 

    SDSS images were used as reference images for the P60 data.

  • 35 

    All ground-based spectra presented in this Letter are released on WISeREP (http://www.weizmann.ac.il/astrophysics/wiserep/; Yaron & Gal-Yam 2012).

  • 36 

    The measured Hγ, Hβ, Hα, [O iii] λ4959,[O iii] λ5007, and [N ii] λ6584 emission-line fluxes are 34.1 ± 3.3, 88.5 ± 3.6, 290.9 ± 3.4, 81.5 ± 3.2, 232.7 ± 4.8, and 10.3 ± 1.0 × 10−17 erg s−1 cm−2 before extinction correction, respectively. Errors include statistical measurement uncertainties.

  • 37 

    The fireball models can at best give an underestimate on the explosion date, since they do not incorporate the photon diffusion time.

Please wait… references are loading.
10.1088/2041-8205/760/2/L33