Jan 6, 2007 - arXiv:astro-ph/0610660v2 6 Jan 2007. The Remarkable Afterglow of GRB 061007: ...... by A. Wilson. ESA SP-604, Volume 2, Noordwijk: ESA ...
arXiv:astro-ph/0610660v2 6 Jan 2007
The Remarkable Afterglow of GRB 061007: Implications for Optical Flashes and GRB Fireballs C. G. Mundell1 , A. Melandri1 , C. Guidorzi1,2,3 , S. Kobayashi1 , I. A. Steele1 , D. Malesani4 , L. Amati5 , P. D’Avanzo3,6 , D. F. Bersier1 , A. Gomboc1,7 , E. Rol8 , M. F. Bode1 , D. Carter1 , C. J. Mottram1 , A. Monfardini1,9 , R. J. Smith1 , S. Malhotra10 , J. Wang11 , N. Bannister8 , P. T. O’Brien8 , N. R. Tanvir8 . (cgm)@astro.livjm.ac.uk ABSTRACT We present a multiwavelength analysis of Swift GRB 061007. The 2-m robotic Faulkes Telescope South (FTS) began observing 137 s after the onset of the γ-ray emission, when the optical counterpart was already decaying from R∼10.3 mag, and continued observing for the next 5.5 hours. These observations begin during the final γ-ray flare and continue through and beyond a long, soft tail of γ-ray emission whose flux shows an underlying simple power-law decay identical to that seen at optical and X-ray wavelengths, with temporal slope α ∼ 1.7 (F ∝ t−α ). This remarkably simple decay in all of these bands is rare for Swift bursts, which 1
Astrophysics Research Institute, Liverpool John Moores University, Twelve Quays House, Birkenhead, CH41 1LD, UK 2
Dipartimento di Fisica, Universit´a di Milano-Bicocca, Piazza delle Scienze 3, 20126 Milano, Italy
3
INAF-Osservatorio Astronomico di Brera, via Bianchi 46, 23807 Merate (LC), Italy
4
International School for Advanced Studies (SISSA-ISAS), via Beirut 2-4, I-34014 Trieste, Italy
5
INAF - IASF Bologna, via P. Gobetti 101, Bologna, Italy
6
Universit`a dell’Insubria, Dipartimento di Fisica e Matematica, via Valleggio 11, I-22100 Como, Italy
7
FMF, University in Ljubljana, Jadranska 19, 1000 Ljubljana, Slovenia.
8
Department of Physics and Astronomy, University of Leicester, University Road, Leicester LE1 7RH,
UK 9
CNRS-CTBT, 25 Avenue des Martyrs, 38000 Grenoble, France
10
School of Earth and Space Exploration, Arizona State University, P.O. Box 871404, Tempe, AZ 852871404 11
Center for Astrophysics, University of Science and Technology of China, Hefei, Anhui 230026, P. R. China
–2– often show much more complex light curves. We suggest the afterglow emission begins as early as 30 − 100 s and is contemporaneous with the on-going variable prompt emission from the central engine, but originates from a physically distinct region dominated by the forward shock. The afterglow continues unabated until at least ∼105 seconds showing no evidence of a break. The observed multiwavelength evolution of GRB 061007 is explained by an expanding fireball whose optical, X-ray and late-time γ-ray emission is dominated by emission from a forward shock with typical synchrotron frequency, νm , that is already below the optical band as early as t=137 s and a cooling frequency, νc , above the X-ray band to at least t=105 s. In contrast, the typical frequency of the reverse shock lies in the radio band at early time. We suggest that the unexpected lack of bright optical flashes from the majority of Swift GRBs may be explained with a low νm originating from small microphysics parameters, ǫe and ǫB . We derive a minimum jet opening angle θ = 4.7◦ from the optical light curves and conclude that GRB 061007 is a secure outlier to spectral energy correlations because no X-ray jet break occurred before t ∼ 106 s. Subject headings: gamma rays: bursts – shock waves – radiation mechanisms: nonthermal – cosmology: observations
1.
Introduction
Gamma-ray bursts (GRBs) are the most instantaneously luminous cosmological objects known, and although their transient nature makes them observationally challenging to follow up at other wavelengths after the initial burst of γ rays, their brightness and large redshift range (0.0085 < z < 6.29) ensure they remain uniquely useful as probes of the early Universe. With the launch of the GRB-dedicated Swift satellite (Gehrels et al. 2004), early-time monitoring of the radiation from GRBs across the electromagnetic spectrum is now possible. Rapid and accurate localisation of GRBs and immediate dissemination of their positions drive ground-based follow-up campaigns with small and large optical and infrared robotic telescopes such as ROTSE, REM, the Liverpool and Faulkes Telescopes, while Swift’s onboard narrow field instruments the X-Ray Telescope (XRT) and the UltraViolet Optical Telescope (UVOT) observe a large fraction of GRBs within minutes of their detection with the Burst Alert Telescope (BAT). Following the detection of a very bright optical counterpart to GRB 990123 by ROTSE (Akerlof et al. 1999) contemporaneous with the γ-ray emission and consistent with models predicting optical emission from a reverse shock in the relativistic outflow (Sari & Piran 1999;
–3– Kobayashi & Sari 2000; Kobayashi 2000; Nakar & Piran 2005), there was strong anticipation of routine detection of optical flashes at early times with UVOT, which can detect optical emission below the detection threshold of ROTSE. However, despite rapid and accurate localisation of γ-ray emission and subsequent follow-up by UVOT, bright optical flashes remain elusive (Roming et al. 2006). In some cases, the high redshift of the burst renders it undetectable due to the lack of UVOT filters redder than V-band (i.e. z>4); in others, intrinsic faintness coupled with local extinction pushes the OT below the detection threshold of the UVOT (Oates et al. 2006). Ground-based observations from large robotic telescopes, optimised for the rapid follow-up of GRB alerts (Guidorzi et al. 2005; Guidorzi et al. 2006a; Dai et al. 2006), are increasing the number of well-sampled, early-time multicolor optical light curves over a wide range of magnitudes (Gomboc et al. 2006; Melandri et al. 2006) and revealing more complex decays than expected from pre-Swift light curves (Guidorzi et al. 2005; Stanek et al. 2006; Dai et al. 2006; Monfardini et al. 2006; Oates et al. 2006). Small departures from a power-law decay (i.e. ∆m.0.5 mag) are interpreted as interactions between the expanding fireball and dense clumps in the circumburst medium (Lazzati et al. 2002; Heyl & Perna 2003; Nakar et al. 2003; Guidorzi et al. 2005) while major rebrightenings suggest late-time energy injections (Kumar & Piran 2000; Sari & M´esz´aros 2000; Zhang & M´esz´aros 2002; Stanek et al. 2006; Wozniak et al. 2006; Monfardini et al. 2006); bright optical flares consistent with reverse shocks, however, are still rare or unconfirmed (Bo¨er et al. 2006; Jel´ınek et al. 2006). Complex light curves are now common at X-ray energies with Swift light curves revealing a canonical shape comprising a “3-segment” structure inconsistent with a single power-law decay, with X-ray flares superimposed in at least 50% of bursts. Such complexity has been a major surprise of the Swift era and has led to a general acceptance of long-lived activity (Burrows et al. 2005; O’Brien et al. 2006), challenging standard central engine models. A significant minority of GRB X-ray light curves show a fairly smooth decay consistent with an afterglow at early time (O’Brien et al. 2006) but most of these cases lack simultaneous optical data sufficient to fully test current models. Here we present multicolor BVRi′ light curves of the very bright optical counterpart to GRB 061007 obtained automatically by the 2-m robotic telescope, the Faulkes Telescope South (FTS) at Siding Spring, Australia which began observing 137 s after the onset of the γ-ray emission when the optical counterpart was already decaying from R∼10.3 mag, and continued observing for the next 5.5 hours1 . We show that GRB 061007 is unusual for a Swift burst in that it shows a remarkably simple power-law decay from γ rays to optical wavelengths. Its isotropic equivalent radiated energy is one of the highest ever measured (Eiso ∼ 1054 erg) and its peak νFν energy, Epeak , is well-determined. It therefore provides a useful test 1
Observations obtained as part of the RoboNet-1.0 project: http://www.astro.livjm.ac.uk/RoboNet/
–4– of current spectral energy correlations without the complication of highly-structured light curves (Amati et al. 2002; Ghirlanda et al. 2004; Liang & Zhang 2005; Firmani et al. 2006; Amati 2006). The multiwavelength observations and analyses are presented in Section 2; results and discussion are presented in Section 3, where a shock model is proposed to explain the observed multiwavelength evolution and broad-band spectral energy distribution. We end by considering implications for cosmological spectral energy correlations and present our conclusions in Section 4. Throughout, 1-σ errors are given unless otherwise stated.
2.
Observations and Analysis
On 2006 Oct 07, 10:08:08 UT2 Swift-BAT detected GRB 061007 (trigger 232683) at coordinates R.A.(J2000)=03:05:16, Dec(J2000)=−50:29:15 (3 arc-min 90% containment radius). The prompt emission from GRB 061007 was also detected at MeV energies by Konus/Wind (20 keV−10 MeV) (Golenetskii et al. 2006), Suzaku/WAM (50 keV−5 Mev) (Yamaoka et al. 2006) and RHESSI (20 keV−10 MeV) (Wigger et al. 2006), whose detections provide +29 estimates of the peak energy Epeak of 399+19 −18 keV (Golenetskii et al. 2006), 561−27 keV (Yamaoka et al. 2006) and 391+58 −52 keV (Wigger et al. 2006) respectively. The Swift-BAT γ-ray light curve shows three significant flat-topped peaks with substantial sub-structure, and a small fourth peak (from t0+75s) that shows long exponential decay and faint emission detectable to t∼ 900 s (Markwardt et al. 2006; Schady et al. 2006c). The fluence, Sγ =2.49×10−4 erg cm−2 , measured by Konus-Wind, corresponds to an isotropic equivalent energy Eiso ∼1054 erg (Golenetskii et al. 2006) for the spectroscopic redshift z=1.261 (Osip et al. 2006; Jakobsson et al. 2006). The Swift XRT and UVOT began observing 80 s after the BAT trigger time and identified a very bright X-ray/optical counterpart with a rapid decay rate αx =1.6±0.1 (F ∝ t−α ) (Schady et al. 2006a; Vetere et al. 2006; Schady et al. 2006b). Robotically-triggered photometric observations with ground-based telescopes ROTSE and the FTS began at 26 and 137 s respectively, each identifying the optical counterpart with brightnesses 13.6 mag (unfiltered) (Rykoff et al. 2006) and 10.3 mag (R-band) (Mundell et al. 2006; Bersier et al. 2006) respectively at R.A.(J2000)=03:05:19.56, Dec(J2000)=−50:30:2.6 (±1.′′ 2). After taking the initial 3 × 10 s exposure R-band images, the FTS continued to observe for the next 5.5 hours, cycling through BRi′ filters with V imaging interspersed and using increasing exposure times, until finishing with a 300-s V-band image at UT=15:39:25.322. Initial photometric calibration was performed with GAIA (Optimal PSF Photometry Tool) 2
t=0 throughout the paper.
–5– relative to all the catalogued stars (USNO B1 and NOMAD) in the 4.′ 6×4.′ 6 field of view and the bright star near the OT (RA=03:05:23.5, Dec=-50:30:16.6) was used for PSF determination. A 5 × 120 s R-band image was also acquired using FORS1 on the VLT at 2006 Oct 8.04702 UT (0.6247 days after the burst); the observing conditions were photometric and calibration was performed relative to Landolt standard stars. The FTS R-band were adjusted to the VLT photometry to provide a well-calibrated R-band light curve. Later-time R-band images were acquired with the IMACS instrument on the Magellan-I Baade telescope at Las Campanas Observatory on Oct 9 and 10 and calibrated with respect to the VLT image. Finally, magnitudes are corrected for Galactic extinction: from EB−V = 0.020 (Schlegel et al. 1998), and AV = RV · EB−V = 0.06 mag, with RV = 3.1, we evaluated the extinction in the other filters following Cardelli et al. (1989): AB = 0.08 mag, AR = 0.05 mag and Ai′ = 0.04 mag. Fluxes were derived following Fukugita et al. (1996) and calibration uncertainties for the remaining FTS filters, which do not affect the overall light curve properties, were taken into account when deriving the broad-band spectral energy distribution. Results are summarized in Table 1 and plotted in Figure 1 along with the mask-tagged background subtracted BAT light curve, which was produced by the standard BAT pipeline (HEADAS v 6.1.1). The BAT light curve was masked with the ground-refined coordinates provided by the BAT team (Markwardt et al. 2006). BAT spectra have been corrected for systematics with the “batphasyserr” tool and the original 80 energy channels have been grouped to have 3-σ significant grouped channels. We processed XRT data from 200 s to 400 s after the trigger adopting standard screening and we extracted the XRT 0.3–10 keV spectrum: data were in Window Timing (WT) mode and the average count rate was below 100 counts s−1 so that no correction for pileup was required (Romano et al. 2006). The brightness of the burst and the similarity of the decay rate in all bands (see Figure 1 and Table 2) allows the construction of an early-time broad-band spectral energy distribution (SED), shown in Figure 2. The joint XRT/BAT 0.3–150 keV spectrum from 200 s to 400 s is well fit with a power law with significant absorption in excess of the Galactic, NH (Gal) = 1.75 × 1020 cm−2 (Kalberla et al. 2005), such that βX−γ = 1.01 ± 0.03 (F ∝ ν −β ) and restframe NH = (5.8±0.4)×1021 cm−2 (90% confidence level). As shown in Figure 2 and Table 3, when β is free to vary, the best fit is obtained for βopt−X−γ = 1.02 ± 0.05 and AV =0.48±0.19 mag with an SMC-like extinction profile (Pei 1992).
3.
Results and Discussion
Figure 1 shows the FTS BVRi′ and background-subtracted Swift-BAT light curves of GRB 061007; a single power-law decay with α = 1.72 ± 0.10 (F ∝ t−α ) provides a good
–6– fit to all optical light curves for t < 40 min, resulting in small residuals in this time range and emphasizing the presence of an additional broad, multicolor bump over the time range 40 < t < 350 min (see Table 2). Late-time VLT R-band photometry confirms that the underlying power-law decay continues, showing no evidence for a break before 1000 min (0.624 d). This is consistent with the continuation of the simple power-law decay in the XRT light curve (Schady et al. 2006c) . Similarly, the marginal Magellan detection extends the continuation of the unbroken power-law decay to t = 2495 min (1.73 d), and, combined with the final upper limit R> 23.3 mag at 4087 min (2.84 d - Table 1), constrains the brightness of the host galaxy. A similar power law (α = 1.7 ± 0.1) also provides a good fit to the decay of the long, soft tail of the γ-ray emission for t > 120 s, although slight modulation of the power law due to the intrinsic variability is evident, as illustrated in the expanded view of the light curves of the first 10 min (Figure 1 right panel; Table 2). A similar deviation from the R-band power law may be present at ∼2.3 min (see residuals Figure 1), coincident in time with a small γ-ray flare, but the effect is marginal. The broad-band spectral energy distribution shown in Figure 2 confirms the overall common power-law decay in all bands from γ to optical and is well fitted by an absorbed power law with optical extinction AV ∼ 0.48 ± 0.19 mag with an SMC profile. Comparing this with the NH derived from the X-ray spectrum, we find a deficit of optical extinction, as already found in most GRBs (Stratta et al. 2004; Kann et al. 2006; Galama & Wijers 2001). In principle this could be explained in a number of ways (see e.g. the case of GRB 051111, Guidorzi et al. 2006b): 1) an extinction profile different from those typical of MW, SMC, LMC (e.g. Stratta et al. 2004; Savaglio & Fall 2004); 2) a population of dust grains skewed towards big sizes, possibly as a result of dust destruction due to the GRB itself (e.g. Maiolino et al. 2001; Chen et al. 2006); 3) a significant presence of molecular gas (1999); 4) overabundance of some alpha metals responsible for absorption in X-rays (e.g., GRB 050401, Watson et al. 2006). Below, we explain the observed multiwavelength evolution of GRB 061007 at t > 70 s with an expanding fireball whose optical, X-ray and γ-ray emission is dominated by synchrotron emission from a forward shock whose cooling frequency lies above the X-ray band and typical synchrotron frequency lies just below the optical band; in contrast, the typical frequency of the reverse shock lies in the radio band at early time, with the typical frequency of the forward shock entering the radio band at ∼2.5 days after the burst.
–7– 3.1.
Forward Shock Emission from a Decelerating Fireball
The optical and X-ray afterglow light curves are described by the same power law α ∼ 1.7 to more than 105 sec after the burst; this indicates that these photons are radiated from the same forward shock, and optical and X-ray bands are in the same domain of a synchrotron spectrum3 . In principle a bump of a forward shock (blast wave) light curve, as seen for t> 40 min (Figure 1), can be produced by inhomogeneity of the ambient medium or by energy injection into a blast wave. However, energy injection causes a transition from a blast wave to another with a larger energy, and it leads to a shift of the afterglow decay baseline after the bump (at this late time the evolution of a blast wave is adiabatic). Since the optical flux comes back to the pre-bump decay line, inhomogeneity of ambient medium is favored to explain the bump around 100 min in the optical. The X-ray band (and the optical band) should be below the synchrotron cooling frequency νc , otherwise the inhomogeneity in the ambient medium does not produce an afterglow bump (Kumar 2000). For the synchrotron shock model, the decay index and spectral index at frequencies νm < ν < νc are α = 3(p − 1)/4 and β = (p − 1)/2, respectively (Sari, Piran and Narayan 1998). With p ∼ 3 these give reasonable fits to the observed indices. We note that if a burst occurred in a wind medium with density n ∝ R−s , the decay index in this spectral domain νm < ν < νc is given by α = [2s + 3(p − 1)(4 − s)]/(16 − 4s) ∼ 1.7 for p = 3 and s = 1 (Monfardini et al. 2006). The additional parameter s allows a better fits to the decay index, but still a rather large value of p ∼ 3 is required to explain the observed β. Since the optical and X-ray light curves show no break during the observations, the typical synchrotron frequency νm ∝ t−3/2 should be already below the optical band νm < νopt ∼ 5 × 1014 Hz before our earliest optical observation was made at t = 137 sec, while the cooling frequency νc ∝ t−1/2 should be still above X-ray band νX < νc ∼ 1018 Hz at t ∼ 105 sec. Using the results in Sari, Piran and Narayan (1998), we get 1/2
1/2
νm (t) ∼ 1.0 × 1021 ζ 1/2 ǫB ǫ2e E54 t−3/2 Hz, m νc (t) ∼ 6.8 ×
−3/2 −1/2 1012 ζ −1/2 ǫB E54 n−1 t−1/2 m
(1) Hz
(2)
where ζ = (1 + z)/2.26, E54 = E/1054 ergs, tm is the observer’s time in units of min, n is the density in cm−3 , νm is proportional to (p − 2)/(p − 1) and we have assumed p = 3. The requirements for the break frequencies at t = 137 s and at t = 105 s give constraints on the 3
Even if we assume reverse shock emission, the observed bands should be in the same spectral domain νm < νopt , νX < νc , because reverse shock emission above the cooling frequency vanishes after the shock crossing. In this case we obtain the theoretical value α − 3β/2 ∼ 1 which is not consistent with observation α − 3β/2 ∼ 0
–8– microphysics parameters, ǫB < 3.0 × 10−5 n−2/3 (ζE54 )−1/3 −1/4 ǫB −2 ǫe < 1.7 × 10 (ζE54 )−1/4 3.0 × 10−5
(3) (4)
These energy partition values are somewhat small but are not inconsistent with values derived for other bursts (e.g. Panaitescu & Kumar 2002).
Onset of the Afterglow and Lack of Optical Flash ROTSE detected a 13.6-mag optical counterpart at 26 s after the GRB trigger (Rykoff et al. 2006). This magnitude is dimmer than the extrapolation of our observations. It indicates that the afterglow peaked between 26 s and 137 s, and very likely around the end of the main pulses of the prompt emission tp ∼ 100 s. Using this peak time, we can estimate the Lorentz factor of the fireball at that time (Sari and Piran 1999). It is Γ∼
3(1 + z)3 E 32πnmp c5 t3p
1/8
tp ∼ 230 102 sec
−3/8
−1/8 ζ 3E54 /n
(5)
It is interesting that the tail of the γ-ray emission is described by the same power-law decay α ∼ 1.7 as the optical and X-ray emission (Figure 1); indeed the flux at the dip in the prompt phase at ∼40 s is consistent with the underlying power-law line. If the afterglow started at ∼ 30 s, the Lorentz factor is ∼ 360. No jet break within 105 s 4 gives a lower limit −1/8 of the GRB jet opening angle θ0 (rad) > 0.07 n1/8 E54 ζ −3/8 (t/105 sec)3/8 . There is no reverse shock emission component even though the optical observations started right after the prompt emission phase. We can explain the lack of optical flash naturally in the standard model if the typical frequency of the forward shock emission is lower than optical band νm < νopt at the onset of afterglow (the peak time tp ). At the peak time (the shock crossing time), the forward and reverse shocked regions have the same Lorentz factor and internal energy density, the cooling frequency of the reverse shock is equal to that of the forward shock νc,r (tp ) ∼ νc (tp ). The matter density in the reverse shocked region is much higher than in the forward shock region, and it makes the electron 4
Another possible scenario is that a jet break happened immediately after the prompt phase at t < 100 sec. This scenario requires a very small electron index to explain the decay indices of the optical and X-ray light curves.
–9– temperature lower. The typical frequency of the reverse shock is much lower than that of the forward shock νm,r (tp ) ≪ νm (tp ) < νopt (e.g. Kobayashi & Zhang 2003). The shock regions have the same amount of shock energies at tp , and their luminosities (∼ νFν at νc ) are comparable at tp . Considering that the two emission components have the same spectral index, at tp they contribute equally to the flux at any observed frequency between νm,f and νc . After the shock crossing t > tp the reverse shock emission decays more rapidly as ∼ t−2 , so the contribution from the reverse shock should be masked by the forward shock component in the decay phase. Forward shock emission peaks when a fireball decelerates or when the typical frequency passes through the observed band. If an optical peak is not associated with a spectral change from ν 1/3 to ν −(p−1)/2 , it could be because of the deceleration of a fireball, and the typical frequency could be located below the observed band at that time. In such a case, the optical flash should always be absent. In recent years, ground-based robotic telescopes reported the lack of optical flashes, except for a few cases. The lack might be explained with low νm originating from rather small microphysics parameters, ǫe and ǫB . GRB 060117 and GRB 061007 might be clear example cases.
3.2.
Implications for Late-Time Radio Emission
After the deceleration, the emission frequency of each electron in the shocked ejecta (reverse shock region) drops quickly with time as νe ∝ t−73/48 . Both the typical frequency and cooling frequency drop with this scaling, because after a reverse shock has crossed the ejecta (deceleration), no new electrons are injected and all electrons cool by adiabatic expansion only. The peak power decays as Fν,max ∝ t−47/48 . The flux below or above the typical frequency νm,r evolves as (Sari & Piran 1999; Kobayashi & Sari 2000, their Eq. 3), Fν ∝
t−17/36 ν < νm,r , −(73p+21)/96 −5/2 t ∝t (p = 3) ν > νm,r .
(6)
The optical afterglow of GRB 061007 was very bright. If the forward shock typical frequency is close to the optical band at the deceleration time νm (tp ) ∼ 5 × 1014 Hz, extrapolating our observational results to the deceleration time tp ∼ 100 s, we obtain the peak flux of the forward shock emission Fν,max ∼ 420 mJy. The typical frequency of the reverse shock emission is lower by a factor of ∼ Γ2 (Kobayashi & Zhang 2002), and it is in the radio band νm,r (tp ) ∼ νm /Γ2 ∼ 9.5(Γ/230)−2 GHz. At the deceleration time, the peak flux of the reverse shock emission is larger by a factor of Γ ∼ 230 than that of the forward shock emission. Since the reverse shock emission in the radio band initially decays as t−17/36 until
– 10 – the typical frequency passage, then as t−5/2 , the reverse shock radio emission at 1 day is about ∼ 8.6 (ν/4.8GHz)−1 (t/1day)−5/2 µJy. The typical frequency of the forward shock emission νm ∝ t−3/2 comes to radio band ν = 4.8 GHz about 2.5 day after the burst. Before the passage, the flux increases as t1/2 , and then decays as t−1.7 . The peak value is expected to be ∼ 420 mJy. At low frequencies and early times, self absorption plays an important role and significantly reduces the afterglow flux. A simple rough estimate of the maximum flux is the emission from a black body with the shock temperature (e.g. Sari & Piran 1999; Kobayashi & Sari 2000). The temperature is given by the random energy of the typical electron kB T ∼ me c2 γm ∼ ǫe mp c2 Γ/2 for p=3. If the observed frequency is above the typical frequency, the electrons radiating into the observed frequency have an energy higher than the typical energy, me c2 γm , by a factor of (ν/νm )1/2 . Following the notation of Sari & Piran (1999), the upper limit of blackbody emission is FνBB = π(R⊥ /d)2 ΓSν
(7)
where d = dL /(1 + z)1/2 and Sν = (2ν 2 /c2 )kB T . Therefore FνBB
2 R⊥ ∼ π(1 + z)mp ν ǫe Γ , dL ν 2 t 1/2 ǫe,−2 (E54 /n)1/2 mJy ∼ 6 4.8GHz 1day 2
2
(8) (9)
where ǫe,−2 = ǫe /10−2 and R⊥ ∼ 4.6Γct is the observed size of the fireball, assuming z = 1.26 and dL = 2.7×1028 cm. Note that other emission estimates use only scalings (normalizations are given by the observations), the blackbody upper limit is more model dependent. The results from ATCA (van der Horst et al. 2006a,b) give even lower limits about 50µJy level. Their formal flux measurements are −25 ± 45µ Jy (1.00-1.24 days after the burst) and −1 ± 43µJy (5.03-5.28 days after) at 4.8 GHz. Although the blackbody limit is possibly less reliable, the radio non-detection might imply a high-density environment n ≫ 1. The peak flux of the forward shock emission is given by (Sari, Piran & Narayan 1998), Fν,max (z = 1.26) ∼ 20
ǫB 5 × 10−5
1/2
E54 n1/2 mJy
(10)
As can be seen from Equations 3 and 10, the dependence of Fν,max on density is weak (∼ n1/6 ), and thus the bright afterglow cannot be easily explained by a high density environment.
– 11 – 3.3.
GRB 061007 and Cosmological Spectral–Energy Correlations
Swift GRB 061007 is one of the most luminous GRBs ever detected and thus, with its high-quality γ-ray data, well-characterised light curves and well-determined peak energy, it is an ideal object to test the cosmological spectral–energy correlations, namely the Epeak -Eiso (Amati et al. 2002), Epeak -Eγ (Ghirlanda et al. 2004) and Epeak -Eiso-tb (Liang & Zhang 2005) correlations. As described in Section 2 the prompt γ-ray emission at Mev energies was detected by three additional satellites which provide estimates of the peak energy. While the peak energies measured by Konus/Wind and RHESSI are comparable, that measured by Suzaku/WAM is significantly larger; in addition to possible systematics, this may be explained by the fact that this instrument measured the brightest part of the GRB, as suggested by their low duration estimate (T90∼50 s vs 90 s measured by Konus/Wind)5 . We therefore use the spectral information from Konus/Wind and RHESSI, conservatively assuming an Epeak range including both 90% c.l. intervals provided by the two instruments, i.e. Epeak =394±55 keV, corresponding to a cosmological rest-frame peak energy of 890±124 keV. The isotropic–equivalent radiated energy is Eiso =(1.0±0.1)×1054 erg in the 1-10000 keV cosmological rest–frame, assuming H0 = 70 km/s/Mpc, ΩΛ = 0.3 and Ωm = 0.7 and a Band spectral shape (Band et al. 1993) with parameters and fluence provided by Konus/Wind (Golenetskii at al. 2006). Following Amati (2006) the Epeak -Eiso correlation predicts a corresponding Epeak = 907 keV, fully consistent with the measured value for GRB 061007. Thus, GRB 061007 further confirms the validity of the Epeak -Eiso correlation in the very high radiated energy regime and since this correlation should only hold for prompt emission, confirms that GRB 061007 has prompt emission typical of other Swift bursts on the correlation, despite its unusually simple afterglow emission. Finally, the test of the Epeak -Eγ and Epeak -Eiso -tb correlations, which is of particular importance given their possible use for the estimate of cosmological parameters (e.g. Firmani et al. 2005) requires the detection of a late break in the afterglow light curves. As can be seen in Figure 1, our optical light curves show no evidence of a break up to ∼ 150 ks and the XRT light curve continues unbroken to at least 106 s (Schady et al. 2006c), providing a firm lower limit to the jet break time. Combining our estimates of Eiso and the lower limit to the jet opening angle of ∼ 4.7◦ inferred from the optical light curve, we derive a collimationcorrected energy Eγ > 3.4×1051 erg; this limit would predict Epeak > 1274 keV in order to be consistent with the Epeak -Eγ correlation (Nava et al. 2006), a value somewhat higher 5
The lower bound of the energy band of WAM is somewhat higher than that of Konus/Wind and RHESSI, which may result in an overestimate of Epeak (Amati 2006).
– 12 – than the measured value of 890±124 (90% c.l.). By accounting for the logarithmic dispersion of this correlation (Nava et al. 2006), the deviation of the measured value from this lower limit is ∼1.6 σ. Better consistency is found with the Epeak -Eiso -tb correlation, which, when adopting the parameters reported by Liang & Zhang (2005) predicts a lower limit to Epeak of 986 keV. Using the lower limit to the jet break time t ∼ 106 s suggested by the X-ray light curve, we conclude that GRB 061007 becomes a secure outlier to both Epeak -Eγ and Epeak -Eiso-tb correlations, deviating by more than 3 σ. X-ray light curves for other GRBs suggest the existence of a population of objects similar to GRB 061007 that are inconsistent with these correlations (Sato et al. 2006; Willingdale et al. in prep.), although optical observations were not available for these objects. GRB 061007 might therefore transpire to be the norm rather than an anomaly.
4.
Conclusions
We have presented a multiwavelength study of the very bright GRB 061007 based on optical observations that cover the period from 137 s to 3 days after the burst and Swift XRT and BAT data. We conclude that the afterglow commences as early as 30 − 100 s after the onset of the GRB and the optical, X-ray and underlying γ-ray emission, which are described by the same power-law decay α ∼ 1.7 after this time, represent photons that are radiated from the same forward shock, with optical and X-ray bands in the same domain of a synchrotron spectrum. In contrast, the typical frequency of the reverse shock emission is already in the radio domain at early time, explaining the lack of a bright early-time optical flash and late time radio flare. We note the similarity of the light curves of GRB 061007 to that of GRB 060117, which was only observed in a single optical filter (Jel´ınek et al. 2006) and suggest it too may represent the onset of an early afterglow similar to GRB 061007. Finally, we highlight the simplicity of the early-onset afterglow of GRB 061007, which with its well-behaved light curves, typical circumburst environment and well-determined peak γray energy, make it a useful candidate for inclusion in current cosmological correlations. We find that this event is fully consistent within 1 σ with the Epeak -Eiso correlation and also potentially consistent within 1.6 σ with the Epeak -Eγ and Epeak -Eiso -tb correlations using a lower limit to the jet break time tb , and thus to Eγ , derived from our data. We note that since the X-ray afterglow continued without a jet break until t > 106 s, this event is firm outlier to the Epeak -Eγ and Epeak -Eiso-tb correlations, deviating from these correlations by more than ∼ 3 σ. CGM acknowledges financial support from the Royal Society. AXM, MFB, ER, PTO, acknowledge funding from the Particle Physics and Astronomy Research Council (PPARC).
– 13 – RoboNet1.0 and Swift are supported by PPARC. The Faulkes Telescopes, now owned by Las Cumbres Observatory, are operated with support from the Dill Faulkes Educational Trust. This work is partly based on data taken with ESO telescopes at VLT under program 078.D0752. We thank the kind assistance of the observing staff at Paranal. We thank J.E. Rhoads for assistance in acquiring the Magellan data, L. Nava for input on the spectral correlations and the anonymous referee for useful comments and suggestions.
REFERENCES Akerlof, C., et al. 1999, Nature, 398, 400 Amati L. et al., 2002, A&A, 390, 81 Amati, L. 2006, MNRAS, 372, 233 Arabadjis, J.S., & Bregman, J.N. 1999 ApJ, 510, 806 Band, D. et al. 1993 ApJ, 413, 281 Bersier et al. 2006, GCN Circ. 5709 Bo¨er, M., Atteia, J.L., Damerdji, Y., Gendre, B., Klotz, A. & Stratta, G. 2006 ApJ, 623, L71 Burrows, D. N. et al. 2006 Proceedings of the ”The X-ray Universe 2005”, 26-30 September 2005, El Escorial, Madrid, Spain. Organised by European Space Agency (ESA). Ed. by A. Wilson. ESA SP-604, Volume 2, Noordwijk: ESA Publications Division, 2006, p.877 Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989 ApJ, 345, 245 Chen, S.L., Li, A.G., & Wei, D.M. 2006 ApJ, 647, L13 Dai, X., Halpern, J.P., Morgan, N.D., Armstrong, E., Mirabal, N., Haislip, J.B., Reichart, D.E. & Stanek, K.Z. 2006, astro-ph/0609269 Firmani C., Ghisellini G., Avila-Reese V., Ghirlanda G., 2006, MNRAS, 370, 185 Fukugita, M., et al. 1996, AJ, 111 Number 4, 1748 Galama, T.J. & Wijers, A.M.J. 2001 ApJ, 549, L209 Gehrels, N., et al. 2004 ApJ, 611, 1005
– 14 – Ghirlanda G., Ghisellini G., Lazzati D. 2004, ApJ, 616, 331 Golenetskii, S. et al. 2006, GCN Circ. 5722 Gomboc et al. 2006, to appear in ”Swift and GRBs:Unveiling the Relativistic Universe”, Il Nuovo Cimento, in press Guidorzi, C., et al. 2005 ApJ, 630 Issue 2, L121 Guidorzi, C., et al. 2006, PASP, 118 N.840, 288 Guidorzi, C. et al. 2006, A&A, submitted Heyl, J.S. & Perna, R. 2003 ApJ, 586, L13 Jakobsson, P. et al. 2006, GCN Circ. 5716 Jel´ınek, M. et al. 2006, A&A, 454, L119 Kalberla, P.M.W., Burton, W.B., Hartmann, D., Amal, E.M., Bajaja, E., Morras, R., & P¨oppel, W.G.L., 2005, A&A, 440, 775 Kann, D.A., Klose, S. & Zeh, A. 2006 ApJ, 641, 993 Kobayashi, S. & Sari, R. 2000 ApJ, 542, 819 Kobayashi, S. 2000 ApJ, 545, 807 Kobayashi, S. & Zhang, B. 2003 ApJ, 582, L75 Kumar, P. 2000 ApJ, 538, L125 Kumar, P. & Piran, T. 2000 ApJ, 532, 286 Lazzati, D., Rossi, E., Covino, S., Ghisellini, G., & Malesani, D. 2002, A&A, 396, L5 Liang E. & Zhang B., 2005, ApJ, 633, 611 Maiolino, R., Marconi, A., & Oliva, E. 2001, A&A, 365, 37 Markwardt, C. 2006, GCN Circ. 5713 Melandri, A. et al. 2006, to appear in ”The Multicolored Landscape of Compact Objects, AIP Conf. Proc, in preparation Monfardini, A. et al. 2006, ApJ, 648, 1125
– 15 – Mundell, C. et al. 2006, GCN Circ. 5708 Nakar, E., Piran, T., Granot, J. 2003, New Astronomy, 8, 495 Nakar, E. & Piran, T. 2005, ApJ, 619, L147 Nava, L., Ghisellini, G., Ghirlanda, G., Tavecchio, F., & Firmani, C. 2006, A&A, 450, 471 Oates, S., et al. 2006, MNRAS, 372, 327 O’Brien, P. et al. 2006, ApJ, 647, 1213 Osip, D. et al. 2006, GCN Circ. 5715 Panaitescu, A. and Kumar, P. 2002 ApJ, 571, 779 Pei, Y.C. 1992 ApJ, 395, 130 Romano, P. et al., 2006, A&A, 456, 917 Roming, P. et al., 2006 ApJ, in press (astro-ph/0509273) Rykoff et al. 2006, GCN Circ. 5706 Sari, R. & Piran, T. 1999, A&AS, 138, 537 Sari, R. & M´esz´aros, P. 2000 ApJ, 535,L33 Sari, R., Piran, T. & Narayan,R. 1998 ApJ, 497, L17 Sato, G. et al. 2006 ApJ, accepted Savaglio, S. & Fall, S.M. 2004 ApJ, 614, 293 Schady et al. 2006a, GCN Circ. 5707 Schady et al. 2006b, GCN Circ. 5719 Schady et al. 2006c, MNRAS, submitted (astro-ph/0611081) Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998 ApJ, 500, 525 Stanek, K. Z., Dai, X. et al. 2006, astro-ph/0602495 Stratta, G., Fiore, F., Antonelli, L.A., Piro, L. & de Pasquale, M. 2004 ApJ, 608, 846 van der Horst, A.J., Rol, E. et al. 2006a, GCN Circ. 5720
– 16 – van der Horst, A.J., Rol, E. et al. 2006b, GCN Circ. 5726 Vetere, L. et al. 2006, GCN Circ. 5712 Watson, D. et al. 2006, ApJ, in press (astro-ph/0510368) Wigger, C. et al. 2006, GCN Circ. 5725 Wozniak, P. R., Vestrand, J.A., Wren, J.A., White, R.R., Evans, S.M., Casperson, D. 2006 ApJ, 642, L99 Yamaoka, K. et al. 2006, GCN Circ. 5724 Zhang, B. & M´esz´aros, P. 2002 ApJ, 566, 712
This preprint was prepared with the AAS LATEX macros v5.2.
– 17 –
Table 1. FTS, VLT and Magellan optical observations of GRB 061007. Telescope
∆t (min)
Filter
Exposure Time (s)
FTS
2.28 2.63 3.00 4.63 5.52 6.38 7.50 8.80 12.25 16.73 18.42 20.18 21.93 24.63 27.38 30.15 33.78 37.55 41.45 44.08 46.83 49.63 53.28 57.07 60.85 68.40 70.50 71.23 72.13 73.17
RC RC RC B V SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC B V SDSS-i′ B
10.0 10.0 10.0 10.0 10.0 10.0 30.0 30.0 30.0 60.0 60.0 60.0 120.0 120.0 120.0 180.0 180.0 180.0 120.0 120.0 120.0 180.0 180.0 180.0 240.0 30.0 10.0 10.0 10.0 30.0
Mag ± Err†
10.34 10.68 10.96 12.81 12.52 11.67 13.56 12.85 13.03 15.18 14.28 13.90 15.64 14.78 14.47 16.24 15.36 15.04 16.70 15.71 15.31 16.92 15.98 15.59 17.27 16.25 17.48 17.16 15.94 17.66
± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±
0.11 0.12 0.11 0.20 0.20 0.04 0.10 0.10 0.04 0.09 0.10 0.11 0.07 0.08 0.10 0.09 0.08 0.09 0.09 0.08 0.10 0.10 0.07 0.09 0.11 0.09 0.11 0.21 0.11 0.13
– 18 –
Table 1—Continued Telescope
∆t (min)
Filter
Exposure Time (s)
74.40 75.62 76.90 78.53 80.32 82.08 84.77 87.52 90.30 93.95 97.73 101.55 104.20 106.93 109.68 113.43 117.20 120.98 125.62 145.55 150.92 160.12 171.98 191.77 197.15 203.75 209.13 215.67 242.75 249.33
RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC SDSS-i′ B RC RC RC V B SDSS-i′ SDSS-i′ RC RC V B SDSS-i′
30.0 30.0 60.0 60.0 60.0 120.0 120.0 120.0 180.0 180.0 180.0 120.0 120.0 120.0 180.0 180.0 180.0 240.0 240.0 300.0 300.0 450.0 450.0 300.0 300.0 300.0 150.0 600.0 300.0 300.0
Mag ± Err†
16.48 16.06 17.66 16.59 16.15 17.78 16.72 16.30 17.92 16.95 16.51 18.18 17.03 16.63 18.36 17.26 16.87 18.67 17.50 17.55 17.65 18.38 19.25 17.57 17.63 18.19 18.20 18.90 19.62 18.17
± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ± ±
0.09 0.10 0.20 0.09 0.10 0.18 0.09 0.09 0.16 0.09 0.10 0.17 0.08 0.10 0.14 0.08 0.11 0.17 0.15 0.11 0.10 0.18 0.15 0.13 0.12 0.11 0.11 0.17 0.20 0.13
– 19 –
Table 1—Continued Telescope
VLT Magellan
†
∆t (min)
Filter
Exposure Time (s)
Mag ± Err†
254.72 261.28 273.20 295.43 307.42 319.37 331.28 899.58 2494.80 4086.64
SDSS-i′ RC V B SDSS-i′ RC V RC RC RC
300.0 600.0 600.0 600.0 600.0 300.0 300.0 600.0 240.0 240.0
18.17 ± 0.12 18.95 ± 0.13 19.38 ± 0.20 20.09 ± 0.26 18.39 ± 0.15 19.20 ± 0.25 19.99 ± 0.30 21.48 ± 0.03 23.65 ± 0.50 >23.30
Quoted errors for BVi′ filters are statistical. Absolute calibration requires inclusion of the systematic error 0.25 mag due to the uncertainty in the USNO B1 photometric calibration. This error has been fully included when deriving the SED (See Figure 2 and Table 3).
– 20 –
Table 2: Results of power-law fits to optical and γ-ray light curves as shown in Figure 1. Band Time Range For Fit Temporal Power-Law Decay (α) B t < 40 min 1.72±0.05 V − Fixed to αR R t < 40 min 1.72±0.01 ′ i t < 40 min 1.75±0.06 BAT(15 − 350 keV) t > 40 s 1.7±0.1
Table 3: Best-fitting parameters derived from optical-X-γ spectral energy distribution. SMC‡ LMC MW SMC, βXγ fixed β 1.02±0.05 1.03±0.05 1.03±0.05 1.01 AV (mag) 0.48±0.19 0.54±0.20 0.54±0.20 0.29±0.06 2 χ /d.o.f 195/174 195/174 195/174 197/175
‡
See Figure 2
– 21 –
Fig. 1.— Left: FTS BVRi′ optical and Swift γ-ray light curves of GRB 061007. Optical light curves are fitted with a single power-law fits for t< 40 min and the background-subtracted γ-ray light curve is fitted with a similar power law to data at t > 120 s. The BAT count rate (counts per sec per fully illuminated detector for an equivalent on-axis source) and optical magnitudes are given on the left and right vertical axes respectively. The lower panel shows the optical residuals after removal of the corresponding power law (see Table 2), highlighting the presence of a broad bump in all filters for t > 40 min. Late time R-band photometry from VLT and Magellan images are included at t∼1000 and t∼2500 min respectively and are consistent with extrapolation of the power-law decay. Right: An expanded view of the optical and γ-ray light curves at early time (1