Off-Axis Afterglow Light Curves from High-Resolution Hydrodynamical Jet Simulations

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.Comment: 4 pages, 3 figures. First of two contributions to proceedings GRB2010 Maryland conference. Editors: McEnery, Racusin and Gehrels. The data from this paper is publicly available from http://cosmo.nyu.edu/afterglowlibrary

An on-line library of afterglow light curves

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.Comment: 4 pages, 2 figures. Second of two contributions to proceedings GRB2010 Maryland conference. Editors: McEnery, Racusin and Gehrels. The data from this paper is publicly available from http://cosmo.nyu.edu/afterglowlibrary

Fitting afterglows with multi-dimensional simulations

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 scaleinvariant 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

Gamma-ray burst afterglow theory

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.Comment: Invited review for 7th Huntsville Gamma-Ray Burst Symposium, GRB 2013: paper 24 in eConf Proceedings C130414

GRB afterglow blast wave encountering sudden circumburst density change produces no flares

Afterglows of gamma-ray bursts are observed to produce light curveswith 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 relativisticblast wave with a circumburst density transition. In this paper, we modelthis type of interaction computationally in one and two dimensions, usinga relativistic hydrodynamics code with adaptive mesh refinement calledram, and analytically in one dimension. We simulate a blast wave travelingin a stellar wind environment that encounters a sudden change indensity, followed by a homogeneous medium, and compute the observedradiation using a synchrotron model. We show that flares are not observablefor an encounter with a sudden density increase, such as a windtermination shock, nor for an encounter with a sudden density decrease.Furthermore, by extending our analysis to two dimensions, we are able toresolve the spreading, collimation, and edge effects of the blast wave as itencounters the change in circumburst medium. In all cases considered inthis paper, we find that a flare will not be observed for any of the densitychanges studied

From blast wave to observation

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.Comment: conference proceedings 6th Huntsville symposium, 20-23 October 2008. Editors: C.A. Meegan, N. Gehrels, and C. Kouvelioto