Gamma-ray Burst Blast Waves Research Articles

Ultrahigh-energy cosmic-ray signature in GRB 221009A

The brightest long gamma-ray burst (GRB) detected so far by the Swift-BAT and Fermi-GBM telescopes, GRB 221009A, provides an unprecedented opportunity for understanding the high-energy processes in extreme transient phenomena. We find that the conventional leptonic models for the afterglow emission from this source, synchrotron and synchrotron-self-Compton, have difficulties explaining the observation of ≳10 TeV γ rays (as high as 18 TeV) by the LHAASO detector. We modeled the γ-ray spectrum estimated in the energy range 0.1–1 GeV by the Fermi-LAT detector. The flux predicted by our leptonic models is severely attenuated at > 1 TeV due to γγ pair production with extragalactic background light, and hence an additional component is required at ≳10 TeV. Ultrahigh-energy cosmic rays can be accelerated in the GRB blast wave, and their propagation induces an electromagnetic cascade in the extragalactic medium. The line-of-sight component of this flux can explain the emission at ≳10 TeV detected by LHAASO, which requires a fraction of the GRB blast wave energy to be in ultrahigh-energy cosmic rays. This could be an indication of ultrahigh-energy cosmic-ray acceleration in GRBs.

Constraints on electron acceleration in gamma-ray bursts afterglows from radio peaks

ABSTRACT Studies of gamma-ray bursts (GRBs) and their multiwavelength afterglows have led to insights in electron acceleration and emission properties from relativistic, high-energy astrophysical sources. Broad-band modelling across the electromagnetic spectrum has been the primary means of investigating the physics behind these sources, although independent diagnostic tools have been developed to inform and corroborate assumptions made in particle acceleration simulations and broad-band studies. We present a methodology to constrain three physical parameters related to electron acceleration in GRB blast waves: the fraction of shock energy in electrons, ϵe; the fraction of electrons that gets accelerated into a power-law distribution of energies, ξe; and the minimum Lorentz factor of the accelerated electrons, γm. These parameters are constrained by observations of the peaks in radio afterglow light curves and spectral energy distributions. From a sample of 49 radio afterglows, we are able to find narrow distributions for these parameters, hinting at possible universality of the blast wave microphysics, although observational bias could play a role in this. Using radio peaks and considerations related to the prompt gamma-ray emission efficiency, we constrain the allowed parameter ranges for both ϵe and ξe to within about one order of magnitude, 0.01 ≲ ϵe ≲ 0.2 and 0.1 ≲ ξe ≲ 1. Such stringent constraints are inaccessible for ξe from broad-band studies due to model degeneracies.

Rapid-response radio observations of short GRB 181123B with the Australia Telescope Compact Array

ABSTRACT We introduce the Australia Telescope Compact Array (ATCA) rapid-response mode by presenting the first successful trigger on the short-duration gamma-ray burst (GRB) 181123B. Early-time radio observations of short GRBs may provide vital insights into the radio afterglow properties of Advanced LIGO- and Virgo-detected gravitational wave events, which will in turn inform follow-up strategies to search for counterparts within their large positional uncertainties. The ATCA was on target within 12.6 h post-burst, when the source had risen above the horizon. While no radio afterglow was detected during the 8.3 h observation, we obtained force-fitted flux densities of 7 ± 12 and $15 \pm 11\, \mu$Jy at 5.5 and 9 GHz, respectively. Afterglow modelling of GRB 181123B showed that the addition of the ATCA force-fitted radio flux densities to the Swift X-ray Telescope detections provided more stringent constraints on the fraction of thermal energy in the electrons (log $\epsilon _e = -0.75^{+0.39}_{-0.40}$ rather than log $\epsilon _e = -1.13^{+0.82}_{-1.2}$ derived without the inclusion of the ATCA values), which is consistent with the range of typical ϵe derived from GRB afterglow modelling. This allowed us to predict that the forward shock may have peaked in the radio band ∼10 d post-burst, producing detectable radio emission ≳3–4 d post-burst. Overall, we demonstrate the potential for extremely rapid radio follow-up of transients and the importance of triggered radio observations for constraining GRB blast wave properties, regardless of whether there is a detection, via the inclusion of force-fitted radio flux densities in afterglow modelling efforts.

The Second Plateau in X-Ray Afterglow Providing Additional Evidence for Rapidly Spinning Magnetars as the GRB Central Engine

Abstract Evidence for the central engine of gamma-ray bursts (GRBs) has been collected in the Neil Gehrels Swift data. For instance, some GRBs show an internal X-ray plateau followed by very steep decay, which is difficult interpret within the framework of a black hole (BH) central engine, but is consistent within a rapidly spinning magnetar engine picture. The very steep decay at the end of the plateau suggests a sudden cessation of the central engine, which is explained as the collapse of a supramassive magnetar into a BH when it spins down. Here we propose that some additional evidence, such as a second X-ray plateau feature, would show up if the fallback accretion could activate the newborn BH and sufficient energy could be transferred from the newborn BH to the GRB blast wave. With a systematic data analysis for all long GRBs, we find three candidates in the Swift sample, i.e., GRBs 070802, 090111, and 120213A, whose X-ray afterglow lightcurves contain two plateaus, with the first one being an internal plateau. We find that in a fairly loose and reasonable parameter space, the second X-ray plateau data for all 3 GRBs could be interpreted with our proposed model. Future observations are likely to discover similar events, which could offer more information on the properties of the magnetar, as well as the newborn BH.

Gamma-ray burst afterglow blast waves

The various stages of baryonic gamma-ray burst (GRB) afterglow blast waves are reviewed. These are responsible for the afterglow emission from which much of our understanding of gamma-ray bursts derives. Initially, the blast waves are confined to the dense medium surrounding the burster (stellar envelope or dense wind), giving rise to a jet-cocoon structure. A massive ejecta is released and potentially fed by ongoing energy release from the burster and a forward–reverse shock system is set up between ejecta and ambient density. Ultimately the blast wave spreads sideways and slows down, and the dominant afterglow emission shifts from X-rays down to radio. Over the past years significant progress has been made both observationally and theoretically/numerically in our understanding of these blast waves, unique in the universe due to their often incredibly high initial Lorentz factors of 100–1000. The recent discovery of a short gamma-ray burst counterpart to a gravitational wave detection (GW 170817) brings the promise of a completely new avenue to explore and constrain the dynamics of gamma-ray burst blast waves.

Very High Energy Neutrinos from nearby long GRB Afterglows

Long duration Gamma Ray Bursts (GRBs) are well-motivated sources of Ultra High Energy Cosmic Rays (UHECRs) and neutrinos. During the afterglow phase these particles can be produced as a result of acceleration and interaction there in. We have modeled afterglow spectra and light curves from synchrotron cooling of accelerated electrons. We have fitted data of 17 long GRBs detected within redshift 0.5 in case of the GRB blastwave evolving in a wind and constant density interstellar medium. The afterglow photons can interact with the shock accelerated protons to produce very high energy neutrinos. We have calculated the neutrino flux for photo-pion interactions for all these GRBs. As IceCube have been detecting very high energy neutrinos for the last four years and a larger future extension called Gen 2 is planned, this calculation will help in understanding more about GRB neutrino production. Calculation of flux and estimation of events for Northern Hemisphere GRBs are done for the upcoming neutrino observatory KM3NeT.

Electron and ion acceleration in relativistic shocks with applications to GRB afterglows

We have modeled the simultaneous first-order Fermi shock acceleration of protons, electrons, and helium nuclei by relativistic shocks. By parameterizing the particle diffusion, our steady-state Monte Carlo simulation allows us to follow particles from particle injection at nonthermal thermal energies to above PeV energies, including the nonlinear smoothing of the shock structure due to cosmic-ray (CR) backpressure. We observe the mass-to-charge (A/Z) enhancement effect believed to occur in efficient Fermi acceleration in non-relativistic shocks and we parameterize the transfer of ion energy to electrons seen in particle-in-cell (PIC) simulations. For a given set of environmental and model parameters, the Monte Carlo simulation determines the absolute normalization of the particle distributions and the resulting synchrotron, inverse-Compton, and pion-decay emission in a largely self-consistent manner. The simulation is flexible and can be readily used with a wide range of parameters typical of gamma-ray burst (GRB) afterglows. We describe some preliminary results for photon emission from shocks of different Lorentz factors and outline how the Monte Carlo simulation can be generalized and coupled to hydrodynamic simulations of GRB blast waves. We assume Bohm diffusion for simplicity but emphasize that the nonlinear effects we describe stem mainly from an extended shock precursor where higher energy particles diffuse further upstream. Quantitative differences will occur with different diffusion models, particularly for the maximum CR energy and photon emission, but these nonlinear effects should be qualitatively similar as long as the scattering mean free path is an increasing function of momentum.

PeV-EeV neutrinos from GRB blast waves in IceCube and future neutrino telescopes

Ultrahigh-energy cosmic rays (UHECRs), if accelerated in the gamma-ray burst (GRB) blast wave, are expected to produce PeV-EeV neutrinos by interacting with long-lived GRB afterglow photons. Detailed spectral and temporal properties of the flux of these neutrinos depend on the GRB blast wave evolution scenario, but can last for days to years time scale in contrast to the seconds to minutes time scale for ``burst'' neutrino flux contemporaneous with the prompt gamma-ray emission and which has been constrained by IceCube in the $\ensuremath{\sim}50\text{ }\text{ }\mathrm{TeV}\ensuremath{-}2\text{ }\text{ }\mathrm{PeV}$ range. We compute expected neutrino events in IceCube in the PeV-EeV range from the blast wave of long-duration GRBs, both for the diffuse flux and for individual GRBs in the nearby universe. We show that IceCube will be able to detect the diffuse GRB blast wave neutrino flux after 5 years in operation, and will be able to distinguish it from the cosmogenic neutrino flux arising from the Greisen-Zatsepin-Kuzmin process in case the ultrahigh-energy cosmic rays are heavy nuclei. We also show that EeV neutrinos from the blast wave of an individual GRB can be detected with long-term monitoring by a future high-energy extension of IceCube for redshift up to $z\ensuremath{\sim}0.5$.

Detailed afterglow modelling and host galaxy properties of the dark GRB 111215A

Gamma-ray burst (GRB) 111215A was bright at X-ray and radio frequencies, but not detected in the optical or near-infrared (nIR) down to deep limits. We have observed the GRB afterglow with the Westerbork Synthesis Radio Telescope and Arcminute Microkelvin Imager at radio frequencies, with the William Herschel Telescope and Nordic Optical Telescope in the nIR/optical, and with the Chandra X-ray Observatory. We have combined our data with the Swift X-Ray Telescope monitoring, and radio and millimeter observations from the literature to perform broadband modeling, and determined the macro- and microphysical parameters of the GRB blast wave. By combining the broadband modeling results with our nIR upper limits we have put constraints on the extinction in the host galaxy. This is consistent with the optical extinction we have derived from the excess X-ray absorption, and higher than in other dark bursts for which similar modeling work has been performed. We also present deep imaging of the host galaxy with the Keck I telescope, Spitzer Space Telescope, and Hubble Space Telescope (HST), which resulted in a well-constrained photometric redshift, giving credence to the tentative spectroscopic redshift we obtained with the Keck II telescope, and estimates for the stellar mass and star formation rate of the host. Finally, our high resolution HST images of the host galaxy show that the GRB afterglow position is offset from the brightest regions of the host galaxy, in contrast to studies of optically bright GRBs.

ON THE ORIGIN OF GeV EMISSION IN GAMMA-RAY BURSTS

The most common progenitors of gamma-ray bursts (GRBs) are massive stars with strong stellar winds. We show that the GRB blast wave in the wind should emit a bright GeV flash. It is produced by inverse Compton cooling of the thermal plasma behind the forward shock wave. The main part of the flash is shaped by scattering of the prompt MeV radiation (emitted at smaller radii) which streams through the external blast wave. The inverse-Compton flash is bright due to the huge e+- enrichment of the external medium. At late times, the blast wave switches to normal synchrotron-self-Compton cooling. The mechanism is demonstrated by a detailed transfer simulation. The observed prompt MeV radiation is taken as an input of the simulation; we use GRB 080916C as an example. The result reproduces the GeV flash observed by the Fermi telescope. It explains the delayed onset, the steep rise, the peak flux, the time of the peak, the long smooth decline, and the spectral slope of GeV emission. The wind density required to reproduce all these features is typical of Wolf-Rayet stars. Our simulation predicts strong TeV emission 1 min after the burst trigger; then a cutoff in the observed high-energy spectrum is expected from absorption by extragalactic background light. In addition, a bright optical counterpart of the GeV flash is predicted for plausible values of the magnetic field; such a double (optical+GeV) flash has been observed in GRB 130427A.

The plateau phase of gamma-ray burst afterglows in the thick-shell scenario

We present analytic calculations of synchrotron radiation from the forward and the reverse shock of gamma-ray burst blast waves, in the thick-shell scenario (i.e. when the reverse shock is relativistic). We show that this scenario can naturally account for the plateau phase, observed early in the afterglows of about half the bursts detected by Swift. We generalize our approach to include power-law luminosity of the central engine and show that when radiation from both regions (forward and reverse shock) is taken into account, a wide range of possibilities emerge, including chromatic and achromatic breaks, frequency-dependent spectral evolution during the injection break and widely varying decay indices in different bands. For both the forward and the reverse shock, we derive formulas for the spectral parameters and the observed flux in different power-law segments of the spectrum, as a function of observer time. We explore the Fb-tb relation (between the observed time of the end of the plateau phase and the flux at that point) in the framework of the presented model and show that model predictions favour the reverse shock as the dominant source of emission in both optical and X-rays. As case studies, we present simultaneous fits to X-ray and optical/IR afterglow data of GRB 080928 and GRB 090423. We identify the end of the plateau phase with the cessation of energy injection and infer the corresponding upper limits to central-engine activity, which are about 1 h for the former and 1.5 h for the latter. We conclude that smooth energy injection through the reverse shock is a plausible explanation for the plateau phase of gamma-ray burst afterglows. During that phase, radiation from the reverse shock is likely to be important, or even dominant, and should be taken into account when fitting model parameters to observations.

Long-lived PeV–EeV neutrinos from gamma-ray burst blastwave

Long duration gamma-ray bursts are powerful sources that can accelerate particles to ultrahigh energies. Acceleration of protons in the forward shock of the highly relativistic gamma-ray burst (GRB) blastwave allows PeV-EeV neutrino production by photopion interactions of ultrahigh energy protons with x-ray to optical photons of the GRB afterglow emission. Four different blastwave evolution scenarios are considered: adiabatic and fully radiative blastwaves in a constant density circumburst medium and in a wind environment with the particle density in the wind decreasing inversely proportional to the square of the radius from the center of the burst. The duration of the neutrino flux depends on the evolution of the blastwave and can last up to a day in the case of an adiabatic blastwave in a constant density medium. Neutrino fluxes from the three other blastwave evolution scenarios are also calculated. Diffuse neutrino fluxes calculated using the observed rate of long-duration GRBs are consistent with the recent IceCube upper limit on the prompt GRB neutrino flux below PeV. The diffuse neutrino flux needed to explain the two neutrino events at PeV energies recently detected by IceCube can partially come from the presented GRB blastwave diffuse fluxes. Future observations by IceCube and upcoming huge radio Askaryan experiments will be able to probe the flux models presented here or constrain the GRB blastwave properties.

OBSERVATIONS OF PARTICLE ACCELERATION IN THE BLAST WAVES OF GAMMA-RAY BURSTS

The electron energy distribution index, p, is a fundamental parameter of the process by which electrons at shock fronts are accelerated to relativistic speeds and by which they radiate, via synchrotron emission. This acceleration process is applicable to a myriad of astronomical jet sources such as AGN, X-ray binaries and gamma-ray bursts (GRBs), as well as to particle acceleration in the solar wind and supernovae, and the acceleration of cosmic rays. The accurate measurement of the distribution of p is of fundamental importance to differentiate between the possible theories of electron acceleration at any relativistic shock front; there is division as to whether p has a universal value or whether it has a distribution, and if so, what that distribution is. Here one such source of synchrotron emission is examined: the blast waves of GRB afterglows observed by the Swift satellite. Within the framework of the GRB blast wave model, the constraints placed on the distribution of p by the observed X-ray spectral and temporal indices are examined and the distribution parametrized. The results show that the observed distribution of spectral indices is inconsistent with an underlying distribution of p composed of a single discrete value but consistent with a Gaussian distribution centred at p = 2.4 and having a width of 0.6. This finding disagrees with theoretical work that argues for a single, universal value of p, but also demonstrates that the width of the distribution is not as wide as has been suggested by some authors.

Gamma-ray burst afterglow light curves from realistic density profiles

The afterglow emission that follows gamma-ray bursts (GRBs) contains valuable information about the circumburst medium and, therefore, about the GRB progenitor. Theoretical studies of GRB blast waves, however, are often limited to simple density profiles for the external medium (mostly constant density and power-law R^{-k} ones). We argue that a large fraction of long-duration GRBs should take place in massive stellar clusters where the circumburst medium is much more complicated. As a case study, we simulate the propagation of a GRB blast wave in a medium shaped by the collision of the winds of O and Wolf-Rayet stars, the typical distance of which is d /sim 0.1 - 1 pc. Assuming a spherical blast wave, the afterglow light curve shows a flattening followed by a shallow break on a timescale from hours up to a week after the burst, which is a result of the propagation of the blast wave through the shocked wind region. If the blast wave is collimated, the jet break may, in some cases, become very pronounced with the post break decline of the light curve as steep as t-5. Inverse Compton scattering of ultra-violet photons from the nearby star off energetic electrons in the blast wave leads to a bright \simGeV afterglow flare that may be detectable by Fermi.

A SEMI-ANALYTIC FORMULATION FOR RELATIVISTIC BLAST WAVES WITH A LONG-LIVED REVERSE SHOCK

This paper performs a semi-analytic study of relativistic blast waves in the context of gamma-ray bursts (GRBs). Although commonly used in a wide range of analytical and numerical studies, the equation of state (EOS) with a constant adiabatic index is a poor approximation for relativistic hydrodynamics. Adopting a more realistic EOS with a variable adiabatic index, we present a simple form of jump conditions for relativistic hydrodynamical shocks. Then we describe in detail our technique of modeling a very general class of GRB blast waves with a long-lived reverse shock. Our technique admits an arbitrary radial stratification of the ejecta and ambient medium. We use two different methods to find dynamics of the blast wave: (1) customary pressure balance across the blast wave and (2) the "mechanical model". Using a simple example model, we demonstrate that the two methods yield significantly different dynamical evolutions of the blast wave. We show that the pressure balance does not satisfy the energy conservation for an adiabatic blast wave while the mechanical model does. We also compare two sets of afterglow light curves obtained with the two different methods.

Neutrino emission from a gamma-ray burst afterglow shock during an inner supernova shock breakout

The observations of a nearby low-luminosity gamma-ray burst (GRB) 060218 associated with supernova SN 2006aj may imply an interesting astronomical picture where a supernova shock breakout locates behind a relativistic GRB jet. Based on this picture, we study neutrino emission for early afterglows of GRB 060218-like GRBs, where neutrinos are expected to be produced from photopion interactions in a GRB blast wave that propagates into a dense wind. Relativistic protons for the interactions are accelerated by an external shock, while target photons are basically provided by the incoming thermal emission from the shock breakout and its inverse-Compton scattered component. Because of a high estimated event rate of low-luminosity GRBs, we would have more opportunities to detect afterglow neutrinos from a single nearby GRB event of this type by IceCube. Such a possible detection could provide evidence for the picture described above.

GRB blastwaves through wind-shaped circumburst media

Context:A significant fraction of progenitors for long gamma-ray bursts (GRBs) are believed to be massive stars. The investigation of long GRBs therefore requires modeling the propagation of ultra-relativistic blastwaves through the circumburst medium surrounding massive stars. We simulate the expansion of an isotropic, adiabatic relativistic fireball into the wind-shaped medium around a massive GRB progenitor. The circumburst medium is composed of a realistically stratified stellar wind zone up to its termination shock, followed by a region of shocked wind characterized by a constant density. Aims: We followed the evolution of the blastwave through all its stages, including the extremely rapid acceleration up to a Lorentz factor 75 flow, its deceleration by interaction with stellar wind, its passage of the wind termination shock, until its propagation through shocked wind. Methods: We used the adaptive mesh refinement versatile advection code to follow the evolution of the fireball. Results: We show that the acceleration from purely thermal to ultra-relativistic kinetic regimes is abrupt and produces an internally structured blastwave. We resolved the structure of this ultra-relativistic shell in all stages, thanks to the adaptive mesh. We comment on the dynamical roles played by forward and reverse shock pairs in the phase of interaction with the free stellar wind and clearly identify the complex shock-dominated structure created when the shell crosses the terminal shock. Conclusion: We show that in our model where the terminal shock is taken relatively close to the massive star, the phase of self-similar deceleration of Blandford-McKee type can only be produced in the constant density, shocked wind zone.

GRB 030329: 3 years of radio afterglow monitoring

Radio observations of gamma-ray burst (GRB) afterglows are essential for our understanding of the physics of relativistic blast waves, as they enable us to follow the evolution of GRB explosions much longer than the afterglows in any other wave band. We have performed a 3-year monitoring campaign of GRB 030329 with the Westerbork Synthesis Radio Telescopes and the Giant Metrewave Radio Telescope. Our observations, combined with observations at other wavelengths, have allowed us to determine the GRB blast wave physical parameters, such as the total burst energy and the ambient medium density, as well as to investigate the jet nature of the relativistic outflow. Further, by modelling the late-time radio light curve of GRB 030329, we predict that the Low-Frequency Array (30-240 MHz) will be able to observe afterglows of similar GRBs, and constrain the physics of the blast wave during its non-relativistic phase.

Ultra-high energy cosmic rays, cascade gamma rays, and high-energy neutrinos from gamma-ray bursts

Gamma-ray bursts (GRBs) are sources of energetic, highly variable fluxes of γ rays, which demonstrates that they are powerful particle accelerators. Besides relativistic electrons, GRBs should also accelerate high-energy hadrons, some of which could escape cooling to produce ultra-high energy cosmic rays (UHECRs). Acceleration of high-energy hadrons in GRB blast waves will be established if high-energy neutrinos produced through photopion interactions in the blast wave are detected from GRBs. Limitations on the energy in non-thermal hadrons and the number of expected neutrinos are imposed by the fluxes from pair-photon cascades initiated in the same processes that produce neutrinos. Only the most powerful bursts at fluence levels ≳ 3× 10−4 erg cm−2 offer a realistic prospect for detection of ≫TeV neutrinos. Detection of high-energy neutrinos is likely if GRB blast waves have large baryon loads and Doppler factors ≲ 200. Cascade γ rays will accompany neutrino production and might already have been detected as anomalous emission components in the spectra of some GRBs. Prospects for detection of GRBs in the Milky Way are also considered.

Acceleration mechanisms 1: shock acceleration

Acceleration by shock waves is one of the favorite models envisaged to explain the origin of the highest energy cosmic rays. This view was developed in the framework of various astrophysical acceleration sites: hot spots of radio galaxies, clusters of galaxies, jets powered by supermassive black holes, gamma ray burst blast waves, etc. These mechanisms along with their limitations and specific signatures for each type are reviewed. To cite this article: M. Ostrowski, C. R. Physique 5 (2004).