Issue |
A&A
Volume 495, Number 2, February IV 2009
|
|
---|---|---|
Page(s) | 505 - 512 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361:200811072 | |
Published online | 14 January 2009 |
HESS observations of
-ray bursts in 2003-2007
F. Aharonian1,13 - A. G. Akhperjanian2 - U. Barres de Almeida8,
- A. R. Bazer-Bachi3 - B. Behera14 - W. Benbow1 - K. Bernlöhr1,5 - C. Boisson6 - A. Bochow1 - V. Borrel3 - I. Braun1 - E. Brion7 - J. Brucker16 - P. Brun7 - R. Bühler1 - T. Bulik24 - I. Büsching9 - T. Boutelier17 - S. Carrigan1 - P. M. Chadwick8 - A. Charbonnier19 - R. C. G. Chaves1 - A. Cheesebrough8 - L.-M. Chounet10 - A. C. Clapson1 - G. Coignet11 - M. Dalton5 - B. Degrange10
- C. Deil1 - H. J. Dickinson8 - A. Djannati-Ataï12 - W. Domainko1 - L. O'C. Drury13 - F. Dubois11 - G. Dubus17 - J. Dyks24 - M. Dyrda28 - K. Egberts1 - D. Emmanoulopoulos14 - P. Espigat12 - C. Farnier15 - F. Feinstein15 - A. Fiasson15 - A. Förster1 - G. Fontaine10 - M. Füßling5 - S. Gabici13 - Y. A. Gallant15 - L. Gérard12 - B. Giebels10 - J. F. Glicenstein7 - B. Glück16 - P. Goret7 - C. Hadjichristidis8 - D. Hauser14 - M. Hauser14 - S. Heinz16 - G. Heinzelmann4 - G. Henri17 - G. Hermann1 - J. A. Hinton25 - A. Hoffmann18 - W. Hofmann1 - M. Holleran9 - S. Hoppe1 - D. Horns4 - A. Jacholkowska19 - O. C. de Jager9 - I. Jung16 - K. Katarzynski27 - S. Kaufmann14 - E. Kendziorra18 - M. Kerschhaggl5 - D. Khangulyan1 - B. Khélifi10 - D. Keogh8 - Nu. Komin7 - K. Kosack1 - G. Lamanna11 - J.-P. Lenain6 - T. Lohse5 - V. Marandon12 - J. M. Martin6 - O. Martineau-Huynh19 - A. Marcowith15 - D. Maurin19 - T. J. L. McComb8 - M. C. Medina6 - R. Moderski24 - E. Moulin7 - M. Naumann-Godo10 - M. de Naurois19 - D. Nedbal20 - D. Nekrassov1 - J. Niemiec28 - S. J. Nolan8 - S. Ohm1 - J.-F. Olive3 - E. de Oña Wilhelmi12,29 - K. J. Orford8 - J. L. Osborne8 - M. Ostrowski23 - M. Panter1 - G. Pedaletti14 - G. Pelletier17 - P.-O. Petrucci17 - S. Pita12 - G. Pühlhofer14 - M. Punch12 - A. Quirrenbach14 - B. C. Raubenheimer9 - M. Raue1,29 - S. M. Rayner8 - M. Renaud1 - F. Rieger1,29 - J. Ripken4 - L. Rob20 - S. Rosier-Lees11 - G. Rowell26 - B. Rudak24 - C. B. Rulten8 - J. Ruppel21 - V. Sahakian2 - A. Santangelo18 - R. Schlickeiser21 - F. M. Schöck16 - R. Schröder21 - U. Schwanke5 - S. Schwarzburg18 - S. Schwemmer14 - A. Shalchi21 - J. L. Skilton25 - H. Sol6 - D. Spangler8 -
.
Stawarz23 - R. Steenkamp22 - C. Stegmann16 - G. Superina10 - P. H. Tam14 - J.-P. Tavernet19 - R. Terrier12 - O. Tibolla14 - C. van Eldik1 - G. Vasileiadis15 - C. Venter9 - L. Venter6 - J. P. Vialle11 - P. Vincent19 - M. Vivier7 - H. J. Völk1 - F. Volpe10,29 - S. J. Wagner14 - M. Ward8 - A. A. Zdziarski24 - A. Zech6
1 -
Max-Planck-Institut für Kernphysik, PO Box 103980, 69029
Heidelberg, Germany
2 -
Yerevan Physics Institute, 2 Alikhanian Brothers St., 375036 Yerevan,
Armenia
3 -
Centre d'Étude Spatiale des Rayonnements, CNRS/UPS, 9 Av. du Colonel Roche, BP
4346, 31029 Toulouse Cedex 4, France
4 -
Universität Hamburg, Institut für Experimentalphysik, Luruper Chaussee
149, 22761 Hamburg, Germany
5 -
Institut für Physik, Humboldt-Universität zu Berlin, Newtonstr. 15,
12489 Berlin, Germany
6 -
LUTH, Observatoire de Paris, CNRS, Université Paris Diderot, 5 place Jules Janssen, 92190 Meudon,
France
7 -
IRFU/DSM/CEA, CE Saclay, 91191
Gif-sur-Yvette Cedex, France
8 -
University of Durham, Department of Physics, South Road, Durham DH1 3LE,
UK
9 -
Unit for Space Physics, North-West University, Potchefstroom 2520,
South Africa
10 -
Laboratoire Leprince-Ringuet, École Polytechnique, CNRS/ IN2P3,
91128 Palaiseau, France
11 -
Laboratoire d'Annecy-le-Vieux de Physique des Particules,CNRS/IN2P3,
9 chemin de Bellevue, BP 110, 74941 Annecy-le-Vieux Cedex, France
12 -
Astroparticule et Cosmologie (APC), CNRS, Université Paris 7 Denis Diderot,
10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 13, France (UMR 7164, CNRS, Université Paris VII, CEA, Observatoire de Paris)
13 -
Dublin Institute for Advanced Studies, 5 Merrion Square, Dublin 2,
Ireland
14 -
Landessternwarte, Universität Heidelberg, Königstuhl, 69117 Heidelberg, Germany
15 -
Laboratoire de Physique Théorique et Astroparticules, CNRS/IN2P3,
Université Montpellier II, CC 70, Place Eugène Bataillon, 34095
Montpellier Cedex 5, France
16 -
Universität Erlangen-Nürnberg, Physikalisches Institut, Erwin-Rommel-Str. 1,
91058 Erlangen, Germany
17 -
Laboratoire d'Astrophysique de Grenoble, INSU/CNRS, Université Joseph Fourier, BP
53, 38041 Grenoble Cedex 9, France
18 -
Institut für Astronomie und Astrophysik, Universität Tübingen,
Sand 1, 72076 Tübingen, Germany
19 -
LPNHE, Université Pierre et Marie Curie Paris 6, Université Denis Diderot
Paris 7, CNRS/IN2P3, 4 place Jussieu, 75252 Paris Cedex 5, France
20 -
Institute of Particle and Nuclear Physics, Charles University,
V Holesovickach 2, 180 00 Prague 8, Czech Republic
21 -
Institut für Theoretische Physik, Lehrstuhl IV: Weltraum und
Astrophysik,
Ruhr-Universität Bochum, 44780 Bochum, Germany
22 -
University of Namibia, Private Bag 13301, Windhoek, Namibia
23 -
Obserwatorium Astronomiczne, Uniwersytet Jagiellonski, ul. Orla 171,
30-244 Kraków, Poland
24 -
Nicolaus Copernicus Astronomical Center, ul. Bartycka 18, 00-716 Warsaw,
Poland
25 -
School of Physics & Astronomy, University of Leeds, Leeds LS2 9JT, UK
26 -
School of Chemistry & Physics,
University of Adelaide, Adelaide 5005, Australia
27 -
Torun Centre for Astronomy, Nicolaus Copernicus University, ul.
Gagarina 11, 87-100 Torun, Poland
28 -
Instytut Fizyki Jadrowej PAN, ul. Radzikowskiego 152, 31-342 Kraków,
Poland
29 -
European Associated Laboratory for Gamma-Ray Astronomy, jointly
supported by CNRS and MPG
Received 2 October 2008 / Accepted 22 December 2008
Abstract
Aims. Very-high-energy (VHE; 100 GeV)
-rays are expected from
-ray bursts (GRBs) in some scenarios. Exploring this photon energy regime is necessary for understanding the energetics and properties of GRBs.
Methods. GRBs have been one of the prime targets for the HESS experiment, which makes use of four Imaging Atmospheric Cherenkov Telescopes (IACTs) to detect VHE -rays. Dedicated observations of 32 GRB positions were made in the years 2003-2007 and a search for VHE
-ray counterparts of these GRBs was made. Depending on the visibility and observing conditions, the observations mostly start minutes to hours after the burst and typically last two hours.
Results. Results from observations of 22 GRB positions are presented. Evidence of a VHE signal was found neither in observations of any individual GRBs nor from stacking data from subsets of GRBs with higher expected VHE flux according to a model-independent ranking scheme. Upper limits for the VHE -ray flux from the GRB positions were derived. For those GRBs with measured redshifts, differential upper limits at the energy threshold after correcting for absorption due to extra-galactic background light are also presented.
Key words: gamma rays: bursts - gamma rays: observations
1 Introduction
Gamma-ray bursts (GRBs) are the most energetic events in
the -ray regime. Depending on their duration (e.g. T90), GRBs are categorized into long GRBs (T90>2 s) and short GRBs (T90<2 s). First detected in late 1960s (Klebesadel et al. 1973), GRBs remained mysterious for three decades.
Breakthroughs in understanding GRBs came only after
the discovery of longer-wavelength afterglows
with the launch of BeppoSAX in 1997 (van Paradijs et al. 2000). Multi-wavelength (MWL) observations have proved to be crucial in our understanding of GRBs, and provide valuable information about their physical
properties. These MWL afterglow observations are generally explained by synchrotron emission from shocked electrons in the relativistic fireball model (Zhang & Mészáros 2004; Piran 1999). A plateau phase is revealed in many of the Swift/XRT light curves, the origin of which is still not clear (Zhang et al. 2006). Observations of GRBs at energies >10 GeV may test some of the ideas that have been suggested to explain the X-ray observations (Fan et al. 2008).
In the framework of the relativistic fireball model, photons with energies up to 10 TeV or higher are expected from the GRB afterglow phase (Zhang & Mészáros 2004; Fan & Piran 2008). Possible leptonic radiation mechanisms include forward-shocked electrons up-scattering self-emitted synchrotron
photons (SSC processes; Fan et al. 2008; Dermer et al. 2000; Zhang & Mészáros 2001) or photons from other shocked
regions (Wang et al. 2001). Physical parameters, such as the ambient density of the surrounding
material (n), magnetic field equipartition fraction (
), and bulk
Lorentz factor (
)
of the outflow, may be constrained by observations at these energies
(Pe'er & Waxman 2005; Wang et al. 2001).
A possible additional contribution to VHE emission
relates to the X-ray flare phenomenon. X-ray flares are
found in more than 50% of the Swift GRBs during the afterglow phase (Chincarini et al. 2007). The
energy fluence of some of them (e.g. GRB 050502B) is comparable to that of the
prompt emission. Most of them
are clustered at 102-103 s after the GRB (see Fig. 2
in Chincarini et al. 2007), while late X-ray flares (>104 s) are also
observed; when these happen they can cause an increase in the X-ray flux of
an order of magnitude or more over the power-law temporal
decay (Curran et al. 2008). The cause of X-ray flares is still a subject of
debate, but corresponding VHE
-ray flares from inverse-Compton (IC) processes
are predicted (Wang et al. 2006; Galli & Piro 2007; Fan et al. 2008). The accompanying external-Compton flare may be
weak if the flare originated behind the
external shock, e.g. from prolonged central engine activity (Fan et al. 2008). However, in the
external shock model, the expected SSC flare is very strong at GeV energies
and can be readily detected using a VHE instrument
with an energy threshold of
100 GeV (Galli & Piro 2008), such as
the HESS array, for a typical GRB at
.
Therefore, VHE
-ray data taken during an X-ray flare may help for distinguishing the internal/external shock origin of the X-ray flares, and may be used as a diagnostic tool for the late central engine
activity.
Waxman & Bahcall (2000) and Murase et al. (2008) suggest that GRBs may be sources of ultra-high-energy cosmic rays (UHECRs). In this case, -decays from proton-
interaction may generate VHE emission. The VHE
-ray emission produced from such a hadronic component is generally expected to
decay more slowly than the leptonic sub-MeV radiation (Böttcher & Dermer 1998). Dermer (2007)
suggests a combined leptonic/hadronic scenario to explain the
rapidly-decaying phase and plateau phase seen in many of the Swift/XRT
light curves. This model can be tested with VHE observations taken minutes to hours after the burst.
Most searches for VHE -rays from GRBs have obtained negative results
(Atkins et al. 2005; Connaughton et al. 1997). There may be indications of excess photon
events from some observations, but these results are not conclusive
(Amenomori et al. 1996; Poirier et al. 2003; Atkins et al. 2000; Padilla et al. 1998). Currently, the most
sensitive detectors in the VHE
-ray regime are IACTs. Horan et al. (2007) presented upper limits from 7 GRBs observed with the Whipple Telescope during the pre-Swift era. Upper limits for 9 GRBs with redshifts that were either unknown or >3.5 were also reported by the MAGIC collaboration
(Albert et al. 2007). In general,
these limits do not violate a power-law extrapolation of the keV spectra
obtained with satellite-based instruments. However, most GRBs are now believed to originate at cosmological distances, therefore absorption of VHE
-rays by the EBL (Nikishov 1962) must be considered when interpreting these limits.
In this paper, observations of 22 -ray bursts made with
HESS during the years 2003-2007 are reported. They represent the largest sample of GRB afterglow observations made by an IACT array and result in the most stringent upper limits obtained in the VHE band. The prompt phase of GRB 060602B was observed serendipitously with HESS The results of
observations before, during, and after this burst are presented in Aharonian et al. (2009).
2 The HESS experiment and GRB observation strategy
The HESS array is a
system of four 13 m-diameter IACTs located at 1800 m above sea level in the
Khomas Highland of Namibia (
S,
E). Each of the four telescopes is located at a
corner of a square with a side length of 120 m. This configuration was
optimized for maximum sensitivity to
100 GeV photons. The effective
collection area increases from
10
at 100 GeV to more than
at 1 TeV for observations at a zenith angle (ZA) of 20
.
The
system has a point source sensitivity above 100 GeV of
(
of the flux from the
Crab nebula) for a
detection in a 2 h observation. Each
HESS camera consists of 960 photomultiplier tubes (PMTs), which in total
provide a field of view (FoV) of
5
.
This relatively large
FoV allows for the simultaneous determination of the background events from
off-source positions, so that no dedicated off run is
needed (Aharonian et al. 2006c). The slew rate of the array is
100
per
minute, enabling it to point to any sky position within
2 min. The
HESS array is currently the only IACT array in the Southern Hemisphere
used for an active GRB observing
programme
.
The trigger system of the HESS array is described in Funk et al. (2004). The stereoscopic technique is used, i.e. a coincidence of at least two telescopes triggering within a window of (normally) 80 nanoseconds is required. This largely rejects background events caused by local muons that trigger only a single telescope.
The observations reported here were obtained over the period March 2003 to October 2007. The observations of two GRBs in 2003 were made using two telescopes while the system was under construction. Before 2003 July, each of the two telescopes took data separately. Stereo analysis was then performed on the data, which requires coincidence of events to be determined offline using GPS time stamps. After the installation of the central trigger system in 2003 July, the stereo multiplicity requirement was determined on-line. All observations since 2004 have made use of the completed four-telescope array and the stereo technique (Aharonian et al. 2006c).
Most of the data were taken in 28 min runs using wobble mode, i.e. the
GRB position is placed at an offset,
,
of
or
(in RA and Dec) relative to the centre of the camera FoV during observations.
Onboard GRB triggers distributed by the Swift
satellite, as well as triggers from INTEGRAL and HETE-II confirmed by
ground-based analysis, are followed by HESS observations. Upon the reception of a GCN notice from one of
these satellites (with appropriate indications
that the source is a genuine GRB),
the burst position is observed if ZA
45
(to ensure a
reasonably low-energy threshold) during HESS dark
time
. An
automated program is running on site to keep the shift crew alerted of any new
detected GRBs in real time. Depending on the observational constraints and the
measured redshifts of the GRBs reported through GCN circulars
, observations of the burst positions are started up to
24 h after the burst time,
typically with an exposure time of
120 min in wobble
mode. The remarkably nearby, bright GRB 030329 was an exceptional case. It was not observed
until 11.5 days after the burst because of poor weather, which
prohibited observation any earlier.
3 The GRB observations
Thirty-two GRBs were observed with HESS during the period from March 2003 to October 2007. After applying a set of data-quality criteria that rejects observation runs with non-optimal weather conditions and hardware status, 22 GRB observations were selected for analysis and are described in this section.
Table 1: Properties of GRBs observed with HESS from March 2003 to October 2007.
Table 2: HESS observations of GRBs from March 2003 to October 2007.
3.1 Properties of the GRBs
For each burst, the observational properties as obtained from the triggering satellite are
shown in Table 1. These include trigger number, energy band,
fluence in that energy band, and the duration of the burst
(T90). Whenever there were follow-up observations in the X-ray, optical, or
radio bands, whether a detection has occurred (denoted by a tick )
or not (denoted by a cross
)
is also shown. If no observation at a given wavelength was
reported, a dot (.) is shown. The reported redshifts (z) of 10 GRBs
are also presented, of which 6 are lower than one. Two observed bursts,
GRB 070209 and GRB 070724A, are short
GRBs while the rest are long GRBs. The population of short GRBs has a
redshift distribution (Berger et al. 2007) significantly less than that of the long
GRBs (Jakobsson et al. 2006). Therefore, on average they are likely to suffer from a lower level of EBL absorption.
X-ray flares were detected from three of the GRBs in the HESS sample. They occurred at 273 s
after the burst for GRB 050726, 284 s for GRB 050801, and
s for GRB 070429A (Curran et al. 2008). Unfortunately, the flares occurred
outside the time windows of the HESS observations.
3.2 HESS observations
For each burst, the start time,
,
of the HESS observations
after the burst is shown in Table 2. Since an observing strategy to start observing the burst position up to
24 h after the burst time is applied,
the mean
is of the
order of 10 h. The (good-quality) exposure time of the observations using
telescopes for each burst is included. The mean ZA of the
observations is also presented.
3.3 The ranking scheme
As mentioned in the introduction, there is no lack of models predicting VHE emission from GRBs. However, the evolution of the possible VHE -ray emission with time is model-dependent. To give an empirical, model-independent estimate of the relative expected
VHE flux of each GRB (which also depends on
), it is assumed that: (1) the relative VHE signal scales
as the energy released in the prompt emission, taken as a typical energy
measure of a GRB. Hence
where
is the fluence in the Swift/BAT band. For bursts not triggered by BAT, the measured fluence is extrapolated into this energy band; (2) the possible VHE signal fades as
time goes on, as observed in longer wavelength (e.g. X-ray) data. In particular, the VHE
flux follows the average decay of the X-ray flux and therefore
where t denotes the time after the burst and 1.3 is the average X-ray afterglow late-time power-law decay index (Nousek et al. 2006). Since in most cases the exposure time of
the observations is much shorter than
(the start time of the
corresponding HESS observations after the trigger), the expected
flux at
can be used as
a measure of the strength of the VHE signal, and therefore of the relative
possibility of detecting a VHE signal from that GRB. By setting t to
,
we have
The rank of each GRB according to Eq. (1) is shown in the last column in Table 1. Note that redshift information (available for only a few GRBs), and thus the corresponding EBL absorption, is not taken into account in the ranking scheme.
4 Data analysis
Calibration of data, event reconstruction and rejection of the cosmic-ray background (i.e. -ray event selection criteria) were performed as described in Aharonian et al. (2006c), which employs the techniques described by Hillas (1996).
Gamma-like events were then taken from a circular region (on-source) of radius
centred at the burst position given in
Table 1. The background was
estimated using the reflected-region background model as described
in Berge et al. (2007), in which the number of background events in the on-source
region (
)
is
estimated from
off-source regions located at the same
as the on-source region during the same observation. The number of
-like events is given by
where
is the total number of events detected in the on-source
region and
the normalization factor.
Independent analyses of various GRBs using different methods and background estimates (Berge et al. 2007) yielded consistent results.
4.1 Analysis technique
Two sets of analysis cuts were applied
to search for a VHE -ray signal from observational data taken with three or four telescopes. These are ``standard'' cuts (Aharonian et al. 2006c) and ``soft'' cuts
(the latter have lower energy thresholds, as described in Aharonian et al. 2006a).
For standard (soft) cuts,
(
). While standard cuts are optimized for a source with a power-law spectrum of photon index
,
soft cuts are optimized for a source with a steep spectrum (
), and have better sensitivity at lower energies. Since EBL absorption is less severe for lower energy photons, the soft-cut analysis is useful in searching for VHE
-rays from GRBs which are at cosmological distances. For example, the photon indices of two blazars PKS 2005-489 (Aharonian et al. 2005) and PG 1553+113 (Aharonian et al. 2008) were measured to be
.
An exception to this analysis scheme is GRB 030329. As the central trigger system had yet to be installed when this observation was made, a slightly different analysis technique was used. The description of the image and analysis cuts used for the data from GRB 030329 can be found in Aharonian et al. (2005). For GRB 030821, only the standard-cut analysis (for two-telescope data) was performed (see Sect. 5.4).
The positional error circle of most GRBs, with the exceptions of GRB 030821, GRB 050209, and
GRB 070209, is small compared to the HESS point spread function (PSF). The 68%
-ray containment radius,
,
of the HESS PSF can be as small
as
3
,
depending on the ZA and
of the observations, and the analysis cuts applied. The 68% containment radius,
,
of the observations of GRB 050209 and
GRB 070209 is about 9
using standard-cut analysis
, slightly larger
than the corresponding error circles. Therefore, point-source analyses were
performed for all GRBs except GRB 030821, the error box of which is much bigger than the HESS PSF (see Sect. 5.4 for its treatment).
4.2 Energy threshold
The energy threshold,
,
is conventionally defined as the peak in
the differential
-ray rate versus energy
curve of a fictitious source with photon index
(Konopelko et al. 1999). This curve is a
convolution of the effective area with the expected energy spectrum of the
source as seen on Earth. Such energy thresholds, obtained by the
standard-cut analysis and the soft-cut analysis for each GRB observation, are
shown in Table 2, assuming
.
The energy threshold depends on the ZA of the observations and the analysis used.
The larger the ZA, the higher is the energy threshold. Moreover, soft-cut
analysis gives a lower value of
than that of standard-cut
analysis. Note that
-ray photons with energies below
can
be detected by the telescopes.
4.3 Optical efficiency of the instrument
The data presented were also corrected for the long-term changes in the optical efficiency of the instrument. The optical efficiency has decreased over a period of a few years. This has changed the effective area and energy threshold of the instrument. Specifically, the energy threshold has increased with time. Using images of local muons in the FoV, this effect in the calculation of flux upper limits is corrected (cf. Aharonian et al. 2006c).
5 Results
No evidence of a significant excess of VHE -ray events from any of the
GRB positions given in Table 1 during the period covered by the
HESS observations was found. The number of on-source (
)
and off-source events (
), normalization factor (
),
excess, and statistical significance
of the excess in standard deviations (
)
are given for each of the 21 GRBs in
Table 2. The results for GRB 030821 are given in
Sect. 5.4. Figure 1 shows the distribution of
the significance obtained from the soft-cut analysis of the
observations of each of the 21 GRBs. A Gaussian distribution with mean zero
and standard deviation one, which is expected in the case of no
detection, is shown for comparison. The distribution of the statistical
significance is consistent with this Gaussian distribution. Thus no
significant signal was found from any of the individual GRBs. A search for serendipitous source discoveries in
the HESS FoV during observations of the GRBs also resulted in no significant
detection. The 99.9% confidence level (c.l.) flux upper limits (above
)
have been calculated using the method
of Feldman & Cousins (1998) for both standard cuts (assuming
)
and soft cuts (assuming
), and are included in
Table 2. The limits are as observed on Earth, i.e. the
EBL absorption factor was not taken into account. The systematic error on a
HESS integral flux measurement is estimated to be
20%, and it was not
included in the calculation of the upper limits.
For those GRBs with reported
redshifts, the effect of the EBL on the HESS limits can be
estimated. Using the EBL model P0.45 described in Aharonian et al. (2006b), differential upper
limits (again assuming )
at the energy threshold were calculated from the integral
upper limits obtained using soft-cut analysis. These upper limits, as
well as those calculated without taking the EBL into account, are shown in Table 3.
![]() |
Figure 1: Distribution of the statistical significance ( histogram) as derived from the observations of 21 GRBs using soft-cut analysis. The mean is -0.4 and the standard deviation is 1.4. Each entry corresponds to one GRB. The solid line is a Gaussian function with mean zero and standard deviation unity. |
Open with DEXTER |
Table 3: Differential flux upper limits at the energy thresholds from the HESS observations of GRBs with reported redshifts.
Table 4: Combined significance of 3 subsets of GRBs selected based on the requirements listed in Sect. 5.1.
5.1 Stacking analysis
Although no significant excess was found from any individual GRB, co-adding the excess events from the observations of a number of GRBs may reveal a signal that is too weak to be seen in the data from one GRB, provided that the PSFs of the HESS observations are bigger than the error box of the GRB positions (which is the case, see Sect. 4.1). Firstly, stacking of all GRBs (except GRB 030821, which has a high positional uncertainty) in the sample was performed. This yielded a total of -157 excess events and a statistical significance of -1.98 using the soft-cut analysis. Use of standard cuts produced a similar result (see Table 4). Secondly, combining the significance of the results from three selected subsets extracted from the whole sample was performed. The a priori selection criteria were to choose those GRBs with a higher expected VHE flux or a lower level of EBL absorption. The following requirements were used to select three subsets:
- Sample A: the first 10 in the ranking described in Sect. 3.3;
- Sample B: all GRBs with a measured redshift z<1;
- Sample C: all GRBs with a soft-cut energy threshold lower than 300 GeV and with either a measured redshift z<1 or with an unknown redshift.
5.2 Temporal analysis
As possible VHE radiation from GRBs is expected to vary with time, a temporal analysis to search for deviation from zero excess in the observed data was performed. Soft-cut analysis was used for all GRBs (except GRB 030329)
since this analysis has a lower energy
threshold and a better acceptance of -rays and cosmic rays and therefore
increases the statistics. The
-like excess events were binned in 10-min time intervals
for each GRB data set and were compared to the assumption of no excess
throughout the observed period.
The
/d.o.f. value and the
corresponding probability are shown in Table 2 for each
GRB. Within the whole sample, the lowest probability that the
hypothesis that the excess was zero throughout the observation period is
correct is
(for GRB 071003) and no significant
deviation from zero within any of the GRB temporal
data was found. Standard-cut analysis produced consistent results.
5.3 GRB 070621: observations of a GRB with the fastest reaction and the longest exposure time
GRB 070621 is the highest-ranked GRB in the sample
(Sect. 3.3), i.e. it has the highest relative expected VHE flux at
the start time of the observations. The duration of the Swift burst was
s, thus clearly classifying the burst as a long GRB. The fluence in
the 15-150 keV band was
erg cm-2. The XRT
light curve is represented by an initial rapidly-decaying phase and a shallow phase, with the transition happening around t0 + 380s where t0 denotes the trigger
time (Sbarufatti et al. 2007). Despite extensive optical monitoring, no fading
optical counterpart was found. The HESS observations started at t0+420 s
and lasted for
5 h, largely coincident with the X-ray shallow
phase. These observations were both the most prompt and the longest among those presented. Figure 2 shows the 99.9% HESS energy flux upper limits
above 200 GeV (using soft-cut analysis), together with
the XRT results (Evans et al. 2007). As seen, the limits for
this period are at levels comparable to the X-ray energy flux during the same
period. Unfortunately the lack of redshift information for
this burst prevents further interpretation of the limits.
![]() |
Figure 2: The 99.9% confidence level energy flux upper limits (in red) at energies >200 GeV derived from HESS observations at the position of GRB 070621. The ends of the horizontal lines indicate the start and end times of the observations from which the upper limits were derived. The XRT energy flux in the 0.3-10 keV band is shown in black for comparison (Evans et al. 2007). |
Open with DEXTER |
5.4 GRB 030821: observations of a GRB with a high positional uncertainty
Some GRBs, such as GRB 030821, have a high uncertainty in position; with a relatively large camera
FoV (5
), the HESS telescopes are able to cover the whole
positional error box of such GRBs.
Observations of GRB 030821 started 18 h after the burst
and lasted for a live-time of 55.5 min, with a mean ZA of 28
.
The observations were taken when the array was under construction and only two telescopes were operating, resulting in an energy threshold of 260 GeV. The GRB has a relatively high uncertainty in
position as determined from IPN (the third Interplanetary Network)
triangulation (Hurley et al. 2003), and its error box is bigger than the
PSF of HESS However, because of the relatively large FoV of the camera, the whole
error box, and thus the possible GRB position, is within the HESS FoV. The sky excess map overlaid with the error box is shown in
Fig. 3. As can be seen, there is no significant excess at any position within the error
box. The sky region with the largest number of peak excess events is located in the south-eastern part of the error box. Using a point-source analysis centred at this peak, a flux upper limit (above 260 GeV) of
cm-2 s-1 was derived. Since an upper limit derived for any location in the error box with fewer excess events is lower than this value
, it may be regarded as a conservative upper limit of the VHE flux associated with GRB 030821 during the period of the HESS observations.
![]() |
Figure 3:
The |
Open with DEXTER |
6 Discussion
The upper limits presented in this paper are among the most stringent
ever derived from VHE -ray observations of GRBs during the
afterglow period. In fact, the 99.9% confidence level limits (in energy flux) are at levels comparable to the X-ray energy flux as observed by Swift/XRT
during the same period (see, e.g. Fig. 2).
Unless most of the GRBs are located at high redshifts and thus their VHE flux is severely
absorbed by the EBL (this possibility is discussed below), one expects detection of the predicted VHE component with
energy flux levels comparable to those in X-rays in some
scenarios (Pe'er & Waxman 2005; Fan et al. 2008; Dermer et al. 2000; Wang et al. 2001; Zhang & Mészáros 2001).
On the other hand, the unknown redshifts of many of the GRBs in the sample
(including GRB 070621, the highest-ranking, which is discussed in Sect. 3.3) complicate the physical interpretation of the data,
because EBL absorption at VHE energies is severe for a GRB with z>1. The
mean and median redshift of the 10 GRBs with reported redshifts is 1.3
and 0.7, respectively. If the 12 GRBs without redshift have the same redshift distribution, one would expect 40% of them (
5 GRBs) to have z<0.5. In this case, the EBL absorption may not preclude the detection of the predicted VHE
-rays for the GRB sample
presented here
.
There is no reported X-ray flare during the HESS observational time windows, therefore no conclusion on whether or not X-ray flares are accompanied by VHE flares, as well as the origin of X-ray flares, can be drawn. If UHECRs are generated in nearby GRB sources, as suggested by some authors, a detectable VHE flux is expected from nearby GRBs. Therefore, although the unknown redshifts of a significant fraction of GRBs in our sample (12 out of 22) and the uncertainty in the modeled VHE temporal evolution are surely in play, the results presented here do not indicate (but also not exclude) that GRBs are dominant sources of UHECRs.
7 Outlook
The data from our sample of 22 GRBs do not provide any evidence for a strong VHE -ray
component from GRBs during the afterglow phase. EBL absorption can explain the lack of detection in our sample.
However, this does not exclude a
population of GRBs that exhibit a strong VHE component. While the EGRET
experiment did not detect MeV-GeV photons from most BATSE GRBs in its FoV, some strong bursts
(e.g. GRB 940217) have proved to emit delayed emission,
1.5
hours after the burst, at energies as high as
20 GeV (Hurley et al. 1994). With Fermi's observations of GRBs having started in mid-2008, it is likely that our knowledge of the high-energy emission of GRBs will be improved in the near future.
The future prospects for detection at VHE energies rely on the likelihood of observing a GRB with low redshift (e.g. z<0.5) early enough. In the cases where there is no detection, sensitive and early upper limits on the intrinsic VHE luminosity of these nearby GRBs will still improve our understanding of the radiation mechanisms of GRBs. The existence of a distinct population of low-luminosity (LL) GRBs was suggested based on the high detection rate of low-redshift LL GRBs such as GRB 980425 and GRB 060218 (e.g. Soderberg et al. 2004). Due to their proximity, they are good targets for VHE observations. Since they are sub-energetic compared to other GRBs, they may be accompanied by a lower VHE luminosity. On the other hand, if most radiation are emitted at high energies, the detection probability would be much higher.
Franceschini et al. (2008) claimed a very small -ray opacity due to EBL absorption. The optical depth is about a factor of three less than the one we used, depending on the energies (Aharonian et al. 2006b). Therefore, on-going GRB observations with HESS, as well as other ground-based VHE detectors, are crucial to test this model.
8 Conclusions
During 5 years of operation (2003-2007), 32 GRBs were observed during the afterglow phase using the HESS experiment. Those 22 GRBs with high-quality data were analysed and the results presented in this paper. Depending on the visibility and observing conditions, the start time of the observations varied from minutes to hours after the burst.
There is no evidence of VHE emission from any individual
GRB during the period covered by the HESS observations, nor from stacking
analysis using the whole sample and a priori selected sub-sets of GRBs. Fine-binned temporal data revealed no short-term variability from any observation and no indication of VHE signal from any of these time bins was found. Upper limits of VHE -ray flux during the observations from the GRBs were
derived. These 99.9% confidence level energy flux upper limits are
at levels comparable to the contemporary X-ray energy flux. For those GRBs with reported redshifts,
differential upper limits at the energy threshold after correcting for EBL absorption
are presented.
HESS phase II will have an energy threshold of about 30 GeV. With much less
absorption by the EBL at such low energies, it is hoped that the HESS experiment will
enable the detection of VHE -ray counterparts of
GRBs.
Acknowledgements
The support of the Namibian authorities and of the University of Namibia in facilitating the construction and operation of HESS is gratefully acknowledged, as is the support by the German Ministry for Education and Research (BMBF), the Max Planck Society, the French Ministry for Research, the CNRS-IN2P3 and the Astroparticle Interdisciplinary Programme of the CNRS, the U.K. Science and Technology Facilities Council (STFC), the IPNP of the Charles University, the Polish Ministry of Science and Higher Education, the South African Department of Science and Technology and National Research Foundation, and by the University of Namibia. We appreciate the excellent work of the technical support staff in Berlin, Durham, Hamburg, Heidelberg, Palaiseau, Paris, Saclay, and in Namibia in the construction and operation of the equipment. P.H. Tam acknowledges support from IMPRS-HD. This work has made use of the GCN Notices and Circulars provided by NASA's Goddard Space Flight Center, as well as data supplied by the UK Swift Science Data Centre at the University of Leicester.
References
- Aharonian, F. A., Akhperjanian, A. G., Aye, K.-M., et al. (HESS collaboration) 2005, A&A, 436, L17 (In the text)
- Aharonian, F. A., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. (HESS collaboration) 2006a, A&A, 448, L19 (In the text)
- Aharonian, F. A., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. (HESS collaboration) 2006b, Nature, 440, 1018 (In the text)
- Aharonian, F. A., Akhperjanian, A. G., Bazer-Bachi, A. R., et al. (HESS collaboration) 2006c, A&A, 457, 899 (In the text)
- Aharonian, F. A., Akhperjanian, A. G., Barres de Almeida, U., et al. (HESS collaboration) 2008, A&A, 477, 481 (In the text)
- Aharonian, F. A., Akhperjanian, A. G., Barres de Almeida, U., et al. (HESS collaboration) 2009, ApJ, 690, 1068 (In the text)
- Albert, J., Aliu, E., Anderhub, H., et al. (MAGIC collaboration) 2007, ApJ, 667, 358 (In the text)
- Atkins, R., Benbow, W., Berley, D., et al. 2000, ApJ, 533, L119 [NASA ADS] [CrossRef]
- Atkins, R., Benbow, W., Berley, D., et al. 2005, ApJ, 630, 996 [NASA ADS] [CrossRef]
- Amenomori, M., Cao, Z., Dai, B. Z., et al. 1996, A&A, 311, 919 [NASA ADS]
- Baring, M. G. 2006, ApJ, 650, 1004 [NASA ADS] [CrossRef]
- Berge, D., Funk, S., & Hinton, J. A. 2007, A&A, 466, 1219 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Berger, E., & Fox, D. 2007, GCN Circ., 6101 [NASA ADS]
- Berger, E., & Gladders, M. 2006, GCN Circ., 5170 [NASA ADS]
- Berger, E., Fox, D. B., Price, P. A., et al. 2007, ApJ, 664, 1000 [NASA ADS] [CrossRef] (In the text)
- Böttcher, M., & Dermer, C. D. 1998, ApJ, 499, L131 [NASA ADS] [CrossRef] (In the text)
- Butler, N. R., Kocevski, D., Bloom, J. S., et al. 2007, ApJ, 671, 656 [NASA ADS] [CrossRef]
- Chincarini, G., Moretti, A., Romano, P., et al. 2007, ApJ, 671, 1903 [NASA ADS] [CrossRef] (In the text)
- Connaughton, V., Akerlof, C. W., Barthelmy, S., et al. 1997, ApJ, 479, 859 [NASA ADS] [CrossRef]
- Cucchiara, A., Fox, D. B., & Cenko, S. B., et al. 2007, GCN Circ., 6665 [NASA ADS]
- Curran, P. A., Starling, R. L. C., O'Brien, P. T., et al. 2008, A&A, 487, 533 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Dermer, C. D. 2007, ApJ, 664, 384 [NASA ADS] [CrossRef] (In the text)
- Dermer, C. D., Chiang, J., & Mitman, K. E. 2000, ApJ, 537, 785 [NASA ADS] [CrossRef]
- Evans, P. A., Beardmore, A. P., & Page, K. L., et al. 2007, A&A, 469, 379 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Fan, Y.-Z., & Piran, T. 2008, Frontiers of Physics in China, 3, 306 [NASA ADS] [CrossRef]
- Fan, Y.-Z., Piran, T., Narayan, R., & Wei, D.-M. 2008, MNRAS, 384, 1483 [NASA ADS] [CrossRef] (In the text)
- Feldman, G. J., & Cousins, R. D. 1998, Phys. Rev. D., 57, 3873 [CrossRef] (In the text)
- Franceschini, A., Rodighiero, G., & Vaccari, M. 2008, A&A, 487, 837 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Funk, S., Hermann, G., Hinton, J. A. et al. 2004, Astropart. Phys., 22, 285 [NASA ADS] [CrossRef] (In the text)
- Fynbo, J. P. U., Thoene, C. C., & Malesani, D., et al. 2007, GCN Circ., 6759 [NASA ADS]
- Galli, A., & Piro, L. 2007, A&A, 475, 421 [NASA ADS] [CrossRef] [EDP Sciences]
- Galli, A., & Piro, L. 2008, A&A, 489, 1073 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Hillas, A. M. 1996, Space Sci. Rev., 75, 17 (In the text)
- Horan, D., Atkins, R. W., Badran, H. M., et al. 2007, ApJ, 655, 396 [NASA ADS] [CrossRef] (In the text)
- Hurley, K., Dingus, B. L., Mukherjee, R., et al. 1994, Nature, 372, 652 [NASA ADS] [CrossRef] (In the text)
- Hurley, K., Cline, T., & Mazets, E., et al. 2003, GCN Circ., 2359 [NASA ADS] (In the text)
- Jakobsson, P., Levan, A., Fynbo, J. P. U., et al. 2006, A&A, 447, 897 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Klebesadel, R. W., Strong, I. B., & Olson, R. A. 1973, ApJ, 182, L85 [NASA ADS] [CrossRef] (In the text)
- Konopelko, A., Hemberger, M., Aharonian, F., et al. (HEGRA Collaboration) 1999, Astropart. Phys., 10, 275 (In the text)
- Li, T.-P., & Ma, Y.-Q. 1983, ApJ, 272, 317 [NASA ADS] [CrossRef]
- Malesani, D., Jakobsson, P., Fynbo, J. P. U., Hjorth, J., & Vreeswijk, P. M. 2007, GCN Circ., 6651 [NASA ADS]
- Murase, K., Ioka, K., Nagataki, S., & Nakamura, T. 2008, Phys. Rev. D, 78, 023005 [NASA ADS] [CrossRef] (In the text)
- Nikishov, A. I. 1962, Sov. Phys. JETP, 14, 393 (In the text)
- Nousek, J. A., Kouveliotou, C., Grupe, D., et al. 2006, ApJ, 642, 389 [NASA ADS] [CrossRef] (In the text)
- Ofek, E. O., Cenko, S. B., Gal-Yam, A., et al. 2007, ApJ, 662, 1129 [NASA ADS] [CrossRef]
- Padilla, L., Funk, B., Krawczynski, H., et al. 1998, A&A, 337, 43 [NASA ADS]
- Palmer, D., Cummings, J., Stamatikos, M., Markwardt, C., & Sakamoto, T. 2006, GCN Circ., 5076 [NASA ADS]
- van Paradijs, J., Kouveliotou, C., & Wijers, R. A. M. J., 2000, ARA&A, 38, 379 [NASA ADS] [CrossRef] (In the text)
- de Pasquale, M., Oates, S. R., Page, M. J., et al. 2007, MNRAS, 377, 1638 [NASA ADS] [CrossRef]
- Pe'er, A., & Waxman, E. 2005, ApJ, 633, 1018 [NASA ADS] [CrossRef]
- Pélangeon, A., Atteia, J.-L., Lamb, D. Q., & Ricker, G. R. 2006, AIPC, 836, 149 [NASA ADS]
- Perley, D. A., Li, W., Chornock, R., et al. 2008, ApJ, 688, 470 [NASA ADS] [CrossRef]
- Piran, T. 1999, Phys. Rep., 314, 575 [NASA ADS] [CrossRef]
- Poirier, J., D'Andrea, C., Fragile, P. C., et al. 2003, Phys. Rev. D, 67, 042001 [NASA ADS] [CrossRef]
- Razzaque, S., Mészáros, P., & Zhang, B. 2004, ApJ, 613, 1072 [NASA ADS] [CrossRef]
- Sakamoto, T., Lamb, D. Q., Kawai, N., et al. 2005, ApJ, 629, 311 [NASA ADS] [CrossRef]
- Sakamoto, T., Barthelmy, S. D., Barbier, L., et al. 2008, ApJS, 175, 179 [NASA ADS] [CrossRef]
- Sbarufatti, B., Barthelmy, S. D., Holland, S. T., et al. 2007, GCN Rep., 67 [NASA ADS].2 (In the text)
- Shirasaki, Y., Yoshida, A., Kawai, N., et al. 2008, PASJ, 60, 919 [NASA ADS]
- Soderberg, A. M., Kulkarni, S. R., Berger, E., et al. 2004, Nature, 430, 648 [NASA ADS] [CrossRef] (In the text)
- Soderberg, A. M., Kulkarni, S. R., Price, P. A., et al. 2006, ApJ, 636, 391 [NASA ADS] [CrossRef]
- Stamatikos, M., Barbier, L., Barthelmy, S. D., et al. 2007, GCN Circ., 6711 [NASA ADS]
- Stanek, K. Z., Matheson, T., Garnavich, P. M., et al. 2003, ApJ, 591, L17 [NASA ADS] [CrossRef]
- Wang, X. Y., Dai, Z. G., & Lu, T. 2001, ApJ, 556, 1010 [NASA ADS] [CrossRef] (In the text)
- Wang, X. Y., Li, Z., & Mészáros, P. 2006, ApJ, 641, L89 [NASA ADS] [CrossRef]
- Waxman, E., & Bahcall, J. N. 2000, ApJ, 541, 707 [NASA ADS] [CrossRef] (In the text)
- Zhang, B., & Mészáros, P. 2001, ApJ, 559, 110 [NASA ADS] [CrossRef]
- Zhang, B., & Mészáros, P. 2004, Int. J. Mod. Phys. A., 19, 2385 [NASA ADS] [CrossRef]
- Zhang, B., Fan, Y.-Z., Dyks, J., et al. 2006, ApJ, 642, 354 [NASA ADS] [CrossRef] (In the text)
Footnotes
- ...
- Supported by CAPES Foundation, Ministry of Education of Brazil.
- ... array
- http://www.mpi-hd.mpg.de/hfm/HESS/HESShtml
- ...
programme
- http://www.lsw.uni-heidelberg.de/projects/hess/HESS/grbs.phtml
- ... GCN
- The Gamma ray bursts Coordinates Network, http://gcn.gsfc.nasa.gov/
- ... indications
- Which include, e.g., burst position incompatible with known sources, and a high signal-to-noise ratio of the burst.
- ...
time
- HESS observations are taken in darkness and when the moon is below the horizon. The fraction of HESS dark time is about 0.2.
- ... circulars
- http://gcn.gsfc.nasa.gov/gcn3_archive.html
- ... cuts
- ``Soft'' cuts were called ``spectrum'' cuts in Aharonian et al. (2006a).
- ... analysis
-
is larger using soft-cut analysis.
- ... significance
- Calculated by Eq. (17) in Li & Ma (1983).
- ... value
- A larger excess implies a higher value of the upper limit, since the integrated exposure, that depends on ZA and
of the observations, is largely the same over the whole error box.
- ... here
- The optical depth of EBL absorption for a
100 GeV photon is
3 at z=1, according to the P0.45 model demonstrated in Aharonian et al. (2006b).
All Tables
Table 1: Properties of GRBs observed with HESS from March 2003 to October 2007.
Table 2: HESS observations of GRBs from March 2003 to October 2007.
Table 3: Differential flux upper limits at the energy thresholds from the HESS observations of GRBs with reported redshifts.
Table 4: Combined significance of 3 subsets of GRBs selected based on the requirements listed in Sect. 5.1.
All Figures
![]() |
Figure 1: Distribution of the statistical significance ( histogram) as derived from the observations of 21 GRBs using soft-cut analysis. The mean is -0.4 and the standard deviation is 1.4. Each entry corresponds to one GRB. The solid line is a Gaussian function with mean zero and standard deviation unity. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: The 99.9% confidence level energy flux upper limits (in red) at energies >200 GeV derived from HESS observations at the position of GRB 070621. The ends of the horizontal lines indicate the start and end times of the observations from which the upper limits were derived. The XRT energy flux in the 0.3-10 keV band is shown in black for comparison (Evans et al. 2007). |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
The |
Open with DEXTER | |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.