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.

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.

GRB 160821B is a short duration gamma-ray burst (GRB) detected and localized by the Neil Gehrels Swift Observatory in the outskirts of a spiral galaxy at z=0.1613, at a projected physical offset of 16 kpc from the galaxy's center. We present X-ray, optical/nIR and radio observations of its counterpart and model them with two distinct components of emission: a standard afterglow, arising from the interaction of the relativistic jet with the surrounding medium, and a kilonova, powered by the radioactive decay of the sub-relativistic ejecta. Broadband modeling of the afterglow data reveals a weak reverse shock propagating backward into the jet, and a likely jet-break at 3.5 d. This is consistent with a structured jet seen slightly off-axis while expanding into a low-density medium. Analysis of the kilonova properties suggests a rapid evolution toward red colors, similar to AT2017gfo, and a low nIR luminosity, possibly due to the presence of a long-lived neutron star. The global properties of the environment, the inferred low mass (M_ej < 0.006 Msun) and velocities (v > 0.05 c) of lanthanide-rich ejecta are consistent with a binary neutron star merger progenitor.

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.

Most violent and energetic processes in our universe, including mergers of compact objects, explosions of massive stars and extreme accretion events, produce copious amounts of X-rays. X-ray follow-up is an efficient tool for identifying transients because (1) X-rays can quickly localize transients with large error circles, and (2) X-rays reveal the nature of transients that may not have unique signatures at other wavelengths. In this white paper, we identify key science questions about several extragalactic multi-messenger and multi-wavelength transients, and demonstrate how X-ray follow-up helps answer these questions

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.

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.

Newtonian simulations have demonstrated that accretion onto binary black holes produces accretion disks around each black hole ("minidisks"), fed by gas streams flowing through the circumbinary cavity from the surrounding circumbinary disk. We study the dynamics and radiation of an individual black hole minidisk using two-dimensional hydrodynamical simulations performed with a new general relativistic version of the moving mesh code Disco. We introduce a co-moving energy variable which enables highly accurate integration of these high Mach number flows. Tidally induced spiral shock waves are excited in the disk and propagate through the ISCO providing a Reynolds stress which causes efficient accretion by purely hydrodynamic means and producing a radiative signature brighter in hard X-rays than the Novikov-Thorne model. Disk cooling is provided by a local blackbody prescription that allows the disk to evolve self-consistently to a temperature profile where hydrodynamic heating is balanced by radiative cooling. We find that the spiral shock structure is in agreement with the relativistic dispersion relation for tightly-wound linear waves. We measure the shock induced dissipation and find outward angular momentum transport corresponding to an effective alpha parameter of order 0.01. We perform ray-tracing image calculations from the simulations to produce theoretical minidisk spectra and viewing angle dependent images for comparison with observations.

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

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.