A&A 457, 857-864 (2006)
DOI: 10.1051/0004-6361:20065530

The GRB 060218/SN 2006aj event in the context of other gamma-ray burst supernovae[*],[*]

P. Ferrero1 - D. A. Kann1 - A. Zeh2 - S. Klose1 - E. Pian2 - E. Palazzi3 - N. Masetti3 - D. H. Hartmann4 - J. Sollerman5 - J. Deng6 - A. V. Filippenko7 - J. Greiner8 - M. A. Hughes9 - P. Mazzali2,10 - W. Li7 - E. Rol11 - R. J. Smith12 - N. R. Tanvir9,11


1 - Thüringer Landessternwarte Tautenburg, 07778 Tautenburg, Germany
2 - Istituto Nazionale di Astrofisica-OATs, 34131 Trieste, Italy
3 - INAF, Istituto di Astrofisica Spaziale e Fisica Cosmica, Sez. di Bologna, 40129 Bologna, Italy
4 - Clemson University, Department of Physics and Astronomy, Clemson, SC 29634-0978, USA
5 - Dark Cosmology Center, Niels Bohr Institute, Copenhagen University, 2100 Copenhagen, Denmark
6 - National Astronomical Observatories, CAS, Chaoyang District, Beijing 100012, PR China
7 - Department of Astronomy, University of California, Berkeley, CA 94720-3411, USA
8 - Max-Planck-Institut für extraterrestische Physik, 85741 Garching, Germany
9 - Centre for Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, AL10 9AB, UK
10 - Max-Planck Institut für Astrophysik, 85748 Garching, Germany
11 - Department of Physics and Astronomy, University of Leicester, Leicester, LE1 7RH, UK
12 - Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Birkenhead, CH41 1LD, UK

Received 2 May 2006 / Accepted 13 July 2006

Abstract
The supernova SN 2006aj associated with GRB 060218 is the second-closest GRB-SN observed to date (z=0.033). We present Very Large Telescope, Liverpool Telescope, and Katzman Automatic Imaging Telescope multi-color photometry of SN 2006aj. This supernova is found to be subluminous and rapidly evolving. Its early light curve includes an additional wavelength-dependent component, which can be interpreted as shock break-out. We compare the photometric evolution of multi-band light curves with the corresponding properties of the present sample of more than 10 GRB-SNe with precisely known redshifts. Using host-galaxy extinction measurements, we derive extinction-corrected GRB-SN luminosities and place SN 2006aj in the context of this GRB-selected supernova sample as well as in the context of local stripped-envelope supernovae.

Key words: gamma rays: bursts: X-rays: individuals: GRB 060218 - supernovae: individual: SN 2006aj

1 Introduction

There is now substantial evidence for the association of long-duration GRBs with core-collapse supernovae (SNe of type Ic to be more specific; see Filippenko 1997 for a review of supernova spectral classification). The first hint for such a relation came with the contemporaneous discovery of GRB 980425 and of a local SN (1998bw) in its error circle (Galama et al. 1998). Subsequently, Bloom et al. (1999) reported a possible underlying SN component in the afterglow of GRB 980326. The first spectroscopic identification of a SN Ic superposed on a GRB afterglow (GRB 030329/SN 2003dh; Hjorth et al. 2003; Matheson et al. 2003; Stanek et al. 2003) solidified this association considerably.

The photometric database on GRB-SNe was summarized by Zeh et al. (2004,2005, henceforth Z04, Z05) who, by studying a sample of 21 (Z04) and 29 (Z05) bursts with well established redshifts in a systematic way, found or confirmed a weak SN component for 9 (Z04) and 13 (Z05) sources in those sets, respectively. In particular, all long bursts with redshifts $z\la 0.7$ were found to have a late-time bump in their optical afterglows, suggesting that in fact all long bursts are physically associated with SN explosions. Additional photometric and partly spectroscopic results on recent GRB-SNe with well-known redshift were presented by Malesani et al. (2004; GRB 031203/SN 2003lw), Soderberg et al. (2006) and Stanek et al. (2005) on SNe associated with GRB 040924 and GRB 041006, as well as Della Valle et al. (2006; GRB 050525A/SN 2005nc). A recent review of the status of the supernova - GRB connection is presented by Woosley & Bloom (2006).

GRB 980425/SN 1998bw is still the closest GRB-SN to date, and poses the question whether it is a typical representative of this class, a "standard source'' in some sense. Its comparatively low isotropic energy release in the gamma-ray band (less than 1048 erg, Galama et al. 1998) makes this burst stand out from the cosmological GRB ensemble with known redshifts (cf. Ghirlanda et al. 2004). It is therefore of great interest to find additional events in the nearby universe. Consequently, the discovery of GRB 060218 (Cusumano et al. 2006) (more precisely an X-Ray Flash, XRF) and its identification with an energetic SN by Masetti et al. (2006), now designated SN 2006aj (Soderberg et al. 2006), attracted much attention (Modjaz et al. 2006a; Sollerman et al. 2006; Mirabal et al. 2006; Cobb et al. 2006). Its small distance (z=0.033, Mirabal & Halpern 2006) results in an isotropic equivalent energy of the burst of $5 \times 10^{49}$ erg in the 15.5-154.8 keV band, which again is underenergetic in comparison to the cosmological GRB sample (Modjaz et al. 2006a), whose mean redshift is larger than z=1 (Jakobsson et al. 2006; Bagoly et al. 2006).

Here, we report results from a photometric ESO Very Large Telescope (VLT) campaign on GRB 060218/SN 2006aj covering a time span of nearly 4 weeks, with additional data gathered with the robotic Liverpool Telescope (LT) at La Palma and the robotic Katzman Automatic Imaging Telescope (KAIT; Filippenko et al. 2001; Li et al. 2003a) at Lick Observatory. Results from the spectroscopic VLT campaign are presented in Pian et al. (2006) and modeled in Mazzali et al. (2006b). While previous publications on SN 2006aj discussed the X-ray and optical light curves in detail, we place the photometric evolution of SN 2006aj in the context of all presently known GRB-SNe with precise redshifts. In particular, we update the fits reported in Z04 and Z05, add several more recent GRB-SNe, and correct the photometry for host-galaxy extinction as derived in Kann et al. (2006, henceforth K06).

2 Observations and data reduction

2.1 VLT data

We observed SN 2006aj both spectroscopically and photometrically with the ESO Very Large Telescope FORS1 and FORS2 instruments. Observations were performed until 26 days after the burst (Table 1). After this point in time, the SN location was no longer observable due to high airmass. For the photometry the exposure time was between 30 and 60 seconds. The images, in the Bessell B, V, R, and I filters, were bias-subtracted and flat-fielded with standard reduction procedures using IRAF[*] and final photometry was performed with standard Point Spread Function (PSF) fitting using the DAOPHOT II image data analysis package (Stetson 1987) within MIDAS[*]. Photometric calibration of the images was performed using the standard star 685 of the field SA98 (Landolt 1992), imaged during the same night of our first observational run, which was used to create 16 secondary standards in the field around the SN that were subsequently used for all individual images to draw a linear regression among instrumental and reference magnitudes and calculate the value of the optical transient. The accuracy of the photometry was confirmed by considering the zero point values of several stars in the calibration field in all bands for different nights. In addition, we investigated a possible influence of color terms. We found that they could affect the photometry with an additional error of up to 0.02 mag, which was then added in quadrature to the individual photometric measurement errors.

2.2 Liverpool Telescope data and KAIT data

Further photometric data were obtained with the 2 m robotic Liverpool Telescope on La Palma on several occasions over a period of about 2.5 weeks post-burst. Observations were made in Bessell B, V and SDSS $r^{\prime}$ filters and zero-point reduced by comparison to two field stars (star S1 and star 4 of Modjaz et al. 2006a). The magnitudes were transformed into Johnson B, V and Cousins R-band via equations previously derived for the camera (Steele 2004). The quoted errors reflect the scatter in individual measurements for subexposures at each epoch, and include an estimate of the calibration uncertainty.

Additional data were also obtained with the 0.8 m robotic KAIT at Lick Observatory on four consecutive nights between 3 to 6 days post-burst. Unfiltered observations were obtained in all four nights, and Bessell B, V, R, and I filtered observations were made in the last three nights. The images were bias-subtracted and flat-fielded with standard reduction procedures using IRAF, and the final photometry was performed with PSF-fitting technique in IRAF/DAOPHOT. The secondary standards calibrated with VLT data were used for the final photometric calibration. The unfiltered magnitudes were converted to the R band following the method described in Li et al. (2003b).


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{5530fig1.eps}
\end{figure} Figure 1: The location of SN 2006aj is shown here in a VLT/FORS2 R-band image taken on 2006 March 1, at 00:52 UT (30 s exposure time). The field of view is approximately 2 $\times $ 3 arcmin.
Open with DEXTER

2.3 The magnitudes of the host galaxy

The host galaxy of GRB 060218/SN 2006aj was imaged pre-burst by the Sloan Digital Sky Survey (Cool et al. 2006; Adelman-McCarthy et al. 2006). It was pointed out (Modjaz et al. 2006a; Hicken et al. 2006) that there are offsets between the SDSS calibration and later efforts, rendering the transformed host-galaxy magnitudes of Cool et al. (2006) too bright. The values used for the host galaxy in this paper are $B=20.57\pm0.07$, $V=20.18\pm0.04$, $R=20.03\pm0.03$ and $I=19.58\pm0.06$, not corrected for Galactic extinction. These magnitudes have been evaluated from those reported by Cool et al. (2006), by applying the corrections recommended by Modjaz et al. (2006a) and Hicken et al. (2006). Our $\it BVRI$ magnitudes have been obtained from Cool et al.'s u*b*v*r*i* magnitudes by adopting the conversion of Lupton (2005)[*] and taking into account that the offsets determined by Modjaz et al. (2006) in the $\it BVR$ bands are 0.40, 0.27, 0.20 mag, respectively[*]. In order to determine the offset in the I band, we have considered that the offset determined by Hicken et al. (2006) and Modjaz et al. (2006a) in the i*-filter is 0.15 mag, and we have used the SDSS transformations[*] to derive from this the offset in the i-band. We note that our host-galaxy magnitudes are consistent to those of Sollerman et al. (2006) within the errors. As the host galaxy and SN 2006aj were not separated at the resolution of the images, the host magnitudes have to be subtracted.

2.4 The extinction correction

For the extinction correction, we used the values derived by Guenther et al. (2006) from an analysis of the Na I D2 absorption lines in our Galaxy ( $A_{\rm V}=0.39$ mag) and the host galaxy ( $A_{\rm V}=0.13$ mag) in a VLT UVES spectrum, obtained during this campaign close to SN maximum. We note that this host-galaxy extinction is typical of the lines of sight toward GRBs (K06). Using these extinction and host magnitude values, we derive the magnitudes of the pure SN light component.

3 Results

We followed previous works (Z04, Z05) and modeled the light curve of SN2006aj using SN 1998bw as a template (Galama et al. 1998) that was shifted[*] to the corresponding redshift and scaled in luminosity (in the SN  rest frame) by a factor k and in time evolution by a factor s, while zero host extinction was assumed for SN 1998bw (Patat et al. 2001). In doing so, we find an additional component visible in the early data that makes the light curve systematically brighter than that of SN 1998bw. This component has also been noted by Mirabal et al. (2006). Furthermore, strong spectral evolution is noticeable, with the flux excess being stronger toward shorter wavelengths. The most reasonable explanation for this additional blue component is the light due to shock break-out through a dense progenitor wind (Campana et al. 2006). However, we caution that at least part of this additional component may be the result of an intrinsically different (as compared to SN 1998bw) early SN light curve. SN 2002ap (cf. Appendix) is also overluminous at early times in comparison to SN 1998bw, and has been compared, both photometrically and spectroscopically, with SN 2006aj (e.g. Pian et al. 2006). But the early multi-color evolution of SN 2002ap is achromatic, showing no evidence for strong spectral evolution as seen in SN 2006aj.

We thus used two different methods to fit the SN light curve. In the first fit, we excluded all data earlier than 8.8 days, but fit only the SN 1998bw model light curve to the data, without any additional contributions. The result is shown in Fig. 2. The early excess blue flux is clearly visible in the residuals. The data after 8.8 days are matched very well by the SN 1998bw model. Note that for typical GRB-SNe we also fit only the peak and decaying part of the SN, as the early flux is typically dominated by the GRB afterglow.

In the second fit, we assumed that an additional component decays according to a power-law, with the decay index $\alpha $ ( $F_{\nu}\propto
t^{-\alpha}$) being a free parameter of the fit, and independent for each band (Fig. 3). The derived values of the luminosity ratio k and the stretch factor s as well as the decay index $\alpha $ are given in Table 2. It is apparent that while the results derived from the two methods are similar, they do in fact differ by as much as 10% (in k). The additional component improves the quality of the fit in all bands. The steeper decay at shorter wavelengths is also in agreement with the Swift UVOT photometry of the blue shock break-out (Campana et al. 2006). The strong wavelength dependence of the decay (Table 2) argues against this component being the optical afterglow of the GRB. On the other hand, the very shallow decay in $\it VRI$ indicates that the time evolution of the SED of this additional component is probably more complex than a simple decay according to a pure power-law, as assumed in our analysis.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{5530fig2.eps}
\end{figure} Figure 2: The light curves of SN 2006aj based on our $\it BVRI$ data after correcting for extinction and host flux contribution. The data were fitted using the light curves of SN 1998bw as a template, not considering data from SN 2006aj that was taken before 8.8 days post burst. The residuals $\Delta $m represent observed values minus the fit.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{5530fig3.eps}
\end{figure} Figure 3: The light curves of SN 2006aj fitted by including an additional early component which decays according to a power-law. Solid lines mark the fit, the dashed lines are the power-law component and the dotted lines are the SN component. The time scale is logarithmic to show in detail the earlier data.
Open with DEXTER

Table 2: Fitting results for the light curves of SN 2006aj attained by either excluding early data (top) or including an additional power-law component decaying with index $\alpha $ (see text for details).

Table 3: Luminosity ratio k and stretch factor s for GRB-SNe with known redshift normalized to SN 1998bw before and after correcting for host extinction. In cases where no extinction could be derived, only the non-corrected values are given. In all cases the SN fits were derived in the observer R-band frame. This list is complete up to the end of 2005, with the exception of GRB 040924 (Soderberg et al. 2006). These results supersede those presented in Z05. Differences compared to Z05 are basically due to the inclusion of new observational data or a revision of the used value for the Galactic extinction along the line of sight. Furthermore, note that we use different cosmological parameters than those in Z04 and Z05.

The small value of the stretch factor s means that SN 2006aj evolves much faster than SN 1998bw (see also Pian et al. 2006; Modjaz et al. 2006a; Sollerman et al. 2006; Mirabal et al. 2006; Cobb et al. 2006; Soderberg et al. 2006). Furthermore, SN 2006aj has a slightly different color than SN 1998bw, being brighter in the B and I band than in the V and R band. The k values we obtain are also in accordance with the values found by the authors cited above. In all bands, SN 2006aj is about 30% less luminous than SN 1998bw. From the light curve data, independent of any fit, we find the following peak absolute magnitudes of the SN: $M_B=-18.29\pm0.05$, $M_{V}=-18.76\pm0.05$, $M_R=-18.89\pm0.05$, $M_I=-18.95\pm0.10$.

In order to place SN 2006aj in the context of the presently known GRB-SNe, we started with the results presented in Z04 and Z05, updated them in cases where new data were available, and corrected the k value using host-galaxy extinctions derived in K06, including uncertainties. We added two further GRB-SNe to the sample. The result for XRF 020903 is based on data from Bersier et al. (2006). We found that the SN of this event is also superimposed on an underlying power-law afterglow component. Unfortunately, the afterglow data are too sparse to derive any conclusions on the extinction in the host galaxy. For GRB 031203, we used the extinction derived by Mazzali et al. (2006a). In the case of GRB 050525A, we employed the data presented by Della Valle et al. (2006) and applied the extinction found by Blustin et al. (2006). In those cases where no extinction could be found in K06, we only present the observed values and note that these are lower limits on the luminosity of the SNe (Table 3). On the other hand, the line of sight extinction is often low and the true values are thus not expected to be much higher.

4 Discussion

4.1 SN 2006aj in the context of known GRB-SNe

SN 2006aj is a rare nearby GRB-SN, with a distance of only about 140 Mpc. Its luminosity is about 70% of that of SN 1998bw, confirming that GRB-selected supernovae may not qualify for the label "standard candle''. Given the UVES data of SN 2006aj, taken close to SN maximum light (Guenther et al. 2006; Wiersema et al. 2006), it is unlikely that the deficiency in luminosity is the result of dust extinction in the host galaxy of SN 2006aj.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{5530fig4.eps}
\end{figure} Figure 4: The distribution of the luminosity factor k in the R band with redshift z. There is no apparent correlation with redshift, even though we caution that all the k values refer to the R band in the observer frame, corresponding to different wavelengths in the host frame. Filled symbols mark the k values that have been corrected for host extinction, open symbols represent supernovae for which this correction could not be applied. A star marks the value derived from the R-band light curve of SN 2006aj by removing data before 8.8 days from the fit. SN 1998bw is at k=1, by definition.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=16cm,clip]{5530fig5.eps}
\end{figure} Figure 5: Luminosity factor k versus stretch factor s in the R band for all GRB-SNe in our sample (Table 3). The symbols are identical to those used in Fig. 4. SN 1998bw defines k=1, s=1. Note that in the case of SN 2006aj the values for $\it BVRI$ are plotted. Furthermore, we plot the type Ic supernova SN 1994I and the type Ic broad-lined SN 2002ap, neither of which are associated with GRBs. Both are fainter than any GRB-SN in our sample for which the host extinction has been corrected.
Open with DEXTER

By the end of 2005, there were 14 optical afterglows[*] with known redshifts that showed evidence for extra light in the R band at late times (GRBs 970228, 990712, 991208, 000911, 010921, 011121, 020405, 021211, 030329, 031203, 040924[*], 041006, 050525A, XRF 020903). For eight afterglows with a known SN bump (GRBs 991208, 000911, 010921, 011121, 020405, 021211, 030329, 041006) K06 were able to derive a host extinction value. We used these to determine the extinction-corrected luminosities of these SNe in the observer frame (Fig. 4). It turns out that only two (030329, 031203) or perhaps three (including 991208) of the 14 SNe were actually more luminous than SN 1998bw with high significance (but see, e.g., Deng et al. 2005, for a spectroscopic modeling of SN 2003dh that results in k<1). Remarkably, most extinction corrected SNe cluster around 0.6<k<1.5. We find that the four SNe not corrected for extinction (GRBs 970228, 990712, 020903, 021211) are typically fainter, implying that a correction for host galaxy extinction will probably also shift them into this range. While one might worry about the fact that for each individual GRB-SN, due to the different redshifts, the k factor refers to a different wavelength region, SN 2006aj does not contradict our assumption that k is not strongly dependent on wavelength (Fig. 5).

We conclude that the present data, even though some have large uncertainties, indicate that the width of the GRB-SN luminosity function is at least 2 mag, comparable to what is known for the other types of SNe (cf. Richardson et al. 2002; Richardson et al. 2006, henceforth RBB06). In particular, there is no evidence that the luminosity function evolves with redshift: The width of the luminosity function for z<0.2 is comparable to the width at $z\approx 0.7$. It is also much narrower than the distribution of intrinsic afterglow luminosities (K06).

4.2 SN 2006aj in the context of local stripped-envelope supernovae

It is also interesting to compare the light curve properties of SN 2006aj and other GRB-SNe with well-studied local stripped-envelope (i.e. types Ib, Ic, IIb) supernovae, as is shown in Fig. 6. Distance moduli and absolute magnitudes for these SNe were taken from RBB06 (who use H0=60 km s-1 Mpc-1) and have been transformed to the world model used here (H0=71 km s-1 Mpc-1) by adding $-5 \log 71/60 =-
0.365$ mag to the former and +0.365 mag to the latter values. Similarly, for the world model used here SN 1998bw is 0.19 mag less luminous than given in Galama et al. (1998), who use H0=65 km s-1 Mpc-1, i.e. k=1 corresponds to $M_{\rm V} = - 19.16$. The absolute visual magnitudes of the GRB-SNe in our sample were then calculated according to $M_{\rm V}^{\rm GRB-SN} - M_{\rm V}^{\rm 98bw} = - 2.5 \log k$, assuming that k is independent of wavelength. Figure 6 shows the result obtained in this way. Most notable is that the ensemble of GRB-SNe is at the bright end of the luminosities of local stripped-envelope supernovae. In other words, the present data indicate that type Ic supernovae with associated (detected) GRBs are on average more luminous in the optical bands than those without detected GRBs. In particular, SN 2006aj is no exception from this rule. We caution, however, that our assumption of wavelength-independence of k could be an oversimplification. The larger the redshift, the more uncertain is the absolute V-band magnitude of a GRB-SN derived in this way. This uncertainty is not included in the error bars plotted in Fig. 6.

Finally, we note that an observational bias might affect the interpretation of Fig. 6: the larger the redshift, the more difficult it is to observe less luminous stripped-envelope supernovae.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{5530fig6.eps}
\end{figure} Figure 6: The absolute V-band magnitude $M_{\rm V}$ of stripped-envelope supernovae versus their distance modulus. This is an expanded version of Figure 1 in RBB06, and the values for SN of type IIb, Ib and Ic have been taken from that work but transformed according to the world model used here. The individual GRB-SNe have the same numbers as in Fig. 7. We have also labeled two local type Ic SNe we have included in our study (Fig. 5; Appendix). GRB-SNe without known host extinction are not included. The dashed line at $M_{\rm V}=-19.5$ denotes the typical absolute magnitude of type Ia SNe (the "ridgeline''). The slanted dotted line denotes a constant visual magnitude $m_{\rm V}=16$ in the case of no extinction.
Open with DEXTER

RBB06 give host extinction values for all 27 events in their sample. Therefore, we investigated whether host extinction of GRB-SNe is different from those of local stripped-envelope supernovae. No substantial differences are apparent (Fig. 7). While it is interesting that SN 2006aj is less affected by host extinction compared to the local sample of type Ic SNe, the current sample is too small to draw reliable conclusions from this finding.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{5530fig7.eps}
\end{figure} Figure 7: The visual host-galaxy extinction $A_{\rm V}$ (in the host frame) of stripped-envelope supernovae versus their distance modulus. The values for SNe of type IIb, Ib and Ic have been taken from RBB06, while the data for the GRB-SNe are from K06. We labeled individual GRB-SNe and the two type Ic SNe we have also studied. GRB-SNe without known host extinction were not included.
Open with DEXTER

4.3 The stretch factor

Finally, we consider the statistics of the stretch factor s, which is in some sense more reliable since it is not affected by the extinction issue. First of all, among all known GRB-SNe, SN 2006aj has the smallest s value, i.e. it is the fastest GRB-SN ever seen (Pian et al. 2006; Mazzali et al. 2006b; Modjaz et al. 2006a). Interestingly, about half of the GRB-SNe have s<1, i.e., they are evolving faster than SN 1998bw. In contrast, the evolution of the SN associated with GRB 041006 was very slow (Fig. 5) with high significance (Stanek et al. 2005). Basically only two SNe (000911, 970228) occupy the (k<1, s>1) region, i.e., they are slow and subluminous, contrary to the general trend, but both have large error bars. Stanek et al. (2005) noted that for GRB-SNe a relation may exist between light curve shape and luminosity, similar to the one established for type Ia SNe (see Fig. 3 in Stanek et al. 2005). Figure 5 neither supports nor contradicts the existence of such a potential relation. While a trend of rising k with rising s seems to be visible, a fit to the data does not support this trend with reasonable statistical significance.

4.4 Shortcomings of the procedure

The procedure applied here has shortcomings that are basically related to cosmological effects on one hand and to data quality on the other.

As long as one is only concerned with broad-band photometry, which is the case in this study, the observed light curves of GRB-SNe usually refer to different wavelength bands in their host frames. However for SN 1998bw it has been shown that the light curve shape, the time of the peak flux, and the peak flux itself are a function of wavelength (as is the case for other SN types, too). Unfortunately, the available photometric data base is in most cases restricted to R-band observations (in the observer frame). In order to be able to compare light curves of GRB-SNe that occurred at various redshifts with light curves of SN 1998bw, the simplest and in some sense only useful approach is then to assume that the luminosity ratio k and the stretch factor s, which are both normalized to SN 1998bw, are independent of wavelength. In other words, we assume that the SED of all GRB-SNe (in their host frame) is the same at all times.

In addition, the data base used for our light curve fits is usually weighted to data obtained past the peak time of a GRB-SN under consideration. The reason is that before the peak time usually the afterglow dominates the light of the optical transient, while after the peak time usually the SN light dominates. This problem is difficult to overcome. Therefore, for basically all cosmologically remote SNe we have no information about the details of the SN light curves at early times up to several days after the corresponding burst. In particular the stretch factor s is then mainly affected by the late-time behaviour of the SN light curves.

It is therefore by no means clear if a relation between s and k indeed exists. It is possible that such a relation is hidden by the relatively large error bars for individual (k,s) values of the GRB-SNe on the one hand (Fig. 5) and by the basic assumptions that went into the procedure we applied here on the other hand. More spectral data from GRB-SNe might finally solve this issue. However, progress made in this regard during the last years was only modest, at best. Most Swift detected GRBs are at such a high redshift that no SN spectroscopy can be performed with current telescopes within a reasonable amount of observing time. Therefore, the photometric approach utilized here derives some of its justification from these spectroscopic limitations.

5 Conclusions

SN 2006aj is the fastest evolving and one of the least luminous GRB-SN discovered so far, being only about 70% as luminous as SN 1998bw. This places it at the faint end of the GRB-SN luminosity distribution which so far, after extinction correction, covers the range from about 0.6 to 2 times the luminosity of SN 1998bw. Placing SN 2006aj in the context of the luminosities of local Type Ic supernovae without detected GRB (RBB06) reveals, however, that SN 2006aj is still at the bright end of their luminosity distribution. SN 2006aj thus follows the general "rule'' that GRB-SNe tend to be more luminous than (local) Type Ic SNe without detected GRB.

GRBs 980425 and 060218 suggest that the prompt emission properties are not correlated with the optical properties of the associated SNe. GRB 980425 had a duration of 31 s (Pian et al. 1999), while GRB 060218 lasted for more than 2000 s (Campana et al. 2006), and the corresponding isotropic equivalent energy release in the gamma-ray band was different by a factor of 10. The corresponding SN luminosities, however, differ by only 30%. Assuming that differences in instrumental sensitivities (BeppoSAX versus Swift) do not play a role here, this supports the notion that the properties of the GRB and the associated SN are to a large extent independent of each other.

Acknowledgements

We thank the staff at ESO/Paranal, in particular A. Bik, C. Dumas, O. Hainaut, E. Jehin, C. Ledoux, O. Marco, C. Melo, D. Naef, L. Schmidtobreick for performing the observations. F.P., D.A.K. and S.K. acknowledge support by DFG grant Kl 766/13-2 and by the German Academic Exchange Service (DAAD) under grant No. D/05/54048. A.V.F. is supported by NASA/Swift grants NNG05GF35G and NNG06GI86G. KAIT and its ongoing research were made possible by generous donations from Sun Microsystems, Inc., the Hewlett-Packard Company, AutoScope Corporation, Lick Observatory, the National Science Foundation, the University of California, the Sylvia & Jim Katzman Foundation, and the TABASGO Foundation. This work has benefitted from collaboration within the EU FP6 Research Training Network HPRN-CT-2002-00294 "Gamma-Ray Bursts: an Enigma and a Tool.'' We thank the anonymous referee for a prompt and constructive report.

Appendix A: SN 1994I and SN 2002ap

Following the suggestion of the referee, we consider here two local type Ic supernovae that were not physically related to a detected GRB but are often used for comparison with GRB-SNe, and place them in the k,s context used in the present paper. While the fit appears to be worse, it should be noted that when fitting a GRB-SN in all cases the fit is practically always weighted to data points lying after the light curve maximum when the SN light dominates the afterglow light.

SN 1994I in the Whirlpool Galaxy M 51 (Puckett et al. 1994) counts as the "prototypical'' type Ic supernova (Filippenko et al. 1995), although it is faster than any other (RBB06). For the analysis of the light curve we used R-band data from Yokoo et al. (1994), Schmidt & Kirschner (1994), Lee et al. (1995) and Richmond et al. (1996) (with late data taken only from Schmidt & Kirschner and Richmond et al.). SN 1994I occurred near the bulge of M 51 and was highly reddened: Richmond et al. estimate $E_{(B-V)}=0.45\pm0.16$ with a Milky Way extinction curve. The explosion date of the supernova is not well known. From light curve modeling, Iwamoto et al. (1994) derive an explosion date of JD 2449433. These authors equate this with a B-band maximum 12 days after the explosion. A comparison with the data from Richmond et al. (1996) shows, however, that JD 2449437 corresponds to 12 days before the B maximum. Using the latter, the early light curve becomes strongly subluminous in comparison to SN 1998bw, has to be excluded from the fit and $s=0.69\pm0.01$ is found from the location of the maximum alone. On the other hand, if we take March 31.0 UT (JD 2 449 442.5) as the explosion date, we find that up to an age of 25 days the early light curve is well fitted by the SN 1998bw light curve. Assuming a distance modulus of 29.60 mag (Richmond et al. 1996), corresponding to a distance d=8.34 Mpc for the world model used here, we find, after correcting for the host-galaxy extinction (Richmond et al. 1996), k=0.39+0.17-0.12. Thus SN 1994I is fainter than any GRB-SN when we account for host extinction.


  \begin{figure}
\par\includegraphics[width=8.2cm,clip]{5530fig8.eps}
\end{figure} Figure A.1: The two local type Ic supernovae SN 1994I and SN 2002ap fitted with the model light curve of SN 1998bw at the respective redshifts. In both cases, the light curve of SN 1998bw does not fit the light curves of the other SNe completely. The best agreement is obtained for SN 2002ap, where data after ten days up to 60 days are fitted well by the SN 1998bw light curve. For SN 1994I data before and after the maximum had to be excluded from the fit.

SN 2002ap, which exploded in the outer spiral arms of M 74 (Nakano et al. 2002), is the closest and best followed-up broad-lined SN Ic to date. A search through data of all IPN satellites did not reveal any GRB that could be associated with the SN (Hurley et al. 2002); still, the explosion date was well constrained, as the supernova was discovered very rapidly (Mazzali et al. 2002). For the analysis of the light curve we collected data from the following sources: Gal-Yam et al. (2002), Borisov et al. (2002), Cook et al. (2002), Nishihara et al. (2002), Yoshii et al. (2003), Pandey et al. (2003), Doroshenko et al. (2003), Foley et al. (2003), and Vinkó et al. (2004). In addition we made use of the data collected at SNWEB[*]. Based on these data SN 2002ap is overluminous at early times in comparison to SN 1998bw. The host-galaxy extinction is very small, $E_{(B-V)}=0.020\pm0.0008$ (Takada-Hidai et al. 2002). Correcting for this and the Galactic extinction, and assuming a distance modulus of 29.32 mag (Foley et al. 2003), which corresponds to a distance d=7.3 Mpc for the cosmological model used here, we find $k=0.219\pm0.0002$ and $s=0.747\pm0.0006$. Thus, the stretch factor s for SN 2002ap is very close to what we find for SN 2006aj, but the SN is much fainter.

References

 

  
Online Material

Table 1: The photometry of SN 2006aj from the Very Large Telescope (VLT), the Liverpool Telescope (LT) and the Katzman Automatic Imaging Telescope (KAIT), including the acquisition images for the VLT spectroscopy (Pian et al. 2006). The date is the UT exposure mid-time. Column t-t0 refers to the time in days after the burst trigger at t0 = 3:34:30 UT of Feb. 18, 2006 (Cusumano et al. 2006). The third column gives measured magnitudes which are corrected neither for extinction nor flux from the host galaxy. The fourth column gives the pure supernova magnitudes after correcting for Galactic extinction, host-galaxy flux and host-galaxy extinction. Exposure times are in seconds.



Copyright ESO 2006