GRB 131108A is a bright long Gamma-Ray Burst (GRB) detected by the Large Area Telescope and the Gamma-ray Burst Monitor on board the \textit{Fermi Gamma-ray Space Telescope}. Dedicated temporal and spectral analyses reveal three $\gamma$-ray flares dominating above 100 MeV, which are not directly related to the prompt emission in the GBM band (10 keV--10 MeV). The high-energy light curve of GRB 131108A (100 MeV -- 10 GeV) shows an unusual evolution: a steep decay, followed by three flares with an underlying emission, and then a long-lasting decay phase. The detailed analysis of the $\gamma$-ray flares finds that the three flares are 6 -- 20 times brighter than the underlying emission and are similar to each other. The fluence of each flare, (1.6 $\sim$ 2.0) $\times$ 10$^{-6}$ erg cm$^{-2}$, is comparable to that of emission during the steep decay phase, 1.7 $\times$ 10$^{-6}$ erg cm$^{-2}$. The total fluence from three $\gamma$-ray flares is 5.3 $\times$ 10$^{-6}$ erg cm$^{-2}$. The three $\gamma$-ray flares show properties similar to the usual X-ray flares that are sharp flux increases, occurring in $\sim$ 50\% of afterglows, in some cases well after the prompt emission. Also, the temporal and spectral indices during the early steep decay phase and the decaying phase of each flare show the consistency with a relation of the curvature effect ($\hat{\alpha}$ = 2 + $\hat{\beta}$), which is the first observational evidence of the high-latitude emission in the GeV energy band.

The characteristics of the spectral evolution of the prompt emission of gamma-ray bursts (GRBs), which are closely related to the radiation mechanism (synchrotron or photosphere), are still an unsolved subject. Here, by performing the detailed time-resolved spectral fitting of GRB 131231A, which has a very bright and well-defined single pulse, some interesting spectral evolution features have been found. (i) Both the low-energy spectral index $\alpha$ and the peak energy $E_{\rm p}$ exhibit the "flux-tracking" pattern ("double-tracking" characteristics). (ii) The parameter relations, i.e., $F$ (the energy flux)-$\alpha$, $F$-$E_{\rm p}$, and $E_{\rm p}$-$\alpha$, along with the analogous Yonetoku $E_{\rm p}$-$L_{\gamma,\rm iso}$ relation for the different time-resolved spectra, show strong monotonous (positive) correlations, both in the rising and the decaying phases. (iii) The values of $\alpha$ do not exceed the synchrotron limit ($\alpha$= -2/3) in all slices across the pulse, favoring the synchrotron origin. We argue that the one-zone synchrotron emission model with the emitter streaming away at a large distance from the central engine can explain all of these special spectral evolution characteristics.

The synchrotron external shock model predicts the evolution of the spectral ($\beta$) and temporal ($\alpha$) indices during the gamma-ray burst (GRB) afterglow for different environmental density profiles, electron spectral indices, electron cooling regimes, and regions of the spectrum. We study the relationship between $\alpha$ and $\beta$, the so-called "closure relations" with GRBs detected by \textit{Fermi} Large Area Telescope (\textit{Fermi}-LAT) from 2008 August to 2018 August. The spectral and temporal indices for the > 100 MeV emission from the \textit{Fermi}-LAT as determined in the Second Fermi-LAT Gamma-ray Burst Catalog (2FLGC) are used in this work. We select GRBs whose spectral and temporal indices are well constrained (58 long-duration GRBs and 1 short-duration GRBs) and classify each GRB into the best-matched relation. As a result, we found that a number of GRBs require a very small fraction of the total energy density contained in the magnetic field ($\epsilon_{B}$ $\lesssim$ 10$^{-7}$). The estimated mean and standard deviation of electron spectral index $\mathit{p}$ are 2.40 and 0.44, respectively. The GRBs satisfying a closure relation of the slow cooling tend to have a softer $\mathit{p}$ value compared to those of the fast cooling. Moreover, the Kolmogorov--Smirnov test of the two $\mathit{p}$ distributions from the fast and slow coolings rejects a hypothesis that the two distributions are drawn from the single reference distribution with a significance of 3.2 $\sigma$. Lastly, the uniform density medium is preferred over the medium that decreases like the inverse of distance squared for long-duration GRBs.

GRB 160709A is one of the few bright short gamma-ray bursts detected by both the Gamma-ray Burst Monitor and the Large Area Telescope on board the $Fermi$ $Gamma$-$ray$ $Space$ $Telescope$. The $\gamma$-ray prompt emission of GRB 160709A is adequately fitted by combinations of three distinct components: (i) a nonthermal component described by a power law (PL) with a high-energy exponential cutoff, (ii) a thermal component modeled with a Planck function, and (iii) a second nonthermal component shaped by an additional PL crossing the whole $\gamma$-ray spectrum. While the thermal component dominates during $\sim$ 0.12 s of the main emission episode of GRB 160709A with an unusually high temperature of $\sim$ 340 keV, the nonthermal components dominate in the early and late time. The thermal component is consistent with the photospheric emission resulting in the following parameters: the size of the central engine, $R_{0}$ = $3.8 \substack{+5.9 \\ -1.8}$ $\times 10^{8}$ cm, the size of the photosphere, R$_{ph}$ = $7.4 \substack{+0.8 \\ -1.2}$ $\times 10^{10}$ cm, and a bulk Lorentz factor, $\Gamma$ = $728 \substack{+75 \\ -93}$ assuming a redshift of 1. The slope of the additional PL spectrum stays unchanged throughout the burst duration; however, its flux decreases continuously as a function of time. A standard external shock model has been tested for the additional PL component using the relation between the temporal and spectral indices (the closure relation). Each set of spectral and temporal indices from two energy bands (200 keV--40 MeV and 100 MeV--10 GeV) satisfies a distinct closure relation. From the closure relation test we derived the index for the electron spectral distribution, $p$ = 2.5 $\pm$ 0.1. The interaction of the jet with the interstellar environment is preferred over the interaction with the wind medium.

The prompt emission phase of gamma-ray bursts (GRBs) exhibits two distinct patterns of the peak-energy ($E_p$) evolution, i.e., time-resolved spectral analyses of $\nu F_{\nu}$ spectra of broad pulses reveal (1) "hard-to-soft" and (2) "flux-tracking" patterns of $E_p$ evolution in time, the physical origin of which still remains not well understood. We show here that these two patterns can be successfully reproduced within a simple physical model invoking synchrotron radiation in a bulk-accelerating emission region. We show further that the evolution patterns of the peak energy have, in fact, direct connections to the existence of two different (positive or negative) types of spectral lags, seen in the broad pulses. In particular, we predict that (1) only the positive type of spectral lags is possible for the hard-to-soft evolution of the peak energy, (2) both the positive and negative type of spectral lags can occur in the case of flux-tracking pattern of the peak energy, (3) for the flux-tracking pattern, the peak location of the flux light curve slightly lags behind the peak of the $E_p$ evolution with time if the spectral lags are positive, and (4) in the case of flux-tracking pattern, double-peaked broad pulses can appear in the light curves, the shape of which is energy-dependent.

Applying our recently-developed generalized version of the high-latitude emission theory to the observations of X-ray flares in gamma-ray bursts (GRBs), we present here clear observational evidence that the X-ray flare emission region is undergoing rapid bulk acceleration as the photons are emitted. We show that both the observed X-ray flare light curves and the photon index evolution curves can be simultaneously reproduced within a simple physical model invoking synchrotron radiation in an accelerating emission region far from the GRB central engine. Such an acceleration process demands an additional energy dissipation source other than kinetic energy, which points towards a significant Poynting-flux in the emission region of X-ray flares. As the X-ray flares are believed to share a similar physical mechanism as the GRB prompt emission, our finding here hints that the GRB prompt emission jets may also carry a significant Poynting-flux in their emitting region.

When emission in a conical relativistic jet ceases abruptly (or decays sharply), the observed decay light curve is controlled by the high-latitude "curvature effect". Recently, Uhm & Zhang found that the decay slopes of three GRB X-ray flares are steeper than what the standard model predicts. This requires bulk acceleration of the emission region, which is consistent with a Poynting-flux-dominated outflow. In this paper, we systematically analyze a sample of 85 bright X-ray flares detected in 63 Swift GRBs, and investigate the relationship between the temporal decay index $\alpha$ and spectral index $\beta$ during the steep decay phase of these flares. The $\alpha$ value depends on the choice of the zero time point $t_0$. We adopt two methods. "Method I" takes $t_0^I$ as the first rising data point of each flare, and is the most conservative approach. We find that at 99.9% condifence level 56/85 flares have decay slopes steeper than the simplest curvature effect prediction, and therefore, are in the acceleration regime. "Method II" extrapolates the rising light curve of each flare backwards until the flux density is three orders of magnitude lower than the peak flux density, and defines the corresponding time as the time zero point (t_0^II). We find that 74/85 flares fall into the acceleration regime at 99.9% condifence level. This suggests that bulk acceleration is common, may be even ubiquitous among X-ray flares, pointing towards a Poynting-flux-dominated jet composition for these events.

Despite decades of investigations, the physical mechanism that powers the bright prompt $\gamma$-ray emission from gamma-ray bursts (GRBs) is still not identified. One important observational clue that remains not properly interpreted so far is the existence of time lags of broad light curve pulses in different energy bands, named "spectral lags". Here we show that the traditional view invoking the high-latitude emission "curvature effect" of a relativistic jet cannot account for spectral lags. Rather, the observed spectral lags demand the sweep of a spectral peak across the observing energy band in a specific manner. The duration of the broad pulses and inferred typical Lorentz factor of GRBs require that the emission region is in an optically thin emission region far from the GRB central engine. We construct a simple physical model invoking synchrotron radiation from a rapidly expanding outflow. We show that the observed spectral lags appear naturally in our model light-curves given that (1) the gamma-ray photon spectrum is curved (as observed), (2) the magnetic field strength in the emitting region decreases with radius as the region expands in space, and (3) the emission region itself undergoes rapid bulk acceleration as the prompt $\gamma$-rays are produced. These requirements are consistent with a Poynting-flux-dominated jet abruptly dissipating magnetic energy at a large distance from the engine.

We consider a relativistic spherical shell and calculate its spectral flux as received by a distant observer. Using two different methods, we derive a simple analytical expression of the observed spectral flux and show that the well-known relation $\hat \alpha = 2+\hat \beta$ (between temporal index $\hat \alpha$ and spectral index $\hat \beta$) of the high-latitude emission is achieved naturally in our derivation but holds only when the shell moves with a constant Lorentz factor $\Gamma$. Presenting numerical models where the shell is under acceleration or deceleration, we show that the simple $\hat \alpha = 2+\hat \beta$ relation is indeed deviated as long as $\Gamma$ is not constant. For the models under acceleration, we find that the light curves produced purely by the high-latitude emission decay initially much steeper than the constant $\Gamma$ case and gradually resume the $\hat \alpha = 2+\hat \beta$ relation in about one and half orders of magnitude in observer time. For the models under deceleration, the trend is opposite. The light curves made purely by the high-latitude emission decay initially shallower than the constant $\Gamma$ case and gradually resume the relation $\hat \alpha = 2+\hat \beta$ in a similar order of magnitude in observer time. We also show that how fast the Lorentz factor $\Gamma$ of the shell increases or decreases is the main ingredient determining the initial steepness or shallowness of the light curves.

We perform a time-resolved spectral analysis of GRB 130606B within the framework of a fast-cooling synchrotron radiation model with magnetic field strength in the emission region decaying with time, as proposed by Uhm & Zhang. The data from all time intervals can be successfully fit by the model. The same data can be equally well fit by the empirical Band function with typical parameter values. Our results, which involve only minimal physical assumptions, offer one natural solution to the origin of the observed GRB spectra and imply that, at least some, if not all, Band-like GRB spectra with typical Band parameter values can indeed be explained by synchrotron radiation.

We present a wide dataset of gamma-ray, X-ray, UVOIR, and radio observations of the Swift GRB100814A. At the end of the slow decline phase of the X-ray and optical afterglow, this burst shows a sudden and prominent rebrightening in the optical band only, followed by a fast decay in both bands. The optical rebrightening also shows chromatic evolution. Such a puzzling behaviour cannot be explained by a single component model. We discuss other possible interpretations, and we find that a model that incorporates a long-lived reverse shock and forward shock fits the temporal and spectral properties of GRB100814A the best.

We investigate the dynamics and afterglow light curves of gamma-ray burst (GRB) blast waves that encounter various density structures (such as bumps, voids, or steps) in the surrounding ambient medium. We present and explain the characteristic response features that each type of density structures in the medium leaves on the forward shock (FS) and reverse shock (RS) dynamics, for blast waves with either a long-lived or short-lived RS. We show that, when the ambient medium density drops, the blast waves exhibit in some cases a period of an actual acceleration (even during their deceleration stage), due to adiabatic cooling of blast waves. Comparing numerical examples that have different shapes of bumps or voids, we propose a number of consistency tests that correct modeling of blast waves needs to satisfy. Our model results successfully pass these tests. Employing a Lagrangian description of blast waves, we perform a sophisticated calculation of afterglow emission. We show that, as a response to density structures in the ambient medium, the RS light curves produce more significant variations than the FS light curves. Some observed features (such as re-brightenings, dips, or slow wiggles) can be more easily explained within the RS model. We also discuss on the origin of these different features imprinted on the FS and RS light curves.

Synchrotron radiation of relativistic electrons is an important radiation mechanism in many astrophysical sources. In the sources where the synchrotron cooling time scale $t_c$ is shorter than the dynamical time scale $t_{dyn}$, electrons are cooled down below the minimum injection energy. It has been believed that such "fast cooling" electrons have an energy distribution $dN_e /d\gamma_e \propto \gamma_e^{-2}$, and their synchrotron radiation flux density has a spectral shape $F_\nu \propto \nu^{-1/2}$. On the other hand, in a transient expanding astrophysical source, such as a gamma-ray burst (GRB), the magnetic field strength in the emission region continuously decreases with radius. Here we study such a system, and find that in a certain parameter regime, the fast cooling electrons can have a harder energy spectrum, and the standard $d N_e / d \gamma_e \propto \gamma_e^{-2}$ spectrum is achieved only in the deep fast cooling regime when $t_c \ll t_{dyn}$. We apply this new physical regime to GRBs, and suggest that the GRB prompt emission spectra whose low-energy photon index $\alpha$ has a typical value -1 could be due to synchrotron radiation in this moderately fast cooling regime.

Although the widely-used analytical afterglow model of gamma-ray bursts (GRBs) predicts a sharp cooling break $\nu_c$ in its afterglow spectrum, the GRB observations so far rarely show clear evidence for a cooling break in their spectra or its corresponding temporal break in their light curves. Employing a Lagrangian description of the blast wave, we conduct a sophisticated calculation of the afterglow emission. We precisely follow the cooling history of non-thermal electrons accelerated into each Lagrangian shell. We show that a detailed calculation of afterglow spectra does not in fact give rise to a sharp cooling break at $\nu_c$. Instead, it displays a very mild and smooth transition, which occurs gradually over a few orders of magnitude in energy or frequency. The main source of this slow transition is that different mini-shells have different evolution histories of the comoving magnetic field strength $B$, so that deriving the current value of $\nu_c$ of each mini-shell requires an integration of its cooling rate over the time elapsed since its creation. We present the time evolution of optical and X-ray spectral indices to demonstrate the slow transition of spectral regimes, and discuss the implications of our result in interpreting GRB afterglow data.

We perform a detailed study on the dynamics of a relativistic blast wave with the presence of a long-lived reverse shock (RS). Although a short-lived RS has been widely considered, the RS is believed to be long-lived as a consequence of a stratification expected on the ejecta Lorentz factors. The existence of a long-lived RS makes the forward shock (FS) dynamics to deviate from a self-similar Blandford-McKee solution. Employing the "mechanical model" that correctly incorporates the energy conservation, we present an accurate solution for both the FS and RS dynamics. We conduct a sophisticated calculation of the afterglow emission. Adopting a Lagrangian description of the blast wave, we keep track of an adiabatic evolution of numerous shells between the FS and RS. An evolution of the electron spectrum is also followed individually for every shell. We then find the FS and RS light curves by integrating over the entire FS and RS shocked regions, respectively. In particular, we make use of an analytic expression for observed spectral flux, which we derive here in terms of an observed frequency and observer time. Exploring a total of 20 different ejecta stratifications, we explain in detail how a stratified ejecta affects its blast wave dynamics and afterglow light curves. We show that, while the FS light curves are not sensitive to the ejecta stratifications, the RS light curves exhibit much richer features, including steep declines, plateaus, bumps, re-brightenings, and a variety of temporal decay indices. These distinctive RS features may be observable if the RS has higher values of the micophysics parameters than the FS. We discuss possible applications of our results in understanding the GRB afterglow data.

We present the gamma-ray, X-ray, optical and radio data for GRB100814A. At the end of the slow decline phase of the X-ray and optical afterglow, a sudden and prominent rebrightening in the optical band occurs followed by a fast decay in both bands. This optical rebrightening is accompanied by possible chromatic variations. We discuss possible interpretations, such as double component scenarios and internal dissipation mechanism, with their virtues and drawbacks. We also compare GRB100814A with other Swift bursts that show optical rebrightenings with similar properties.

A few long gamma-ray bursts such as GRB 050421 show no afterglow emission beyond the usual initial steep decay phase. It has been suggested that these events correspond to "naked" bursts that occur in a very low density environment. We reconsider this possibility in the context of various scenarios for the origin of the afterglow. In the standard model where the afterglow results from the forward shock as well as in the alternative model where the afterglow comes from the reverse shock, we aim to obtain constraints on the density of the environment, the microphysics parameters, or the Lorentz factor of the ejecta, which are imposed by the absence of a detected afterglow. For the two models we compute the afterglow evolution for different values of the external density (uniform or wind medium) and various burst parameters. We then compare our results to the Swift data of GRB 050421, which is the best example of a long burst without afterglow. In the standard model we show that consistency with the data imposes that the external density does not exceed 1E-5 cm-3 or that the microphysics parameters are very small with epsilon_e <~ 1E-2 and epsilon_B <~ 1E-4. If the afterglow is caused by the reverse shock, we find that its contribution can be strongly reduced if the central source has mainly emitted fast-moving material (with less than 10 - 30 % of the kinetic energy at Gamma<100 and was located in a dense environment. The two considered scenarios therefore lead to opposite constraints on the circumburst medium. The high-density environment, favored by the reverse shock model, better corresponds to what is expected if the burst progenitor was a massive star.

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.

The curvature of a relativistic blast wave implies that its emission arrives to observers with a spread in time. This effect is believed to wash out fast variability in the lightcurves of GRB afterglows. We note that the spreading effect is reduced if emission is anisotropic in the rest-frame of the blast wave (i.e. if emission is limb-brightened or limb-darkened). In particular, synchrotron emission is almost certainly anisotropic, and may be strongly anisotropic, depending on details of electron acceleration in the blast wave. Anisotropic afterglows can display fast and strong variability at high frequencies (above the 'fast-cooling' frequency). This may explain the existence of bizarre features in the X-ray afterglows of GRBs, such as sudden drops and flares. We also note that a moderate anisotropy can significantly delay the 'jet break' in the lightcurve, which makes it harder to detect.

Hundreds of gamma-ray burst (GRB) UV-optical light curves have been measured since the discovery of optical afterglows, however, even after nearly 5 years of operation of the SWIFT observatory, only a handful of measurements have been made soon (within a minute) after the gamma ray signal. This lack of early observations fails to address burst physics at the short time scales associated with burst events and progenitors. Because of this lack of sub-minute data, the characteristics of the UV-optical light curve of short-hard type GRB and rapid-rising GRB, which may account for ~30% of all GRB, remain practically unknown. We have developed methods for reaching the sub-minute and the sub-second timescales in a small spacecraft observatory appropriate for launch on a microsatellite. Rather than slewing the entire spacecraft to aim the UV-optical instrument at the GRB position, we use rapidly moving mirrors to redirect our optical beam. Our collaboration has produced a unique MEMS (microelectromechanical systems) micromirror array which can point and settle on target in only 1 ms. This technology is proven, flying successfully as the MTEL (MEMS Telescope for Extreme Lightning) on the Tatiana-2 Spacecraft in September 2009 and as the KAMTEL on the International Space Station in April 2008. The sub-minute measurements of the UV-optical emission of dozens of GRB each year will result in a more rigorous test of current internal shock models, probe the extremes of bulk Lorentz factors, and provide the first early and detailed measurements of fast-rise and short type GRB UV-optical light curves.

The standard model of afterglow production by the forward shock wave is not supported by recent observations. We propose a model in which the forward shock is invisible and afterglow is emitted by a long-lived reverse shock in the burst ejecta. It explains observed optical and X-ray light curves, including the plateau at 10^3-10^4 s with a peculiar chromatic break, and the second break that was previously associated with a beaming angle of the explosion. The plateau forms following a temporary drop of the reverse-shock pressure much below the forward-shock pressure. A simplest formalism that can describe such blast waves is the ``mechanical'' model (Beloborodov, Uhm 2006); we use it in our calculations.

Relativistic blast waves can be described by a mechanical model. In this model, the "blast" -- the compressed gas between the forward and reverse shocks -- is viewed as one hot body. Equations governing its dynamics are derived from conservation of mass, energy, and momentum. Simple analytical solutions are obtained in the two limiting cases of ultra-relativistic and non-relativistic reverse shock. Equations are derived for the general explosion problem.

We propose a simple mechanical model for relativistic explosions with both forward and reverse shocks, which allows one to do fast calculations of GRB afterglow. The blast wave in the model is governed by pressures P_F and P_R at the forward and reverse shocks. We show that the simplest assumption P_F=P_R is in general inconsistent with energy conservation law. The model is applied to GRBs with non-uniform ejecta. Such "stratified fireballs" are likely to emerge with a monotonic velocity profile after an internal-shock stage. We calculate the early afterglow emission expected from stratified fireballs.