We outline the science prospects for gamma-ray bursts (GRBs) with the Cherenkov Telescope Array (CTA), the next-generation ground-based gamma-ray observatory operating at energies above few tens of GeV. With its low energy threshold, large effective
area and rapid slewing capabilities, CTA will be able to measure the spectra and variability of GRBs at multi-GeV energies with unprecedented photon statistics, and thereby break new ground in elucidating the physics of GRBs, which is still poorly understood. Such measurements will also provide crucial diagnostics of ultra-high-energy cosmic ray and neutrino production in GRBs, advance observational cosmology by probing the high-redshift extragalactic background light and intergalactic magnetic fields, and contribute to fundamental physics by testing Lorentz invariance violation with high precision. Aiming to quantify these goals, we present some simulated observations of GRB spectra and light curves, together with estimates of their detection rates with CTA. Although the expected detection rate is modest, of order a few GRBs per year, hundreds or more high-energy photons per burst may be attainable once they are detected. We also address various issues related to following up alerts from satellites and other facilities with CTA, as well as follow-up observations at other wavelengths. The possibility of discovering and observing GRBs from their onset including short GRBs during a wide-field survey mode is also briefly discussed.
We report on the development of Mezcal-SRHD, a new adaptive mesh refinement, special relativistic hydrodynamics (SRHD) code, developed with the aim of studying the highly relativistic flows in Gamma-Ray Burst sources. The SRHD equations are solved us
ing finite volume conservative solvers. The correct implementation of the algorithms is verified by one-dimensional (1D) shock tube and multidimensional tests. The code is then applied to study the propagation of 1D spherical impulsive blast waves expanding in a stratified medium with $rho propto r^{-k}$, bridging between the relativistic and Newtonian phases, as well as to a two-dimensional (2D) cylindrically symmetric impulsive jet propagating in a constant density medium. It is shown that the deceleration to non-relativistic speeds in one-dimension occurs on scales significantly larger than the Sedov length. This transition is further delayed with respect to the Sedov length as the degree of stratification of the ambient medium is increased. This result, together with the scaling of position, Lorentz factor and the shock velocity as a function of time and shock radius, is explained here using a simple analytical model based on energy conservation. The method used for calculating the afterglow radiation by post-processing the results of the simulations is described in detail. The light curves computed using the results of 1D numerical simulations during the relativistic stage correctly reproduce those calculated assuming the self-similar Blandford-McKee solution for the evolution of the flow. The jet dynamics from our 2D simulations and the resulting afterglow lightcurves, including the jet break, are in good agreement with those presented in previous works. Finally, we show how the details of the dynamics critically depend on properly resolving the structure of the relativistic flow.
