The binary neutron star (BNS) merger GW170817 was the first astrophysical source detected in gravitational waves and multi-wavelength electromagnetic radiation. The almost simultaneous observation of a pulse of gamma-rays proved that BNS mergers are associated with at least some short gamma-ray bursts (GRBs). However, the gamma-ray pulse was faint, casting doubts on the association of BNS mergers with the luminous, highly relativistic outflows of canonical short GRBs. Here we show that structured jets with a relativistic, energetic core surrounded by slower and less energetic wings produce afterglow emission that brightens characteristically with time, as recently seen in the afterglow of GW170817. Initially, we only see the relatively slow material moving towards us. As time passes, larger and larger sections of the outflow become visible, increasing the luminosity of the afterglow. The late appearance and increasing brightness of the multi-wavelength afterglow of GW170817 allow us to constrain the geometry of its ejecta and thus reveal the presence of an off-axis jet pointing about 30 degrees away from Earth. Our results confirm a single origin for BNS mergers and short GRBs: GW170817 produced a structured outflow with a highly relativistic core and a canonical short GRB. We did not see the bright burst because it was beamed away from Earth. However, approximately one in 20 mergers detected in gravitational waves will be accompanied by a bright, canonical short GRB.

We present the results of numerical simulations of the prompt emission of short-duration gamma-ray bursts. We consider emission from the relativistic jet, the mildly relativistic cocoon, and the non-relativistic shocked ambient material. We find that the cocoon material is confined between off-axis angles 15<theta<45 degrees and gives origin to X-ray transients with a duration of a few to ~10 seconds, delayed by a few seconds from the time of the merger. We also discuss the distance at which such transients can be detected, finding that it depends sensitively on the assumptions that are made about the radiation spectrum. Purely thermal cocoon transients are detectable only out to a few Mpc, Comptonized transients can instead be detected by the FERMI GBM out to several tens of Mpc.

We carry out a series of simulations of G2-type clouds interacting with the black hole at the galactic center, to determine why no large changes in the luminosity of Sgr A* were seen, and to determine the nature of G2. We measure the accretion rate from the gas cloud onto Sgr A* for a range of simulation parameters, such as cloud structure, background structure, background density, grid resolution, and accretion radius. For a broad range of parameters, the amount of cloud material accreted is small relative to the amount of background material accreted. The total accretion rate is not significantly effected for at least 30 years after periapsis. We find that reproducing observations of G2 likely requires two components for the object: an extended, cold gas cloud responsible for the Br-gamma emission, and a compact core or dusty stellar object dominating the bolometric luminosity. In simulations, the bolometric and X-ray luminosity have a peak lasting from about 1 year before to 1 year after periapsis, a feature not detected in observations. Our simulated Br-gamma emission is largely consistent with observations leading up to periapsis, with a slight increase in luminosity and a large increase in the FWHM of the line velocity. All emission from a gaseous component of G2 should fade rapidly after periapsis and be undetectable after 1 year, due to shock heating and expansion of the cloud. Any remaining emission should be from the compact component of G2.

We present calculations of the wide angle emission of short-duration gamma-ray bursts from compact binary merger progenitors. Such events are expected to be localized by their gravitational wave emission, fairly irrespective of the orientation of the angular momentum vector of the system, along which the gamma-ray burst outflow is expected to propagate. We show that both the prompt and afterglow emission are dim and challenging to detect for observers lying outside of the cone within which the relativistic outflow is propagating. If the jet initially propagates through a baryon contaminated region surrounding the merger site, however, a hot cocoon forms around it. The cocoon subsequently expands quasi-isotropically producing its own prompt emission and external shock powered afterglow. We show that the cocoon prompt emission is detectable by Swift BAT and Fermi GBM, We also show that the cocoon afterglow peaks a few hours to a few days after the burst and is detectable for up to a few weeks at all wavelengths. The timing and brightness of the transient are however uncertain due to their dependence on unknown quantities such as the density of the ambient medium surrounding the merger site, the cocoon energy, and the cocoon Lorentz factor. For a significant fraction of the gravitationally-detected neutron-star-binary mergers, the cocoon afterglow could possibly be the only identifiable electromagnetic counterpart, at least at radio and X-ray frequencies.

We model the possible afterglow of the Fermi GBM event associated with LIGO detection GW150914, under the assumption that the gamma-ray are produced by a short GRB-like relativistic outflow. We model GW150914-GBM as both a weak, on-axis short GRB and normal short GRB seen far off axis. Given the large uncertainty in the position of GW150914, we determine that the best chance of finding the afterglow is with ASKAP or possibly the MWA, with the flux from an off-axis short GRB reaching 0.2 - 4 mJy (0.12 - 16 mJy) at 150 MHz (863.5 MHz) by 1 - 12 months after the initial event. At low frequencies, the source would evolve from a hard to soft spectrum over several months. The radio afterglow would be detectable for several months to years after it peaks, meaning the afterglow may still be detectable and increasing in brightness NOW (mid-July 2016). With a localization from the MWA or ASKAP, the afterglow would be detectable at higher radio frequencies with the ATCA and in X-rays with Chandra or XMM.

Within the context of exoplanetary atmospheres, we present a comprehensive linear analysis of forced, damped, magnetized shallow water systems, exploring the effects of dimensionality, geometry (Cartesian, pseudo-spherical and spherical), rotation, magnetic tension and hydrodynamic and magnetic sources of friction. Across a broad range of conditions, we find that the key governing equation for atmospheres and quantum harmonic oscillators are identical, even when forcing (stellar irradiation), sources of friction (molecular viscosity, Rayleigh drag and magnetic drag) and magnetic tension are included. The global atmospheric structure is largely controlled by a single, key parameter that involves the Rossby and Prandtl numbers. This near-universality breaks down when either molecular viscosity or magnetic drag acts non-uniformly across latitude or a poloidal magnetic field is present, suggesting that these effects will introduce qualitative changes to the familiar chevron-shaped feature witnessed in simulations of atmospheric circulation. We also find that hydrodynamic and magnetic sources of friction have dissimilar phase signatures and affect the flow in fundamentally different ways, implying that using Rayleigh drag to mimic magnetic drag is inaccurate. We exhaustively lay down the theoretical formalism (dispersion relations, governing equations and time-dependent wave solutions) for a broad suite of models. In all situations, we derive the steady state of an atmosphere, which is relevant to interpreting infrared phase and eclipse maps of exoplanetary atmospheres. We elucidate a pinching effect that confines the atmospheric structure to be near the equator. Our suite of analytical models may be used to decisively develop physical intuition and as a reference point for three-dimensional, magnetohydrodynamic (MHD) simulations of atmospheric circulation.

Galaxy clusters contain much more metal per star, typically 3 times as much, than is produced in normal galaxies. We set out to determine what changes are needed to the stellar mass function and supernovae rates to account for this excess metal. In particular, we vary the Type Ia supernovae rate, IMF slope, upper and lower mass cutoffs, and the merger rate of massive stars. We then use existing simulation results for metal production from AGB stars, Type Ia SNe and core-collapse SNe to calculate the total amount of each element produced per solar mass of star formation. For models with very massive stars, we also include metal production from pair-instability supernovae (PISNe).We find that including PISNe makes it much easier to increase the amount of metal produced per stellar mass. Therefore a separate population of high-mass stars is not needed to produce the high amounts of metal found in galaxy clusters. We also find that including at least some PISNe increases the abundance of intermediate-mass elements relative to both oxygen and iron, consistent with observations of ICM abundances.

Low Mach number, high beta fast mode shocks can occur in the magnetic reconnection outflows of solar flares. These shocks, which occur above flare loop tops, may provide the electron energization responsible for some of the observed hard X-rays and contemporaneous radio emission. Here we present new 2D particle-in-cell simulations of low Mach number/high beta quasi-perpendicular shocks. The simulations show that electrons above a certain energy threshold experience shock-drift-acceleration. The transition energy between the thermal and non-thermal spectrum and the spectral index from the simulations are consistent with some of the X-ray spectra from RHESSI in the energy regime of $E\lesssim 40\sim 100$ keV. Plasma instabilities associated with the shock structure such as the modified-two-stream and the electron whistler instabilities are identified using numerical solutions of the kinetic dispersion relations. We also show that the results from PIC simulations with reduced ion/electron mass ratio can be scaled to those with the realistic mass ratio.

Astrophysical shocks are often studied in the high Mach number limit but weakly compressive fast shocks can occur in magnetic reconnection outflows and are considered to be a site of particle energization in solar flares. Here we study the microphysics of such perpendicular, low Mach number collisionless shocks using two-dimensional particle-in-cell (PIC) simulations with a reduced ion/electron mass ratio and employ a moving wall boundary method for initial generation of the shock. This moving wall method allows for more control of the shock speed, smaller simulation box sizes, and longer simulation times than the commonly used fixed wall, reflection method of shock formation. Our results, which are independent of the shock formation method, reveal the prevalence shock drift acceleration (SDA) of both electron and ions in a purely perpendicular shock with Alfv\'en Mach number $M_A=6.8$ and ratio of thermal to magnetic pressure $\beta=8$. We determine the respective minimum energies required for electrons and ions to incur SDA. We derive a theoretical electron distribution via SDA that compares to the simulation results. We also show that a modified two-stream instability due to the incoming and reflecting ions in the shock transition region acts as the mechanism to generate collisionless plasma turbulence that sustains the shock.

Magnetic reconnection is commonly perceived to drive flow and particle acceleration in flares of solar, stellar, and astrophysical disk coronae but the relative roles of different acceleration mecha- nisms in a given reconnection environment are not well understood. We show via direct numerical simulations that reconnection outflows produce weak fast shocks, when conditions for fast recon- nection are met and the outflows encounter an obstacle. The associated compression ratios lead to a Fermi acceleration particle spectrum that is significantly steeper than the strong fast shocks commonly studied, but consistent with the demands of solar flares. While this is not the only acceleration mechanism operating in a reconnection environment, it is plausibly a ubiquitous one.

We consider radiation emitted by the jitter mechanism in a Blandford-McKee self-similar blastwave. We assume the magnetic field configuration throughout the whole blastwave meets the condition for the emission of jitter radiation and we compute the ensuing images, light curves and spectra. The calculations are performed for both a uniform and a wind environment. We compare our jitter results to synchrotron results. We show that jitter radiation produces slightly different spectra than synchrotron, in particular between the self-absorption and the peak frequency, where the jitter spectrum is flat, while the synchrotron spectrum grows as \nu^{1/3}. The spectral difference is reflected in the early decay slope of the light curves. We conclude that jitter and synchrotron afterglows can be distinguished from each other with good quality observations. However, it is unlikely that the difference can explain the peculiar behavior of several recent observations, such as flat X-ray slopes and uncorrelated optical and X-ray behavior.

Accretion disks in which angular momentum transport is dominated by the magnetorotational instability (MRI) can also possess additional, purely hydrodynamic, drivers of turbulence. Even when the hydrodynamic processes, on their own, generate negligible levels of transport, they may still affect the evolution of the disk via their influence on the MRI. Here, we study the interaction between the MRI and hydrodynamic turbulence using local MRI simulations that include hydrodynamic forcing. As expected, we find that hydrodynamic forcing is generally negligible if it yields a saturated kinetic energy density that is small compared to the value generated by the MRI. For stronger hydrodynamic forcing levels, we find that hydrodynamic turbulence modifies transport, with the effect varying depending upon the spatial scale of hydrodynamic driving. Large scale forcing boosts transport by an amount that is approximately linear in the forcing strength, and leaves the character of the MRI (for example the ratio between Maxwell and Reynolds stresses) unchanged, up to the point at which the forced turbulence is an order of magnitude stronger than that generated by the MRI. Low amplitude small scale forcing may modestly suppress the MRI. We conclude that the impact of hydrodynamic turbulence on the MRI is generically ignorable in cases, such as convection, where the additional turbulence arises due to the accretion energy liberated by the MRI itself. Hydrodynamic turbulence may affect (and either enhance or suppress) the MRI if it is both strong, and driven by independent mechanisms such as self-gravity, supernovae, or solid-gas interactions in multiphase protoplanetary disks.

Relativistic electrons moving into a highly tangled magnetic field emit jitter radiation. We present a detailed computation of the jitter radiation spectrum, including self-absorption, for electrons inside Weibel-like shock generated magnetic fields. We apply our results to the case of the prompt and afterglow emission of gamma-ray bursts. We show that jitter emission can reproduce most of the observed features with some important differences with respect to standard synchrotron, especially in the frequency range between the self-absorption and the peak frequency. We discuss the similarities and differences between jitter and synchrotron and discuss experiments that can disentangle the two mechanisms.

The standard model of GRB afterglows assumes that the radiation observed as a delayed emission is of synchrotron origin, which requires the shock magnetic field to be relatively homogeneous on small scales. An alternative mechanism -- jitter radiation, which traditionally has been applied to the prompt emission -- substitutes synchrotron when the magnetic field is tangled on a microscopic scale. Such fields are produced at relativistic shocks by the Weibel instability. Here we explore the possibility that small-scale fields populate afterglow shocks. We derive the spectrum of jitter radiation under the afterglow conditions. We also derive the afterglow lightcurves for the ISM and Wind profiles of the ambient density. Jitter self-absorption is calculated here for the first time. We find that jitter radiation can produce afterglows similar to synchrotron-generated ones, but with some important differences. We compare the predictions of the two emission mechanisms. By fitting observational data to the synchrotron and jitter afterglow lightcurves, it can be possible to discriminate between the small-scale vs large-scale magnetic field models in afterglow shocks.