A&A 456, 917-927 (2006)
DOI: 10.1051/0004-6361:20065071
P. Romano1 - S. Campana1 - G. Chincarini1,2 - J. Cummings3,4 - G. Cusumano5 - S. T. Holland3,6 - V. Mangano5 - T. Mineo5 - K. L. Page7 - V. Pal'shin8 - E. Rol7 - T. Sakamoto3,4 - B. Zhang9 - R. Aptekar8 - S. Barbier3 - S. Barthelmy3 - A. P. Beardmore7 - P. Boyd3 - D. N. Burrows10 - M. Capalbi11 - E. E. Fenimore12 - D. Frederiks8 - N. Gehrels3 - P. Giommi11 - M. R. Goad7 - O. Godet7 - S. Golenetskii8 - D. Guetta13 - J. A. Kennea10 - V. La Parola5 - D. Malesani14 - F. Marshall3 - A. Moretti1 - J. A. Nousek10 - P. T. O'Brien7 - J. P. Osborne7 - M. Perri11 - G. Tagliaferri1
1 - INAF-Osservatorio Astronomico di Brera, via E. Bianchi 46, 23807 Merate (LC), Italy
2 - Università degli Studi di Milano, Bicocca, Piazza delle Scienze 3, 20126, Milano, Italy
3 - NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA
4 - National Research Council, 2101 Constitution Avenue, NW, TJ2114, Washington, DC 20418, USA
5 - INAF-Istituto di Astrofisica Spaziale e Fisica Cosmica Sezione di Palermo,
via U. La Malfa 153, 90146 Palermo, Italy
6 - Universities Space Research Association, 10211 Wincopin Circle, Suite 500, Columbia, MD, 21044-3432, USA
7 - Department of Physics & Astronomy, University of Leicester, LE1 7RH, UK
8 - Ioffe Physico-Technical Institute, 26 Polytekhnicheskaya, St. Petersburg 194021, Russian Federation
9 - Department of Physics, University of Nevada, Las Vegas, NV 89154-4002, USA
10 - Department of Astronomy & Astrophysics, Pennsylvania State University,
525 Davey Lab, University Park, PA 16802, USA
11 - ASI Science Data Center, via G. Galilei, 00044 Frascati (Roma), Italy
12 - Los Alamos National Laboratory, MS B244, NM 87545, USA
13 - INAF-Osservatorio Astronomico di Roma, via di Frascati 33, 00040 Monteporzio Catone, Italy
14 - International School for Advanced Studies (SISSA-ISAS), via Beirut 2-4, 34014 Trieste, Italy
Received 22 February 2006 / Accepted 13 June 2006
Abstract
We present observations of GRB 060124, the first event for which
both the prompt and the afterglow emission could be observed simultaneously
and in their entirety by the three Swift instruments.
Indeed, Swift-BAT triggered on a precursor 570 s before the main burst peak,
and this allowed Swift to repoint the narrow field instruments to the burst position
350 s before the main burst occurred.
GRB 060124 also triggered Konus-Wind, which observed the prompt emission in a harder gamma-ray band (up to 2 MeV). Thanks to these exceptional circumstances, the temporal and spectral properties
of the prompt emission can be studied in the optical, X-ray and gamma-ray ranges.
While the X-ray emission (0.2-10 keV) clearly tracks the
gamma-ray burst, the optical component follows a different pattern,
likely indicating a different origin, possibly the onset of external shocks.
The prompt GRB spectrum shows significant spectral evolution,
with both the peak energy and the spectral index varying. As observed in
several long GRBs, significant lags are measured between the hard- and
low-energy components, showing that this behaviour extends over 3 decades in energy.
The GRB peaks are also much broader at soft energies.
This is related to the temporal evolution of the spectrum, and can be
accounted for by assuming that the electron spectral index softened with
time. The burst energy (
1053 erg)
and average peak energy (
keV)
make GRB 060124 consistent with the Amati relation.
The X-ray afterglow is characterized by a decay which presents a break at
s.
Key words: gamma rays: bursts - X-rays: bursts - X-rays: individuals: GRB 060124
The Swift Gamma-Ray Burst Explorer (Gehrels et al. 2004) was successfully launched on 2004 Nov. 20.
Its payload includes one wide-field instrument, the gamma-ray Burst Alert Telescope (BAT; Barthelmy et al. 2005a, 15-350 keV energy band), and two narrow-field instruments (NFIs),
the X-Ray Telescope (XRT; Burrows et al. 2005a, 0.2-10 keV)
and the Ultraviolet/Optical Telescope (UVOT; Roming et al. 2005a, 1700-6500 Å).
The BAT detects the bursts, calculates their positions
to
accuracy and triggers an autonomous
slew of the observatory to point the two narrow-field instruments,
typically within 100 s from the burst onset.
The XRT can provide
positions,
while the UVOT further refines the afterglow localization to
.
GRB 060124 is the first event for which the three Swift instruments
have a clear detection of both the prompt and the afterglow emission.
Indeed, Swift-BAT triggered on a precursor on 2006-01-24 at 15:54:52 UT,
570 s before the main burst peak. This allowed Swift to immediately repoint the NFIs
and acquire a pointing towards the burst
350 s before the burst occurred.
The burst, which had a highly structured burst profile, comprises
three major peaks following the precursor and had one of longest
total durations (even excluding the precursor)
recorded by either BATSE or Swift.
GRB 060124 also triggered Konus-Wind (Aptekar et al. 1995) 559.4 s after the BAT trigger (Golenetskii et al. 2006). The Konus light curve confirmed the presence of both the precursor and the three peaks of prompt emission. The main peak of GRB 060124 was also bright enough to trigger the FREGATE instrument aboard HETE-II (HETE trigger 4012, Lamb et al. 2006) 557.7 s after the BAT trigger.
The prompt emission of GRB 060124 was observed simultaneously by XRT
with exceptional signal-to-noise (S/N) and was detected by UVOT at
V=16.96
0.08 (T+183 s) and V=16.79
0.04 (T+633 s).
This fact makes it an exceptional test case to study prompt emission models, since
this is the very first case that the burst proper could be observed with an X-ray CCD with high spatial resolution imaging down to 0.2 keV. There have been a handful of previous optical detections while the burst prompt emission was still active, due to the increasing
number of robotic followup telescopes. GRB 990123 was observed by ROTSE (Akerlof et al. 1999), and its light curve showed a flare which rapidly decayed, uncorrelated with the
gamma-ray emission. This behaviour was generally interpreted as emission
from the external reverse shock (Mészáros & Rees 1999; Sari & Piran 1999b).
A reverse shock was also likely observed in GRB 050904 (Boër et al. 2006).
A different behaviour was observed for GRB 041219A
(Vestrand et al. 2005; Blake et al. 2005). In this case, the optical emission showed a temporal
variability correlated with the gamma-ray emission, suggesting a common
origin for the two components. The prompt afterglows of GRB 050319
(Wozniak et al. 2005; Quimby et al. 2006) and GRB 050401 (Rykoff et al. 2005) showed a decay matching the late-time one, without the need for extra components. Last,
GRB 050801 (Rykoff et al. 2006) showed a period of
plateaux up to a few minutes after the burst, excluding contribution from
the internal shocks, but possibly arguing for energisation from the
central engine correlated with the X-ray emission.
GRB 060124 was extensively followed from ground-based observatories
with detections starting from one hour after the trigger.
The optical counterpart was identified by the Tautenburg telescope (Kann 2006).
Further observations were carried out with the ARIES 1.04 m (Misra 2006),
the HCT 2 m (Bhatt et al. 2006), the CrAO 2.6 m (Rumyantsev et al. 2006) and
the Asiago 1.82 m (Masetti et al. 2006). Limits were given by ART (Torii 2006),
MASTER (Lipunov et al. 2006), UoM 30 cm (Sonoda et al. 2006),
TAROT (Klotz et al. 2006), and Loiano-1.52 m (Greco et al. 2006).
An absorption redshift of z=2.297
0.001 was determined based on the
CIV
,
1550 doublet
(Mirabal & Halpern 2006; Cenko et al. 2006; Prochaska et al. 2006).
![]() |
Figure 1:
Top: mask-tagged (i.e., background subtracted) BAT light curve of GRB 060124 in the 15-350 keV energy band. The dot-long dash vertical line denotes the Konus trigger time
![]() |
Open with DEXTER |
This paper is organized as follows. In Sect. 2 we describe our observations and data reduction; in Sect. 3 we describe our instrument-based spatial, timing and spectral data analysis; in Sects. 4 and 5 we describe our detailed,
multi-wavelength timing and spectral analysis of the prompt and afterglow emission.
In Sect. 6 we discuss our findings.
Finally, in Sect. 7 we summarise our findings and conclusions.
Throughout this paper the quoted uncertainties are given at 90% confidence level
for one interesting parameter (i.e.,
)
unless otherwise stated.
Times are referred to the BAT trigger T0, t=T-T0, unless otherwise specified.
The decay and spectral indices are parameterized as
,
where
(erg cm-2 s-1 Hz-1) is the
monochromatic flux as a function of time t and frequency
;
we also use
as the photon index,
(ph keV-1 cm-2 s-1).
We adopt a standard cosmology model with H0 = 70 km s-1 Mpc-1,
,
.
At 15:54:52 UT, Swift-BAT triggered and localized GRB 060124
(trigger 178 750, Fenimore et al. 2006; Holland et al. 2006).
The on-board location was distributed 21 s after the trigger, and its position was
RA(J2000
,
Dec(J2000
,
with an
uncertainty of 3
.
The partial coding was 62%.
The spacecraft started to slew to the burst position at t=19 s, and settled at
t=89 s. Table 1 reports the log of the Swift-BAT
observations that were used for this work.
Since BAT was triggered by the precursor, the event-by-event data
which are available from
only include the precursor and do not include the main burst.
The temporal analysis around the main peak was performed using the
mask-tagged 4 channel light curves (15-25, 25-50, 50-100, and
100-350 keV bands) generated by the flight software. For the
spectral analysis, the detector plane histogram (DPH) data
(80-channel spectral data with time-intervals shown in the
bottom panel of Fig. 1) were used.
The BAT mask-weighted light curve (Fig. 1) shows a precursor from t=-3 to t=13 s, with a peak count rate of 1700 counts s-1 (15-350 keV)
at the trigger time, then three major peaks: t= 520-550 s, t=560-580 s
(with the largest flux), and t= 690-710 s. The inset of Fig. 1 shows the 16-s precursor. We note that this is the longest interval (
500 s)
recorded between a precursor and the main GRB event (Lazzati 2005).
![]() |
Figure 2:
Konus light curves of GRB 060124 in the
18-70 keV a), 70-300 keV b), 300-1160 keV c) and
18-1160 keV d) energy ranges at 2.944 s time resolution.
The background levels are marked with dotted horizontal lines.
The vertical lines denote the BAT (T0, dot-short dash) and
Konus (
![]() |
Open with DEXTER |
The most intense part of GRB 060124 triggered Konus-Wind
on 2006-01-24 at
16:04:13.9, 559.4 s after the BAT trigger
(taking into account the
3-s propagation delay from Swift to Wind).
It was detected by the S2 detector which observes the North ecliptic
hemisphere; the incident angle was
.
The count rates are recorded by Konus in three energy ranges:
G1 (18-70 keV), G2 (70-300 keV), and G3 (300-1160 keV).
The early (trigger) part of the time history is recorded from
0.512 s to
229.632 s with time resolution ranging from 2 to 256 ms.
The data before
0.512 s are collected in the waiting mode with
2.944 s time resolution. Time history recorded in the three energy
ranges can be considered as 3-channel spectra.
The full Konus time history of GRB 060124 is reported in Fig. 2,
which shows a weak (4-
and 2.5-
in the 18-70 keV and 70-300 keV
energy bands, respectively) precursor corresponding to the Swift-BAT trigger time,
no statistically significant emission in any Konus energy band
from t=20 to 340 s, resuming emission at t=340 s,
the main double-peaked pulse in the interval t=550-590 s,
and weaker pulses up to t=800 s.
The XRT observations of GRB 060124 started on 2006-01-24 15:56:36 UT, only 104 s after the BAT trigger, and lasted until 2006-02-23 23:32:57 UT. Table 1 reports the log of the Swift-XRT observations that were used for this work,
which include a total net exposure time of 867 s in Windowed Timing (WT) mode and
247 ks in Photon Counting (PC) mode, spread over a
31 d baseline.
The monitoring is organized in 24 sequences. The first (000) was performed as an automated target (AT) and with XRT in auto state. Therefore, in sequence 000 the XRT took an initial 2.5 s image (IM at t=104 s) followed by WT frames (at t=112 s) for the remainder of the first snapshot (continuous pointing at the target).
During most of the first sequence, the count rate was sufficiently high for XRT
to remain in WT mode, the exception being 15 s recorded in PC mode
(for a description of read-out modes, see Hill et al. 2004).
The latter were used to provide a first localization of the afterglow
(Mangano et al. 2006).
The following observations were performed in both PC and WT modes, but the WT data
have a very low S/N and were therefore not used here.
Table 1 reports the log of the Swift-UVOT observations that were used for this work. The first UVOT observation was taken in IMAGE mode, because the burst interrupted a Target-of-Opportunity (ToO) pointing which had required this observing mode. Therefore, the UVOT did not execute the standard AT observing sequence, which uses EVENT mode to provide high time resolution. By contrast, the IMAGE mode data were collected in long integrations. After the first observation, the normal follow-up mode was used. Figure 3 shows the observations for which a detection was obtained (the upper limits are not shown because they do not usefully constrain the behaviour).
![]() |
Figure 3: UVOT light curves. The V and B magnitudes are shown as (red) filled circles and (blue) empty squares, respectively. |
Open with DEXTER |
The BAT data were analysed using the standard BAT analysis
software distributed within FTOOLS v6.0.4. The BAT ground-calculated position, based on data from t=-40 to +19 s, is RA(J2000
,
Dec(J2000
,
with an uncertainty of
(radius, statistical and systematic errors).
Mask-tagged BAT light curves were created in the standard 4 energy bands,
15-25, 25-50, 50-100, 100-350 keV (Fig. 4)
and in the total band 15-350 keV (Fig. 1)
at 3.2-s time resolution. The BAT spectra were extracted over the interval
which includes the precursor (A in Fig. 1,
t=-1.5 to 13.5 s), and in the other 5 time intervals.
Response matrices were generated with the task batdrmgen
using the latest spectral redistribution matrices.
For our spectral fitting (XSPEC v11.3.2) we considered the 14-150 keV energy range and applied an energy-dependent systematic error
vector
.
The Konus data were processed using standard Konus analysis tools.
From
to
s, 64 spectra in 101 channels were
accumulated on timescales varying from 64 ms near the trigger (first four spectra) to 8.2 s afterwards. We also derived the Konus 3-channel spectra and the corresponding response
for several time intervals prior to
,
as well as for the total
burst, for the joint fitting with BAT and XRT.
T90 and T50 (the timescales over which 90% and 50% of the burst fluence is measured) were calculated from the BAT light curves with the task battblocks, using data binned to 1.6 s resolution in the time range t=19.2-1194.4 s (thus excluding the precursor) and are reported in Table 2. These were also calculated for Konus data with Konus-specific tasks. We note, that the Konus T90 and T50 in the G3 energy band are affected by low S/N effects and therefore are probably underestimated.
Although the individual peak durations are
comparable to the duration of other long GRBs, the total duration
(excluding the precursor) is among
the longest recorded by either BATSE (Paciesas et al. 1999) or Swift.
In the 15-150 keV band the fluences of the precursor emission,
the main emission and the following peaks are
10-7,
10-5 and
10-6 erg cm-2, respectively.
![]() |
Figure 4: UVOT (V), XRT, BAT and Konus light curves. The XRT light curves are corrected for pile-up. The count rates have been normalized to the peak of each light curve and offset vertically for clarity. |
Open with DEXTER |
The XRT data were first processed by the Swift Data Center at NASA/GSFC into
Level 1 products (event lists). Then
they were further processed with the XRTDAS (v1.7.1) software package, written by the
ASI Science Data Center (ASDC) and distributed within
FTOOLS to produce the final cleaned event lists.
In particular, we ran the task xrtpipeline (v0.9.9) applying calibration and
standard filtering and screening criteria.
An on-board event threshold of 0.2 keV
was applied to the central pixel of each
event, which has been proven to reduce most of the background due to either
the bright Earth limb or the CCD dark current (which depends on the CCD temperature).
For our analysis we selected XRT grades 0-12 and 0-2 for PC and WT data,
respectively (according to Swift nomenclature; Burrows et al. 2005a).
The X-ray counterpart was detected at the position
RA(J2000
,
Dec(J2000
,
with an estimated uncertainty of
.
This position was determined using the xrtcentroid task (v0.2.7)
on the sequence 002, which is not affected by pile-up, and it takes
into account the correction for the misalignment
between the telescope and the satellite optical axis.
It is
from the refined BAT position (Sect. 3.1).
Table 2: Multi-wavelength T90 and T50.
During sequence 000 the count rate of the burst was high enough
to cause pile-up in the WT mode data. Therefore, to account for this effect, the WT data were
extracted in a rectangular 40
20-pixel region
with a 4
20-pixel region excluded from its centre.
The size of the exclusion region was determined following the
procedure illustrated in the Appendix.
To account for the background, WT events were also extracted within a rectangular box (40
20 pixels) far from background sources. WT events were extracted in three energy bands,
0.2-10 keV (total), 0.2-1 keV (soft, S) and 1-10 keV (hard, H).
We dynamically subtracted their respective backgrounds.
Figure 5 shows the soft and hard
(background-subtracted and corrected for Point-Spread Function,
PSF, losses) light curves, as well as the ratio H/S, with the
BAT trigger as the origin of time. We also created WT light curves in the 1-4 and 4-10 keV bands.
These are shown in Fig. 4.
T90 and T50 were calculated from the XRT/WT light curves and
are reported in Table 2.
During most of sequence 001 the PC mode data were piled-up, as well.
Therefore, we extracted the source events in an annulus with a 35-pixel
outer radius (83
)
and a 3-pixel inner radius.
These values were derived by comparing the observed and nominal PSF
(Vaughan et al. 2006).
Given the decreasing brightness of the source, for the PC data
collected in sequence 002 and the following ones
we extracted the spectra in circular regions of decreasing radius (35 to 10 pixels) to maximise the S/N. PC background data were also extracted in a source-free
circular region (radius 42 pixels).
Ancillary response files were generated with the task xrtmkarf within FTOOLS, and account for different extraction regions and PSF corrections. We used the latest spectral redistribution matrices in the Calibration Database (CALDB 2.3) maintained by HEASARC.
![]() |
Figure 5: X-ray WT mode background-subtracted light curves of GRB 060124. a) Soft band ( S, 0.2-1 keV). b) Hard band ( H, 1-10 keV). c) Ratio of hard to soft count rates. |
Open with DEXTER |
The UVOT located the afterglow of GRB 060124 at RA(J2000
,
Dec(J2000
,
with an internal accuracy of
and an estimated absolute
astrometric uncertainty of
.
This is
from the XRT coordinates (Sect. 3.2). To avoid contamination from a V = 15.4 source located
from GRB 060124, we performed aperture photometry using a circular
aperture with a radius of
centred on the optical afterglow.
A sky annulus of width
and inner radius
was
used since this annulus encloses a large sample of sky pixels, and it avoids
the bright sources near the optical afterglow. Aperture corrections
were performed to convert the
photometry aperture into the
standard aperture radii used to define the photometric zero points.
Six isolated stars were used to compute aperture corrections for each
exposure.
The rms scatter in the aperture corrections is typically 0.02 mag.
The instrumental magnitudes were transformed to Vega magnitudes using the
photometric zero points in the Swift-UVOT Calibration Database (CALDB 2.3). Colour terms have not been applied, but preliminary calibrations suggest that they are negligible.
The adopted photometric zero points were
ZPV= 17.88
0.09,
ZPB = 19.16
0.12,
0.23, and
0.02. Table 3 reports the photometry.
The UVOT fluxes reported in Table 3
are monochromatic fluxes computed at the
central wavelength of the appropriate filter.
Table 3: UVOT multicolour data.
Figure 3 shows the Swift-UVOT V and B light
curves for all observations where a detection of the afterglow was achieved.
The first two V-band points, at 183 and 633 s, differ
from each other by 0.2 mag, which corresponds to a 2.3- variation.
The UVOT light curve shows a remarkable amount of flickering
(also see Sect. 4.2).
We note that since these observations were performed in IMAGE mode,
finer time resolution is not available for these data. Therefore,
the actual, possibly variable behaviour of the optical light curve
during the prompt phase remains unknown.
We can estimate that if the prompt emission in the optical
had occurred only in a time interval comparable to the one of the X-ray
(and with no detectable emission outside, i.e. the effective exposure
time is shorter),
then the UVOT measurement would have been 1.7 mag brighter.
Even with the (ad hoc) estimate for this shorter interval the
optical flux would not be larger by the factor
100 by which the X- and gamma-ray
flux increase in the same amount of time.
In order for the optical emission to
fully track the hard energy band behaviour, a timescale of emission
much shorter (
10 s) than that observed in the X-ray would have been required.
![]() |
Figure 6:
Shift of the emission peak with energy band. The XRT 0.2-1 keV light curve is compared with the BAT 15-100 keV light curve. The lag between the peaks is ![]() |
Open with DEXTER |
The prompt phase of GRB 060124 was observed by UVOT, XRT, BAT and Konus, as shown in Fig. 4, which illustrates the prompt emission in different energy bands. A few observational features can be noticed.
First, the relative importance of the three main peaks (as can be represented by the count rate ratios) varies through the prompt phase, in the sense that the third peak flux increases, relative to the second, as the energy band becomes softer. For example, the third peak in the 0.1-1 keV light curve is actually stronger than the second one, which indicates that strong spectral evolution is taking place (also, see Fig. 5c). This has often been observed before (e.g., Ford et al. 1995).
Second, the peak times of the emission shift with the energy band,
as is shown in detail in Fig. 6,
which shows the 0.2-1 keV XRT light curve compared with the BAT 15-100 keV one.
In this extreme case, the lag between the peaks is 10 s.
The emission starts off simultaneously in all bands (albeit with the shift described above), although it has a considerably lower S/N in the BAT and Konus bands (see Figs. 1 and 2 where the total light curves are shown). Furthermore, the emission becomes less spiky as the energy of the band decreases. This is more quantitatively described by the XRT and BAT T90 and T50 values which are reported in Table 2 and illustrated in Fig. 7. While the T90s are within 12% (fractional), the T50s vary by as much as 42%. This means that the peaks are much broader in the softer bands.
The X-ray peaks can be formally modelled with 12 Gaussians. When modelled with power laws, the prompt emission peaks show very high rise slopes, ranging from 5 to 12 if the time is referred to the BAT trigger, of the order of 1-2 if the initial time is set to the first point of the flare. We note that there is evidence of continued flaring after the end of the first snapshot (see Figs. 5 and 8), and our T90s and T50s may therefore be underestimated.
The afterglow phase was observed by XRT and UVOT. Since the XRT data are characterized by detection over a longer time and with a finer sampling, a joint fit is strongly biased in favour of the X-ray data. Therefore, we fit the XRT and UVOT data separately.
The XRT/PC portion of the light curve (t> 104 s)
was first fit with a simple power-law model using the BAT trigger as origin of time
(we perform a fit with a different origin,
,
below).
The fit yields a decay index
0.02 and
(94 degrees of freedom, d.o.f.).
A fit with a broken power law
for
and
for
,
where
is the time of the break, yields
0.04 and
0.06,
and a break at
105 s
(
,
92 d.o.f.).
This is a significant improvement over the simple power law
(null hypothesis probability 1.9
10-6).
![]() |
Figure 7:
Variation of T90 (black filled circles) and T50
(red empty squares) with energy. The T90 and T50 error bars correspond to the light curve binning (with the exception of the points at ![]() |
Open with DEXTER |
We note that the XRT/PC light curve is affected by flickering,
so that the residuals to the fits are dominated by rapid variations.
Their intensity decreases as the time increases, a behaviour
which has been observed before
and which is due to both the logarithmic binning in our light curves and
to the power-law decrease of the source counts.
We tested the effect of
flickering by removing the points in the light curve which were
at least 2
above the model being tested for.
Although the normalization of both the power-law and the
broken power-law fits decrease, the slopes are consistent with
the previous values within the uncertainties, and
drops down to
1.2.
It is also interesting to fit the light curves considering an initial time
different from the BAT trigger time, i.e. the beginning of the main peak.
To this end, we performed a fit of the PC data,
using a broken power-law model with an initial time fixed to the Konus trigger time
s. We obtain
,
,
and break at
105 s
(
,
92 d.o.f.),
which is consistent with what we found for the fits which adopted
the BAT trigger as the starting time.
The late-time UVOT light curve (
s)
is characterized by
0.15 in the V band
and
0.23 in the B band.
A simultaneous fit of V and B band data with the
constraint of a single slope yields
0.06 (68% error,
,
7 d.o.f.). We infer an observed value of the colour B-V=0.47
0.08 and an optical spectral index
0.3
(the latter was derived using an extinction correction of
mag).
Figure 8 shows the homogenized Swift-UVOT light curve
(the wavelength dependency of the UVOT bands was taken into account)
superimposed (by scaling) on the XRT flux light curve.
![]() |
Figure 8: XRT (filled circles) light curve in units of erg s-1 cm-2 and UVOT (empty squares) flux density (erg s-1 cm-2 Å-1) scaled to match the X-ray. The solid line is the best fit model to the XRT afterglow data. The arrow marks the break position. |
Open with DEXTER |
Table 4: Multi-wavelength spectral fits of the prompt phase of GRB 060124.
Table 5: Multi-wavelength energetics of the prompt phase of GRB 060124.
As a way to characterize the fundamental burst parameters, we
used the full Swift and Konus energy ranges.
To this end we considered the 6 time intervals
described in Sect. 3.1 and extracted spectra
from BAT in interval A, from XRT and BAT data in intervals 1 and 2,
and XRT, BAT and Konus in intervals 3-5. We fit all the spectra extracted in the same time interval simultaneously. We considered an absorbed power law model, with two absorption components,
the Galactic one
1020 cm-2,
and one at the GRB redshift (z=2.297, Cenko et al. 2006).
was free to vary between 7.65
1020and 9.90
1020 cm-2.
The lower limit was derived from the Galactic extinction estimate
along the GRB line of sight
mag (Schlegel et al. 1998), with a conversion into Hydrogen column
cm-2 (Predehl & Schmitt 1995).
The second limit is the upper limit of neutral Hydrogen maps
(Dickey & Lockman 1990), calculated with the nh tool within FTOOLS.
We note that the Leiden/Argentine/Bonn (LAB) Survey of Galactic HI (Kalberla et al. 2005)
provides
1020 cm-2, which is within the considered range. We considered different emission models, ranging from a simple power-law,
,
a cutoff power-law,
and the GRBM model (Band et al. 1993),
for
,
for
,
where E0 is the cutoff energy, and
and
are the low- and high-energy photon indices. The fits were performed in the 0.5-10, 14-150, 20-2000 keV energy ranges for XRT/WT, BAT and Konus, respectively.
Our results for the best fitting model in each time interval
are reported in Table 4
and illustrated in Fig. 9.
At this time resolution, the photon index generally
increases (softens) throughout the prompt phase (Fig. 9c).
The precursor is significantly softer than the main burst prompt emission.
However, we note that in the
X-ray range, where a finer time sampling is available for spectroscopy, the
hardness ratio tracks the count rate on a fine scale (see Fig. 5c).
Therefore, the emission is softer when its intensity is lower
(also, see Sect. 5.2).
There is marginal evidence for variations in the intrinsic column density,
with a 3-
decrease between the third and fourth data point
(Fig. 9b).
We note that the break energy lies outside the XRT band even in the
last bin (Fig. 9d).
For comparison, we considered the Konus spectrum around the Konus trigger (t=567.3-575.5) as representative of the most intense part of the burst. Fitting
this spectrum with a GRBM model, we obtain
,
and
E0=307-110+190 (
,
61 d.o.f.).
We derive
keV.
![]() |
Figure 9:
Multi-wavelength time-resolved spectral analysis.
The models and fit results are reported in
Table 4.
a) BAT 15-350 keV light curve with the 6 time intervals chosen
for the analysis marked;
b)
![]() ![]() ![]() |
Open with DEXTER |
We also fitted the total spectra
(t=301.2-811.2 s, last line in Table 4).
Based on this last fit, we derive an observed and rest-frame peak energy
of
keV and
keV, respectively.
We also calculated integrated quantities,
such as the fluence
,
the peak photon and energy flux
in several energy bands.
These are reported in Table 5.
The derived value of the isotropic energy, for a time integration
of 510 s in the source 1-104 keV energy range, is
1053 erg.
With these values, this burst is consistent with the
-
relation (Amati et al. 2002).
The rest-frame isotropic luminosity in the 30-10 000 keV range is
1053 erg s-1.
This value is also consistent with the
-
relation (Yonetoku et al. 2004).
To study the spectral evolution during the prompt phase in more detail than possible with BAT or Konus data, we created 23 time intervals that covered the prompt phase, and extracted XRT/WT spectra in the regions defined above. We fit them in the 0.5-10 keV band, with an absorbed power law model, with the two absorption components described in Sect. 5.1. The results are shown in Fig. 10. As Fig. 5c also shows, there is strong spectral evolution during the prompt phase. In particular, the photon index has a minimum (i.e. the spectrum becomes harder) at the peaks of emission. Also, we note that the photon index achieves a plateau following each peak in the light curve, with each plateau becoming successively softer. During this time interval, the photon index increases by about 1, reaching a value consistent with the later afterglow spectrum by the end of the last peak.
The only spectroscopic information on the pure afterglow comes from the XRT/PC data.
We performed a fit of the PC data in the 0.5-10 keV band,
which were rebinned with a minimum of 20 counts per energy bin
to allow fitting within XSPEC.
We considered an absorbed power law model, with the two absorption components
described in Sect. 5.1.
The fits for each sequence yield consistent
photon indices and column density at the GRB redshift, as reported in
Table 6. We also fitted all data sets together and obtained a photon index
0.06,
1020 cm-2, and a column density at the GRB redshift
1022 cm-2 (
;
216 d.o.f., see Fig. 11).
We confirm the presence of an excess of
with respect to the
Galactic hydrogen column along the line of sight on the order of 1022 cm-2.
![]() |
Figure 10: Time-resolved spectroscopy of the prompt emission of GRB 060124 with XRT/WT data. a) Normalized XRT count rate. b) Time evolution of the photon index. The vertical dashed lines mark the position of the peaks. |
Open with DEXTER |
![]() |
Figure 11: Spectroscopy of the afterglow emission of GRB 060124 with XRT/PC data. Top: PC data fit with an absorbed power law model. The normalizations are relative to sequence 001, and are 0.29 (002), 0.09 (003), 0.05 (004), and 0.02 (005). Bottom: ratio of the best-fit models to the data. |
Open with DEXTER |
One of the main results of Swift is the identification of a steep-decay phase which is observed in the early X-ray light curves of many GRBs, followed by a shallower phase (O'Brien et al. 2006; Nousek et al. 2006; Tagliaferri et al. 2005). Quite often an erratic flaring behaviour is observed superimposed on these phases (Burrows et al. 2005b; Romano et al. 2006; Falcone et al. 2006). The most accredited explanation for the steep phase is that we are seeing high-latitude emission due to the termination of the central engine activity (Kumar & Panaitescu 2000; Zhang et al. 2006), while flares are probably due to the reactivation of the central engine. Here, for the first time, we have the chance to follow in the soft X-ray and in the optical/UV both the prompt phase and its transition to the afterglow phase.
![]() |
Figure 12: Spectral energy distribution of GRB 060124 during the prompt phase. SEDs (1) and (2) are derived from simultaneous fits to XRT and BAT data fits, while SEDs (3)-(5) from XRT, BAT and Konus fits (see Table 4). |
Open with DEXTER |
Table 6: Spectral fits of the afterglow phase of GRB 060124.
The panchromatic observation of GRB 060124 during the prompt phase
provides the opportunity to study the broadband
emission of the GRB prompt emission in great detail. Time-resolved
spectra were analysed for 5 temporal segments
(Fig. 12), each spanning nearly 4 decades in
frequency (from the XRT to the BAT or Konus band). The
flux levels at even lower frequencies (the UVOT band) are also
measured, but at a much coarser temporal resolution.
Whenever the joint BAT-Konus data are available (bin 3-5), a peak energy
can be computed. The time-resolved spectral fits reveal an evolution of
to lower energies.
The origin of
has been extensively discussed in the literature
(Ryde 2005; Pe'er et al. 2005; Rees & Mészáros 2005; Tavani 1996; Zhang & Mészáros 2002; Mészáros et al. 1994; Pe'er & Waxman 2004; Pilla & Loeb 1998).
The suggestions range from the optically-thin characteristic synchrotron
emission frequency, to the optically-thick Comptonised spectral peak, to
the thermal peak of the fireball photosphere. In particular, there
have been questions about whether the synchrotron model can account
for the observed spectrum (Lloyd & Petrosian 2000; Ghisellini et al. 2000), since the
standard synchrotron model requires fast cooling so that the photon
index below the
is expected to be
1.5. This is
because for internal emission models (e.g. the internal shock model or
the internal magnetic dissipation model), the magnetic fields in the
dissipation region are so high that the electron cooling time scale is
much shorter than the dynamical time scale. The measured
photon index in the bins 2-5 are all close to this value, suggesting
that the fast-cooling synchrotron model is a plausible mechanism at
least for this burst. Lacking a temporally resolved UVOT light curve,
one cannot derive a reliable time-resolved SED extending to the UVOT
band. Nonetheless, connecting the time-averaged UVOT flux level to the
flux level in the XRT band requires a further spectral break between
the XRT and the UVOT band, which may be consistent with a cooling
break or a synchrotron self-absorption break.
Both the spectral peak energy
and the spectral indices
evolve with time. This is likely
due to the complicated interplay between particle heating and cooling
as well as the emission contributions from many emission units
(e.g. internal shocks). One interesting property is that the second
main peak (bin 5,
keV) is significantly softer than the
first main peak (bin 4,
keV). Within the internal
dissipation models, one has
(Zhang & Mészáros 2002),
where L is the wind luminosity of the GRB outflow, r is the emission radius,
is the bulk Lorentz factor, and
is the time variability from the central
engine. The second part of proportionality is valid specifically for
the internal shock model. The 15-150 keV fluence of the second main peak
is about one order of
magnitude lower than that of the first main peak
.
Given the similar durations of the
two episodes, the luminosity of the second epoch is lower by
a factor of
10. This would account for a factor of 3 reduction
of the
.
The data then require a larger emission radius
(by a factor of
3). In the internal shock scenario, this
requires that the bulk Lorentz factor is somewhat higher. This may be
consistent with the idea that the environment of the burst becomes
cleaner with time as the initial jet has loaded most of the baryons
along the path.
This would be the first piece of evidence of an increase of the
Lorentz factor during the prompt phase.
The fact that the pulses in the softer band tend to be broader is related to the temporal evolution of the spectrum. If the electron spectral index does not evolve with time (as is usually assumed in the afterglow modeling), the rising/falling time scales of the light curves should be similar in different energy bands. However, if during the acceleration processes the electron spectrum hardens with time, one can get a wider, soft wing with respect to a narrow, hard spike. The temporal decay slopes after the peaks are shallower than 2. This suggests that the decay is not due to the adiabatic cooling of the emission region. Rather, electron heating may be still going on but with a reduced rate. If one assumes that the electron spectral index softens with time during the decay, one can account for a longer soft tail with respect to the narrow hard spike. This is supported by our finding that the photon indices are systematically smaller near the peaks (Fig. 10).
The super-long duration of the burst suggests that the GRB central engine is active for hundreds of seconds. This has become a general feature of both long and short GRBs in view of the commonly detected X-ray flares hundreds of seconds after the burst trigger, which are generally interpreted as late central engine activities (Burrows et al. 2005b; Zhang et al. 2006; Romano et al. 2006; Falcone et al. 2006; Barthelmy et al. 2005b; Campana et al. 2006; Liang et al. 2006). The second main peak (excluding the precursor) would have been categorized as an X-ray flare had it been even softer. The "precursor'' emission leads the main peak by about 500 s. This is inconsistent with precursor models that involve an erupting fireball from a massive stellar envelope, which predicts a precursor leading the main episode by several seconds (Ramirez-Ruiz et al. 2002). A straightforward possibility is that the precursor is also due to early central engine activity. Suggestions to produce episodic, erratic central engine activity include fragmentation of the collapsing star (King et al. 2005) or the accretion disk (Perna et al. 2006), a magnetic barrier modulated accretion flow (Proga & Zhang 2006), or the magnetic activity of a massive neutron star central engine (Dai et al. 2006). An interesting question for these models would be how to generate a weak precursor 570/(1+z) =173 s before the main episode.
The origin of the optical emission during the
prompt phase is uncertain from the available data.
However, a tracking behaviour would require timescales of optical emission
in the order of 10 s, which is short, compared with GRB 990123 (Akerlof et al. 1999).
Lacking a clear tracking behaviour with the prompt gamma-ray and X-ray light curves
suggests that it may come from a different emission site, e.g. the
reverse shock (Sari & Piran 1999a; Mészáros & Rees 1999,1997).
Nonetheless, the level of the prompt optical emission clearly rules out the
presence of a prominent reverse shock component as discovered in GRB 990123
(Akerlof et al. 1999) and GRB 021211
(Fox et al. 2003; Li et al. 2003), which generally requires that the magnetic fields in the reverse shock
region are significantly stronger than that in the forward shock region
(Fan et al. 2002; Kumar & Panaitescu 2003; Zhang et al. 2003), and that
the fireball has comparable energy in baryonic materials
and magnetic fields (i.e.
,
Zhang & Kobayashi 2005).
The contribution of the reverse shock could be less significant when
is low (baryonic) or very high (Poynting-flux-dominated, Zhang & Kobayashi 2005). In particular, various
physical processes could lead to the suppression of the reverse shock
emission even if the fireball is baryonic
(e.g. Kobayashi et al. 2005; Kobayashi 2000; Nakar & Piran 2004).
As a result, although the insignificant reverse shock emission is consistent with a Poynting
flux dominated GRB model (e.g. Lyutikov & Blandford 2003; Roming et al. 2005b),
the baryonic GRB model is not disfavoured by the data.
The early part of the afterglow (
)
can be explained well in
terms of the standard model (Sari et al. 1998). In this case, the
optical and X-ray data constrain the cooling frequency
between
the optical and X-ray bands. In fact, the X-ray and optical spectral indices
0.10 and
0.3 differ by
0.5, and constrain
1016 Hz. The
temporal behaviour is fully consistent with this interpretation: both the
optical and X-ray decay slopes are consistent with the model prediction:
0.10 (valid for
)
and
0.4 (valid for
). The latter equation is valid only for a homogeneous external medium (the wind model would require a much steeper
). The electron distribution index is
0.15, which is in agreement with the shock theory.
The X-ray decay slope after the break is too shallow for the break being due to
a standard, sideways-expanding jet, in which case
would
be expected (Rhoads 1999). The break magnitude
0.08 may be consistent with a jet which propagates in a wind environment
suffering no sideways expansion (for which
). However, the
wind model is ruled out by the early-time optical data. A solution which
satisfies all the avaliable constraints is a structured jet having a homogeneous core surrounded by power-law wings with energy profile
(Panaitescu 2005). The
break magnitude is dependent upon the index q (and on the location of
). With the observed values, we infer
.
The break
magnitude below the cooling frequency (i.e. in the optical) is slightly
different:
.
This difference, however, is quite small and
may be easily missed.
It is interesting to compare the expected break time according to the
-
-
relation (Ghirlanda et al. 2005; Liang & Zhang 2005), where
is the temporal break identified in optical data.
According to this relation one would expect a break at
105 s, for
1053 erg and
keV,
which is broadly consistent with the observed break time 105 s.
According to the Amati relation (Amati 2006; Amati et al. 2002), the
expected intrinsic spectral break is
keV.
This is generally consistent with the measured value
keV.
We have presented observations of GRB 060124, the first event for which
both the prompt and the afterglow emission could be observed simultaneously
and in their entirety by the three Swift instruments.
Since BAT triggered on a precursor 570 s before the main burst peak,
Swift repointed the narrow field instruments to the burst position
350 s before the main burst occurred.
GRB 060124 also triggered Konus-Wind, which observed the prompt
emission up to
MeV energies.
Thanks to these exceptional circumstances, we could report on the temporal
and spectral properties of the prompt emission in the optical, X-ray and
gamma-ray ranges, as well as on the optical and X-ray properties of the
afterglow.
We find that while the X-ray emission (0.2-10 keV) clearly tracks the
gamma-ray burst, the optical component follows a different pattern,
likely indicating a different origin, possibly the onset of external
shocks.
The prompt GRB spectrum shows significant spectral evolution,
with both the peak energy and the spectral index varying. As observed in
several long GRBs, significant lags are measured between the hard- and
low-energy components, showing that this behaviour extends over 3
decades in energy. The GRB peaks are shown to be much broader at soft energies.
This is related to the temporal evolution of the spectrum, and can be
accounted for by assuming the electron spectral index softened with
time. The burst energy (
1053 erg)
and average rest-frame peak energy
keV)
make GRB 010624 consistent with the Amati relation.
The X-ray afterglow is characterized by a decay presenting a break at
s, possibly due to the crossing of the cooling frequency
.
Acknowledgements
We thank the referee for a swift and through report and D. A. Kann for useful discussions. This work is supported at OAB by ASI grant I/R/039/04, at Penn State by NASA contract NAS5-00136 and at the University of Leicester by PPARC. We gratefully acknowledge the contributions of dozens of members of the XRT and UVOT team at OAB, PSU, UL, GSFC, ASDC, and MSSL and our subcontractors, who helped make this instruments possible. The Konus-Wind experiment is supported by the Russian Space Agency contract and RFBR grant 06-02-16070.
Table 1: Observation log of GRB 060124.
![]() |
Figure A.1: Pile-up-driven spectral variations. XRT/WT photon index as a function of observed source count rate and size of the excluded region at the centre of the PSF (Sect. A). The photon indices are shifted along the ordinate direction for clarity. |
The pile-up effect is a well known feature of CCD detectors (Ballet 1999).
It occurs when the input count rate is so high that
more than one X-ray photon is collected in a given 3
3 pixel
region defining an X-ray event in a single frame and the charges produced by the two separate
events are recorded as one. This effect causes groups of relatively soft
X-ray photons to be detected as one single high-energy photon.
The observed spectrum of an astrophysical source, therefore, is
distorted from its true form, with a general hardening of the slope
(decrease of the photon index).
Another feature of CCD detectors is that the charge produced by a cosmic ray is usually spread over many pixels, as opposed to charge produced by X-ray photons. For this reason, a classification of the possible shapes of the charge distribution is defined associating to each of them a grade number: low grades correspond to confined charge distributions and high grades to extended charge distributions. Cosmic rays are then recognized and rejected on the basis of the event grade number. Pile-up also occurs when two X-rays are collected in neighbouring pixels in one frame. In this case groups of low grade events can be recorded as a single high grade event. Then, when pile-up occurs, the observed grade distribution of the events will show a deficit of low grade events and an excess of high grade events as compared with event grade distribution of fainter sources. Since higher grade events are rejected as cosmic rays, the pile-up effect causes an apparent loss of flux, particularly from the centre of the PSF where the count rate intensity is higher.
The standard way to correct for the pile-up effect is to determine the extension of the PSF core affected by pile-up and exclude from the analysis all the events that fell within that region. Based on the three observable effects we have described, three different methods can be formulated to check for the presence of pile-up and estimate the inner PSF region to exclude. They involve a study of the distortion of 1) the radial intensity profile; 2) the spectrum; 3) the grade distribution. The first method was discussed in detail by Vaughan et al. (2006), with particular attention to Swift-XRT observations in PC mode, and it is widely applied as standard correction technique for PC data. The present analysis will focus on the following two methods as favourite pile-up diagnostic tools for Swift-XRT WT data.
Ground testing and calibration of the XRT integrated system
carried out in 2002 Sep. 23-Oct. 4 at the Panter laboratory
of the Max-Planck-Institut fur extraterrestrische Physik showed
that WT mode data are affected by pile-up for point source
input intensities above 250 counts s-1.
To determine the extent to which pile-up distorts the spectral shape,
we selected 5 time intervals during which the observed source
count rate was
<100, 100-200, 200-300, 300-400, and >400 counts s-1.
For each count rate range, spectra were extracted from
rectangular 40
20-pixel regions
(40 pixels long along the image strip and 20 pixels wide)
with a region of increasing size (0
20-15
20 pixels)
excluded from its centre. The spectra were fit with an absorbed power law with
fixed to the one obtained from the fit to the spectrum
extracted with the exclusion of 15 pixels (which is assumed to be
unaffected by piled-up).
The spectra had 1433-2529 counts, depending on intensity cut and
exclusion region. Figure A.1 illustrates the pile-up-driven spectral
changes. It demonstrates that in the range 0-100 counts s-1 the
source is not piled-up, since no point deviates from the horizontal line,
which represents the best fit to the photon index values obtained in the
10-15 pixel exclusion range. For this intensity range, there is no need
to exclude any pixels from the extraction region.
The source is only affected by moderate pile-up for 100-300 counts s-1;
in this case it is necessary to exclude 1 pixel from
the extraction region. Finally, two and four pixels must be excluded for 300-400 counts s-1)
and >400 counts s-1, respectively.
To determine the extent to which pile-up distorts the grades, we studied the grade distribution for the same 5 time intervals and the same exclusion regions as for the spectral study. We derive analogous conclusions on the number of pixels to exclude from the extraction regions as a function of observed source count rate. The definitions of XRT grades for timing modes are different from the ones for PC mode (Burrows et al. 2005a). Grades in WT span from 0 to 15, and good events are likely associated to grades from 0 to 2 only. If we assume that a good approximation of the expected grade distribution with no pile-up is obtained excluding 10 central pixels, the distributions we obtained show a deficit at grade zero (associated to an increase at higher grades) that becomes negligible if a region of at least 4-pixels is excluded.