-
PDF
- Split View
-
Views
-
Cite
Cite
J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, K. Berger, E. Bernardini, A. Biland, O. Blanch, R. K. Bock, A. Boller, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, D. Carreto Fidalgo, P. Colin, E. Colombo, J. L. Contreras, J. Cortina, L. Cossio, S. Covino, P. Da Vela, F. Dazzi, A. De Angelis, G. De Caneva, B. De Lotto, C. Delgado Mendez, M. Doert, A. Domínguez, D. Dominis Prester, D. Dorner, M. Doro, D. Eisenacher, D. Elsaesser, E. Farina, D. Ferenc, M. V. Fonseca, L. Font, K. Frantzen, C. Fruck, R. J. García López, M. Garczarczyk, D. Garrido Terrats, M. Gaug, G. Giavitto, N. Godinović, A. González Munoz, S. R. Gozzini, A. Hadamek, D. Hadasch, A. Herrero, J. Hose, D. Hrupec, W. Idec, V. Kadenius, M. L. Knoetig, T. Krähenbühl, J. Krause, J. Kushida, A. La Barbera, D. Lelas, N. Lewandowska, E. Lindfors, S. Lombardi, R. López-Coto, M. López, A. López-Oramas, E. Lorenz, I. Lozano, M. Makariev, K. Mallot, G. Maneva, N. Mankuzhiyil, K. Mannheim, L. Maraschi, B. Marcote, M. Mariotti, M. Martínez, J. Masbou, D. Mazin, U. Menzel, M. Meucci, J. M. Miranda, R. Mirzoyan, J. Moldón, A. Moralejo, P. Munar-Adrover, D. Nakajima, A. Niedzwiecki, K. Nilsson, N. Nowak, R. Orito, A. Overkemping, S. Paiano, M. Palatiello, D. Paneque, R. Paoletti, J. M. Paredes, S. Partini, M. Persic, F. Prada, P. G. Prada Moroni, E. Prandini, S. Preziuso, I. Puljak, I. Reichardt, R. Reinthal, W. Rhode, M. Ribó, J. Rico, J. Rodriguez Garcia, S. Rügamer, A. Saggion, K. Saito, T. Saito, M. Salvati, K. Satalecka, V. Scalzotto, V. Scapin, C. Schultz, T. Schweizer, S. N. Shore, A. Sillanpää, J. Sitarek, I. Snidaric, D. Sobczynska, F. Spanier, V. Stamatescu, A. Stamerra, J. Storz, S. Sun, T. Surić, L. Takalo, F. Tavecchio, P. Temnikov, T. Terzić, D. Tescaro, M. Teshima, J. Thaele, O. Tibolla, D. F. Torres, T. Toyama, A. Treves, M. Uellenbeck, P. Vogler, R. M. Wagner, Q. Weitzel, F. Zandanel, R. Zanin, A. Bouvier, M. Hayashida, H. Tajima, F. Longo, MAGIC upper limits on the GRB 090102 afterglow, Monthly Notices of the Royal Astronomical Society, Volume 437, Issue 4, 01 February 2014, Pages 3103–3111, https://doi.org/10.1093/mnras/stt2041
- Share Icon Share
Abstract
Indications of a GeV component in the emission from gamma-ray bursts (GRBs) are known since the Energetic Gamma-Ray Experiment Telescope observations during the 1990s and they have been confirmed by the data of the Fermi satellite. These results have, however, shown that our understanding of GRB physics is still unsatisfactory. The new generation of Cherenkov observatories and in particular the MAGIC telescope, allow for the first time the possibility to extend the measurement of GRBs from several tens up to hundreds of GeV energy range. Both leptonic and hadronic processes have been suggested to explain the possible GeV/TeV counterpart of GRBs. Observations with ground-based telescopes of very high energy (VHE) photons (E > 30 GeV) from these sources are going to play a key role in discriminating among the different proposed emission mechanisms, which are barely distinguishable at lower energies. MAGIC telescope observations of the GRB 090102 (z = 1.547) field and Fermi Large Area Telescope data in the same time interval are analysed to derive upper limits of the GeV/TeV emission. We compare these results to the expected emissions evaluated for different processes in the framework of a relativistic blastwave model for the afterglow. Simultaneous upper limits with Fermi and a Cherenkov telescope have been derived for this GRB observation. The results we obtained are compatible with the expected emission although the difficulties in predicting the HE and VHE emission for the afterglow of this event makes it difficult to draw firmer conclusions. Nonetheless, MAGIC sensitivity in the energy range of overlap with space-based instruments (above about 40 GeV) is about one order of magnitude better with respect to Fermi. This makes evident the constraining power of ground-based observations and shows that the MAGIC telescope has reached the required performance to make possible GRB multiwavelength studies in the VHE range.
INTRODUCTION
Since the discovery of gamma-ray bursts (GRBs) in the late 1960s (Klebesadel, Strong & Olson 1973), these energetic and mysterious phenomena have been targets of large observational efforts. The discovery of their afterglow in late 1990s (Costa et al. 1997; Van Paradijs et al. 1997) provided a great boost in GRB studies at all wavelengths. The wealth of available information put severe constraints on the various families of interpretative scenarios, showing an unexpected richness and complexity of possible behaviours (see e.g. Gehrels, Ramirez-Ruiz & Fox 2009). The first observations at MeV–GeV energies with the Energetic Gamma-Ray Experiment Telescope (EGRET) on board the Compton Gamma-Ray Observatory (Hurley et al. 1994; Dingus 1995), showed that the high energy (HE: 1 MeV–30 GeV) and very high energy range (VHE: 30 GeV–30 TeV) can be powerful diagnostic tools for the emission processes and physical conditions of GRBs. The launch of Fermi (Band et al. 2009), with its Large Area Telescope (LAT; Atwood et al. 2009b), showed that, at least for the brightest events, GeV emission from GRBs is a relatively common phenomenon (Granot et al. 2010). However, a satisfactory interpretative framework of the GeV emission is still lacking. In this context, ground-based imaging atmospheric Cherenkov telescopes (IACTs), such as MAGIC,1 H.E.S.S.2 and VERITAS,3 despite the reduced duty cycle of ground-based facilities, provide access to the ∼100 GeV to TeV energy interval for GRB observations. Furthermore, the energy range down to ∼80 GeV, which was accessible almost exclusively with space-based instruments, has been opened to ground-based observations by the MAGIC observatory (Aliu et al. 2008; Schweizer et al. 2010). Together with the multiwavelength coverage provided by the LAT instrument, this makes possible the complete coverage of the 1–100 GeV energy range with the advantage, in the VHE domain, of an increase of ∼2–3 order of magnitude in the sensitivity relative to space-based instruments. Moreover, the low-energy trigger threshold of MAGIC makes less relevant the effect of the source distance. The flux above ∼100 GeV is attenuated by pair production with the lower energetic (optical/IR) photons of the diffuse Extragalactic Background Light (EBL; Nikishov 1962; Gould & Schreder 1966). The resulting cosmic opacity to VHE gamma-rays heavily affects Cherenkov observations, especially for GRBs which are sources with an average redshift slightly larger than 2 (Fynbo et al. 2009). Therefore, the higher the redshift, the lower the likelihood of detection at a given energy [i.e. about a 90 per cent of flux reduction at 100 GeV for a z = 2 source following Domínguez et al. (2011a)]. In addition, the transient and unpredictable nature of GRBs makes it difficult for large ground-based instruments such as IACTs to point them rapidly enough to catch the prompt emission and the early afterglow phases, when these sources are expected to be observable at high energies (see e.g. Covino et al. 2009a, for a discussion about IACTs perspectives for GRB observations). MAGIC has the advantage, compared to the other IACTs, in its low-energy sensitivity and pointing speed (e.g. Garczarczyk et al. 2009). Several attempts to observe GRB emission have been discussed (Albert et al. 2006, 2007; Aleksič et al. 2010). In all cases, only upper limits (ULs) have been derived. Similar results have also been reported by other IACTs (Tam et al. 2006; Aharonian et al. 2009; Acciari et al. 2011). As discussed above, the two most limiting factors are the high redshift of the source and the delay of the observation.
In this paper, we report and discuss the MAGIC observation of GRB 090102, a GRB at a redshift about 1.5 observed at low zenith angle and good weather conditions. These observations permitted data-taking with an energy threshold of about 30 GeV. However, no gamma-ray signal was detected, and hence only ULs could be derived.
Section 2 gives general information about GRB 090102. In Sections 3 and 4, we discuss the MAGIC and LAT data sample, respectively, while Sections 5 and 6 introduce and develop the interpretative scenario. In Section 7, we evaluate the effect of the EBL absorption on the lowest energy bins allowed by our observation and finally, we discuss our results in a general theoretical scenario in the last section.
We assume a Λ cold dark matter cosmology with Ωm = 0.27, |$\Omega _\Lambda = 0.73$| and h0 = 0.71. At the redshift of the GRB, the proper distance is ∼4.5 Gpc (∼1.38 × 1028 cm). Throughout this paper the convention Qx = Q/10x has been adopted in CGS units.
GRB 090102
GRB 090102 was detected and located by the Swift satellite (Gehrels et al. 2004) on 2009 January 2 at 02:55:45 ut (Mangano et al. 2009b) and also by the Fermi Gamma-ray Burst Monitor (GBM) detector. The prompt light curve was structured in four partially overlapping peaks (Sakamoto et al. 2009) for a total T90 of 27.0 ± 2.0 s. Since the burst was also detected by Konus Wind (Golenetskii et al. 2009) and Integral (Mangano et al. 2009a), it has been possible to obtain a very good reconstruction of the prompt emission spectral parameters. The time-averaged spectrum can be modelled with the classical Band function (Band et al. 1993) with peak energy Epeak = 451|$^{ \ +73}_{ \ -58}$| keV and a total fluence in the 20 keV–2 MeV range of 3.09|$^{ \ +0.29}_{ \ -0.25} \times 10^5$| erg cm−2 (Golenetskii et al. 2009). Early optical follow-up measurements were performed by many groups like TAROT (Klotz et al. 2009) at T0+40.8 s, the REM robotic telescope at T0+53 s (Covino, D'Avanzo & Antonelli 2009b) and GROND telescope (Afonso et al. 2009) at T0+2.5 h. Optical spectroscopy was rapidly obtained with the NOT telescope by De Ugarte Postigo et al. (2009a). They found evidence of several absorption metal lines, including Fe ii, Mg ii, Mg i, Al ii, Al iii and C iv, at a common redshift of z = 1.547. The resulting isotropic energy value Eiso = 5.75 × 1053 erg and the rest-frame peak energy Epeak = 1149|$^{ \ +186}_{ \ -148}$| keV are in good agreement with the Amati relation (Amati, Frontera & Guidorzi 2009). The multiwavelength light curve is shown in Fig. 1 in which, data in the R and H band correspond to a rest-frame UV and optical emission, respectively. X-ray data are the unbinned Swift X-Ray Telescope (XRT) and Burst Alert Telescope (BAT) data in the 0.5–10 keV. According to Gendre et al. (2010), it is very difficult to model the whole afterglow in a standard scenario (see the next section). Moreover, it showed a distinct behaviour in the optical and in X-rays. The X-ray light curve showed an uninterrupted decay from about 400 s from the GRB onset up to 5 × 106 s, when Swift ceased observations of the event. The optical light curve, monitored from several tens of seconds to slightly more than a day from the T0, showed a steep-to-shallow behaviour with a break at about 1 ks. Before the break, the optical flux decay index is α1 = 1.50 ± 0.06 while the index becomes α2 = 0.97 ± 0.03 after the break, steeper and flatter, respectively, when compared to the simultaneous X-ray emission. This behaviour strongly resembles that showed by GRB 061126 (Gomboc et al. 2008; Perley et al. 2008) and GRB 060908 (Covino et al. 2010).
![Light curve for GRB 090102. Data in the R band (red points) were taken from the TAROT, REM, NOT, GROND and Palomar telescopes. The early near-infrared H band (blue points) are from the observations of the REM telescope. All magnitudes are expressed in the Vega system. X-ray data are the unbinned Swift/XRT and BAT data in the 0.5–10 keV (green and magenta point, respectively). The MAGIC observation window is also plotted. R and H data from Gendre et al. (2010). XRT and BAT data retrieved from http://www.swift.ac.uk/burst_analyser (Evans et al. 2010).](https://cdn.statically.io/img/oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/mnras/437/4/10.1093_mnras_stt2041/4/m_stt2041fig1.jpeg?Expires=1723624766&Signature=Xe0K3Xc5-IUbIp85d3SSrv1s4pZIXICo6j9vmWjHKPhB0aPfyT~MwUKgpywd9Rth24lzPNN3Lu20x4cIMGBUY1ZemGS8stYfw~qDJvsqAUM6BNqJTMeJIwjmKXO-ztXqJP7FuU4v4xfLrcnag1ErDIejhWrU~dSQYXvgy78UJOorCacigMJ84sFwvwWQBDdRJ~-NBdsZelNOLbdKJNBPt4LRMINjKIMbBdtUYbiUgNZjEsukcxGqqNTNeL~7L10S6lgP9anusTPpO4Pr6WXQUH~eQFGu49kVRhJRkDAEHAucxvfdNvyxCIj86S9goyOSrmRFKI-iRg-4C3HtsQuGmw__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Light curve for GRB 090102. Data in the R band (red points) were taken from the TAROT, REM, NOT, GROND and Palomar telescopes. The early near-infrared H band (blue points) are from the observations of the REM telescope. All magnitudes are expressed in the Vega system. X-ray data are the unbinned Swift/XRT and BAT data in the 0.5–10 keV (green and magenta point, respectively). The MAGIC observation window is also plotted. R and H data from Gendre et al. (2010). XRT and BAT data retrieved from http://www.swift.ac.uk/burst_analyser (Evans et al. 2010).
Other observations were performed at later time by GROND (T0+2.5 h; Afonso et al. 2009), Palomar (T0+50 min; Cenko, Rau & Salvato 2009) and IAC80 (T0+19.2 h; De Ugarte Postigo, Blanco & Castro-Tirando 2009b) telescopes while during the following days, the NOT (Malesani et al. 2009) and HST (Levan et al. 2009) provided the detection of the host galaxy. In the radio energy band, the VLA (Chandra & Frail 2009) and the Westerbork Synthesis Radio Telescope (Van der Horst, Wijers & Kamble 2009) performed follow-up observation at 8.46 and 4.9 GHz with no afterglow detection and UL evaluation. A detailed discussion of the follow-up observations for this burst can be found in Gendre et al. (2010).
MAGIC FOLLOW-UP OBSERVATION AND ANALYSIS
The MAGIC telescope located at Roque de los Muchachos (28| $_{.}^{\circ}$|75 N, 17| $_{.}^{\circ}$|89 W, La Palma, Canary Islands) performed a follow-up measurement of GRB 090102. The data presented in this paper were taken when MAGIC was operating as a single telescope. The MAGIC telescope was autonomously repointed and started the observations at T0+255 s, following the GRB alert from Fermi-GBM. Later on, the shift crew operating the telescope realized that the GBM coordinates (RA: 08h35m06s; Dec.: 37°16′48′′) differed from the BAT coordinates (RA: 08h33m02s; Dec.: 33°05′29′′) by more than 4°. Consequently, the telescope was repointed to the BAT coordinates and re-started observations by T0+1161 s. After this burst, the alert system was modified to cope with this situation. First data runs were taken at very low zenith angles from 5° reaching 52° at the end of data taking at 06:54:01 ut after 13 149 s of observation. MAGIC ULs above 80 GeV have already been published for this GRB (Gaug et al. 2009a), while results and scientific discussion about a subsequent dedicated analysis focused in the low-energy band (Gaug et al. 2009b) will be presented here. To ensure the lowest energy threshold, only data taken with zenith distance <25°, corresponding to the first 5919 s of observation (data subsample up to 04:53:32 ut) have been taken into account during this analysis. By employing the MAGIC-1 sum trigger system (Aliu et al. 2008), an analysis threshold of around 30 GeV is achieved, which is evaluated from Monte Carlo (MC) simulations. In order to accurately estimate the background from hadronic atmospheric showers, an OFF data sample was taken one night later with the telescope pointing close to the burst location and in the same observational conditions and instrument setup. Data were analysed using the MAGIC Analysis and Reconstruction Software (mars; Albert et al. 2008; Aliu et al. 2009a) and processed using the standard Hillas parameters (Hillas 1985). Gamma/hadron separation and energy estimation were performed using a multidimensional classification method (Random Forest; Breiman 2001) while arrival directions of the gamma photons is reconstructed using the disp algorithm (Fomin et al. 1994). The alpha parameter is then used to evaluate the significance of the signal in six energy bins. In spite of the low-energy analysis threshold, no significant excess of gamma-ray photons have been detected from a position consistent with GRB 090102. Differential ULs assuming a power-law gamma-ray spectrum with spectral index of Γ = −2.5 and using the method of Rolke, López & Conrad (2005) were evaluated with a 95 per cent confidence level (CL) and 30 per cent estimation of systematic uncertainties and are reported in Table 1 and Fig. 2.
![SSC modelled emission during the afterglow of GRB 090102. Blue triangles are 95 per cent CL ULs derived by MAGIC for low-energy (LE) analysis. The relatively more constraining UL in the 50–80 GeV is due to a negative significance energy bin. For comparison, the regular energy range MAGIC ULs (Gaug et al. 2009a) are also reported in light grey. The red triangles report the Fermi-LAT 95 per cent CL ULs. The purple and black curves depict the expected energy flux according to the GRB afterglow model described in Sections 6 and 5. Physical parameters are ϵe = 0.1, ϵB = 0.01, E52 = 4.5 and T = T0 + 4 ks at a redshift z = 1.547. The shaded region shows the uncertainty in the EBL absorption, as prescribed in Domínguez et al. (2011a).](https://cdn.statically.io/img/oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/mnras/437/4/10.1093_mnras_stt2041/4/m_stt2041fig2.jpeg?Expires=1723624766&Signature=AEaHZkEMJjBcpU2xq~S~X5BkwknRvoN39k6G8iSIsO5BMakfU10SrI2bqyB44KXLBMUVbW-fnENMfx3Kx-G9H6i~VP-Tj22lGQ1EUGMn0CxoCGmb-8Qyc9ssCkW2X5~SZhYNlnEgfCG6rSVQQXdYTcIKmgbecQ35fY595HAlcurkw9DQAoNG7G8Bco04yUxOmhhGPjCGx-TJ~u1dY6yQ1~xgi6nRtO5i7HQIpTmk6pq0l7XbKEWALxF4YV3buBurgJkQ-g7tRMGHoj~3XawycE5fMTbbNvYIf~N8HuK7IqTXL1Kb6Hn734Ue0mp-NUP2J~P8idNjhW~xqGzJ4LFzew__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
SSC modelled emission during the afterglow of GRB 090102. Blue triangles are 95 per cent CL ULs derived by MAGIC for low-energy (LE) analysis. The relatively more constraining UL in the 50–80 GeV is due to a negative significance energy bin. For comparison, the regular energy range MAGIC ULs (Gaug et al. 2009a) are also reported in light grey. The red triangles report the Fermi-LAT 95 per cent CL ULs. The purple and black curves depict the expected energy flux according to the GRB afterglow model described in Sections 6 and 5. Physical parameters are ϵe = 0.1, ϵB = 0.01, E52 = 4.5 and T = T0 + 4 ks at a redshift z = 1.547. The shaded region shows the uncertainty in the EBL absorption, as prescribed in Domínguez et al. (2011a).
MAGIC-I 95 per cent confidence level ULs for the afterglow emission of GRB 090102. The values correspond to the first 5919 s of observation from 03:14:52 to 04:53:32 ut. α Bins central energy was evaluated applying all analysis cuts to MC simulations. β Statistical significance of the excess events observed by MAGIC.
E bin . | 〈E〉 α . | σ β . | Average flux limits . |
---|---|---|---|
(GeV) . | (GeV) . | . | (erg cm−2 s−1) . |
25–50 | 43.9 | 0.83 | 8.7 × 10−10 |
50–80 | 57.3 | −0.30 | 1.5 × 10−10 |
80–125 | 90.2 | 1.09 | 3.1 × 10−10 |
125–175 | 137.2 | 0.51 | 2.2 × 10−10 |
175–300 | 209.4 | 0.90 | 1.6 × 10−10 |
300–1000 | 437.6 | −0.48 | 0.3 × 10−10 |
E bin . | 〈E〉 α . | σ β . | Average flux limits . |
---|---|---|---|
(GeV) . | (GeV) . | . | (erg cm−2 s−1) . |
25–50 | 43.9 | 0.83 | 8.7 × 10−10 |
50–80 | 57.3 | −0.30 | 1.5 × 10−10 |
80–125 | 90.2 | 1.09 | 3.1 × 10−10 |
125–175 | 137.2 | 0.51 | 2.2 × 10−10 |
175–300 | 209.4 | 0.90 | 1.6 × 10−10 |
300–1000 | 437.6 | −0.48 | 0.3 × 10−10 |
MAGIC-I 95 per cent confidence level ULs for the afterglow emission of GRB 090102. The values correspond to the first 5919 s of observation from 03:14:52 to 04:53:32 ut. α Bins central energy was evaluated applying all analysis cuts to MC simulations. β Statistical significance of the excess events observed by MAGIC.
E bin . | 〈E〉 α . | σ β . | Average flux limits . |
---|---|---|---|
(GeV) . | (GeV) . | . | (erg cm−2 s−1) . |
25–50 | 43.9 | 0.83 | 8.7 × 10−10 |
50–80 | 57.3 | −0.30 | 1.5 × 10−10 |
80–125 | 90.2 | 1.09 | 3.1 × 10−10 |
125–175 | 137.2 | 0.51 | 2.2 × 10−10 |
175–300 | 209.4 | 0.90 | 1.6 × 10−10 |
300–1000 | 437.6 | −0.48 | 0.3 × 10−10 |
E bin . | 〈E〉 α . | σ β . | Average flux limits . |
---|---|---|---|
(GeV) . | (GeV) . | . | (erg cm−2 s−1) . |
25–50 | 43.9 | 0.83 | 8.7 × 10−10 |
50–80 | 57.3 | −0.30 | 1.5 × 10−10 |
80–125 | 90.2 | 1.09 | 3.1 × 10−10 |
125–175 | 137.2 | 0.51 | 2.2 × 10−10 |
175–300 | 209.4 | 0.90 | 1.6 × 10−10 |
300–1000 | 437.6 | −0.48 | 0.3 × 10−10 |
Spectrum break energies for the different considered processes. The value of the expected SSC emission in the first MAGIC energy bin (∼40 GeV) is also shown with and without considering the EBL absorption. We refer to Zhang & Mésźaros (2001) for the numerical results presented in this paper.
Synchrotron (e) . | SSC . | Synchrotron (p) . |
---|---|---|
Em ≃ 0.6 eV | |$E_{\rm m}^{\rm ssc} \simeq 1.1$| MeV | |$E_{\rm m}^p \simeq 10^{-8}$|eV |
Ec ≃ 4.1 eV | |$E_{\rm c}^{\rm {\rm ssc}} \simeq 47$| MeV | |$E_{\rm c}^p \simeq 140$| TeV |
Emax ≃ 207 MeV | EKN ≃ 60 TeV | |$E_{\rm max}^p \simeq$| 1.7 MeV |
≈5 × 10−11 | ≈1.1 × 10−10 | ≈4 × 10−17 |
– | 4.3 × 10−11 (3.4 × 10−11) | – |
Synchrotron (e) . | SSC . | Synchrotron (p) . |
---|---|---|
Em ≃ 0.6 eV | |$E_{\rm m}^{\rm ssc} \simeq 1.1$| MeV | |$E_{\rm m}^p \simeq 10^{-8}$|eV |
Ec ≃ 4.1 eV | |$E_{\rm c}^{\rm {\rm ssc}} \simeq 47$| MeV | |$E_{\rm c}^p \simeq 140$| TeV |
Emax ≃ 207 MeV | EKN ≃ 60 TeV | |$E_{\rm max}^p \simeq$| 1.7 MeV |
≈5 × 10−11 | ≈1.1 × 10−10 | ≈4 × 10−17 |
– | 4.3 × 10−11 (3.4 × 10−11) | – |
Spectrum break energies for the different considered processes. The value of the expected SSC emission in the first MAGIC energy bin (∼40 GeV) is also shown with and without considering the EBL absorption. We refer to Zhang & Mésźaros (2001) for the numerical results presented in this paper.
Synchrotron (e) . | SSC . | Synchrotron (p) . |
---|---|---|
Em ≃ 0.6 eV | |$E_{\rm m}^{\rm ssc} \simeq 1.1$| MeV | |$E_{\rm m}^p \simeq 10^{-8}$|eV |
Ec ≃ 4.1 eV | |$E_{\rm c}^{\rm {\rm ssc}} \simeq 47$| MeV | |$E_{\rm c}^p \simeq 140$| TeV |
Emax ≃ 207 MeV | EKN ≃ 60 TeV | |$E_{\rm max}^p \simeq$| 1.7 MeV |
≈5 × 10−11 | ≈1.1 × 10−10 | ≈4 × 10−17 |
– | 4.3 × 10−11 (3.4 × 10−11) | – |
Synchrotron (e) . | SSC . | Synchrotron (p) . |
---|---|---|
Em ≃ 0.6 eV | |$E_{\rm m}^{\rm ssc} \simeq 1.1$| MeV | |$E_{\rm m}^p \simeq 10^{-8}$|eV |
Ec ≃ 4.1 eV | |$E_{\rm c}^{\rm {\rm ssc}} \simeq 47$| MeV | |$E_{\rm c}^p \simeq 140$| TeV |
Emax ≃ 207 MeV | EKN ≃ 60 TeV | |$E_{\rm max}^p \simeq$| 1.7 MeV |
≈5 × 10−11 | ≈1.1 × 10−10 | ≈4 × 10−17 |
– | 4.3 × 10−11 (3.4 × 10−11) | – |
LAT OBSERVATION AND ANALYSIS
The Fermi observatory is operating in a sky survey mode and the Swift localization of GRB 090102 was observable by the LAT instrument approximately 3300 s after trigger and remained within the LAT field of view (Θboresight ≲ 60°) for a duration of ∼2300 s. We analysed the Fermi-LAT data using the Science Tools 09-30-01 with Pass7V6 ‘Source’ event class. We used the publicly available models for the Galactic and isotropic diffuse emissions, |$gal\_2yearp7v6\_trim\_v0.fits$| and |$iso\_p7v6source.txt$|, that can be retrieved from the Fermi Science Support Center.4 No significant excess was found in this observation, so we computed ULs in three different energy bands: [0.1–1 GeV], [1–10 GeV] and [10–100 GeV]. We first fit the broad energy range (from 0.1 to 100 GeV) using the unbinned likelihood analysis, which was then used to constrain the background model. Then we froze the normalizations of the isotropic and Galactic diffuse templates, and independently fit the source in the three different energy bands, using the unbinned profile likelihood method to derive 95 per cent LAT ULs. The following UL values were derived for the [0.1–1 GeV], [1–10 GeV], [10–100 GeV] energy ranges, respectively: 2.73 × 10−10, 4.58 × 10−10, 3.45 × 10−9 erg cm−2 s−1 and are depicted in Fig. 2. These ULs are more constraining than the ones reported in Inoue et al. (2013). The reason for that is the usage of P7V6 ‘Source’ instead of P6V3 ‘Diffuse’, and also the usage of a different procedure to parametrize the diffuse background in the three differential energy bins. Even if observed with a considerable time delay, the achieved energy threshold of MAGIC permits a better overlap with LAT in the GeV range when compared with previous results on GRB by MAGIC and other IACTs. Thus, it has been possible to derive simultaneous ULs with a complete coverage of the energy range from 0.1 GeV up to TeV using MAGIC and Fermi-LAT. Furthermore, it is worth stressing that, in the energy range where the two instruments overlap (range [25–100 GeV]), the ULs derived by MAGIC are about one order of magnitude lower than those from Fermi-LAT.
THE LOW-ENERGY SCENARIO
MODELING THE VHE EMISSION
Basically, the new spectral feature has the same shape of the underlying synchrotron component with a new break in the spectrum (EKN) due to the decreasing of the IC cross-section with energy (Fragile et al. 2004). However, this cut-off is found to be above few tens of TeV in our case, securing that MAGIC ULs stay below this limit. The relevant break energies for the assumed model are summarized in Table 2.
EBL ATTENUATION
Gamma-ray absorption by pair production with EBL plays a key role in VHE astronomy since it significantly limits the IACTs capability in detecting sources at redshift z > 1. The optical depth τ is strictly connected to the light content of the Universe and the source distance. In the past years, several EBL models have been proposed providing a wide range of values for τ from 1 up to 6 for a z ∼ 1 source at 100 GeV (see e.g. Kneiske et al. 2004; Stecker, Malkan & Scully 2006; Stecker & Scully 2008), which gives an attenuation in the expected flux ranging between 1/3 and to more than 1/100. However, the more recent EBL models (Franceschini, Rodighiero & Vaccari 2008; Domínguez et al. 2011a), although based on different assumptions, are converging to stable results. Within this context, Domínguez et al. (2011a) have used real data on the evolution of the galaxy population taken from the All-wavelength Extended Groth Strip International Survey (AEGIS) catalogue to evaluate EBL intensity for a wide range of redshift. The reliability of the results have been tested on the three most distant object observed by MAGIC (Domínguez et al. 2011b). Moreover, the EBL intensity evaluated using this model matches the minimum level allowed by galaxy counts which leads to the highest transparency of the universe to VHE gamma-rays. We used the model of Domínguez et al. (2011a) to evaluate the EBL absorption obtaining a value for τ of |$0.218^{ +0.075}_{ -0.041}$| at about 40 GeV. This gives an attenuation of the flux at the same energy of the order of ∼20 per cent, a value that does not significantly compromise detection capability of MAGIC. However, the optical depth increases quickly with energy reaching the values of ∼1.5 and ∼14.4 at 100 and 500 GeV, respectively, and this makes necessary to lower the energy threshold of the observation. In the case of GRB 090102, MAGIC shows its capability to perform observation at very low energy limiting the gamma absorption even for moderate redshift sources.
DISCUSSION
Although only ULs have been obtained, the possibility of having simultaneous observations with Fermi-LAT in the energy range 0.1–100 GeV and MAGIC in the energy range that starts at 25 GeV (hence overlapping with LAT) make the GRB 090102 a good case study in spite of its relatively high redshift. However, it has to be remarked that GRB 090102 can be considered as a common GRB in terms of both energetics and redshifts. Higher expected fluxes can be foreseen in the case of more energetic events that are not so rare accordingly to Fermi results. We have used the equations of relativistic shock model in order to predict, in a reliable way, the expected VHE emission in the LAT and MAGIC energy range. From numerical results, it is evident that for the chosen parameters and at our observation time, leptonic components are the dominant mechanisms from the low to the VHE. Following equation (6), the cooling frequency for protons is usually located well above the VHE range (>TeV) and this make the process potentially interesting for MAGIC observations. However, to make the two emissions comparable in the low-energy regime, and the proton component to dominate the leptonic one at high energies, a fine tuning in the parameters choice is needed implying |$\frac{\epsilon _{\rm e}}{\epsilon _{\rm p}} \approx \frac{m_{\rm e}}{m_{\rm p}} \approx 10^{-3}$|. Similar results can be obtained for the IC component. In both cases, however, a higher total energy release of ≈1055 erg and a circumburst density medium of ≈100 cm−3 are needed to maintain the low-energy flux at the observed level. This makes the possibility of observing the hadronic emission component with the MAGIC telescope unrealistic, at least for a canonical model. A sketch of the scenario described above is shown in Fig. 3.
![Modeled SED in hadronic-dominated scenario for GRB afterglow. Used parameters are E52 = 103, Tobs = T0 + 4 ks, ϵe = 10−3, ϵB = 0.01, n = 100 cm−3.](https://cdn.statically.io/img/oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/mnras/437/4/10.1093_mnras_stt2041/4/m_stt2041fig3.jpeg?Expires=1723624766&Signature=hmiltieSfOnbMbutvlLrvqWRU8~j7OQHJYdYvOtg9jdQFOtfxxqzZiNS6nYSvCNXuu9QXudIL70e4MH1HVXiPXa6NZzOBPkC7Af42FgWZoQ0PVtSrBeIDTe5GkSY6klmA7ktzyEvQs48bBV2nN~fznFQ54jJi9Tgdn11gNpX1LR2EdXvLBSISjidaRqzeWdUuh2lL-BHTXnC6xl0segj-mremPMHIFIxwDKELXsViZx6MvC7EUVxF0bizy8~kgvg9724akVp-MhE1Q7Tw6mfxk5pH6mw~hw~ke0s0YMC0GJQCCHlYvkUsKjUyUl45gY2S02nYfH~lYaf0cC0R0X8Cg__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Modeled SED in hadronic-dominated scenario for GRB afterglow. Used parameters are E52 = 103, Tobs = T0 + 4 ks, ϵe = 10−3, ϵB = 0.01, n = 100 cm−3.
![The expected SSC emission at different time. T0 + 0.8 ks (purple curve), T0 + 2 ks (yellow curve), T0 + 10 ks (orange curve). In the latter case, the corresponding VERITAS 99 per cent ULs on ∼5000 s of observation have been plotted (Acciari et al. 2011) for comparison.](https://cdn.statically.io/img/oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/mnras/437/4/10.1093_mnras_stt2041/4/m_stt2041fig4.jpeg?Expires=1723624766&Signature=yGrcDXVRJEZMmb3GSW4oVLsX4kPh7nywHhzcyC7M0aLR6BlRzsOtppr2R7a0GEpfd36i5ygymlZHMKVpm9EXjZRWOQw74dot1qkA0FyQ-M4tfbcjUWD36P3DtKqaBHrPeZTrEW7YOmN8ux7XjRKmQnopa6KFNG1fwP5jrNNiqVWW0dyNpRrzqUtE9GMNJkpfaBZWEuroWfu9pSMcbhBXfUh7mOiA5PjxNSjKsjaFdDiqLm4x2Q5GF4nYCdWaGfYIIjOy7-EWJ01IPgO~hyehC6srjWq2dD1at3So15KiqW-oi5izFeRmy3lnEZlHJoL7mg-jI1ZfMF3-qYpMBLyC5Q__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
The expected SSC emission at different time. T0 + 0.8 ks (purple curve), T0 + 2 ks (yellow curve), T0 + 10 ks (orange curve). In the latter case, the corresponding VERITAS 99 per cent ULs on ∼5000 s of observation have been plotted (Acciari et al. 2011) for comparison.
FUTURE PROSPECTS
Catching VHE signal from GRBs is one of the primary target of the MAGIC telescope and future IACTs like the Cherenkov Telescope Array (CTA). Our estimates show that for this particular GRB, MAGIC follow-up observations made within the first 1–2 min from the trigger time would have the potential to detect the VHE component or at least to evaluate constraining ULs (see Fig. 5). This demonstrates both the capabilities of the system and the necessity of a fast-response observations.
![Integrated σLiMa significance as a function of observation time in the 63–158 GeV for the MAGIC stereo system in the case of a GRB event similar to the one reported in this paper. Significance curve evolutions are showed for different starting of observation times after GRB onset: 180 s (green), 600 s (blue), 1100 s (red). In this latter case (inner plot), the foreseen performance for CTA assuming the preliminary sensitivity achieved with MC simulations in the same energy range is showed. The coloured areas shows the assumed 50 per cent systematic errors in the effective area evaluation as explained in Lombardi, Carosi & Antonelli (2013).](https://cdn.statically.io/img/oup.silverchair-cdn.com/oup/backfile/Content_public/Journal/mnras/437/4/10.1093_mnras_stt2041/4/m_stt2041fig5.jpeg?Expires=1723624766&Signature=S45iE42PpleKle1CK6Onfq-WggXauAaZTnS~z6KKd7FeAUQKcLuw35WMtQ6I2rzTg7wlihPUqCQF4DTsZrVKIdsJInwxCDcf1mZ57g-OZW9FQTn-sVrm6HR2xs6bl-91Dh2BwoeV2z5Pyc4p5LcS853~nZla3~a2QRHbht9h9NXKH-7WJafnUhtOpH~VFLmCRbT3W5MNC0-KJT7PRnuHIHVVukBUxIsiJJgcqPKV9YEh9f6VqGUq9wM0YTwoC5zadKEXglKj9TCrpLbE6j9zlUnzdDLTE5Y85GElLZE7c6kLLqQYY6BcNtvXqo7Q833yHj2J1YkrE9umNZ5wJKGHQQ__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Integrated σLiMa significance as a function of observation time in the 63–158 GeV for the MAGIC stereo system in the case of a GRB event similar to the one reported in this paper. Significance curve evolutions are showed for different starting of observation times after GRB onset: 180 s (green), 600 s (blue), 1100 s (red). In this latter case (inner plot), the foreseen performance for CTA assuming the preliminary sensitivity achieved with MC simulations in the same energy range is showed. The coloured areas shows the assumed 50 per cent systematic errors in the effective area evaluation as explained in Lombardi, Carosi & Antonelli (2013).
As GeV emission is found to be relatively common in Fermi GRBs (see e.g Abdo et al. 2009c), the unique opportunity of having simultaneous follow-up with LAT and the MAGIC telescope will make accessible the end of the electromagnetic spectrum of GRBs and will have an important role in constraining different emission mechanisms and the space parameters. Moreover, the recent technical improvement of the MAGIC stereo system (Aleksič et al. 2012) will bring an improvement in the instrument sensitivity in its low-energy range. The steeper decay of the flux makes in any case difficult late time (>200 s) detections for such moderate high-redshift event (see Fig. 5). On the other hand, such a time-scale is well within the pointing capabilities of the present generation of IACTs (e.g. MAGIC) that are able to perform follow-up measurements within few hundreds of seconds. Basing on the preliminary sensitivity of the future CTA,6 a detection will instead be possible, within the assumed model, even on later time (>1000 s) and higher redshift events.
We acknowledge an anonymous referee for useful comments.
We would like to thank the Instituto de Astrofísica de Canarias for the excellent working conditions at the Observatorio del Roque de los Muchachos in La Palma. The support of the German BMBF and MPG, the Italian INFN, the Swiss National Fund SNF and the Spanish MICINN is gratefully acknowledged. This work was also supported by the CPAN CSD2007-00042 and MultiDark CSD2009-00064 projects of the Spanish Consolider-Ingenio 2010 programme, by grant DO02-353 of the Bulgarian NSF, by grant 127740 of the Academy of Finland, by the DFG Cluster of Excellence ‘Origin and Structure of the Universe’, by the DFG Collaborative Research Centers SFB823/C4 and SFB876/C3 and by the Polish MNiSzW grant 745/N-HESS-MAGIC/2010/0.
The Fermi-LAT Collaboration acknowledges support from a number of agencies and institutes for both development and the operation of the LAT as well as scientific data analysis. These include NASA and DOE in the United States; CEA/Irfu and IN2P3/CNRS in France; ASI and INFN in Italy; MEXT, KEK and JAXA in Japan; and the K. A. Wallenberg Foundation, the Swedish Research Council and the National Space Board in Sweden. Additional support from INAF in Italy and CNES in France for science analysis during the operations phase is also gratefully acknowledged. This research is partially supported by NASA through the Fermi Guest Investigator Grants NNX09AT92G and NNX10AP22G. This work made use of data supplied by the UK Swift Science Data Centre at the University of Leicester.
Present address: Ecole polytechnique fédérale de Lausanne (EPFL), Lausanne, Switzerland.
Present address: Finnish Centre for Astronomy with ESO (FINCA), University of Turku, Finland.
Present address: GRAPPA Institute, University of Amsterdam, NL-1098 XH Amsterdam, the Netherlands.
GRBs show their phenomenology mainly in the X-ray and soft γ-ray energy band (1 keV–1 MeV). To avoid confusion with the Fermi-LAT and IACT operational energy range (>20 MeV and >25 GeV, respectively), we will refer to the former as a ‘low-energy’ range.