A&A 450, 59-68 (2006)
DOI: 10.1051/0004-6361:20054172
P. Romano1 - A. Moretti1 - P. L. Banat1 - D. N. Burrows2 - S. Campana1 - G. Chincarini1 - S. Covino1 - D. Malesani3 - G. Tagliaferri1 - S. Kobayashi2,4 - B. Zhang5 - A. D. Falcone2 - L. Angelini6,7 - S. Barthelmy6 - A. P. Beardmore8 - M. Capalbi9 - G. Cusumano10 - P. Giommi11 - M. R. Goad8 - O. Godet8 - D. Grupe2 - J. E. Hill6,12 - J. A. Kennea2 - V. La Parola10 - V. Mangano10 - P. Mészáros2 - D. C. Morris2 - J. A. Nousek2 - P. T. O'Brien8 - J. P. Osborne8 - A. Parsons6 - M. Perri9 - C. Pagani2,1 - K. L. Page8 - A. A. Wells8 - N. Gehrels6
1 - INAF - Osservatorio Astronomico di Brera, via E. Bianchi 46, 23807 Merate (LC), Italy
2 - Department of Astronomy & Astrophysics, Pennsylvania State University,
525 Davey Lab, University Park, PA 16802, USA
3 - International School for Advanced Studies (SISSA-ISAS), via Beirut 2-4, 34014 Trieste, Italy
4 - Center for Gravitational Wave Physics, Pennsylvania State University,
104 Davey Lab, University Park, PA 16802, USA
5 - Department of Physics, University of Nevada, Las Vegas, NV 89154-4002, USA
6 - NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA
7 - Department of Physics and Astronomy, The Johns Hopkins University,
3400 North Charles Street, Baltimore MD 21218 USA
8 - X-Ray and Observational Astronomy Group, Department of Physics & Astronomy, University of Leicester, LE1 7RH, UK
9 - ASI Science Data Center, via G. Galilei, 00044 Frascati, Roma, Italy
10 - INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica Sezione di Palermo,
via U. La Malfa 153, 90146 Palermo, Italy
11 - Agenzia Spaziale Italiana, Unità Osservazione dell'Universo, Viale Liegi 26, 00198 Roma, Italy
12 - Universities Space Research Association, 10211 Wincopin Circle, Suite 500, Columbia, MD, 21044-3432, USA
Received 7 September 2005 / Accepted 3 January 2006
Abstract
We present observations of XRF 050406, the first burst detected by Swift
showing a flare in its X-ray light curve.
During this flare, which peaks at
s after the BAT trigger,
a flux variation of
in
a very short time
was observed.
Its measured fluence in the 0.2-10 keV band was
erg cm-2,
which corresponds to 1-15% of the prompt fluence.
We present indications of spectral variations during the flare.
We argue that the producing mechanism is late internal shocks,
which implies that the central engine is still active at 210 s,
though with a reduced power with respect to the prompt emission.
The X-ray light curve flattens to a very shallow slope with decay index
of
0.5 after
4400 s, which also supports continued central
engine activity at late times.
This burst is classified as an
X-ray flash, with a relatively low fluence (
10-7 erg cm-2 in
the 15-350 keV band,
erg),
a soft spectrum (photon index 2.65), no significant
flux above
50 keV and a peak energy
keV.
XRF 050406 is one of the first examples of a well-studied X-ray light curve of an XRF.
We show that the main afterglow characteristics are qualitatively similar
to those of normal GRBs.
In particular, X-ray flares superimposed on a power-law light curve have now been
seen in both XRFs and GRBs. This indicates that a similar mechanism may be at work
for both kinds of events.
Key words: gamma rays: bursts - X-rays: bursts - X-rays: individuals: XRF 050406
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 (15-350 keV)
Burst Alert Telescope (BAT; Barthelmy et al. 2005a),
and two narrow-field instruments,
the X-Ray Telescope (XRT; Burrows et al. 2005a)
and the Ultraviolet/Optical Telescope (UVOT; Roming et al. 2005).
BAT detects the bursts, calculates their position
to
accuracy and triggers an autonomous
slew of the observatory to point the two narrow-field instruments.
The XRT, which operates in the 0.2-10 keV energy range,
can provide
positions,
while the UVOT, which operates in the 1700-6000 Å wavelength range,
can further refine the afterglow localization to
.
With its unique fast re-pointing capabilities Swift set out to investigate
the very early phases of gamma-ray burst (GRB) afterglows,
beginning as early as one minute after the BAT trigger.
During the initial activation and calibration phases, which ended on 2005 Apr. 5,
Swift discovered 25 GRBs.
The narrow-field instruments were re-pointed towards seven of them within a
few hundred seconds, and such is the case for GRB 050406.
On 2005 Apr. 6 at 15:58:48.40 UT, BAT triggered on GRB 050406
(trigger 113872; Parsons et al. 2005), and located it at
RA(J2000
,
Dec(J2000
,
with an uncertainty of 3 arcmin (95% containment; Krimm et al. 2005).
The derived value for the time during which 90% of the burst fluence is observed
was
s in the 15-350 keV band.
In the 15-25 keV band the light curve peak had a fast-rise,
exponential decay (FRED) profile, while in the 25-50 keV band, the shape was more
symmetric, with the peak starting
2 s earlier (Krimm et al. 2005).
Both the peak and time-averaged spectra were well fit by a simple power-law
with a time-averaged spectrum photon index of
(90% confidence; Krimm et al. 2005).
The fluence in the 15-350 keV band was
erg cm-2.
The gamma-ray characteristics of this burst, i.e.
the softness of the observed spectrum and
the absence of significant emission above
50 keV,
classify GRB 050406 as an X-ray flash (XRF; Heise et al. 2001).
From now on, we shall therefore refer to this event as XRF 050406.
Swift executed a prompt slew.
The XRT imaged the BAT field only 84 s after the trigger
but no bright X-ray source could be detected within the field of view.
However, a refined on-ground analysis revealed a previously uncatalogued
X-ray source (Cusumano et al. 2005a; Capalbi et al. 2005).
From the very first examination of the down-linked data it was clear that the
afterglow of this burst was peculiar. Indeed, after an initial decay, the
X-ray count rate began rising, peaking at 220 s, and subsequently
decaying again (Capalbi et al. 2005).
Ground-based observations started as soon as the burst discovery was reported via the
GCN network.
The Magellan/Clay telescope imaged the XRT error circle
with LDSS-3 in the R and i bands and found a single faint source (
mag,
7.8 h after the burst) located at
RA(J2000
,
Dec(J2000
with an uncertainty of
in each coordinate (Berger et al. 2005b,a).
Similarly to XRT, UVOT also imaged the field at the end of the slew
(starting from
88 s after the trigger) and though it failed to detect the afterglow on-board
(Landsman et al. 2005),
subsequent on-ground analysis revealed a source within the XRT error circle at the
4.3- (19.0 mag), 3.0- and 2.5-
detection levels in the U, B and V bands,
respectively.
The UVOT position was RA(J2000
,
Dec(J2000
,
consistent with the
Magellan one.
By the time the second UVOT observation (1.3 h later) was performed,
the source was not detected in the U band, confirming it as the
afterglow of XRF 050406.
Schady et al. (2006) obtained
an estimate of
from fitting the broad band spectrum (combined UVOT and XRT data).
In this paper we present observations of the first Swift burst
where a flare is clearly detected in its X-ray light curve,
during which the source count rate increased by a factor of 6.
This feature had never been observed before in Swift data, and
had rarely been observed before in any X-ray afterglow (Piro et al. 2005).
This paper is organized as follows.
In Sect. 2 we describe our observations and data reduction;
in Sect. 3 we describe our spatial, timing and spectral data analysis;
in Sect. 4 we discuss our findings.
Finally, in Sect. 5 we summarize our 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.
The decay and spectral indices are parameterized as follows,
,
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).
Table 1: Observation log of XRF 050406.
Table 2: Spectral fit results.
Table 1 reports the log of the observations that were used for this work. The BAT data were analyzed using the standard BAT analysis software distributed within FTOOLS v6.0. The burst is detected in the first two standard bands (15-25 and 25-50 keV) while virtually no flux is observed above 50 keV. We find
The BAT spectra were extracted over the full time interval over which the burst was
detected (
), in the interval covering the 1-s peak
,
and for the
and
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 15-150 keV energy range.
All spectra are well fit with a simple power law with photon index
(see details in Table 2).
There is no evidence of a spectral break within the BAT energy range, thus
constraining the peak energy
keV.
The indices are steeper (softer) although consistent with the ones reported by
Krimm et al. (2005), due to the different energy ranges used for the spectral fitting.
No significant improvements are found using either a cutoff power-law or a Band model
(Band et al. 1993).
The 1-s peak photon flux was
erg cm-2 s-1 (15-350 keV band), while
the fluence was
erg cm-2 (15-350 keV band).
This fluence corresponds to an isotropic-equivalent energy
erg
(in the rest frame 52-1204 keV) assuming
(Schady et al. 2006).
![]() |
Figure 1:
XRT image of XRF 050406,
obtained from the total ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In order to cover the dynamic range and rapid variability expected from GRB afterglows and to provide rapid-response, automated observations, XRT was designed to support different readout modes that optimize the collected information as the flux of the burst diminishes. The XRT supports four major readout modes, one imaging (IM), two timing, Piled-up/Low-rate Photodiode (PuPD and LrPD) and Windowed Timing (WT), and one Photon-Counting (PC). A detailed description of XRT modes can be found in Hill et al. (2004). In the nominal operating state the mode switching is based on the source flux and is fully automated (auto state) to minimize pile-up in the data.
The XRT observations of XRF 050406 started on 2005 Apr. 6 at 16:00:12 UT,
only 84 s after the trigger, and ended on 2005 Apr. 22, thus summing up a
total net exposure (in PC mode) of 163 ks spread over a
16 d baseline.
The monitoring is organized in 9 observations (000, 001, 002, 005, 006, 008, 009, 010, 011)
and 183 snapshots (continuous pointings at the target).
This was the first burst to occur after the formal end of the calibration phase (2005 Apr. 5), and
the first (000) observation was performed as an automated target (AT) with XRT in auto state.
Therefore, during observation 000 the automated mode switching made XRT take an initial 2.5 s
image (IM at t=84 s), immediately followed by one PuPD (t=90 s) and one LrPD (t=91 s) frame.
Then at t=92 s a series of 5 WT frames was taken until the on-board measured count rate
was low enough for XRT to switch to PC mode (t=99 s).
After this, XRT repeatedly switched between WT and PC modes because of an increased
background level (see below).
Since the signal-to-noise (S/N) in these late WT frames is low, we did not include
them in our analysis (Table 1).
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.4.0) software package written by the
Agenzia Spaziale Italiana (ASI) Science Data Center and distributed within
FTOOLS v6.0 to produce the final cleaned event lists. We ran the
task xrtpipeline (v0.8.8) applying calibration and standard filtering and screening criteria,
i.e., we cut out temporal intervals during which the CCD temperature was higher
than
C, and we removed hot and flickering pixels. These are present
because the CCD is operating at a temperature higher than the design
temperature of -100
C due to a failure in the active cooling system.
An on-board event threshold of
0.2 keV (un-reconstructed pulse-height PHAS[1]>80)
was also applied to the central pixel, 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).
These two sources of background are the main reason for the switching between PC and WT mode
even when the source count rate is below 1 counts s-1.
Throughout the monitoring campaign the CCD temperature was <-50 C,
with the exception of part of observations 002 and 005,
where it became as high as -43.5 and -45
C, respectively;
those data were therefore screened out.
For our analysis we further selected XRT grades 0-12 and 0-2 for PC and WT data, respectively
(according to Swift nomenclature; Burrows et al. 2005a).
Figure 1 shows the 163 ks XRT image accumulated in PC mode in the
0.2-10 keV energy band.
We detected two previously uncatalogued sources within 1 arcmin of the optical
burst coordinates.
The brightest uncatalogued source, which we identified as the fading X-ray
counterpart of the burst, is present in the first four XRT snapshots.
The source is piled-up during the initial 500 s of PC data.
Therefore, to obtain an unbiased position, we rely on the
remainder of the PC data in the first observation, which has a
net exposure of 49.8 ks.
We used the xrtcentroid task (v0.2.7) and found that the
afterglow position is RA(J2000
,
Dec(J2000
.
We estimate its uncertainty to be 4
2 (90% confidence level).
This position takes into account the correction for the misalignment
between the telescope and the satellite optical axis (Moretti et al. 2006).
Figure 1 shows the XRT error circle, as well as the
3
BAT error circle (Krimm et al. 2005; 95% containment) and the
optical counterpart coordinates determined by Magellan (Berger et al. 2005a)
and by UVOT (Rol et al. 2005).
The XRT coordinates are 23
from
the BAT ones, and 1
6 and 2
8 from the Magellan and UVOT ones, respectively.
XRF 050406 was detected (XIMAGE v4.3) in the first four snapshots individually,
but not from the 5th on.
The second source, S2, is located at RA(J2000
,
Dec(J2000
and has a constant rate
(
counts s-1 throughout the observation
campaign.
![]() |
Figure 2:
X-ray light curve of the XRF 050406 afterglow in the 0.2-10 keV energy band.
The curve is background-subtracted and the time is referred to the BAT trigger, 2005 Apr. 06 at 15:58:48.4 UT (Parsons et al. 2005). The last point after 106 s is a 3-![]() ![]() |
Open with DEXTER |
Table 3: Light curve fit resultsa.
During the first 500 s of the XRT observation the intensity of the afterglow was high enough
to cause pile-up in the PC mode data. To account for this effect we extracted
the source events within an annulus with a 30-pixel outer radius (71
)
and a 2-pixel inner radius.
These values were derived by comparing the observed and nominal PSF.
For the PC data collected after the first 500 s, the entire circular region (30-pixel radius)
was used, instead.
In both cases we further disregarded data within a circular region centred on the serendipitous
source S2 (which lies within the 30-pixel PC source extraction region)
with a 7.17 pixel (17
)
radius.
The WT data were extracted in a rectangular region 31 pixels long along the image strip
(and 20 pixels wide), which excludes the data from the source S2.
The selected extraction regions correspond to
69% (piled-up PC),
95% (non piled-up PC),
and
94% (WT) of the XRT PSF.
To account for the background, data were also extracted in PC mode within
a circular region (radius
pixels)
and in WT mode within a rectangular box (
pixels),
in locations far from background sources.
The mean PC background in the 0.2-10 keV band was found to be constant throughout the observations
and, normalized to the PC source extraction region, it had a value of
counts s-1.
The mean WT background in the same energy band and normalized to the WT source region was
counts s-1.
Figure 2 shows the background-subtracted light curve extracted in the
0.2-10 keV energy band, with the BAT trigger as origin of time.
We considered WT data for the first snapshot of the first observation,
and PC data for all 9 available observations (see Table 1).
During the initial phases of the afterglow evolution (
s)
we binned the source counts with a minimum of 30 counts per time bin, and dynamically subtracted the
normalized background counts in each bin. The PC mode data were corrected for the effects of pile-up.
We note that, by keeping to the minimum number of counts per time bin criterion,
we created several bins during the first snapshot, but subsequently needed to merge data belonging
to snapshots 1 and 2 (point at
4 ks), then from snapshots 3 and 4 (point at
20 ks),
and later on from snapshots 5 through 8 (point at
35 ks).
Afterwards, we used XIMAGE with the option SOSTA, which
calculates vignetting- and PSF-corrected count rates within a specified box,
and the background in a user-specified region.
To ensure uniformity with the early light curve, the background was estimated in the same
region as the one used for the initial part of the light curve.
We thus obtained a signal-to-noise ratio
(the only exception being the point at
33 ks which has
). The last point is a 3-
upper limit.
This latter method is preferred for the construction of the late part of the light curve
since it better accounts for the background in a low-counts regime.
We note, however, that extracting the light curve in the same 30-pixel source region up to the
end of the last observation, we obtained fully consistent results, albeit with a noisier
light curve.
We also note that the residual contribution of the serendipitous source S2 within the source
extraction region is
of the S2 counts, which corresponds to
counts s-1.
Therefore, S2 only makes a marginal contribution to the afterglow light curve, which amounts to <20% of the last point.
The light curve clearly shows a complex behaviour, with a power law decay underlying
a remarkable flare which peaks at 210 s after the BAT trigger
(see Fig. 2, inset).
To fit the light curve we used the BAT trigger as reference time
and we only considered spectroscopic-mode data obtained while XRT was pointing,
thus excluding the early LrPD, the LrPD upper limit and the IM point.
Further excluding the data taken during the flare (180 s <t<300 s),
a fit with a simple power law yields
(12 degrees of freedom, d.o.f.),
which is unacceptable.
A fit with a broken power law
for
and
for
,
where
is the time of the break, yields
and
,
and a break at
4200 s after the
BAT trigger.
This latter model yields a good fit (
,
10 d.o.f.), a significant improvement
over the simple power law (null hypothesis probability
,
equivalent to 3.2
),
but some of the parameters are not well constrained.
Alternatively, a fit with two smoothly joined power laws
yields
(10 d.o.f.) with similar values for the inferred parameters.
A summary of the fits to the light curve can be found in Table 3.
As a reference, the 0.2-10 keV unabsorbed
flux at
is
erg cm-2 s-1(we adopted a count rate to unabsorbed flux conversion factor of
erg cm-2 count-1,
obtained from the best fit models derived in Sect. 3.3)
and the luminosity in the 0.7-34.4 keV band is
erg s-1.
![]() |
Figure 3: WT background-subtracted light curves. a): Total band ( T, 0.2-10 keV). b): Soft band ( S, 0.2-1 keV). c): Hard band ( H, 1-10 keV). d): Ratio of hard to soft count rates. |
Open with DEXTER |
During the flare a rebrightening of the source by a factor of 6 in flux was observed
between
s and the peak at
210 s.
Both the rising and the falling part of the flare had very steep slopes that,
when fit with a simple power law, yield
and
.
When the underlying power-law afterglow is subtracted,
the fit yields
and
and the peak is at
s from the BAT trigger.
In all cases the errors are dominated by the uncertainty in the placement of the flare boundaries.
The flare can also be characterised, as a simple parametric description,
as a Gaussian line. A combined broken power law and Gaussian model fit yields a
peak at
211.1+5.4-4.4 s
(
61.4+1.6-1.3 s in the rest frame)
and a width
17.9+12.3-4.6 s
(
,
17 d.o.f.).
In this case the ratio of the characteristic time-scale and the peak time is
or 0.20, when using the Gaussian width or its
FWHM (
42.2+29.0-10.8 s), respectively. In either case,
,
which puts severe constraints on the
emission mechanisms that can produce the flare.
We shall address this issue in the discussion section.
Integration of the Gaussian best-fitting function yields
an estimate of the fluence of the flare,
erg cm-2,
corresponding to an energy of
erg.
The large error reflects the uncertainty on the actual model used
for the integration of the flare.
We also extracted events from the first snapshot WT data in two more energy bands,
0.2-1 keV (soft, S) and 1-10 keV (hard, H), as well as the total band,
0.2-10 keV. We used the same regions as the ones described above,
a constant time binning of 30 s and dynamically subtracted their respective backgrounds.
Figure 3 shows the three background-subtracted light curves,
as well as the ratio H/S.
Indeed, during the rising portion of the flare the hard band flux increases by a factor of 6 while the soft band flux only increases slightly,
so that the spectrum of the flare starts off harder than the underlying afterglow,
and then evolves into a softer state as its flux decreases;
this can be seen in the following time bin, when the soft band flux peaks with a flare to
pre-flare flux ratio of
3.5.
This yields an indication of spectral evolution during the flare
as a
3-
excess over a constant fit to H/S.
It should be noted that this behaviour is reminiscent of that observed in the
prompt emission (Ford et al. 1995),
with the harder band peak preceding the softer band peak.
At
s a second faint bump is observed.
Its significance is not high, since it is detected as a 2-
excess over the underlying afterglow. Similar late-time bumps
have been observed in other Swift-detected GRBs
(e.g. GRB 050502B; Falcone et al. 2006).
The afterglow of XRF 050406 was very faint, hence it is not possible to perform
time-resolved spectroscopy to distinguish the spectral properties of the
afterglow proper from the ones of the flare observed in the light curve.
Therefore, we proceeded as follows.
Spectra of the source and background were extracted in the regions described in
Sect. 3.2 from the first observation (000) event files.
PC and WT spectra were extracted during the first 500 s
of the PC observation (see Table 2 for times referred to
T0), when PC data are piled-up and when the flare is observed in the light curve.
We also extracted PC spectra after the first 500 s during the
first 4 snapshots. For the latter we used a circular region with a 10 pixel radius
(corresponding to
80% of the XRT PSF) to minimize the background
and to be able to use the Cash statistics (Cash 1979).
Ancillary response files (ARF) were generated with the task xrtmkarf
within FTOOLS v6.0 using the latest ARF distribution (v003).
These ARFs account for different extraction regions and PSF corrections.
We used the latest spectral redistribution matrices (v007).
The energy ranges adopted for spectral fitting were 0.5-10 keV
and 0.2-10 keV for WT and PC, respectively.
We first performed a fit with an absorbed (wabs in XSPEC) power law
to the WT data (166 counts),
which were rebinned with a minimum of 20 counts per energy bin to
allow
fitting within XSPEC.
The hydrogen column was initially kept as a free parameter,
and then frozen to the Galactic value (
cm-2, Dickey & Lockman 1990)
when the fit yielded a value lower than (although consistent with) the Galactic one.
The fit was good,
for 6 d.o.f.,
and yielded
.
We then performed a fit with the same model to the remainder
of the PC data during snapshots 1 through 4
using Cash statistics which is more appropriate given the low number of counts
(21 un-binned counts) and calculated the goodness of the fit via 104 Montecarlo
simulations. The fit was good and yielded consistent results.
We also performed simultaneous fits to the WT and PC (60 counts) spectra extracted during the first
500 s (using
statistics) and of the PC data alone (using Cash statistics),
also obtaining consistent results.
Table 2 summarizes the results of the fits.
We note that, given the current goodness of the XRT calibration
(5% systematic uncertainty for all observing modes and grade selections in the
0.5-10 keV range; e.g., Romano et al. 2005), an excess of
cannot be excluded and we find a 3-
upper limit to the
total (Galactic plus intrinsic) hydrogen column along the line of sight of
cm-2.
We can therefore conclude that, during the first 600 s after the burst,
which include the X-ray flare observed in the light curve, the mean photon index is
,
and that the photon index does not vary after the end of the flare.
However, we do have clues regarding the presence of spectral evolution during the flare
coming from the hardness ratio analysis (Sect. 3.2),
even though the statistics are not high enough to show it in the spectral analysis.
As we will discuss later (Sect. 4.2), other afterglows with
larger amplitude X-ray flares demonstrate a strong spectral evolution
of the flares.
The duration of this burst (
s
in the 15-350 keV band) places this burst in the short
tail of the long GRB population (Kouveliotou et al. 1993).
Its fluence is relatively low (
erg cm-2 in the
15-350 keV band) but not unusually faint.
The gamma-ray characteristics of this burst are consistent with a classification
as an X-ray flash (Heise et al. 2001), or as an "X-ray rich GRB'' (XRR).
The softness of the observed spectrum, which is well fit in the 15-150 keV band with a simple power law with photon index
,
and with no significant emission observed above
50 keV,
implies that the peak energy is below the BAT bandpass (
keV).
The operational definition of XRFs/XRRs (e.g. Lamb et al. 2004) is
of a fast transient X-ray source characterized by a softness ratio
for an XRF and
for an XRR.
Extrapolation of the BAT spectrum, with the assumption of
keV, yields
,
which classifies this burst as an XRF.
However, a break in the spectrum may well be present in the 2-15 keV band.
In the most conservative case, i.e. assuming no flux below 15 keV, this event
would be an XRR GRB, with
.
The isotropic-equivalent gamma-ray energy of this event is
erg (Sect. 2.1),
and this effectively puts XRF 050406 in the low-energy tail of GRB energies
(Bloom et al. 2003b). Assuming that the Amati relation (Amati et al. 2002) holds,
we can infer a rest-frame
keV, which corresponds to
an observer-frame
keV.
This value is consistent with the nondetection of
in the BAT energy range.
To date, X-ray afterglows of XRFs have been detected in just a few cases (XRF 011030, XRF 020427: Levan et al. 2005; Bloom et al. 2003a; XRF 030723: Butler et al. 2004; XRF 040701: Fox 2004; XRF 050315: Vaughan et al. 2006). This is one of the first examples of a well-studied X-ray light curve of an XRF. Its main characteristics are not qualitatively different from those of normal GRBs (Chincarini et al. 2005; Nousek et al. 2006). As observations accumulate, it is becoming clear that these two classes of phenomena share many properties, and both have afterglows with similar characteristics (Sakamoto et al. 2005). This is a clue that both types of events may have a common origin and is supported by recent evidence that some XRFs are associated with supernovae (Fynbo et al. 2004; Bersier et al. 2005; Soderberg et al. 2004).
The general behaviour of the afterglow of XRF 050406 is a typical one.
The observed X-ray photon index (
)
is common among X-ray afterglows (De Pasquale et al. 2005; Chincarini et al. 2005).
The light curve shows a break from a relatively steep decay (
)
to a flatter one (
). Its overall shape is similar to the one
typically observed by the XRT (Chincarini et al. 2005; Nousek et al. 2006), even though
the initial slope is less steep than average.
The most striking characteristic of this burst is the strong flare in its X-ray
light curve, a feature which had never been detected by Swift before and had been previously
observed in very few GRBs
(GRB 970508, Piro et al. 1999; GRB 011121 and GRB 011211, Piro et al. 2005).
The fluence of the flare is
erg cm-2 in the
0.2-10 keV band, which amounts to
14% of the observed (15-350 keV band)
prompt fluence.
A better estimate of the flare-to-prompt energy ratio would require
the knowledge of the prompt spectral energy distribution (SED).
Since the actual peak energy of the prompt SED is unknown (
keV),
the extrapolation of the BAT fluence to the XRT band is highly uncertain.
For plausible values of
,
the flare to prompt fluence ratio
is in the 1-10% range.
The observed rebrightening is by a factor of 6 in flux, presents a peak at
s
and takes place on a very short timescale, with a
ratio of the characteristic time-scale and the peak time
.
Both the rising and the falling parts of the afterglow-subtracted flare had very steep slopes,
and
,
assuming the burst trigger as the time origin.
According to the standard relativistic fireball model, the prompt emission is
caused by internal shocks within the expanding fireball, while the afterglow is
produced by the fireball shocking the external medium (external shocks, Zhang & Mészáros 2004; Piran 1999).
Available models to explain flares include refreshed shocks (Rees & Mészáros 1998),
external shocks with a clumpy medium (Lazzati et al. 2002) and
angular inhomogeneities in the outflow (Nakar et al. 2003; Fenimore et al. 1999).
However, it can be argued (Burrows et al. 2005b, Zhang et al. 2005,
Nousek et al. 2006) that such models cannot produce the observed large
flux variations
in such short timescales
(Ioka et al. 2005).
Similarly, none of the above mechanisms would explain the steep slopes observed in the flare.
External reverse shocks, created when the fireball slows down because of the
interaction with the external medium, are expected to emerge at optical and radio
wavelengths, hence synchrotron self-Compton (SSC) must be invoked to produce emission
in the X-ray band. This would require carefully balanced conditions (Kobayashi et al. 2005).
Piro et al. (2005) suggested that the X-ray flares observed in GRB 011121 and GRB 011211 were due to the onset of the afterglow. The steep slopes and the short timescale variability can only be accounted for within the thick shell scenario (Sari & Piran 1999). Galli & Piro (2005) successfully modeled XRF 011030 using this model. In this scenario, the emission before and after the flare is due to different processes (prompt tail and afterglow, respectively), hence a discontinuity in the light curve is generally expected underlying the flare. This is not the case for XRF 050406, where the same component describes the X-ray emission both before and after the flare. Even if a fine-tuning may explain this particular event, the lack of a light curve break is common to a large fraction of the flares observed by Swift (Burrows et al. 2005c). Therefore, while the explanation of flares in terms of the afterglow onset is attractive, it is unlikely to be applicable to the vast majority of the X-ray flares seen by XRT.
A promising mechanism to produce the flare is late internal shocks
(Zhang et al. 2005; Perna et al. 2006; Fan & Wei 2005; King et al. 2005),
which implies that the central engine is still
active at t=213 s, even though the prompt emission ended after s.
The late-time activity in this case must have a reduced power with respect
to the prompt emission, as the relative fluences indicate.
Such a mechanism would naturally explain the steep rise and decay slopes.
Second, the energy required to power the flare would be much lower than in the
other scenarios (Zhang et al. 2005).
The indications of spectral evolution throughout the flare further support
this interpretation. The flare appears to be harder than the underlying afterglow,
which suggests a distinct origin for this emission.
Furthermore, there are indications of spectral evolution, which shows the typical
hard-to-soft pattern. Such a behaviour
is commonly observed in the prompt emission spikes of GRBs (e.g. Ford et al. 1995),
which are produced in internal shocks.
Further evidence of late engine activity comes from both the flat part of the
light curve (
,
see Sect. 4.3)
and possibly by the presence of the late-time bump
observed at
s.
Following the discovery of a flare in the afterglow of XRF 050406,
initially reported by Burrows et al. (2005b), many others were identified:
GRB 050502B (Falcone et al. 2006), GRB 050724 (Barthelmy et al. 2005b)
and GRB 050904 (Cusumano et al. 2005c), just to mention a few.
At the time of writing (2005 Oct.), 50% of the bursts detected by XRT which were immediately re-pointed towards showed flares,
making flaring quite a common behaviour.
Furthermore, all the characteristics of the XRF 050406 flares
have now been observed in most flaring GRBs
(see Burrows et al. 2005c for a recent review).
For example, highly significant spectral evolution throughout the flare
has been reported in GRB 050502B (which was the brightest observed so far)
and GRB 050724.
In several cases the flares present large amplitudes and occur on short timescales.
Furthermore, several flares are often observed in the same event,
at times ranging from
100 s to
104-105 s after the burst.
Finally, in most cases the afterglow is clearly present before the onset
of the flare, and has consistent decay slope and flux levels with after the flare.
The present case shows that flares are present both in XRFs and in GRBs.
Since flares are likely tied to the central engine activity, this finding
further supports the idea that a similar mechanism is at work for both kind
of events (Fan & Wei 2005).
The prompt reaction of Swift has allowed us to observe the X-ray
light curves of GRB afterglows starting from a few tens of seconds
after the burst explosion.
In most cases the X-ray light curves are characterized by an initial
steep decay (up to 500 s) followed by a shallow decay, and then by a
steeper decay with a second break normally occurring at a few thousand
seconds later (Chincarini et al. 2005; Nousek et al. 2006).
The early steep decay seen in the X-ray light curve can be explained as
the tail of the prompt emission (however, see Panaitescu et al. 2005).
The few cases where the XRT light curve
lies well above the extrapolation of the prompt emission into the X-ray band
can be explained either by a strong spectral
evolution or by an X-ray flare with the maximum
located before the XRT observation (Tagliaferri et al. 2005).
There are other instances where the first steep decay is not observed at all
(e.g. Campana et al. 2005).
In the case of XRF 050406, however, the initial slope is shallower than the
steep values
observed in other early afterglows
(Tagliaferri et al. 2005). Moreover, the curvature relation
(Dermer 2004; Kumar & Panaitescu 2000) is not
satisfied, even after taking into account the effects pointed out by
Zhang et al. (2005) that would alter such relation.
Therefore, we also investigate whether the initial decline seen in XRF 050406
is consistent with afterglow emission.
Comparison of spectral indices and temporal decay slopes with
theoretical relativistic fireball models (e.g. Table 2 in Zhang et al. 2005)
indicates that the first decay index
and
energy index
rule out fast cooling models
(for which the injection frequency
exceeds the cooling frequency
)
for
.
For
,
the
closure relation is
satisfied within the errors and an electron power-law distribution
index
is obtained.
The same relation holds for the slow cooling regime (where
)
for
(both wind and ISM). In this case a consistent solution is
also found for
,
although with a large
.
The ISM environment is favoured on the basis of a better satisfied closure relation.
In conclusion, the spectral indices and temporal decay slopes of the first part of the
X-ray curve can be interpreted in terms of relativistic fireball models,
even though the large uncertainties associated with the slopes do not allow us to choose
among the available models.
An alternative explanation for the initial XRT emission is the presence of an additional flare which started before the beginning of the XRT observation, and of which we only see the decaying part. The superposition of two (and possibly more, fainter) flares would then mimic the initial steep power law decay. However, this interpretation seems less likely since recent Swift observations of X-ray flares within the first several hundred seconds of the prompt emission all had temporal decay indices much steeper than the observed XRF 050406 pre-flare index.
At
s the XRT light curve breaks to
.
Such a flat decay cannot be explained in terms of the standard afterglow model.
The only possibility would be to observe, in the fast cooling regime, the segment with
(where
is expected, marginally consistent
with the observed value). However, the fast cooling regime is expected to end much
earlier. To maintain the observed decay unbroken up to
106 s,
large values of the equipartition parameters
and
or of the Compton parameter would be required.
We consider this possibility quite unlikely.
Another possibility is that the angular energy profile of the fireball is not trivial
(a structured jet), so that emission coming from the (brighter) wings of the jet may
increase the observed flux as the fireball Lorentz factor decreases (Panaitescu et al. 2005).
An interesting explanation for the shallow-decay phase is injection of new
energy into the fireball through refreshed shocks (Sari & Mészáros 2000; Zhang & Mészáros 2001).
For this to happen, the energy
release inside refreshed shocks must be sizeable, since the whole fireball dynamics
has to be modified. Assuming an energy injection rate
,
we find
q in the range 0 to 0.5 depending on the model details (Zhang et al. 2005).
In this model, the initial part of the XRT afterglow light curve can be due to standard
afterglow emission only if the fireball evolution is not influenced at these stages.
Indeed, the energy supply provided by refreshed shocks is steadily growing,
and at the beginning it cannot alter the fireball dynamics.
In this case, the break would identify the time when the new, injected
energy is comparable to the fireball energy.
On the contrary, if the first XRT phase were due to late engine activity,
then the energy injection could have begun much earlier
and its emission would have been masked.
Integration of the light curve from the onset of the flat slope phase yields
erg cm-2, where
is the time at which the shallow phase ends, for
which we can only set a lower limit.
We note that this depends weakly on the onset time of the shallow phase, therefore
the calculated fluence is correct in both presented scenarios.
For comparison, the amount of energy released during the steep phase of the light
curve (excluding the flare) is
erg cm-2.
We note that the shallow phase lasts a considerable time.
Zhang et al. (2005) propose three explanations for the energy injection mechanism.
In the impulsive case (Sari & Mészáros 2000),
the central engine ejects material with a wide distribution of Lorentz factors.
In this case, slower moving shells will catch the fireball at a later time.
We can estimate the minimum Lorentz factor as
,
where
erg is the isotropic-equivalent energy,
and n0 is the external medium particle density in units of cm-3.
This implies that the acceleration process works from ultra- to mildly-relativistic
velocities.
Within the putative Poyinting flux scenario (Zhang & Kobayashi 2005),
the energy supply is provided by the transfer of
magnetic energy to the fireball, and the time at which the injection stops is related to
the ratio
of the electromagnetic to baryonic kinetic energy.
If this scenario is correct, we can infer a lower limit
of
,
for q=0.5-0, where
is the start time of injection.
Therefore, after the end of the energy transfer phase, the energy of the blast-wave would
be increased by a comparable factor.
In the third scenario (the prolonged energy output by the central engine, Zhang & Mészáros 2001),
the end of the injection phase is simply the end of the engine activity.
In this case, this activity produces a large amount of energy, particularly so since
the radiative efficiency may be lower during the late afterglow than during the prompt
emission, as is generally the case.
This was previously noticed by Nousek et al. (2006) in a sample of several Swift GRBs.
The monitoring of XRF 050406 was discontinued 22 days after the trigger.
By then, the source was no longer detectable and only a
3-
upper limit could be drawn at
counts s-1.
In order for the afterglow energy not to diverge, a further, late break is necessary.
One interesting possibility is that this may be due to seeing the edge of the jet.
A steepening in the light curve is expected when the fireball Lorentz factor becomes
comparable to the inverse of the jet half-opening angle. Such a late break is not
unexpected for an XRF. The few XRFs with known redshift
(Fynbo et al. 2004; Bersier et al. 2005; Soderberg et al. 2004)
have a very low isotropic-energy release, and this may be at least in part accommodated
if they have very wide jets.
This picture is consistent with the result found by
Frail et al. (2001; see also Ghirlanda et al. 2004),
who found that low-energy GRBs tend to have wider opening angles.
Using the standard formalism (Sari et al. 1999; Rhoads 1999), the jet half-opening angle is
deg,
where
s is the jet break time and
is the
burst radiative efficiency.
Therefore, using our lower limit on the jet break time
s,
we can infer a lower limit on the jet half-opening angle of 16 deg.
This value is at the high end of the distribution of jet angles (Bloom et al. 2003b).
XRF 050406 is classified as an X-ray flash, with fluence
erg cm-2 (15-350 keV),
a soft spectrum (
), no significant
flux above
50 keV and a peak energy
keV.
Its main characteristics are however not qualitatively different from those of
normal GRBs. As observations accumulate, it becomes clear that these two
classes of phenomena share many properties, and both have afterglows with similar
characteristics. This is a clue that both events may have a common origin.
XRF 050406 is the first Swift-detected burst that showed a flare in its
X-ray light curve, a feature now found in 50% of the XRT afterglows.
The flare peaked at
210 s after the BAT trigger
(
61 s in the rest frame).
The best fit of the afterglow decay is obtained with a broken power law with
,
,
and a break at
4400 s after the BAT trigger.
The mean photon index is
.
During the X-ray flare a flux variation of
in
a timescale
is observed, and its measured fluence
in the 0.2-10 keV band is
erg cm-2
[
erg],
which corresponds to 1-15% of the prompt fluence.
We argued that the flare-producing mechanism is late internal shocks,
which implies that the central engine is still active at
s,
though with a reduced power with respect to the prompt emission.
We showed possible indications of spectral variations during the flare, and a
flattening of the X-ray light curve after
s in
support of continued central engine activity at late times.
Since XRF 050406 was observed, flares have been detected by XRT in both X-ray flashes and normal GRBs, indicating that flares are linked to some common properties of both kinds of bursts, and probably tied to their central engine.
Acknowledgements
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 instrument possible.