Among the potential gravitational wave (GW) sources for the upcoming space-based interferometer LISA are extreme- or intermediate-mass ratio inspirals (EMRI/IMRIs). These events involve the coalescence of supermassive black holes in the mass range $10^5M_{\odot}\lesssim M \lesssim 10^7M_{\odot}$ with companion BHs of much lower masses. A subset of E/IMRIs are expected to occur in the accretion discs of active galactic nuclei (AGN). Previous work has shown that torques exerted by the disc can interfere with the inspiral and cause a phase shift in the GW waveform detectable by LISA. Here we use a suite of two-dimensional hydrodynamical simulations with the moving-mesh code DISCO to present a systematic study of disc torques. We measure torques on an inspiraling BH and compute the corresponding waveform deviations as a function of the binary mass ratio $q\equiv M_2/M_1$, the disc viscosity ($\alpha$), and gas temperature (or equivalently Mach number; $\mathcal{M}$). We find that the absolute value of the gas torques is within an order of magnitude of previously determined planetary migration torques, but their precise value and sign depends non-trivially on the combination of these parameters, the inspiral rate, and the accretion rate onto the satellite BH. The gas imprint is generally detectable by LISA for binaries embedded in AGN discs with surface densities above $\Sigma_0\ge10^{4-6} \rm \, g \, cm^{-2}$, depending on $q$, $\alpha$ and $\mathcal{M}$. We find that the deviations are most pronounced in discs with higher viscosities, and for E/IMRIs detected at frequencies where LISA is most sensitive. Torques in colder discs exhibit a noticeable dependence on the GW-driven inspiral rate as well as strong fluctuations at late stages of the inspiral. Our results further suggest that LISA may be able to place constraints on AGN disc parameters and the physics of disc-satellite interaction.

Numerical studies of gas accretion onto supermassive black hole binaries (SMBHBs) have generally been limited to conditions where the circumbinary disk (CBD) is 10-100 times thicker than expected for disks in active galactic nuclei (AGN). This discrepancy arises from technical limitations, and also from publication bias toward replicating fiducial numerical models. Here we present the first systematic study of how the binary's orbital evolution varies with disk scale height. We report three key results: (1) Binary orbital evolution switches from outspiralling for warm disks (aspect ratio ~0.1), to inspiralling for more realistic cooler, thinner disks at a critical aspect ratio ~0.04, corresponding to orbital Mach number ~25. (2) The net torque on the binary arises from a competition between positive torque from gas orbiting close to the black holes, and negative torque from the inner edge of the CBD, which is denser for thinner disks. This leads to increasingly negative net torques on the binary for increasingly thin disks. (3) The accretion rate is modestly suppressed with increasing Mach number. We discuss how our results may influence modeling of the nano-Hz gravitational wave background, as well as estimates of the LISA merger event rate.

Using numerical hydrodynamics calculations and a novel method for densely sampling parameter space, we measure the accretion and torque on a binary system from a circumbinary disk. In agreement with previous studies, we find that the net torque on the binary is positive for mass ratios close to unity, and that accretion always drives the binary towards equal mass. Accretion variability depends sensitively on the numerical sink prescription, but the torque and relative accretion onto each component do not depend on the sink timescale. Positive torque and highly variable accretion occurs only for mass ratios greater than around $0.05$. This means that for mass ratios below $0.05$, the binary would migrate inward until the secondary accreted sufficient mass, after which it would execute a U-turn and migrate outward. We explore a range of viscosities, from $\alpha = 0.03$ to $\alpha = 0.15$, and find that this outward torque is proportional to the viscous torque, simply proportional to viscosity in this range. Dependence of accretion and torque on mass ratio is explored in detail, densely sampling mass ratios between $0.01$ and unity. For mass ratio $q > 0.6$, accretion variability is found to exhibit a distinct sawtooth pattern, typically with a five-orbit cycle that provides a "smoking gun" prediction for variable quasars observed over long periods, as a potential means to confirm the presence of a binary.

The coalescence of a compact object with a $10^{4}-10^{7} {\rm M_\odot}$ supermassive black hole (SMBH) produces mHz gravitational waves (GWs) detectable by the future Laser Interferometer Space Antenna (LISA). If such an inspiral occurs in the accretion disc of an active galactic nucleus (AGN), the gas torques imprint a small deviation in the GW waveform. Here we present two-dimensional hydrodynamical simulations with the moving-mesh code DISCO of a BH inspiraling at the GW rate in a binary system with a mass ratio $q\!=\!M_2/M_1\!=\!10^{-3}$, embedded in an accretion disc. We assume a locally isothermal equation of state for the gas (with Mach number $\mathcal{M}=20$) and implement a standard $\alpha$-prescription for its viscosity (with $\alpha = 0.03$). We find disc torques on the binary that are weaker than in previous semi-analytic toy models, and are in the opposite direction: the gas disc slows down, rather than speeds up the inspiral. We compute the resulting deviations in the GW waveform, which scale linearly with the mass of the disc. The SNR of these deviations accumulates mostly at high frequencies, and becomes detectable in a 5-year LISA observation if the total phase shift exceeds a few radians. We find that this occurs if the disc surface density exceeds $\Sigma_0 \gtrsim 10^{2-3}\rm g\,cm^{-2}$, as may be the case in thin discs with near-Eddington accretion rates. Since the characteristic imprint on the GW signal is strongly dependent on disc parameters, a LISA detection of an intermediate mass ratio inspiral would probe the physics of AGN discs and migration.

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

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

We systematically investigate the outflow structure of GW170817 in comparison with a sample of 27 cosmological short GRBs by modelling their afterglow light curves. We find that cosmological short GRBs share the same outflow structures with GW170817, relativistic structured jets. The jet opening angle of GW170817 is $6.3^{+1.1}_{-0.6}{}^{\circ}$, which is consistent with that of cosmological short GRBs ($\theta_0 = 6.9^\circ \pm 2.3^{\circ}$). Our analysis indicates that GW170817 is viewed off-axis ($\theta_{\rm obs} = 30^{+7}_{-4}{}^{\circ}$), while cosmological short GRBs are viewed on-axis ($\theta_{\rm obs} \lesssim \theta_0$). The exceptional properties of the GW170817 afterglow can be explained by the difference in observation angle alone. We demonstrate that the light curves of the GW170817 afterglow, if viewed on-axis, are consistent with those of cosmological short GRBs. Other properties of GW170817, such as Lorentz factor $\Gamma \approx 150$, spectral index $p \approx 2.15$, isotropic equivalent energy $E_{\rm iso} \approx 8 \times 10^{52}$ erg and interstellar medium density $n_0 \approx 10^{-2}$ proton cm$^{-3}$, fit well within the ranges of those of cosmological short GRBs. The similarity between the GW170817 outflow structure and those of cosmological short GRBs indicates that cosmological short GRBs are likely neutron star mergers.

We present full-time-domain, moving-mesh, relativistic hydrodynamic simulations of jets launched from the center of a massive progenitor star and compute the resulting synchrotron light curves for observers at a range of viewing angles. We follow jet evolution from ignition inside the stellar center, propagation in the stellar envelope and breakout from the stellar surface, then through the coasting and deceleration phases. The jet compresses into a thin shell, sweeps up the circumstellar medium, and eventually enters the Newtonian phase. The jets naturally develop angular and radial structure due to hydro-dynamical interaction with surrounding gas. The calculated synchrotron light curves cover the observed temporal range of prompt to late afterglow phases of long gamma-ray bursts (LGRBs). The on-axis light curves exhibit an early emission pulse originating in shock-heated stellar material, followed by a shallow decay and a later steeper decay. The off-axis light curves rise earlier than previously expected for top-hat jet models -- on a time scale of seconds to minutes after jet breakout, and decay afterwards. Sometimes the off-axis light curves have later re-brightening components that can be contemporaneous with SNe Ic-bl emission. Our calculations may shed light on the structure of GRB outflows in the afterglow stage. The off-axis light curves from full-time-domain simulations advocate new light curve templates for the search of off-axis/orphan afterglows.

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.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

We present the first extensive radio to gamma-ray observations of a fast-rising blue optical transient (FBOT), AT2018cow, over its first ~100 days. AT2018cow rose over a few days to a peak luminosity $L_{pk}\sim4\times 10^{44}$ erg/s exceeding those of superluminous supernovae (SNe), before declining as $\propto t^{-2}$. Initial spectra at $\lesssim 15$ days were mostly featureless and indicated large expansion velocities v~0.1c and temperatures reaching 30000 K. Later spectra revealed a persistent optically-thick photosphere and the emergence of H and He emission features with v~sim 4000 km/s with no evidence for ejecta cooling. Our broad-band monitoring revealed a hard X-ray spectral component at $E\ge 10$ keV, in addition to luminous and highly variable soft X-rays, with properties unprecedented among astronomical transients. An abrupt change in the X-ray decay rate and variability appears to accompany the change in optical spectral properties. AT2018cow showed bright radio emission consistent with the interaction of a blastwave with $v_{sh}$~0.1c with a dense environment ($\dot M\sim10^{-3}-10^{-4}\,M_{\odot}yr^{-1}$ for $v_w=1000$ km\s). While these properties exclude Ni-powered transients, our multi-wavelength analysis instead indicates that AT2018cow harbored a "central engine", either a compact object (magnetar or black hole) or an embedded internal shock produced by interaction with a compact, dense circumstellar medium. The engine released $\sim10^{50}-10^{51.5}$ erg over $\sim10^3-10^5$ s and resides within low-mass fast-moving material with equatorial-polar density asymmetry ($M_{ej,fast}\lesssim0.3\,\rm{M_{\odot}}$). Successful SNe from low-mass H-rich stars (like electron-capture SNe) or failed explosions from blue supergiants satisfy these constraints. Intermediate-mass black-holes are disfavored by the large environmental density probed by the radio observations.

The multi-wavelength non-thermal emission from the binary neutron star (BNS) merger GW170817 has raised a heated debate concerning the post-merger outflow structure. Both a relativistic structured jet viewed off-axis and a mildly relativistic quasi-spherical outflow can explain the observational data of GW170817 up to ~260 days. We utilize a physically motivated analytic two-parameter model called the "boosted fireball" for the outflow structure after it has expanded far from the merger site and has entered the self-similar coasting phase. This model consists of a family of outflows with a structure varying smoothly between a highly collimated ultra-relativistic jet and an isotropic fireball. We simulate the dynamical evolution, starting with "boosted fireball" initial conditions, of 240 outflows using the moving-mesh relativistic hydrodynamics code JET to follow their evolution through the afterglow phase. We compute nearly 2,000,000 synchrotron spectra from the hydrodynamic simulations using the standard synchrotron radiation model. By making use of scaling relations in the hydrodynamic and radiation equations, we develop a synthetic light curve generator with an efficient sampling speed. This allows us to fit the observational data by performing Markov-Chain Monte Carlo (MCMC) analysis in a 8-dimensional parameter space, consisting of hydrodynamic parameters, radiation parameters and observational parameters. Our results favor the relativistic structured jet, with a jet opening angle ~5 deg and Lorentz factor ~175, viewed from an off-axis angle of 27(+9-3) deg. Due to parameter degeneracies, we find broad distributions for the explosion energy E_0, the circumburst density n_0, epsilon_e and epsilon_B. The combination of a high n_0 and a low epsilon_B can also produce a good fit, indicating that an extremely low n_0 may not be required for GW170817.

We present synthetic radio images of the GRB170817A afterglow, computed from moving-mesh hydrodynamic simulations of binary neutron star merger outflows. Having expanded for nearly a year, the merger remnant is expected to subtend roughly 5 milli-arcseconds (mas) on the sky, potentially resolvable by very long baseline radio imaging techniques. Any observations revealing the radio centroid to be offset from the line-of-site to the merger would be the smoking gun of a jetted outflow. However, our results indicate that a measurement of the centroid position alone cannot independently determine whether that jet escaped successfully from the merger debris cloud, or was "choked," yielding a quasi-spherical explosion. We find that in both scenarios, the centroid exhibits superluminal proper motion away from the merger site at roughly 4 - 10 mas per day for at least the first 300 days. We argue that a successful strategy for differentiating among the explosion models will need to include multiple observations over the coming months - years. In particular, we find the time at which the centroid attains its maximum offset, and begins heading back toward the merger site, is considerably later if the jet was choked. Detecting a reversal of the centroid trajectory earlier than 600 days would uniquely identify a successful jet. Our results indicate the source might be resolved using VLBI radio observing techniques with roughly 1 mas resolution starting at roughly 400 days post-merger, and that the the angular extent of a successful jet is significantly smaller than that of a choked jet (4.5 versus 7 mas respectively).

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

We present numerical simulations of energetic flows propagating through the debris cloud of a binary neutron star (BNS) merger. Starting from the scale of the central engine, we use a moving-mesh hydrodynamics code to simulate the complete dynamical evolution of the produced relativistic jets. We compute synchrotron emission directly from the simulations and present multi-band light curves of the early (sub-day) through late (weeks to years) afterglow stages. Our work systematically compares two distinct models for the central engine, referred to as the narrow and wide engine scenario, which is associated with a successful structured jet and a quasi-isotropic explosion respectively. Both engine models naturally evolve angular and radial structure through hydrodynamical interaction with the merger debris cloud. They both also result in a relativistic blast wave capable of producing the observed multi-band afterglow data. However, we find that the narrow and wide engine scenario might be differentiated by a new emission component that we refer to as a merger flash. This component is a consequence of applying the synchrotron radiation model to the shocked optically thin merger cloud. Such modeling is appropriate if injection of non-thermal electrons is sustained in the breakout relativistic shell, for example by internal shocks or magnetic reconnection. The rapidly declining signature may be detectable for future BNS mergers during the first minutes to day following the GW chirp. Furthermore, its non-detection for the GRB170817A event may disfavor the wide, quasi-isotropic explosion model.

THESEUS is a space mission concept aimed at exploiting Gamma-Ray Bursts for investigating the early Universe and at providing a substantial advancement of multi-messenger and time-domain astrophysics. These goals will be achieved through a unique combination of instruments allowing GRB and X-ray transient detection over a broad field of view (more than 1sr) with 0.5-1 arcmin localization, an energy band extending from several MeV down to 0.3 keV and high sensitivity to transient sources in the soft X-ray domain, as well as on-board prompt (few minutes) follow-up with a 0.7 m class IR telescope with both imaging and spectroscopic capabilities. THESEUS will be perfectly suited for addressing the main open issues in cosmology such as, e.g., star formation rate and metallicity evolution of the inter-stellar and intra-galactic medium up to redshift $\sim$10, signatures of Pop III stars, sources and physics of re-ionization, and the faint end of the galaxy luminosity function. In addition, it will provide unprecedented capability to monitor the X-ray variable sky, thus detecting, localizing, and identifying the electromagnetic counterparts to sources of gravitational radiation, which may be routinely detected in the late '20s / early '30s by next generation facilities like aLIGO/ aVirgo, eLISA, KAGRA, and Einstein Telescope. THESEUS will also provide powerful synergies with the next generation of multi-wavelength observatories (e.g., LSST, ELT, SKA, CTA, ATHENA).

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.

We present the results of 2D, moving-mesh, viscous hydrodynamical simulations of an accretion disc around a merging supermassive black hole binary (SMBHB). The simulation is pseudo-Newtonian, with the BHs modeled as point masses with a Paczynski-Wiita potential, and includes viscous heating, shock heating, and radiative cooling. We follow the gravitational inspiral of an equal-mass binary with a component mass M_bh=10^6 M_sun from an initial separation of 60r_g (where r_g=GM_bh/c^2 is the gravitational radius) to the merger. We find that a central, low-density cavity forms around the binary, as in previous work, but that the BHs capture gas from the circumbinary disc and accrete efficiently via their own minidiscs, well after their inspiral outpaces the viscous evolution of the disc. The system remains luminous, displaying strong periodicity at twice the binary orbital frequency throughout the entire inspiral process, all the way to the merger. In the soft X-ray band, the thermal emission is dominated by the inner edge of the circumbinary disc with especially clear periodicity in the early inspiral. By comparison, harder X-ray emission is dominated by the minidiscs, and the light curve is initially more noisy but develops a clear periodicity in the late inspiral stage. This variability pattern should help identify the EM counterparts of SMBHBs detected by the space-based gravitational-wave detector LISA.

Unusually high velocities (< ~0.1c) and correspondingly high kinetic energies have been observed in a subset of Type Ic supernovae (so-called "broad-lined Ic" supernovae; SNe Ic-BL), prompting a search for a central engine model capable of generating such energetic explosions. A clue to the explosion mechanism may lie in the fact that all supernovae that accompany long-duration gamma-ray bursts belong to the SN Ic-BL class. Using a combination of two-dimensional relativistic hydrodynamics and radiation transport calculations, we demonstrate that the central engine responsible for long gamma-ray bursts can also trigger a SN Ic-BL. We find that a reasonable gamma-ray burst engine injected into a stripped Wolf-Rayet progenitor produces a relativistic jet with energy ~10^51 ergs, as well as a SN whose synthetic light curves and spectra are fully consistent with observed SNe Ic-BL during the photospheric phase. As a result of the jet's asymmetric energy injection, the SN spectra and light curves depend on viewing angle. The impact of viewing angle on the spectrum is particularly pronounced at early times, while the viewing angle dependence for the light curves (~10% variation in bolometric luminosity) persists throughout the photospheric phase.

Gaseous circumbinary accretion discs provide a promising mechanism to facilitate the mergers of supermassive black holes (SMBHs) in galactic nuclei. We measure the torques exerted on accreting SMBH binaries, using 2D, isothermal, moving-mesh, viscous hydrodynamical simulations of circumbinary accretion discs. Our computational domain includes the entire inner region of the circumbinary disk with the individual black holes (BHs) included as point masses on the grid and a sink prescription to model accretion onto each BH. The BHs each acquire their own well-resolved accretion discs ("minidiscs"). We explore a range of mass removal rates for the sink prescription removing gas from the central regions of the minidiscs. We find that the torque exerted on the binary is primarily gravitational, and dominated by the gas orbiting close behind and ahead of the individual BHs. The torques from the distorted circumbinary disc farther out and from the direct accretion of angular momentum are subdominant. The torques are sensitive to the sink prescription: slower sinks result in more gas accumulating near the BHs and more negative torques, driving the binary to merger more rapidly. For faster sinks, the torques are less negative and eventually turn positive (for unphysically fast sinks). When the minidiscs are modeled as standard alpha discs, our results are insensitive to the choice of sink radius. Scaling the simulations to a binary orbital period tbin = 1yr and background disc accretion rate Mdot = 0.3MEdd in Eddington units, the binary inspirals on a timescale of 3X10^6 years, irrespective of the SMBH masses. For binaries with total mass <10^7Msun, this is shorter than the inspiral time due to gravitational wave (GW) emission alone, implying that gas discs will have a significant impact on the SMBH binary population and can affect the GW signal for Pulsar Timing Arrays.

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

We study circumbinary accretion discs in the framework of the restricted three-body problem (R3Bp) and via numerically solving the height-integrated equations of viscous hydrodynamics. Varying the mass ratio of the binary, we find a pronounced change in the behaviour of the disc near mass ratio $q \equiv M_s/M_p \sim 0.04$. For mass ratios above $q=0.04$, solutions for the hydrodynamic flow transition from steady, to strongly-fluctuating; a narrow annular gap in the surface density around the secondary's orbit changes to a hollow central cavity; and a spatial symmetry is lost, resulting in a lopsided disc. This phase transition is coincident with the mass ratio above which stable orbits do not exist around the L4 and L5 equilibrium points of the R3B problem. Using the DISCO code, we find that for thin discs, for which a gap or cavity can remain open, the mass ratio of the transition is relatively insensitive to disc viscosity and pressure. The $q=0.04$ transition has relevance for the evolution of massive black hole binary+disc systems at the centers of galactic nuclei, as well as for young stellar binaries and possibly planets around brown dwarfs.

The merger of binary neutron stars (BNSs) can lead to large amplifications of the magnetic field due to the development of turbulence and instabilities in the fluid, such as the Kelvin-Helmholtz shear instability, which drive small-scale dynamo activity. In order to properly resolve such instabilities and obtain the correct magnetic field amplification, one would need to employ resolutions that are currently unfeasible in global general relativistic magnetohydrodynamic (GRMHD) simulations of BNS mergers. Here, we present a subgrid model that allows global simulations to take into account the small-scale amplification of the magnetic field which is caused by the development of turbulence during BNS mergers. Assuming dynamo saturation, we show that magnetar-level fields ($\sim 10^{16}\,{\rm G}$) can be easily reached, and should therefore be expected from the merger of magnetized BNSs. The total magnetic energy can reach values up to $\sim 10^{51}\,{\rm erg}$ and the post-merger remnant can therefore emit strong electromagnetic signals and possibly produce short gamma-ray bursts.

Graham et al. (2015) have detected a 5.2 year periodic optical variability of the quasar PG 1302-102 at redshift $z=0.3$, which they interpret as the redshifted orbital period $(1+z)t_{\rm bin}$ of a putative supermassive black hole binary (SMBHB). Here we consider the implications of a $3-8$ times shorter orbital period, suggested by hydrodynamical simulations of circumbinary discs (CBDs) with nearly equal--mass SMBHBs ($q\equiv M_2/M_1\gtrsim 0.3$). With the corresponding $2-4$ times tighter binary separation, PG 1302 would be undergoing gravitational wave dominated inspiral, and serve as a proof that the BHs can be fueled and produce bright emission even in this late stage of the merger. The expected fraction of binaries with the shorter $t_{\rm bin}$, among bright quasars, would be reduced by 1-2 orders of magnitude, compared to the 5.2 year period, in better agreement with the rarity of candidates reported by Graham et al. (2015). Finally, shorter periods would imply higher binary speeds, possibly imprinting periodicity on the light curves from relativistic beaming, as well as measurable relativistic effects on the Fe K$\alpha$ line. The CBD model predicts additional periodic variability on timescales of $t_{\rm bin}$ and $\approx 0.5 t_{\rm bin}$, as well as periodic variation of broad line widths and offsets relative to the narrow lines, which are consistent with the observations. Future observations will be able to test these predictions and hence the binary+CBD hypothesis for PG 1302.

We use two-dimensional relativistic hydrodynamic numerical calculations to show that highly collimated relativistic jets can be produced in neutron star merger models of short-duration gamma ray bursts without the need for a highly directed engine or a large net magnetic flux. Even a hydrodynamic engine generating a very wide sustained outflow on small scales can in principle produce a highly collimated relativistic jet, facilitated by a dense surrounding medium which provides a cocoon surrounding the jet core. An oblate geometry to the surrounding gas significantly enhances the collimation process. Previous numerical simulations have shown that the merger of two neutron stars produces an oblate, expanding cloud of dynamical ejecta. We show that this gas can efficiently collimate the central engine power much as the surrounding star does in long-duration GRB models. For typical short-duration GRB central engine parameters, we find jets with opening angles of order 10 degrees in which a large fraction of the total outflow power of the central engine resides in highly relativistic material. These results predict large differences in the opening angles of outflows from binary neutron star mergers versus neutron star-black hole mergers.

We demonstrate that the steep decay and long plateau in the early phases of gamma ray burst (GRB) X-ray afterglows are naturally produced in the collapsar model, by a means ultimately related to the dynamics of relativistic jet propagation through a massive star. We present two-dimensional axisymmetric hydrodynamical simulations which start from a collapsar engine and evolve all the way through the late afterglow phase. The resultant outflow includes a jet core which is highly relativistic after breaking out of the star, but becomes baryon-loaded after colliding with a massive outer shell, corresponding to mass from the stellar atmosphere of the progenitor star which became trapped in front of the jet core at breakout. The prompt emission produced before or during this collision would then have the signature of a high Lorentz factor jet, but the afterglow is produced by the amalgamated post-collision ejecta which has more inertia than the original highly relativistic jet core and thus has a delayed deceleration. This naturally explains the early light curve behavior discovered by Swift, including a steep decay and a long plateau, without invoking late-time energy injection from the central engine. The numerical simulation is performed continuously from engine to afterglow, covering a dynamic range of over ten orders of magnitude in radius. Light curves calculated from the numerical output demonstrate that this mechanism reproduces basic features seen in early afterglow data. Initial steep decays are produced by internal shocks, and the plateau corresponds to the coasting phase of the outflow.

Gamma Ray Bursts (GRBs) are characterized by ultra-relativistic outflows, while supernovae are generally characterized by non-relativistic ejecta. GRB afterglows decelerate rapidly usually within days, because their low-mass ejecta rapidly sweep up a comparatively larger mass of circumstellar material. However supernovae, with heavy ejecta, can be in nearly free expansion for centuries. Supernovae were thought to have non-relativistic outflows except for few relativistic ones accompanied by GRBs. This clear division was blurred by SN 2009bb, the first supernova with a relativistic outflow without an observed GRB. Yet the ejecta from SN 2009bb was baryon loaded, and in nearly-free expansion for a year, unlike GRBs. We report the first supernova discovered without a GRB, but with rapidly decelerating mildly relativistic ejecta, SN 2012ap. We discovered a bright and rapidly evolving radio counterpart driven by the circumstellar interaction of the relativistic ejecta. However, we did not find any coincident GRB with an isotropic fluence of more than a sixth of the fluence from GRB 980425. This shows for the first time that central engines in type Ic supernovae, even without an observed GRB, can produce both relativistic and rapidly decelerating outflows like GRBs.

We describe a technique for improving the performance of hydrodynamics codes which employ a moving Voronoi mesh. Currently, such codes are susceptible to high-frequency noise produced by rapid adjustments in the grid topology on the smallest scales. The treatment for this grid noise is simple; instead of moving the mesh-generating marker points with the local fluid velocity, this velocity field is smoothed on small scales, so that neighboring marker points generally have similar velocities. We demonstrate significant improvement gained by this adjustment in several code tests relevant to the physics which moving-mesh codes are designed to capture.

We present a sample of 27 GRBs with detailed Swift light curves supplemented by late time Chandra observations. To answer the missing jet-break problem in general, we develop a numerical simulation based model which can be directly fit to the data using Monte Carlo methods. Our numerical model takes into account all the factors that can shape a jet-break: (i) lateral expansion (ii) edge effects and (iii) off-axis effects. Our results provide improved fits to the light curves and constraints on physical parameters. More importantly, our results suggest that off-axis effects are important and must be included in interpretations of GRB jet breaks.

We present the results of 2D, moving mesh, viscous hydrodynamical simulations of accretion onto merging supermassive black hole (SMBH) binaries. We include viscous heating, shock heating, and radiative cooling, and simulate the transition from the "pre-decoupling" epoch, where the inspiral timescale is longer than the viscous timescale, to the "post-decoupling" epoch, where the inspiral timescale is shorter than the viscous timescale. We find that there is no abrupt halt to the accretion at decoupling, but rather the accretion shows a slow decay, with significant accretion well after the expected decoupling. Moreover, we find that the luminosity in X-rays is significantly higher prior to the merger, as orbital energy from the SMBH binary is converted to heat via strong shocks inside the cavity, and radiated away. Following the merger, the cavity refills viscously and the accretion rate relaxes to the Shakura-Sunyaev value, while the X-ray luminosity drops as the shocks quickly dissipate.

Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk.

Afterglow jets are Rayleigh-Taylor unstable and therefore turbulent during the early part of their deceleration. There are also several processes which actively cool the jet. In this letter, we demonstrate that if cooling significantly increases the compressibility of the flow, the turbulence collides with the forward shock, destabilizing and corrugating it. In this case, the forward shock is turbulent enough to produce the magnetic fields responsible for synchrotron emission via small scale turbulent dynamo. We calculate light curves assuming the magnetic field is in energy equipartition with the turbulent kinetic energy and discover that dynamic magnetic fields are well-approximated by a constant magnetic-to- thermal energy ratio of 1%, though there is a sizeable delay in the time of peak flux as the magnetic field turns on only after the turbulence has activated. The reverse shock is found to be significantly more magnetized than the forward shock, with a magnetic-to-thermal energy ratio of order 10%. This work motivates future Rayleigh-Taylor calculations using more physical cooling models.

The kinetic energy of supernovae (SNe) accompanied by gamma-ray bursts (GRBs) tends to cluster near E52 erg, with 2.E52 erg an upper limit to which no compelling exceptions are found (assuming a certain degree of asphericity), and it is always significantly larger than the intrinsic energy of the GRB themselves (corrected for jet collimation). This energy is strikingly similar to the maximum rotational energy of a neutron star rotating with period 1 ms. It is therefore proposed that all GRBs associated with luminous SNe are produced by magnetars. GRBs that result from black hole formation (collapsars) may not produce luminous SNe. X-ray Flashes (XRFs), which are associated with less energetic SNe, are produced by neutron stars with weaker magnetic field or lower spin.

We present the results of a calculation of the thermal spectrum from a 2D, moving mesh, high-accuracy, viscous hydrodynamical simulation of an accreting supermassive black hole binary. We include viscous heating, shock heating, and radiative cooling, evolving for longer than a viscous time so that we reach a quasi-steady accretion state. In agreement with previous work, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We also find that emission from the shock-heated mini-disks and accretion streams prevents any deficit in high-energy emission that may be expected inside the circumbinary cavity, and instead leads to a characteristic brightening of the spectrum beginning in soft X-rays.

We use the ScaleFit package to perform Markov chain Monte Carlo light curve fitting on a large sample of Swift-XRT gamma-ray burst afterglows. The ScaleFit model uses scaling relations in the hydrodynamic and radiation equations to compute synthetic light curves directly from a set of high resolution two-dimensional relativistic blast wave simulations. The data sample consists of all Swift-XRT afterglows from 2005 to 2012 with sufficient coverage and a known redshift, 188 bursts in total. We find the jet half opening angle varies widely but is commonly less than 0.1 radians. The distribution of the electron spectral index is also broad, with a median at 2.23. This approach allows, for the first time, for the off-axis observer angle to be inferred directly from the light curve. We find the observer angle to have a median value of 0.57 of the jet opening angle over our sample, which has profound consequences for the predicted rate of observed jet breaks and effects the beaming corrected total energies of gamma-ray bursts.

The presence of fallback disks around young neutron stars has been invoked over the years to explain a large variety of phenomena. Here we perform a numerical investigation of the formation of such disks during a supernova explosion, considering both neutron star (NS) and black hole (BH) remnants. Using the public code MESA, we compute the angular momentum distribution of the pre-supernova material, for stars with initial masses M in the range 13 - 40 Msun, initial surface rotational velocities vsurf between 25% and 75% of the critical velocity, and for metallicities Z of 1%, 10% and 100% of the solar value. These pre SN models are exploded with energies E varying between 10^{50} - 3x10^{52} ergs, and the amount of fallback material is computed. We find that, if magnetic torques play an important role in angular momentum transport, then fallback disks around NSs, even for low-metallicity main sequence stars, are not an outcome of SN explosions. Formation of such disks around young NSs can only happen under the condition of negligible magnetic torques and a fine-tuned explosion energy. For those stars which leave behind BH remnants, disk formation is ubiquitous if magnetic fields do not play a strong role; however, unlike the NS case, even with strong magnetic coupling in the interior, a disk can form in a large region of the {Z,M,vsurf,E} parameter space. Together with the compact, hyperaccreting fallback disks widely discussed in the literature, we identify regions in the above parameter space which lead to extended, long-lived disks around BHs. We find that the physical conditions in these disks may be conducive to planet formation, hence leading to the possible existence of planets orbiting black holes.

A near-equal mass binary black hole can clear a central cavity in a circumbinary accretion disk; however, previous works have revealed accretion streams entering this cavity. Here we use 2D hydrodynamical simulations to study the accretion streams and their periodic behavior. In particular, we perform a suite of simulations, covering different binary mass ratios $q=M_2/M_1$ in the range $0.01 \leq q \leq 1$. In each case, we follow the system for several thousand binary orbits, until it relaxes to a stable accretion pattern. We find the following results: (i) while the binary is efficient in maintaining a low-density cavity, the time-averaged mass accretion rate into the cavity, through narrow coherent accretion streams, is suppressed by at most a factor of $\sim 5$, compared to a disk with a single BH with the same mass; (ii) the largest suppression occurs for $q\approx 0.05$; binaries whose mass ratios are either lower or higher both suppress accretion less significantly; (iii) for $q \gtrsim 0.05$, the accretion rate is strongly modulated by the binary, and depending on the precise value of $q$, the power spectrum of the accretion rate shows either one, two, or three distinct periods; and (v) for $q \lesssim 0.05$, the accretion rate becomes steady, with no time-variations. Most binaries produced in galactic mergers are expected to have $q\gtrsim 0.05$. If the luminosity of these binaries tracks their accretion rate, then a periodogram of their light-curve could help in their identification, and to constrain their mass ratio and disk properties.

We numerically calculate the growth and saturation of the Rayleigh-Taylor instability caused by the deceleration of relativistic outflows with Lorentz factor {\Gamma} = 10, 30, and 100. The instability generates turbulence whose scale exhibits strong dependence on Lorentz factor, as only modes with angular size smaller than 1/{\Gamma} can grow. We develop a simple diagnostic to measure the kinetic energy in turbulent fluctuations, and calculate a ratio of turbulent kinetic energy to thermal energy of .03 in the region affected by the instability. Although our numerical calculation does not include magnetic fields, we argue that small scale turbulent dynamo amplifies magnetic fields to nearly this same fraction, giving a ratio of magnetic to thermal energy of ~ .01, to within a factor of two. The instability completely disrupts the contact discontinuity between the ejecta and the swept up circumburst medium. The reverse shock is stable, but is impacted by the Rayleigh-Taylor instability, which strengthens the reverse shock and pushes it away from the forward shock. The forward shock front is unaffected by the instability, but Rayleigh-Taylor fingers can penetrate of order 10% of the way into the energetic region behind the shock during the two-shock phase of the explosion. We calculate afterglow emission from the explosion and find the reverse shock emission peaks at a later time due to its reduced Lorentz factor and modified density and pressure at the shock front. These calculations are performed using a novel numerical technique that includes a moving computational grid. The moving grid is essential as it maintains contact discontinuities to high precision and can easily evolve flows with extremely large Lorentz factors.

We present a model for relativistic jets which generates a particular angular distribution of Lorentz factor and energy per solid angle. We consider a fireball with specific internal energy E/M launched with bulk Lorentz factor \gamma_B. This "boosted fireball" model is motivated by the phenomenology of collapsar jets, but is applicable to a wide variety of relativistic flows. In its center-of-momentum frame the fireball expands isotropically, converting its internal energy into radially expanding flow with asymptotic Lorentz factor \eta_0 ~ E/M. In the lab frame the flow is beamed, expanding with Lorentz factor \Gamma = 2 \eta_0 \gamma_B in the direction of its initial bulk motion and with characteristic opening angle \theta_0 ~ 1/\gamma_B. The flow is jet-like with \Gamma \theta_0 ~ 2 \eta_0 such that jets with \Gamma > 1/\theta_0 are naturally produced. The choice \eta_0 ~ \gamma_B ~ 10 yields a jet with \Gamma ~ 200 on-axis and angular structure characterized by opening angle \theta_0 ~ 0.1 of relevance for cosmological GRBs, while \gamma_B >~ 1 may be relevant for low-luminosity GRBs. The model produces a family of outflows, of relevance for different relativistic phenomena with structures completely determined by \eta_0 and \gamma_B. We calculate the energy per unit solid angle for the model and use it to compute light curves for comparison with the widely used top-hat model. The jet break in the boosted fireball light curve is greatly subdued when compared to the top-hat model because the edge of the jet is smoother than for a top-hat. This may explain missing jet breaks in afterglow light curves.

We present the results of 2D hydrodynamical simulations of circumbinary disk accretion using the finite-volume code DISCO. This code solves the 2D viscous Navier-Stokes equations on a high-resolution moving mesh which shears with the fluid flow, greatly reducing advection errors in comparison with a fixed grid. We perform a series of simulations for binary mass ratios in the range 0.026 < q < 1.0, each lasting longer than a viscous time so that we reach a quasi-steady accretion state. In each case, we find that gas is efficiently stripped from the inner edge of the circumbinary disk and enters the cavity along accretion streams, which feed persistent "mini-disks" surrounding each black hole. We find that for q > 0.1, the binary excites eccentricity in the inner region of the circumbinary disk, creating an overdense lump which gives rise to enhanced periodicity in the accretion rate. The dependence of the periodicity on mass ratio may provide a method for observationally inferring mass ratios from measurements of the accretion rate. We also find that for all mass ratios studied, the magnitude of the accretion onto the secondary is sufficient to drive the binary toward larger mass ratio. This suggests a mechanism for biasing mass ratio distributions toward equal mass.

Afterglows of gamma-ray bursts are observed to produce light curves with the flux following power law evolution in time. However, recent observations reveal bright flares at times on the order of minutes to days. One proposed explanation for these flares is the interaction of a relativistic blast wave with a circumburst density transition. In this paper, we model this type of interaction computationally in one and two dimensions, using a relativistic hydrodynamics code with adaptive mesh refinement called ram, and analytically in one dimension. We simulate a blast wave traveling in a stellar wind environment that encounters a sudden change in density, followed by a homogeneous medium, and compute the observed radiation using a synchrotron model. We show that flares are not observable for an encounter with a sudden density increase, such as a wind termination shock, nor for an encounter with a sudden density decrease. Furthermore, by extending our analysis to two dimensions, we are able to resolve the spreading, collimation, and edge effects of the blast wave as it encounters the change in circumburst medium. In all cases considered in this paper, we find that a flare will not be observed for any of the density changes studied.

We analyze the equations of relativistic magnetized plasma dynamics in the limiting case that electromagnetic stress-energy is dominant over pressure and rest mass energy density. The naive formulation of these equations is shown to be not hyperbolic. Modifying the equations by terms that vanish for all physical solutions, we obtain a symmetric hyperbolic evolution system, which should exhibit improved numerical behavior.

We present preliminary data fit results of synthetic light curves computed from numerical afterglow blast wave simulations. Our technique uses Markov chain Monte Carlo (MCMC) in a new data analysis tool, ScaleFit. Scaling relations in both the hydrodynamics and radiation equations allow light curves to be parameterized by a small set of scale-invariant characteristic quantities. These quantities have been calculated and tabulated from high resolution two-dimensional hydrodynamic simulations. Producing a light curve from the characteristics takes only a millisecond, allowing for the use of MCMC data fitting techniques which can require millions of iterations. ScaleFit is a portable, lightweight, python package which performs this analysis on afterglow light curves. Using the set of Swift-XRT light curves from 2011 & 2012 with known redshifts, we find ScaleFit can measure the jet opening angle, observer angle, and spectral index of most afterglows. Globally we find gamma-ray burst afterglows tend to be observed off axis, at a significant fraction of the jet opening angle.

Afterglows of gamma-ray bursts are observed to produce light curves with the flux following power law evolution in time. However, recent observations reveal bright flares at times on the order of minutes to days. One proposed explanation for these flares is the interaction of a relativistic blast wave with a circumburst density transition. In this paper, we model this type of interaction computationally in one and two dimensions, using a relativistic hydrodynamics code with adaptive mesh refinement called RAM, and analytically in one dimension. We simulate a blast wave traveling in a stellar wind environment that encounters a sudden change in density, followed by a homogeneous medium, and compute the observed radiation using a synchrotron model. We show that flares are not observable for an encounter with a sudden density increase, such as a wind termination shock, nor for an encounter with a sudden density decrease. Furthermore, by extending our analysis to two dimensions, we are able to resolve the spreading, collimation, and edge effects of the blast wave as it encounters the change in circumburst medium. In all cases considered in this paper, we find that a flare will not be observed for any of the density changes studied.

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

By numerically integrating the compressible Navier-Stokes equations in two dimensions, we calculate the criterion for gap formation by a very low mass (q ~10^{-4}) protoplanet on a fixed orbit in a thin viscous disk. In contrast with some previously proposed gap-opening criteria, we find that a planet can open a gap even if the Hill radius is smaller than the disk scale height. Moreover, in the low-viscosity limit, we find no minimum mass necessary to open a gap for a planet held on a fixed orbit. In particular, a Neptune-mass planet will open a gap in a minimum mass solar nebula with suitably low viscosity (\alpha <10^{-4}). We find that the mass threshold scales as the square root of viscosity in the low mass regime. This is because the gap width for critical planet masses in this regime is a fixed multiple of the scale height, not of the Hill radius of the planet.

GRB afterglow jets have been notoriously difficult to resolve numerically using 2D hydrodynamical simulations due to high outflow Lorentz factors. By performing simulations in a boosted frame, it is possible to calculate light curves from numerically computed flows in sufficient detail to accurately quantify the shape of the jet break and the post-break steepening of the light curve. We study jet breaks for jets with opening angles of theta_0 = 0.05, 0.1 and 0.2 radians decelerating in a constant density medium, observed at angles (theta_obs) ranging from on-axis to the jet edge. We present a single set of scale-invariant functions describing the time evolution of synchrotron spectral break frequencies and peak flux, depending only on theta_0 and theta_obs, sufficient to reconstruct light curves for arbitrary explosion energy, circumburst density and synchrotron slope p. We compare our light curves to fit functions using sharp breaks and smooth transitions and confirm that the measured jet break time can be delayed due to theta_obs. We find that the difference in temporal indices across the jet break is larger than theoretically anticipated leading to post-break slopes of roughly 0.25 - 1.3 p, sufficiently steep to not easily reconcile with post-break slopes measured for the Swift sample, suggesting that most Swift GRBs do not explode in homogeneous media or that jet breaks are hidden by additional physics such as prolonged energy injection or viewing angle effects. Light curves and spectral break and peak flux evolution functions will be made publicly available on-line at http://cosmo.nyu.edu/afterglowlibrary .

The simultaneous detection of electromagnetic and gravitational wave emission from merging neutron star binaries would aid greatly in their discovery and interpretation. By studying turbulent amplification of magnetic fields in local high-resolution simulations of neutron star merger conditions, we demonstrate that magnetar-level (~10^16) G fields are present throughout the merger duration. We find that the small-scale turbulent dynamo converts 60% of the randomized kinetic energy into magnetic fields on a merger time scale. Since turbulent magnetic energy dissipates through reconnection events which accelerate relativistic electrons, turbulence may facilitate the conversion of orbital kinetic energy into radiation. If 10^-4 of the ~ 10^53 erg of orbital kinetic available gets processed through reconnection, and creates radiation in the 15-150 keV band, then the fluence at 200 Mpc would be 10^-7 erg/cm^2, potentially rendering most merging neutron stars in the advanced LIGO and Virgo detection volumes detectable by Swift BAT.

We present X-ray, UV/optical, and radio observations of the stripped-envelope, core-collapse supernova (SN) 2011ei, one of the least luminous SNe IIb or Ib observed to date. Our observations begin with a discovery within 1 day of explosion and span several months afterward. Early optical spectra exhibit broad, Type II-like hydrogen Balmer profiles that subside rapidly and are replaced by Type Ib-like He-rich features on the timescale of one week. High-cadence monitoring of this transition suggests that absorption attributable to a high velocity (> 12,000 km/s) H-rich shell is not rare in Type Ib events. Radio observations imply a shock velocity of v = 0.13c and a progenitor star mass-loss rate of 1.4 x 10^{-5} Msun yr^{-1} (assuming wind velocity v_w=10^3 km/s). This is consistent with independent constraints from deep X-ray observations with Swift-XRT and Chandra. Overall, the multi-wavelength properties of SN 2011ei are consistent with the explosion of a lower-mass (3-4 Msun), compact (R* <= 1x10^{11} cm), He core star. The star retained a thin hydrogen envelope at the time of explosion, and was embedded in an inhomogeneous circumstellar wind suggestive of modest episodic mass-loss. We conclude that SN 2011ei's rapid spectral metamorphosis is indicative of time-dependent classifications that bias estimates of explosion rates for Type IIb and Ib objects, and that important information about a progenitor star's evolutionary state and mass-loss immediately prior to SN explosion can be inferred from timely multi-wavelength observations.

High-resolution numerical simulations are utilized to examine isotropic turbulence in a compressible fluid when long wavelength velocity fluctuations approach light speed. Spectral analysis reveals an inertial sub-range of relativistic motions with a broadly 5/3 index. The use of generalized Lorentz-covariant structure functions based on the four-velocity is proposed. These structure functions extend the She-Leveque model for intermittency into the relativistic regime.

We calculate the global quasi-steady state of a thin disk perturbed by a low-mass protoplanet orbiting at a fixed radius using extremely high-resolution numerical integrations of Euler's equations in two dimensions. The calculations are carried out using a moving computational domain, which greatly reduces advection errors and allows for much longer time-steps than a fixed grid. We calculate the angular momentum flux and the torque density as a function of radius and compare them with analytical predictions. We discuss the quasi-steady state after 100 orbits and the prospects for gap formation by low mass planets.

We present a powerful new tool for fitting broadband gamma-ray burst afterglow data, which can be used to determine the burst explosion parameters and the synchrotron radiation parameters. By making use of scale invariance between relativistic jets of different energies and different circumburst medium densities, and by capturing the output of high-resolution two-dimensional relativistic hydrodynamical (RHD) jet simulations in a concise summary, the jet dynamics are generated quickly. Our method calculates the full light curves and spectra using linear radiative transfer sufficiently fast to allow for a direct iterative fit of RHD simulations to the data. The fit properly accounts for jet features that so far have not been successfully modeled analytically, such as jet decollimation, inhomogeneity along the shock front and the transitory phase between the early time relativistic and late time non-relativistic outflow. As a first application of the model we simultaneously fit the radio, X-ray and optical data of GRB 990510. We not only find noticeable differences between our findings for the explosion and radiation parameters and those of earlier authors, but also an improved model fit when we include the observer angle in the data fit. The fit method will be made freely available on request and on-line at http://cosmo.nyu.edu/afterglowlibrary. In addition to data fitting, the software tools can also be used to quickly generate a light curve or spectrum for arbitrary observer position, jet and radiation parameters.

We discuss jet dynamics for narrow and wide gamma-ray burst (GRB) afterglow jets and the observational implications of numerical simulations of relativistic jets in two dimensions. We confirm earlier numerical results that sideways expansion of relativistic jets during the bulk of the afterglow emission phase is logarithmic in time and find that this also applies to narrow jets with half opening angle of 0.05 radians. As a result, afterglow jets remain highly nonspherical until after they have become nonrelativistic. Although sideways expansion steepens the afterglow light curve after the jet break, the jet edges becoming visible dominates the jet break, which means that the jet break is sensitive to the observer angle even for narrow jets. Failure to take the observer angle into account can lead to an overestimation of the jet energy by up to a factor 4. This weakens the challenge posed to the magneter energy limit by extreme events such as GRB090926A. Late time radio calorimetry based on a spherical nonrelativistic outflow model remains relevant when the observer is approximately on-axis and where differences of a few in flux level between the model and the simulation are acceptable. However, this does not imply sphericity of the outflow and therefore does not translate to high observer angles relevant to orphan afterglows. For more accurate calorimetry and in order to model significant late time features such as the rise of the counterjet, detailed jet simulations remain indispensable.

We demonstrate that gamma-ray burst afterglow spectra and light curves can be calculated for arbitrary explosion and radiation parameters by scaling the peak flux and the critical frequencies connecting different spectral regimes. Only one baseline calculation needs to be done for each jet opening angle and observer angle. These calculations are done numerically using high-resolution relativistic hydrodynamical afterglow blast wave simulations which include the two-dimensional dynamical features of expanding and decelerating afterglow blast waves. Any light curve can then be generated by applying scaling relations to the baseline calculations. As a result, it is now possible to fully fit for the shape of the jet break, e.g. at early time X-ray and optical frequencies. In addition, late-time radio calorimetry can be improved since the general shape of the transition into the Sedov-Taylor regime is now known for arbitrary explosion parameters so the exact moment when the Sedov-Taylor asymptote is reached in the light curve is no longer relevant. When calculating the baselines, we find that the synchrotron critical frequency and the cooling break frequency are strongly affected by the jet break. The synchrotron break temporal slope quickly drops to the steep late time Sedov-Taylor slope, while the cooling break first steepens then rises to meet the level of its shallow late time asymptote.

We have generalized a method for the numerical solution of hyperbolic systems of equations using a dynamic Voronoi tessellation of the computational domain. The Voronoi tessellation is used to generate moving computational meshes for the solution of multi-dimensional systems of conservation laws in finite-volume form. The mesh generating points are free to move with arbitrary velocity, with the choice of zero velocity resulting in an Eulerian formulation. Moving the points at the local fluid velocity makes the formulation effectively Lagrangian. We have written the TESS code to solve the equations of compressible hydrodynamics and magnetohydrodynamics for both relativistic and non-relativistic fluids on a dynamic Voronoi mesh. When run in Lagrangian mode, TESS is significantly less diffusive than fixed mesh codes and thus preserves contact discontinuities to high precision while also accurately capturing strong shock waves. TESS is written for Cartesian, spherical and cylindrical coordinates and is modular so that auxilliary physics solvers are readily integrated into the TESS framework and so that the TESS framework can be readily adapted to solve general systems of equations. We present results from a series of test problems to demonstrate the performance of TESS and to highlight some of the advantages of the dynamic tessellation method for solving challenging problems in astrophysical fluid dynamics.

Models for GRB outflows invoke turbulence in relativistically hot magnetized fluids. In order to investigate these conditions we have performed high-resolution three-dimensional numerical simulations of relativistic magneto-hydrodynamical (RMHD) turbulence. We find that magnetic energy is amplified to several percent of the total energy density by turbulent twisting and folding of magnetic field lines. Values of epsilon_B near 1% are thus naturally expected. We study the dependence of saturated magnetic field energy fraction as a function of Mach number and relativistic temperature. We then present power spectra of the turbulent kinetic and magnetic energies. We also present solenoidal (curl-like) and dilatational (divergence-like) power spectra of kinetic energy. We propose that relativistic effects introduce novel couplings between these spectral components. The case we explore in most detail is for equal amounts of thermal and rest mass energy, corresponding to conditions after collisions of shells with relative Lorentz factors of several. These conditions are relevant in models for internal shocks, for the late afterglow phase, for cocoon material along the edge of a relativistic jet as it propagates through a star, as well neutron stars merging with each other and with black hole companions. We find that relativistic turbulence decays extremely quickly, on a sound crossing time of an eddy. Models invoking sustained relativistic turbulence to explain variability in GRB prompt emission are thus strongly disfavored unless a persistant driving of the turbulence is maintained for the duration of the prompt emission.

A wide variety of astrophysical phenomena involve the flow of turbulent magnetized gas with relativistic velocity or energy density. Examples include gamma-ray bursts, active galactic nuclei, pulsars, magnetars, micro-quasars, merging neutron stars, X-ray binaries, some supernovae, and the early universe. In order to elucidate the basic properties of the relativistic magnetohydrodynamical (RMHD) turbulence present in these systems, we present results from numerical simulations of fully developed driven turbulence in a relativistically warm, weakly magnetized and mildly compressible ideal fluid. We have evolved the RMHD equations for many dynamical times on a uniform grid with 1024^3 zones using a high order Godunov code. We observe the growth of magnetic energy from a seed field through saturation at about 1% of the total fluid energy. We compute the power spectrum of velocity and density-weighted velocity and conclude that the inertial scaling is consistent with a slope of -5/3. We compute the longitudinal and transverse velocity structure functions of order p up to 11, and discuss their possible deviation from the expected scaling for non-relativistic media. We also compute the scale-dependent distortion of coherent velocity structures with respect to the local magnetic field, finding a weaker scale dependence than is expected for incompressible non-relativistic flows with a strong mean field.

Self similarity allows for analytic or semi-analytic solutions to many hydrodynamics problems. Most of these solutions are one dimensional. Using linear perturbation theory, expanded around such a one-dimensional solution, we find self-similar hydrodynamic solutions that are two- or three-dimensional. Since the deviation from a one-dimensional solution is small, we call these slightly two-dimensional and slightly three-dimensional self-similar solutions, respectively. As an example, we treat strong spherical explosions of the second type. A strong explosion propagates into an ideal gas with negligible temperature and density profile of the form rho(r,theta,phi)=r^{-omega}[1+sigma*F(theta,phi)], where omega>3 and sigma << 1. Analytical solutions are obtained by expanding the arbitrary function F(theta,phi) in spherical harmonics. We compare our results with two dimensional numerical simulations, and find good agreement.

We present results for a large number of gamma-ray burst (GRB) afterglow light curve calculations, done by combining high resolution two-dimensional relativistic hydrodynamics simulations using RAM with a synchrotron radiation code. Results were obtained for jet energies, circumburst medium densities and jet angles typical for short and underluminous GRBs, different observer angles and observer frequencies from low radio (75 MHz) to X-ray (1.5 keV). We summarize the light curves through smooth power law fits with up to three breaks, covering jet breaks for small observer angles, the rising phase for large observer angles and the rise and decay of the counterjet. All light curve data are publicly available via http://cosmo.nyu.edu/afterglowlibrary . The data can be used for model fits to observational data and as an aid for predicting observations by future telescopes such as LOFAR or SKA and will benefit the study of neutron star mergers using different channels, such as gravitational wave observations with LIGO or Virgo. For small observer angles, we find jet break times that vary significantly between frequencies, with the break time in the radio substantially postponed. Increasing the observer angle also postpones the measured jet break time. The rising phase of the light curve for large observer angle has a complex shape that can not always be summarized by a simple power law. Except for very large observer angles, the counter jet is a distinct feature in the light curve, although in practice the signal will be exceedingly difficult to observe by then.

Numerical jet simulations serve a valuable role in calculating gamma-ray burst afterglow emission beyond analytical approximations. Here we present the results of high resolution 2D simulations of decelerating relativistic jets performed using the RAM adaptive mesh refinement relativistic hydrodynamics code. We have applied a separate synchrotron radiation code to the simulation results in order to calculate light curves at frequencies varying from radio to X-ray for observers at various angles from the jet axis. We provide a confirmation from radio light curves from simulations rather than from a simplified jet model for earlier results in the literature finding that only a very small number of local Ibc supernovae can possibly harbor an orphan afterglow. Also, recent studies have noted an unexpected lack of observed jet breaks in the Swift sample. Using a jet simulation with physical parameters representative for an average Swift sample burst, such as a jet half opening angle of 0.1 rad and a source redshift of z = 2.23, we have created synthetic light curves at 1.5 keV with artificial errors while accounting for Swift instrument biases as well. A large set of these light curves have been generated and analyzed using a Monte Carlo approach. Single and broken power law fits are compared. We find that for increasing observer angle, the jet break quickly becomes hard to detect. This holds true even when the observer remains well within the jet opening angle. We find that the odds that a Swift light curve from a randomly oriented 0.1 radians jet at z = 2.23 will exhibit a jet break at the 3 sigma level are only 12 percent. The observer angle therefore provides a natural explanation for the lack of perceived jet breaks in the Swift sample.

Numerical studies of gamma-ray burst afterglow jets reveal significant qualitative differences with simplified analytical models. We present an on-line library of synthetic afterglow light curves and broadband spectra for use in interpreting observational data. Light curves have been calculated for various physics settings such as explosion energy and circumburst structure, as well as differing jet parameters and observer angle and redshift. Calculations gave been done for observer frequencies ranging from low radio to X-ray and for observer times from hours to decades after the burst. The light curves have been calculated from high-resolution 2D hydrodynamical simulations performed with the RAM adaptive-mesh refinement code and a detailed synchrotron radiation code. The library will contain both generic afterglow simulations as well as specific case studies and will be freely accessible at http://cosmo.nyu.edu/afterglowlibrary . The synthetic light curves can be used as a check on the accuracy of physical parameters derived from analytical model fits to afterglow data, to quantitatively explore the consequences of varying parameters such as observer angle and for accurate predictions of future telescope data.

Starting as highly relativistic collimated jets, gamma-ray burst outflows gradually decelerate and become non-relativistic spherical blast waves. Although detailed analytical solutions describing the afterglow emission received by an on-axis observer during both the early and late phases of the outflow evolution exist, a calculation of the received flux during the intermediate phase and for an off-axis observer requires either a more simplified analytical model or direct numerical simulations of the outflow dynamics. In this paper we present light curves for off-axis observers covering the long-term evolution of the blast wave calculated from a high resolution two-dimensional relativistic hydrodynamics simulation using a synchrotron radiation model. We compare our results to earlier analytical work and calculate the consequence of the observer angle with respect to the jet axis both for the detection of orphan afterglows and for jet break fits to the observational data. We find that observable jet breaks can be delayed for up to several weeks for off-axis observers, potentially leading to overestimation of the beaming corrected total energy. When using our off-axis light curves to create synthetic Swift X-ray data, we find that jet breaks are likely to remain hidden in the data. We also confirm earlier results in the literature finding that only a very small number of local Type Ibc supernovae can harbor an orphan afterglow.

While there is mounting evidence that long Gamma-Ray Bursts (GRBs) are associated with the collapse of massive stars, the detailed structure of their pre-supernova stage is still debatable. Particularly uncertain is the degree of mixing among shells of different composition, and hence the role of magnetic torques and convection in transporting angular momentum. Here we show that early-time afterglow observations with the Swift satellite place constraints on the allowed GRB pre-supernova models. In particular, they argue against pre-supernova models in which different elemental shells are unmixed. These types of models would produce energy injections into the GRB engine on timescales between several hundreds of seconds to a few hours. Flaring activity has {\em not} been observed in a large fraction of well-monitored long GRBs. Therefore, if the progenitors of long GRBs have common properties, then the lack of flares indicates that the massive stars which produce GRBs are mostly well mixed, as expected in low-metallicity, rapidly rotating massive stars.

Finding electromagnetic (EM) counterparts of future gravitational wave (GW) sources would bring rich scientific benefits. A promising possibility, in the case of the coalescence of a super-massive black hole binary (SMBHB), is that prompt emission from merger-induced disturbances in a supersonic circumbinary disk may be detectable. We follow the post-merger evolution of a thin, zero-viscosity circumbinary gas disk with two-dimensional simulations, using the hydrodynamic code FLASH. We analyze perturbations arising from the 530 km/s recoil of a 10^6 M_sun binary, oriented in the plane of the disk, assuming either an adiabatic or a pseudo-isothermal equation of state for the gas. We find that a single-armed spiral shock wave forms and propagates outward, sweeping up about 20% of the mass of the disk. The morphology and evolution of the perturbations agrees well with those of caustics predicted to occur in a collisionless disk. Assuming that the disk radiates nearly instantaneously to maintain a constant temperature, we estimate the amount of dissipation and corresponding post-merger light-curve. The luminosity rises steadily on the time-scale of months, and reaches few times 10^{43} erg/s, corresponding to about 10% of the Eddington luminosity of the central SMBHB. We also analyze the case in which gravitational wave emission results in a 5% mass loss in the merger remnant. The mass-loss reduces the shock overdensities and the overall luminosity of the disk by 15-20%, without any other major effects on the spiral shock pattern.

Direct multi-dimensional numerical simulation is the most reliable approach for calculating the fluid dynamics and observational signatures of relativistic jets in gamma-ray bursts (GRBs). We present a two-dimensional relativistic hydrodynamic simulation of a GRB outflow during the afterglow phase, which uses the fifth-order weighted essentially non-oscillatory scheme and adaptive mesh refinement. Initially, the jet has a Lorentz factor of 20. We have followed its evolution up to 150 years. Using the hydrodynamic data, we calculate synchrotron radiation based upon standard afterglow models and compare our results with previous analytic work. We find that the sideways expansion of a relativistic GRB jet is a very slow process and previous analytic works have overestimated its rate. In our computed lightcurves, a very sharp jet break is seen and the post-break lightcurves are steeper than analytic predictions. We find that the jet break in GRB afterglow lightcurves is mainly caused by the missing flux when the edge of the jet is observed. The outflow becomes nonrelativistic at the end of the Blandford-McKee phase. But it is still highly nonspherical, and it takes a rather long time for it to become a spherical Sedov-von Neumann-Taylor blast wave. We find that the late-time afterglows become increasingly flatter over time. But we disagree with the common notion that there is a sudden flattening in lightcurves due to the transition into the Sedov-von Neumann-Taylor solution. We have also found that there is a bump in lightcurves at very late times ($\sim 1000$ days) due to radiation from the counter jet. We speculate that such a counter jet bump might have already been observed in GRB 980703.

Magnetic field strengths inferred for relativistic outflows including gamma-ray bursts (GRB) and active galactic nuclei (AGN) are larger than naively expected by orders of magnitude. We present three-dimensional relativistic magnetohydrodynamics (MHD) simulations demonstrating amplification and saturation of magnetic field by a macroscopic turbulent dynamo triggered by the Kelvin-Helmholtz shear instability. We find rapid growth of electromagnetic energy due to the stretching and folding of field lines in the turbulent velocity field resulting from non-linear development of the instability. Using conditions relevant for GRB internal shocks and late phases of GRB afterglow, we obtain amplification of the electromagnetic energy fraction to $\epsilon_B \sim 5 \times 10^{-3}$. This value decays slowly after the shear is dissipated and appears to be largely independent of the initial field strength. The conditions required for operation of the dynamo are the presence of velocity shear and some seed magnetization both of which are expected to be commonplace. We also find that the turbulent kinetic energy spectrum for the case studied obeys Kolmogorov's 5/3 law and that the electromagnetic energy spectrum is essentially flat with the bulk of the electromagnetic energy at small scales.

This paper reports on the results of a numerical investigation designed to address how the initially anisotropic appearance of a GRB remnant is modified by the character of the circumburst medium and by the possible presence of an accompanying supernova (SN). Axisymmetric hydrodynamical calculations of light, impulsive jets propagating in both uniform and inhomogeneous external media are presented, which show that the resulting dynamics of their remnants since the onset of the non-relativistic phase is different from the standard self-similar solutions. Because massive star progenitors are expected to have their close-in surroundings modified by the progenitor winds, we consider both free winds and shocked winds as possible external media for GRB remnant evolution. Abundant confirmation is provided here of the important notion that the morphology and visibility of GRB remnants are determined largely by their circumstellar environments. For this reason, their detectability is highly biased in favor of those with massive star progenitors; although, in this class of models, the beamed component may be difficult to identify because the GRB ejecta is eventually swept up by the accompanying SN. The number density of asymmetric GRB remnants in the local Universe could be, however, far larger if they expand in a tenuous interstellar medium, as expected for some short GRB progenitor models. In these sources, the late size of the observable, asymmetric remnant could extend over a wide, possibly resolvable angle and may be easier to constrain directly.

We present a two-dimensional grid-based hydrodynamic simulation of a thin, viscous, locally-isothermal corotating disk orbiting an equal-mass Newtonian binary point mass on a fixed circular orbit. We study the structure of the disk after multiple viscous times. The binary maintains a central hole in the viscously-relaxed disk with radius equal to about twice the binary semimajor axis. Disk surface density within the hole is reduced by orders of magnitude relative to the density in the disk bulk. The inner truncation of the disk resembles the clearing of a gap in a protoplanetary disk. An initially circular disk becomes elliptical and then eccentric. Disturbances in the disk contain a component that is stationary in the rotating frame in which the binary is at rest; this component is a two-armed spiral density wave. We measure the distribution of the binary torque in the disk and find that the strongest positive torque is exerted inside the central low-density hole. We make connection with the linear theory of disk forcing at outer Lindblad resonances (OLRs) and find that the measured torque density distribution is consistent with forcing at the 3:2 (m=2) OLR, well within the central hole. We also measure the time dependence of the rate at which gas accretes across the hole and find quasi-periodic structure. We discuss implications for variability and detection of active galactic nuclei containing a binary massive black hole.

We derive constraints that must be satisfied by the sources of ~10^{15} to ~10^{18} eV cosmic rays, under the assumption that the sources are Galactic. We show that while these constraints are not satisfied by ordinary supernovae (SNe), which are believed to be the sources of <10^{15} eV cosmic rays, they may be satisfied by the recently discovered class of trans-relativistic supernovae (TRSNe), which were observed in association with gamma-ray bursts. We define TRSNe as SNe that deposit a large fraction, f_R>10^{-2}, of their kinetic energy in mildly relativistic, \gamma\beta>1, ejecta. The high velocity ejecta enable particle acceleration to ~10^{18} eV, and the large value of f_R (compared to f_R~10^{-7} for ordinary SNe) ensures that if TRSNe produce the observed ~10^{18} eV cosmic ray flux, they do not overproduce the flux at lower energies. This, combined with the estimated rate and energy production of TRSNe, imply that Galactic TRSNe may be the sources of cosmic rays with energies up to ~10^{18}eV .

It has been suggested that magnetic fields play a dynamically-important role in core-collapse explosions of massive stars. In particular, they may be important in the collapsar scenario for gamma-ray bursts (GRB), where the central engine is a hyper-accreting black hole or a millisecond magnetar. The present paper is focussed on the magnetar scenario, with a specific emphasis on the interaction of the magnetar magnetosphere with the infalling stellar envelope. First, the ``Pulsar-in-a-Cavity'' problem is introduced as a paradigm for a magnetar inside a collapsing star. The basic set-up of this fundamental plasma-physics problem is described, outlining its main features, and simple estimates are derived for the evolution of the magnetic field. In the context of a collapsing star, it is proposed that, at first, the ram pressure of the infalling plasma acts to confine the magnetosphere, enabling a gradual build-up of the magnetic pressure. At some point, the growing magnetic pressure overtakes the (decreasing) ram pressure of the gas, resulting in a magnetically-driven explosion. The explosion should be highly anisotropic, as the hoop-stress of the toroidal field, confined by the surrounding stellar matter, collimates the magnetically-dominated outflow into two beamed magnetic-tower jets. This creates a clean narrow channel for the escape of energy from the central engine through the star, as required for GRBs. In addition, the delayed onset of the collimated-explosion phase can explain the production of large quantities of Nickel-56, as suggested by the GRB-Supernova connection. Finally, the prospects for numerical simulations of this scenario are discussed.

We consider a newly-born millisecond magnetar, focusing on its interaction with the dense stellar plasma in which it is initially embedded. We argue that the confining pressure and inertia of the surrounding plasma acts to collimate the magnetar's Poynting-flux-dominated outflow into tightly beamed jets and increases its magnetic luminosity. We propose this process as an essential ingredient in the magnetar model for gamma-ray burst and asymmetric supernova central engines. We introduce the ``pulsar-in-a-cavity'' as an important model problem representing a magnetized rotating neutron star inside a collapsing star. We describe its essential properties and derive simple estimates for the evolution of the magnetic field and the resulting spin-down power. We find that the infalling stellar mantle confines the magnetosphere, enabling a gradual build-up of the toroidal magnetic field due to continuous twisting. The growing magnetic pressure eventually becomes dominant, resulting in a magnetically-driven explosion. The initial phase of the explosion is quasi-isotropic, potentially exposing a sufficient amount of material to $^{56}$Ni-producing temperatures to result in a bright supernova. However, if significant expansion of the star occurs prior to the explosion, then very little $^{56}$Ni is produced and no supernova is expected. In either case, hoop stress subsequently collimates the magnetically-dominated outflow, leading to the formation of a magnetic tower. After the star explodes, the decrease in bounding pressure causes the magnetic outflow to become less beamed. However, episodes of late fallback can reform the beamed outflow, which may be responsible for late X-ray flares.

Geometrical shock dynamics, also called CCW theory, yields approximate equations for shock propagation in which only the conditions at the shock appear explicitly; the post-shock flow is presumed approximately uniform and enters implicitly via a Riemann invariant. The nonrelativistic theory, formulated by G. B. Whitham and others, matches many experimental results surprisingly well. Motivated by astrophysical applications, we adapt the theory to ultra-relativistic shocks advancing into an ideal fluid whose pressure is negligible ahead of the shock, but one third of its proper energy density behind the shock. Exact results are recovered for some self-similar cylindrical and spherical shocks with power-law pre-shock density profiles. Comparison is made with numerical solutions of the full hydrodynamic equations. We review relativistic vorticity and circulation. In an ultrarelativistic ideal fluid, circulation can be defined so that it changes only at shocks, notwithstanding entropy gradients in smooth parts of the flow.

We present detailed optical, near-infrared, and radio observations of the X-ray flash 050416a obtained with Palomar and Siding Springs Observatories as well as HST and the VLA, placing this event among the best-studied X-ray flashes to date. In addition, we present an optical spectrum from Keck LRIS from which we measure the redshift of the burst, z=0.6528. At this redshift the isotropic-equivalent prompt energy release was about 10^51 erg, and using a standard afterglow synchrotron model we find that the blastwave kinetic energy is a factor of 10 larger, E_K,iso ~ 10^52 erg. The lack of an observed jet break to t ~ 20 days indicates that the opening angle is larger than 7 deg and the total beaming-corrected relativistic energy is larger than 10^50 erg. We further show that the burst produced a strong radio flare at t ~ 40 days accompanied by an observed flattening in the X-ray band which we attribute to an abrupt circumburst density jump or an episode of energy injection (either from a refreshed shock or off-axis ejecta). Late-time observations with HST show evidence for an associated supernova with peak optical luminosity roughly comparable to that of SN 1998bw. Next, we show that the host galaxy of XRF 050416a is actively forming stars at a rate of at least 2 M_solar per year with a luminosity of L_B ~ 0.5L* and metallicity of Z ~ 0.2-0.8 Z_solar. Finally, we discuss the nature of XRF 050416a in the context of short-hard gamma-ray bursts and under the framework of off-axis and dirty fireball models for X-ray flashes.

We propose a magnetic mechanism for the collimated explosion of a massive star relevant for GRBs, XRFs and asymmetric supernovae. We apply Lynden-Bell's magnetic tower scenario to the interior of a massive rotating star after the core has collapsed to form a black hole with an accretion disk or a millisecond magnetar acting as a central engine. We solve the force-free Grad-Shafranov equation to calculate the magnetic structure and growth of a tower embedded in a stellar environment. The pressure of the toroidal magnetic field, continuously generated by differential rotation of the central engine, drives a rapid expansion which becomes vertically collimated after lateral force balance with the surrounding gas pressure is reached. The collimation naturally occurs because hoop stress concentrates magnetic field toward the rotation axis and inhibits lateral expansion. This leads to the growth of a self-collimated magnetic tower. When embedded in a massive star, the supersonic expansion of the tower drives a strong bow shock behind which an over-pressured cocoon forms. The cocoon confines the tower by supplying collimating pressure and provides stabilization against disruption due to MHD instabilities. Because the tower consists of closed field lines starting and ending on the central engine, mixing of baryons from the cocoon into the tower is suppressed. The channel cleared by the growing tower is thus plausibly free of baryons and allows the escape of magnetic energy from the central engine through the star. While propagating down the stellar density gradient, the tower accelerates and becomes relativistic. During the expansion, fast collisionless reconnection becomes possible resulting in dissipation of magnetic energy which may be responsible for GRB prompt emission.

We have developed a new computer code, RAM, to solve the conservative equations of special relativistic hydrodynamics (SRHD) using adaptive mesh refinement (AMR) on parallel computers. We have implemented a characteristic-wise, finite difference, weighted essentially non-oscillatory (WENO) scheme using the full characteristic decomposition of the SRHD equations to achieve fifth-order accuracy in space. For time integration we use the method of lines with a third-order total variation diminishing (TVD) Runge-Kutta scheme. We have also implemented fourth and fifth order Runge-Kutta time integration schemes for comparison. The implementation of AMR and parallelization is based on the FLASH code. RAM is modular and includes the capability to easily swap hydrodynamics solvers, reconstruction methods and physics modules. In addition to WENO we have implemented a finite volume module with the piecewise parabolic method (PPM) for reconstruction and the modified Marquina approximate Riemann solver to work with TVD Runge-Kutta time integration. We examine the difficulty of accurately simulating shear flows in numerical relativistic hydrodynamics codes. We show that under-resolved simulations of simple test problems with transverse velocity components produce incorrect results and demonstrate the ability of RAM to correctly solve these problems. RAM has been tested in one, two and three dimensions and in Cartesian, cylindrical and spherical coordinates. We have demonstrated fifth-order accuracy for WENO in one and two dimensions and performed detailed comparison with other schemes for which we show significantly lower convergence rates. Extensive testing is presented demonstrating the ability of RAM to address challenging open questions in relativistic astrophysics.

We present the results from a {\it Hubble Space Telescope/ACS} study of the supernovae associated with gamma-ray bursts 040924 ($z=0.86$) and 041006 ($z=0.71$). We find evidence that both GRBs were associated with a SN 1998bw-like supernova dimmed by $\sim 1.5$ and $\sim 0.3$ magnitudes, respectively, making GRB 040924 the faintest GRB-associated SN ever detected. We study the luminosity dispersion in GRB/XRF-associated SNe and compare to local Type Ibc supernovae from the literature. We find significant overlap between the two samples, suggesting that GRB/XRF-associated SNe are not necessarily more luminous nor produce more $^{56}$Ni than local SNe. Based on the current (limited) datasets, we find that the two samples may share a similar $^{56}$Ni production mechanism.

The discovery of long-lasting (~100 s) X-ray flares following short gamma-ray bursts initially called into question whether they were truly classical short-hard bursts. Opinion over the last few years has coalesced around the view that the short-hard bursts arise from the merger of pairs of neutron stars, or a neutron star merging with a stellar-mass black hole. The natural timescales associated with these processes, however, essentially preclude an X-ray flare lasting ~100 s. Here we show that an interaction between the GRB outflow and a non-compact stellar companion at a distance of ~a light-minute provides a natural explanation for the flares. In the model, the burst is triggered by the collapse of a neutron star after accreting matter from the companion. This is reminiscent of type Ia supernovae, where there is a wide distribution of delay times between formation and explosion, leading to an association with both star-forming galaxies and old ellipticals.

The final chapter in the long-standing mystery of the gamma-ray bursts (GRBs) centres on the origin of the short-hard class, suspected on theoretical grounds to result from the coalescence of neutron star or black hole binary systems. Numerous searches for the afterglows of short-hard bursts have been made, galvanized by the revolution in our understanding of long-duration GRBs that followed the discovery in 1997 of their broadband (X-ray, optical, and radio) afterglow emission. Here we present the discovery of the X-ray afterglow of a short-hard burst whose accurate position allows us to unambiguously associate it with a star-forming galaxy at redshift z=0.160, and whose optical lightcurve definitively excludes a supernova association. Together with results from three other recent short-hard bursts, this suggests that short-hard bursts release much less energy than the long-duration GRBs. Models requiring young stellar populations, such as magnetars and collapsars, are ruled out, while coalescing degenerate binaries remain the most promising progenitor candidates.

We present the results from an Hubble Space Telescope/ACS search for supernovae associated with X-ray flashes 020903, 040701, 040812 and 040916. We find strong evidence that XRF 020903 (z=0.25) was associated with a SN 1998bw-like supernova and confirm this using optical spectroscopy at t ~ 25 days. We find no evidence, however, for SN 1998bw-like supernovae associated with the other three events. In the case of XRF 040701 (z=0.21), we rule out even a faint supernova similar to SN 2002ap, using template light-curves for several local Type Ic supernovae. For the two cases in which the redshift is not known, XRFs 040812 and 040916, we derive robust redshift limits assuming they were accompanied by supernovae similar to SN 1998bw and compare these limits with photometric redshift constraints provided by their host galaxies. We supplement this analysis with results for three additional events (XRFs 011030, 020427 and 030723) and discuss the observed diversity of supernovae associated with X-ray flashes and gamma-ray bursts. We conclude that XRF-SNe exist, but can be significantly fainter than SN 1998bw, possibly consistent with the observed spread in local Type Ibc supernovae.

We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.

We present results from axisymmetric, time-dependent magnetohydrodynamic (MHD) simulations of the collapsar model for gamma-ray bursts. We begin the simulations after the 1.7 MSUN iron core of a 25 MSUN presupernova star has collapsed and study the ensuing accretion of the 7 MSUN helium envelope onto the central black hole formed by the collapsed iron core. We consider a spherically symmetric progenitor model, but with spherical symmetry broken by the introduction of a small, latitude-dependent angular momentum and a weak radial magnetic field. Our MHD simulations include a realistic equation of state, neutrino cooling, photodisintegration of helium, and resistive heating. Our main conclusion is that, within the collapsar model, MHD effects alone are able to launch, accelerate and sustain a strong polar outflow. We also find that the outflow is Poynting flux-dominated, and note that this provides favorable initial conditions for the subsequent production of a baryon-poor fireball.

Winds blown from accretion disks formed inside massive rotating stars may result in stellar explosions observable as Type Ibc and Type II supernovae. A key feature of the winds is their ability to produce the radioactive Nickel-56 necessary to power a supernova light curve. The wind strength depends on accretion disk cooling by neutrino emission and photo-disintegration of bound nuclei. These cooling processes depend on the angular momentum of the stellar progenitor via the virial temperature at the Kepler radius where the disk forms. The production of an observable supernova counterpart to a Gamma-Ray Burst (GRB) may therefore depend on the angular momentum of the stellar progenitor. Stars with low angular momentum may produce a GRB without making an observable supernova. Stars with large angular momentum may make extremely bright and energetic supernovae like SN 1998bw. Stars with an intermediate range of angular momentum may simultaneously produce a supernova and a GRB.

We examine the propagation of 2-dimensional relativistic jets through the stellar progenitor in the collapsar model for gamma-ray bursts. In agreement with previous studies, we find that relativistic jets are collimated by their passage through the stellar mantle. Interaction of these jets with the star and their own cocoons also causes mixing that sporadically decelerates the flow. We speculate that this mixing instability is chiefly responsible for the variable Lorentz factor needed in the internal shock model and for the complex light curves seen in many GRBs. In all cases studied, the jet is shocked deep inside the star following a brief period of adiabatic expansion. The jet that finally emerges from the star thus has a moderate Lorentz factor, modulated by mixing, and a very large internal energy. In a second series of calculations, we follow the escape of that sort of jet. Because of the large ratio of internal to kinetic energy in both the jet and its cocoon, the opening angle of the final jet is significantly greater than at breakout. A small amount of material emerges at large angles, but with a Lorentz factor still sufficiently large to make a weak GRB. This leads us to propose a "unified model" in which a variety of high energy transients, ranging from x-ray flashes to "classic" GRBs, may be seen depending upon the angle at which a standard collapsar is observed. We also speculate that the breakout of a relativistic jet and its collision with the stellar wind will produce a brief transient with properties similar to the class of "short-hard" GRBs. Implications of our calculations for GRB light curves, the luminosity-variability relation, and the GRB-supernova association are also discussed. (Abridged)

Within the framework of the internal shock model, we study the luminosity and the variability in gamma-ray bursts from collimated fireballs. In particular we pay attention to the role of the photosphere due to $e^\pm$ pairs produced by internal shock synchrotron photons. It is shown that the observed Cepheid-like relationship between the luminosity and the variability can be interpreted as a correlation between the opening angle of the fireball jet and the mass included at the explosion with a standard energy output. We also show that such a correlation can be a natural consequence of the collapsar model. Using a multiple-shell model, we numerically calculate the temporal profiles of gamma-ray bursts. Collimated jets, in which the typical Lorentz factors are higher than in wide jets, can produce more variable temporal profiles due to smaller angular spreading time scales at the photosphere radius. Our simulations quantitatively reproduce the observed correlation.

Recent observations suggest that long-duration gamma-ray bursts and their afterglows are produced by highly relativistic jets emitted in core-collapse explosions. As the jet makes its way out of the stellar mantle, a bow shock runs ahead and a strong thermal precursor is produced as the shock breaks out. Such erupting fireballs produce a very bright gamma-ray precursor as they interact with the thermal break-out emission. The prompt gamma-ray emission propagates ahead of the fireball before it becomes optically thin, leading to electron-positron pair loading and radiative acceleration of the external medium. The detection of such precursors would offer the possibility of diagnosing not only the radius of the stellar progenitor and the initial Lorentz factor of the collimated fireball, but also the density of the external environment.

A variety of stellar explosions powered by black hole accretion are discussed. All involve the failure of neutrino energy deposition to launch a strong supernova explosion. A key quantity which determines the type of high energy transient produced is the ratio of the engine operation time, $\rm t_{engine}$, to the time for the explosion to break out of the stellar surface, $\rm t_{bo}$. Stars with sufficient angular momentum produce collapsars -- black holes accreting rapidly through a disk -- in their centers. Collapsars can occur in stars with a wide range of radii depending on the amount of pre-collapse mass loss. The stellar radius and jet properties determine the degree of beaming of the explosion. In some cases the stellar envelope serves to focus the explosion to narrow beaming angles. The baryon loading of various models for classical GRBs formed in massive stars is examined and the consequences are explored. For $\rm t_{engine} > t_{bo}$, highly relativistic outflow is possible and classical GRBs accompanied by supernovae can be produced. In other cases hyper-energetic, asymmetric supernovae are produced. Longer GRBs ($t \gtaprx 100 \s$) can be produced by fallback following a weak neutrino-driven supernova explosion.

Using a collapsar progenitor model of MacFadyen & Woosley we have simulated the propagation of an axisymmetric jet through a collapsing rotating massive star with the GENESIS multi-dimensional relativistic hydrodynamic code. The jet forms as a consequence of an assumed (constant or variable) energy deposition in the range $10^{50}$ erg s$^{-1}$ to $10^{51}$ erg s$^{-1}$ within a $30^{\circ}$ cone around the rotation axis. The jet flow is strongly beamed ($\la$ few degrees), spatially inhomogeneous, and time dependent. The jet reaches the surface of the stellar progenitor ($R_{\ast} = 2.98 \times 10^{10} $cm) intact. At breakout the maximum Lorentz factor of the jet flow is 33. After breakout the jet accelerates into the circumstellar medium, whose density is assumed to decrease exponentially and then being constant $\rho_{\rm ext} = 10^{-5}$ gcm$^{-3}$. Outside the star the flow begins to expand also laterally ($v \sim c$), but the beam remains very well collimated. At a distance of $2.54 R_{\ast}$, where the simulation ends, the Lorentz factor has increased to 44.

We have studied the relativistic beamed outflow proposed to occur in the collapsar model of gamma-ray bursts. A jet forms as a consequence of an assumed energy deposition of $\sim 10^{50}- 10^{51}$ erg/s within a $30^{\circ}$ cone around the rotation axis of the progenitor star. The generated jet flow is strongly beamed ($\la$ few degrees) and reaches the surface of the stellar progenitor (r $\approx 3 10^{10} $cm) intact. At break-out the maximum Lorentz factor of the jet flow is about 33. Simulations have been performed with the GENESIS multi-dimensional relativistic hydrodynamic code.

We consider the explosion of supernovae and the possible production of a variety of high energy transients by delayed black hole formation in massive stars endowed with rotation. Following the launch of a ``successful'' shock by the usual neutrino powered mechanism, the inner layers of the star move outwards, but lack adequate energy to eject all the matter exterior to the neutron star. Over a period of minutes to hours a variable amount of mass, about 0.1 to 5 solar masses falls back into the collapsed remnant, often turning it into a black hole and establishing an accretion disk. The accretion rate, about 0.001 to 0.01 solar masses per second, is inadequate to produce a jet mediated by neutrino annihilation, but similar to that invoked in magnetohydrodynamic (MHD) models for gamma-ray bursts (GRBs). We thus consider the effect of jets formed by ``fallback'' in stars that are already in the process of exploding. We justify a parameterization of the jet power as a constant times the mass accretion rate, $\epsilon \dot {\rm M} c^2$, and explore the consequences of $\epsilon$ = 0.001 and 0.01. In supergiants, shock breakout produces bright x-ray transients that might be a diagnostic of the model. Jets produced by fallback should be more frequent than those made by the prompt formation of a black hole and may power the most common form of gamma-ray transient in the universe, although not the most common form seen so far by BATSE. Those are still attributed to prompt black hole formation, but it may be that the diverse energies observed for GRBs so far reflect chiefly the variable collimation of the jet inside the star and a consequently highly variable fraction of relativistic ejecta. Indeed, these events may all have a common total energy near 10$^{52}$ erg.

A diverse range of phenomena is possible when a black hole experiences very rapid accretion from a disk due to the incomplete explosion of a massive presupernova star endowed with rotation. In the most extreme case, the outgoing shock fails promptly in a rotating helium star, a black hole and an accretion disk form, and a strong gamma-ray burst (GRB) results. However, there may also be more frequently realized cases where the black hole forms after a delay of from several tens of seconds to several hours as approximately 0.1 to 5 solar masses falls back into the collapsed remnant following a mildly successful supernova explosion. There, the same MHD mechanisms frequently invoked to produce GRBs would also produce jets in stars already in the process of exploding. The presupernova star could be a Wolf-Rayet star or a red or blue supergiant. Depending upon its initial pressure, the collimation of the jet may also vary since ``hot'' jets will tend to diverge and share their energy with the rest of the star. From these situations, one expects diverse outcomes ranging from GRBs with a large range of energies and durations, to asymmetric, energetic supernovae with weak GRBs. SN 1998bw may have been the explosion of a star in which fall back produced a black hole and a less collimated jet than in the case of prompt black hole formation.

Using a two-dimensional hydrodynamics code (PROMETHEUS), we study the continued evolution of rotating massive helium stars whose iron core collapse does not produce a successful outgoing shock, but instead forms a black hole. We study the formation of a disk, the associated flow patterns, and the accretion rate for disk viscosity parameter, alpha ~ 0.001 and 0.1. For the standard 14 solar mass model the average accretion rate for 15 s is 0.07 solar masses per second and the total energy deposited along the rotational axes by neutrino annihilation is (1 - 14) x 10**51 erg, depending upon the evolution of the Kerr parameter and uncertain neutrino efficiencies. Simulated deposition of this energy in the polar regions results in strong relativistic outflow - jets beamed to about 1.5% of the sky. The jets remain highly focused, and are capable of penetrating the star in 5 - 10 s. After the jet breaks through the surface of the star, highly relativistic flow can commence. Because of the sensitivity of the mass ejection and jets to accretion rate, angular momentum, and disk viscosity, and the variation of observational consequences with viewing angle, a large range of outcomes is possible ranging from bright GRBs like GRB 971214 to faint GRB-supernovae like SN 1998bw. X-ray precursors are also possible as the jet first breaks out of the star. While only a small fraction of supernovae make GRBs, we predict that all GRBs longer than a few seconds will make supernovae similar to SN 1998bw. However, hard, energetic GRBs shorter than a few seconds will be difficult to make in this model.