Massive galaxy clusters undergo strong evolution from z~1.6 to z~0.5, with overdense environments at high-z characterized by abundant dust-obscured star formation and stellar mass growth which rapidly give way to widespread quenching. Data spanning the near- to far-infrared (IR) spectrum can directly trace this transformation; however, such studies have largely been limited to the massive galaxy end of cluster populations. In this work, we present ``total light" stacking techniques spanning 3.4-500{\mu}m aimed at revealing the total cluster IR emission, including low mass members and potential intracluster dust. We detail our procedures for WISE, Spitzer, and Herschel imaging, including corrections to recover the total stacked emission in the case of high fractions of detected galaxies. We apply our stacking techniques to 232 well-studied massive (log M200/Msun~13.8) clusters across multiple z bins, recovering extended cluster emission at all wavelengths, typically at >5sigma. We measure the averaged near- to far-IR radial profiles and SEDs, quantifying the total stellar and dust content. The near-IR radial profiles are well described by an NFW model with a high (c~7) concentration parameter. Dust emission is similarly concentrated, albeit suppressed at small radii (r<0.2Mpc). The measured SEDs lack warm dust, consistent with the colder SEDs expected for low mass galaxies. We derive total stellar masses consistent with the theoretical Mhalo-M_star relation and specific-star formation rates that evolve strongly with redshift, echoing that of massive (log Mstar/Msun>10) cluster galaxies. Separating out the massive galaxy population reveals that the majority of cluster far-IR emission (~70-80%) is provided by the low mass constituents, which differs from field galaxies. This effect may be a combination of mass-dependent quenching and excess dust in low mass cluster galaxies.

Star formation depends critically on cooling mechanisms in the interstellar medium (ISM); however, thermal properties of gas in galaxies at the peak epoch of star formation (z ~ 2) remain poorly understood. A limiting factor in understanding the multiphase ISM is the lack of multiple tracers detected in the same galaxies, such as Polycyclic Aromatic Hydrocarbon (PAH) emission, a tracer of a critical photoelectric heating mechanism in interstellar gas, and [C II] 158{\mu}m fine-structure emission, a principal coolant. We present ALMA Band 9 observations targeting [C II] in six z ~ 2 star-forming galaxies with strong Spitzer IRS detections of PAH emission. All six galaxies are detected in dust continuum and marginally resolved. We compare the properties of PAH and [C II] emission, and constrain their relationship as a function of total infrared luminosity (LIR) and IR surface density. [C II] emission is detected in one galaxy at high signal-to-noise (34{\sigma}), and we place a secure upper limit on a second source. The rest of our sample are not detected in [C II] likely due to redshift uncertainties and narrow ALMA bandpass windows. Our results are consistent with the deficit in [C II]/LIR and PAH/LIR observed in the literature. However, the ratio of [C II] to PAH emission at z ~ 2 is possibly much lower than what is observed in nearby dusty star-forming galaxies. This could be the result of enhanced cooling via [O I] at high-z, hotter gas and dust temperatures, and/or a reduction in the photoelectric efficiency, in which the coupling between interstellar radiation and gas heating is diminished.

As some of the only Lyman continuum (LyC) emitters at z~0, Green Pea (GP) galaxies are possible analogs of the sources that reionized the universe. We present HST COS spectra of 13 of the most highly ionized GPs, with [O III]/[O II]=6-35, and investigate correlations between Ly-alpha, galaxy properties, and low-ionization UV lines. Galaxies with high [O III]/[O II] have higher H-alpha equivalent widths (EWs), and high intrinsic Ly-alpha production may explain the prevalence of high Ly-alpha EWs among GPs. While Ly-alpha escape fraction is closely linked to low gas covering fractions, implying a clumpy gas geometry, narrow Ly-alpha velocity peak separation (delta_v,LyA) correlates with the ionization state, suggesting a density-bounded geometry. We therefore suggest that delta_v,LyA may trace the residual transparency of low-column-density pathways. Metallicity is associated with both [O III]/[O II] and delta_v,LyA. This trend may result from catastrophic cooling around low-metallicity star clusters, which generates a compact geometry of dense clouds within a low-density inter-clump medium. We find that the relative strength of low-ionization UV emission to absorption correlates with Ly-alpha emission strength and is related to Ly-alpha profile shape. However, as expected for optically thin objects, the GPs with the lowest delta_v,LyA show both weak low-ionization emission and weak absorption. The strengths of the low-ionization absorption and emission lines in a stacked spectrum do not correspond to any individual spectrum. Galaxies with high [O III]/[O II] contain a high fraction of LyC emitter candidates, but [O III]/[O II] alone is an insufficient diagnostic of LyC escape.

Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems.

Mechanisms regulating the escape of Ly$\alpha$ photons and ionizing radiation remain poorly understood. To study these processes we analyze VLA 21cm observations of one Green Pea (GP), J160810+352809 (hereafter J1608), and HST COS spectra of 17 GP galaxies at $z<0.2$. All are highly ionized: J1608 has the highest [O III] $\lambda5007$/[O II] $\lambda3727$ for star-forming galaxies in SDSS, and the 17 GPs have [O III]/[O II] $\geq6.6$. We set an upper limit on J1608's HI mass of $\log M_{HI}/M_\odot=8.4$, near or below average compared to similar mass dwarf galaxies. In the COS sample, eight GPs show Ly$\alpha$ absorption components, six of which also have Ly$\alpha$ emission. The HI column densities derived from Ly$\alpha$ absorption are high, $\log N_{HI}/$cm$^{-2}=19-21$, well above the LyC optically thick limit. Using low-ionization absorption lines, we measure covering fractions ($f_{\mbox{cov}}$) of $0.1-1$, and find that $f_{\mbox{cov}}$ strongly anti-correlates with Ly$\alpha$ escape fraction. Low covering fractions may facilitate Ly$\alpha$ and LyC escape through dense neutral regions. GPs with $f_{\mbox{cov}}\sim1$ all have low neutral gas velocities, while GPs with lower $f_{\mbox{cov}}=0.2-0.6$ have a larger range of velocities. Conventional mechanical feedback may help establish low $f_{\mbox{cov}}$ in some cases, whereas other processes may be important for GPs with low velocities. Finally, we compare $f_{\mbox{cov}}$ with proposed indicators of LyC escape. Ionizing photon escape likely depends on a combination of neutral gas geometry and kinematics, complicating the use of emission-line diagnostics for identifying LyC emitters.

The classical, relativistic thin-disk theory of Novikov and Thorne (NT) predicts a maximum accretion efficiency of 40% for an optically thick, radiatively efficient accretion disk around a maximally spinning black hole (BH). However, when a strong magnetic field is introduced to numerical simulations of thin disks, large deviations in efficiency are observed, in part due to mass and energy carried by jets and winds launched by the disk or BH spin. The total efficiency of accretion can be significantly enhanced beyond that predicted by NT but it has remained unclear how the radiative component is affected. In order to study the effect of a dynamically relevant large-scale magnetic field on radiatively efficient accretion, we have performed numerical 3D general relativistic - radiative - magnetohydroynamic (GRRMHD) simulations of a disk with scale height to radius ratio of $H/R~0.1$ around a moderately spinning BH (a=0.5) using the code HARMRAD. Our simulations are fully global and allow us to measure the jet, wind, and radiative properties of a magnetically arrested disk (MAD) that is kept thin via self-consistent transport of energy by radiation using the M1 closure scheme. Our fiducial disk is MAD out to a radius of ~16R_g and the majority of the total ~13\% efficiency of the accretion flow is carried by a magnetically driven wind. We find that the radiative efficiency is slightly suppressed compared to NT, contrary to prior MAD GRMHD simulations with an ad hoc cooling function, but it is unclear how much of the radiation and thermal energy trapped in the outflows could ultimately escape.

In accretion disks with large-scale ordered magnetic fields, the magnetorotational instability (MRI) is marginally suppressed, so other processes may drive angular momentum transport leading to accretion. Accretion could then be driven by large-scale magnetic fields via magnetic braking, but large-scale magnetic flux can build-up onto the black hole and within the disk leading to a magnetically-arrested disk (MAD). Such a MAD state is unstable to the magnetic Rayleigh-Taylor (RT) instability, which itself leads to vigorous turbulence and the emergence of low-density highly-magnetized bubbles. This instability was studied in a thin (ratio of half-height H to radius R, $H/R \approx 0.1$) MAD simulation, where it has a more dramatic effect on the dynamics of the disk than for thicker disks. We find that the low-density bubbles created by the magnetic RT instability decrease the stress (leading to angular momentum transport) in the disk rather than increasing magnetic torques. Indeed, we find that the dominant component of the stress is due to turbulent magnetic fields, despite the suppression of the axisymmetric MRI and the dominant presence of large-scale magnetic fields. This suggests that the magnetic RT instability plays a significant role in driving angular momentum transport in MADs.

In the past few years wide-field optical and UV transient surveys as well as X-ray telescopes have allowed us to identify a few dozen candidate tidal disruption events (TDEs). While in theory the physical processes in TDEs are expected to be ubiquitous, a few distinct classes of TDEs have been observed. Some TDEs radiate mainly in NUV/optical while others produce prominent X-rays. Moreover, relativistic jets have been observed in only a handful of TDEs. This diversity might be related to the details of the super-Eddington accretion and emission physics relevant to TDE disks. In this Letter, we utilize novel three-dimensional general relativistic radiation magnetohydrodynamics simulations to study the super-Eddington compact disk phase expected in TDEs. Consistent with previous studies, geometrically thick disks, wide-angle optically-thick fast outflows and relativistic jets are produced. The outflow density and velocity depend sensitively on the inclination angle, and hence so does the reprocessing of emission produced from the inner disk. We then use Monte-Carlo radiative transfer to calculate the reprocessed spectra and find that that the observed ratio of optical to X-ray fluxes increases with increasing inclination angle. This naturally leads to a unified model for different classes of TDEs in which the spectral properties of the TDE depend mainly on the viewing-angle of the observer with respect to the orientation of the disk.

It is widely believed that relativistic jets in X-ray binaries and active-galactic nuclei are powered by the rotational energy of black holes. This idea is supported by general-relativistic magnetohydrodynamic (GRMHD) simulations of accreting black holes, which demonstrate efficient energy extraction via the Blandford-Znajek mechanism. However, due to uncertainties in the physics of mass-loading, and the failure of GRMHD numerical schemes in the highly-magnetized funnel region, the matter content of the jet remains poorly constrained. We investigate the observational signatures of mass-loading in the funnel by performing general-relativistic radiative transfer calculations on a range of 3D GRMHD simulations of accreting black holes. We find significant observational differences between cases in which the funnel is empty and cases where the funnel is filled with plasma, particularly in the optical and X-ray bands. In the context of Sgr A*, current spectral data constrains the jet filling only if the black hole is rapidly rotating with $a\gtrsim0.9$. In this case, the limits on the infrared flux disfavour a strong contribution from material in the funnel. We comment on the implications of our models for interpreting future Event Horizon Telescope observations. We also scale our models to stellar-mass black holes, and discuss their applicability to the low-luminosity state in X-ray binaries.

Blazars show variability on timescales ranging from minutes to years, the former being comparable to and in some cases even shorter than the light-crossing time of the central black hole. The observed gamma-ray lightcurves can be described by a power-law power density spectrum (PDS), with a similar index for both BL Lacs and flat-spectrum radio quasars. We show that this variability can be produced by turbulence in relativistic jets launched by magnetically arrested accretion flows (MADs). We perform radiative transport calculations on the turbulent, highly-magnetized jet launching region of a MAD with a rapidly rotating supermassive black hole. The resulting synchrotron and synchrotron self-Compton emission, originating from close to the black hole horizon, is highly variable. This variability is characterized by PDS which is remarkably similar to the observed power-law spectrum at frequencies less than a few per day. Furthermore, turbulence in the jet launching region naturally produces fluctuations in the plasma on scales much smaller than the horizon radius. We speculate that similar turbulent processes, operating in the jet at large radii (and therefore high bulk Lorentz factor), are responsible for blazar variability over many decades in frequency, including on minute timescales.

The James Webb Space Telescope's Medium Resolution Spectrometer (MRS), will offer nearly 2 orders of magnitude improvement in sensitivity and >3X improvement in spectral resolution over our previous space-based mid-IR spectrometer, the Spitzer IRS. In this paper, we make predictions for spectroscopic pointed observations and serendipitous detections with the MRS. Specifically, pointed observations of Herschel sources require only a few minutes on source integration for detections of several star-forming and active galactic nucleus lines, out to z$=$3 and beyond. But the same data will also include tens of serendipitous 0$\lesssim$z$\lesssim$4 galaxies per field with infrared luminosities ranging $\sim10^6-10^{13}$L$_{\odot}$. In particular, for the first time and for free we will be able to explore the $L_{IR}<10^{9}L_{\odot}$ regime out to $z\sim3$. We estimate that with $\sim$100 such fields, statistics of these detections will be sufficient to constrain the evolution of the low-$L$ end of the infrared luminosity function, and hence the star formation rate function. The above conclusions hold for a wide range in potential low-$L$ end of the IR luminosity function, and accounting for the PAH deficit in low-$L$, low-metallicity galaxies.

We present a numerical method which evolves a two-temperature, magnetized, radiative, accretion flow around a black hole, within the framework of general relativistic radiation magnetohydrodynamics. As implemented in the code KORAL, the gas consists of two sub-components -- ions and electrons -- which share the same dynamics but experience independent, relativistically consistent, thermodynamical evolution. The electrons and ions are heated independently according to a standard prescription from the literature for magnetohydrodynamical turbulent dissipation. Energy exchange between the particle species via Coulomb collisions is included. In addition, electrons gain and lose energy and momentum by absorbing and emitting synchrotron and bremsstrahlung radiation, and through Compton scattering. All evolution equations are handled within a fully covariant framework in the relativistic fixed-metric spacetime of the black hole. Numerical results are presented for five models of low luminosity black hole accretion. In the case of a model with a mass accretion rate $\dot{M}\sim10^{-8} \dot M_{\rm Edd}$, we find that radiation has a negligible effect on either the dynamics or the thermodynamics of the accreting gas. In contrast, a model with a larger $\dot{M}\sim 4\times 10^{-4} \dot M_{\rm Edd}$ behaves very differently. The accreting gas is much cooler and the flow is geometrically less thick, though it is not quite a thin accretion disk.

The possible Fermi detection of an electromagnetic counterpart to the double black hole merger GW150914 has inspired many theoretical models, some of which propose that the holes spiraled together inside a massive star. However, we show that the heat produced by the dynamical friction on such black hole orbits can exceed the stellar binding energy by a large factor, which means that this heat could destroy the star and thus make it difficult for enough gas to be near the holes at merger to produce detectable photons. These considerations must be taken into account when models are proposed for electromagnetic counterparts to the coalescence of two stellar-mass black holes. We find that only when the two black holes form very close to the center can the star avoid destruction. In that case, dynamical friction can make the black holes coalesce faster than they would in vacuum, which leads to a modification of the gravitational waveform that is potentially observable by advanced LIGO.

We present an extension to the general relativistic radiation magnetohydrodynamic code HARMRAD to account for emission and absorption by thermal cyclo-synchrotron, double Compton, bremsstrahlung, low-temperature OPAL opacities as well as Thomson and Compton scattering. We approximate the radiation field as a Bose-Einstein distribution and evolve it using the radiation number-energy-momentum conservation equations in order to track photon hardening. We perform various simulations to study how these extensions affect the radiative properties of magnetically-arrested disks accreting at Eddington to super-Eddington rates. We find that double Compton dominates bremsstrahlung in the disk within a radius of $r\sim 15r_g$ (gravitational radii) at a hundred times the Eddington accretion rate, and within smaller radii at lower accretion rates. Double Compton and cyclo-synchrotron regulate radiation and gas temperatures in the corona, while cyclo-synchrotron regulates temperatures in the jet. Interestingly, as the accretion rate drops to Eddington, an optically thin corona develops whose gas temperature of $T\sim 10^9$K is $\sim 100$ times higher than the disk's black body temperature. Our results show the importance of double Compton and synchrotron in super-Eddington disks, magnetized coronae, and jets.

Observations of jets in X-ray binaries show a correlation between radio power and black hole spin. This correlation, if confirmed, points towards the idea that relativistic jets may be powered by the rotational energy of black holes. In order to examine this further, we perform general-relativistic radiative transport calculations on magnetically arrested accretion flows, which are known to produce powerful jets via the Blandford-Znajek (BZ) mechanism. We find that the X-ray and gamma-ray emission strongly depend on spin and inclination angle. Surprisingly, the high-energy power does not show the same dependence on spin as the BZ jet power, but instead can be understood as a redshift effect. In particular, photons observed perpendicular to the spin axis suffer little net redshift until originating from close to the horizon. Such observers see deeper into the hot, dense, highly-magnetized inner disk region. This effect is largest for rapidly rotating black holes due to a combination of frame dragging and decreasing horizon radius. While the X-ray emission is dominated by the near horizon region, the near-infrared radiation originates at larger radii. Therefore, the ratio of X-ray to near-infrared power is an observational signature of black hole spin.

The radiative and jet efficiencies of thin magnetized accretion disks around black holes (BHs) are affected by BH spin and the presence of a magnetic field that, when strong, could lead to large deviations from Novikov-Thorne (NT) thin disk theory. To seek the maximum deviations, we perform general relativistic magnetohydrodynamic (GRMHD) simulations of radiatively efficient thin (half-height $H$ to radius $R$ of $H/R\approx 0.10$) disks around moderately rotating BHs with $a/M=0.5$. First, our simulations, each evolved for more than $70,000r_g/c$ (gravitational radius $r_g$ and speed of light $c$), show that large-scale magnetic field readily accretes inward even through our thin disk and builds-up to the magnetically-arrested disk (MAD) state. Second, our simulations of thin MADs show the disk achieves a radiative efficiency of $\eta_{\rm r}\approx 15\%$ (after estimating photon capture), which is about twice the NT value of $\eta_{\rm r}\sim 8\%$ for $a/M=0.5$ and gives the same luminosity as a NT disk with $a/M\approx 0.9$. Compared to prior simulations with $\lesssim 10\%$ deviations, our result of an $\approx 80\%$ deviation sets a new benchmark. Building on prior work, we are now able to complete an important scaling law which suggest that observed jet quenching in the high-soft state in BH X-ray binaries is consistent with an ever-present MAD state with a weak yet sustained jet.

Jets are observed as radio emission in active galactic nuclei and during the low/hard state in X-ray binaries (XRBs), but their contribution at higher frequencies has been uncertain. We study the dynamics of jets in XRBs using the general-relativistic magnetohydrodynamic code HARM. We calculate the high-energy spectra and variability properties using a general-relativistic radiative transport code based on grmonty. We find the following signatures of jet emission (i) a significant gamma-ray peak above ~10^22 Hz, (ii) a break in the optical/UV spectrum, with a change in luminosity from L ~ nu^0 to L ~ nu, followed by another break at higher frequencies where the spectrum roughly returns to L ~ nu^0, and (iii) a pronounced synchrotron peak near or below ~10^14 Hz indicates that a significant fraction of any observed X-ray emission originates in the jet. We investigate the variability during a large-scale magnetic field inversion in which the Blandford-Znajek (BZ) jet is quenched and a new transient hot reconnecting plasmoid is launched by the reconnecting field. The ratio of the gamma-rays to X-rays changes from L_gamma / L_X > 1 in the BZ jet to L_gamma / L_X < 1 during the launching of the transient plasmoid.

Magnetic fields are believed to drive accretion and relativistic jets in black hole accretion systems, but the magnetic-field structure that controls these phenomena remains uncertain. We perform general relativistic (GR) polarized radiative transfer of time-dependent three-dimensional GR magnetohydrodynamical (MHD) simulations to model thermal synchrotron emission from the Galactic Center source Sagittarius A$^\ast$ (Sgr A*). We compare our results to new polarimetry measurements by the Event Horizon Telescope (EHT) and show how polarization in the visibility (Fourier) domain distinguishes and constrains accretion flow models with different magnetic field structures. These include models with small-scale fields in disks driven by the magnetorotational instability (MRI) as well as models with large-scale ordered fields in magnetically-arrested disks (MAD). We also consider different electron temperature and jet mass-loading prescriptions that control the brightness of the disk, funnel-wall jet, and Blandford-Znajek-driven funnel jet. Our comparisons between the simulations and observations favor models with ordered magnetic fields near the black hole event horizon in Sgr A*, although both disk- and jet-dominated emission can satisfactorily explain most of the current EHT data. We show that stronger model constraints should be possible with upcoming circular polarization and higher frequency ($349~{\rm GHz}$) measurements.

Accretion discs and black holes (BHs) have angular momenta that are generally misaligned with respect to each other, which can lead to warps in the discs and bends in any jets produced. We consider a disc that is misaligned at large radii and torqued by Lense-Thirring (LT) precession and a Blandford-Znajek (BZ) jet torque. We consider a variety of disc states that include radiatively inefficient thick discs, radiatively efficient thin discs, and super-Eddington accretion discs. The magnetic field strength of the BZ jet is chosen as either from standard equipartition arguments or from magnetically arrested disc (MAD) simulations. We show that standard thin accretion discs can reach spin-disc alignment out to large radii long before LT would play a role, as caused by the slow infall time that gives even a weak BZ jet time to align the disc. We show that geometrically thick radiatively inefficient discs and super-Eddington discs in the MAD state reach disc-spin alignment near the black hole when density profiles are shallow as in magnetohydrodynamical simulations, while the BZ jet aligns discs with steep density profiles (as in advection-dominated accretion flows) with the BH spin out to larger radii. Our results imply that the BZ jet torque should affect the cosmological evolution of BH spin magnitude and direction, BH spin measurements in active galactic nuclei and X-ray binaries, and the interpretations for Event Horizon Telescope observations of discs or jets in strong-field gravity regimes.

Near a black hole, differential rotation of a magnetized accretion disk is thought to produce an instability that amplifies weak magnetic fields, driving accretion and outflow. These magnetic fields would naturally give rise to the observed synchrotron emission in galaxy cores and to the formation of relativistic jets, but no observations to date have been able to resolve the expected horizon-scale magnetic-field structure. We report interferometric observations at 1.3-millimeter wavelength that spatially resolve the linearly polarized emission from the Galactic Center supermassive black hole, Sagittarius A*. We have found evidence for partially ordered fields near the event horizon, on scales of ~6 Schwarzschild radii, and we have detected and localized the intra-hour variability associated with these fields.

This paper describes the Polarization Spectroscopic Telescope Array (PolSTAR), a mission proposed to NASA's 2014 Small Explorer (SMEX) announcement of opportunity. PolSTAR measures the linear polarization of 3-50 keV (requirement; goal: 2.5-70 keV) X-rays probing the behavior of matter, radiation and the very fabric of spacetime under the extreme conditions close to the event horizons of black holes, as well as in and around magnetars and neutron stars. The PolSTAR design is based on the technology developed for the Nuclear Spectroscopic Telescope Array (NuSTAR) mission launched in June 2012. In particular, it uses the same X-ray optics, extendable telescope boom, optical bench, and CdZnTe detectors as NuSTAR. The mission has the sensitivity to measure ~1% linear polarization fractions for X-ray sources with fluxes down to ~5 mCrab. This paper describes the PolSTAR design as well as the science drivers and the potential science return.

One of the puzzles associated with tidal disruption event candidates (TDEs) is that there is a dichotomy between the color temperatures of ${\rm few}\times 10^4$~K for TDEs discovered with optical and UV telescopes, and the color temperatures of ${\rm few}\times 10^5 - 10^6$~K for TDEs discovered with X-ray satellites. Here we propose that high-temperature TDEs are produced when the tidal debris of a disrupted star self-intersects relatively close to the supermassive black hole, in contrast to the more distant self-intersection that leads to lower color temperatures. In particular, we note from simple ballistic considerations that greater apsidal precession in an orbit is the key to closer self-intersection. Thus larger values of $\beta$, the ratio of the tidal radius to the pericenter distance of the initial orbit, are more likely to lead to higher temperatures of more compact disks which are super-Eddington and geometrically and optically thick. For a given star and $\beta$, apsidal precession also increases for larger black hole masses, but larger black hole masses imply a lower temperature at the Eddington luminosity. Thus the expected dependence of the temperature on the mass of the black hole is non-monotonic. We find that in order to produce a soft X-ray temperature TDE, a deep plunging stellar orbit with $\beta> 3$ is needed and a black hole mass of $\lesssim 5\times 10^6 M_\odot$ is favored. Although observations of TDEs are comparatively scarce and are likely dominated by selection effects, it is encouraging that both expectations are consistent with current data.

The radiative efficiency of super-Eddington accreting black holes (BHs) is explored for magnetically-arrested disks (MADs), where magnetic flux builds-up to saturation near the BH. Our three-dimensional general relativistic radiation magnetohydrodynamic (GRRMHD) simulation of a spinning BH (spin $a/M=0.8$) accreting at $\sim 50$ times Eddington shows a total efficiency $\sim 50\%$ when time-averaged and total efficiency $\gtrsim 100\%$ in moments. Magnetic compression by the magnetic flux near the rotating BH leads to a thin disk, whose radiation escapes via advection by a magnetized wind and via transport through a low-density channel created by a Blandford-Znajek (BZ) jet. The BZ efficiency is sub-optimal due to inertial loading of field lines by optically thick radiation, leading to BZ efficiency $\sim 40\%$ on the horizon and BZ efficiency $\sim 5\%$ by $r\sim 400r_g$ (gravitational radii) via absorption by the wind. Importantly, radiation escapes at $r\sim 400r_g$ with efficiency $\eta\approx 15\%$ (luminosity $L\sim 50L_{\rm Edd}$), similar to $\eta\approx 12\%$ for a Novikov-Thorne thin disk and beyond $\eta\lesssim 1\%$ seen in prior GRRMHD simulations or slim disk theory. Our simulations show how BH spin, magnetic field, and jet mass-loading affect the radiative and jet efficiencies of super-Eddington accretion.

We present a sub-grid model that emulates the magnetic dynamo operating in magnetized accretion disks. We have implemented this model in the general relativisic radiation magnetohydrodynamic (GRRMHD) code \koral, using results from local shearing sheet simulations of the magnetorotational instability to fix the parameters of the dynamo. With the inclusion of this dynamo, we are able to run 2D axisymmetric GRRMHD simulations of accretion disks for arbitrarily long times. The simulated disks exhibit sustained turbulence, with the poloidal and toroidal magnetic field components driven towards a state similar to that seen in 3D studies. Using this dynamo code, we present a set of long-duration global simulations of super-Eddington, optically-thick disks around non-spinning and spinning black holes. Super-Eddington disks around non-rotating black holes exhibit a surprisingly large efficiency, $\eta\approx0.04$, independent of the accretion rate, where we measure efficiency in terms of the total energy output, both radiation and mechanical, flowing out to infinity. Super-Eddington disks around spinning black holes are even more efficient, and appear to extract black hole rotational energy through a process similar to the Blandford-Znajek mechanism. All the simulated models are characterized by highly super-Eddington radiative fluxes collimated along the rotation axis. We also present a set of simulations that were designed to have Eddington or slightly sub-Eddington accretion rates ($\dot{M} \lesssim 2\dot M_{\rm Edd}$). None of these models reached a steady state. Instead, the disks collapsed as a result of runaway cooling, presumably because of a thermal instability.

Black hole (BH) accretion flows and jets are dynamic hot relativistic magnetized plasma flows whose radiative opacity can significantly affect flow structure and behavior. We describe a numerical scheme, tests, and an astrophysically relevant application using the M1 radiation closure within a new three-dimensional (3D) general relativistic (GR) radiation (R) magnetohydrodynamics (MHD) massively parallel code called HARMRAD. Our 3D GRRMHD simulation of super-Eddington accretion (about $20$ times Eddington) onto a rapidly rotating BH (dimensionless spin $j=0.9375$) shows sustained non-axisymmemtric disk turbulence, a persistent electromagnetic jet driven by the Blandford-Znajek effect, and a total radiative output consistently near the Eddington rate. The total accretion efficiency is of order $20\%$, the large-scale electromagnetic jet efficiency is of order $10\%$, and the total radiative efficiency that reaches large distances remains low at only order $1\%$. However, the radiation jet and the electromagnetic jet both emerge from a geometrically beamed polar region, with super-Eddington isotropic equivalent luminosities. Such simulations with HARMRAD can enlighten the role of BH spin vs.\ disks in launching jets, help determine the origin of spectral and temporal states in x-ray binaries, help understand how tidal disruption events (TDEs) work, provide an accurate horizon-scale flow structure for M87 and other active galactic nuclei (AGN), and isolate whether AGN feedback is driven by radiation or by an electromagnetic, thermal, or kinetic wind/jet. For example, the low radiative efficiency and weak BH spin-down rate from our simulation suggest that BH growth over cosmological times to billions of solar masses by redshifts of $z\sim 6-8$ is achievable even with rapidly rotating BHs and ten solar mass BH seeds.

Magnetically arrested accretion discs (MADs), where the magnetic pressure in the inner disc is dynamically important, provide an alternative mechanism for regulating accretion to what is commonly assumed in black hole systems. We show that a global magnetic field inversion in the MAD state can destroy the jet, significantly increase the accretion rate, and move the effective inner disc edge in to the marginally stable orbit. Reconnection of the MAD field in the inner radii launches a new type of transient outflow containing hot plasma generated by magnetic dissipation. This transient outflow can be as powerful as the steady magnetically-dominated Blandford-Znajek jet in the MAD state. The field inversion qualitatively describes many of the observational features associated with the high luminosity hard to soft state transition in black hole X-ray binaries: the jet line, the transient ballistic jet, and the drop in rms variability. These results demonstrate that the magnetic field configuration can influence the accretion state directly, and hence the magnetic field structure is an important second parameter in explaining observations of accreting black holes across the mass and luminosity scales.

A new general relativistic radiation magnetohydrodynamical code KORAL, is described, which employs the M1 scheme to close the radiation moment equations. The code has been successfully verified against a number of tests. Axisymmetric simulations of super-critical magnetized accretion on a non-rotating black hole (a=0.0) and a spinning black hole (a=0.9) are presented. The accretion rates in the two models are \dot M = 100-200 \dot M_Edd. These first general relativistic simulations of super-critical black hole accretion are potentially relevant to tidal disruption events and hyper-accreting supermassive black holes in the early universe. Both simulated models are optically and geometrically thick, and have funnels through which energy escapes in the form of relativistic gas, Poynting flux and radiative flux. The jet is significantly more powerful in the a=0.9 run. The net energy outflow rate in the two runs correspond to efficiencies of 5% (a=0) and 33% (a=0.9), as measured with respect to the mass accretion rate at the black hole. These efficiencies agree well with those measured in previous simulations of non-radiative geometrically thick disks. Furthermore, in the a=0.9 run, the outflow power appears to originate in the spinning black hole, suggesting that the associated physics is again similar in non-radiative and super-critical accretion flows. While the two simulations are efficient in terms of total energy outflow, both runs are radiatively inefficient. Their luminosities are only \sim 1-10 L_Edd, which corresponds to a radiative efficiency \sim 0.1%. Interestingly, most of the radiative luminosity emerges through the funnels, which subtend a very small solid angle. Therefore, measured in terms of a local radiative flux, the emitted radiation is highly super-Eddington.

High-frequency quasi-periodic oscillations (QPOs) appear in general-relativistic magnetohydrodynamic simulations of magnetically choked accretion flows around rapidly rotating black holes (BHs). We perform polarized radiative transfer calculations with ASTRORAY code to explore the manifestations of these QPOs for SgrA*. We construct a simulation-based model of a radiatively inefficient accretion flow and find model parameters by fitting the mean polarized source spectrum. The simulated QPOs have a total sub-mm flux amplitude up to 5% and a linearly polarized flux amplitude up to 2%. The oscillations reach high levels of significance 10-30\sigma and high quality factors Q~5. The oscillation period T~100M~35min corresponds to the rotation period of the BH magnetosphere that produces a trailing spiral in resolved disk images. The total flux signal is significant over noise for all tested frequencies 87GHz, 230GHz, and 857GHz and inclination angles 10deg, 37deg, and 80deg. The non-detection in the 230GHz Sub-Millimeter Array light curve is consistent with a low signal level and a low sampling rate. The presence of sub-mm QPOs in SgrA* will be better tested with the Atacama Large Millimeter Array.

We briefly summarize the method of simulating Sgr A* polarized sub-mm spectra from the accretion flow and fitting the observed spectrum. The dynamical flow model is based on three-dimensional general relativistic magneto hydrodynamic simulations. Fully self-consistent radiative transfer of polarized cyclo-synchrotron emission is performed. We compile a mean sub-mm spectrum of Sgr A* and fit it with the mean simulated spectra. We estimate the ranges of inclination angle theta=42-75deg, mass accretion rate Mdot=(1.4-7.0)*10^{-8}Msun/yr, and electron temperature Te=(3-4)*10^{10}K at 6M. We discuss multiple caveats in dynamical modeling, which must be resolved to make further progress.

Accreting black holes (BHs) produce intense radiation and powerful relativistic jets, which are affected by the BH's spin magnitude and direction. While thin disks might align with the BH spin axis via the Bardeen-Petterson effect, this does not apply to jet systems with thick disks. We used fully three-dimensional general relativistic magnetohydrodynamical simulations to study accreting BHs with various BH spin vectors and disk thicknesses with magnetic flux reaching saturation. Our simulations reveal a "magneto-spin alignment" mechanism that causes magnetized disks and jets to align with the BH spin near BHs and further away to reorient with the outer disk. This mechanism has implications for the evolution of BH mass and spin, BH feedback on host galaxies, and resolved BH images for SgrA* and M87.

We obtain estimates of Sgr A* accretion flow and black hole parameters by fitting polarized sub-mm observations with spectra computed using three-dimensional (3D) general relativistic (GR) magnetohydrodynamical (MHD) (GRMHD) simulations. Observations are compiled from averages over many epochs from reports in 29 papers for estimating the mean fluxes Fnu, linear polarization (LP) fractions, circular polarization (CP) fractions, and electric vector position angles (EVPAs). GRMHD simulations are computed with dimensionless spins a_*=0,0.5,0.7,0.9,0.98 over a 20,000M time interval. We perform fully self-consistent GR polarized radiative transfer using our new code to explore the effects of spin a_*, inclination angle \theta, position angle (PA), accretion rate Mdot, and electron temperature Te (Te is reported for radius 6M). By fitting the mean sub-mm fluxes and LP/CP fractions, we obtain estimates for these model parameters and determine the physical effects that could produce polarization signatures. Our best bet model has a_*=0.5, \theta=75deg, PA=115deg, Mdot=4.6*10^{-8}M_Sun/year, and Te=3.1*10^10K at 6M. The sub-mm CP is mainly produced by Faraday conversion as modified by Faraday rotation, and the emission region size at 230GHz is consistent with the VLBI size of 37microas. Across all spins, model parameters are in the ranges \theta=42deg-75deg, Mdot=(1.4-7.0)*10^{-8}M_Sun/year, and Te=(3-4)*10^10K. Polarization is found both to help differentiate models and to introduce new observational constraints on the effects of the magnetic field that might not be fit by accretion models so-far considered.

Black hole (BH) accretion flows and jets are qualitatively affected by the presence of ordered magnetic fields. We study fully three-dimensional global general relativistic magnetohydrodynamic (MHD) simulations of radially extended and thick (height $H$ to cylindrical radius $R$ ratio of $|H/R|\sim 0.2--1$) accretion flows around BHs with various dimensionless spins ($a/M$, with BH mass $M$) and with initially toroidally-dominated ($\phi$-directed) and poloidally-dominated ($R-z$ directed) magnetic fields. Firstly, for toroidal field models and BHs with high enough $|a/M|$, coherent large-scale (i.e. $\gg H$) dipolar poloidal magnetic flux patches emerge, thread the BH, and generate transient relativistic jets. Secondly, for poloidal field models, poloidal magnetic flux readily accretes through the disk from large radii and builds-up to a natural saturation point near the BH. For sufficiently high $|a/M|$ or low $|H/R|$ the polar magnetic field compresses the inflow into a geometrically thin highly non-axisymmetric "magnetically choked accretion flow" (MCAF) within which the standard linear magneto-rotational instability is suppressed. The condition of a highly-magnetized state over most of the horizon is optimal for the Blandford-Znajek mechanism that generates persistent relativistic jets with $\gtrsim 100$% efficiency for $|a/M|\gtrsim 0.9$. A magnetic Rayleigh-Taylor and Kelvin-Helmholtz unstable magnetospheric interface forms between the compressed inflow and bulging jet magnetosphere, which drives a new jet-disk quasi-periodic oscillation (JD-QPO) mechanism. The high-frequency QPO has spherical harmonic $|m|=1$ mode period of $\tau\sim 70GM/c^3$ for $a/M\sim 0.9$ with coherence quality factors $Q\gtrsim 10$. [abridged]

The outflow efficiency (eta) from black hole (BH) accretion disc systems is known to depend upon both the BH spin (a) and the amount of large-scale magnetic flux threading the BH and disc. Semi-analytical flux-trapping models suggest retrograde BHs should trap much more large-scale magnetic flux near the BH leading to much higher eta than for prograde BHs. We self-consistently determine the amount of large-scale magnetic flux trapped by rapidly spinning (a = -0.9 and 0.9) BHs using global 3D time-dependent non-radiative general relativistic magnetohydrodynamic simulations of thick (h/r ~ 0.3-0.6) discs. We find that BH-trapped flux builds up until it is strong enough to disrupt the inner accretion disc. Contrary to prior flux-trapping models, which do not include the back-reaction of magnetic flux on the disc, our simulations show prograde BHs trap more magnetic flux, leading to about 3 times higher eta than retrograde BHs for |a| = 0.9. Both spin orientations can produce highly efficient jets, eta ~ 100%, with increasing eta for increasing disc thickness. The similarity of eta for prograde and retrograde BHs makes it challenging to infer the sign of BH spin based on jet energetics alone.

Recent advances in general relativistic magnetohydrodynamic modeling of jets offer unprecedented insights into the inner workings of accreting black holes that power the jets in active galactic nuclei (AGN) and other accretion systems. I will present the results of recent studies that determine spin-dependence of jet power and discuss the implications for the AGN radio loud/quiet dichotomy and recent observations of high jet power in a number of AGN.

The supermassive black hole candidate at the center of M87 drives an ultra-relativistic jet visible on kiloparsec scales, and its large mass and relative proximity allow for event horizon scale imaging with very long baseline interferometry at millimeter wavelengths (mm-VLBI). Recently, relativistic magneto-hydrodynamic (MHD) simulations of black hole accretion flows have proven capable of launching magnetically-dominated jets. We construct time-dependent disc/jet models of the innermost portion of the M87 nucleus by performing relativistic radiative transfer calculations from one such simulation. We identify two types of models, jet-dominated or disc/jet, that can explain the spectral properties of M87, and use them to make predictions for current and future mm-VLBI observations. The Gaussian source size for the favored sky orientation and inclination from observations of the large-scale jet is 33-44 microarcseconds (~4-6 Schwarzschild radii) on current mm-VLBI telescopes, very similar to existing observations of Sgr A*. The black hole shadow, direct evidence of an event horizon, should be visible in future measurements using baselines between Hawaii and Mexico. Both models exhibit variability at millimeter wavelengths with factor of ~2 amplitudes on year timescales. For the low inclination of M87, the counter-jet dominates the event horizon scale millimeter wavelength emission from the jet-forming region.

Ongoing millimeter VLBI observations with the Event Horizon Telescope allow unprecedented study of the innermost portion of black hole accretion flows. Interpreting the observations requires relativistic, time-dependent physical modeling. We discuss the comparison of radiative transfer calculations from general relativistic MHD simulations of Sagittarius A* and M87 with current and future mm-VLBI observations. This comparison allows estimates of the viewing geometry and physical conditions of the Sgr A* accretion flow. The viewing geometry for M87 is already constrained from observations of its large-scale jet, but, unlike Sgr A*, there is no consensus for its millimeter emission geometry or electron population. Despite this uncertainty, as long as the emission region is compact, robust predictions for the size of its jet launching region can be made. For both sources, the black hole shadow may be detected with future observations including ALMA and/or the LMT, which would constitute the first direct evidence for a black hole event horizon.

We present an analytical solution for thin disk accretion onto a Kerr black hole that extends the standard Novikov-Thorne alpha-disk in three ways: (i) it incorporates nonzero stresses at the inner edge of the disk, (ii) it extends into the plunging region, and (iii) it uses a corrected vertical gravity formula. The free parameters of the model are unchanged. Nonzero boundary stresses are included by replacing the Novikov-Thorne no torque boundary condition with the less strict requirement that the fluid velocity at the innermost stable circular orbit is the sound speed, which numerical models show to be the correct behavior for luminosities below ~30% Eddington. We assume the disk is thin so we can ignore advection. Boundary stresses scale as alpha*h and advection terms scale as h^2 (where h is the disk opening angle (h=H/r)), so the model is self-consistent when h < alpha. We compare our solution with slim disk models and general relativistic magnetohydrodynamic disk simulations. The model may improve the accuracy of black hole spin measurements.

The "no hair" theorem, a key result in General Relativity, states that an isolated black hole is defined by only three parameters: mass, angular momentum, and electric charge; this asymptotic state is reached on a light-crossing time scale. We find that the "no hair" theorem is not formally applicable for black holes formed from collapse of a rotating neutron star. Rotating neutron stars can self-produce particles via vacuum breakdown forming a highly conducting plasma magnetosphere such that magnetic field lines are effectively "frozen-in" the star both before and during collapse. In the limit of no resistivity, this introduces a topological constraint which prohibits the magnetic field from sliding off the newly-formed event horizon. As a result, during collapse of a neutron star into a black hole, the latter conserves the number of magnetic flux tubes $N_B = e \Phi_\infty /(\pi c \hbar)$, where $\Phi_\infty \approx 2 \pi^2 B_{NS} R_{NS}^3 /(P_{\rm NS} c)$ is the initial magnetic flux through the hemispheres of the progenitor and out to infinity. We test this theoretical result via three-dimensional general relativistic plasma simulations of rotating black holes that start with a neutron star dipole magnetic field with no currents initially present outside the event horizon. The black hole's magnetosphere subsequently relaxes to the split monopole magnetic field geometry with self-generated currents outside the event horizon. The dissipation of the resulting equatorial current sheet leads to a slow loss of the anchored flux tubes, a process that balds the black hole on long resistive time scales rather than the short light-crossing time scales expected from the vacuum "no-hair" theorem.

Prompt gamma-ray burst (GRB) emission requires some mechanism to dissipate an ultrarelativistic jet. Internal shocks or some form of electromagnetic dissipation are candidate mechanisms. Any mechanism needs to answer basic questions, such as what is the origin of variability, what radius does dissipation occur at, and how does efficient prompt emission occur. These mechanisms also need to be consistent with how ultrarelativistic jets form and stay baryon pure despite turbulence and electromagnetic reconnection near the compact object and despite stellar entrainment within the collapsar model. We use the latest magnetohydrodynamical models of ultrarelativistic jets to explore some of these questions in the context of electromagnetic dissipation due to the slow collisional and fast collisionless reconnection mechanisms, as often associated with Sweet-Parker and Petschek reconnection, respectively. For a highly magnetized ultrarelativistic jet and typical collapsar parameters, we find that significant electromagnetic dissipation may be avoided until it proceeds catastrophically near the jet photosphere at large radii ($r\sim 10^{13}--10^{14}{\rm cm}$), by which the jet obtains a high Lorentz factor ($\gamma\sim 100--1000), has a luminosity of $L_j \sim 10^{50}--10^{51}\ergs$, has observer variability timescales of order 1s (ranging from 0.001-10s), achieves $\gamma\theta_j\sim 10--20 (for opening half-angle $\theta_j$) and so is able to produce jet breaks, and has comparable energy available for both prompt and afterglow emission. This reconnection switch mechanism allows for highly efficient conversion of electromagnetic energy into prompt emission and associates the observed prompt GRB pulse temporal structure with dissipation timescales of some number of reconnecting current sheets embedded in the jet.[abridged]

We describe global, 3D, time-dependent, non-radiative, general-relativistic, magnetohydrodynamic simulations of accreting black holes (BHs). The simulations are designed to transport a large amount of magnetic flux to the center, more than the accreting gas can force into the BH. The excess magnetic flux remains outside the BH, impedes accretion, and leads to a magnetically arrested disc. We find powerful outflows. For a BH with spin parameter a = 0.5, the efficiency with which the accretion system generates outflowing energy in jets and winds is eta ~ 30%. For a = 0.99, we find eta ~ 140%, which means that more energy flows out of the BH than flows in. The only way this can happen is by extracting spin energy from the BH. Thus the a = 0.99 simulation represents an unambiguous demonstration, within an astrophysically plausible scenario, of the extraction of net energy from a spinning BH via the Penrose-Blandford-Znajek mechanism. We suggest that magnetically arrested accretion might explain observations of active galactic nuclei with apparent eta ~ few x 100%.

The X-ray spectra of accretion discs of eight stellar-mass black holes have been analyzed to date using the thermal continuum fitting method, and the spectral fits have been used to estimate the spin parameters of the black holes. However, the underlying model used in this method of estimating spin is the general relativistic thin-disc model of Novikov & Thorne, which is only valid for razor-thin discs. We therefore expect errors in the measured values of spin due to inadequacies in the theoretical model. We investigate this issue by computing spectra of numerically calculated models of thin accretion discs around black holes, obtained via three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulations. We apply the continuum fitting method to these computed spectra to estimate the black hole spins and check how closely the values match the actual spin used in the GRMHD simulations. We find that the error in the dimensionless spin parameter is up to about 0.2 for a non-spinning black hole, depending on the inclination. For black holes with spins of 0.7, 0.9 and 0.98, the errors are up to about 0.1, 0.03 and 0.01 respectively. These errors are comparable to or smaller than those arising from current levels of observational uncertainty. Furthermore, we estimate that the GRMHD simulated discs from which these error estimates are obtained correspond to effective disc luminosities of about 0.4-0.7 Eddington, and that the errors will be smaller for discs with luminosities of 0.3 Eddington or less, which are used in the continuum-fitting method. We thus conclude that use of the Novikov-Thorne thin-disc model does not presently limit the accuracy of the continuum-fitting method of measuring black hole spin.

Magnetic reconnection, a fundamental plasma process associated with a rapid dissipation of magnetic energy, is believed to power many disruptive phenomena in laboratory plasma devices, the Earth magnetosphere, and the solar corona. Traditional reconnection research, geared towards these rather tenuous environments, has justifiably ignored the effects of radiation on the reconnection process. However, in many reconnecting systems in high-energy astrophysics (e.g., accretion-disk coronae, relativistic jets, magnetar flares) and, potentially, in powerful laser plasma and z-pinch experiments, the energy density is so high that radiation, in particular radiative cooling, may start to play an important role. This observation motivates the development of a theory of high-energy-density radiative magnetic reconnection. As a first step towards this goal, we present in this paper a simple Sweet--Parker-like theory of non-relativistic resistive-MHD reconnection with strong radiative cooling. First, we show how, in the absence of a guide magnetic field, intense cooling leads to a strong compression of the plasma in the reconnection layer, resulting in a higher reconnection rate. The compression ratio and the layer temperature are determined by the balance between ohmic heating and radiative cooling. The lower temperature in the radiatively-cooled layer leads to a higher Spitzer resistivity and hence to an extra enhancement of the reconnection rate. We then apply our general theory to several specific astrophysically important radiative processes (bremsstrahlung, cyclotron, and inverse-Compton) in the optically thin regime, for both the zero- and strong-guide-field cases. We derive specific expressions for key reconnection parameters, including the reconnection rate. We also discuss the limitations and conditions for applicability of our theory.

For the first time it has become possible to compare global 3D general relativistic magnetohydrodynamic (GRMHD) jet formation simulations directly to very-long baseline interferometric multi-frequency polarization observations of the pc-scale structure of active galactic nucleus (AGN) jets. Unlike the jet emission, which requires post hoc modeling of the non-thermal electrons, the Faraday rotation measures (RMs) depend primarily upon simulated quantities and thus provide a robust way in which to confront simulations with observations. We compute RM distributions of 3D global GRMHD jet formation simulations, with which we explore the dependence upon model and observational parameters, emphasizing the signatures of structures generic to the theory of MHD jets. With typical parameters, we find that it is possible to reproduce the observed magnitudes and many of the structures found in AGN jet RMs, including the presence of transverse RM gradients. In our simulations the RMs are generated within a smooth extension of the jet itself, containing ordered toroidally dominated magnetic fields. This results in a particular bilateral morphology that is unlikely to arise due to Faraday rotation in distant foreground clouds. However, critical to efforts to probe the Faraday screen will be resolving the transverse jet structure. Therefore, the RMs of radio cores may not be reliable indicators of the properties of the rotating medium. Finally, we are able to constrain the particle content of the jet, finding that at pc-scales AGN jets are electromagnetically dominated, with roughly 2% of the comoving energy in nonthermal leptons and much less in baryons.

The standard general relativistic model of a razor-thin accretion disk around a black hole, developed by Novikov & Thorne (NT) in 1973, assumes the shear stress vanishes at the radius of the innermost stable circular orbit (ISCO) and that, outside the ISCO, the shear stress is produced by an effective turbulent viscosity. However, astrophysical accretion disks are not razor-thin, it is uncertain whether the shear stress necessarily vanishes at the ISCO, and the magnetic field, which is thought to drive turbulence in disks, may contain large-scale structures that do not behave like a simple local scalar viscosity. We describe three-dimensional general relativistic magnetohydrodynamic simulations of accretion disks around black holes with a range of spin parameters, and we use the simulations to assess the validity of the NT model. Our fiducial initial magnetic field consists of multiple (alternating polarity) poloidal field loops whose shape is roughly isotropic in the disk in order to match the isotropic turbulence expected in the poloidal plane. For a thin disk with an aspect ratio |h/r| ~ 0.07 around a non-spinning black hole, we find a decrease in the accreted specific angular momentum of 2.9% relative to the NT model and an excess luminosity from inside the ISCO of 3.5%. The deviations in the case of spinning black holes are of the same order. In addition, the deviations decrease with decreasing |h/r|. We therefore conclude that magnetized thin accretion disks in x-ray binaries in the thermal/high-soft spectral state ought to be well-described by the NT model, especially at luminosities below 30% of Eddington where we expect a very small disk thickness |h/r| <~ 0.05. We also discuss how the stress and the luminosity inside the ISCO depend on the assumed initial magnetic field geometry. (abridged)

Recent high resolution observations of the Galactic center black hole allow for direct comparison with accretion disk simulations. We compare two-temperature synchrotron emission models from three dimensional, general relativistic magnetohydrodynamic simulations to millimeter observations of Sgr A*. Fits to very long baseline interferometry and spectral index measurements disfavor the monochromatic face-on black hole shadow models from our previous work. Inclination angles \le 20 degrees are ruled out to 3 \sigma. We estimate the inclination and position angles of the black hole, as well as the electron temperature of the accretion flow and the accretion rate, to be i=50+35-15 degrees, \xi=-23+97-22 degrees, T_e=(5.4 +/- 3.0)x10^10 K and Mdot=(5+15-2)x10^-9 M_sun / yr respectively, with 90% confidence. The black hole shadow is unobscured in all best fit models, and may be detected by observations on baselines between Chile and California, Arizona or Mexico at 1.3mm or .87mm either through direct sampling of the visibility amplitude or using closure phase information. Millimeter flaring behavior consistent with the observations is present in all viable models, and is caused by magnetic turbulence in the inner radii of the accretion flow. The variability at optically thin frequencies is strongly correlated with that in the accretion rate. The simulations provide a universal picture of the 1.3mm emission region as a small region near the midplane in the inner radii of the accretion flow, which is roughly isothermal and has \nu/\nu_c ~ 1-20, where \nu_c is the critical frequency for thermal synchrotron emission.

Achromatic breaks in afterglow light curves of gamma-ray bursts (GRBs) arise naturally if the product of the jet's Lorentz factor \gamma and opening angle \Theta_j satisfies (\gamma \Theta_j) >> 1 at the onset of the afterglow phase, i.e., soon after the conclusion of the prompt emission. Magnetohydrodynamic (MHD) simulations of collimated GRB jets generally give (\gamma \Theta_j) <~ 1, suggesting that MHD models may be inconsistent with jet breaks. We work within the collapsar paradigm and use axisymmetric relativistic MHD simulations to explore the effect of a finite stellar envelope on the structure of the jet. Our idealized models treat the jet-envelope interface as a collimating rigid wall, which opens up outside the star to mimic loss of collimation. We find that the onset of deconfinement causes a burst of acceleration accompanied by a slight increase in the opening angle. In our fiducial model with a stellar radius equal to 10^4.5 times that of the central compact object, the jet achieves an asymptotic Lorentz factor \gamma ~ 500 far outside the star and an asymptotic opening angle \Theta_j ~ 0.04 rad ~ 2 deg, giving (\gamma \Theta_j) ~ 20. These values are consistent with observations of typical long-duration GRBs, and explain the occurrence of jet breaks. We provide approximate analytic solutions that describe the numerical results well.

Radio loud active galactic nuclei (AGN) are on average 1000 times brighter in the radio band compared to radio quiet AGN. We investigate whether this radio loud/quiet dichotomy can be due to differences in the spin of the central black holes that power the radio-emitting jets. Using general relativistic magnetohydrodynamic simulations, we construct steady state axisymmetric numerical models for a wide range of black hole spins (dimensionless spin parameter 0.1 <= a <= 0.9999 and a variety of jet geometries. We assume that the total magnetic flux through the black hole horizon at radius r_H(a) is held constant. If the black hole is surrounded by a thin accretion disk, we find that the total black hole power output depends approximately quadratically on the angular frequency of the hole, P \propto \Omega_H^2 \propto (a/r_H)^2. We conclude that, in this scenario, differences in the black hole spin can produce power variations of only a few tens at most. However, if the disk is thick such that the jet subtends a narrow solid angle around the polar axis, then the power dependence becomes much steeper, P \propto \Omega_H^4 or even \propto \Omega_H^6. Power variations of 1000 are then possible for realistic black hole spin distributions. We derive an analytic solution that accurately reproduces the steeper scaling of jet power with \Omega_H, and we provide a numerical fitting formula that reproduces all our simulation results. We discuss other physical effects that might contribute to the observed radio loud/quiet dichotomy of AGN.

The connection between collimation and acceleration of magnetized relativistic jets is discussed. The focus is on recent numerical simulations which shed light on some longstanding problems.

We have developed time-dependent models of FU Ori accretion outbursts to explore the physical properties of protostellar disks. Our two-dimensional, axisymmetric models incorporate full vertical structure with a new treatment of the radiative boundary condition for the disk photosphere. We find that FU Ori-type outbursts can be explained by a slow accumulation of matter due to gravitational instability. Eventually this triggers the magnetorotational instability, which leads to rapid accretion. The thermal instability is triggered in the inner disk but this instability is not necessary for the outburst. An accurate disk vertical structure, including convection, is important for understanding the outburst behavior. Large convective eddies develop during the high state in the inner disk. The models are in agreement with Spitzer IRS spectra and also with peak accretion rates and decay timescales of observed outbursts, though some objects show faster rise timescale. We also propose that convection may account for the observed mild-supersonic turbulence and the short-timescale variations of FU Orionis objects.

Unconfined relativistic outflows from rotating, magnetized compact objects are often well-modeled by assuming the field geometry is approximately a split-monopole at large radii. Earlier work has indicated that such an unconfined flow has an inefficient conversion of magnetic energy to kinetic energy. This has led to the conclusion that ideal magnetohydrodynamical (MHD) processes fail to explain observations of, e.g., the Crab pulsar wind at large radii where energy conversion appears efficient. In addition, as a model for astrophysical jets, the monopole field geometry has been abandoned in favor of externally confined jets since the latter appeared to be generically more efficient jet accelerators. We perform time-dependent axisymmetric relativistic MHD simulations in order to find steady state solutions for a wind from a compact object endowed with a monopole field geometry. Our simulations follow the outflow for 10 orders of magnitude in distance from the compact object, which is large enough to study both the initial "acceleration zone" of the magnetized wind as well as the asymptotic "coasting zone." We obtain the surprising result that acceleration is actually {\it efficient} in the polar region, which develops a jet despite not being confined by an external medium. Our models contain jets that have sufficient energy to account for moderately energetic long and short gamma-ray burst (GRB) events (~10^{51}--10^{52} erg), collimate into narrow opening angles (opening half-angle \theta_j \approx 0.03 rad), become matter-dominated at large radii (electromagnetic energy flux per unit matter energy flux \sigma<1), and move at ultrarelativistic Lorentz factors (\gamma_j ~ 200 for our fiducial model). (abridged)

Rotating magnetized compact objects and their accretion discs can generate strong toroidal magnetic fields driving highly magnetized plasmas into relativistic jets. Of significant concern, however, has been that a strong toroidal field in the jet should be highly unstable to the non-axisymmetric helical kink (screw) $m=1$ mode leading to rapid disruption. In addition, a recent concern has been that the jet formation process itself may be unstable due to the accretion of non-dipolar magnetic fields. We describe large-scale fully three-dimensional global general relativistic magnetohydrodynamic simulations of rapidly rotating, accreting black holes producing jets. We study both the stability of the jet as it propagates and the stability of the jet formation process during accretion of dipolar and quadrupolar fields. For our dipolar model, despite strong non-axisymmetric disc turbulence, the jet reaches Lorentz factors of $\Gamma\sim 10$ with opening half-angle $\theta_j\sim 5^\circ$ at $10^3$ gravitational radii without significant disruption or dissipation with only mild substructure dominated by the $m=1$ mode. On the contrary, our quadrupolar model does not produce a steady relativistic ($\Gamma\gtrsim 3$) jet due to mass-loading of the polar regions caused by unstable polar fields. Thus, if produced, relativistic jets are roughly stable structures and may reach up to an external shock with strong magnetic fields. We discuss the astrophysical implications of the accreted magnetic geometry playing such a significant role in relativistic jet formation, and we outline avenues for future work.

We describe three-dimensional general relativistic magnetohydrodynamic simulations of a geometrically thin accretion disk around a non-spinning black hole. The disk has a thickness $h/r\sim0.05-0.1$ over the radial range $(2-20)GM/c^2$. In steady state, the specific angular momentum profile of the inflowing magnetized gas deviates by less than 2% from that of the standard thin disk model of Novikov & Thorne (1973). Also, the magnetic torque at the radius of the innermost stable circular orbit (ISCO) is only $\sim2%$ of the inward flux of angular momentum at this radius. Both results indicate that magnetic coupling across the ISCO is relatively unimportant for geometrically thin disks.

Long-duration gamma-ray bursts (GRBs) require an engine capable of driving a jet of plasma to ultrarelativistic bulk Lorentz factors of up to several hundred and into narrow opening angles of a few degrees. We use global axisymmetric stationary solutions of magnetically-dominated (force-free) ultrarelativistic jets to test whether the popular magnetic-driving paradigm can generate the required Lorentz factors and opening angles. Our global solutions are obtained via time-dependent relativistic ideal magnetodynamical numerical simulations which follow the jet from the central engine to beyond six orders of magnitude in radius. Our model is primarily motivated by the collapsar model, in which a jet is produced by a spinning black hole or neutron star and then propagates through a massive stellar envelope. We find that the size of the presupernova progenitor star and the radial profile of pressure inside the star determine the terminal Lorentz factor and opening angle of the jet. At the radius where the jet breaks out of the star, our well-motivated fiducial model generates a Lorentz factor $\gamma\sim 400$ and a half-opening angle $\theta_j\sim 2^\circ$, consistent with observations of many long-duration GRBs. Other models with slightly different parameters give $\gamma$ in the range 100 to 5000 and $\theta_j$ from $0.1^\circ$ to $10^\circ$, thus reproducing the range of properties inferred for GRB jets. A potentially observable feature of some of our solutions is that the maximum Poynting flux in the jet is found at $\theta \sim \theta_j$ with the jet power concentrated in a hollow cone, while the maximum in the Lorentz factor occurs at an angle $\theta$ substantially smaller than $\theta_j$ also in a hollow cone. [abridged]

Active galactic nuclei, x-ray binaries, pulsars, and gamma-ray bursts are all believed to be powered by compact objects surrounded by relativistic plasma flows driving phenomena such as accretion, winds, and jets. These flows are often accurately modelled by the relativistic magnetohydrodynamics (MHD) approximation. Time-dependent numerical MHD simulations have proven to be especially insightful, but one regime that remains difficult to simulate is when the energy scales (kinetic, thermal, magnetic) within the plasma become disparate. We develop a numerical scheme that significantly improves the accuracy and robustness of the solution in this regime. We use a modified form of the WENO method to construct a finite-volume general relativistic hydrodynamics code called WHAM that converges at fifth order. We avoid (1) field-by-field decomposition by adaptively reducing down to 2-point stencils near discontinuities for a more accurate treatment of shocks, and (2) excessive reduction to low order stencils, as in the standard WENO formalism, by maintaining high order accuracy in smooth monotonic flows. Our scheme performs the proper surface integral of the fluxes, converts cell averaged conserved quantities to point conserved quantities before performing the reconstruction step, and correctly averages all source terms. We demonstrate that the scheme is robust in strong shocks, very accurate in smooth flows, and maintains accuracy even when the energy scales in the flow are highly disparate.

The role of the equation of state for a perfectly conducting, relativistic magnetized fluid is the main subject of this work. The ideal constant $\Gamma$-law equation of state, commonly adopted in a wide range of astrophysical applications, is compared with a more realistic equation of state that better approximates the single-specie relativistic gas. The paper focus on three different topics. First, the influence of a more realistic equation of state on the propagation of fast magneto-sonic shocks is investigated. This calls into question the validity of the constant $\Gamma$-law equation of state in problems where the temperature of the gas substantially changes across hydromagnetic waves. Second, we present a new inversion scheme to recover primitive variables (such as rest-mass density and pressure) from conservative ones that allows for a general equation of state and avoids catastrophic numerical cancellations in the non-relativistic and ultrarelativistic limits. Finally, selected numerical tests of astrophysical relevance (including magnetized accretion flows around Kerr black holes) are compared using different equations of state. Our main conclusion is that the choice of a realistic equation of state can considerably bear upon the solution when transitions from cold to hot gas (or viceversa) are present. Under these circumstances, a polytropic equation of state can significantly endanger the solution.

The expulsion of axisymmetric magnetic field from the event horizons of rapidly rotating black holes has been seen as an astrophysically important effect that may significantly reduce or even nullify the efficiency of the Blandford-Znajek mechanism of powering the relativistic jets in Active Galactic Nuclei and Gamma Ray Bursts. However, this Meissner-like effect is seen in vacuum solutions of black hole electrodynamics whereas the Blandford-Znajek mechanism is concerned with plasma-filled magnetospheres. In this paper we argue that conductivity dramatically changes the properties of axisymmetric electromagnetic solutions -- even for a maximally rotating Kerr black hole the magnetic field is pulled inside the event horizon. Moreover, the conditions resulting in outgoing Poynting flux in the Blandford-Znajek mechanism exist not on the event horizon but everywhere within the black hole ergosphere. Thus, the Meissner effect is unlikely to be of interest in astrophysics of black holes, at least not in the way this has been suggested so far. These conclusions are supported by the results of time-dependent numerical simulations with three different computer codes. The test problems involve black holes with the rotation parameter ranging from $a=0.999$ to $a=1$. The pure electrodynamic simulations deal with the structure of conductive magnetospheres of black holes placed in a uniform-at-infinity magnetic field (Wald's problem) and the magnetohydrodynamic simulations are used to study the magnetospheres arising in the problem of disk accretion.

In paper I, we showed that time-dependent general relativistic magnetohydrodynamic (GRMHD) numerical models of accretion disks, although being highly turbulent, have surprisingly simple electromagnetic properties. In particular, the toroidal current density in the disk takes the form $dI_\phi/dr \propto r^{-5/4}$. Guided by this simplicity, we use a time-dependent general relativistic force-free electrodynamics (GRFFE) code to study an idealized problem in which the accretion disk is replaced by an infinitely thin rotating equatorial current sheet. We consider both an $r^{-5/4}$ current profile and an $r^{-1}$ profile, the latter corresponding to the paraboloidal model of Blandford & Znajek (1977). The force-free magnetosphere we obtain with the $r^{-5/4}$ current sheet matches remarkably well to the Poynting-dominated jet seen in GRMHD numerical models. By comparing to the non-rotating force-free model studied in paper I, rotation is seen to lead to mild decollimation of the jet suggesting that hoop-stress forces nearly cancel centrifugal forces. In order to study the process that generates the corona and disk wind and destroys the ordered field in the corona in GRMHD numerical models, the force-free field with the $r^{-5/4}$ current distribution is embedded in an accretion disk and followed in a GRMHD simulation. Reconnection and magnetic stresses contribute to a magnetized, thermal wind without the aid of an ordered field threading the disk.

General relativistic numerical simulations of magnetized accretion flows around black holes show a disordered electromagnetic structure in the disk and corona and a highly relativistic, Poynting-dominated funnel jet in the polar regions. The polar jet is nearly consistent with the stationary paraboloidal Blandford-Znajek model of an organized field threading the polar regions of a rotating black hole. How can a disordered accretion disk and corona lead to an ordered jet? We show that the polar jet is associated with a strikingly simple angular-integrated toroidal current distribution $dI_\phi/dr \propto r^{-5/4}$, where $I_\phi(r)$ is the toroidal current enclosed inside radius $r$. We demonstrate that the poloidal magnetic field in the simulated jet agrees well with the force-free field solution for a non-rotating thin disk with an $r^{-5/4}$ toroidal current, suggesting rotation leads to negligible self-collimation. We find that the polar field is confined/collimated by the corona. The electromagnetic field in the disk also scales as $r^{-5/4}$, which is consistent with some Newtonian accretion models that assume rough equipartition between magnetic and gas pressure. However, the agreement is accidental since toward the black hole the magnetic pressure increases faster than the gas pressure. This field dominance near the black hole is associated with magnetic stresses that imply a large effective viscosity parameter $\alpha\sim 1$, whereas the typically assumed value of $\alpha\sim 0.1$ holds far from the black hole.[abridged]

We consider a self-similar force-free wind flowing out of an infinitely thin disk located in the equatorial plane. On the disk plane, we assume that the magnetic stream function $P$ scales as $P\propto R^\nu$, where $R$ is the cylindrical radius. We also assume that the azimuthal velocity in the disk is constant: $v_\phi = Mc$, where $M<1$ is a constant. For each choice of the parameters $\nu$ and $M$, we find an infinite number of solutions that are physically well-behaved and have fluid velocity $\leq c$ throughout the domain of interest. Among these solutions, we show via physical arguments and time-dependent numerical simulations that the minimum-torque solution, i.e., the solution with the smallest amount of toroidal field, is the one picked by a real system. For $\nu \geq 1$, the Lorentz factor of the outflow increases along a field line as $\gamma \approx M(z/\Rfp)^{(2-\nu)/2} \approx R/R_{\rm A}$, where $\Rfp$ is the radius of the foot-point of the field line on the disk and $R_{\rm A}=\Rfp/M$ is the cylindrical radius at which the field line crosses the Alfven surface or the light cylinder. For $\nu < 1$, the Lorentz factor follows the same scaling for $z/\Rfp < M^{-1/(1-\nu)}$, but at larger distances it grows more slowly: $\gamma \approx (z/\Rfp)^{\nu/2}$. For either regime of $\nu$, the dependence of $\gamma$ on $M$ shows that the rotation of the disk plays a strong role in jet acceleration. On the other hand, the poloidal shape of a field line is given by $z/\Rfp \approx (R/\Rfp)^{2/(2-\nu)}$ and is independent of $M$. Thus rotation has neither a collimating nor a decollimating effect on field lines, suggesting that relativistic astrophysical jets are not collimated by the rotational winding up of the magnetic field.

The formation and large-scale propagation of Poynting-dominated jets produced by accreting, rapidly rotating black hole systems are studied by numerically integrating the general relativistic magnetohydrodynamic equations of motion to follow the self-consistent interaction between accretion disks and black holes. This study extends previous similar work by studying jets till $t\approx 10^4GM/c^3$ out to $r\approx 10^4GM/c^2$, by which the jet is super- fast magnetosonic and moves at a lab-frame bulk Lorentz factor of $\Gamma\sim 10$ with a maximum terminal Lorentz factor of $\Gamma_\infty\lesssim 10^3$. The radial structure of the Poynting-dominated jet is piece-wise self-similar, and fits to flow quantities along the field line are provided. Beyond the \alf surface at $r\sim 10$--$100GM/c^2$, the jet becomes marginally unstable to (at least) current-driven instabilities. Such instabilities drive shocks in the jet that limit the efficiency of magnetic acceleration and collimation. These instabilities also induce jet substructure with $3\lesssim\Gamma\lesssim 15$. The jet is shown to only marginally satisfy the necessary and sufficient conditions for kink instability, so this may explain how astrophysical jets can extend to large distances without completely disrupting. At large distance, the jet angular structure is Gaussian-like (or uniform within the core with sharp exponential wings) with a half-opening angle of $\approx 5^\circ$ and there is an extended component out to $\approx 27^\circ$. Unlike in some hydrodynamic simulations, the environment is found to play a negligible role in jet structure, acceleration, and collimation as long as the ambient pressure of the surrounding medium is small compared to the magnetic pressure in the jet.

The luminosity and structure of neutron star magnetospheres are crucial to our understanding of pulsar and plerion emission. A solution found using the force-free approximation would be an interesting standard with which any model with more physics could be compared. Prior quasi-analytic force-free solutions may not be stable, while prior time-dependent magnetohydrodynamic models used unphysical model parameters. We use a time-dependent relativistic force-free electrodynamics code with no free parameters to find a unique stationary solution for the axisymmetric rotating pulsar magnetosphere in a Minkowski space-time in the case of no surface currents on the star. The solution is similar to the force-free quasi-analytic solution of \citet{cont99} and the numerical magnetohydrodynamic solution of \citet{kom05}. The magnetosphere structure and the usefulness of the classical y-point in the general dissipative regime are discussed. The pulsar luminosity is found to be $L \approx 0.99\pm 0.01 \mu^2\Omega_\star^4/c^3$ for a dipole moment $\mu$ and stellar angular frequency $\Omega_\star$.

The force-free limit of magnetohydrodynamics (MHD) is often a reasonable approximation to model black hole and neutron star magnetospheres. We describe a general relativistic force-free (GRFFE) formulation that allows general relativistic magnetohydrodynamic (GRMHD) codes to directly evolve the GRFFE equations of motion. Established, accurate, and well-tested conservative GRMHD codes can simply add a new inversion piece of code to their existing code, while continuing to use all the already-developed facilities present in their GRMHD code. We show how to enforce the $\mathbf{E}\cdot\mathbf{B}=0$ constraint and energy conservation, and we introduce a simplified general model of the dissipation of the electric field to enforce the $B^2-E^2>0$ constraint. We also introduce a simplified yet general method to resolve current sheets, without much reconnection, over many dynamical times. This formulation is incorporated into an existing GRMHD code (HARM), which is demonstrated to give accurate and robust GRFFE results for Minkowski and black hole space-times.

Conservative numerical schemes for general relativistic magnetohydrodynamics (GRMHD) require a method for transforming between ``conserved'' variables such as momentum and energy density and ``primitive'' variables such as rest-mass density, internal energy, and components of the four-velocity. The forward transformation (primitive to conserved) has a closed-form solution, but the inverse transformation (conserved to primitive) requires the solution of a set of five nonlinear equations. Here we discuss the mathematical properties of the inverse transformation and present six numerical methods for performing the inversion. The first method solves the full set of five nonlinear equations directly using a Newton-Raphson scheme and a guess from the previous timestep. The other methods reduce the five nonlinear equations to either one or two nonlinear equations that are solved numerically. Comparisons between the methods are made using a survey over phase space, a two-dimensional explosion problem, and a general relativistic MHD accretion disk simulation. The run-time of the methods is also examined. Code implementing the schemes is available for download on the web.

A rotating black hole probably powers ultrarelativistic jets in gamma-ray bursts, relativistic jets from some active galactic nuclei, and jets from some black hole x-ray binaries. Prior estimates of the power output of a black hole have assumed an infinitely thin disk, a magnetic field based upon a slowly rotating black hole, and have not self-consistently determined the geometry or magnitude of the magnetic field for a realistic accretion disk. We provide useful formulae for the total and jet Blandford-Znajek (BZ) power and efficiency as determined self-consistently from general relativistic magnetohydrodynamic numerical models. Of all jet mechanisms, we suggest that only the BZ mechanism is able to produce an ultrarelativistic jet.

Two types of relativistic jets are suggested to form near accreting black holes: a potentially ultrarelativistic Poynting-dominated jet and a Poynting-baryon jet. One source of jet matter is electron-positron pair production, which is driven by neutrino annihilation in GRBs and photon annihilation in AGN and x-ray binaries. GRB Poynting-dominated jets are also loaded by electron-proton pairs by the collisional cascade of Fick-diffused free neutrons. We show that, for the collapsar model, the neutrino-driven enthalpy flux (classic fireball model) is probably dominated by the Blandford-Znajek energy flux, which predicts a jet Lorentz factor of $\Gamma\sim 100-1000$. We show that radiatively inefficient AGN, such as M87, are synchrotron-cooling limited to $\Gamma\sim 2-10$. Radiatively efficient x-ray binaries, such as GRS1915+105, are Compton-drag limited to $\Gamma \lesssim 2$, but the jet may be destroyed by Compton drag. However, the Poynting-baryon jet is a collimated outflow with $\Gamma \sim 1-3$. The jet from radiatively efficient systems, such as microquasar GRS1915+105, may instead be a Poynting-baryon jet that is only relativistic when the disk is geometrically thick. In a companion paper, general relativistic hydromagnetic simulations of black hole accretion with pair creation are used to simulate jet formation in GRBs, AGN, and x-ray binaries.

In a companion theory paper, we presented a unified model of jet formation. We suggested that primarily two types of relativistic jets form near accreting black holes: a potentially ultrarelativistic Poynting-dominated jet and a Poynting-baryon jet. We showed that, for the collapsar model, the neutrino-driven enthalpy flux (classic fireball model) is probably dominated by the Blandford-Znajek energy flux, which predicts a jet Lorentz factor of $\Gamma\sim 100-1000$. We showed that radiatively inefficient AGN, such as M87, are synchrotron-cooling limited to $\Gamma\sim 2-10$. Radiatively efficient x-ray binaries, such as GRS1915+105, are Compton-drag limited to $\Gamma \lesssim 2$, but the jet may be destroyed by Compton drag. However, the Poynting-baryon jet is a collimated outflow with $\Gamma \sim 1-3$. Here we present general relativistic hydromagnetic simulations of black hole accretion with pair creation used to simulate jet formation in GRBs, AGN, and x-ray binaries. Our collapsar model shows the development of a patchy ``magnetic fireball'' with typically $\Gamma\sim 100-1000$ and a Gaussian structure. Temporal variability of the jet is dominated by toroidal field instabilities for $\gtrsim 10^2$ gravitational radii. A broader Poynting-baryon jet with $\Gamma\sim 1.5$ could contribute to a supernova.

Some active galactic nuclei, microquasars, and gamma ray bursts may be powered by the electromagnetic braking of a rapidly rotating black hole. We investigate this possibility via axisymmetric numerical simulations of a black hole surrounded by a magnetized plasma. The plasma is described by the equations of general relativistic magnetohydrodynamics, and the effects of radiation are neglected. The evolution is followed for $2000 G M/c^3$, and the computational domain extends from inside the event horizon to typically $40 G M/c^2$. We compare our results to two analytic steady state models, including the force-free magnetosphere of Blandford & Znajek. Along the way we present a self-contained rederivation of the Blandford-Znajek model in Kerr-Schild (horizon penetrating) coordinates. We find that (1) low density polar regions of the numerical models agree well with the Blandford-Znajek model; (2) many of our models have an outward Poynting flux on the horizon in the Kerr-Schild frame; (3) none of our models have a net outward energy flux on the horizon; and (4) one of our models, in which the initial disk has net magnetic flux, shows a net outward angular momentum flux on the horizon. We conclude with a discussion of the limitations of our model, astrophysical implications, and problems to be addressed by future numerical experiments.

We consider astrophysical maser radiation that is created in the presence of mildly supersonic, magnetohydrodynamic (MHD) turbulence. The focus is on the OH masers for which the magnetic field is strong enough that the separations of the Zeeman components are greater than the spectral linebreadths. A longstanding puzzle has been the absence of the Zeeman $\pi $ components and the high circular polarization in the observed spectra of these masers. We first argue that the elongation of eddies along the field that has recently been recognized in MHD turbulence will enhance the optical depth parallel to the magnetic field in comparison with that perpendicular to the magnetic field. We then simulate maser emission with a numerical model of MHD turbulence to demonstrate quantitatively how the intensities of the linearly polarized $\pi $ components are suppressed and the intensities of the nearly circularly polarized $\sigma $ components are enhanced. This effect is also generic in the sense that most spectral lines in MHD turbulence with Mach number $M\sim 1$ should have larger optical depths parallel to the magnetic field than perpendicular. The effect is reduced considerably when $M<1$. The simulations also demonstrate that the velocity and magnetic field variations due to the turbulence can (but do not necessarily) cause one of the $\sigma $ components to be much more intense than the other, as is often observed for mainline OH masers.

We consider a subset of the physical processes that determine the spin j = a/M of astrophysical black holes. These include: (1) Initial conditions. Recent models suggest that the collapse of supermassive stars are likely to produce black holes with j ~ 0.7. (2) Major mergers. The outcome of a nearly equal mass black hole-black hole merger is not yet known, but we review the current best guesses and analytic bounds. (3) Minor mergers. We recover the result of Blandford & Hughes that accretion of small companions with isotropically distributed orbital angular momenta results in spindown, with j ~ M^{-7/3}. (4) Accretion. We present new results from fully relativistic magnetohydrodynamic accretion simulations. These show that, at least for one sequence of flow models, spin equilibrium (dj/dt = 0) is reached for j ~ 0.9, far less than the canonical value 0.998 of Thorne that was derived in the absence of MHD effects. This equilibrium value may not apply to all accretion flows, particularly thin disks. Nevertheless, it opens the possibility that black holes that have grown primarily through accretion are not maximally rotating.

We describe a conservative, shock-capturing scheme for evolving the equations of general relativistic magnetohydrodynamics. The fluxes are calculated using the Harten, Lax, and van Leer scheme. A variant of constrained transport, proposed earlier by T\'oth, is used to maintain a divergence free magnetic field. Only the covariant form of the metric in a coordinate basis is required to specify the geometry. We describe code performance on a full suite of test problems in both special and general relativity. On smooth flows we show that it converges at second order. We conclude by showing some results from the evolution of a magnetized torus near a rotating black hole.

We report on a numerical study of viscous fluid accretion onto a black hole. The flow is axisymmetric and uses a pseudo-Newtonian potential to model relativistic effects near the event horizon. The numerical method is a variant of the ZEUS code. As a test of our numerical scheme, we are able to reproduce results from earlier, similar work by Igumenshchev and Abramowicz and Stone et al. We consider models in which mass is injected onto the grid as well as models in which an initial equilibrium torus is accreted. In each model we measure three ``eigenvalues'' of the flow: the accretion rate of mass, angular momentum, and energy. We find that the eigenvalues are sensitive to r_{in}, the location of the inner radial boundary. Only when the flow is always supersonic on the inner boundary are the eigenvalues insensitive to small changes in r_{in}. We also report on the sensitivity of the results to other numerical parameters.