Abstract Blazar TXS 0506+056 is the main candidate for a coincident neutrino and gamma-ray flare event. In this paper, we compose a detailed kinetic lepto-hadronic emission model capable of producing a photon and neutrino spectrum given a set of parameters. Our model includes a range of large-scale geometries and both dynamical and steady-state injection models for electrons and protons. We link this model with a Markov Chain Monte Carlo sampler to obtain a powerful statistical tool that allows us to both fit the Spectral Energy Distribution and study the probability density functions and correlations of the parameters. Assuming a fiducial neutrino flux, we demonstrate how multi-messenger observations can be modelled jointly in a Bayesian framework. We find the best parameters for each of the variants of the model tested and report on their cross-correlations. Additionally, we confirm that reproducing the neutrino flux of TXS 0506+056 requires an extreme proton to electron ratio either in the local acceleration process or from external injection.

Recent observations with the Chandra X-ray telescope continue to detect X-ray emission from the transient GW170817. In a total exposure of 96.6 ks, performed between March 9 and March 16 2020 (935 d to 942 d after the merger), a total of 8 photons are measured at the source position, corresponding to a significance of about 5 sigma. Radio monitoring with the Australian Telescope Compact Array (ATCA) shows instead that the source has faded below our detection threshold (<33 uJy, 3 sigma). By assuming a constant spectral index beta=0.585, we derive an unabsorbed X-ray flux of approximately 1.4E-15 erg/cm^2/s, higher than earlier predictions, yet still consistent with a simple structured jet model. We discuss possible scenarios that could account for prolonged emission in X-rays. The current dataset appears consistent both with energy injection by a long-lived central engine and with the onset of a kilonova afterglow, arising from the interaction of the sub-relativistic merger ejecta with the surrounding medium. Long-term monitoring of this source will be essential to test these different models.

Late activity of the central engine is often invoked in order to explain the flares observed in the early X-ray afterglow of gamma-ray bursts, either in the form of an active neutron star remnant or (fall-back) accretion onto a black hole. However, these scenarios are not always plausible, in particular when flares are delayed to very late times after the burst. Recently, a new scenario was proposed that suggests X-ray flares can be the result of the passing of a long-lived reverse shock through a stratified ejecta, with the advantage that it does not require late-time engine activity. In this work, we numerically demonstrate this scenario to be physically plausible, by performing onedimensional simulations of ejecta dynamics and emission using our novel moving-mesh relativistic hydrodynamics code. Improved efficiency and precision over previous work enables the exploration of a broader range of setups. We can introduce a more physically realistic description of the circumburst medium mass density. We can also locally trace the cooling of electrons when computing the broadband emission from these setups. We show that the synchrotron cooling timescale can dominate the flare decay time if the stratification in the ejecta is constrained to a localised angular region inside the jet, with size corresponding to the relativistic causal connection angle, and that it corresponds to values reported in observations. We demonstrate that this scenario can produce a large range of observed flare times, suggesting a connection between flares and initial ejection dynamics rather than with late-time remnant activity.

The afterglow emission from gamma-ray bursts (GRBs) is a valuable source of information to understand the physics of these energetic explosions. The blast wave model has become the standard to describe the evolution of the afterglow emission over time and frequency. Thanks to recent developments in the theory of afterglows and numerical simulations of relativistic outflows, we are able to model the afterglow emission with realistic dynamics and radiative processes. Although the models agree with observations remarkably well, the afterglow emission still contains additional physics, instrumental systematics, and propagation effects which make the modelling of these events challenging. In this work, we present a new approach to modelling GRB afterglows, using Gaussian processes (GPs) to take into account systematics in the afterglow data. We show that, using this new approach, it is possible to obtain more reliable estimates of the explosion and microphysical parameters of GRBs. We present fit results for 5 long GRBs and find a preliminary correlation between the isotropic energetics and opening angles of GRBs, which confirms the idea of a common energy reservoir for the kinetic energy of long GRBs.

We present comprehensive multiwavelength radio to X-ray observations of GRB 181201A spanning from $\approx150$ s to $\approx163$ days after the burst, comprising the first joint ALMA-VLA-GMRT observations of a gamma-ray burst (GRB) afterglow. The radio and mm-band data reveal a distinct signature at $\approx3.9$ days, which we interpret as reverse shock (RS) emission. Our observations present the first time that a single radio-frequency spectral energy distribution can be decomposed directly into RS and forward shock (FS) components. We perform detailed modeling of the full multiwavelength data set, using Markov Chain Monte Carlo sampling to construct the joint posterior density function of the underlying physical parameters describing the RS and FS synchrotron emission. We uncover and account for all degeneracies in the model parameters. The joint RS-FS modeling reveals a weakly magnetized ($\sigma\approx3\times10^{-3}$), mildly relativistic RS, from which we derive an initial bulk Lorentz factor of $\Gamma_0\approx103$ for the GRB jet. Our results support the hypothesis that low-density environments are conducive to the observability of RS emission. We compare our observations to other events with strong RS detections, and find a likely observational bias selecting for longer lasting, non-relativistic reverse shocks. We present and begin to address new challenges in modeling posed by the present generation of comprehensive, multi-frequency data sets.

We present ALMA 97.5 GHz total intensity and linear polarization observations of the mm-band afterglow of GRB 190114C spanning 2.2 to 5.2 hours after the burst. We detect linear polarization at the $\approx 5\,\sigma$ level, decreasing from $\Pi=(0.87\pm0.13)\%$ to $(0.60\pm0.19)\%$, and evolving in polarization position angle from $(10\pm5)^\circ$ to $(-44\pm12)^\circ$ during the course of the observations. This represents the first detection of polarized millimeter emission in a $\gamma$-ray burst. We show that the optical and X-ray observations between $0.03$ days and $\sim0.3$ days are consistent with a fast cooling forward shock expanding into a wind environment. However, the optical observations at $\lesssim0.03$ days, as well as the radio and millimeter observations arise from a separate component, which we interpret as emission from the reverse-shocked ejecta. Using the measured linear polarization, we constrain the coherence scale of tangled magnetic fields in the ejecta to an angular size of $\theta_{\rm B} \approx10^{-3}$ radian, while the rotation of the polarization angle rules out the presence of large scale, ordered axisymmetric magnetic fields, and in particular a large scale toroidal field, in the jet.

Gamma-ray bursts (GRBs) associated with gravitational wave events are, and will likely continue to be, viewed at a larger inclination than GRBs without gravitational wave detections. As demonstrated by the afterglow of GW170817, this requires an extension of the common GRB afterglow models which typically assume emission from an on-axis top hat jet. We present a characterization of the afterglows arising from structured jets, providing a framework covering both successful and choked jets. We compute new closure relations for structured jets and compare with the established relations for energy injection and refreshed shock models. The temporal slope before the jet break is found to be a simple function of the ratio between the viewing angle and effective opening angle of the jet. A numerical model to calculate synthetic light curves and spectra is publicly available as the open source Python package afterglowpy.

We present the results of our year-long afterglow monitoring of GW170817, the first binary neutron star (NS) merger detected by advanced LIGO and advanced Virgo. New observations with the Australian Telescope Compact Array (ATCA) and the Chandra X-ray Telescope were used to constrain its late-time behavior. The broadband emission, from radio to X-rays, is well-described by a simple power-law spectrum with index ~0.585 at all epochs. After an initial shallow rise ~t^0.9, the afterglow displayed a smooth turn-over, reaching a peak X-ray luminosity of ~5e39 erg/s at 160 d, and has now entered a phase of rapid decline ~t^(-2). The latest temporal trend challenges most models of choked jet/cocoon systems, and is instead consistent with the emergence of a relativistic structured jet seen at an angle of ~22 deg from its axis. Within such model, the properties of the explosion (such as its blastwave energy E_K~2E50 erg, jet width theta_c~4 deg, and ambient density n~3E-3 cm^(-3)) fit well within the range of properties of cosmological short GRBs.

We investigate the prompt emission and the afterglow properties of short duration gamma-ray burst (sGRB) 130603B and another eight sGRB events during 2012-2015, observed by several multi-wavelength facilities including the GTC 10.4m telescope. Prompt emission high energy data of the events were obtained by INTEGRAL/SPI/ACS, Swift/BAT and Fermi/GBM satellites. The prompt emission data by INTEGRAL in the energy range of 0.1-10 MeV for sGRB 130603B, sGRB 140606A, sGRB 140930B, sGRB 141212A and sGRB 151228A do not show any signature of the extended emission or precursor activity and their spectral and temporal properties are similar to those seen in case of other short bursts. For sGRB130603B, our new afterglow photometric data constraints the pre jet-break temporal decay due to denser temporal coverage. For sGRB 130603B, the afterglow light curve, containing both our new as well as previously published photometric data is broadly consistent with the ISM afterglow model. Modeling of the host galaxies of sGRB 130603B and sGRB 141212A using the LePHARE software supports a scenario in which the environment of the burst is undergoing moderate star formation activity. From the inclusion of our late-time data for 8 other sGRBs we are able to; place tight constraints on the non-detection of the afterglow, host galaxy or any underlying kilonova emission. Our late-time afterglow observations of the sGRB 170817A/GW170817 are also discussed and compared with the sub-set of sGRBs.

We investigate the afterglow of GRB 140713A, a gamma-ray burst (GRB) that was detected and relatively well-sampled at X-ray and radio wavelengths, but was not present at optical and near-infrared wavelengths, despite searches to deep limits. We present the emission spectrum of the likely host galaxy at $z = 0.935$ ruling out a high-redshift explanation for the absence of the optical flux detection. Modelling the GRB multi-wavelength afterglow using the radiative transfer hydrodynamics code BOXFIT provides constraints on physical parameters of the GRB jet and its environment, for instance a relatively wide jet opening angle and an electron energy distribution slope $p$ below 2. Most importantly, the model predicts an optical flux about two orders of magnitude above the observed limits. We calculated that the required host extinction to explain the observed limits in the $r$, $i$ and $z$ bands was $A^{\rm host}_{V} > 3.2$ mag, equivalent to $E(B-V)^{\rm host} > 1.0$ mag. From the X-ray absorption we derive that the GRB host extinction is $A^{\rm host}_{\rm V} = 11.6^{+7.5}_{-5.3}$ mag, equivalent to $E(B-V)^{\rm host} = 3.7^{+2.4}_{-1.7}$ mag, which is consistent with the extinction required from our BOXFIT derived fluxes. We conclude that the origin of the optical darkness is a high level of extinction in the line of sight to the GRB, most likely within the GRB host galaxy.

Multi-messenger observations of GW170817 have not conclusively established whether the merger remnant is a black hole (BH) or a neutron star (NS). We show that a long-lived magnetized NS with a poloidal field $B\approx 10^{12}$G is fully consistent with the electromagnetic dataset, when spin down losses are dominated by gravitational wave (GW) emission. The required ellipticity $\epsilon\gtrsim 10^{-5}$ can result from a toroidal magnetic field component much stronger than the poloidal component, a configuration expected from a NS newly formed from a merger. Abrupt magnetic dissipation of the toroidal component can lead to the appearance of X-ray flares, analogous to the one observed in gamma-ray burst (GRB) afterglows. In the X-ray afterglow of GW170817 we identify a low-significance ($\gtrsim 3\sigma$) temporal feature at 155 d, consistent with a sudden reactivation of the central NS. Energy injection from the NS spin down into the relativistic shock is negligible, and the underlying continuum is fully accounted for by a structured jet seen off-axis. Whereas radio and optical observations probe the interaction of this jet with the surrounding medium, observations at X-ray wavelengths, performed with adequate sampling, open a privileged window on to the merger remnant.

The recent discovery of a faint gamma-ray burst (GRB) coincident with the gravitational wave (GW) event GW 170817 revealed the existence of a population of low-luminosity short duration gamma-ray transients produced by neutron star mergers in the nearby Universe. These events could be routinely detected by existing gamma-ray monitors, yet previous observations failed to identify them without the aid of GW triggers. Here we show that GRB150101B was an analogue of GRB170817A located at a cosmological distance. GRB 150101B was a faint short duration GRB characterized by a bright optical counterpart and a long-lived X-ray afterglow. These properties are unusual for standard short GRBs and are instead consistent with an explosion viewed off-axis: the optical light is produced by a luminous kilonova component, while the observed X-rays trace the GRB afterglow viewed at an angle of ~13 degrees. Our findings suggest that these properties could be common among future electromagnetic counterparts of GW sources.

GW170817 was not merely an absolute breakthrough in gravitational wave astrophysics and a first in multi-messenger astronomy. The quality and diversity of the electro-magnetic counterpart emission is staggering on its own as well, including unprecedented kilonova spectra and a broadband off-axis gamma-ray burst afterglow that has progressed along a trajectory of rise and decay and by now has even been measured using very large baseline interferometry. For these proceedings, I will summarize the points for discussion that I presented during the workshop regarding off-axis short gamma-ray bursts and their (un-)successful jets and their emission. Given that developments are currently moving very fast in the field, I also touch on some results that have appeared in the literature following the Vulcano meeting.

We present our broadband study of GW170817 from radio to hard X-rays, including NuSTAR and Chandra observations up to 165 days after the merger, and a multi-messenger analysis including LIGO constraints. The data are compared with predictions from a wide range of models, providing the first detailed comparison between non-trivial cocoon and jet models. Homogeneous and power-law shaped jets, as well as simple cocoon models are ruled out by the data, while both a Gaussian shaped jet and a cocoon with energy injection can describe the current dataset for a reasonable range of physical parameters, consistent with the typical values derived from short GRB afterglows. We propose that these models can be unambiguously discriminated by future observations measuring the post-peak behaviour, with slope -1.0 for the cocoon and -2.5 for the jet model.

Deriving physical parameters from gamma-ray burst afterglow observations remains a challenge, even now, 20 years after the discovery of afterglows. The main reason for the lack of progress is that the peak of the synchrotron emission is in the sub-mm range, thus requiring radio observations in conjunction with X-ray/optical/near-infrared data in order to measure the corresponding spectral slopes and consequently remove the ambiguity wrt. slow vs. fast cooling and the ordering of the characteristic frequencies. We observed GRB 151027B, the 1000th Swift-detected GRB, with GROND in the optical-NIR, ALMA in the sub-millimeter, ATCA in the radio band, and combine this with public Swift-XRT X-ray data. While some observations at crucial times only return upper limits or surprising features, the fireball model is narrowly constrained by our data set, and allows us to draw a consistent picture with a fully-determined parameter set. Surprisingly, we find rapid, large-amplitude flux density variations in the radio band which are extreme not only for GRBs, but generally for any radio source. We interpret these as scintillation effects, though the extreme nature requires either the scattering screen to be at much smaller distance than usually assumed, multiple screens, or a combination of the two.

The various stages of baryonic gamma-ray burst 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.

A long-standing paradigm in astrophysics is that collisions- or mergers- of two neutron stars (NSs) form highly relativistic and collimated outflows (jets) powering gamma-ray bursts (GRBs) of short (< 2 s) duration. However, the observational support for this model is only indirect. A hitherto outstanding prediction is that gravitational wave (GW) events from such mergers should be associated with GRBs, and that a majority of these GRBs should be off-axis, that is, they should point away from the Earth. Here we report the discovery of the X-ray counterpart associated with the GW event GW170817. While the electromagnetic counterpart at optical and infrared frequencies is dominated by the radioactive glow from freshly synthesized r-process material in the merger ejecta, known as kilonova, observations at X-ray and, later, radio frequencies exhibit the behavior of a short GRB viewed off-axis. Our detection of X-ray emission at a location coincident with the kilonova transient provides the missing observational link between short GRBs and GWs from NS mergers, and gives independent confirmation of the collimated nature of the GRB emission.

Short-duration GRBs with extended emission form a subclass of short GRBs, comprising about 15% of the short-duration sample. Afterglow detections of short GRBs are also rare (about 30%) due to their smaller luminosity. We present a multi-band data set of the short burst with extended emission GRB 150424A, comprising of GROND observations, complemented with data from Swift/UVOT, Swift/XRT, HST, Keck/LRIS and data points from the literature. The GRB 150424A afterglow shows an extended plateau phase, lasting about 8hrs. The analysis of this unique GRB afterglow might shed light on the understanding of afterglow plateau emission, the nature of which is still under debate. We present a phenomenological analysis by applying fireball closure relations, and interpret the findings in the context of the fireball model. We discuss the plausibility of a magnetar as a central engine, being responsible for additional and prolonged energy injection into the fireball. We find convincing evidence for energy injection into the afterglow of GRB 150424A. We find that a magnetar spin down as source for a prolonged energy injection requires that at least 4% of the spin-down energy is converted to radiation.

We explain the results of Yu et al. (2015b) of the novel sharpness angle measurement to a large number of spectra obtained from the Fermi gamma-ray burst monitor. The sharpness angle is compared to the values obtained from various representative emission models: blackbody, single-electron synchrotron, synchrotron emission from a Maxwellian or power-law electron distribution. It is found that more than 91% of the high temporally and spectrally resolved spectra are inconsistent with any kind of optically thin synchrotron emission model alone. It is also found that the limiting case, a single temperature Maxwellian synchrotron function, can only contribute up to 58+23 -18% of the peak flux. These results show that even the sharpest but non-realistic case, the single-electron synchrotron function, cannot explain a large fraction of the observed spectra. Since any combination of physically possible synchrotron spectra added together will always further broaden the spectrum, emission mechanisms other than optically thin synchrotron radiation are likely required in a full explanation of the spectral peaks or breaks of the GRB prompt emission phase.

Aims. The aim of the study is to constrain the physics of gamma-ray bursts (GRBs) by analysing the multi-wavelength afterglow data set of GRB 121024A that covers the full range from radio to X-rays. Methods. Using multi-epoch broad-band observations of the GRB 121024A afterglow, we measured the three characteristic break frequencies of the synchrotron spectrum. We used six epochs of combined XRT and GROND data to constrain the temporal slopes, the dust extinction, the X-ray absorption, and the spectral slope with high accuracy. Two more epochs of combined data from XRT, GROND, APEX, CARMA, and EVLA were used to set constraints on the break frequencies and therefore on the micro-physical and dynamical parameters. Results. The XRT and GROND light curves show a simultaneous and achromatic break at around 49 ks. As a result, the crossing of the synchrotron cooling break is no suitable explanation for the break in the light curve. The multi wavelength data allow us to test two plausible scenarios explaining the break: a jet break, and the end of energy injection. The jet-break scenario requires a hard electron spectrum, a very low cooling break frequency, and a non-spreading jet. The energy injection avoids these problems, but requires $\epsilon_e > 1 (k = 2)$, spherical outflow, and $\epsilon_B < 10^{-9}$. Conclusions. In light of the extreme microphysical parameters required by the energy-injection model, we favour a jet-break scenario where $\nu_m < \nu_{sa}$ to explain the observations. This scenario gives physically meaningful microphysical parameters, and it also naturally explains the reported detection of linear and circular polarisation.

We aim to obtain high-quality time-resolved spectral fits of gamma-ray bursts (GRBs) observed by the Gamma-ray Burst Monitor (GBM) on board the Fermi Gamma-ray Space Telescope. We perform time-resolved spectral analysis with high temporal and spectral resolution of the brightest bursts observed by Fermi GBM in its first 4 years of mission. We present the complete catalog containing 1,491 spectra from 81 bursts with high spectral and temporal resolution. Distributions of parameters, statistics of the parameter populations, parameter-parameter and parameter-uncertainty correlations, and their exact values are obtained and presented as main results in this catalog. We report a criterion that is robust enough to automatically distinguish between different spectral evolutionary trends between bursts. We also search for plausible blackbody emission components and find that only 3 bursts (36 spectra in total) show evidence of a pure Planck function. It is observed that the averaged time-resolved low-energy power-law index and peak energy are slightly harder than the time-integrated values. Time-resolved spectroscopic results should be used when interpreting physics from the observed spectra, instead of the time-integrated results.

We aim to obtain a measure of the curvature of time-resolved spectra that can be compared directly to theory. This tests the ability of models such as synchrotron emission to explain the peaks or breaks of GBM prompt emission spectra. We take the burst sample from the official Fermi GBM GRB time-resolved spectral catalog. We re-fit all spectra with a measured peak or break energy in the catalog best-fit models in various energy ranges, which cover the curvature around the spectral peak or break, resulting in a total of 1,113 spectra being analysed. We compute the sharpness angles under the peak or break of the triangle constructed under the model fit curves and compare to the values obtained from various representative emission models: blackbody, single-electron synchrotron, synchrotron emission from a Maxwellian or power-law electron distribution. We find that 35% of the time-resolved spectra are inconsistent with the single-electron synchrotron function, and 91% are inconsistent with the Maxwellian synchrotron function. The single temperature, single emission time and location blackbody function is found to be sharper than all the spectra. No general evolutionary trend of the sharpness angle is observed, neither per burst nor for the whole population. It is found that the limiting case, a single temperature Maxwellian synchrotron function, can only contribute up to $58^{+23}_{-18}$% of the peak flux. Our results show that even the sharpest but non-realistic case, the single-electron synchrotron function, cannot explain a large fraction of the observed GRB prompt spectra. Because of the fact that any combination of physically possible synchrotron spectra added together will always further broaden the spectrum, emission mechanisms other than optically thin synchrotron radiation are likely required in a full explanation of the spectral peaks or breaks of the GRB prompt emission phase.

A new class of ultra-long duration (>10,000 s) gamma-ray bursts has recently been suggested. They may originate in the explosion of stars with much larger radii than normal long gamma-ray bursts or in the tidal disruptions of a star. No clear supernova had yet been associated with an ultra-long gamma-ray burst. Here we report that a supernova (2011kl) was associated with the ultra-long duration burst 111209A, at z=0.677. This supernova is more than 3 times more luminous than type Ic supernovae associated with long gamma-ray bursts, and its spectrum is distinctly different. The continuum slope resembles those of super-luminous supernovae, but extends farther down into the rest-frame ultra-violet implying a low metal content. The light curve evolves much more rapidly than super-luminous supernovae. The combination of high luminosity and low metal-line opacity cannot be reconciled with typical type Ic supernovae, but can be reproduced by a model where extra energy is injected by a strongly magnetized neutron star (a magnetar), which has also been proposed as the explanation for super-luminous supernovae.

Advances in our numerical and theoretical understanding of gamma-ray burst afterglow processes allow us to construct models capable of dealing with complex relativistic jet dynamics and non-thermal emission, that can be compared directly to data from instruments such as Swift. Because afterglow blast waves and power law spectra are intrinsically scale-invariant under changes of explosion energy and medium density, templates can be generated from large-scale hydrodynamics simulations. This allows for iterative template-based model fitting using the physical model parameters (quantifying the properties of the burster, emission and observer) directly as fit variables. Here I review how such an approach to afterglow analysis works in practice, paying special attention to the underlying model assumptions, possibilities, caveats and limitations of this type of analysis. Because some model parameters can be degenerate in certain regions of parameter space, or unconstrained if data in a limited number of a bands is available, a Bayesian approach is a natural fit. The main features of the standard afterglow model are reviewed in detail.

We present a sample of 27 GRBs with detailed Swift light curves supplemented by late time Chandra observations. To answer the missing jet-break problem in general, we develop a numerical simulation based model which can be directly fit to the data using Monte Carlo methods. Our numerical model takes into account all the factors that can shape a jet-break: (i) lateral expansion (ii) edge effects and (iii) off-axis effects. Our results provide improved fits to the light curves and constraints on physical parameters. More importantly, our results suggest that off-axis effects are important and must be included in interpretations of GRB jet breaks.

We study the time-resolved spectra of eight GRBs observed by Fermi GBM in its first five years of mission, with 1 keV - 1 MeV fluence $f>1.0\times10^{-4}$ erg cm$^{-2}$ and signal-to-noise level $\text{S/N}\geq10.0$ above 900 keV. We aim to constrain in detail the spectral properties of GRB prompt emission on a time-resolved basis and to discuss the theoretical implications of the fitting results in the context of various prompt emission models. We perform time-resolved spectral analysis using a variable temporal binning technique according to optimal S/N criteria, resulting in a total of 299 time-resolved spectra. We fit the Band function to all spectra and obtain the distributions for the low-energy power-law index $\alpha$, the high-energy power-law index $\beta$, the peak energy in the observed $\nu F_\nu$ spectrum $E_\text{p}$, and the difference between the low- and high-energy power-law indices $\Delta s=\alpha-\beta$. Using the distributions of $\Delta s$ and $\beta$, the electron population index $p$ is found to be consistent with the "moderately fast" scenario which fast- and slow-cooling scenarios cannot be distinguished. We also apply a physically motivated synchrotron model, which is a triple power-law with constrained power-law indices and a blackbody component, to test for consistency with a synchrotron origin for the prompt emission and obtain the distributions for the two break energies $E_\text{b,1}$ and $E_\text{b,2}$, the middle segment power-law index $\beta$, and the Planck function temperature $kT$. A synchrotron model is found consistent with the majority of time-resolved spectra for these eight energetic Fermi GBM bursts with good high-energy photon statistics, as long as both the cooling and injection break are included and the leftmost spectral slope is lifted either by inclusion of a thermal component or when an evolving magnetic field is accounted for.

A number of correlations between observables have been found to exist for gamma-ray burst (GRB) afterglows, linking ejecta energy to prompt and afterglow energy release and linking early stage optical and X-ray luminosity to the end times of these stages. Here, these correlations are compared to thick and thin shell models for GRB afterglows. In the thick shell model, the time evolution of the underlying relativistic blast wave is still influenced by the original ejecta, while in the thin shell model most energy in the explosion has been transferred to the external medium. It is shown here that the observed correlations rule out basic thin shell models but not the basic thick shell model. In the thick shell case, both forward shock and reverse shock dominated outflows are shown to be consistent with the correlations, using randomly generated samples of thick shell model afterglows.

Prompt or early optical emission in gamma-ray bursts is notoriously difficult to measure, and observations of the dozen cases show a large variety of properties. Yet, such early emission promises to help us achieve a better understanding of the GRB emission process(es). We performed dedicated observations of the ultra-long duration (T90 about 7000 s) GRB 130925A in the optical/near-infrared with the 7-channel "Gamma-Ray Burst Optical and Near-infrared Detector" (GROND) at the 2.2m MPG/ESO telescope. We detect an optical/NIR flare with an amplitude of nearly 2 mag which is delayed with respect to the keV--MeV prompt emission by about 300--400 s. The decay time of this flare is shorter than the duration of the flare (500 s) or its delay. While we cannot offer a straightforward explanation, we discuss the implications of the flare properties and suggest ways toward understanding it.

A sufficiently powerful astrophysical source with power law luminosity in time will give rise to a self-similar relativistic blast wave with a reverse shock traveling into the ejecta and a forward shock moving into the surrounding medium. Once energy injection ceases and the last energy is delivered to the shock front, the blast wave will transit into another self-similar stage depending only on the total amount of energy injected. I describe the effect of limited duration energy injection into environments with density depending on radius as a power law, emphasizing optical / X-ray Gamma-ray Burst afterglows as applications. The blast wave during injection is treated analytically, the transition following last energy injection with one-dimensional simulations. Flux equations for synchrotron emission from the forward and reverse shock regions are provided. The reverse shock emission can easily dominate, especially with different magnetizations for both regions. Reverse shock emission is shown to support both the reported X-ray and optical correlations between afterglow plateau duration and end time flux, independently of the luminosity power law slope. The model is demonstrated by application to bursts 120521A and 090515, and can accommodate their steep post-plateau light curve slopes.

We use the ScaleFit package to perform Markov chain Monte Carlo light curve fitting on a large sample of Swift-XRT gamma-ray burst afterglows. The ScaleFit model uses scaling relations in the hydrodynamic and radiation equations to compute synthetic light curves directly from a set of high resolution two-dimensional relativistic blast wave simulations. The data sample consists of all Swift-XRT afterglows from 2005 to 2012 with sufficient coverage and a known redshift, 188 bursts in total. We find the jet half opening angle varies widely but is commonly less than 0.1 radians. The distribution of the electron spectral index is also broad, with a median at 2.23. This approach allows, for the first time, for the off-axis observer angle to be inferred directly from the light curve. We find the observer angle to have a median value of 0.57 of the jet opening angle over our sample, which has profound consequences for the predicted rate of observed jet breaks and effects the beaming corrected total energies of gamma-ray bursts.

It is by now fairly well established that gamma-ray burst afterglows result from initially relativistic outflows interacting with the medium surrounding the burster and emitting non-thermal radiation ranging from radio to X-rays. However, beyond that, many big and small questions remain about afterglows, with the accumulating amount of observational data at the various frequencies raising as many questions as they answer. In this review I highlight a number of current theoretical issues and how they fit or do not fit within our basic theoretical framework. In addition to theoretical progress I will also emphasize the increasing role and usefulness of numerical studies of afterglow blast waves and their radiation.

Afterglows of gamma-ray bursts are observed to produce light curves with the flux following power law evolution in time. However, recent observations reveal bright flares at times on the order of minutes to days. One proposed explanation for these flares is the interaction of a relativistic blast wave with a circumburst density transition. In this paper, we model this type of interaction computationally in one and two dimensions, using a relativistic hydrodynamics code with adaptive mesh refinement called ram, and analytically in one dimension. We simulate a blast wave traveling in a stellar wind environment that encounters a sudden change in density, followed by a homogeneous medium, and compute the observed radiation using a synchrotron model. We show that flares are not observable for an encounter with a sudden density increase, such as a wind termination shock, nor for an encounter with a sudden density decrease. Furthermore, by extending our analysis to two dimensions, we are able to resolve the spreading, collimation, and edge effects of the blast wave as it encounters the change in circumburst medium. In all cases considered in this paper, we find that a flare will not be observed for any of the density changes studied.

We present preliminary data fit results of synthetic light curves computed from numerical afterglow blast wave simulations. Our technique uses Markov chain Monte Carlo (MCMC) in a new data analysis tool, ScaleFit. Scaling relations in both the hydrodynamics and radiation equations allow light curves to be parameterized by a small set of scale-invariant characteristic quantities. These quantities have been calculated and tabulated from high resolution two-dimensional hydrodynamic simulations. Producing a light curve from the characteristics takes only a millisecond, allowing for the use of MCMC data fitting techniques which can require millions of iterations. ScaleFit is a portable, lightweight, python package which performs this analysis on afterglow light curves. Using the set of Swift-XRT light curves from 2011 & 2012 with known redshifts, we find ScaleFit can measure the jet opening angle, observer angle, and spectral index of most afterglows. Globally we find gamma-ray burst afterglows tend to be observed off axis, at a significant fraction of the jet opening angle.

Afterglows of gamma-ray bursts are observed to produce light curves with the flux following power law evolution in time. However, recent observations reveal bright flares at times on the order of minutes to days. One proposed explanation for these flares is the interaction of a relativistic blast wave with a circumburst density transition. In this paper, we model this type of interaction computationally in one and two dimensions, using a relativistic hydrodynamics code with adaptive mesh refinement called RAM, and analytically in one dimension. We simulate a blast wave traveling in a stellar wind environment that encounters a sudden change in density, followed by a homogeneous medium, and compute the observed radiation using a synchrotron model. We show that flares are not observable for an encounter with a sudden density increase, such as a wind termination shock, nor for an encounter with a sudden density decrease. Furthermore, by extending our analysis to two dimensions, we are able to resolve the spreading, collimation, and edge effects of the blast wave as it encounters the change in circumburst medium. In all cases considered in this paper, we find that a flare will not be observed for any of the density changes studied.

GRB afterglow jets have been notoriously difficult to resolve numerically using 2D hydrodynamical simulations due to high outflow Lorentz factors. By performing simulations in a boosted frame, it is possible to calculate light curves from numerically computed flows in sufficient detail to accurately quantify the shape of the jet break and the post-break steepening of the light curve. We study jet breaks for jets with opening angles of theta_0 = 0.05, 0.1 and 0.2 radians decelerating in a constant density medium, observed at angles (theta_obs) ranging from on-axis to the jet edge. We present a single set of scale-invariant functions describing the time evolution of synchrotron spectral break frequencies and peak flux, depending only on theta_0 and theta_obs, sufficient to reconstruct light curves for arbitrary explosion energy, circumburst density and synchrotron slope p. We compare our light curves to fit functions using sharp breaks and smooth transitions and confirm that the measured jet break time can be delayed due to theta_obs. We find that the difference in temporal indices across the jet break is larger than theoretically anticipated leading to post-break slopes of roughly 0.25 - 1.3 p, sufficiently steep to not easily reconcile with post-break slopes measured for the Swift sample, suggesting that most Swift GRBs do not explode in homogeneous media or that jet breaks are hidden by additional physics such as prolonged energy injection or viewing angle effects. Light curves and spectral break and peak flux evolution functions will be made publicly available on-line at http://cosmo.nyu.edu/afterglowlibrary .

We present results of model fits to afterglow data sets of GRB970508, GRB980703 and GRB070125, characterized by long and broadband coverage. The model assumes synchrotron radiation (including self-absorption) from a spherical adiabatic blast wave and consists of analytic flux prescriptions based on numerical results. For the first time it combines the accuracy of hydrodynamic simulations through different stages of the outflow dynamics with the flexibility of simple heuristic formulas. The prescriptions are especially geared towards accurate description of the dynamical transition of the outflow from relativistic to Newtonian velocities in an arbitrary power-law density environment. We show that the spherical model can accurately describe the data only in the case of GRB970508, for which we find a circumburst medium density consistent with a stellar wind. We investigate in detail the implied spectra and physical parameters of that burst. For the microphysics we show evidence for equipartition between the fraction of energy density carried by relativistic electrons and magnetic field. We also find that for the blast wave to be adiabatic, the fraction of electrons accelerated at the shock has to be smaller than 1. We present best-fit parameters for the afterglows of all three bursts, including uncertainties in the parameters of GRB970508, and compare the inferred values to those obtained by different authors.

We present analytic flux prescriptions for broadband spectra of self-absorbed and optically thin synchrotron radiation from gamma-ray burst afterglows, based on one-dimensional relativistic hydrodynamic simulations. By treating the evolution of critical spectrum parameters as a power-law break between the ultrarelativistic and non-relativistic asymptotic solutions, we generalize the prescriptions to any observer time. Our aim is to provide a set of formulas that constitutes a useful tool for accurate fitting of model-parameters to observational data, regardless of the dynamical phase of the outflow. The applicability range is not confined to gamma-ray burst afterglows, but includes all spherical outflows (also jets before the jet-break) that produce synchrotron radiation as they adiabatically decelerate in a cold, power-law medium. We test the accuracy of the prescriptions and show that numerical evidence suggests that typical relative errors in the derivation of physical quantities are about 10 per cent. A software implementation of the presented flux prescriptions combined with a fitting code is freely available on request and on-line. Together they can be used in order to directly fit model parameters to data.

We present a powerful new tool for fitting broadband gamma-ray burst afterglow data, which can be used to determine the burst explosion parameters and the synchrotron radiation parameters. By making use of scale invariance between relativistic jets of different energies and different circumburst medium densities, and by capturing the output of high-resolution two-dimensional relativistic hydrodynamical (RHD) jet simulations in a concise summary, the jet dynamics are generated quickly. Our method calculates the full light curves and spectra using linear radiative transfer sufficiently fast to allow for a direct iterative fit of RHD simulations to the data. The fit properly accounts for jet features that so far have not been successfully modeled analytically, such as jet decollimation, inhomogeneity along the shock front and the transitory phase between the early time relativistic and late time non-relativistic outflow. As a first application of the model we simultaneously fit the radio, X-ray and optical data of GRB 990510. We not only find noticeable differences between our findings for the explosion and radiation parameters and those of earlier authors, but also an improved model fit when we include the observer angle in the data fit. The fit method will be made freely available on request and on-line at http://cosmo.nyu.edu/afterglowlibrary. In addition to data fitting, the software tools can also be used to quickly generate a light curve or spectrum for arbitrary observer position, jet and radiation parameters.

We discuss jet dynamics for narrow and wide gamma-ray burst (GRB) afterglow jets and the observational implications of numerical simulations of relativistic jets in two dimensions. We confirm earlier numerical results that sideways expansion of relativistic jets during the bulk of the afterglow emission phase is logarithmic in time and find that this also applies to narrow jets with half opening angle of 0.05 radians. As a result, afterglow jets remain highly nonspherical until after they have become nonrelativistic. Although sideways expansion steepens the afterglow light curve after the jet break, the jet edges becoming visible dominates the jet break, which means that the jet break is sensitive to the observer angle even for narrow jets. Failure to take the observer angle into account can lead to an overestimation of the jet energy by up to a factor 4. This weakens the challenge posed to the magneter energy limit by extreme events such as GRB090926A. Late time radio calorimetry based on a spherical nonrelativistic outflow model remains relevant when the observer is approximately on-axis and where differences of a few in flux level between the model and the simulation are acceptable. However, this does not imply sphericity of the outflow and therefore does not translate to high observer angles relevant to orphan afterglows. For more accurate calorimetry and in order to model significant late time features such as the rise of the counterjet, detailed jet simulations remain indispensable.

We demonstrate that gamma-ray burst afterglow spectra and light curves can be calculated for arbitrary explosion and radiation parameters by scaling the peak flux and the critical frequencies connecting different spectral regimes. Only one baseline calculation needs to be done for each jet opening angle and observer angle. These calculations are done numerically using high-resolution relativistic hydrodynamical afterglow blast wave simulations which include the two-dimensional dynamical features of expanding and decelerating afterglow blast waves. Any light curve can then be generated by applying scaling relations to the baseline calculations. As a result, it is now possible to fully fit for the shape of the jet break, e.g. at early time X-ray and optical frequencies. In addition, late-time radio calorimetry can be improved since the general shape of the transition into the Sedov-Taylor regime is now known for arbitrary explosion parameters so the exact moment when the Sedov-Taylor asymptote is reached in the light curve is no longer relevant. When calculating the baselines, we find that the synchrotron critical frequency and the cooling break frequency are strongly affected by the jet break. The synchrotron break temporal slope quickly drops to the steep late time Sedov-Taylor slope, while the cooling break first steepens then rises to meet the level of its shallow late time asymptote.

We present results for a large number of gamma-ray burst (GRB) afterglow light curve calculations, done by combining high resolution two-dimensional relativistic hydrodynamics simulations using RAM with a synchrotron radiation code. Results were obtained for jet energies, circumburst medium densities and jet angles typical for short and underluminous GRBs, different observer angles and observer frequencies from low radio (75 MHz) to X-ray (1.5 keV). We summarize the light curves through smooth power law fits with up to three breaks, covering jet breaks for small observer angles, the rising phase for large observer angles and the rise and decay of the counterjet. All light curve data are publicly available via http://cosmo.nyu.edu/afterglowlibrary . The data can be used for model fits to observational data and as an aid for predicting observations by future telescopes such as LOFAR or SKA and will benefit the study of neutron star mergers using different channels, such as gravitational wave observations with LIGO or Virgo. For small observer angles, we find jet break times that vary significantly between frequencies, with the break time in the radio substantially postponed. Increasing the observer angle also postpones the measured jet break time. The rising phase of the light curve for large observer angle has a complex shape that can not always be summarized by a simple power law. Except for very large observer angles, the counter jet is a distinct feature in the light curve, although in practice the signal will be exceedingly difficult to observe by then.

Strong optical and radio flares often appear in the afterglow phase of Gamma-Ray Bursts (GRBs). It has been proposed that colliding ultra-relativistic shells can produce these flares. Such consecutive shells can be formed due to the variability in the central source of a GRB. We perform high resolution 1D numerical simulations of late collisions between two ultra-relativistic shells in order to explore these events. We examine the case where a cold uniform shell collides with a self-similar Blandford and McKee shell in a constant density environment and consider cases with different Lorentz factor and energy for the uniform shell. We produce the corresponding on-axis light curves and emission images for the afterglow phase and examine the occurrence of optical and radio flares assuming a spherical explosion and a hard-edged jet scenario. For our simulations we use the Adaptive Mesh Refinement version of the Versatile Advection Code (AMRVAC) coupled to a linear radiative transfer code to calculate synchrotron emission. We find steeply rising flare like behavior for small jet opening angles and more gradual rebrightenings for large opening angles. Synchrotron self-absorption is found to strongly influence the onset and shape of the radio flare.

Numerical studies of gamma-ray burst afterglow jets reveal significant qualitative differences with simplified analytical models. We present an on-line library of synthetic afterglow light curves and broadband spectra for use in interpreting observational data. Light curves have been calculated for various physics settings such as explosion energy and circumburst structure, as well as differing jet parameters and observer angle and redshift. Calculations gave been done for observer frequencies ranging from low radio to X-ray and for observer times from hours to decades after the burst. The light curves have been calculated from high-resolution 2D hydrodynamical simulations performed with the RAM adaptive-mesh refinement code and a detailed synchrotron radiation code. The library will contain both generic afterglow simulations as well as specific case studies and will be freely accessible at http://cosmo.nyu.edu/afterglowlibrary . The synthetic light curves can be used as a check on the accuracy of physical parameters derived from analytical model fits to afterglow data, to quantitatively explore the consequences of varying parameters such as observer angle and for accurate predictions of future telescope data.

Numerical jet simulations serve a valuable role in calculating gamma-ray burst afterglow emission beyond analytical approximations. Here we present the results of high resolution 2D simulations of decelerating relativistic jets performed using the RAM adaptive mesh refinement relativistic hydrodynamics code. We have applied a separate synchrotron radiation code to the simulation results in order to calculate light curves at frequencies varying from radio to X-ray for observers at various angles from the jet axis. We provide a confirmation from radio light curves from simulations rather than from a simplified jet model for earlier results in the literature finding that only a very small number of local Ibc supernovae can possibly harbor an orphan afterglow. Also, recent studies have noted an unexpected lack of observed jet breaks in the Swift sample. Using a jet simulation with physical parameters representative for an average Swift sample burst, such as a jet half opening angle of 0.1 rad and a source redshift of z = 2.23, we have created synthetic light curves at 1.5 keV with artificial errors while accounting for Swift instrument biases as well. A large set of these light curves have been generated and analyzed using a Monte Carlo approach. Single and broken power law fits are compared. We find that for increasing observer angle, the jet break quickly becomes hard to detect. This holds true even when the observer remains well within the jet opening angle. We find that the odds that a Swift light curve from a randomly oriented 0.1 radians jet at z = 2.23 will exhibit a jet break at the 3 sigma level are only 12 percent. The observer angle therefore provides a natural explanation for the lack of perceived jet breaks in the Swift sample.

We present a study of the intermediate regime between ultra-relativistic and nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of spherically symmetric blast waves is numerically calculated using the AMRVAC adaptive mesh refinement code. Spectra and light curves are calculated using a separate radiation code that, for the first time, links a parametrisation of the microphysics of shock acceleration, synchrotron self-absorption and electron cooling to a high-performance hydrodynamics simulation. For the dynamics we find that the transition to the nonrelativistic regime generally occurs later than expected, that the Sedov-Taylor solution overpredicts the late time blast wave radius and that the analytical formula for the blast wave velocity from Huang (1999) overpredicts the late time velocity by a factor 4/3. For the radiation we find that the flux may differ up to an order of magnitude depending on the equation of state that is used for the fluid and that the counterjet leads to a clear rebrightening at late times for hard-edged jets. Simulating GRB030329 using predictions for its physical parameters from the literature leads to spectra and light curves that may differ significantly from the actual data, emphasizing the need for very accurate modelling. Predicted light curves at low radio frequencies for a hard-edged jet model of GRB030329 with opening angle 22 degrees show typically two distinct peaks, due to the combined effect of jet break, non relativistic break and counterjet.

Starting as highly relativistic collimated jets, gamma-ray burst outflows gradually decelerate and become non-relativistic spherical blast waves. Although detailed analytical solutions describing the afterglow emission received by an on-axis observer during both the early and late phases of the outflow evolution exist, a calculation of the received flux during the intermediate phase and for an off-axis observer requires either a more simplified analytical model or direct numerical simulations of the outflow dynamics. In this paper we present light curves for off-axis observers covering the long-term evolution of the blast wave calculated from a high resolution two-dimensional relativistic hydrodynamics simulation using a synchrotron radiation model. We compare our results to earlier analytical work and calculate the consequence of the observer angle with respect to the jet axis both for the detection of orphan afterglows and for jet break fits to the observational data. We find that observable jet breaks can be delayed for up to several weeks for off-axis observers, potentially leading to overestimation of the beaming corrected total energy. When using our off-axis light curves to create synthetic Swift X-ray data, we find that jet breaks are likely to remain hidden in the data. We also confirm earlier results in the literature finding that only a very small number of local Type Ibc supernovae can harbor an orphan afterglow.

The conventional derivation of the gamma-ray burst afterglow jet break time uses only the blast wave fluid Lorentz factor and therefore leads to an achromatic break. We show that in general gamma-ray burst afterglow jet breaks are chromatic across the self-absorption break. Depending on circumstances, the radio jet break may be postponed significantly. Using high-accuracy adaptive mesh fluid simulations in one dimension, coupled to a detailed synchrotron radiation code, we demonstrate that this is true even for the standard fireball model and hard-edged jets. We confirm these effects with a simulation in two dimensions. The frequency dependence of the jet break is a result of the angle dependence of the emission, the changing optical depth in the self-absorbed regime and the shape of the synchrotron spectrum in general. In the optically thin case the conventional analysis systematically overestimates the jet break time, leading to inferred opening angles that are underestimated by a factor 1.32 and explosion energies that are underestimated by a factor 1.73, for explosions in a homogeneous environment.

We present a study of the intermediate regime between ultra-relativistic and nonrelativistic flow for gamma-ray burst afterglows. The hydrodynamics of spherically symmetric blast waves is numerically calculated using the AMRVAC adaptive mesh refinement code. Spectra and light curves are calculated using a separate radiation code that, for the first time, links a parametrisation of the microphysics of shock acceleration, synchrotron self-absorption and electron cooling to a high-performance hydrodynamics simulation.

Gamma-ray burst afterglow flares and rebrightenings of the optical and X-ray light curve have been attributed to both late time inner engine activity and density changes in the medium surrounding the burster. To test the latter, we study the encounter between the relativistic blast wave from a gamma-ray burster and a stellar wind termination shock. The blast wave is simulated using a high performance adaptive mesh relativistic hydrodynamics code, AMRVAC, and the synchrotron emission is analyzed in detail with a separate radiation code. We find no bump in the resulting light curve, not even for very high density jumps. Furthermore, by analyzing the contributions from the different shock wave regions we are able to establish that it is essential to resolve the blast wave structure in order to make qualitatively correct predictions on the observed output and that the contribution from the reverse shock region will not stand out, even when the magnetic field is increased in this region by repeated shocks. This study resolves a controversy in recent literature.

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. We also show the light curves resulting from a relativistic blast wave encountering a wind termination shock. From high resolution calculations we find that the observed transition from a stellar wind type light curve to an interstellar medium type light curve is smooth and without short-time transitory features.

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