Spectroscopy of the γ-ray burst GRB 021004: a structured jet ploughing through a massive stellar wind

We present spectra of the afterglow of the γ-ray burst GRB 021004 taken with the ISIS spectrograph on the William Herschel Telescope (WHT) and with the Focal Reducer/Low Dispersion Spectrograph 1 (FORS1) on the Very Large Telescope (VLT) at three epochs spanning 0.49–6.62 d after the burst. We observe strong absorption probably coming from the host galaxy, alongside absorption in H  i , Si  iv and C  iv with blueshifts of up to 2900 km s⁻¹ from the explosion centre, which we assume originates close to the progenitor. We find no significant variability of these spectral features. We investigate the origin of the outflowing material and evaluate various possible progenitor models. The most plausible explanation is that these result in the fossil stellar wind of a highly evolved Wolf–Rayet (WR) star. However, ionization from the burst itself prevents the existence of H  i , Si  iv and C  iv close to the afterglow surface where the fast stellar wind should dominate, and large amounts of blueshifted hydrogen are not expected in a WR star wind. We propose that the WR star wind is enriched by a hydrogen-rich companion, and that the GRB has a structured jet geometry in which the γ-rays emerge in a small opening angle within the wider opening angle of the cone of the afterglow. This scenario is able to explain both the spectral-line features and the irregular light curve of this afterglow.     

A unified treatment of the gamma-ray burst 021211 and its afterglow

The gamma-ray burst (GRB) 021211 had a simple light curve, containing only one peak and the expected Poisson fluctuations. Such a burst may be attributed to an external shock, offering the best chance for a unified understanding of the gamma-ray burst and afterglow emissions. We analyse the properties of the prompt (burst) and delayed (afterglow) emissions of GRB 021211 within the fireball model. Consistency between the optical emission during the first 11 min (which, presumably, comes from the reverse shock heating of the ejecta) and the later afterglow emission (arising from the forward shock) requires that, at the onset of deceleration (∼2 s), the energy density in the magnetic field in the ejecta, expressed as a fraction of the equipartition value  (ɛ _B )  , is larger than in the forward shock at 11 min by a factor of approximately 10³. We find that synchrotron radiation from the forward shock can account for the gamma-ray emission of GRB 021211; to explain the observed GRB peak flux requires that, at 2 s,  ɛ _B   in the forward shock is larger by a factor 100 than at 11 min. These results suggest that the magnetic field in the reverse shock and early forward shock is a frozen-in field originating in the explosion and that most of the energy in the explosion was initially stored in the magnetic field. We can rule out the possibility that the ejecta from the burst for GRB 021211 contained more than 10 electron–positron pairs per proton.     

Shock breakout in Type Ibc supernovae and application to GRB 060218/SN 2006aj

Recently, a soft blackbody component was observed in the early X-ray afterglow of GRB 060218, which was interpreted as shock breakout from the thick wind of the progenitor Wolf–Rayet (WR) star of the underlying Type Ic supernova 2006aj. In this paper, we present a simple model for computing the characteristic quantities (including energy, temperature and time duration) for the transient event from the shock breakout in Type Ibc supernovae produced by the core-collapse of WR stars surrounded by dense winds. In contrast to the case of a star without a strong wind, the shock breakout occurs in the wind region rather than inside the star, caused by the large optical depth in the wind. We find that, for the case of a WR star with a dense wind, the total energy of the radiation generated by the supernova shock breakout is larger than that in the case of the same star without a wind by a factor of >10. The temperature can be either hotter or colder, depending on the wind parameters. The time duration is larger caused by the increase in the effective radius of the star due to the presence of a thick wind. Then, we apply the model to GRB 060218/SN 2006aj. We show that, to explain both the temperature and the total energy of the blackbody component observed in GRB 060218 by the shock breakout, the progenitor WR star has to have an unrealistically large core radius (the radius at optical depth of 20), larger than 100 R_⊙. In spite of this disappointing result, our model is expected to have important applications to the observations on Type Ibc supernovae in which the detection of shock breakout will provide important clues to the progenitors of Type Ibc supernovae.     

Gamma-ray bursts from stellar remnants: probing the Universe at high redshift

A gamma-ray burst (GRB) releases an amount of energy similar to that of a supernova explosion, which combined with its rapid variability suggests an origin related to neutron stars or black holes. Since these compact stellar remnants form from the most massive stars not long after their birth, GRBs should trace the star formation rate in the Universe; we show that the GRB flux distribution is consistent with this. Because of the strong evolution of the star formation rate with redshift, it follows that the dimmest known bursts have z  ∼ 6, much above the value usually quoted and beyond the most distant quasars. This explains the absence of bright galaxies in well-studied GRB error boxes. The increased distances imply a peak luminosity of 8.3 × 10⁵¹ erg s⁻¹ and a rate density of 0.025 per million years per galaxy. These values are 20 times higher and 150 times lower, respectively, than are implied by fits with non-evolving GRB rates. This means either that GRBs are caused by a much rarer phenomenon than mergers of binary neutron stars, or that their gamma-ray emission is often invisible to us due to beaming. Precise burst locations from optical transients will discriminate between the various models for GRBs from stellar deaths, because the distance between progenitor birth place and burst varies greatly among them. The dimmest GRBs are then the most distant known objects, and may probe the Universe at an age when the first stars were forming.     

The circumstellar environment of Wolf–Rayet stars and gamma-ray burst afterglows

We study the evolution of the circumstellar medium of massive stars. We pay particular attention to Wolf-Rayet stars that are thought to be the progenitors of some long gamma-ray bursts (GRBs). We detail the mass-loss rates we use in our stellar evolution models and how we estimate the stellar wind speeds during different phases. With these details we simulate the interactions between the wind and the interstellar medium to predict the circumstellar environment around the stars at the time of core-collapse. We then investigate how the structure of the environment might affect the GRB afterglow. We find that when the afterglow jet encounters the free-wind/stalled-wind interface, rebrightening occurs and a bump is seen in the afterglow light curve. However, our predicted positions of this interface are too distant from the site of the GRB to reach while the afterglow remains observable. The values of the final wind density,   A _*  , from our stellar models are of the same order (≲1) as some of the values inferred from observed afterglow light curves. We do not reproduce the lowest   A _*  values below 0.5 inferred from afterglow observations. For these cases, we suggest that the progenitors could have been a WO-type Wolf–Rayet (WR) star or a very low-metallicity star. Finally, we turn our attention to the matter of stellar wind material producing absorption lines in the afterglow spectra. We discuss the observational signatures of two WR stellar types, WC and WO, in the afterglow light curve and spectra. We also indicate how it may be possible to constrain the initial mass and metallicity of a GRB progenitor by using the inferred wind density and wind velocity.     

A possible energy mechanism for cosmological γ-ray bursts

We suggest that an extreme Kerr black hole with a mass ∼10⁶ M_⊙, a dimensionless angular momentum     and a marginally stable orbital radius     located in a normal galaxy may produce a γ -ray burst (GRB) by capturing and disrupting a star. During the capture period, a transient accretion disc is formed and a strong transient magnetic field ∼     lasting for     may be produced at the inner boundary of the accretion disc. A large amount of rotational energy of the black hole is extracted and released in an ultrarelativistic jet with a bulk Lorentz factor Γ larger than 10³ via the Blandford–Znajek process. The relativistic jet energy can be converted into γ -radiation via an internal shock mechanism. The GRB duration should be the same as the lifetime of the strong transient magnetic field. The maximum number of sub-bursts is estimated to be     because the disc material is likely to break into pieces with a size about the thickness of the disc h at the cusp     The shortest risetime of the burst estimated from this model is ∼     The model GRB density rate is also estimated.     

The extreme, red afterglow of GRB 060923A: distance or dust?

Gamma-ray bursts (GRBs) are powerful probes of the early Universe, but locating and identifying very distant GRBs remain challenging. We report here the discovery of the K -band afterglow of Swift GRB 060923A, imaged within the first hour post-burst, and the faintest so far found. It was not detected in any bluer bands to deep limits, making it a candidate very high- z burst  ( z ≳ 11)  . However, our later-time optical imaging and spectroscopy reveal a faint galaxy coincident with the GRB position which, if it is the host, implies a more moderate redshift (most likely   z ≲ 2.8  ) and therefore that dust is the likely cause of the very red-afterglow colour. This being the case, it is one of the few instances so far found of a GRB afterglow with high-dust extinction.     

Massive stars exploding in a He-rich circumstellar medium – III. SN 2006jc: infrared echoes from new and old dust in the progenitor CSM

We present near- (NIR) and mid-infrared (MIR) photometric data of the Type Ibn supernova (SN) 2006jc obtained with the United Kingdom Infrared Telescope (UKIRT), the Gemini North Telescope and the Spitzer Space Telescope between days 86 and 493 post-explosion. We find that the IR behaviour of SN 2006jc can be explained as a combination of IR echoes from two manifestations of circumstellar material. The bulk of the NIR emission arises from an IR echo from newly condensed dust in a cool dense shell (CDS) produced by the interaction of the ejecta outward shock with a dense shell of circumstellar material ejected by the progenitor in a luminous blue variable (LBV)-like outburst about two years prior to the SN explosion. The CDS dust mass reaches a modest  3.0 × 10⁻⁴ M_⊙  by day 230. While dust condensation within a CDS formed behind the ejecta inward shock has been proposed before for one event (SN 1998S), SN 2006jc is the first one showing evidence for dust condensation in a CDS formed behind the ejecta outward shock in the circumstellar material. At later epochs, a substantial and growing contribution to the IR fluxes arises from an IR echo from pre-existing dust in the progenitor wind. The mass of the pre-existing circumstellar medium (CSM) dust is at least  ∼8 × 10⁻³ M_⊙  . This paper therefore adds to the evidence that mass-loss from the progenitors of core-collapse SNe could be a major source of dust in the Universe. However, yet again, we see no direct evidence that the explosion of an SN produces anything other than a very modest amount of dust.     

Gamma-ray burst afterglows: effects of radiative corrections and non-uniformity of the surrounding medium

The afterglow of a gamma-ray burst (GRB) is commonly thought to be the result of continuous deceleration of a relativistically expanding fireball in the surrounding medium. Assuming that the expansion of the fireball is adiabatic and that the density of the medium is a power-law function of shock radius, i.e. n _ext ∝  R ^{− k} , we study the effects of the first-order radiative correction and the non-uniformity of the medium on a GRB afterglow analytically. We first derive a new relation among the observed time, the shock radius and the Lorentz factor of the fireball: t _⊕ =  R /4(4− k ) γ²c, and also derive a new relation among the comoving time, the shock radius and the Lorentz factor of the fireball: t _co = 2 R /(5− k ) γc. We next study the evolution of the fireball by using the analytic solution of Blandford &38; McKee. The radiation losses may not significantly influence this evolution. We further derive new scaling laws both between the X-ray flux and observed time and between the optical flux and observed time. We use these scaling laws to discuss the afterglows of GRB 970228 and GRB 970616, and find that if the spectral index of the electron distribution is p  = 2.5, implied from the spectra of GRBs, the X-ray afterglow of GRB 970616 is well fitted by assuming k  = 2.     

Millimetric properties of gamma-ray burst host galaxies

We present millimetre (mm) and submillimetre (submm) photometry of a sample of five host galaxies of gamma-ray bursts (GRBs), obtained using the Max Planck Millimetre Bolometer (MAMBO2) array and Submillimetre Common-User Bolometer Array (SCUBA). These observations were obtained as part of an ongoing project to investigate the status of GRBs as indicators of star formation. Our targets include two of the most unusual GRB host galaxies, selected as likely candidate submm galaxies: the extremely red  ( R − K ≈ 5)  host of GRB 030115, and the extremely faint  ( R > 29.5)  host of GRB 020124. Neither of these galaxies is detected, but the deep upper limits for GRB 030115 impose constraints on its spectral energy distribution, requiring a warmer dust temperature than is commonly adopted for submillimetre galaxies (SMGs).
As a framework for interpreting these data, and for predicting the results of forthcoming submm surveys of Swift -derived host samples, we model the expected flux and redshift distributions based on luminosity functions of both submm galaxies and GRBs, assuming a direct proportionality between the GRB rate density and the global star formation rate density. We derive the effects of possible sources of uncertainty in these assumptions, including (1) introducing an anticorrelation between GRB rate and the global average metallicity, and (2) varying the dust temperature.     

Magnetized relativistic jets and long-duration GRBs from magnetar spin-down during core-collapse supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spin-down of a newly formed millisecond,   B ∼ 10¹⁵ G  , magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spin-down powers  (∼10⁵¹–10⁵² erg s⁻¹)  , the magnetar wind is superfast at almost all latitudes, while for lower spin-down powers  (∼10⁵⁰ erg s⁻¹)  , the wind is subfast but still super-Alfvénic. In all cases, the rates at which the neutron star loses mass, angular momentum and energy are very similar to the corresponding free wind values (≲30 per cent differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated  (∼5–10°)  relativistic jet out along the rotation axis of the star. Nearly all of the spin-down power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.     

A unified picture for the γ-ray and prompt optical emissions of GRB 990123

The prompt optical emission of GRB 990123 was uncorrelated to the γ-ray light curve and exhibited temporal properties similar to those of the steeply decaying, early X-ray emission observed by Swift at the end of many bursts. These facts suggest that the optical counterpart of GRB 990123 was the large-angle emission released during (the second pulse of) the burst. If the optical and γ-ray emissions of GRB 990123 have, indeed, the same origin then their properties require that (i) the optical counterpart was synchrotron emission and γ-rays arose from inverse-Compton scatterings (the 'synchrotron self-Compton model'), (ii) the peak energy of the optical-synchrotron component was at ∼20 eV and (iii) the burst emission was produced by a relativistic outflow moving at Lorentz factor  ≳450  and at a radius  ≳10¹⁵  cm, which is comparable to the outflow deceleration radius. Because the spectrum of GRB 990123 was optically thin above 2 keV, the magnetic field behind the shock must have decayed on a length-scale of  ≲1  per cent  of the thickness of the shocked gas, which corresponds to  10⁶–10⁷  plasma skin depths. Consistency of the optical counterpart decay rate and its spectral slope (or that of the burst, if they represent different spectral components) with the expectations for the large-angle burst emission represents the most direct test of the unifying picture proposed here for GRB 990123.     

Was GRB 990123 a unique optical flash?

GRB 990123 was a long, complex gamma-ray burst accompanied by an extremely bright optical flash. We find different constraints on the bulk Lorentz of this burst to be consistent with the speculation that the optical light is emission from the reverse shock component of the external shock. Motivated by this currently favoured idea, we compute the prompt reverse shock emission to be expected for bursts in which multiwavelength observations allow the physical parameters to be constrained. We find that for reasonable assumptions about the velocity of source expansion, a strong optical flash  m_V≈9  was expected from the reverse shocks, which were usually found to be mildly relativistic. The best observational prospects for detecting these prompt flashes are highlighted, along with the possible reasons for the absence of optical prompt detections in ongoing observations.     

Ly+ and ultraviolet emission from high-redshift gamma-ray burst hosts: to what extent do gamma-ray bursts trace star formation?

We report the result of a search for Lyα emission from the host galaxies of the gamma-ray bursts  (GRBs) 030226 ( z = 1.986), 021004 ( z = 2.335)  and  020124 ( z = 3.198)  . We find that the host galaxy of GRB 021004 is an extended (around 8 kpc) strong Lyα emitter with a rest-frame equivalent width (EW) of 68⁺¹²₋₁₁Å, and a star formation rate of  10.6 ± 2.0 M_⊙ yr⁻¹  . We do not detect the hosts of GRB 030226 and GRB 020124, but the upper limits on their Lyα fluxes do not rule out large rest-frame EWs. In the fields of GRB 021004 and GRB 030226 we find seven and five other galaxies, respectively, with excess emission in the narrow-band filter. These galaxies are candidate Lyα-emitting galaxies in the environment of the host galaxies. We have also compiled a list of all   z ≳ 2  GRB hosts, and demonstrate that a scenario where they trace star formation in an unbiased way is compatible with current observational constraints. Fitting the   z = 3  luminosity function (LF) under this assumption results in a characteristic luminosity of   R *= 24.6  and a faint-end slope of  α=−1.55  , consistent with the LF measured for Lyman-break galaxies.     

GRB 080319B: evidence for relativistic turbulence, not internal shocks

We show that the excellent optical and gamma-ray data available for GRB 080319B rule out the internal shock model for the prompt emission. The data instead point to a model in which the observed radiation was produced close to the deceleration radius  (∼10¹⁷ cm)  by a turbulent source with random Lorentz factors of ∼10 in the comoving frame. The optical radiation was produced by synchrotron emission from relativistic electrons, and the gamma-rays by inverse-Compton scattering of the synchrotron photons. The gamma-ray emission originated both in eddies and in an inter-eddy medium, whereas the optical radiation was mostly from the latter. Therefore, the gamma-ray emission was highly variable whereas the optical was much less variable. The model explains all the observed features in the prompt optical and gamma-ray data of GRB 080319B. We are unable to determine with confidence whether the energy of the explosion was carried outwards primarily by particles (kinetic energy) or magnetic fields. Consequently, we cannot tell whether the turbulent medium was located in the reverse shock (we can rule out the forward shock) or in a Poynting-dominated jet.     

A Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars

We introduce a new Rigid-Field Hydrodynamics approach to modelling the magnetospheres of massive stars in the limit of very strong magnetic fields. Treating the field lines as effectively rigid, we develop hydrodynamical equations describing the one-dimensional flow along each, subject to pressure, radiative, gravitational and centrifugal forces. We solve these equations numerically for a large ensemble of field lines to build up a three-dimensional time-dependent simulation of a model star with parameters similar to the archetypal Bp star σ Ori E. Since the flow along each field line can be solved independently of other field lines, the computational cost of this approach is a fraction of an equivalent magnetohydrodynamical treatment.
The simulations confirm many of the predictions of previous analytical and numerical studies. Collisions between wind streams from opposing magnetic hemispheres lead to strong shock heating. The post-shock plasma cools initially via X-ray emission, and eventually accumulates into a warped, rigidly rotating disc defined by the locus of minima of the effective (gravitational plus centrifugal) potential. However, a number of novel results also emerge. For field lines extending far from the star, the rapid area divergence enhances the radiative acceleration of the wind, resulting in high shock velocities (up to  ∼3000 km s⁻¹  ) and hard X-rays. Moreover, the release of centrifugal potential energy continues to heat the wind plasma after the shocks, up to temperatures around twice those achieved at the shocks themselves. Finally, in some circumstances the cool plasma in the accumulating disc can oscillate about its equilibrium position, possibly due to radiative cooling instabilities in the adjacent post-shock regions.     

Gamma-ray burst afterglow scaling coefficients for general density profiles

Gamma-ray burst (GRB) afterglows are well described by synchrotron emission originating from the interaction between a relativistic blast wave and the external medium surrounding the GRB progenitor. We introduce a code to reconstruct spectra and light curves from arbitrary fluid configurations, making it especially suited to study the effects of fluid flows beyond those that can be described using analytical approximations. As a check and first application of our code, we use it to fit the scaling coefficients of theoretical models of afterglow spectra. We extend earlier results of other authors to general circumburst density profiles. We rederive the physical parameters of GRB 970508 and compare with other authors.     

The Properties of Supernova 1997cy Associated with GRB 970514

The extraordinary SN 1997cy associated with GRB 970514 has been observed photometrically and spectroscopically for nearly 2 yr. At the time of discovery, SN 1997cy was the brightest supernova (SN) ever observed (MV相似文献   

Precursor activity in bright, long BATSE gamma-ray bursts

We study a sample of bright, long Burst and Transient Source Experiment (BATSE) gamma-ray burst (GRB) light curves in the 200 s before the detection of the GRB prompt emission. We find that in a sizable fraction of cases (∼20 per cent) there is evidence of emission above the background coming from the same direction as the GRB. This emission is characterized by a softer spectrum with respect to the main spectrum and contains a small fraction (0.1–1 per cent) of the total event counts. The precursors have typical delays of several tens of seconds extending (in few cases) up to 200 s (the limit of the investigated period). Their spectra are typically non-thermal power-law spectra, aside from a few cases. Such long delays and the non-thermal origin of their spectra are hard to reconcile with any model for the progenitor.