The afterglow emission from gamma-ray bursts (GRBs) is believed to originate from a relativistic blast wave driven into the circumburst medium. Although the afterglow emission from radio up to X-ray frequencies is thought to originate from synchrotron radiation emitted by relativistic, non-thermal electrons accelerated by the blast wave, the origin of the emission at high energies (HE; $\gtrsim$~GeV) remains uncertain. The recent detection of sub-TeV emission from GRB~190114C by MAGIC raises further debate on what powers the very high-energy (VHE; $\gtrsim 300$GeV) emission. Here, we explore the inverse Compton scenario as a candidate for the HE and VHE emissions, considering two sources of seed photons for scattering: synchrotron photons from the blast wave (synchrotron self-Compton or SSC) and isotropic photon fields external to the blast wave (external Compton). For each case, we compute the multi-wavelength afterglow spectra and light curves. We find that SSC will dominate particle cooling and the GeV emission, unless a dense ambient infrared photon field, typical of star-forming regions, is present. Additionally, considering the extragalactic background light attenuation, we discuss the detectability of VHE afterglows by existing and future gamma-ray instruments for a wide range of model parameters. Studying GRB~190114C, we find that its afterglow emission in the \fermi-LAT band is synchrotron-dominated.The late-time \fermi-LAT measurement (i.e., $t\sim 10^4$~s), and the MAGIC observation also set an upper limit on the energy density of a putative external infrared photon field (i.e. $\lesssim 3\times 10^{-9}\,{\rm erg\,cm^{-3}}$), making the inverse Compton dominant in the sub-TeV energies.

The central engine in long gamma-ray bursts (GRBs) is thought to be a compact object produced by the core collapse of massive stars, but its exact nature (black hole or millisecond magnetar) is still debatable. Although the central engine of GRB collapsars is hidden to direct observation, its properties may be imprinted on the accompanying electromagnetic signals. We aim to decipher the generic properties of central engines that are consistent with prompt observations of long GRBs detected by the Burst Alert Telescope (BAT) on board the Neil Gehrels Swift Observatory. Adopting a generic model for the central engine, in which the engine power and activity timescale are independent of each other, we perform Monte Carlo simulations of long GRBs produced by jets that successfully breakout from the star. Our simulations consider the dependence of the jet breakout timescale on the engine luminosity and the effects of the detector's flux threshold. The two-dimensional (2D) distribution of simulated detectable bursts in the gamma-ray luminosity versus gamma-ray duration plane is consistent with the observed one for a range of parameter values describing the central engine. The intrinsic 2D distribution of simulated collapsar GRBs peaks at lower gamma-ray luminosities and longer durations than the observed one, a prediction that can be tested in the future with more sensitive detectors. Black-hole accretors, whose power and activity time are set by the large-scale magnetic flux through the progenitor star and stellar structure, respectively, are compatible with the properties of the central engine inferred by our model.

The joint detection of GW~170817/GRB 170817 confirmed the long-standing theory that binary neutron star mergers produce short gamma-ray burst (sGRB) jets that can successfully break out of the surrounding ejecta. At the same time, the association with a kilonova provided unprecedented information regarding the physical properties (such as masses and velocities) of the different ejecta constituents. Combining this knowledge with the observed luminosities and durations of cosmological sGRBs detected by the Burst Alert Telescope (BAT) onboard the Neil Gehrels Swift Observatory, we revisit the breakout conditions of sGRB jets. Assuming self-collimation of sGRB jets does not play a critical role, we find that the time interval between the binary merger and the launching of a typical sGRB jet is $\lesssim0.1$~s. We also show that for a fraction of at least $\sim 30\%$ of sGRBs, the usually adopted assumption of static ejecta is inconsistent with observations, even if the polar ejecta mass is an order of magnitude smaller than the one in GRB 170817. Our results disfavour magnetar central engines for powering cosmological sGRBs, limit the amount of energy deposited in the cocoon prior to breakout, and suggest that the observed delay of $\sim 1.$7~s in GW 170817 /GRB 170817 between the gravitational wave and $\gamma$-ray signals is likely dominated by the propagation time of the jet to the $\gamma$-ray production site.

Periodicities observed in two Fast Radio Burst (FRB) sources (16 days in FRB 180916.J0158+65 and 160 days in FRB 121102) are consistent with that of tight, stellar mass binary systems. In the case of FRB 180916.J0158+65 the primary is an early OB-type star with mass loss rate $\dot{M} \sim 10^{-8}- 10^{-7} M_\odot$ yr$^{-1}$, and the secondary a neutron star. The observed periodicity is not intrinsic to the FRB's source, but is due to the orbital phase-dependent modulation of the absorption conditions in the massive star's wind. The observed relatively narrow FRB activity window implies that the primary's wind dynamically dominates that of the pulsar, $\eta = L_{sd}/(\dot{M} v_w c) \leq 1$, where $L_{sd} $ is pulsar spin-down, $\dot{M}$ is the primary's wind mass loss rate and $v_w$ is its velocity. The condition $\eta \leq 1$ requires mildly powerful pulsar with $L_{sd} \lesssim 10^{37}$ erg $s^{-1}$. The observations are consistent with magnetically-powered radio emission originating in the magnetospheres of strongly magnetized neutron stars, the classical magnetars.

Relativistic jets are highly collimated plasma outflows emerging from accreting black holes. They are launched with a significant amount of magnetic energy, which can be dissipated to accelerate nonthermal particles and give rise to electromagnetic radiation at larger scales. Kink instabilities can be an efficient mechanism to trigger dissipation of jet magnetic energy. While previous works have studied the conditions required for the growth of kink instabilities in relativistic jets, the radiation signatures of these instabilities have not been investigated in detail. In this paper, we aim to self-consistently study radiation and polarization signatures from kink instabilities in relativistic jets. We combine large-scale relativistic magnetohydrodynamic (RMHD) simulations with polarized radiation transfer of a magnetized jet, which emerges from the central engine and propagates through the surrounding medium. We observe that a localized region at the central spine of the jet exhibits the strongest kink instabilities, which we identify as the jet emission region. Very interestingly, we find quasi-periodic oscillation (QPO) signatures in the light curve from the emission region. Additionally, the polarization degree appears to be anti-correlated to flares in the light curves. Our analyses show that these QPO signatures are intrinsically driven by kink instabilities, where the period of the QPOs is associated to the kink growth time scale. The latter corresponds to weeks to months QPOs in blazars. The polarization signatures offer unique diagnostics for QPOs driven by kink instabilities.

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.

Blazar emission models based on magnetic reconnection succeed in reproducing many observed spectral and temporal features, including the short-duration luminous flaring events. Plasmoids, a self-consistent by-product of the tearing instability in the reconnection layer, can be the main source of blazar emission. Kinetic simulations of relativistic reconnection have demonstrated that plasmoids are characterized by rough energy equipartition between their radiating particles and magnetic fields. This is the main reason behind the apparent shortcoming of plasmoid-dominated emission models to explain the observed Compton ratios of BL Lac objects. Here, we demonstrate that the radiative interactions among plasmoids, which have been neglected so far, can assist in alleviating this contradiction. We show that photons emitted by large, slow-moving plasmoids can be a potentially important source of soft photons to be then up-scattered, via inverse Compton, by small fast-moving, neighboring plasmoids. This inter-plasmoid Compton scattering process can naturally occur throughout the reconnection layer, imprinting itself as an increase in the observed Compton ratios from those short and luminous plasmoid-powered flares within BL Lac sources, while maintaining energy equipartition between radiating particles and magnetic fields.

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.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

Magnetic reconnection is often invoked to explain the non-thermal radiation of relativistic outflows, including jets of active galactic nuclei (AGN). Motivated by the largely unknown plasma composition of AGN jets, we study reconnection in the unexplored regime of electron-positron-proton (pair-proton) plasmas with large-scale two-dimensional particle-in-cell simulations. We cover a wide range of pair multiplicities (lepton-to-proton number ratio $\kappa=1-199$) for different values of the all-species plasma magnetization ($\sigma=1,3$ and 10) and electron temperature ($\Theta_e\equiv kT_e/m_ec^2=0.1-100$). We focus on the dependence of the post-reconnection energy partition and lepton energy spectra on the hot pair plasma magnetization $\sigma_{e,h}$ (i.e., the ratio of magnetic to pair enthalpy densities). We find that the post-reconnection energy is shared roughly equally between magnetic fields, pairs, and protons for $\sigma_{e,h}\gtrsim 3$. We empirically find that the mean lepton Lorentz factor in the post-reconnection region depends on $\sigma, \Theta_e$, and $\sigma_{e,h}$ as $\langle \gamma_e-1\rangle \approx \sqrt{\sigma}(1+4\Theta_e)\left(1+\sigma_{e,h}/30\right)$, for $\sigma\ge1$. The high-energy part of the post-reconnection lepton energy distributions can be described by a power law, whose slope is mainly controlled by $\sigma_{e,h}$ for $\kappa \gtrsim 3-6$, with harder power laws obtained for higher magnetizations. We finally show that reconnection in pair-proton plasmas with multiplicities $\kappa \sim 1-20$, magnetizations $\sigma \sim 1-10$, and temperatures $\Theta_e \sim 1-10$ results in particle power law slopes and average electron Lorentz factors that are consistent with those inferred in leptonic models of AGN jet emission.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

The possible association of the blazar TXS 0506+056 with a high-energy neutrino detected by IceCube holds the tantalizing potential to answer three astrophysical questions: 1. Where do high-energy neutrinos originate? 2. Where are cosmic rays produced and accelerated? 3. What radiation mechanisms produce the high-energy {\gamma}-rays in blazars? The MeV gamma-ray band holds the key to these questions, because it is an excellent proxy for photo-hadronic processes in blazar jets, which also produce neutrino counterparts. Variability in MeV gamma-rays sheds light on the physical conditions and mechanisms that take place in the particle acceleration sites in blazar jets. In addition, hadronic blazar models also predict a high level of polarization fraction in the MeV band, which can unambiguously distinguish the radiation mechanism. Future MeV missions with a large field of view, high sensitivity, and polarization capabilities will play a central role in multi-messenger astronomy, since pointed, high-resolution telescopes will follow neutrino alerts only when triggered by an all-sky instrument.

The constantly improving sensitivity of ground-based and space-borne observatories has made possible the detection of high-energy emission (X-rays and gamma-rays) from several thousands of extragalactic sources. Enormous progress has been made in measuring the continuum flux enabling us to perform imaging, spectral and timing studies. An important remaining challenge for high-energy astronomy is measuring polarization. The capability to measure polarization is being realized currently at X-ray energies (e.g. with IXPE), and sensitive gamma-ray telescopes capable of measuring polarization, such as AMEGO, AdEPT, e-ASTROGAM, etc., are being developed. These future gamma-ray telescopes will probe the radiation mechanisms and magnetic fields of relativistic jets from active galactic nuclei at spatial scales much smaller than the angular resolution achieved with continuum observations of the instrument. In this white paper, we discuss the scientific potentials of high-energy polarimetry, especially gamma-ray polarimetry, including the theoretical implications, and observational technology advances being made. In particular, we will explore the primary scientific opportunities and wealth of information expected from synergy of multi-wavelength polarimetry that will be brought to multi-messenger astronomy.

Active galactic nuclei (AGN) with relativistic jets, powered by gas accretion onto their central supermassive black hole (SMBH), are unique laboratories for studying the physics of matter and elementary particles in extreme conditions that cannot be realized on Earth. For a long time since the discovery of AGN, photons were the only way to probe the underlying physical processes. The recent discovery of a very high energy neutrino, IceCube-170922A, coincident with a flaring blazar, TXS 0506+056, provides the first evidence that AGN jets are multi-messenger sources; they are capable of accelerating hadrons to very high energies, while producing non-thermal EM radiation and high-energy neutrinos. This new era of multi-messenger astronomy, which will mature in the next decade, offers us the unprecedented opportunity to combine more than one messenger to solve some long-standing puzzles of AGN jet physics: How do jets dissipate their energy to accelerate particles? What is the jet total kinetic power? Where and how do jets produce high-energy emission and neutrinos? What physical mechanisms drive the particle acceleration?

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.

The characteristic two-component blazar spectral energy distribution (SED) can be of either leptonic and/or hadronic origins. The potential association of the high-energy neutrino event IceCube-170922A with the flaring blazar TXS~0506+056 indicates that hadronic processes may operate in a blazar jet. Despite multi-wavelength follow-ups of the event and extensive theoretical modelings, the radiation mechanisms and the underlying magnetic field strength and configuration remain poorly understood. In this paper, we consider generic leptonic and hadronic blazar spectral models with distinct magnetic field strengths and radiation mechanisms. We analytically reproduce the SEDs and the neutrino flux of hadronic models, and predict their X-ray to $\gamma$-ray polarization degrees. Furthermore, by performing relativistic magnetohydrodynamic (RMHD) simulations taking into account the polarization-dependent radiation transfer, we study the time-dependent multi-wavelength polarization variability of the proton synchrotron model under a shock scenario. Our results suggest that the high-energy polarization degree and the neutrino flux can be jointly used to pinpoint the leptonic and/or hadronic blazar radiation mechanisms in the X-ray and $\gamma$-ray bands, and to infer the magnetic field strength in the emission region. Additionally, the temporal multi-wavelength polarization signatures in the proton synchrotron model shed light on the jet energy composition and the dynamical importance of magnetic fields in the blazar emission region. Future multi-wavelength polarimetry facilities such as {\it IXPE} and {\it AMEGO} together with neutrino telescopes such as {\it IceCube} can provide unprecedented observational constraints to probe the blazar radiation mechanisms and jet dynamics.

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.

The joint detection of gravitational waves (GWs) and $\gamma$-rays from a binary neutron-star (NS) merger provided a unique view of off-axis GRBs and an independent measurement of the NS merger rate. Comparing the observations of GRB170817 with those of the regular population of short GRBs (sGRBs), we show that an order unity fraction of NS mergers result in sGRB jets that breakout of the surrounding ejecta. We argue that the luminosity function of sGRBs, peaking at $\approx 2\times 10^{52}\, \mbox{erg s}^{-1}$, is likely an intrinsic property of the sGRB central engine and that sGRB jets are typically narrow with opening angles $\theta_0 \approx 0.1$. We perform Monte Carlo simulations to examine models for the structure and efficiency of the prompt emission in off axis sGRBs. We find that only a small fraction ($\sim 0.01-0.1$) of NS mergers detectable by LIGO/VIRGO in GWs is expected to be also detected in prompt $\gamma$-rays and that GW170817-like events are very rare. For a NS merger rate of $\sim 1500$ Gpc$^{-3}$ yr$^{-1}$, as inferred from GW170817, we expect within the next decade up to $\sim 12$ joint detections with off-axis GRBs for structured jet models and just $\sim 1$ for quasi-spherical cocoon models where $\gamma$-rays are the result of shock breakout. Given several joint detections and the rates of their discoveries, the different structure models can be distinguished. In addition the existence of a cocoon with a reservoir of thermal energy may be observed directly in the UV, given a sufficiently rapid localisation of the GW source.

The multi-wavelength spectral and temporal variability observed in blazars set tight constraints on current theoretical emission models. Here, we investigate the relativistic magnetic reconnection process as a source of blazar emission in which quasi-spherical plasmoids, containing relativistic particles and magnetic fields, are associated with the emission sites in blazar jets. By coupling recent two-dimensional particle-in-cell simulations of relativistic reconnection with a time-dependent radiative transfer code, we compute the non-thermal emission from a chain of plasmoids formed during a reconnection event. The derived photon spectra display characteristic features observed in both BL Lac sources and flat spectrum radio quasars, with the distinction made by varying the strength of the external photon fields, the jet magnetization, and the number of pairs per proton contained within. Light curves produced from reconnection events are composed of many fast and powerful flares that appear on excess of a slower evolving envelope produced by the cumulative emission of medium-sized plasmoids. The observed variability is highly dependent upon the orientation of the reconnection layer with respect to the blazar jet axis and to the observer. Our model provides a physically motivated framework for explaining the multi-timescale blazar variability across the entire electromagnetic spectrum.

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.

Gamma-ray bursts are the most energetic electromagnetic sources in the Universe. Their prompt gamma-ray radiation corresponds to an energy release of 1E42-1E47J. Fifty years after their discovery and several dedicated space-based instruments, the physical origin of this emission is still unknown. Synchrotron emission has been one of the early contenders but was criticized because spectral fits of empirical models (such as a smoothly-connected broken power law or a cut-off power law) suggest too hard a slope of the low-energy power law, violating the so-called synchrotron line-of-death. We perform time-resolved gamma-ray spectroscopy of single-peaked GRBs as measured with Fermi/GBM. We demonstrate that idealized synchrotron emission, when properly incorporating time-dependent cooling of the electrons, is capable of fitting ~95% of all these GBM spectra. The comparison with spectral fit results based on previous empirical models demonstrates that the past exclusion of synchrotron radiation as an emission mechanism derived via the line-of-death was misleading. Our analysis probes the physics of these ultra-relativistic outflows and the related microphysical processes, and for the first time provides estimates of magnetic field strength and Lorentz factors of the emitting region directly from spectral fits. Our modeling of the Fermi/GBM observations provides evidence that GRBs are produced by moderately magnetized jets in which relativistic mini-jets emit optically-thin synchrotron radiation at large emission radii.

Frequency-dependent brightness fluctuations of radio sources, the so-called extreme scattering events (ESEs), have been observed over the last three decades. They are caused by Galactic plasma structures whose geometry and origin are still poorly understood. In this paper, we construct axisymmentric two-dimensional (2D) column density profiles for the plasma lens and explore the resulting ESEs for both point-like and extended sources. A quantity that becomes relevant is the impact parameter $b$, namely the distance of the observer's path from the lens symmetry axis. We demonstrate its effects on the shape of ESE light curves and use it for a phenomenological classification of ESEs into four main types. Three of them are unique outcomes of the 2D model and do not show a characteristic U-shaped dip in the light curve, which has been traditionally used as an identification means of ESEs. We apply our model to five well-studied ESEs and show that elongated plasma tubes or quasi-spherical clouds are favoured over plasma sheets for four of them, while the remaining one is compatible with both lens geometries.

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.

Relativistic magnetic reconnection events can widely exist in magnetized plasmas in astrophysical systems. During this process, oppositely directed magnetic field lines reconnect and release magnetic energy, efficiently accelerating nonthermal particles. However, so far there is little clear observational signatures of relativistic magnetic reconnection events in astrophysical systems. Blazars are relativistic magnetized plasma outflows from supermassive black holes. Their multi-wavelength flares may be powered by relativistic magnetic reconnection. The highly variable radiation and polarization signatures are well covered by multi-wavelength observation campaigns, making them ideal targets to examine the magnetic reconnection model. Recent observations have found that several blazar flares are accompanied by optical polarization angle swings which may be of as large amplitude as $> 180^{\circ}$, challenging existing theoretical models. In this paper, we present integrated particle-in-cell (PIC) and polarized radiation transfer simulations of magnetic reconnection events. We find that plasmoid coalescences in the reconnection layer can give rise to highly variable light curves, low and fluctuating polarization degree, and rotating polarization angle. In particular, large-amplitude polarization angle swings, similar to those observed during blazar flares, can be a unique signature of relativistic magnetic reconnection events.

Black-hole driven relativistic astrophysical jets, such as blazars and gamma-ray bursts (GRBs), are powerful sources of electromagnetic radiation. Their emission is powered by some energy dissipation and particle acceleration mechanism operating in-situ at large distances from the black hole. We propose that the formation of the dissipative structures in the jet is controlled by the time variability of the accretion disc. We argue that the open magnetic field lines through the black hole, which drive a strongly magnetized jet, may have their polarity reversing over time-scales related to the growth of the magneto-rotational dynamo in the disc. Consequently, the jet is characterized by an alternating toroidal field polarity along its propagation axis, i.e., it is a "striped jet". Magnetic reconnection in the current sheets that form between the stripes dissipates the alternating-field energy and thus powers further jet acceleration. Here, we consider a jet with a broad distribution of stripe widths $l>l_{\rm min}$, above a dominant scale $l_{\rm min}$. We find that the bulk acceleration of the jet, driven by the annihilation of the stripes, is very gradual. The dissipation rate peaks at a distance $z_{\rm peak}\sim 10^6 R_g\, (\Gamma_\infty/30)^2\, (l_{\rm min}/1000R_g)$, where $R_{g}$ is the black-hole's gravititional radius and $\Gamma_{\infty}$ the jet's asymptotic Lorentz factor, and exhibits a very broad plateau extending by $\sim 4-5$ orders of magnitude in distance. The prolonged energy dissipation accounts for the flat-to-inverted long-wavelength spectra commonly observed in jets. The model can also account for the broad range of flaring timescales of blazars and the fact that their bulk acceleration appears to continue out to ~100 pc scales. In GRB jets, the model predicts comparable power for the photopsheric and Thomson-thin emission components.

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.

The recent discovery of gravitational waves from GW170817, associated with a short Gamma-Ray Burst (GRB) at a distance of $40$Mpc, has demonstrated that short GRBs can occur locally and at a reasonable rate. Furthermore, gravitational waves enable us to detect close by GRBs, even when we are observing at latitudes far from the jet's axis. We consider here Compton echoes, the scattered light from the prompt and afterglow emission. Compton echoes, an as yet undetected counterpart of GRBs, peak in X-rays and maintain a roughly constant flux for hundreds to thousands of years after the burst. Though too faint to be detected in typical cosmological GRBs, a fraction of close by bursts with a sufficiently large energy output in X-rays, and for which the surrounding medium is sufficiently dense, may indeed be observed in this way. The detection of a Compton echo could provide unique insight into the burst properties and the environment's density structure. In particular, it could potentially determine whether or not there was a successful jet that broke through the compact binary merger ejecta. We discuss here the properties and expectations from Compton echoes and suggest methods for detectability.

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.

Rapidly spinning, strongly magnetized proto-neutron stars ("millisecond proto-magnetars") are candidate central engines of long-duration gamma-ray bursts (GRB), superluminous supernovae (SLSNe), and binary neutron star mergers. Magnetar birth may be accompanied by the fall-back of stellar debris, lasting for seconds or longer following the explosion. Accretion alters the magnetar evolution by (1) providing an additional source of rotational energy (or a potential sink, if the propeller mechanism operates); (2) enhancing the spin-down luminosity above the dipole rate by compressing the magnetosphere and expanding the polar cap region of open magnetic field lines; (3) supplying an additional accretion-powered neutrino luminosity that sustains the wind baryon-loading, even after the magnetar's internal neutrino luminosity has subsided. The more complex evolution of the jet power and magnetization of an accreting magnetar more readily accounts for the high 56Ni yields GRB SNe and irregular time evolution of some GRB light curves (e.g.~bursts with precursors followed by a long quiescent interval before the main emission episode). Additional baryon-loading from accretion-powered neutrino irradiation of the magnetar polar cap lengthens the timeframe over which the jet magnetization is in the requisite range sigma ~< 1e3 for efficient gamma-ray emission, thereby accommodating GRBs with ultra-long durations. Though accretion does not significantly raise the maximum energy budget from the limit of <~ few 1e52 ergs for an isolated magnetar, it greatly expands the range of magnetic field strengths and birth spin periods capable of powering GRB jets, reducing the differences between the magnetar properties normally invoked to explain GRBs versus SLSNe.

Previous studies have considered synchrotron as the emission mechanism for prompt Gamma-Ray Bursts (GRBs). These works have shown that the electrons must cool on a timescale comparable to the dynamic time at the source in order to satisfy spectral constraints while maintaining high radiative efficiency. We focus on conditions where synchrotron cooling is balanced by a continuous source of heating, and in which these constraints are naturally satisfied. Assuming that a majority of the electrons in the emitting region are contributing to the observed peak, we find that the energy per electron has to be $E\gtrsim 20$ GeV and that the Lorentz factor of the emitting material has to be very large $10^3\lesssim \Gamma_{\rm em} \lesssim 10^4$, well in excess of the bulk Lorentz factor of the jet inferred from GRB afterglows. A number of independent constraints then indicate that the emitters must be moving relativistically, with $\Gamma'\approx 10$, relative to the bulk frame of the jet and that the jet must be highly magnetized upstream of the emission region, $\sigma_{\rm up}\gtrsim 30$. The emission radius is also strongly constrained in this model to $R\gtrsim 10^{16}$cm. These values are consistent with magnetic jet models where the dissipation is driven by magnetic reconnection that takes place far away from the base of the jet.

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.

The recent coincident detection of gravitational waves (GW) from a binary neutron star merger with aLIGO/Virgo and short-lived gamma-ray emission with Fermi/GBM (called GW 170817) is a milestone for the establishment of multi-messenger astronomy. Merging neutron stars (NS) represent the standard scenario for short-duration (< 2 sec) gamma-ray bursts (GRBs) which are produced in a collimated, relativistically expanding jet with an opening angle of a few degrees and a bulk Lorentz factor of 300-1000. While the present aLIGO detection is consistent with predictions, the measured faint gamma-ray emission from GW 170817A, if associated to the merger event at a distance of 40 Mpc, is about 1000x less luminous than known short-duration GRBs (sGRBs). Hence, the presence of this sGRB in the local Universe is either a very rare event, or points to a dramatic ignorance of the emission properties of sGRBs outside their narrow jets. Here we show that the majority of previously detected faint sGRBs are local, at redshift smaller than 0.1, seen off-axis. In contrast, the brighter sGRBs are seen on-axis, and therefore out to larger distances, consistent with the measured redshift distribution. Examining the observer-frame parameter space of all Fermi/GBM sGRBs shows that the sGRB associated with GW 170817A is extreme in its combination of flux, spectral softness and temporal structure. We identify a group of similar GRBs, one of which has been associated to a bright galaxy at 75 Mpc. We incorporate off-axis emission in the estimate of the rates of sGRBs, and predict that the majority of future GW-detections of NS-NS mergers will be accompanied by faint gamma-ray emission, contrary to previous thinking. The much more frequent off-axis emission of sGRBs also implies a much higher deadly rate of gamma-rays for extraterrestrial life in the Universe.

PSR J2032+4127 is a young and rapidly rotating pulsar on a highly eccentric orbit around the high-mass Be star MT91 213. X-ray monitoring of the binary system over a $\sim$4000 day period with Swift has revealed an increase of the X-ray luminosity which we attribute to the synchrotron emission of the shocked pulsar wind. We use Swift X-ray observations to infer a clumpy stellar wind with $r^{-2}$ density profile and constrain the Lorentz factor of the pulsar wind to $10^5<\gamma_w<10^6$. We investigate the effects of an axisymmetric stellar wind with polar gradient on the X-ray emission. Comparison of the X-ray light curve hundreds of days before and after the periastron can be used to explore the polar structure of the wind.

We examine millisecond magnetars as central engines of Gamma Ray Bursts' (GRB) prompt emission. Using the proto-magnetar wind model of Metzger et al. 2011, we estimate the temporal evolution of the magnetization and power injection at the base of the GRB jet and apply these to different prompt emission models to make predictions for the GRB energetics, spectra and lightcurves. We investigate both shock and magnetic reconnection models for the particle acceleration, as well as the effects of energy dissipation across optically thick and thin regions of the jet. The magnetization at the base of the jet, $\sigma_0$, is the main parameter driving the GRB evolution in the magnetar model and the emission is typically released for $100\lesssim \sigma_0 \lesssim 3000$. Given the rapid increase in $\sigma_0$ as the proto-magnetar cools and its neutrino-driven mass loss subsides, the GRB duration is typically limited to $\lesssim 100$ s. This low baryon loading at late times challenges magnetar models for ultra-long GRBs, though black hole models likely run into similar difficulties without substantial entrainment from the jet walls. The maximum radiated gamma-ray energy is $\lesssim 5 \times 10^{51}$erg, significantly less than the magnetar's total initial rotational energy and in strong tension with the high end of the observed GRB energy distribution. However, the gradual magnetic dissipation model (Beniamini & Giannios 2017) applied to a magnetar central engine, naturally explains several key observables of typical GRBs, including energetics, durations, stable peak energies, spectral slopes and a hard to soft evolution during the burst.

Plasmoids, quasi-spherical regions of plasma containing magnetic fields and high-energy particles, are a self-consistent by-product of the reconnection process in the relativistic regime. Recent two-dimensional particle-in-cell (PIC) simulations have shown that plasmoids can undergo a variety of processes (e.g. mergers, bulk acceleration, growth, and advection) within the reconnection layer. We developed a Monte Carlo (MC) code, benchmarked with the recent PIC simulations, to examine the effects of these processes on the steady-state size and momentum distributions of the plasmoid chain. The differential plasmoid size distribution is shown to be a power law, $N(w)\propto w^{-\chi}$, ranging from a few plasma skin depths to $\sim 0.1$ of the reconnection layer's length. We demonstrate numerically and analytically that the power law slope $\chi$ is linearly dependent upon the ratio of the plasmoid acceleration and growth rates and that it slightly decreases (from $\sim 2$ to $\sim 1.3$) with increasing plasma magnetization (from 3 to 50). We perform a detailed comparison of our results with those of recent PIC simulations and briefly discuss the astrophysical implications of our findings through the representative case of flaring events from blazar jets.

Jets in long-duration $\gamma$-ray bursts (GRBs) have to drill through the collapsing star in order to break out of it and produce the $\gamma$-ray signal while the central engine is still active. If the breakout time is shorter for more powerful engines, then the jet-collapsar interaction acts as a filter of less luminous jets. We show that the observed broken power-law GRB luminosity function is a natural outcome of this process. For a theoretically motivated breakout time that scales with jet luminosity as $L^{-\chi}$ with $\chi\sim 1/3-1/2$, we show that the shape of the $\gamma$-ray duration distribution can be uniquely determined by the GRB luminosity function and matches the observed one. This analysis has also interesting implications about the supernova-central engine connection. We show that not only successful jets can deposit sufficient energy in the stellar envelope to power the GRB-associated supernovae, but also failed jets may operate in all Type Ib/c supernovae.

Relativistic jets in active galactic nuclei (AGN) convert as much as half of their energy into radiation. To explore the poorly understood processes that are responsible for this conversion, we carry out fully 3D magnetohydrodynamic (MHD) simulations of relativistic magnetized jets. Unlike the standard approach of injecting the jets at large radii, our simulated jets self-consistently form at the source and propagate and accelerate outward for several orders of magnitude in distance before they interact with the ambient medium. We find that this interaction can trigger strong energy dissipation of two kinds inside the jets, depending on the properties of the ambient medium. Those jets that form in a new outburst and drill a fresh hole through the ambient medium fall victim to a 3D magnetic kink instability and dissipate their energy primarily through magnetic reconnection in the current sheets formed by the instability. On the other hand, those jets that form during repeated cycles of AGN activity and escape through a pre-existing hole in the ambient medium maintain their stability and dissipate their energy primarily at MHD recollimation shocks. In both cases the dissipation region can be associated with a change in the density profile of the ambient gas. The Bondi radius in AGN jets serves as such a location.

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

We considered a model for the prompt phase of Gamma-Ray Burst (GRB) emission arising from a magnetized jet undergoing gradual energy dissipation due to magnetic reconnection. The dissipated magnetic energy is translated to bulk kinetic energy and to acceleration of particles. The energy in these particles is released via synchrotron radiation as they gyrate around the strong magnetic fields in the jet. At small radii, the optical depth is large, and the radiation is reprocessed through Comptonization into a narrow, strongly peaked, component. At larger distances the optical depth becomes small and radiation escapes the jet with a non-thermal distribution. The obtained spectra typically peak around $\approx 300$keV (as observed) and with spectral indices below and above the peak that are, for a broad range of the model parameters, close to the observed values. The small radius of dissipation causes the emission to become self absorbed at a few keV and can sufficiently suppress the optical and X-ray fluxes within the limits required by observations (Beniamini & Piran 2014).

A small fraction of Tidal Disruption Events (TDE) produce relativistic jets, evidenced by their non-thermal X-ray spectra and transient radio emission. Here we present milliarcsecond-resolution imaging results on TDE J1644+5734 with the European VLBI Network (EVN). These provide a strong astrometric constraint on the average apparent jet velocity <0.27, that constrains the intrinsic jet velocity for a given viewing angle.

Radiatively inefficient accretion flow models have been shown to accurately account for the spectrum and luminosity observed from Sgr A* in the X-ray regime down to mm wavelengths. However, observations at a few GHz cannot be explained by thermal electrons alone but require the presence of an additional non-thermal particle population. Here, we propose a model for the origin of such a population in the accretion flow via means of a pulsar orbiting the supermassive black hole in our Galaxy. Interactions between the relativistic pulsar wind with the disc lead to the formation of a bow shock in the wind. During the pulsar's transit through the accretion disc, relativistic pairs, accelerated at the shock front, are injected into the disc. The radio-emitting particles are long-lived and remain within the disc long after the pulsar's transit. Periodic pulsar transits through the disc result in regular injection episodes of non-thermal particles. We show that for a pulsar with spin-down luminosity $L_{\rm sd} \sim 3\times 10^{35}$ erg s$^{-1}$ and a wind Lorentz factor of $\gamma_{\rm w} \sim 10^4$ a quasi-steady synchrotron emission is established with luminosities in the $1-10$ GHz range comparable to the observed one.

Dozens of stellar tidal disruption events (TDEs) have been identified at optical, UV and X-ray wavelengths. A small fraction of these, most notably Swift J1644+57, produce radio synchrotron emission, consistent with a powerful, relativistic jet shocking the surrounding circumnuclear gas. The dearth of similar non-thermal radio emission in the majority of TDEs may imply that powerful jet formation is intrinsically rare, or that the conditions in galactic nuclei are typically unfavorable for producing a detectable signal. Here we explore the latter possibility by constraining the radial profile of the gas density encountered by a TDE jet using a one-dimensional model for the circumnuclear medium which includes mass and energy input from a stellar population. Near the jet Sedov radius of 10$^{18}$ cm, we find gas densities in the range of $n_{18} \sim$ 0.1$-$1000 cm$^{-3}$ across a wide range of plausible star formation histories. Using one- and two-dimensional relativistic hydrodynamical simulations, we calculate the synchrotron radio light curves of TDE jets (as viewed both on and off-axis) across the allowed range of density profiles. We find that bright radio emission would be produced across the plausible range of nuclear gas densities by jets as powerful as Swift J1644+57, and we quantify the relationship between the radio luminosity and jet energy. We use existing radio detections and upper limits to constrain the energy distribution of TDE jets. Radio follow up observations several months to several years after the TDE candidate will strongly constrain the energetics of any relativistic flow.

Ap Librae is one out of a handful of low-frequency peaked blazars to be detected at TeV $\gamma$-rays and the only one with an identified X-ray jet. Combined observations of Fermi-LAT at high energies (HE) and of H.E.S.S at very high energies (VHE) revealed a striking spectral property of Ap Librae; the presence of a broad high-energy component that extends more than nine orders of magnitude in energy and is, therefore, hard to be explained by the usual single-zone synchrotron self-Compton model. We show that the superposition of different emission components related to photohadronic interactions can explain the $\gamma$-ray emission of Ap Librae without invoking external radiation fields. We present two indicative model fits to the spectral energy distribution of Ap Librae where the VHE emission is assumed to originate from a compact, sub-pc scale region of the jet. A robust prediction of our model is VHE flux variability on timescales similar to those observed at X-rays and HE $\gamma$-rays, which can be further used to distinguish between a sub-pc or kpc scale origin of the TeV emission. We thus calculate the expected variability signatures at X-rays, HE and VHE $\gamma$-rays and show that quasi-simultaneous flares are expected, with larger amplitude flares appearing at $\gamma$-rays. We assess the detectability of VHE variability from Ap Librae with CTA, next generation of IACTs. We show that $\sim$hr timescale variability at $E_{\gamma}>0.1$ TeV could be detectable at high significance with shorter exposure times than current Cherenkov telescopes.

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.

A large number of supernova remnants (SNRs) in our Galaxy and galaxies nearby have been resolved in various radio bands. This radio emission is thought to be produced via synchrotron emission from electrons accelerated by the shock that the supernova ejecta drives into the external medium. Here we consider the sample of radio SNRs in the Magellanic Clouds. Given the size and radio flux of a SNR, we seek to constrain the fraction of shocked fluid energy in non-thermal electrons ($\epsilon_e$) and magnetic field ($\epsilon_B$), and find $\epsilon_e \epsilon_B \sim 10^{-3}$. These estimates do not depend on the largely uncertain values of the external density and the age of the SNR. We develop a Monte Carlo scheme that reproduces the observed distribution of radio fluxes and sizes of the population of radio SNRs in the Magellanic Clouds. This simple model provides a framework that could potentially be applied to other galaxies with complete radio SNRs samples.

Powerful flares from blazars with short ($\sim$ min) variability timescales are challenging for current models of blazar emission. Here, we present a physically motivated ab initio model for blazar flares based on the results of recent particle-in-cell (PIC) simulations of relativistic magnetic reconnection. PIC simulations demonstrate that quasi-spherical plasmoids filled with high-energy particles and magnetic fields are a self-consistent by-product of the reconnection process. By coupling our PIC-based results (i.e., plasmoid growth, acceleration profile, particle and magnetic content) with a kinetic equation for the evolution of the electron distribution function we demonstrate that relativistic reconnection in blazar jets can produce powerful flares whose temporal and spectral properties are consistent with the observations. In particular, our model predicts correlated synchrotron and synchrotron self-Compton flares of duration of several hours--days powered by the largest and slowest moving plasmoids that form in the reconnection layer. Smaller and faster plasmoids produce flares of sub-hour duration with higher peak luminosities than those powered by the largest plasmoids. Yet, the observed fluence in both types of flares is similar. Multiple flares with a range of flux-doubling timescales (minutes to several hours) observed over a longer period of flaring activity (days or longer) may be used as a probe of the reconnection layer's orientation and the jet's magnetization. Our model shows that blazar flares are naturally expected as a result of magnetic reconnection in a magnetically-dominated jet.

The first-known tidal disruption event (TDE) with strong evidence for a relativistic jet -- based on extensive multi-wavelength campaigns -- is Swift J1644+5734. In order to directly measure the apparent speed of the radio jet, we performed very long baseline interferometry (VLBI) observations with the European VLBI network (EVN) at 5 GHz. Our observing strategy was to identify a very nearby and compact radio source with the real-time e-EVN, and then utilise this source as a stationary astrometry reference point in the later five deep EVN observations. With respect to the in-beam source FIRST J1644+5736, we have achieved a statistical astrometric precision about 12 micro-arcsecond (68 % confidence level) per epoch. This is one of the best phase-referencing measurements available to date. No proper motion has been detected in the Swift J1644+5734 radio ejecta. We conclude that the apparent average ejection speed between 2012.2 and 2015.2 was less than 0.3c with a confidence level of 99 %. This tight limit is direct observational evidence for either a very small viewing angle or a strong jet deceleration due to interactions with a dense circum-nuclear medium, in agreement with some recent theoretical studies.

Blobs, or quasi-spherical emission regions containing relativistic particles and magnetic fields, are often assumed ad hoc in emission models of relativistic astrophysical jets, yet their physical origin is still not well understood. Here, we employ a suite of large-scale two-dimensional particle-in-cell simulations in electron-positron plasmas to demonstrate that relativistic magnetic reconnection can naturally account for the formation of quasi-spherical plasmoids filled with high-energy particles and magnetic fields. Our simulations extend to unprecedentedly long temporal and spatial scales, so we can capture the asymptotic physics independently of the initial setup. We characterize the properties of the plasmoids that are continuously generated as a self-consistent by-product of the reconnection process: they are in rough energy equipartition between particles and magnetic fields; the upper energy cutoff of the plasmoid particle spectrum is proportional to the plasmoid width w, corresponding to a Larmor radius ~0.2 w; the plasmoids grow in size at ~0.1 of the speed of light, with most of the growth happening while they are still non-relativistic (first they grow); their growth is suppressed once they get accelerated to relativistic speeds by the field line tension, up to the Alfven speed (then they go). The largest plasmoids, whose typical recurrence interval is ~2.5 L/c, reach a characteristic size w ~ 0.2 L independently of the system length L, they have nearly isotropic particle distributions and they contain the highest energy particles, whose Larmor radius is ~0.03 L. The latter can be regarded as the Hillas criterion for relativistic reconnection. We briefly discuss the implications of our results for the high-energy emission from relativistic jets and pulsar winds.

Sgr A* is an ideal target to study low-luminosity accreting systems. It has been recently proposed that properties of the accretion flow around Sgr A* can be probed through its interactions with the stellar wind of nearby massive stars belonging to the S-cluster. When a star intercepts the accretion disk, the ram and thermal pressures of the disk terminate the stellar wind leading to the formation of a bow shock structure. Here, a semi-analytical model is constructed which describes the geometry of the termination shock formed in the wind. With the employment of numerical hydrodynamic simulations, this model is both verified and extended to a region prone to Kelvin-Helmholtz instabilities. Because the characteristic wind and stellar velocities are in $\sim10^{8}$ cm s$^{-1}$ range, the shocked wind may produce detectable X-rays via thermal bremsstrahlung emission. The application of this model to the pericenter passage of S2, the brightest member of the S-cluster, shows that the shocked wind produces roughly a month long X-ray flare with a peak luminosity of $L\approx 4 \times 10^{33}$ erg s$^{-1}$ for a stellar mass-loss rate, disk number density, and thermal pressure strength of $\dot{M}_{\rm w}= 10^{-7} M_\odot\, {\rm yr}^{-1}$, $n_{\rm d} = 10^{5}$ cm$^{-3}$, and $\alpha=0.1$, respectively. This peak luminosity is comparable to the quiescent X-ray emission detected from Sgr A* and is within the detection capabilities of current X-ray observatories. Its detection could constrain the density and thickness of the disk at a distance of $\sim 3000$ gravitational radii from the supermassive black hole.

In this work we focus on the technical details of the numerical simulations of the non- thermal transient Swift J1644+57, whose emission is probably produced by a two-component jet powered by a tidal disruption event. In this context we provide details of the coupling between the relativistic hydrodynamic simulations and the radiative transfer code. First, we consider the technical demands of one-dimensional simulations of a fast relativistic jet, and show to what extent (for the same physical parameters of the model) do the computed light curves depend on the numerical parameters of the different codes employed. In the second part we explain the difficulties of computing light curves from axisymmetric two dimensonal simulations and discuss a procedure that yields an acceptable tradeoff between the computational cost and the quality of the results.

Despite intensive searches, the only pulsar within 0.1 pc of the central black hole in our Galaxy, Sgr A*, is a radio-loud magnetar. Since magnetars are rare among the Galactic neutron star population, and a large number of massive stars are already known in this region, the Galactic centre (GC) should harbor a large number of neutron stars. Population syntheses suggest several thousand neutron stars may be present in the GC. Many of these could be highly energetic millisecond pulsars which are also proposed to be responsible for the GC gamma-ray excess. We propose that the presence of a neutron star within 0.03~pc from Sgr~A* can be revealed by the shock interactions with the disk around the central black hole. As we demonstrate, these interactions result in observable transient non-thermal X-ray and gamma-ray emission over timescales of months, provided that the spin down luminosity of the neutron star is L_{sd}~10^{35} erg/s. Current limits on the population of normal and millisecond pulsars in the GC region suggest that a number of such pulsars are present with such luminosities.

In order to explain rapid light curve variability without invoking a variable source, several authors have proposed "minijets" that move relativistically relative to the main flow of the jet. Here we consider the possibility that these minijets, instead of being isotropically distributed in the comoving frame of the jet, form primarily perpendicular to the direction of the flow, as the jet dissipates its energy at a large emission radius. This yields two robust features. First, the emission is significantly delayed compared with the isotropic case. This delay allows for the peak of the afterglow emission to appear while the source is still active, in contrast to the simplest isotropic model. Secondly, the flux decline after the source turns off is steeper than the isotropic case. We find that these two features are realized in gamma-ray bursts (GRBs): 1. The peak of most GeV light curves (ascribed to the external shock) appears during the prompt emission phase. 2. Many X-ray light curves exhibit a period of steep decay, which is faster than that predicted by the standard isotropic case. The gamma-ray generation mechanism in GRBs, and possibly in other relativistic flows, may therefore be anisotropic.

The gamma-ray burst (GRB) jet powers the afterglow emission by shocking the surrounding medium, and radio afterglow can now be routinely observed to almost a year after the explosion. Long-duration GRBs are accompanied by supernovae (SNe) that typically contain much more energy than the GRB jet. Here we consider the fact that the SN blast wave will also produce its own afterglow (supernova remnant emission), which will peak at much later time (since it is non-relativistic), when the SN blast wave transitions from a coasting phase to a decelerating Sedov-Taylor phase. We predict that this component will peak generally a few tens of years after the explosion and it will outshine the GRB powered afterglow well-before its peak emission. In the case of GRB 030329, where the external density is constrained by the $\sim 10$-year coverage of the radio GRB afterglow, the radio emission is predicted to start rising over the next decade and to continue to increase for the following decades up to a level of $\sim$ mJy. Detection of the SN-powered radio emission will greatly advance our knowledge of particle acceleration in $ \sim 0.1$c shocks.

We model the non-thermal transient Swift J1644+57 as resulting from a relativistic jet powered by the accretion of a tidally-disrupted star onto a super-massive black hole. Accompanying synchrotron radio emission is produced by the shock interaction between the jet and the dense circumnuclear medium, similar to a gamma-ray burst afterglow. An open mystery, however, is the origin of the late-time radio rebrightening, which occurred well after the peak of the jetted X-ray emission. Here, we systematically explore several proposed explanations for this behavior by means of multi-dimensional hydrodynamic simulations coupled to a self-consistent radiative transfer calculation of the synchrotron emission. Our main conclusion is that the radio afterglow of Swift J1644+57 is not naturally explained by a jet with a one-dimensional top-hat angular structure. However, a more complex angular structure comprised of an ultra-relativistic core (Lorentz factor $\Gamma \sim 10$) surrounded by a slower ($\Gamma \sim $ 2) sheath provides a reasonable fit to the data. Such a geometry could result from the radial structure of the super-Eddington accretion flow or as the result of jet precession. The total kinetic energy of the ejecta that we infer of $\sim$ few $10^{53}\,$erg requires a highly efficient jet launching mechanism. Our jet model providing the best fit to the light curve of the on-axis event Swift J1644+57 is used to predict the radio light curves for off-axis viewing angles. Implications for the presence of relativistic jets from TDEs detected via their thermal disk emission, as well as the prospects for detecting orphan TDE afterglows with upcoming wide-field radio surveys and resolving the jet structure with long baseline interferometry, are discussed.

Observations of gamma-ray-bursts and jets from active galactic nuclei reveal that the jet flow is characterized by a high radiative efficiency and that the dissipative mechanism must be a powerful accelerator of non-thermal particles. Shocks and magnetic reconnection have long been considered as possible candidates for powering the jet emission. Recent progress via fully-kinetic particle-in-cell simulations allows us to revisit this issue on firm physical grounds. We show that shock models are unlikely to account for the jet emission. In fact, when shocks are efficient at dissipating energy, they typically do not accelerate particles far beyond the thermal energy, and vice versa. In contrast, we show that magnetic reconnection can deposit more than 50% of the dissipated energy into non-thermal leptons as long as the energy density of the magnetic field in the bulk flow is larger than the rest mass energy density. The emitting region, i.e., the reconnection downstream, is characterized by a rough energy equipartition between magnetic fields and radiating particles, which naturally accounts for a commonly observed property of blazar jets.

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV gamma-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, both in shape and peak flux density. Assuming that the radio and gamma-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major gamma-ray and radio flares of 2012 using the same framework. The 2012 gamma-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding gamma-ray flares, can be reproduced only for a very specific choice of parameters.

This review discusses the physics of magnetic reconnection, a process in which the magnetic field topology changes and magnetic energy is converted to kinetic energy, in pair plasmas in the relativistic regime. We focus on recent progress in the field driven by theory advances and the maturity of particle-in-cell codes. This work shows that fragmentation instabilities at the current sheet can play a critical role in setting the reconnection speed and affect the resulting particle acceleration, anisotropy, bulk flows, and radiation. Then, we discuss how this novel understanding of relativistic reconnection can be applied to high-energy astrophysical phenomena, with an emphasis on pulsars, pulsar wind nebulae, and active galactic nucleus jets.

We propose a scenario for launching relativistic jets from rotating black holes, in which small-scale magnetic flux loops, sustained by disc turbulence, are forced to inflate and open by differential rotation between the black hole and the accretion flow. This mechanism does not require a large-scale net magnetic flux in the accreting plasma. Estimates suggest that the process could operate effectively in many systems, and particularly naturally and efficiently when the accretion flow is retrograde. We present the results of general-relativistic force-free electrodynamic simulations demonstrating the time evolution of the black hole's magnetosphere, the cyclic formation of jets, and the effect of magnetic reconnection. The jets are highly variable on timescales ~ 10-10^3 r_ g/c, where r_g is the black hole's gravitational radius. The reconnecting current sheets observed in the simulations may be responsible for the hard X-ray emission from accreting black holes.

Long-duration gamma-ray bursts (GRBs) are thought to come from the core-collapse of Wolf-Rayet stars. Whereas their stellar masses $M_*$ have a rather narrow distribution, the population of GRBs is very diverse, with gamma-ray luminosities $L_\gamma$ spanning several orders of magnitude. This suggests the existence of a "hidden" stellar variable whose burst-to-burst variation leads to a spread in $L_\gamma$. Whatever this hidden variable is, its variation should not noticeably affect the shape of GRB lightcurves, which display a constant luminosity (in a time-average sense) followed by a sharp drop at the end of the burst seen with Swift/XRT. We argue that such a hidden variable is progenitor star's large-scale magnetic flux. Shortly after the core collapse, most of stellar magnetic flux accumulates near the black hole (BH) and remains there. The flux extracts BH rotational energy and powers jets of roughly a constant luminosity, $L_j$. However, once BH mass accretion rate $\dot M$ falls below $\sim L_j/c^2$, the flux becomes dynamically important and diffuses outwards, with the jet luminosity set by the rapidly declining mass accretion rate, $L_j\sim \dot M c^2$. This provides a potential explanation for the sharp end of GRBs and the universal shape of their lightcurves. During the GRB, gas infall translates spatial variation of stellar magnetic flux into temporal variation of $L_j$. We make use of the deviations from constancy in $L_j$ to perform stellar magnetic flux "tomography". Using this method, we infer the presence of magnetised tori in the outer layers of progenitor stars for GRB 920513 and GRB 940210.

We report the polarized optical light curve of a gamma-ray burst afterglow obtained using the RoboPol instrument. Observations began 655 seconds after the initial burst of gamma-rays from GRB131030A, and continued uninterrupted for 2 hours. The afterglow displayed a low, constant fractional linear polarization of $p = (2.1 \pm 1.6)\,\%$ throughout, which is similar to the interstellar polarization measured on nearby stars. The optical brightness decay is consistent with a forward-shock propagating in a medium of constant density, and the low polarization fraction indicates a disordered magnetic field in the shock front. This supports the idea that the magnetic field is amplified by plasma instabilities on the shock front. These plasma instabilities produce strong magnetic fields with random directions on scales much smaller than the total observable region of the shock, and the resulting randomly-oriented polarization vectors sum to produce a low net polarization over the total observable region of the shock.

The recent discovery of extragalactic PeV neutrinos opens a new window to the exploration of cosmic-ray accelerators. The observed PeV neutrino flux is close to the Waxman-Bahcall upper bound implying that gamma-ray bursts (GRBs) may be the source of ultra-high energy cosmic rays (UHECRs). Starting with the assumption of the GRB-UHECR connection, we show using both analytical estimates and numerical simulations that the observed neutrinos can originate at the jet as a result of photopion interactions with the following implications: the neutrino spectra are predicted to have a cutoff at energy $\le 10$ PeV; the dissipation responsible for the GRB emission and cosmic-ray acceleration takes place at distances $r_{\rm diss} \approx 3 \times 10^{11}-3 \times 10^{13} \ {\rm cm}$ from the central engine; the Thomson optical depth at the dissipation region is $\tau_{\rm T} \sim 1$; the jet carries a substantial fraction of its energy in the form of Poynting flux at the dissipation region, and has a Lorentz factor $\Gamma \simeq 100-500$. The non-detection of $\sim$PeV neutrinos coincident with GRBs will indicate that GRBs are either poor cosmic accelerators or the dissipation takes place at small optical depths in the jet.

The mechanism behind the launching of gamma-ray-burst (GRB) jets remains debated resulting in large uncertainty over the jet composition. Both magnetohydrodynamical and neutrino annihilation models have been proposed for the energy extraction in a black hole/accretion-disc central engine. In particular, for the extreme accretion rates $\dot M\sim 0.1-1$~ M$_\odot$s$^{-1}$ expected for bursts of duration $T\lesssim 100$~s, the disc can be an efficient neutrino emitter. Neutrino-antineutrino annihilation results in an energy deposition rate at the jet that can, in principle, account for the burst's energetics. Recent discoveries of X-ray flares hours after the burst and of ultra-long GRBs suggest that GRB activity can last for $\sim 10^4$~s or longer. These long-lived events have fluence similar to that of classical GRBs. In view of these findings, we re-evaluate the neutrino annihilation model. We derive the maximum possible energy of a neutrino-powered jet as a function of the burst duration and show that the available energy drops fast for longer bursts. For a standard choice of the parameters, the model falls short by three to four orders of magnitude in explaining the observed energetics of events that last longer than $\sim 10^3$~s.

We explore a one-zone hadronic model that may be able to reproduce $\gamma$-ray burst (GRB) prompt emission with a minimum of free parameters. Assuming only that GRBs are efficient high-energy proton accelerators and without the presence of an ab initio photon field, we investigate the conditions under which the system becomes supercritical, i.e. there is a fast, non-linear transfer of energy from protons to secondary particles initiated by the spontaneous quenching of proton-produced $\gamma$-rays. We first show analytically that the transition to supercriticality occurs whenever the proton injection compactness exceeds a critical value, which favours high proton injection luminosities and a wide range of bulk Lorentz factors. The properties of supercriticality are then studied with a time-dependent numerical code that solves concurrently the coupled equations of proton, photon, electron, neutron and neutrino distributions. For conditions that drive the system deep into the supercriticality we find that the photon spectra obtain a Band-like shape due to Comptonization by cooled pairs and that the energy transfer efficiency from protons to $\gamma$-rays and neutrinos is high reaching $\sim 0.3$. Although some questions concerning its full adaptability to the GRB prompt emission remain open, supercriticality is found to be a promising process in that regard.

The interaction of TeV photons from blazars with the extragalactic background light produces a relativistic beam of electron-positron pairs streaming through the intergalactic medium (IGM). The fate of the beam energy is uncertain. By means of two- and three-dimensional particle-in-cell simulations, we study the non-linear evolution of dilute ultra-relativistic pair beams propagating through the IGM. We explore a wide range of beam Lorentz factors gamma_b>>1 and beam-to-plasma density ratios alpha<<1, so that our results can be extrapolated to the extreme parameters of blazar-induced beams (gamma_b~10^6 and alpha~10^(-15), for the most powerful blazars). For cold beams, we show that the oblique instability governs the early stages of evolution, but its exponential growth terminates - due to self-heating of the beam in the transverse direction - when only a negligible fraction ~(alpha/gamma_b)^(1/3)~10^(-7) of the beam energy has been transferred to the IGM plasma. Further relaxation of the beam proceeds through quasi-longitudinal modes, until the momentum dispersion in the direction of propagation saturates at DeltaP_{b,parallel}/gamma_b m_e c ~ 0.2. This corresponds to a fraction ~10% of the beam energy being ultimately transferred to the IGM plasma, irrespective of gamma_b or alpha. If the initial dispersion in beam momentum satisfies DeltaP_{b,parallel}/gamma_b m_e c > 0.2 (as typically expected for blazar-induced beams), the fraction of beam energy deposited into the IGM is much smaller than ~10%. It follows that at least ~90% of the beam energy is still available to power the GeV emission produced by inverse Compton up-scattering of the Cosmic Microwave Background by the beam pairs.

The afterglow emission from Gamma-Ray Bursts (GRBs) is usually interpreted as synchrotron radiation from relativistic electrons accelerated at the GRB external shock, that decelerates from ultra-relativistic to non-relativistic speeds as it sweeps up the surrounding medium. We investigate the temporal decay of the emission from GRB afterglows at late times, when the bulk of the shock-accelerated electrons are non-relativistic. For a uniform circumburst medium, we show that such "deep Newtonian phase" begins at t_{DN} ~ 3 epsilon_{e,-1}^{5/6} t_{ST}, where t_{ST} marks the transition of the blast wave to the non-relativistic spherically-symmetric Sedov-Taylor solution, and epsilon_e = 0.1 epsilon_{e,-1} quantifies the amount of shock energy transferred to the post-shock electrons. For typical parameters, the deep Newtonian stage starts ~ 0.5-several years after the GRB. The radio flux in this phase decays as F_nu ~ t^{-3(p+1)/10}, with decay slope between -0.9 and -1.2, for a power-law distribution of shock-accelerated electrons with index 2 < p < 3. This is shallower than the commonly assumed scaling F_nu ~ t^{-3(5p-7)/10} (with decay slope between -0.9 and -2.4) derived by Frail et al. (2000), which only applies if the GRB shock is non-relativistic, but the electron distribution still peaks at ultra-relativistic energies (a regime that is relevant for a narrow time interval, and only if t_{DN} > t_{ST}, or epsilon_e > 0.03). We discuss how the deep Newtonian phase can be reliably used for GRB calorimetry, and we comment on the good detection prospects of trans-relativistic blast waves at 0.1-10 GHz with EVLA and LOFAR.

Minute-timescale TeV flares have been observed in several blazars. The fast flaring requires compact regions in the jet that boost their emission towards the observer at an extreme Doppler factor of delta>50. For TeV photons to avoid annihilation in the broad line region of PKS 1222+216, the flares must come from large (pc) scales challenging most models proposed to explain them. Here we elaborate on the magnetic reconnection minijet model for the blazar flaring, focusing on the inherently time-dependent aspects of the process of magnetic reconnection. I argue that, for the physical conditions prevailing in blazar jets, the reconnection layer fragments leading to the formation a large number of plasmoids. Occasionally a plasmoid grows to become a large, "monster" plasmoid. I show that radiation emitted from the reconnection event can account for the observed "envelope" of ~ day-long blazar activity while radiation from monster plasmoids can power the fastest TeV flares. The model is applied to several blazars with observed fast flaring. The inferred distance of the dissipation zone from the black hole and the typical size of the reconnection regions are R_diss~0.3-1 pc and l' \simless 10^{16} cm, respectively. The required magnetization of the jet at this distance is modest: sigma ~a few. Such dissipation distance R_diss and reconnection size l' are expected if the jet contains field structures with size of the order of the black-hole horizon.

How accretion proceeds around the massive black hole in the Galactic center and other highly sub-Eddington accretors remains poorly understood. The orbit of the S2 star in the Galactic center passes through the accretion disk of the massive black hole and any observational signature from such interaction may be used as an accretion probe. Because of its early stellar type, S2 is expected to possess a fairly powerful wind. We show here that the ram pressure of the accretion disk shocks the stellar wind fairly close to the star. The shocked fluid reaches a temperature of ~ 1 keV and cools efficiently through optically thin, thermal bremsstrahlung emission. The radiation from the shocked wind peaks around the epoch of the pericenter passage of the star at a luminosity potentially comparable to the quiescent emission detected from Sgr A*. Detection of shocked wind radiation can constrain the density of the accretion disk at a distance of several thousands of gravitational radii from the black hole.

The unusual transient Swift J1644+57 likely resulted from a collimated relativistic jet powered by accretion onto a massive black hole (BH) following the tidal disruption (TD) of a star. Several mysteries cloud the interpretation of this event: (1) extreme flaring and `plateau' shape of the X-ray/gamma-ray light curve during the first 10 days after the gamma-ray trigger; (2) unexpected rebrightening of the forward shock radio emission months after trigger; (3) no obvious evidence for jet precession, despite misalignment typically expected between the angular momentum of the accretion disk and BH; (4) recent abrupt shut-off in jet X-ray emission after 1.5 years. Here we show that all of these seemingly disparate mysteries are naturally resolved by one assumption: the presence of strong magnetic flux Phi threading the BH. Initially, Phi is weak relative to high fall-back mass accretion rate, Mdot, and the disk and jets precess about the BH axis = our line of sight. As Mdot drops, Phi becomes dynamically important and leads to a magnetically-arrested disk (MAD). MAD naturally aligns disk and jet axis along the BH spin axis, but only after a violent rearrangement phase (jet wobbling). This explains the erratic light curve at early times and the lack of precession at later times. We use our model for Swift J1644+57 to constrain BH and disrupted star properties, finding that a solar-mass main sequence star disrupted by a relatively low mass, M~10^5-10^6 Msun, BH is consistent with the data, while a WD disruption (though still possible) is disfavored. The magnetic flux required to power Swift J1644+57 is too large to be supplied by the star itself, but it could be collected from a quiescent `fossil' accretion disk present in the galactic nucleus prior to the TD. The presence (lack of) of such a fossil disk could be a deciding factor in what TD events are accompanied by powerful jets.[abridged]

The recent transient event Swift J1644+57 has been interpreted as emission from a collimated relativistic jet, powered by the sudden onset of accretion onto a supermassive black hole following the tidal disruption of a star. Here we model the radio-microwave emission as synchrotron radiation produced by the shock interaction between the jet and the gaseous circumnuclear medium (CNM). At early times after the onset of the jet (t < 5-10 days) a reverse shock propagates through and decelerates the ejecta, while at later times the outflow approaches the Blandford-McKee self-similar evolution (possibly modified by additional late energy injection). The achromatic break in the radio light curve of J1644+57 is naturally explained as the transition between these phases. We show that the temporal indices of the pre- and post-break light curve are consistent with those predicted if the CNM has a wind-type radial density profile n ~ 1/r^2. The observed synchrotron frequencies and self-absorbed flux constrain the fraction of the post-shock thermal energy in relativistic electrons epsilon_e ~ 0.03-0.1; the CNM density at 1e18 cm ~ 1-10 1/cm^3; and the initial Lorentz factor Gamma_j ~ 10-20 and opening angle theta_j ~ (0.1-1)Gamma_j^(-1) ~ 0.01-0.1 of the jet. Radio modeling thus provides robust independent evidence for a narrowly collimated outflow. Extending our model to the future evolution of J1644+57, we predict that the radio flux at low frequencies (<~ few GHz) will begin to brighten more rapidly once the characteristic frequency crosses below the self-absorption frequency on a timescale of months (indeed, such a transition may already have begun). Our results demonstrate that relativistic outflows from tidal disruption events provide a unique probe of the conditions in distant, previously inactive galactic nuclei, complementing studies of normal AGN.

We evaluate the prompt observational signatures of the merger between a massive close-in planet (a `hot Jupiter') and its host star, events with an estimated Galactic rate of ~0.1-1/yr. Depending on the ratio of the mean density of the planet rho_p to that of the star rho_star, a merger results in three possible outcomes. If rho_p/rho_star > 5, then the planet directly plunges below the stellar atmosphere before being disrupted by tidal forces. Dissipation of orbital energy creates a hot wake behind the planet, producing a EUV/soft X-ray transient as the planet sinks below the stellar surface. The peak luminosity L_X ~ 1e36 erg/s is achieved weeks-months prior to merger, after which the stellar surface is enshrouded by an outflow. The final inspiral is accompanied by an optical transient powered by the recombination of hydrogen in the outflow, which peaks at L~1e37-38 erg/s on a timescale ~days. If instead rho_planet/rho_star < 5, then Roche Lobe overflow occurs above the stellar surface. For rho_p/rho_star < 1, mass transfer is stable, resulting the planet being accreted on a relatively slow timescale. However, for 1 < rho_p/rho_star < 5, mass transfer may instead be unstable, resulting in the planet being dynamically disrupted into an accretion disk around the star. Super-Eddington outflows from the disk power an optical transient with L~1e37-38 erg/s and characteristic duration ~week-months. The disk itself becomes visible once the accretion rate become sub-Eddington, resulting in a bolometric brightening and spectral shift to the UV. Optical transients from planet merger events may resemble classical novae, but are distinguished by lower ejecta mass and velocity ~100s km/s, and by hard pre- and post-cursor emission, respectively. Promising search strategies include combined optical, UV, and X-ray surveys of nearby massive galaxies with cadences from days to months.

The radiation released at the transparency radius of an ultrarelativistic flow can account for the observed properties of gamma-ray bursts (GRBs) provided that sufficient energy is dissipated in the sub-photospheric region. Here, I investigate how the peak energy of the E*f(E) spectrum and its overall shape depend on the properties of the jet for various "dissipative photospheres". I find that continuous energy release which results in electron heating over a wide range of distances may be the key to explain the GRB emission. In this picture, the peak of the spectrum forms at a Thomson optical depth of several tens. The peak depends mainly on the bulk Lorentz factor Gamma of the flow and can, therefore, be used to determine it. The Gamma is predicted to range from ~10 to 1000 from X-ray flashes to the brightest observed GRBs in agreement with recent observational inferences. The Amati relation can be understood if the brightest bursts are the least baryon loaded ones. Implications from this interpretation of the GRB emission for the central engine are discussed.

The afterglow emission that follows gamma-ray bursts (GRBs) contains valuable information about the circumburst medium and, therefore, about the GRB progenitor. Theoretical studies of GRB blast waves, however, are often limited to simple density profiles for the external medium (mostly constant density and power-law R^{-k} ones). We argue that a large fraction of long-duration GRBs should take place in massive stellar clusters where the circumburst medium is much more complicated. As a case study, we simulate the propagation of a GRB blast wave in a medium shaped by the collision of the winds of O and Wolf-Rayet stars, the typical distance of which is d /sim 0.1 - 1 pc. Assuming a spherical blast wave, the afterglow light curve shows a flattening followed by a shallow break on a timescale from hours up to a week after the burst, which is a result of the propagation of the blast wave through the shocked wind region. If the blast wave is collimated, the jet break may, in some cases, become very pronounced with the post break decline of the light curve as steep as t-5. Inverse Compton scattering of ultra-violet photons from the nearby star off energetic electrons in the blast wave leads to a bright \simGeV afterglow flare that may be detectable by Fermi.

Long duration Gamma-Ray Bursts (GRBs) originate from the core collapse of massive stars, but the identity of the central engine remains elusive. Previous work has shown that rapidly spinning, strongly magnetized proto-neutron stars (`millisecond proto-magnetars') produce outflows with energies, timescales, and magnetizations sigma_0 (maximum Lorentz factor) that are consistent with those required to produce long GRBs. Here we extend this work in order to construct a self-consistent model that directly connects the properties of the central engine to the observed prompt emission. Just after the launch of the supernova shock, a wind heated by neutrinos is driven from the proto-magnetar. The outflow is collimated into a bipolar jet by its interaction with the star. As the magnetar cools, the wind becomes ultra-relativistic and Poynting-flux dominated (sigma_0 >> 1) on a timescale comparable to that required for the jet to clear a cavity through the star. Although the site and mechanism of the prompt emission are debated, we calculate the emission predicted by two models: magnetic dissipation and internal shocks. Our results favor the magnetic dissipation model in part because it predicts a relatively constant `Band' spectral peak energy E_peak with time during the GRB. The jet baryon loading decreases abruptly when the neutron star becomes transparent to neutrinos at t ~ 10-100 seconds. Jets with ultra-high magnetization cannot effectively accelerate and dissipate their energy, suggesting this transition ends the prompt emission and may explain the steep decay phase that follows. We assess several phenomena potentially related to magnetar birth, including low luminosity GRBs, thermal-rich GRBs/X-ray Flashes, very luminous supernovae, and short duration GRBs with extended emission.

Recent measurements by the Pierre Auger Observatory suggest that the composition of ultra-high energy cosmic rays (UHECRs) becomes dominated by heavy nuclei at high energies. However, until now there has been no astrophysical motivation for considering a source highly enriched in heavy elements. Here we demonstrate that the outflows from Gamma-Ray Bursts (GRBs) may indeed be composed primarily of nuclei with masses A ~ 40-200, which are synthesized as hot material expands away from the central engine. In particular, if the jet is magnetically-dominated (rather than a thermally-driven fireball) its low entropy enables heavy elements to form efficiently. Adopting the millisecond proto-magnetar model for the GRB central engine, we show that heavy nuclei are both synthesized in proto-magnetar winds and can in principle be accelerated to energies >1e20 eV in the shocks or regions of magnetic reconnection that are responsible for powering the GRB. Similar results may apply to accretion-powered GRB models if the jet originates from a magnetized disk wind. Depending on the precise distribution of nuclei synthesized, we predict that the average primary mass may continue to increase beyond Fe group elements at the highest energies, possibly reaching the A ~ 90 (Zirconium), A ~ 130 (Tellurium), or even A ~ 195 (Platinum) peaks. Future measurements of the UHECR composition at energies >~ 1e20 eV can thus confirm or constrain our model and, potentially, probe the nature of GRB outflows. The longer attenuation length of ultra-heavy nuclei through the extragalactic background light greatly expands the volume of accesible sources and alleviates the energetic constraints on GRBs as the source of UHECRs.

While gas accretion onto some massive black holes (MBHs) at the centers of galaxies actively powers luminous emission, the vast majority of MBHs are considered dormant. Occasionally, a star passing too near a MBH is torn apart by gravitational forces, leading to a bright panchromatic tidal disruption flare (TDF). While the high-energy transient Swift J164449.3+573451 ("Sw 1644+57") initially displayed none of the theoretically anticipated (nor previously observed) TDF characteristics, we show that the observations (Levan et al. 2011) suggest a sudden accretion event onto a central MBH of mass ~10^6-10^7 solar masses. We find evidence for a mildly relativistic outflow, jet collimation, and a spectrum characterized by synchrotron and inverse Compton processes; this leads to a natural analogy of Sw 1644+57 with a smaller-scale blazar. The phenomenologically novel Sw 1644+57 thus connects the study of TDFs and active galaxies, opening a new vista on disk-jet interactions in BHs and magnetic field generation and transport in accretion systems.

The tidal disruption of a star by a supermassive black hole provides us with a rare glimpse of these otherwise dormant beasts. It has long been predicted that the disruption will be accompanied by a thermal `flare', powered by the accretion of bound stellar debris. Several candidate disruptions have been discovered in this manner at optical, UV and X-ray wavelengths. Here we explore the observational consequences if a modest fraction of the accretion power is channeled into an ultra-relativistic outflow. We show that a relativistic jet decelerates due to its interaction with the interstellar medium at sub-parsec distances from the black hole. Synchrotron radiation from electrons accelerated by the reverse shock powers a bright radio-infrared transient that peaks on a timescale ~1 yr after disruption. Emission from the forward shock may be detectable for several years after the peak. Deep radio follow-up observations of tidal disruption candidates at late times can test for the presence of relativistic ejecta. Upcoming radio transient surveys may independently discover tens to hundreds of tidal disruptions per year, complimenting searches at other wavelengths. Non-thermal emission from tidal disruption probes the physics of jet formation under relatively clean conditions, in which the flow parameters are independently constrained.

Aims: GRB 050502B is well known for the very bright flare displayed in its X-ray light curve. Despite extensive studies, however, the optical light curve has never been discussed and its redshift is unconstrained. Possible correlations between optical and X-ray data are analysed. Methods: Photometric data from TNG in the R and I bands were used to compare the optical afterglow with the X-ray light curve. The HyperZ package and a late time VLT host observation were used to derive redshift estimates. Results: The I-band afterglow decay followed a power-law of index {\alpha} = 2.1 $\pm$ 0.6, after a late break at ~ $1.3 \times 10^5$ s. The R - I color is remarkably red and the broadband spectral index {\beta}_OX = $0.9 \pm 0.1$ is consistent with the X-ray spectral slope {\beta}_X. Although a photometric redshift of z > 4 is the most conservative result to consider, a photometric redshift of z = 5.2 $\pm$ 0.3 is suggested with no extinction in the host, based on which an isotropic energy E_{\gamma},iso = $(3.8 \pm 0.7) \times 10^52$ erg and a jet opening angle {\theta} ~ $3.7^{\circ}$ are subsequently derived. Conclusions: The combined X-ray and optical data suggest an achromatic break, which we interpret as a jet break. The post jet break slope obeys roughly the closure relation for the jet slow cooling model. Because of the afterglow's very red color, in order for the redshift to be low (z < 1), extinction, if present in the host, must be significantly high. Since the optical-to-X-ray index is consistent with the X-ray spectrum, and there is no XRT evidence for excess N_H, GRB 050502B was likely at high redshift.

Long-duration gamma-ray bursts (GRBs) may be powered by the rotational energy of a millisecond magnetar. I argue that the GRB-driving magnetars lie at the high end of the distribution of magnetic field strengths of magnetars. The field of GRB magnetars decays on timescale of hundreds of years and can power SGR-like flares up to ~100 times more powerful than the 2004 event of SGR 1806--20. A few of these flares per year may have been observed by {\it BATSE} and classified as short-duration GRBs. Association of one of these superflares with a nearby d_L\simless 250 Mpc galaxy and the discovery of a, coincident in space, 100-year-old GRB afterglow (observed in the radio) will be the characteristic signature of the magnetar model for GRBs.

Ultra-high energy cosmic rays (UHECRs) may be produced in active galactic nuclei (AGN) or gamma-ray burst (GRB) jets. I argue that magnetic reconnection in jets can accelerate UHECRs rather independently of physical processes in the magnetic dissipation region. First order Fermi acceleration can efficiently take place in the region where the unreconnected (upstream) magnetized fluid converges into the reconnection layer. I find that protons can reach energies up to E~10^{20} eV in GRB and powerful AGN jets while iron nuclei can reach similar energies in AGN jets of more moderate luminosity.

We construct a numerical model of emission from minijets, localized flows driven by magnetic reconnection inside Poynting-flux-dominated jets proposed to explain the ultrafast variability of blazars. The geometrical structure of the model consists of two wedge-like regions of relativistically flowing gas, separated by a stationary shock. The dynamics is based on solutions of relativistic magnetic reconnection with a guide field from Lyubarsky (2005). Electron distributions in each region are chosen to the match the pressure and density of the local plasma. Synchrotron emission from both regions is used to calculate Compton scattering, Compton drag and photon-photon opacity effects, with exact treatment of anisotropy and the Klein-Nishina regime. Radiative effects on plasma are taken into account, including the dependence of pressure on electron radiative losses and adiabatic heating of the flow decelerating under Compton drag. The results are applied to the July 2006 flare in the BL Lac object PKS 2155-304, with the aim of matching TeV flux measurements by H.E.S.S. with models that satisfy the variability constraints, while keeping X-ray emission below simultaneous Chandra observations. We find that models of isolated minijets with a significant guide field overproduce X-ray emission, and that we must take into account the radiative interaction of oppositely-oriented minijets in order to achieve a high enough dominance by Comptonized TeV radiation. We argue that such interactions are likely to occur in a jet where there is substantial internal reconnection, producing a large number of misaligned minijets. Finally, we show that large jet magnetizations are indeed required to satisfy all observational constraints and that the effective Lorentz factor of the minijet plasma has to be larger than 50, in agreement with earlier one-zone estimates.

We use high-resolution relativistic MHD simulations coupled with a radiative transfer code to compute multiwavelength afterglow light curves of magnetized ejecta of gamma-ray bursts interacting with a uniform circumburst medium. The aim of our study is to determine how the magnetization of the ejecta at large distance from the central engine influences the afterglow emission, and to assess whether observations can be reliably used to infer the strength of the magnetic field. We find that, for typical parameters of the ejecta, the emission from the reverse shock peaks for magnetization $\sigma_0 \sim 0.01 - 0.1$ of the flow, and that it is greatly suppressed for higher $\sigma_0$. The emission from the forward shock shows an achromatic break shortly after the end of the burst marking the onset of the self-similar evolution of the blast wave. Fitting the early afterglow of GRB 990123 and 090102 with our numerical models we infer respective magnetizations of $\sigma_0 \sim 0.01$ and $\sigma_0 \sim 0.1$ for these bursts. We argue that the lack of observed reverse shock emission from the majority of the bursts can be understood if $\sigma_0 \simmore 0.1$, since we obtain that the luminosity of the reverse shock decreases significantly for $\sigma_0 \sim 1$. For ejecta with $\sigma_0 \simmore 0.1$ our models predict that there is sufficient energy left in the magnetic field, at least during an interval of ~10 times the burst duration, to produce a substantial emission if the magnetic energy can be dissipated (for instance, due to resistive effects) and radiated away.

We present an updated catalog of 113 X-ray flares detected by Swift in the ~33% of the X-ray afterglows of Gamma-Ray Bursts (GRB). 43 flares have a measured redshift. For the first time the analysis is performed in 4 different X-ray energy bands, allowing us to constrain the evolution of the flare temporal properties with energy. We find that flares are narrower at higher energies: their width follows a power-law relation w~E^{-0.5} reminiscent of the prompt emission. Flares are asymmetric structures, with a decay time which is twice the rise time on average. Both time scales linearly evolve with time, giving rise to a constant rise-to-decay ratio: this implies that both time scales are stretched by the same factor. As a consequence, the flare width linearly evolves with time to larger values: this is a key point that clearly distinguishes the flare from the GRB prompt emission. The flare 0.3-10 keV peak luminosity decreases with time, following a power-law behaviour with large scatter: L_{pk}~ t_{pk}^{-2.7}. When multiple flares are present, a global softening trend is established: each flare is on average softer than the previous one. The 0.3-10 keV isotropic energy distribution is a log-normal peaked at 10^{51} erg, with a possible excess at low energies. The flare average spectral energy distribution (SED) is found to be a power-law with spectral energy index beta~1.1. These results confirmed that the flares are tightly linked to the prompt emission. However, after considering various models we conclude that no model is currently able to account for the entire set of observations.

Recent particle-in-cell simulations suggest that a large fraction of the energy dissipated in a relativistic shock is deposited into a Maxwellian distribution of electrons that is connected to the high-energy power-law tail. Here, we explore the observational implications of such a mixed thermal-nonthermal particle distribution for the afterglow and prompt emission of gamma-ray bursts. When the Maxwellian component dominates the energy budget, the afterglow lightcurves show a very steep decline phase followed by a more shallow decay when the characteristic synchrotron frequency crosses the observed band. The steep decay appears in the X-rays at ~100 sec after the burst and is accompanied by a characteristic hard-soft-hard spectral evolution that has been observed in a large number of early afterglows. If internal shocks produce a similar mixed electron distribution, a bump is expected at the synchrotron peak of the nu*f_nu spectrum.

The jet of the radio galaxy M87 is misaligned, resulting in a Doppler factor delta~1 for emission of plasma moving parallel to the jet. This makes the observed fast TeV flares on timescales of t_v~5R_g/c harder to understand as emission from the jet. In previous work, we have proposed a jets-in-a-jet model for the ultra-fast TeV flares with t_v<<R_g/c seen in Mrk 501 and PKS 2155-304. Here, we show that about half of the minijets beam their emission outside the jet cone. Minijets emitting off the jet axis result in rapidly evolving TeV (and maybe lower energy) flares that can be observed in nearby radio galaxies. The TeV flaring from M87 fits well into this picture, if M87 is a misaligned blazar.

The fast TeV variability of the blazars Mrk 501 and PKS 2155--304 implies a compact emitting region that moves with a bulk Lorentz factor of Gamma_{em}~100 toward the observer. The Lorentz factor is clearly in excess of the jet Lorentz factors Gamma_j\simless 10 measured on sub-pc scales in these sources. We propose that the TeV emission originates from compact emitting regions that move relativistically {\it within} a jet of bulk Gamma_j~10. This can be physically realized in a Poynting flux-dominated jet. We show that if a large fraction of the luminosity of the jet is prone to magnetic dissipation through reconnection, then material outflowing from the reconnection regions can efficiently power the observed TeV flares through synchrotron-self-Compton emission. The model predicts simultaneous far UV/soft X-ray flares.

Bekenstein's Tensor-Vector-Scalar (TeVeS) theory has had considerable success in explaining various phenomena without the need for dark matter. However, it is difficult to observationally discern the differences between TeVeS and predictions made within the Lambda-cold dark matter concordance model. This implies that alternative tests are required that independently verify which theory is correct. For this we turn to the strong-field regime of TeVeS. In particular, we solve the spherically symmetric equations of hydrostatic equilibrium for a perfect fluid with a realistic equation of state to build models of neutron stars in TeVeS. We show that causality within the neutron star is only maintained for certain cosmological values of the scalar field, which allows us to put constraints on this value independently of cosmological observations. We also discuss in detail the internal structure of neutron stars and how each of the free parameters in the theory effects the overall size and mass of the neutron stars. In particular, the radii of neutron stars in TeVeS can significantly differ from those in General Relativity for certain values of the vector field coupling, which allows us to also place extra constraints on this parameter. Finally, we discuss future observations of neutron stars using both the electromagnetic and gravitational wave spectrums that will allow for tests of the appropriate theory of gravity.

(Abridged) We aim to quantitatively understand the dynamical effect and observational signatures of magnetization of the GRB ejecta on the onset of the afterglow. We perform ultrahigh-resolution one-dimensional relativistic MHD simulations of the interaction of a radially expanding, magnetized ejecta with the interstellar medium. The need of ultrahigh numerical resolution derives from the extreme jump conditions in the region of interaction between the ejecta and the circumburst medium. We study the evolution of an ultrarelativistic shell all the way to a the self-similar asymptotic phase. Our simulations show that the complete evolution can be characterized in terms of two parameters, namely, the \xi parameter introduced by Sari & Piran (1995) and the magnetization \sigma_0. We exploit this property by producing numerical models where the shell Lorentz factor is \gamma_0 ~ tens and rescaling the results to arbitrarily large \gamma_0. We find that the reverse shock is typically very weak or absent for ejecta characterized by \sigma_0 >~ 1. The onset of the forward shock emission is strongly affected by the magnetization. On the other hand, the magnetic energy of the shell is transfered to the external medium on a short timescale (~several times the duration of the burst). The later forward shock emission does not contain information for the initial magnetization of the flow. The asymptotic evolution of strongly magnetized shells, after they have suffred a substantial deceleration, resembles that of hydrodynamic shells, i.e., they fully enter in the Blandford-McKee self-similar regime.

The gamma-ray bursts (GRBs) of long duration are very likely to be connected to the death of massive stars. The gamma-ray emission is believed to come from energy released internally in a flow that moves at ultrarelativistic speed. The fast flow drives a shock wave into the external medium leading to the afterglow emission. Most massive stars form in dense clusters, their high luminosity producing a very dense radiation field. Here, I explore the observational consequences of the interaction of the shocked external medium of the burst with the photon field of a nearby O star. I show that inverse Compton scattering of the stellar photons by electrons heated by the shock leads to powerful gamma-ray emission in the ~1-100 GeV range. This emission appears minutes to hours after the burst and can be easily detected by Cherenkov telescopes and probably with the GLAST satellite. This signal may have already been observed in GRB 940217 and can yield important information about the circumburst environment and the extragalactic background light.

I calculate the emission expected from a Poynting-flux-dominated gamma-ray burst (GRB) flow in which energy is dissipated gradually by magnetic reconnection. In this picture, the energy of the radiating particles is determined by heating and cooling balance (slow heating model). Detailed radiative transfer calculations show that, at Thomson optical depths of order of unity, the dominant radiative process is inverse Compton scattering. Synchrotron-self-absorbed emission and inverse Compton dominate in the Thomson thin parts of the flow. The electrons stay thermal throughout the dissipation region because of Coulomb collisions (Thomson thick part of the flow) and exchange of synchrotron photons (Thomson thin part). The resulting spectrum naturally explains the observed sub-MeV break of the GRB emission and the spectral slopes above and below the break. The model predicts that the gamma-ray power-law tail has a high-energy cutoff typically in the ~0.1-1 GeV energy range that should be observable with {\it GLAST}. The model also predicts a prompt emission component in the optical and UV associated with the GeV emission. Observations of the prompt emission of GRB 061121 that cover the energy range from the optical to ~1 MeV are explained by the model.

Recent observations of two black hole candidates (GX 339-4 and J1753.5-0127) in the low-hard state (L_X/L_Edd ~ 0.003-0.05) suggest the presence of a cool accretion disk very close to the innermost stable orbit of the black hole. This runs counter to models of the low-hard state in which the cool disk is truncated at a much larger radius. We study the interaction between a moderately truncated disk and a hot inner flow. Ion-bombardment heats the surface of the disk in the overlap region between a two-temperature advection-dominated accretion flow and standard accretion disk, producing a hot (kT_e ~70 keV) layer on the surface of the cool disk. The hard X-ray flux from this layer heats the inner parts of the underlying cool disk, producing a soft X-ray excess. Together with interstellar absorption these effects mimic the thermal spectrum from a disk extending to the last stable orbit. The results show that soft excesses in the low-hard state are a natural feature of truncated disk models.

Some recent observational results impose significant constraints on all the models that have been proposed to explain the Galactic black-hole X-ray sources in the hard state. In particular, it has been found that during the hard state of Cyg X-1 the power-law photon number spectral index is correlated with the average time lag between hard and soft X-rays. Furthermore, the peak frequencies of the four Lorentzians that fit the observed power spectra are correlated with both the photon index and the time lag. We performed Monte Carlo simulations of Compton upscattering of soft, accretion-disk photons in the jet and computed the time lag between hard and soft photons and the power-law index of the resulting photon number spectra. We demonstrate that our jet model naturally explains the above correlations, with no additional requirements and no additional parameters.

We study the afterglow phases of a GRB through relativistic magnetohydrodynamic simulations. The evolution of a relativistic shell propagating into a homogeneous external medium is followed. We focus on the effect of the magnetization of the ejecta on the initial phases of the ejecta-external medium interaction. In particular we are studying the condition for the existence of a reverse shock into the ejecta, the timescale for the transfer of the energy from the shell to the shocked medium and the resulting multiwavelength light curves. To this end, we have developed a novel scheme to include non-thermal processeses which is coupled to the relativistic magnetohydrodynamic code MRGENESIS in order to compute the non-thermal synchrotron radiation.

The role of magnetic fields in gamma-ray burst (GRB) flows remains controversial. The study of the early afterglow phases and, in particular, of the reverse shock dynamics and associated emission offers a promising probe of the magnetization of the ejecta. In this paper, we derive the conditions for the existence of a reverse shock in arbitrarily magnetized ejecta that decelerate and interact with the circumburst medium. Both constant and wind-like density profiles are considered. We show, in contrast to previous estimates, that ejecta with magnetization larger than unity are not crossed by a reverse shock for a large fraction of the parameter space relevant to GRB flows. Allowing for shell spreading, there is always a relativistic or mildly relativistic reverse shock forming in sigma_o < 0.3 ejecta. From this, we conclude that the paucity of optical flashes, believed to be a distinctive signature of a reverse shock, may be explained by the existence of dynamically important magnetic fields in the ejecta.

We consider gamma-ray burst outflows with a substantial neutron component that are either dominated by thermal energy (fireballs) or by magnetic energy. In the latter case, we focus on the recently introduced `AC' model which relies on magnetic reconnection to accelerate the flow and power the prompt emission. For both the fireball and the AC model, we investigate the dynamical importance of neutrons on the outflow. We study particle creation in inelastic neutron - proton collisions and find that in both models the resulting neutrino emission is too weak to be detectable. The inelastic collisions also produce gamma-rays, which create pairs in interactions with soft photons carried with the flow. In magnetically driven outflows, the energy of these pairs is radiated away as synchrotron emission. The bulk of the emission takes place at a few hundred keV, which makes it difficult to disentangle this signal from the prompt emission. In fireballs, however, pair cascading leads to the emission of gamma-rays with observer energy in the range of 2 - 20 GeV and a fluence well above the GLAST threshold. Therefore this emission can be a useful diagnostic of the nature of the outflow.

We explore the observational appearance of the photosphere of an ultrarelativistic flow with internal dissipation of energy as predicted by the magnetic reconnection model. Previous study of the radiative transfer in the photospheric region has shown that gradual dissipation of energy results in a hot photosphere. There, inverse Compton scattering of the thermal radiation advected with the flow leads to powerful photospheric emission with spectral properties close to those of the observed prompt GRB emission. Here, we build on that study by calculating the spectra for a large range of the characteristics of the flow. An accurate fitting formula is given that provides the photospheric spectral energy distribution in the ~10 keV to ~10 MeV energy range (in the central engine frame) as a function of the basic physical parameters of the flow. It facilitates the direct comparison of the model predictions with observations, including the variability properties of the lightcurves. We verify that the model naturally accounts for the observed clustering in peak energies of the E*f(E) spectrum. In this model, the Amati relation indicates a tendency for the most luminous bursts to have more energy per baryon. If this tendency also holds for individual GRB pulses, the model predicts the observed narrowing of the width of pulses with increasing photon energy.

For accretion rates Mdot~0.1 Msun/s to a few solar mass black hole the inner part of the disk is expected to make a transition from advection dominance to neutrino cooling. This transition is characterized by sharp changes of the disk properties. I argue here that during this transition, a modest increase of the accretion rate leads to powerful enhancement of the Poynting luminosity of the GRB flow and decrease of its baryon loading. These changes of the characteristics of the GRB flow translate into changing gamma-ray spectra from the photosphere of the flow. The photospheric interpretation of the GRB emission explains the observed narrowing of GRB pulses with increasing photon energy and the luminosity-spectral peak relation within and among bursts.

We apply a novel adaptive mesh refinement code, AMRVAC, to numerically investigate the various evolutionary phases in the interaction of a relativistic shell with its surrounding cold Interstellar Medium (ISM). We do this for both 1D isotropic as well as full 2D jetlike fireball models. This is relevant for Gamma Ray Bursts, and we demonstrate that, thanks to the AMR strategy, we resolve the internal structure of the shocked shell-ISM matter, which will leave its imprint on the GRB afterglow. We determine the deceleration from an initial Lorentz factor $\gamma=100$ up to the almost Newtonian $\gamma\sim{\cal O}(2)$ phase of the flow. We present axisymmetric 2D shell evolutions, with the 2D extent characterized by their initial opening angle. In such jetlike GRB models, we discuss the differences with the 1D isotropic GRB equivalents. These are mainly due to thermally induced sideways expansions of both the shocked shell and shocked ISM regions. We found that the propagating 2D ultrarelativistic shell does not accrete all the surrounding medium located within its initial opening angle. Part of this ISM matter gets pushed away laterally and forms a wide bow-shock configuration with swirling flow patterns trailing the thin shell. The resulting shell deceleration is quite different from that found in isotropic GRB models. As long as the lateral shell expansion is merely due to ballistic spreading of the shell, isotropic and 2D models agree perfectly. As thermally induced expansions eventually lead to significantly higher lateral speeds, the 2D shell interacts with comparably more ISM matter and decelerates earlier than its isotropic counterpart.

I explore the observational appearance of the photosphere of an ultrarelativistic flow with internal dissipation of energy (``dissipative'' GRB models). As a case study, I use the magnetic reconnection model (AC model) that makes robust predictions on the energy dissipation rates at different radii in the flow. With analytical and numerical tools for the radiative transfer problem, I show that the flow develops a hot photosphere where inverse Compton scattering leads to highly non-thermal spectra. For a wide range of luminosities and baryon loadings of the flow, this spectrum is very close to the observed prompt GRB emission. Its luminosity ranges from ~3 to 20% of that of the total energy input.

One of the most intriguing discoveries made by the Swift satellite is the flaring activity in about half of the afterglow lightcurves. Flares have been observed on both long and short duration GRBs and on time scales that range from minutes to ~1 day after the prompt emission. The rapid evolution of some flares led to the suggestion that they are caused by late central engine activity. Here, I propose an alternative explanation that does not need reviving of the central engine. Flares can be powered by delayed magnetic dissipation in strongly magnetized (i.e. with initial Poynting to kinetic flux ratio $\simmore 1$) ejecta during its deceleration due to interaction with the external medium. A closer look at the length scales of the dissipation regions shows that magnetic dissipation can give rise to fast evolving and energetic flares. Multiple flares are also expected in the context of the model.

The role of kink instability in magnetically driven jets is explored through numerical one-dimensional steady relativistic MHD calculations. The instability is shown to have enough time to grow and influence the dynamics of Poynting-flux dominated jets. In the case of AGN jets, the flow becomes kinetic flux dominated at distances larger than ~1000 r_g because of the rapid dissipation of Poynting flux. When applied to GRB outflows, the model predicts more gradual Poynting dissipation and moderately magnetized flow at distances of ~10^{16} cm where the deceleration of the ejecta due to its interaction with the external medium is expected. The energy released by the instability can power the compact ``blazar zone'' emission and the prompt emission of GRB outflows with high radiative efficiencies.

Recently, a relativistic gravitation theory has been proposed [J. D. Bekenstein, Phys. Rev. D {\bf 70}, 083509 (2004)] that gives the Modified Newtonian Dynamics (or MOND) in the weak acceleration regime. The theory is based on three dynamic gravitational fields and succeeds in explaining a large part of extragalactic and gravitational lensing phenomenology without invoking dark matter. In this work we consider the strong gravity regime of TeVeS. We study spherically symmetric, static and vacuum spacetimes relevant for a non-rotating black hole or the exterior of a star. Two branches of solutions are identified: in the first the vector field is aligned with the time direction while in the second the vector field has a non-vanishing radial component. We show that in the first branch of solutions the \beta and \gamma PPN coefficients in TeVeS are identical to these of general relativity (GR) while in the second the \beta PPN coefficient differs from unity violating observational determinations of it (for the choice of the free function $F$ of the theory made in Bekenstein's paper). For the first branch of solutions, we derive analytic expressions for the physical metric and discuss their implications. Applying these solutions to the case of black holes, it is shown that they violate causality (since they allow for superluminal propagation of metric, vector and scalar waves) in the vicinity of the event horizon and/or that they are characterized by negative energy density carried by the fields.

We propose a jet model for the low/hard state of Galactic black-hole X-ray sources that can explain the energy spectra from radio to X-rays. The model assumes that i) there is a magnetic field along the axis of the jet; ii) the electron density in the jet drops inversely proportional to distance; and iii) the electrons in the jet follow a power law distribution function. We have performed Monte Carlo simulations of Compton upscattering of soft photons from the accretion disk and have found power-law high-energy spectra with photon-number index in the range 1.5 - 2 and cutoff at a few hundred keV. The spectrum at long wavelengths (radio, infrared, optical) is modeled to come from synchrotron radiation of the energetic electrons in the jet. We find flat to inverted radio spectra that extend from the radio up to about the optical band. For magnetic field strengths of the order of 10^5-10^6G at the base of the jet, the calculated spectra agree well in slope and flux with the observations. Our model has the advantage over other existing models that it also explains many of the existing timing data such as the time lag spectra, the hardening of the power density spectra and the narrowing of the autocorrelation function with increasing photon energy.

We investigate the production of the gamma-ray spectrum of a Poynting-flux dominated GRB ouflow. The very high magnetic field strengths (super-equipartition) in such a flow lead to very efficient synchrotron emission. In contrast with internal shocks, dissipation of magnetic energy by reconnection is gradual and does not produce the spectrum of cooling electrons associated with shock acceleration.We show that a spectrum with a break in the BATSE energy range is produced, instead, if the magnetic dissipation heats a small (\sim 10^{-4}) population of electrons.

We study possible mechanisms of excitation of quasiperiodic oscillations in the accretion flow of black hole accreters in their hard spectral states, in the context of the `truncated disk' model. Quasi-spherical oscillations of the inner ion-supported accretion flow (ISAF) can be excited by the interaction of this hot flow with the cool disk extending outward from it. The fundamental mode of (p-mode) oscillation is most easily excited, and has a frequency near the Kepler frequency at the inner edge of the cool disk. The strongest excitation mechanism is a feedback loop involving cooling of the ISAF by soft photons from the cool disk and heating of the cool disk by the ISAF, while synchrotron emission can be a relatively strong damping effect. Growth times are computed by detailed Comptonization calculations of the interaction of an idealized ISAF with a cool disk. Typical growth times as short as a few dynamical times are found, while amplitudes can be as large as 10%.

we propose a jet model for the low/hard state of black-hole X-ray sources which explains a) the X-ray spectra, b) the timelag spectra, c) the increase in the amplitude (QPO and high frequency) with increasing photon energy, and d) the narrowing of the autocorrelation function with increasing photon energy. The model (in its simplest form) assumes that i) there is a uniform magnetic field along the axis of the jet, ii) the electron density in the jet is inversely proportional to distance and iii) the jet is "hotter" near its center than at its periphery. We have performed Monte Carlo simulations of Compton upscattering of soft photons from the accretion disk and have found power-law high-energy spectra with photon number index in the range 1.5-2, power-law timelags versus Fourier frequency with index ~0.8, and an increase of the rms amplitude of the variability and a narrowing of the autocorrelation function with photon energy as they have been observed in Cygnus X-1.

Most, probably all, accreting binaries that are believed to contain a black-hole emit radio waves when they are in the low/hard state. Whenever this radio emission has been resolved, a jet-like structure has become apparent. We propose that Compton upscattering of low-energy photons in the jet can explain both the energy spectra and the time lags versus Fourier frequency observed in the low/hard state of black-hole systems. The soft photons originate in the inner part of the accretion disk. We have performed Monte Carlo simulations of Compton upscattering in a jet and have found that for a rather wide range of values of the parameters we can obtain power-law high-energy X-ray spectra with photon-number index in the range 1.5 - 2 and power-law time lags versus Fourier frequency with index ~ 0.7. The black-hole source Cyg X-1 in the low/hard state is well described by our model.