A&A 462, L57-L60 (2007)
DOI: 10.1051/0004-6361:20066660

LETTER TO THE EDITOR

Extensive multiband study of the X-ray rich GRB 050408[*]

A likely off-axis event with an intense energy injection[*]

A. de Ugarte Postigo1 - T. A. Fatkhullin2 - G. Jóhannesson3 - J. Gorosabel1 - V. V. Sokolov2 - A. J. Castro-Tirado1 - Yu. Yu. Balega2 - O. I. Spiridonova2 - M. Jelínek1 - S. Guziy1,4 - D. Pérez-Ramírez5 - J. Hjorth6 - P. Laursen6 - D. Bersier7 - S. B. Pandey1,8 - M. Bremer9 - A. Monfardini7 - K. Y. Huang10 - Y. Urata11,12 - W. H. Ip10 - T. Tamagawa11 - D. Kinoshita12 - T. Mizuno13 - Y. Arai13 - H. Yamagishi13 - T. Soyano14 - F. Usui15 - M. Tashiro16 - K. Abe16 - K. Onda16 - Z. Aslan17,18 - I. Khamitov17 - T. Ozisik17 - U. Kiziloglu19 - I. Bikmaev20,21 - N. Sakhibullin20,21 - R. Burenin22 - M. Pavlinsky22 - R. Sunyaev22 - D. Bhattacharya23 - A. P. Kamble23 - C. H. Ishwara Chandra24 - S. A. Trushkin2


1 - Instituto de Astrofísica de Andalucía (IAA-CSIC), Apartado de Correos 3004, 18080 Granada, Spain
2 - Special Astrophysical Observatory, Nizhnij Arkhyz, Zelenchokskaya, Karachaevo-Cherkesia, 369167 Russia
3 - Science Institute, University of Iceland, Dunhaga 3, 107 Reykjavík, Iceland
4 - Nikolaev State University, Nikolska 24, 54030, Nikolaev, Ukraine
5 - Dpto. de Física (EPS), Universidad de Jaén, 23071 Jaén, Spain
6 - Dark Cosmology Centre, Niels Bohr Institute, University of Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen, Denmark
7 - Astrophysics Research Inst., Liverpool John Moores Univ., Twelve Quays House, Egerton Wharf, Birkenhead, CH41 1LD, UK
8 - The UCL Mullard Space Science Laboratory, Holmbury St. Mary, Dorking, Surrey, RH5 6NT, UK
9 - Institut de Radio Astronomie Millimétrique (IRAM), 300 rue de la Piscine, 38406 Saint-Martin d'Hères, France
10 - Institute of Astronomy, National Central University, Chung-Li 32054, Taiwan, PR China
11 - Institute for Physics and Chemical Research (RIKEN), Wako, Saitama 351-0198, Japan
12 - Tokyo Institute of Technology, Ookayama, Meguro, Tokyo 152-8550, Japan
13 - Department of Astronomy and Earth Sciences, Tokyo Gakugei University, Koganei, Tokyo 184, Japan
14 - Kiso Observatory, Institute of Astronomy, University of Tokyo, Mitake-mura, Kiso-gun, Nagano 397-0101, Japan
15 - Japan Aerospace Exploration Agency, Institute of Space and Astronautical Science, Sagamihara, Kanagawa 229-8510, Japan
16 - Saitama University, Sakura-ku, Saitama 338-8570, Japan
17 - TÜBITAK National Observatory, Akdeniz Üniversitesi, 07058, Antalya, Turkey
18 - Akdeniz University, Physics Department, 07058 Antalya, Turkey
19 - Middle East Technical University, Physics Department, Inonu Bulvari, Ankara, 06531, Turkey
20 - Departments of Astronomy, Kazan State University, Kremlevskaya Str., 18, Kazan, 420008, Russia
21 - Academy of Sciences of Tatarstan, Bauman Str., 20, Kazan, 420111, Russia
22 - Space Research Institute (IKI), 84/32 Profsoyuznaya, Moscow, 117997, Russia
23 - Raman Research Institute, Bangalore 560 080, India
24 - National Centre for Radio Astrophysics, Ganeshkhind, Pune, 411007, India

Received 30 October 2006 / Accepted 14 December 2006

Abstract
Aims. Understand the shape and implications of the multiband ligth curve of GRB 050408, an X-ray rich (XRR) burst.
Methods. We present a multiband optical light curve, covering the time from the onset of the $\gamma$-ray event to several months after, when we only detect the host galaxy. Together with X-ray, millimetre and radio observations we compile what, to our knowledge, is the most complete multiband coverage of an XRR burst afterglow to date.
Results. The optical and X-ray light curve is characterised by an early flattening and an intense bump peaking around 6 days after the burst onset. We explain the former by an off-axis viewed jet, in agreement with the predictions made for XRR by some models, and the latter with an energy injection equivalent in intensity to the initial shock. The analysis of the spectral flux distribution reveals an extinction compatible with a low chemical enrichment surrounding the burst. Together with the detection of an underlying starburst host galaxy we can strengthen the link between XRR and classical long-duration bursts.

Key words: gamma rays: bursts - techniques: photometric

1 Introduction

X-ray flashes (XRFs) where first identified by Beppo-SAX (Heise et al. 2001) as those bursts detected by the X-ray camera but not the $\gamma$-ray monitor. Later studies based on the larger sample gathered by HETE-2 (Sakamoto et al. 2005b) gave a more general (and instrument- independent) classification and confirmed the intermediate group of events, the X-ray rich (XRR) class, previously detected by Ginga (Yoshida & Murakami 1994) and Granat/WATCH (Castro-Tirado et al. 1994). It is now known that long-duration GRBs (LGRBs), XRRs and XRFs share the same isotropic distribution in the sky, the same duration range and similar spectrum, with the main difference of having respectively lower observed spectral peak energy  $E_{\rm peak}^{\rm obs}$ in the $\nu F_{\nu}$ spectrum. They seem to form a continuum and thereby, most of the proposed models have tried to explain them as a unified phenomena (see a summary of the different models in Granot et al. 2005).

GRB 050408 was detected by WXM, SXC, and FREGATE aboard HETE-2 (Sakamoto et al. 2005a) at 2005 April 08 16:22:50.93 UT (t0 hereafter). With an observed peak energy of $\sim$20 keV it was classified as an XRR event. The 1.0 m Zeiss and 6.0 m BTA telescopes at the Special Astrophysical Observatory (SAO) in Russia pointed at the position delivered by HETE-2 through the GCN (GRB Circular Network) and detected the optical afterglow (de Ugarte Postigo et al. 2005a) coincident with the possition of the X-ray afterglow detected by Swift/XRT, which began observing 42 min after the burst (Wells et al. 2005). The precise localisation allowed further optical (Covino et al. 2007) and spectroscopic (Berger et al. 2005; Prochaska et al. 2005) observations, this latter ones detemined a redshift of z = 1.236.

In Sect. 2 we present the observations and the reduction methods that have been used for the analysis of the data. Section 3 describes the results that have been obtained, including observations of the host galaxy and modelling of the light curve. Section 4 discusses the implications of the analysis of the light curve.

2 Observations and data reduction

For this work we have compiled over 60 photometric measurements in U, B, V, Rc and Ic bands from 12 telescopes. The images where reduced using standard techniques based on IRAF[*] and JIBARO (de Ugarte Postigo et al. 2005b).

The burst happened during night time in Japan, where a fast follow up was carried out. The very wide field camera, WIDGET was monitoring the field of view of HETE-2 when the event was reported but detected no optical emission before, during or after the gamma-ray emission down to an unfiltered limiting magnitude of 9.7 (all limits given throughout the paper are 3-$\sigma $). The 1.05 m KISO Schmidt telescope pointed to the error box 20 min after the burst but failed to detect the afterglow. Finally, the 1 m LOT telescope observed the field 55 min after the burst, images that later served to confirm the afterglow (Huang et al. 2005b). The discovery of the afterglow was made with the data of the 1.0 m Zeiss and 6.0 m BTA telescopes in Russia (de Ugarte Postigo et al. 2005a), starting 115 min after the burst, when observations became possible from that site.

Further observations were performed from the 1.5 m Russian-Turkish telescope, in TÜBITAK National Observatory, the 4.0 m Blanco telescope in Cerro Tololo, the 4.0 m Mayall telescope in Kitt Peak, the 2.0 m Faulkes Telescope North (FTN) in Haleakala and the 3.5 m Telescopio Nazionale Galileo (TNG) in la Palma. A specially intense multiband campaign was carried out from the 1.54 m Danish telescope in La Silla, where daily observations were obtained during the first 8 days following the burst onset.

Finally, 8 months after the burst, deep observations were made from the 3.5 m telescope at Calar Alto. In these images we detect the host galaxy of the burst in B and Rc bands and impose a limit in Ic band.

Optical photometric calibration is based on the observation of several standard fields (Landolt 1992) using the 1.54 m Danish telescope at La Silla and the 1.5 m telescope at Sierra Nevada Observatory. From these observations we derive 12 secondary standards of different brightnesses. A log with the observations and the calibration stars are given as online material.

Our dataset is completed with several millimetre and radio limits. 6 epochs of millimetre observations were carried out with the 6-antenna Plateau de Bure interferometer (PdB, Guilloteau et al. 1992). No detection was obtained in either of the 1 mm or 3 mm bands, although a 3-$\sigma $ signal was found on the phase center in both observing bands on April 18. Careful re-analysis of the data did not reveal these signals as instrumental artifacts. Based on the extreme spectral slope and the non-detection on April 19, we conclude that this result is either due to a statistical fluctuation or an unusual event of interstellar scattering at high galactic latitude, and not due to a source-intrinsic variation. Data calibration was done using the GILDAS software package[*] using MWC349 as primary flux calibrator and 3C 273 as amplitude and phase calibrator. Further observations where obtained 13 days after the burst at 1.28 GHz from GMRT and at 8.4 GHz from RATAN-600.

3 Results

3.1 Light curve

In order to put together all the radio, millimetre, near infrared (nIR), optical and X-ray data available (including data from Foley et al. 2006 and Soderberg 2005a,b), we have determined the corresponding flux density values for all observations. X-ray afterglow counts, obtained from Foley et al. (2006) have been converted and corrected for hydrogen column extinction using WEBPIMMS[*] taking as spectral model a powerlaw with a slope of $\beta_{\rm X}$ = 1.16 and a column density of $N_{\rm H}$ = 0.25 $\times $ $10^{22}~{\rm cm}^{-1}$ (Nousek et al. 2006; Chincarini et al. 2005). The optical data have been corrected for galactic reddening (using an E(B-V) = 0.026, Schlegel et al. 1998) and intrinsic extinction (see Sect. 3.2). The measured/estimated host galaxy flux has been subtracted from the data to obtain the clean afterglow flux (see Sect. 3.3). The conversion of the optical data to flux density was done using the transformations given by Fukugita et al. (1995) for the optical and by Allen et al. (2000) for the nIR. The resulting light curves are shown in Fig. 1. Note the intense bump, rising at $\sim$3 days and peaking at $\sim$6 days, both in optical and X-rays. These kind of fluctuations have already been detected in the light curves of LGRBs and short-duration bursts (SGRBs) (de Ugarte Postigo et al. 2005c, 2006).

3.2 Study of the optical-nIR SFD

We have constructed the $\it UBVR{\rm c}I{\rm c}JHK$-band Spectral Flux distribution (SFD) of the afterglow 0.6 days after the burst, when near simultaneous optical and nIR observations were available. The $\it UBVR{\rm c}I{\rm c}$-band magnitudes from this work were complemented with the $\it JHK$-band values reported by Foley et al. (2006). Synchronisation to a unique timing is done by assuming a powerlaw with an index of $\alpha$ = 0.7 (F $\propto$ $t^{-\alpha}$), as derived from a linear fit of the nearby multiband data of the afterglow.


  \begin{figure}
\par\includegraphics[width=7.75cm,clip]{6660fig1.eps}\end{figure} Figure 1: Multiwavelength light curve of the GRB 050408 afterglow in the observer frame. The lines show the best fit of a fireball model with one energy injection (at 2.9 days) seen off-axis (see text for details). Our observations are plotted with filled symbols, while the ones obtained from the literature are represented by empty ones, this convention is used for all the figures.
Open with DEXTER


  \begin{figure}
\par\includegraphics[width=8cm,clip]{6660fig2.eps}\end{figure} Figure 2: Spectral flux distribution of the afterglow 0.6 days after the burst onset in the observer frame. The different lines represent results from fitting the data to various extinction laws: Small Magellanic Cloud (SMC), Large Magellanic Cloud (LMC), Milky Way (MW) and No Extinction (NE).
Open with DEXTER

The fluxes are used for fitting an extincted powerlaw ( $F_{\nu} \propto
10^{-0.4A_{\nu}}\nu^{-\beta}$) with 3 different extinction laws: Milky Way (MW), Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) as described by Pei (1992). This allows us to obtain $A_{\rm V}$ and $\beta$ simultaneously. The results of these 3 fits are complemented with an unextincted powerlaw case (NE), see Fig. 2 and Table 1. The best fit to the SFD of the afterglow is obtained when considering a SMC extinction law ( $\chi^{2}/{\rm d.o.f.}$ = 5.0/5). This is consistent with what has been previously found for other LGRB afterglows (Kann et al. 2006).

Table 1: Results of the SFD fitting at 0.6 days for different extinction laws.

3.3 The host galaxy

Several months after the gamma-ray event we revisited the GRB field with the 3.5 m telescope at Calar Alto Observatory in order to search for the host galaxy. Images were obtained in B, Rc and Ic bands, yielding a faint detection in B and Rc and imposing a limit on Ic. We derive galaxy colour indices of $(B-R{\rm c}) = 0.7$ $\pm$ 0.5 and $(R{\rm c}-I{\rm c}) \le 0.73$. These values are corrected for Galactic extinction. We have compared these values to the ones derived from the templates computed by Kinney et al. (1996) for a wide variety of galaxy types. We may conclude that only starburst galaxy templates are consistent with them. The best correlation is obtained with the starbust 2 template, with an intrinsic extinction of E(B-V) = 0.16.

3.4 Modelling of the multiband data

Using the model and methods described by Jóhannesson et al. (2006) we fitted the multiband observations of the afterglow (galaxy subtracted) to a fireball model with energy injections, viewed both on-axis and off-axis (with varying viewing angles). At least one injection is needed in order to account for the bump seen at 6 days which would carry as much energy as the initial shock. Another characteristic of the light curve is a flattening of the early light curve, seen in Rc and X-rays during the first hours of the burst, which has already been reported by Foley et al. (2006). This can be explained either by an early energy injection (single or continous), or an outflow with a low initial Lorentz factor, or as the result of an off-axis viewed burst. Similar early behaviour has already been found in other bursts (Nousek et al. 2006; Zhang et al. 2006).

Our preferred scenario (giving the best fit) describes the burst as a collimated ( $\theta_0=2.7^\circ$) fireball seen off axis ( $\theta_v=1.45 \theta_0$) expanding into a uniform low density environment (n0=0.01 cm-3) with an electron index p=2.03 and having an additional energy injection after 2.9 days with 1.2 times the initial energy. No further injections are needed to explain the light curve with the available amount of data. From the fit we obtain a  $\chi^2/{\rm d.o.f.} = 158.1/93$. The jet break, as defined by Sari et al. (1999), would be expected initially at 1.6 days, or at 3.8 days due to the energy injection. However, due to the effect of the equal arrival time surface, it is further delayed to approximatelly 30 days. Figure 3 shows the radio to X-ray SFD predicted by our model for 3 epochs, together with observational data, host subtraced and corrected for galactic and intrinsic extinction.


  \begin{figure}
\par\includegraphics[width=7.7cm,clip]{6660fig3.eps}\end{figure} Figure 3: Spectral Flux Distribution of the afterglow from radio to X-rays 0.6, 6.0 and 13.0 days after the burst in the observer frame. Several 3-$\sigma $ upper limits from radio and millimetre observations are plotted.
Open with DEXTER

4 Discussion

The optical-nIR SFD shows a clear curvature, implying a need for extinction along the afterglow line of sight in the host galaxy. The only reliable fit ( $\chi / {\rm d.o.f.} = 5.0/5$) is based on a SMC extinction law. This result points towards a low stage of chemical enrichment in the region of the progenitor, as is usually found for LGRBs (Kann et al. 2006). The detection of a starburst host galaxy is also a common feature with most LGRBs (Fruchter et al. 2006), facts that once again favour the hypothesis of shared nature between XRR bursts and LGRBs.

The optical spectral index obtained from this fit ( $\beta_{\rm o} =0.28$ $\pm$ 0.33) and the X-ray one ( $\beta_{\rm X} = 1.14$ $\pm$ 0.19, Nousek et al. 2006) are consistent with a synchrotron spectra in which the cooling break frequency ( $\nu_{\rm c}$) is located between optical and X-rays. A linear fit of the optical ($R{\rm c}$) and X-ray data between 0.1 and 1.0 days (where there is more data available and the light curve seems more stable) returns temporal slope values of $\alpha_{\rm o} = 0.69$ $\pm$ 0.04 and $\alpha_{\rm X} = 0.99$ $\pm$ 0.21. These numbers, together with the optical and X-ray spectral slopes, are consistent with a standard fireball model (Sari et al. 1999) in which a relativistic outflow is expanding in a uniform density environment in the slow cooling regime with an electron power law distribution index of $p \sim 2.0$.

A more complex multirange model, confirms these results and is used to account for the bump that has been detected to peak at $\sim$6 days by allowing for refreshed shocks. This fluctuation is simultaneously observed in optical and X-rays and can be explained by an energy injection of the order of the initial shock. This achromaticity and the simultaneity at both sides of  $\nu_{\rm c}$ rules out other explanations such as a density fluctuation, a dust echo or a supernova bump (which could also be ruled out by amplitude and onset time). Other explanations involving an refreshed energy release such as a double jet (Berger et al. 2003) or a patchy shell (Meszaros et al. 1998) can not be discarded. This injection delays the break, that would be expected for about 1.6 days, to about 30 days and, mainly due to the effect of the equal arrival time surface, transforms it to a very smooth break that expands over a decade in time.

To explain the flattening seen in the earliest points of the light curve we have studied the case of a collimated fireball seen off-axis, as predicted by some unified models (Yamazaki et al. 2002) that simultaneously intend to explain LGRBs, XRR bursts and XRFs by only varying the viewing angle. Our fit accounts reasonably well for the multiband and long scale behaviour of the light curve. However, the fits obtained with an on-axis model with an additional early injection or a low initial gamma factor (dirty fireball) can also interpret the data (although returning worse fits) and can not be ruled out.

Regarding the energetics of the afterglow we find that, with an observed peak energy of $\sim$20 keV and a fluence of $\sim$3.3 $\times $ 10-6 erg cm-2 (2-400 keV) it has an isotropic equivalent energy release in $\gamma$-rays $E_{\gamma,{\rm iso}} \ga
1.3$ $\times $ 1052 erg, at least 6 times greater than the predicted by $E_{\rm peak}-E_{\rm iso}$ relation (Amati 2006).

We encourage polarimetric observations of XRR bursts and XRF events (i.e. Gorosabel et al. 2006) as they will be extremely useful to better understand the physics and geometry of the emission and to discriminate between energy models when explaining the fluctuations seen in the light curves.

Acknowledgements
We acknowledge the generous allocation of observing time by different Time Allocation Committees. This work was partially supported by the Spanish MCyT under programmes AYA2004-01515 and ESP2005-07714-C03-03 (including FEDER funds), RFBR grants 04-02-16300 and 05-02-17744 and grant NSh-784.2006.2. The Dark Cosmology Centre is supported by the Danish National Research Foundation. A.d.U.P. acknowledges support from FPU grant AP2002-0446 from the Spanish MCyT. G.J. acknowledges support from the Icelandic Research Council.

References

 

  
Online Material

Appendix A: Calibration stars

Table A.1: Calibration stars in the field of GRB 050408, as marked in Fig. A.1.

Optical photometric calibration is based on the observation of several standard fields (Landolt 1992) using the 1.54 m Danish at La Silla (Chile) and the 1.5 m telescope at Sierra Nevada Observatory (Spain). From these observations we derive 12 secondary standards of different brightnesses in the field of the GRB. Identification of these stars is shown in Fig. A.1, while their photometric values are displayed in Table A.1.


  \begin{figure}
\par\includegraphics[width=7.4cm,clip]{6660fig4.eps} \end{figure} Figure A.1: Identification chart of GRB 050408. The calibration stars and the afterglow location have been marked. The field of view is 6.0' $\times $ 6.5', North is to the top and East to the left.

Appendix B: Observations

In the following tables we use the HETE-2 onset time t0 = 2005 April 08 16:22:50.93 UT.

Table B.1: Observations in millimetre wavelengths from the Plateau du Bure interferometer (compact 6-antenna configuration on all dates). The errors are based on point-source fits in the UV plane to the phase center.

Table B.2: Observations in radio wavelengths from the Giant Metrewave Radio Telescope (GMRT) and Radio Astronomical Telescope Academy Nauk (RATAN-600).

Table B.3: Optical observations carried out for GRB 050408. The magnitudes are in the Vega system and are not corrected for Galactic reddening.



Copyright ESO 2007