We present the afterglow light curves produced by the deceleration of the non-relativistic ejecta mass in a stratified circumstellar medium with a density profile $n(r)\propto r^{-k}$ with $k=0$, $1$, $1.5$, $2$ and $2.5$. Once the ejecta mass is launched with equivalent kinetic energy parametrized by $E(>\beta)\propto \beta^{-\alpha}$ (where beta is the ejecta velocity) and propagates into the surrounding circumstellar medium, it first moves with constant velocity (the free-coasting phase), and later it decelerates (the Sedov-Taylor expansion). We present the predicted synchrotron and synchrotron-self Compton light curves during the free-coasting phase, and the subsequent Sedov-Taylor expansion. In particular cases, we show the corresponding light curves generated by the deceleration of several ejecta masses with different velocities launched during the coalescence of binary compact objects and the core-collapse of dying massive stars which will contribute at distinct timescales, frequencies, and intensities. We find that before the radiation escapes from the kilonova (KN) and produces a peak in the light curve, the KN afterglow ejecta is very weak and the corresponding ejecta velocity moves as $\propto t^{-\frac{2}{\alpha+1}}$. Finally, using the multi-wavelength observations and upper limits collected by a large campaign of orbiting satellites and ground telescopes, we constrain the parameter space of both the KN afterglow in GW170817 and the possibly generated KN afterglow in S190814bv. Further observations on timescales of years post-merger are needed to derive tighter constraints.

The detection of gravitational waves together with their electromagnetic counterpart, in the gamma-ray burst GRB 170817A, marked a new era of multi-messenger astronomy. Several theoretical models have been proposed to explain the atypical behavior of this event. Recently, it was shown that the multi-wavelength afterglow of GRB 170817A was consistent with a synchrotron forward-shock model when the outflow was viewed off-axis, decelerated in a uniform medium and parametrized through a power-law velocity distribution. Motivated by the upper limits on the very-high-energy emission, and the stratified medium in the close vicinity of a binary neutron star merger proposed to explain the gamma-ray flux in the short GRB 150101B, we extend the mechanism proposed to explain GRB 170817A to a more general scenario deriving the synchrotron self-Compton (SSC) and synchrotron forward-shock model when the off-axis outflow is decelerated in a uniform and stratified circumburst density. As particular cases, we show that the delayed and long-lasting afterglow emission observed in GRB 080503, GRB140903A, GRB 150101B, and GRB 160821B could be interpreted by a similar scenario to the one used to describe GRB 170817A. In addition, we show that the proposed scenario agrees with the MAGIC, Fermi-LAT and H.E.S.S upper limits on gamma-ray emission from GRB 160821B and GRB 170817A.

Early and late multiwavelength observations play an important role in determining the nature of the progenitor, circumburst medium, physical processes and emitting regions associated to the spectral and temporal features of bursts. GRB 180720B is a long and powerful burst detected by a large number of observatories in multiwavelenths that range from radio bands to sub-TeV gamma-rays. The simultaneous multiwavelength observations were presented over multiple periods of time beginning just after the trigger time and extending for more than 30 days. The temporal and spectral analysis of Fermi LAT observations suggests that it presents similar characteristics to other bursts detected by this instrument. Coupled with X-ray and optical observations, the standard external-shock model in a homogeneous medium is favored by this analysis. The X-ray flare is consistent with the synchrotron self-Compton (SSC) model from the reverse-shock region evolving in a thin shell and long-lived LAT, X-ray and optical data with the standard synchrotron forward-shock model. The best-fit parameters derived with the Markov chain Monte Carlo simulations indicate that the outflow is endowed with magnetic fields and that the radio observations are in the self-absorption regime. The SSC forward-shock model with our parameters can explain the LAT photons beyond the synchrotron limit as well as the emission recently reported by the HESS Collaboration.

The production of both gravitational waves and short gamma-ray bursts (sGRBs) is widely associated with the merger of compact objects. Several studies have modelled the evolution of the electromagnetic emission using the synchrotron emission produced by the deceleration of both a relativistic top-hat jet seen off-axis, and a wide-angle quasi-spherical outflow (both using numerical studies). In this study we present an analytical model of the synchrotron and synchrotron self-Compton (SSC) emission for an off-axis top-hat jet and a quasi-spherical outflow. We calculate the light curves obtained from an analytic model in which the synchrotron and SSC emission (in the fast- or slow-cooling regime) of an off-axis top-hat jet and a quasi-spherical outflow are decelerated in either a homogeneous or a wind-like circumburst medium. We show that the synchrotron emission of the quasi-spherical outflow is stronger than that of the off-axis jet during the first $\sim$ 10 - 20 days, and weaker during the next $\gtrsim$ 80 days. Moreover, we show that if the off-axis jet is decelerated in a wind-like medium, then the SSC emission is very likely to be detected. Applying a MCMC code to our model (for synchrotron emission only), we find the best-fit values for the radio, optical and X-ray emission of GRB 170817A which are in accordance with observations. For GRB 170817A, we find using our model that the synchrotron emission generated by the quasi-spherical outflow and off-axis top-hat jet increase as $F_\nu\propto t^{\alpha}$ with $\alpha\lesssim 0.8$ and $\alpha>3$, respectively. Finally, we obtain the correspondent SSC light curves which are in accordance with the very-high-energy gamma-ray upper limits derived with the GeV - TeV observatories.