We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_{\nu}\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_{\nu}\propto \nu^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $\theta_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $\Gamma_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $\theta_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $\epsilon_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $\epsilon_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5\sigma$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (\Gamma\beta)^{-\alpha}$ into the environment, finding that shallow stratification indexes $\alpha\le6$ are disfavored.

We present a revised and complete optical afterglow light curve of the binary neutron star merger GW170817, enabled by deep Hubble Space Telescope (HST) F606W observations at $\approx\!584$ days post-merger, which provide a robust optical template. The light curve spans $\approx 110-362$ days, and is fully consistent with emission from a relativistic structured jet viewed off-axis, as previously indicated by radio and X-ray data. Combined with contemporaneous radio and X-ray observations, we find no spectral evolution, with a weighted average spectral index of $\langle \beta \rangle = -0.583 \pm 0.013$, demonstrating that no synchrotron break frequencies evolve between the radio and X-ray bands over these timescales. We find that an extrapolation of the post-peak temporal slope of GW170817 to the luminosities of cosmological short GRBs matches their observed jet break times, suggesting that their explosion properties are similar, and that the primary difference in GW170817 is viewing angle. Additionally, we place a deep limit on the luminosity and mass of an underlying globular cluster of $L \lesssim 6.7 \times 10^{3}\,L_{\odot}$, or $M \lesssim 1.3 \times 10^{4}\,M_{\odot}$, at least 4 standard deviations below the peak of the globular cluster mass function of the host galaxy, NGC4993. This limit provides a direct and strong constraint that GW170817 did not form and merge in a globular cluster. As highlighted here, HST (and soon JWST) enables critical observations of the optical emission from neutron star merger jets and outflows.

The field of time-domain astrophysics has entered the era of Multi-messenger Astronomy (MMA). One key science goal for the next decade (and beyond) will be to characterize gravitational wave (GW) and neutrino sources using the next generation of Extremely Large Telescopes (ELTs). These studies will have a broad impact across astrophysics, informing our knowledge of the production and enrichment history of the heaviest chemical elements, constrain the dense matter equation of state, provide independent constraints on cosmology, increase our understanding of particle acceleration in shocks and jets, and study the lives of black holes in the universe. Future GW detectors will greatly improve their sensitivity during the coming decade, as will near-infrared telescopes capable of independently finding kilonovae from neutron star mergers. However, the electromagnetic counterparts to high-frequency (LIGO/Virgo band) GW sources will be distant and faint and thus demand ELT capabilities for characterization. ELTs will be important and necessary contributors to an advanced and complete multi-messenger network.

We present a search for late-time rebrightening of radio emission from three supernovae (SNe) with associated gamma-ray bursts (GRBs). It has been previously proposed that the unusually energetic SNe associated with GRBs should enter the Sedov-Taylor phase decades after the stellar explosion, and this SN "remnant" emission will outshine the GRB radio afterglow and be detectable at significant distances. We place deep limits on the radio luminosity of GRB 980425/SN 1998bw, GRB 030329/SN 2003dh and GRB 060218/SN 2006aj, 10-18 years after explosion, with our deepest limit being $L_{\nu}$ $< 4 \times 10^{26}$ erg s$^{-1}$ Hz$^{-1}$ for GRB 980425/SN 1998bw. We put constraints on the density of the surrounding medium for various assumed values of the microphysical parameters related to the magnetic field and synchrotron-emitting electrons. For GRB 060218/SN 2006aj and GRB 980425/SN 1998bw, these density limits have implications for the density profile of the surrounding medium, while the non-detection of GRB 030329/SN 2003dh implies that its afterglow will not be detectable anymore at GHz frequencies.

The neutron star merger, GW170817, was followed by an optical-infrared transient (a kilonova) which indicated that a substantial ejection of mass at trans-relativistic velocities occurred during the merger. Modeling of the kilonova is able to constrain the kinetic energy of the ejecta and its characteristic velocity, but not the high-velocity distribution of the ejecta. Yet, this distribution contains crucial information on the merger dynamics. In this work, we assume a power-law distribution of the form $E(>\beta\Gamma)\propto(\beta\Gamma)^{-\alpha}$ for the energy of the kilonova ejecta and calculate the non-thermal signatures produced by the interaction of the ejecta with the ambient gas. We find that ejecta with minimum velocity $\beta_0\simeq 0.3$ and energy $E\sim 10^{51}$ erg, as inferred from kilonova modeling, has a detectable radio, and possibly X-ray, afterglow for a broad range of parameter space. This afterglow component is expected to dominate the observed emission on a timescale of a few years post merger and peak around a decade later. Its light curve can be used to determine properties of the kilonova ejecta and in particular the ejecta velocity distribution $\alpha$, the minimum velocity $\beta_0$ and its total kinetic energy $E$. We also predict that an afterglow rebrightening, that is associated with the kilonova component, will be accompanied by a shift of the centroid of the radio source towards the initial position of the explosion.

GW170817/GRB170817A has offered unprecedented insight into binary neutron star post-merger systems. Its Prompt and afterglow emission imply the presence of a tightly collimated relativistic jet with a smooth transverse structure. However, it remains unclear whether and how the central engine can produce such structured jets. Here, we utilize 3D GRMHD simulations starting with a black hole surrounded by a magnetized torus with properties typically expected of a post-merger system. We follow the jet, as it is self-consistently launched, from the scale of the compact object out to more than 3 orders of magnitude in distance. We find that this naturally results in a structured jet, which is collimated by the disk wind into a half-opening angle of roughly 10 degrees, its emission can explain features of both the prompt and afterglow emission of GRB170817A for a 30 degree observing angle. Our work is the first to compute the afterglow, in the context of a binary merger, from a relativistic magnetized jet self-consistently generated by an accreting black hole, with the jet's transverse structure determined by the accretion physics and not prescribed at any point.

We present new observations of the binary neutron star merger GW170817 at $\Delta t\approx 220-290$ days post-merger, at radio (Karl G. Jansky Very Large Array; VLA), X-ray (Chandra X-ray Observatory) and optical (Hubble Space Telescope; HST) wavelengths. These observations provide the first evidence for a turnover in the X-ray light curve, mirroring a decline in the radio emission at $\gtrsim5\sigma$ significance. The radio-to-X-ray spectral energy distribution exhibits no evolution into the declining phase. Our full multi-wavelength dataset is consistent with the predicted behavior of our previously published models of a successful structured jet expanding into a low-density circumbinary medium, but pure cocoon models with a choked jet cannot be ruled out. If future observations continue to track our predictions, we expect that the radio and X-ray emission will remain detectable until $\sim 1000$ days post-merger.

We performed calculations of the late radio and X-ray afterglow of GRB/GW170817 in the cocoon-jet paradigm, predicting the appearance of a second peak in the afterglow light curve ~ one-three years after the explosion. The model assumes that the prompt emission and early afterglows originate from a cocoon generated during break-out of the delayed magnetically powered jet. As the jet breaks out from the torus-generated wind, a nearly isotropic mildly relativistic outflow is generated; at the same time the primary jet accelerates to high Lorentz factors and avoids detection. As the fast jet slows down, it should become visible to the off-axis observer. Thus, the model has a clear prediction: the X-ray and radio afterglows should first experience a decay, as the cocoon slows down, followed by a rebrightening when the primary jet starts emitting toward an observer.

We report deep Chandra, HST and VLA observations of the binary neutron star event GW170817 at $t<160$ d after merger. These observations show that GW170817 has been steadily brightening with time and might have now reached its peak, and constrain the emission process as non-thermal synchrotron emission where the cooling frequency $\nu_c$ is above the X-ray band and the synchrotron frequency $\nu_m$ is below the radio band. The very simple power-law spectrum extending for eight orders of magnitude in frequency enables the most precise measurement of the index $p$ of the distribution of non-thermal relativistic electrons $N(\gamma)\propto \gamma^{-p}$ accelerated by a shock launched by a NS-NS merger to date. We find $p=2.17\pm0.01$, which indicates that radiation from ejecta with $\Gamma\sim3-10$ dominates the observed emission. While constraining the nature of the emission process, these observations do \emph{not} constrain the nature of the relativistic ejecta. We employ simulations of explosive outflows launched in NS ejecta clouds to show that the spectral and temporal evolution of the non-thermal emission from GW170817 is consistent with both emission from radially stratified quasi-spherical ejecta traveling at mildly relativistic speeds, \emph{and} emission from off-axis collimated ejecta characterized by a narrow cone of ultra-relativistic material with slower wings extending to larger angles. In the latter scenario, GW170817 harbored a normal SGRB directed away from our line of sight. Observations at $t\le 200$ days are unlikely to settle the debate as in both scenarios the observed emission is effectively dominated by radiation from mildly relativistic material.

Binary neutron star mergers are considered to be the most favorable sources that produce electromagnetic (EM) signals associated with gravitational waves (GWs). These mergers are the likely progenitors of short duration gamma-ray bursts (GRBs). The brief gamma-ray emission (the "prompt GRB" emission) is produced by ultra-relativistic jets, as a result, this emission is strongly beamed over a small solid angle along the jet. It is estimated to be a decade or more before a short GRB jet within the LIGO volume points along our line of sight. For this reason, the study of the prompt signal as an EM counterpart to GW events has been sparse. We argue that for a realistic jet model, one whose luminosity and Lorentz factor vary smoothly with angle, the prompt signal can be detected for a significantly broader range of viewing angles. This can lead to a new type of EM counterpart, an "off-axis" short GRB. Our estimates and simulations show that it is feasible to detect these signals with the aid of the temporal coincidence from a LIGO trigger, even if the observer is substantially misaligned with respect to the jet.

Numerous tidal disruption event (TDE) candidates originating from galactic centres have been detected (e.g., by ${\it Swift}$ and ASASSN). Some of their host galaxies show typical characteristics of a weak active galactic nucleus (AGN), indicative of a pre-existing accretion disc around the supermassive black hole (SMBH). In this work, we develop an analytic model to study how a pre-existing accretion disc affects a TDE. We assume the density of the disc $\rho\propto R^{-\lambda}$, $R$ being the radial distance from the SMBH and $\lambda$ varying between $0.5$ and $1.5$. Interactions between the pre-existing accretion disc and the stream of the tidally disrupted star can stall the stream far from the SMBH, causing a sudden drop in the rate of fallback of gas into the SMBH. These interactions could explain the steep cut-off observed in the light curve of some TDE candidates (e.g., ${\it Swift}$ J1644 and ${\it Swift}$ J2058). With our model, it is possible to use the time of this cut-off to constrain some properties pertaining to the pre-existing accretion disc, such as $\lambda$ and the disc viscosity parameter $\alpha$. We demonstrate this by applying our theory to the TDE candidates ${\it Swift}$ J1644, ${\it Swift}$ J2058 and ASASSN-14li. Our analysis favours a disc profile with $\lambda\sim1$ for viscosity parameters $\alpha\sim0.01-0.1$.

Gamma-ray burst (GRB) afterglows are likely produced in the shock that is driven as the GRB jet interacts with the external medium. Long-duration GRBs are also associated with powerful supernovae (SNe). We consider the optical and radio afterglows of long GRBs for both blasts viewed along the jet axis ("on-axis" afterglows) and misaligned observes ("off-axis" afterglows). Comparing the optical emission from the afterglow with that of the accompanying SN, using SN 1998bw as an archetype, we find that only a few percent of afterglows viewed off-axis are brighter than the SN. For observable optical off-axis afterglows, the viewing angle is at most twice the half-opening angle of the GRB jet. Radio off-axis afterglows should be detected with upcoming radio surveys within a few hundred Mpc. We propose that these surveys will act as "radio triggers", and that dedicated radio facilities should follow-up these sources. Follow-ups can unveil the presence of the radio SN remnant, if present. In addition, they can probe the presence of a mildly relativistic component, either associated with the GRB jet or the SN ejecta, expected in these sources.