A&A 447, 145-156 (2006)
DOI: 10.1051/0004-6361:20053795
P. M. Vreeswijk1,2 - A. Smette1,3,4 - A. S. Fruchter5 - E. Palazzi6 - E. Rol2,7 - R. A. M. J. Wijers2,8 - C. Kouveliotou9 - L. Kaper2 - E. Pian6,10 - N. Masetti6 - F. Frontera11,6 - J. Hjorth12 - J. Gorosabel13 - L. Piro14 - J. P. U. Fynbo12 - P. Jakobsson12 - D. Watson12 - P. T. O'Brien7 - C. Ledoux1
1 - European Southern Observatory, Alonso de Córdova 3107,
Casilla 19001, Santiago 19, Chile
2 -
Astronomical Institute "Anton Pannekoek'', University of Amsterdam &
Center for High Energy Astrophysics, Kruislaan 403, 1098 SJ
Amsterdam, The Netherlands
3 -
Institut d'Astrophysique et de Géophysique, Université de
Liège, Allée du 6 Août 17, 4000 Liège, Belgium
4 -
Chercheur qualifié, FNRS, Belgium
5 -
Space Telescope Science Institute, 3700 San Martin Drive, Baltimore,
MD 21218, USA
6 -
INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica, Sezione di
Bologna, via Gobetti 101, 40129 Bologna, Italy
7 -
Department of Physics and Astronomy, University of Leicester,
University Road, Leicester, LE1 7RH, UK
8 -
Department of Physics and Astronomy, SUNY Stony Brook, NY
11794-3800, USA
9 -
National Space Science Technology Center, NASA/MSFC, XD-12, 320
Sparkman Drive, Huntsville, AL 35805, USA
10 -
INAF - Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11,
34131 Trieste, Italy
11 -
Dipartimento di Fisica, Università di Ferrara, via Paradiso 12,
44100 Ferrara, Italy
12 -
Dark Cosmology Centre, Niels Bohr Institute, University of
Copenhagen, Juliane Maries Vej 30, 2100 Copenhagen Ø, Denmark
13 -
Instituto de Astrofísica de Andalucía, CSIC, Apartado
3004, 18080 Granada, Spain
14 -
INAF - Istituto di Astrofisica Spaziale e Fisica Cosmica, Sezione di
Roma, via Fosso del Cavaliere 100, 00133 Roma, Italy
Received 8 July 2005 / Accepted 6 October 2005
Abstract
We present low-resolution VLT spectroscopy of the afterglow
of the gamma-ray bursts (GRBs) 991216, 011211 and 021211. Our
spectrum of GRB 991216 is the only optical spectrum for this afterglow. It
shows two probable absorption systems at z=0.80 and z=1.02,
where the highest redshift most likely reflects the distance to the
host galaxy. A third system may be detected at z=0.77. HST imaging
of the field, obtained 4 months after the burst, has resulted in the
detection of two amorphous regions of emission, one at the projected
afterglow position, and the other 0
6 away. The spectrum shows
a depression in flux in between 4000 Å and 5500 Å. This could
be the result of a 2175 Å-type extinction feature in the host of
GRB 991216, but at a rather red wavelength of 2360 Å. If this
interpretation is correct, it is the first time the extinction
feature is seen in a GRB afterglow spectrum. It is centered at a
wavelength similar to that of the ultra-violet (UV) bumps inferred
from observations of a few UV-strong, hydrogen-poor stars in the
Galaxy. All significant absorption lines (except for one) detected
in the spectrum of GRB 011211 are identified with lines originating in a
single absorption system at
,
the redshift of the
GRB 011211 host galaxy. We also detect the Ly
absorption line in the
host, to which we fit a neutral hydrogen column density of
log N(H I) =
,
which indicates that it is a
damped Ly
system. Using a curve-of-growth analysis, we estimate
the Si, Fe and Al metallicity at the GRB 011211 redshift to be
[Si/H] =
-0.9+0.6-0.4,
,
and
[Al/H] =
-1.0+0.5-0.3. For GRB 021211, we detect a single emission
line in a spectrum obtained tens of days after the burst, which we
identify as [O II]
3727 at z=1.006. The corresponding
unobscured [O II] star-formation rate is
1.4 M
yr-1.
Key words: gamma rays: bursts - galaxies: abundances - galaxies: distances and redshifts - galaxies: quasars: absorption lines
Table 1: Log of VLT spectroscopic observations.
Spectroscopy of the afterglows of long-duration Gamma-Ray Bursts (GRBs) has been essential for our understanding of the physical mechanism that produces these powerful explosions. Following the discovery of GRB afterglows (Van Paradijs et al. 1997; Costa et al. 1997), the first redshift determination by Metzger et al. (1997) provided conclusive evidence that the origin of the long-duration class of GRBs (Kouveliotou et al. 1993) is cosmological. Obviously, redshifts are required to deduce most meaningful quantities, such as the GRB energetics (Ghirlanda et al. 2004; Bloom et al. 2001; Frail et al. 2001), the brightness of any underlying supernova component (e.g. Bloom et al. 1999), and properties of the GRB host galaxies (e.g. Hogg & Fruchter 1999). So far, optical redshifts have been secured for four dozen GRBs, either through absorption-line spectroscopy of their bright early afterglow (e.g. Metzger et al. 1997), or through the detection of host-galaxy emission lines (e.g. Kulkarni et al. 1999). The average observed GRB redshift is z=1.3(see Jakobsson et al. 2004a; Berger et al. 2005), with the highest being GRB 050904 at z=6.29 (Kawai et al. 2005).
Apart from the vital redshift, spectroscopic afterglow observations
have also shown that GRB hosts are actively star-forming galaxies
(e.g. Djorgovski et al. 1998; Vreeswijk et al. 2001b), that
the afterglow is often situated behind a very large neutral hydrogen
column (e.g. Vreeswijk et al. 2004; Hjorth et al. 2003a) with relatively high
metal column densities and dust depletions compared to damped
Ly
systems observed along QSO sight lines
(Savaglio et al. 2003; Savaglio & Fall 2004), and
have provided evidence for high-velocity outflows (up to 3000 km s-1, Møller et al. 2002; Mirabal et al. 2003; Schaefer et al. 2003),
presumably caused by the wind of the GRB massive-star
progenitor. Finally, spectroscopic monitoring of GRB 030329/SN2003dh
(Stanek et al. 2003; Hjorth et al. 2003b) has revealed a remarkable similarity with the
spectral evolution of GRB 980425/SN1998bw
(Patat et al. 2001; Galama et al. 1998), confirming the
connection between GRBs and supernova explosions. This strongly
suggests that long-duration GRBs are produced by collapsing massive
stars (see Woosley 1993; MacFadyen et al. 2001).
The GRB Afterglow Collaboration at ESO (GRACE) has an on-going program at the European Southern Observatory (ESO) to perform spectroscopic observations of GRB afterglows. In this paper, we present low-resolution Very Large Telescope (VLT) spectroscopy of three GRB afterglows: 991216, 011211 and 021211. For GRB 991216 the spectroscopy is complemented by Hubble Space Telescope (HST) imaging. The organisation of this paper is as follows. After a description of the observations, data reduction and spectral analysis in Sect. 2, each GRB afterglow is presented in a separate section: GRB 991216 in Sect. 3, GRB 011211 in Sect. 4 and GRB 021211 in Sect. 5. Each of these sections starts with an introduction on the afterglow, followed by our results. We briefly conclude in Sect. 6.
Table 1 shows the log of the VLT spectroscopic
observations. All spectra have been reduced in the same manner within
IRAF, using tasks within the kpnoslit package. After overscan
subtraction and flat-fielding, the cosmic rays were removed from the
images using the L.A. Cosmic routine written by
Van Dokkum (2001). The spectra were then optimally extracted
for each 2-dimensional (2D) image separately, with an extraction
aperture width of 2-4
.
The wavelength calibration was applied,
again to each spectrum separately, using an HeNeAr lamp spectrum that
was taken in the morning after the science observations. The formal
error in the wavelength calibration fit (order 5) was roughly
0.3 Å for the 150I grism, and 0.2 Å for the 300V grism. The
individual wavelength-calibrated spectra were averaged, and the
corresponding Poisson error spectra, calculated by the apall
task, were quadratically averaged.
The flux calibration was performed using observations of the standard
HD 49798 for both GRB 991216 and GRB 011211, and LTT 3218 for GRB 021211. These standards
were taken with a 5
slit width during the same night as the
GRB science observations. During the night of observations of GRB 991216,
HD 49798 was also observed with a slit width of 1
,
i.e. the
same width as that of the GRB afterglow spectra. The nights that
GRB 991216 and GRB 021211 were observed were probably photometric, but the night
that GRB 011211 was observed was definitely not.
For the spectra of GRB 011211, we performed a correction for the slit loss,
i.e. the fraction of the surface underneath the spectral profile that
falls outside the slit width. The correction factor was estimated by
fitting a Gaussian profile along the spatial direction (i.e. along the
CCD columns) of the 2D spectra, averaged over 4 pixels in the
dispersion direction (i.e. averaging 4 columns before performing the
fit). The resulting Gaussian full width at half maximum (FWHM) was
then compared to the slit width to obtain the slit loss along the
dispersion axis. The slit loss profile was then fit with a polynomial
to correct the spectra. Note that this does not correct for any
colour-dependent slit losses, but both FORS1 and 2 have a linear
atmospheric dispersion compensator (LADC) in the light path, which
minimises any colour-differential slit losses up to a zenith distance
of 45.
However, several spectra of
GRB 011211 were taken at an airmass above 1.4 (see Table 1), and
therefore these are likely to suffer from such colour-differential
slit losses. For the GRB 991216 spectra a slit-loss correction was not
necessary, as we flux calibrated with the standard that was taken with
the same slit width as the afterglow spectra. This should in principle
result in a correct flux calibration, provided that the seeing was the
same during the afterglow and standard star observations. This was
more or less the case for GRB 991216, with a seeing of 0
6 during the
afterglow observations, and 0
7 during the standard star
observations. For the late-time spectra of GRB 021211, where the afterglow
continuum is absent, we did not attempt a slit-loss correction.
Finally, each spectrum was corrected for Galactic extinction, assuming the EB-V from Schlegel et al. (1998). The EB-V values used are 0.626 for GRB 991216, 0.043 for GRB 011211, and 0.028 for GRB 021211. The flux calibration, slit-loss correction, and dereddening were also applied to the combined wavelength-calibrated error spectrum.
The resulting spectra were analysed within IRAF as well. The equivalent width (EW) of the significant absorption and emission lines were measured using the splot routine in IRAF. By eye, we fit the continuum with a high-order polynomial function. After normalisation, we measure the EWs with the e option. The Poisson error in the EW is calculated from the formal error spectrum calculated within IRAF, using the same pixels that were used to measure the EW. We do not include a contribution in the error from the uncertainty in the location of the continuum.
Following the Burst And Transient Source Experiment (BATSE)
-ray detection of GRB 991216 (see Kippen et al. 1999), the
X-ray afterglow was discovered by the Rossi X-ray Timing
Explorer (RXTE) (Takeshima et al. 1999), and the optical/infrared
afterglow by Uglesich et al. (1999). The afterglow of GRB 991216 is one
of a number of GRB afterglows (as e.g. GRB 990510;
see Harrison et al. 1999; Pian et al. 2001) to show evidence for a
"beaming'' break (Rhoads 1997), suggesting that the
gamma-ray, X-ray, and early optical emission was confined into jets
(Halpern et al. 2000). Although the optical through X-ray data
can be explained with a jet fireball model, inclusion of the radio
data calls for more exotic models (Frail et al. 2000). Chandra X-ray Observatory spectral observations of the X-ray
afterglow show two probable emission features (a 4.7
detection at
keV, and a marginal one at
keV; Piro et al. 2000), that are identified with an
iron line and the recombination continuum at a redshift of
.
These features suggest the presence of 0.01-1
of iron in the vicinity of the burst
(Vietri et al. 2001; Piro et al. 2000). However,
in an independent analysis of the same Chandra data,
Sako et al. (2005) only find a
fluctuation at a slightly
higher energy (3.8 keV), and do not find any feature near 4.4 keV.
![]() |
Figure 1:
Top panel: the optical spectrum of the afterglow of
GRB 991216, taken 1.5 days after the burst. The resolution is approximately
24 Å across the entire spectrum, which corresponds to a resolving
power of 250 at 6000 Å. The spectrum has not been smoothed. The
blue and red spectra (respectively without and with the OG590 order
sorting filter) are connected at 6200 Å. The two strong absorption
lines longward of 6000 Å are telluric. The spectrum flux has been
increased by 10% to match the ![]() ![]() ![]() |
Open with DEXTER |
The spectrum of the afterglow of GRB 991216 is shown in
Fig. 1. The two strong features in the red part
around 6900 Å and 7600 Å are caused by the earth's atmosphere
and are not intrinsic to the GRB host. The apparent lines above
9000 Å are due to bad sky subtraction. In the blue part of the
spectrum, several significant absorption lines are detected. From
inspection of the standard star spectrum at the same wavelengths we
can conclude that the detected lines are not caused by the instrument,
or by the earth's atmosphere. The detected lines cannot be identified
with typical interstellar medium (ISM) absorption features of a system
at a single redshift; at least three systems need to be invoked. An
alternative explanation is that the feature around 5700 Å is
Ly
at z=3.6, and that the features shortward of it are
Ly
forest lines. We have simulated such a spectrum, and the
lines can appear quite similar to the lines in the observed
spectrum. This is a very unlikely explanation, however, since the
break is much too smooth, and the continuum emission around 4000 Å is only about 10-20% lower than the level beyond the supposed
break. We therefore discard the latter possibility, and are left with
the most probable explanation: several absorption systems.
Table 2 lists the identification of the
significant lines in the spectrum of GRB 991216. These significant lines are
indicated with the short solid lines in the lower panel of
Fig. 1. As the spectrum was taken with a very
low-resolution grism (150I), we conservatively use 5
as the
significance threshold. Typical ISM absorption lines that could be
present at low significance, such as Mg I
2852, Ca II K
and H, and other Fe lines, are shown with a short dashed line in the
lower panel of Fig. 1.
On the basis of the strongest features we infer the possible presence
of three absorption-line systems along the line of sight, with the
following redshifts: z=0.77, z=0.80 and z=1.02. The
identification of the z=1.02 system in the blue part of the spectrum
is strengthened by the possible detection of Ca II K (3933 Å)
and H (3968 Å) around 8000 Å, respectively at z=1.021 and
z=1.022. In order to increase the signal-to-noise ratio (S/N) around
the possible Ca lines, we have added the blue spectrum (taken without
the order sorting filter) to the red spectrum (taken with the OG590
order sorting filter), even though the blue spectrum can be
contaminated by its second order above
Å. This second
order appears with twice the resolution and at twice the wavelength of
the first order. However, if the two features around 8000 Å were
due to this contamination, their first order counterparts would be
located at 3976 Å and 4009 Å, where no lines are detected. The
resulting combined spectrum around the Ca lines is also shown in
Fig. 1, offset downward by 30
Jy from the
red-only spectrum for clarity. We measure the following equivalent
widths for the two features in the combined spectrum:
(Ca II K
3933) = (
) Å and
(Ca II H
3968) = (
) Å. The expected ratio of
this doublet on the basis of the oscillator strengths,
,
is roughly two if the
lines are not saturated, but can be unity if they are saturated. The
observed ratio,
,
indicates that these lines are saturated.
Table 2:
Lines detected above 5
in the spectrum of GRB 991216.
A redshift of z=1.02 for GRB 991216 is consistent with the inferred
iron line detection in the Chandra spectrum of
the X-ray afterglow of GRB 991216 (Piro et al. 2000). However, we
note that in an independent analysis of the Chandra data,
Sako et al. (2005) find this iron line to be not significant. The z=1.02system that we detect in our optical spectrum is most likely caused by
the GRB 991216 host-galaxy ISM. In all cases for which both GRB absorption
and emission lines have been detected, the (most distant) absorption
system is at the same redshift as the emission lines from the presumed
host galaxy (e.g. for GRB 980703 and
GRB 021004; Djorgovski et al. 1998; Mirabal et al. 2003). Moreover,
if the detection of Mg I at z=1.02 is real, this suggests a
dense environment at this redshift, probably the ISM of the host
galaxy. The probable detection of the Ca II doublet supports this
hypothesis, since this line is believed to require column densities
that are typically an order of magnitude larger than seen in Mg
II absorbers (Carilli & Van Gorkom 1992; Bowen 1991). We note
that the z=0.80 system shows the strongest Mg II and Fe II
features.
The continuum flux level of optical afterglows, measured from
spectroscopy or broad-band photometry, is generally well described by
a single power law, with a slope of the order
(with
)
(for a high S/N example,
see Savaglio & Fall 2004). This is consistent with the fireball
afterglow theory (e.g. Sari et al. 1998). Any departure
from the pure power law is usually ascribed to extinction caused by
the host-galaxy ISM (see Ramaprakash et al. 1998). The
flux-calibrated spectrum of GRB 991216 contains two features that distinguish
it from "normal'' afterglow spectra: its shallow slope (
over the range 6000-9300 Å), and an apparent depression around
4700 Å.
In order to verify our absolute flux calibration, we first compare the
spectrum with the photometry of Halpern et al. (2000). Their
measurements are shown in Fig. 1; the
horizontal error bars represent the approximate FWHM of the filter
transmission curves. In absolute terms, the spectral flux level is
roughly 10% below that of the photometry; we have scaled the spectrum
in Fig. 1 upward with this amount. Relatively,
however, our flux calibration is in good agreement with the
photometry. There are no published U or B photometry data to
compare with. A cause for concern is that the spectrophotometric
standard that was used for the GRB 991216 flux calibration, HD 49798, has not
been calibrated in the optical regime with an observed spectrum, but
rather with a model extension from the observed 1150-3200 Å range
(see Bohlin & Lindler 1992). The same standard was also observed for GRB 011211
(see Sect. 4.1). Comparison of our GRB 011211
spectrum (see Fig. 3) with the photometry from
Jakobsson et al. (2003) indicates that the absolute flux
calibration of the GRB 011211 spectrum is too high by 70%. Indeed, the night
when this spectrum was taken was definitely not photometric
(alternating thin and thick cloud conditions). However, after scaling
the spectral flux down by this 70%, the relative spectral flux
calibration of our GRB 011211 spectrum is in good agreement with the
photometric calibration. In the B band, the spectral flux is
% higher than the photometry. If this is due to an error in
the model flux calibration of HD 49798, and we were to apply a
correction to the GRB 991216 spectrum based on this difference, it would only
make the GRB 991216 flux depression more pronounced.
Also, during the night that the GRB 991216 afterglow spectra were taken, the
standard HD 49798 was observed with a slit width of both 1
(in
long-slit mode, or LSS) and 5
(in multi-object mode, or MOS;
see the FORS manual
for details). This allows us to check whether the depression could be
due to colour-differential, or wavelength-dependent slit losses.
Flux-calibration of the GRB 991216 spectra with the 5
standard
results in an absolute offset of about 13% with the
1
-flux-calibrated spectra, with a lower flux level for the
5
-calibrated spectra. Relatively, the spectra agree to within
4%, with the depression being more pronounced in the
5
-calibrated spectra, i.e. there is no dependence on wavelength
(e.g. due to colour-differential slit losses) beyond the 4% level.
We note that Fig. 1 shows the
1
-calibrated spectrum. Moreover, we have also performed the
flux calibration using observations of three different 5
standards taken during three different photometric nights, up to one
week before the night of the GRB. The sensitivity curves of all these
standards differ by at most 8% in absolute terms (wavelength range:
4000-6000 Å) with the 5
-standard taken during the
GRB 991216 night. Therefore, we regard it unlikely that the blue flux
depression in the spectrum of GRB 991216 is caused by an error in the flux calibration.
Assuming that the relative flux calibration is acceptable, and that
the redshift of GRB 991216 is z=1.02, we investigate whether the depression
in the blue part of the spectrum can be due to a redshifted 2175 Å absorption feature, as observed in the Galactic extinction curve. With
this aim, we fit the Milky Way (MW) extinction curve of
Pei (1992) to the regions of our spectrum
void of identified absorption lines, and fixing the redshift at
z=1.02. We assume that the intrinsic afterglow spectrum is a single
power law, to obtain the intrinsic spectral slope ()
and
host-galaxy extinction (
)
that best matches the observed spectrum
(i.e. minimum
). For the MW we obtain the following best-fit
values:
,
corresponding to a rest-frame V-band
extinction of
mag. The resulting fit is shown by the
smooth solid line in Fig. 1. Fitting the red
part (6000-9300 Å) of the GRB 991216 spectrum with a single power law model
results in a slope of
.
We have also fit the
Large and Small Magellanic Clouds' (LMC and SMC) extinction curves;
these are shown by the dashed (LMC) and dotted (SMC) curves in
Fig. 1. These MC fits are less successful due
to the 2175 Å feature being less prominent in the LMC as compared
to the MW, and absent in the SMC. Although the MW fit is reasonable,
the 2175 Å feature clearly does not correspond to the minimum
depression in the spectrum. Placing the burst at a redshift of
z=1.19, or alternatively, moving the peak of the bump from 2175 Å to 2360 Å results in the best fit. This fit is shown by the
dash-dotted line in Fig. 1.
Although it seems likely that some form of graphitic carbon (possibly
polycyclic aromatic hydrocarbon - PAH - molecules) is responsible
for the 2175 Å absorption feature in our Galaxy
(Draine 2003), its exact nature is still unclear after
decades of research. Observations along different Galactic sight lines
have shown that although the FWHM of the 2175 Å bump can vary
considerably (in our GRB 991216 fits we did not allow this width to change),
the central wavelength of the bump is very stable
( Å; Draine 2003; Fitzpatrick & Massa 1986).
However, in a few particular cases where the UV component of the star
used to measure the interstellar extinction is strong enough to reveal
extinction by grains in its own circumstellar material, a broad bump
is observed with a central wavelength in between 2300 Å and
2500 Å. Abell 30 (Greenstein 1981) and HD 213985
(Buss et al. 1989) are two examples of such UV-strong,
hydrogen-poor stars. It is tempting to make the connection with the
UV-strong environment of a GRB, which may cause the extinction bump to
be redder than 2175 Å as well. Also, laboratory experiments
(see Blanco et al. 1991,1996), aimed at
testing the physical grain models proposed by
Hecht (1986) and Sorrell (1990),
have shown that, by decreasing the hydrogen content in carbon grains,
the 2175 Å bump not only becomes more pronounced, but also shifts
to higher wavelengths, up to 2600 Å (Mennella et al. 1995).
Beyond the Local Group, the best case for detection of the Galactic
extinction feature is a z=0.83 gravitational lens system, where the
best-fit central wavelength of the bump is found to be only slightly
larger than the Galactic value:
Å (Motta et al. 2002). Therefore, the flux depression in
our spectrum may have been caused by such a red extinction
bump along the line-of-sight to GRB 991216.
The field of GRB 991216 was observed with the Space Telescope Imaging
Spectograph (STIS) approximately 4 months after the burst, on 17 April
2000, starting at 11:36 UT, through the clear (50CCD) and long pass
(LP) filters, each for a total of 4790 s. The pipeline reduced
images were drizzled (see Fruchter & Hook 2002) onto output
images with pixels one-half native scale, or approximately 0
025
on a side. Figure 2 shows the central
square arcseconds of the sum of the 50CCD and LP images.
![]() |
Figure 2:
The
sum of the HST/STIS 50CCD and LP images of the field of GRB 991216. The
position of the early optical transient is marked with a black
circle, with an error radius of 0
![]() ![]() |
Open with DEXTER |
We have projected the optical afterglow (OA) position from an early
VLT image, taken 1.5 days after the burst, to the frame of the
HST drizzled images. Four bright nearby reference stars were
used, and the estimated 1
error in the resulting position is
0
05, corresponding to 2 drizzled pixels. The position and its
error are indicated with a circle in Fig. 2. The
error circle coincides with one of two faint regions of light (A and
B), which are separated by 0
6. The two regions are about equally
bright in the 50CCD as in the LP image, showing that they have similar
colours with most of the light coming out longward of 5500 Å. A and
B are possibly part of the same galaxy, or they could be two systems
that are interacting. Several other host galaxies also show a complex
morphology, e.g. the hosts of GRBs 980613, 990123, 000926, and
011211. The likely detection of the Ca II doublet at z=1.02argues slightly in favour of this merger scenario, since (at low
redshifts) the Ca II lines are detected mainly in disrupted
environments (Carilli & Van Gorkom 1992; Bowen 1991). Another
possibility is that B is actually a foreground system at z=0.80,
corresponding to the absorption-line system detected in the VLT
spectrum, and that A is the host galaxy of GRB 991216 at z=1.02. The
detection of emission lines from A and B would determine which of
these possibilities is correct. Either one of the two other galaxies
that are located roughly 2
away from the OA position could be
responsible for the tentative absorption system at z=0.77. At a
redshift of 0.77, 2
corresponds to about 14 kiloparsec
(assuming H0=65 km s-1 Mpc-1,
,
and
), which is smaller than the typical galaxy halo
size, and so one would indeed expect to detect Mg II in
absorption.
![]() |
Figure 3:
The optical spectrum of the afterglow of GRB 011211, smoothed with
a boxcar the size of 3 pixels, or about 8 Å. The resolution is
approximately 10 Å across the entire spectrum, which corresponds
to a resolving power of 600 at 6000 Å. The spectrum flux level has
been scaled downward by 70% to match the photometric
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Using an aperture radius of 0
2, we measured the following fluxes
in photons per second in the drizzled 50CCD and LP images:
for A in 50CCD,
(B, 50CCD),
(A, LP), and
(B, LP). We also measured the total flux of
A and B combined, with an aperture radius of 1
1:
(A+B, 50CCD), and
(A+B, LP). This aperture size is the
same as used by Halpern et al. (2000), who obtained
for the host galaxy from a Keck II image taken on 4
April 2000, i.e. about two weeks before the HST imaging reported
here. To convert the HST fluxes to standard magnitudes, a spectral
energy distribution has to be assumed for the source. Using synphot within IRAF, combined with the Kinney-Calzetti atlas of
galaxy spectra (Kinney et al. 1996; Calzetti et al. 1994), we
found that the best-fitting galaxy template is the Sc template
(redshifted to z=1.02). The redshifted starburst templates also
provide reasonably consistent magnitudes for both passbands and a
total magnitude corresponding with the value of
of Halpern and colleagues. Adopting the Sc galaxy template, the fluxes
above can be converted to the following magnitudes:
,
and
,
and
colours:
,
,
and
.
We estimate the uncertainty in the
transformation of fluxes to magnitudes to be of the order of 0.2 mag.
We note that these magnitudes (including the one of Halpern et al.)
and colours have not been corrected for the Galactic extinction along
the GRB 991216 sight line (
and AR = 1.7).
The transient afterglow may still be present in these observations,
but the low S/N does not allow an unambiguous identification of the
bright patch at the edge of the galaxy as a point source. We estimate
that any remaining OA is not brighter than R=27.6. Assuming the
single power law decay index,
,
of
Garnavich et al. (2000), the expected magnitude of the afterglow
at the time of our observations is
.
Our observations
therefore tend to confirm the break in the light curve reported by
Halpern et al. (2000). A supernova as bright as SN1998bw at a
redshift of z=1.02 would have R>30 at the epoch of our
observations, and would thus be too faint to be detected.
Following the BeppoSAX localisation of GRB 011211
(Frontera et al. 2002; Gandolfi 2001), the afterglow was
discovered by Grav et al. (2001) (see
also Jakobsson et al. 2004b), and a preliminary redshift of z=2.14was reported by Fruchter et al. (2001), which was confirmed by
Gladders et al. (2001). Rapid optical variability (1 hour
variability time scale at 12 hours after the burst) is detected in the
afterglow lightcurve; these variations could be caused by either
density fluctuations in the external medium
(Jakobsson et al. 2004b; Holland et al. 2002), or by a non-uniform
jet structure (Jakobsson et al. 2004b).
Reeves et al. (2002) obtained an XMM-Newton
spectrum of the X-ray afterglow and detect transient features that
they identify as Mg XI, Si XIV,
S XVI, Ar XVIII, and Ca XX at a redshift of
,
i.e. blue-shifted by about 0.1c with respect to the
optical redshift. The authors suggest this to be due to reprocessing
of the burst flux in a shell expelled by a supernova explosion
preceding the GRB by 4 days. However, Rutledge & Sako (2003)
(but see Reeves et al. 2003) and Sako et al. (2005) claim that
the detection of these metal lines in the X-ray spectrum is not
significant. The host galaxy has
,
with a modest
(
)
SMC-like extinction, and the most likely afterglow
model is that of a jet expanding into a constant mean density
environment (Jakobsson et al. 2003). Fynbo et al. (2003) detect Ly
in emission in the GRB 011211 host galaxy, extending up to 1
to the
North-North-East of the afterglow position.
Figure 3 shows the spectrum of GRB 011211. As was
discussed in Sect. 3.2, the night in which the
GRB 011211 afterglow spectroscopy was performed was not photometric. Comparison
with the photometric calibration of Jakobsson et al. (2003) shows
our calibration to be off by 70%, and we have scaled the spectrum
down by this amount in Fig. 3. The photometric
observations of Jakobsson et al. (2003), scaled in time to the
epoch of the spectrum, are shown by the open circles in
Fig. 3. The horizontal error bars reflect the
FWHM of the filter transmission. Relatively, the spectral and
photometric flux calibration match reasonably well. Fitting a single
power law to our spectrum over the range 4000-8000 Å results in a
spectral slope of
.
Table 3 lists the lines detected with a
significance above 3
in the spectrum. The redshift as
determined from the detected low-ionisation metal lines is:
.
Holland et al. (2002) have reported the
detection of 8 metal absorption lines at a mean redshift of
,
in a spectrum taken with the LDSS-2 imaging
spectrograph at the Magellan 6.5 m telescope. Although the S/N of our
spectra is superior, we are able to confirm the detection of only 4 of
the lines that Holland et al. (2002) identify. Their
Si IV
1368, Cr II
1406, 1422 are
misidentifications, and the line at 6875.2 Å identified as
Fe III
2189 by Holland and colleagues, is a well-known
atmospheric absorption feature. Also, we clearly do not detect the
broad absorption feature that they detect at 4600 Å.
Table 3:
Lines detected above 3
in the spectrum of
GRB 011211.
Figure 4 shows the fit to the Ly
line in the
spectrum of GRB 011211, performed within the LYMAN context in Midas. The
resulting neutral hydrogen H I column density is
log N(H I) =
,
with the line center at
(
) Å, corresponding to a redshift of
.
The fit region is limited to around the core of the
line, 3810-3833 Å. If this is extended to the red, the continuum
around 3850 Å is very poorly fit, and to the blue of 3810 Å the
Ly
forest starts to contaminate the Ly
line.
The H I column density in the host of GRB 011211 is at the low end of
the H I GRB afterglow column densities obtained so far. Of the
7 afterglows to date for which it was possible to detect Ly
in
absorption at the host-galaxy redshift, 6 turn out to be damped Ly
systems (see Vreeswijk et al. 2004; Jensen et al. 2001; Hjorth et al. 2003a)
and therefore have log N(H I) > 20.3. The host of GRB 021004
is the exception, with a column density N(H I)
atoms cm-2.
![]() |
Figure 4:
The
fit (solid continuous line) of the Ly![]() ![]() ![]() ![]() |
Open with DEXTER |
The low resolution of our spectrum makes it difficult to fit an
absorption profile to the detected metal lines to obtain the
corresponding column density and metallicity. However, we can derive a
strict lower limit on these quantities by simply assuming that the
lines are not saturated, even though most lines are in fact
saturated. In this optically thin approximation there is a linear
correspondence between the equivalent width of a line, and its column
density: log
= log
,
where f is the oscillator strength
(see Morton 2003), and the unit of
is
Å. Applying this relation to the detected low-ionisation lines that
are not blended, we obtain: log N(Si II)
14.6,
log N(Fe II)
14.6, and log N(Al II)
13.4. Assuming that these ions are the dominant ionisation state for
the corresponding element, these column density limits can be
translated to the following metallicity limits:
,
,
and
,
where we have used the solar
abundances from Grevesse & Sauval (1998) (see
also Morton 2003), and the H I column density
derived above. Neither Zn II
2026 nor Zn II
2062
are detected, for which we measure:
Å and
Å. These values correspond to
the 2
column density limit: log
and the
metallicity limit:
.
Taking this a step further, the abundances can be estimated through
the curve-of-growth (COG) technique
(see Spitzer 1978). Our COG analysis is very similar
to the one that Savaglio et al. (2003) applied to three GRB host
galaxies. Figure 5 shows the resulting curve of growth
for the absorption lines measured in the afterglow of GRB 011211. The straight
line corresponds to the linear part, where the lines are optically
thin. The bent curve shows the deviation from the optically thin
regime due to saturation of the absorption lines, the so-called flat
part. The degree of bending depends on the amount of Doppler
broadening (b in km s-1) in the medium in which the lines
originate. The best fit curve is obtained by varying the column
density of each element separately (i.e. sliding the data points
horizontally), and finding the value for the Doppler width b that
fits best with the overall shape of all elements together. We have
assumed a Gaussian shape for the line profiles, i.e. we have neglected
damping or natural line broadening, which is important for high column
density absorbers. We only use the low-ionisation lines from
Table 3 which are not obviously blended, i.e.
Si II 1526, Fe II
1608, 2344,
2382, 2586, 2600, and Al II
1670. We also included the
measurement of Fe II
2374, even though this line is detected
at the 2.5
level only; for this reason it is not listed in
Table 3. Minimising
,
we find the
following values for the Doppler width (b) and the
Si II, Fe II and Al II column densities:
b =
70+20-10 km s-1, log N(Si II) =
15.1+0.6-0.3,
log N(Fe II) =
,
and
log N(Al II) =
13.9+0.5-0.3. The errors in these values
are obtained by varying each parameter (while fitting the others), up
to the point where the chi-squared becomes the minimum chi-squared
plus one, i.e. where
+1. The full range
that the parameter can adopt within this chi-squared restriction is
taken as the 1
error range. These metal column densities
correspond to the following metallicities:
[Si/H] =
-0.9+0.6-0.4,
,
and
[Al/H] =
-1.0+0.5-0.3. Again, we have made the common
assumption that these low-ionisation species are the dominant
ionisation states for the corresponding element.
Finally, we do not see the Ly
in emission that was detected by
Fynbo et al. (2003), either in the 2-dimensional spectra or in the
extracted spectra. This is due to the poor S/N of our spectrum around
3800 Å.
![]() |
Figure 5: The curve of growth for the absorption lines measured in the afterglow of GRB 011211. |
Open with DEXTER |
The High Energy Transient Explorer 2 (HETE-2) satellite detected
GRB 021211 (HETE-2 trigger #2493) at 11:18:34 UT on 2002 December 11, as a
bright and X-ray-rich burst (Crew et al. 2003). The optical
afterglow was imaged as early as 65 s after the GRB
(Wozniak et al. 2002), and was reported by
Fox et al. (2003) to the GRB community within an hour of the
burst. The afterglow was fainter than nearly all known afterglows to
that date at an epoch of 1 day after the GRB (see Fig. 2
of Fox et al. 2003). Li et al. (2003) find evidence
for a light-curve break at min, with the decay index
becoming more shallow from -1.82 to -0.82, which they suggest is due
to the emission before the break being dominated by the reverse shock.
Panaitescu & Kumar (2004) calculate the shock microphysical
parameters, and find that the reverse-forward shock scenario provides
a more natural explanation than the wind-bubble scenario for the steep
early decay. Finally, Della Valle et al. (2003b) present
spectroscopic evidence for a supernova component in the late-time
afterglow of GRB 021211.
Our team first reported a redshift of z=0.800 for the putative host
galaxy of GRB 021211, based on the identification of the [O II] and
[O III] emission lines (Vreeswijk et al. 2002b). However, we
later found that, due to saturation of the offset star, the slit
position did in fact not cover the afterglow position
(Vreeswijk et al. 2002a). Instead, the slit mostly contained the
galaxy 1
5 to the North-East of the afterglow, first noted by
McLeod et al. (2002), whose redshift is then z=0.800,
and which is unrelated to the GRB as shown below.
In follow-up spectroscopy of GRB 021211, we positioned the slit to cover both
the position of the early afterglow and that of the galaxy at
z=0.800, as well as a nearby reference star. Using this reference
star and the slit position angle, the exact location of both the
galaxy and the projected position of the early-time afterglow on the
2-dimensional spectrum can be inferred.
Figure 6 shows the spectrum extracted at the
position of the early afterglow; we find one clear emission line at
7478 Å. This line most likely corresponds to [O II] 3727
emission from the host galaxy of GRB 021211 at a redshift of z=1.006(Vreeswijk et al. 2002a). This redshift was later confirmed by
Della Valle et al. (2003a) who besides [O II] also detect
H
4861 and [O III]
4959, 5007. The limited
spectral coverage in the red part of our spectrum does not allow us to
detect these lines. For completeness, we show in
Fig. 7 the wavelength-calibrated spectrum of the
galaxy that is located 1
5 to the North-East of the afterglow,
with a redshift of z=0.800.
![]() |
Figure 6:
The optical spectrum of the probable host galaxy of GRB 021211. The
spectrum has neither been flux calibrated nor smoothed. The
resolution is approximately 11 Å across the entire spectrum,
which corresponds to a resolving power of 550 at 6000 Å. The top
panel shows the S/N spectrum, smoothed with a boxcar with a size
similar to the spectral resolution. Only 1 significant line is
detected, at 7478 Å, which is most likely [O II] ![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
The optical spectrum of the galaxy 1
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
Della Valle et al. (2003b) obtained a total of 8 hours of
spectroscopy with the VLT/FORS2 and the 150I grism and OG590 filter
combination, which covers the optical spectrum above 6000 Å with a
resolution of 19 Å. Della Valle et al. (2003b) focus on the
possible evidence for a supernova component in the late-time afterglow
of GRB 021211, and briefly mention the detection of H
4861 and
[O III]
4959, 5007. Since the host-galaxy metallicity can
be measured from the combination of [O II], H
4861
and [O III]
4959, 5007 (see Kobulnicky et al. 1999),
we decided to reduce and extract these archival spectra as well.
However, we are unable to detect any significant line apart from the
line at 7478 Å mentioned above. This is probably due to a
combination of the very low spectral resolution, which complicates the
subtraction of the numerous sky emission lines in the red, and the low
sensitivity around 10 000 Å, where the [O III] and H
lines are expected.
Figure 6 shows that the continuum level just
redward of the significant line at 7478 Å is similar to that on the
blue side, suggesting that the line is not
Ly
at z=4.15. The most likely identification is [O II]
3727 at a redshift of z=1.006. This line is actually a doublet with
components at 3727.09 Å and 3729.88 Å (vacuum wavelengths), but
our spectral resolution is too low to resolve these. At a redshift of
z=1.006, other prominent emission lines such as H
4861,
and [O III]
4959, 5007 are not covered by our spectral
range. If the significant line were H
at z=0.139, we would
have expected to detect [OII] (at 4246 Å) and the [OIII] lines (at
5649 Å and 5705 Å) in the blue part of the spectrum, for which we
do not find any evidence. And finally, if the line were
[O III]
5007 at z=0.49, we would have expected to detect
[O III]
4959 (at 7406 Å) with roughly one-third of the
strength of [O III]
5007, and [O II] at 5566 Å,
even though this latter detection would have been made slightly
difficult due to the presence of a strong sky line at 5577 Å. The
[O II]
3727 at z=1.006 identification is strengthened by
the observed drop between the 8000-9000 Å continuum and that
blueward of the emission line. This is consistent with the red HST
V-I colour reported by Fruchter et al. (2002), and with the
possibility that this colour is caused by the 4000 Å break.
Assuming that the detected emission line indeed is [O II], we
can estimate the corresponding star-formation rate using the relation
yr-1 from Kennicutt (1998), where
is
in erg s-1. We measure an [OII] flux of
erg s-1 cm-2, where the
error is dominated by the uncertainty in the continuum level. As the
continuum level is close to zero, the measurement of the observed
equivalent width is highly uncertain:
=
90+110-40 Å. At the luminosity distance of GRB 021211 of
cm, the [O II] flux corresponds to a
star-formation rate of
M
yr-1,
which is a typical value of GRB host galaxies (ranging from 0.3
to 24 M
yr-1, e.g. see Table 1 of Vreeswijk et al. 2001a). We note
that this value has not been corrected for host-galaxy extinction, and
is therefore a strict lower limit to the actual star-formation rate.
We have determined the redshift for three GRBs: z=1.02 for GRB 991216,
z=2.142 for GRB 011211 and z=1.006 for GRB 021211. For GRB 991216, we also find evidence
for a foreground absorption system at z=0.80, and possibly another
at z=0.77. We argue that the flux depression from 4000-5500 Å in
the GRB 991216 spectrum is unlikely to be due to an error in the flux
calibration, and find that it could be explained by a host-galaxy
extinction bump similar to the Galactic 2175 Å feature, but at a
wavelength of 2360 Å. In the Galaxy, similar extinction bumps, with
a central wavelength between 2300 Å and 2500 Å, have been observed
toward UV-strong, hydrogen-poor stars. Fitting the Ly
absorption
line that we detect in the afterglow spectrum of GRB 011211, we obtain a
neutral hydrogen column density of log N(H I) =
,
indicating that it is a damped Ly
(DLA) system. This is at the low
end of the H I column densities detected in GRB
afterglows. Using a curve-of-growth analysis, we find the following
metallicity estimates at the redshift of GRB 011211:
[Si/H] =
-0.9+0.6-0.4,
,
and
[Al/H] =
-1.0+0.5-0.3. Assuming that the identification of
[O II]
3727 in the GRB 021211 spectrum is correct, we derive a
star-formation rate of
M
yr-1,
which is a typical value for GRB host galaxies.
Acknowledgements
P.M.V. thanks Thomas Szeifert for an enlightening discussion on the flux-calibration of FORS spectra. The rapid follow-up VLT spectroscopic observations discussed in this paper would not have been possible without the GRB Coordinates Network (GCN), set up by Scott Barthelmy at NASA's GSFC. We are grateful to Jochen Greiner and the University of Texas GRBlog team led by Robert Quimby for their extremely useful GRB webpages. We acknowledge benefits from collaboration within the Research Training Network "Gamma-Ray Bursts: An Enigma and a Tool'', funded by the EU under contract HPRN-CT-2002-00294. P.M.V. was partly supported by the NWO Spinoza grant 08-0 to E. P. J. van den Heuvel. L.K. has been supported by a fellowship of the Royal Academy of Arts & Sciences in the Netherlands. J.G. acknowledges the support of a Ramón y Cajal Fellowship from the Spanish Ministry of Education and Science. The Dark Cosmology Centre is supported by the DNRF.