• We present our current best estimate of the plausible observing scenarios for the Advanced LIGO, Advanced Virgo and KAGRA gravitational-wave detectors over the next several years, with the intention of providing information to facilitate planning for multi-messenger astronomy with gravitational waves. We estimate the sensitivity of the network to transient gravitational-wave signals for the third (O3), fourth (O4) and fifth observing (O5) runs, including the planned upgrades of the Advanced LIGO and Advanced Virgo detectors. We study the capability of the network to determine the sky location of the source for gravitational-wave signals from the inspiral of binary systems of compact objects, that is BNS, NSBH, and BBH systems. The ability to localize the sources is given as a sky-area probability, luminosity distance, and comoving volume. The median sky localization area (90\% credible region) is expected to be a few hundreds of square degrees for all types of binary systems during O3 with the Advanced LIGO and Virgo (HLV) network. The median sky localization area will improve to a few tens of square degrees during O4 with the Advanced LIGO, Virgo, and KAGRA (HLVK) network. We evaluate sensitivity and localization expectations for unmodeled signal searches, including the search for intermediate mass black hole binary mergers.
• ### The influence of circumnuclear environment on the radio emission from TDE jets(1605.08437)

Nov. 4, 2016 astro-ph.HE
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.
• ### Numerical simulations of the jetted tidal disruption event Swift J1644+57(1603.05832)

March 18, 2016 astro-ph.HE
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.
• ### The radio afterglow of Swift J1644+57 reveals a powerful jet with fast core and slow sheath(1501.00361)

April 10, 2015 astro-ph.HE
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.
• ### The influence of the magnetic field on the spectral properties of blazars(1310.5441)

Dec. 2, 2013 astro-ph.HE
We explore the signature imprinted by dynamically relevant magnetic fields on the spectral energy distribution (SED) of blazars. It is assumed that the emission from these sources originates from the collision of cold plasma shells, whose magnetohydrodynamic evolution we compute by numerically solving Riemann problems. We compute the SEDs including the most relevant radiative processes and scan a broad parameter space that encompasses a significant fraction of the commonly accepted values of not directly measurable physical properties. We reproduce the standard double hump SED found in blazar observations for unmagnetized shells, but show that the prototype double hump structure of blazars can also be reproduced if the dynamical source of the radiation field is very ultrarelativistic both, in a kinematically sense (namely, if it has Lorentz factors \gtrsim 50) and regarding its magnetization (e.g., with flow magnetizations \sigma \simeq 0.1). A fair fraction of the {\em blazar sequence} could be explained as a consequence of shell magnetization: negligible magnetization in FSRQs, and moderate or large (and uniform) magnetization in BL Lacs. The predicted photon spectral indices (\gph) in the \gamma-ray band are above the observed values (\Gamma_{\rm ph, obs} \lesssim 2.6 for sources with redshifts 0.4\le z \le 0.6) if the magnetization of the sources is moderate (\sigma \simeq 10^{-2}).
• ### Radiative signature of magnetic fields in internal shocks(1111.6108)

Jan. 4, 2012 astro-ph.HE
Common models of blazars and gamma-ray bursts assume that the plasma underlying the ob- served phenomenology is magnetized to some extent. Within this context, radiative signatures of dissipation of kinetic and conversion of magnetic energy in internal shocks of relativistic magnetized outflows are studied. We model internal shocks as being caused by collisions of homogeneous plasma shells. We compute the flow state after the shell interaction by solving Riemann problems at the contact surface between the colliding shells, and then compute the emission from the resulting shocks. Under the assumption of a constant flow luminosity we find that there is a clear difference between the models where both shells are weakly magne- tized ({\sigma}<\sim0.01) and those where, at least, one shell has a {\sigma}>\sim0.01. We obtain that the radiative efficiency is largest for models in which, regardless of the ordering, one shell is weakly and the other strongly magnetized. Substantial differences between weakly and strongly magne- tized shell collisions are observed in the inverse-Compton part of the spectrum, as well as in the optical, X-ray and 1GeV light curves. We propose a way to distinguish observationally between weakly magnetized from magnetized internal shocks by comparing the maximum frequency of the inverse-Compton and synchrotron part of the spectrum to the ratio of the inverse-Compton and synchrotron fluence. Finally, our results suggest that LBL blazars may correspond to barely magnetized flows, while HBL blazars could correspond to moderately magnetized ones. Indeed, by comparing with actual blazar observations we conclude that the magnetization of typical blazars is {\sigma} <\sim 0.01 for the internal shock model to be valid in these sources.
• ### Efficiency of internal shocks in magnetized relativistic jets(1111.3170)

Nov. 15, 2011 astro-ph.HE
We study the dynamic and radiative efficiency of conversion of kinetic-to-thermal/magnetic energy by internal shocks in relativistic magnetized outflows. A parameter study of a large number of collisions of cylindrical shells is performed. We explore how, while keeping the total flow luminosity constant, the variable fluid magnetization influences the efficiency and find that the interaction of shells in a mildly magnetized jet yields higher dynamic, but lower radiative efficiency than in a non-magnetized flow. A multi-wavelength radiative signature of different shell magnetization is computed assuming that relativistic particles are accelerated at internal shocks.
• ### An unusual stellar death on Christmas Day(1105.3015)

Oct. 3, 2011 astro-ph.HE
Long Gamma-Ray Bursts (GRBs) are the most dramatic examples of massive stellar deaths, usually as- sociated with supernovae (Woosley et al. 2006). They release ultra-relativistic jets producing non-thermal emission through synchrotron radiation as they interact with the surrounding medium (Zhang et al. 2004). Here we report observations of the peculiar GRB 101225A (the "Christmas burst"). Its gamma-ray emission was exceptionally long and followed by a bright X-ray transient with a hot thermal component and an unusual optical counterpart. During the first 10 days, the optical emission evolved as an expanding, cooling blackbody after which an additional component, consistent with a faint supernova, emerged. We determine its distance to 1.6 Gpc by fitting the spectral-energy distribution and light curve of the optical emission with a GRB-supernova template. Deep optical observations may have revealed a faint, unresolved host galaxy. Our proposed progenitor is a helium star-neutron star merger that underwent a common envelope phase expelling its hydrogen envelope. The resulting explosion created a GRB-like jet which gets thermalized by interacting with the dense, previously ejected material and thus creating the observed black-body, until finally the emission from the supernova dominated. An alternative explanation is a minor body falling onto a neutron star in the Galaxy (Campana et al. 2011).
• ### Variable Ly alpha sheds light on the environment surrounding GRB 090426(1101.3488)

Jan. 18, 2011 astro-ph.HE
Long duration gamma-ray bursts are commonly associated with the deaths of massive stars. Spectroscopic studies using the afterglow as a light source provide a unique opportunity to unveil the medium surrounding it, probing the densest region of their galaxies. This material is usually in a low ionisation state and at large distances from the burst site, hence representing the normal interstellar medium in the galaxy. Here we present the case of GRB 090426 at z=2.609, whose optical spectrum indicates an almost fully ionised medium together with a low column density of neutral hydrogen. For the first time, we also observe variations in the Ly alpha absorption line. Photoionisation modeling shows that we are probing material from the vicinity of the burst (~80 pc). The host galaxy is a complex of two luminous interacting galaxies, which might suggest that this burst could have occurred in an isolated star-forming region outside its host galaxy created in the interaction of the two galaxies.
• ### Searching for differences in Swift's intermediate GRBs(1006.4469)

July 15, 2010 astro-ph.CO, astro-ph.HE
Gamma-ray bursts are usually classified through their high-energy emission into short-duration and long-duration bursts, which presumably reflect two different types of progenitors. However, it has been shown on statistical grounds that a third, intermediate population is needed in this classification scheme, although an extensive study of the properties of this class has so far not been done. The large amount of follow-up studies generated during the Swift era allows us to have a suficient sample to attempt a study of this third population through the properties of their prompt emission and their afterglows. Our study is focused on a sample of GRBs observed by Swift during its first four years of operation. The sample contains those bursts with measured redshift since this allows us to derive intrinsic properties. Intermediate bursts are less energetic and have dimmer afterglows than long GRBs, especially when considering the X-ray light curves, which are on average one order of magnitude fainter than long bursts. There is a less significant trend in the redshift distribution that places intermediate bursts closer than long bursts. Except for this, intermediate bursts show similar properties to long bursts. In particular, they follow the Epeak vs. Eiso correlation and have, on average, positive spectral lags with a distribution similar to that of long bursts. Like long GRBs, they normally have an associated supernova, although some intermediate bursts have shown no supernova component. This study shows that intermediate bursts are different from short bursts and, in spite of sharing many properties with long bursts, there are some differences between them as well. We suggest that the physical difference between intermediate and long bursts could be that for the first the ejecta are thin shells while for the latter they are thick shells.
• ### Multiwavelength afterglow light curves from magnetized GRB flows(1004.2720)

April 15, 2010 astro-ph.HE
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.
• ### Detection of the high z GRB 080913 and its implications on progenitors and energy extraction mechanisms(0810.2107)

Dec. 17, 2009 astro-ph
Aims: We present multiwavelength observations of one of the most distant gamma-ray bursts detected so far, GRB 080913. Based on these observations, we consider whether it could be classified as a short-duration GRB and discuss the implications for the progenitor nature and energy extraction mechanisms. Methods: Multiwavelength X-ray, near IR and millimetre observations were made between 20.7 hours and 16.8 days after the event. Results: Whereas a very faint afterglow was seen at the 3.5m CAHA telescope in the nIR, the X-ray afterglow was clearly detected in both Swift and XMM-Newton observations. An upper limit is reported in the mm range. We have modeled the data assuming a collimated $\theta_0$ $\gtrsim$ 3$^\circ$ blast wave with an energy injection at 0.5 days carrying $5\sim 10^{52}$ erg or approximately 12 times the initial energy of the blast wave. We find that GRB 080913 shares many of the gamma-ray diagnostics with the more recent burst GRB 090423 for being classified as short had they ocurred at low redshift. If the progenitor were a compact binary merger, it is likely composed by a NS and BH. The Blandford-Znajek (BZ) mechanism is the preferred one to extract energy from the central, maximally-rotating BH. Both the magnetic field close to the event horizon (B) and the BH mass ($M_{bh}$) are restricted within a relatively narrow range, such that $(B / 3\times 10^{16} \rm{G}) (M_{bh} / 7 M_\odot) \sim 1$. Similar constraints on the central BH hold for collapsar-like progenitor systems if the BZ-mechanism works for the system at hand.
• ### Observational Effects of Anomalous Boundary Layers in Relativistic Jets(0803.2693)

March 18, 2008 astro-ph
Recent theoretical work has pointed out that the transition layer between a jet an the medium surrounding it may be more complex than previously thought. Under physically realizable conditions, the transverse profile of the Lorentz factor in the boundary layer can be non-monotonic, displaying the absolute maximum where the flow is faster than at the jet spine, followed by an steep fall off. Likewise, the rest-mass density, reaches an absolute minimum (coincident with the maximum in Lorentz factor) and then grows until it reaches the external medium value. Such a behavior is in contrast to the standard monotonic decline of the Lorentz factor (from a maximum value at the jet central spine) and the corresponding increase of the rest-mass density (from the minimum reached at the jet core). We study the emission properties of the aforementioned anomalous shear layer structures in kiloparsec-scale jets aiming to show observable differences with respect to conventional monotonic and smooth boundary layers.
• ### The Influence of Helical Magnetic Fields in the Dynamics and Emission of Relativistic Jets(0802.3826)

Feb. 26, 2008 astro-ph
We present numerical relativistic magnetohydrodynamic and emission simulations aimed to study the role played by the magnetic field in the dynamics and emission of relativistic jets in Active Galactic Nuclei. We focus our analysis on the study of the emission from recollimation shocks since they may provide an interpretation for the stationary components seen at parsec-scales in multiple sources. We show that the relative brightness of the knots associated with the recollimation shocks decreases with increasing jet magnetization, suggesting that jets presenting stationary components may have a relatively weak magnetization, with magnetic fields of the order of equipartition or below.
• ### On the existence of a reverse shock in magnetized GRB ejecta(0711.1980)

Nov. 30, 2007 astro-ph
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.
• ### Relativistic Outflows in Gamma-Ray Bursts(astro-ph/0701187)

Jan. 8, 2007 astro-ph
The possibility that gamma-ray bursts (GRBs) were not isotropic emissions was devised theoretically as a way to ameliorate the huge energetic budget implied by the standard fireball model for these powerful phenomena. However, the mechanism by which after the quasy-isotropic release of a few $10^{50}$erg yields a collimated ejection of plasma could not be satisfactory explained analytically. The reason being that the collimation of an outflow by its progenitor system depends on a very complex and non-linear dynamics. That has made necessary the use of numerical simulations in order to shed some light on the viability of some likely progenitors of GRBs. In this contribution I will review the most relevant features shown by these numerical simulations and how they have been used to validate the collapsar model (for long GRBs) and the model involving the merger of compact binaries (for short GRBs).
• ### High--Resolution 3D Simulations of Relativistic Jets(astro-ph/9906428)

June 26, 1999 astro-ph
We have performed high-resolution 3D simulations of relativistic jets with beam flow Lorentz factors up to 7, a spatial resolution of 8 cells per beam radius, and for up to 75 normalized time units to study the morphology and dynamics of 3D relativistic jets. Our simulations show that the coherent fast backflows found in axisymmetric models are not present in 3D models. We further find that when the jet is exposed to non-axisymmetric perturbations, (i) it does not display the strong perturbations found for 3D classical hydrodynamic and MHD jets (at least during the period of time covered by our simulations), and (ii) it does propagate according to the 1D estimate. Small 3D effects in the relativistic beam give rise to a lumpy distribution of apparent speeds like that observed in M87. The beam is surrounded by a boundary layer of high specific internal energy. The properties of this layer are briefly discussed.
• ### GENESIS: A high-resolution code for 3D relativistic hydrodynamics(astro-ph/9903352)

March 23, 1999 astro-ph
The main features of a three dimensional, high-resolution special relativistic hydro code based on relativistic Riemann solvers are described. The capabilities and performance of the code are discussed. In particular, we present the results of extensive test calculations which demonstrate that the code can accurately and efficiently handle strong shocks in three spatial dimensions. Results of the performance of the code on single and multi-processor machines are given. Simulations (in double precision) with 7x10^6 computational cells require less than 1 Gb of RAM memory and ~7x10^-5 CPU seconds per zone and time step (on a SCI Cray-Origin 2000 with a R10000 processor). Currently, a version of the numerical code is under development, which is suited for massively parallel computers with distributed memory architecture (like, e.g., Cray T3E).