We consider the effects of black hole tilt on accretion disk alignment, studying three initial black hole tilts,6 degrees, 12 degrees and 24 degrees, with both magnetohydrodynamic and (inviscid) hydrodynamic evolution. In a number of ways, but not all, the dynamics are homologous in the sense that the alignment fronts resulting from different initial tilts are very similar when analyzed in terms of the fraction of the initial tilt angle. Even when the initial misalignment is 24 degrees, which, for the sound speed studied, is 4 vertical scale heights at the disk fiducial radius, the surface density remains a smooth function of radius; i.e., we find no examples in which the disk inner aligned and outer misaligned regions separate, or "break".

We present a series of simulations in both pure hydrodynamics (HD) and magnetohydrodynamics (MHD) exploring the degree to which alignment of disks subjected to external precessional torques (e.g., as in the `Bardeen-Petterson' effect) is dependent upon the disk sound speed c_s. Across the range of sound speeds examined, we find that the influence of the sound speed can be encapsulated in a simple "lumped-parameter" model proposed by Sorathia et al. (2013a). In this model, alignment fronts propagate outward at a speed ~0.2 rOmega_precess(r), where Omega_precess is the local test-particle precession frequency. Meanwhile, transonic radial motions transport angular momentum both inward and outward at a rate that may be described roughly in terms of an orientation diffusion model with diffusion coefficient ~2c_s^2/Omega, for local orbital frequency Omega. The competition between the two leads, in isothermal disks, to a stationary position for the alignment front at a radius proportional to c_s^(-4/5). For alignment to happen at all, the disk must either be turbulent due to the magnetorotational instability in MHD, or, in HD, it must be cool enough for the bending waves driven by disk warp to be nonlinear at their launch point. Contrary to long-standing predictions, warp propagation in MHD disks is diffusive independent of the parameter c_s/(alpha v_orb$, for orbital speed v_orb and ratio of stress to pressure of alpha. In purely HD disks, i.e., those with no internal stresses other than bulk viscosity, warmer disks align weakly or not at all; cooler disks align qualitatively similarly to MHD disks.

Magnetic fields are fundamental to the dynamics of both accretion disks and the jets that they often drive. We review the basic physics of these phenomena, the past and current efforts to model them numerically with an emphasis on the jet-disk connection, and the observational constraints on the role of magnetic fields in the jets of active galaxies on all scales.

Using only physical mechanisms, i.e., 3D MHD with no phenomenological viscosity, we have simulated the dynamics of a moderately thin accretion disk subject to torques whose radial scaling mimics those produced by lowest-order post-Newtonian gravitomagnetism. In this simulation, we have shown how, in the presence of MHD turbulence, a time-steady transition can be achieved between an inner disk region aligned with the equatorial plane of the central mass's spin and an outer region orbiting in a different plane. The position of the equilibrium orientation transition is determined by a balance between gravitomagnetic torque and warp-induced inward mixing of misaligned angular momentum from the outer disk. If the mixing is interpreted in terms of diffusive transport, the implied diffusion coefficient is ~(0.6--0.8)c_s^2/Omega for sound speed c_s and orbital frequency Omega. This calibration permits estimation of the orientation transition's equilibrium location given the central mass, its spin parameter, and the disk's surface density and scaleheight profiles. However, the alignment front overshoots before settling into an equilibrium, signaling that a diffusive model does not fully represent the time-dependent properties of alignment fronts under these conditions. Because the precessional torque on the disk at the alignment front is always comparable to the rate at which misaligned angular momentum is brought inward to the front by warp-driven radial motions, no break forms between the inner and outer portions of the disk in our simulation. Our results also raise questions about the applicability to MHD warped disks of the traditional distinction between "bending wave" and "diffusive" regimes.

When matter orbits around a central mass obliquely with respect to the mass's spin axis, the Lense-Thirring effect causes it to precess at a rate declining sharply with radius. Ever since the work of Bardeen & Petterson (1975), it has been expected that when a fluid fills an orbiting disk, the orbital angular momentum at small radii should then align with the mass's spin. Nearly all previous work has studied this alignment under the assumption that a phenomenological "viscosity" isotropically degrades fluid shears in accretion disks, even though it is now understood that internal stress in flat disks is due to anisotropic MHD turbulence. In this paper we report a pair of matched simulations, one in MHD and one in pure (non-viscous) HD in order to clarify the specific mechanisms of alignment. As in the previous work, we find that disk warps induce radial flows that mix angular momentum of different orientation; however, we also show that the speeds of these flows are generically transonic and are only very weakly influenced by internal stresses other than pressure. In particular, MHD turbulence does not act in a manner consistent with an isotropic viscosity. When MHD effects are present, the disk aligns, first at small radii and then at large; alignment is only partial in the HD case. We identify the specific angular momentum transport mechanisms causing alignment and show how MHD effects permit them to operate more efficiently. Lastly, we relate the speed at which an alignment front propagates outward (in the MHD case) to the rate at which Lense-Thirring torques deliver angular momentum at smaller radii.

Global disk simulations provide a powerful tool for investigating accretion and the underlying magnetohydrodynamic turbulence driven by the magneto-rotational instability (MRI). Using them to predict accurately quantities such as stress, accretion rate, and surface brightness profile requires that purely numerical effects, arising from both resolution and algorithm, be understood and controlled. We use the flux-conservative Athena code to conduct a series of experiments on disks having a variety of magnetic topologies to determine what constitutes adequate resolution. We develop and apply several resolution metrics: Qz and Qphi, the ratio of the grid zone size to the characteristic MRI wavelength, alpha_mag, the ratio of the Maxwell stress to the magnetic pressure, and the ratio of radial to toroidal magnetic field energy. For the initial conditions considered here, adequate resolution is characterized by Qz > 15, Qphi > 20, alpha_mag = 0.45, and a field energy ratio of 0.2. These values are associated with > 35 zones per scaleheight, a result consistent with shearing box simulations. Numerical algorithm is also important. Use of the HLLE flux solver or second-order interpolation can significantly degrade the effective resolution compared to the HLLD flux solver and third-order interpolation. Resolution at this standard can be achieved only with large numbers of grid zones, arranged in a fashion that matches the symmetries of the problem and the scientific goals of the simulation.

Orbiting disks may exhibit bends due to a misalignment between the angular momentum of the inner and outer regions of the disk. We begin a systematic simulational inquiry into the physics of warped disks with the simplest case: the relaxation of an unforced warp under pure fluid dynamics, i.e. with no internal stresses other than Reynolds stress. We focus on the nonlinear regime in which the bend rate is large compared to the disk aspect ratio. When warps are nonlinear, strong radial pressure gradients drive transonic radial motions along the disk's top and bottom surfaces that efficiently mix angular momentum. The resulting nonlinear decay rate of the warp increases with the warp rate and the warp width, but, at least in the parameter regime studied here, is independent of the sound speed. The characteristic magnitude of the associated angular momentum fluxes likewise increases with both the local warp rate and the radial range over which the warp extends; it also increases with increasing sound speed, but more slowly than linearly. The angular momentum fluxes respond to the warp rate after a delay that scales with the square-root of the time for sound waves to cross the radial extent of the warp. These behaviors are at variance with a number of the assumptions commonly used in analytic models to describe linear warp dynamics.

We present the first three-dimensional magnetohydrodynamic (MHD) simulations of a circumbinary disk surrounding an equal mass binary. The binary maintains a fixed circular orbit of separation $a$. As in previous hydrodynamical simulations, strong torques by the binary can maintain a gap of radius $\simeq 2a$. Streams curve inward from $r \simeq 2a$ toward the binary; some of their mass passes through the inner boundary, while the remainder swings back out to the disk. However, we also find that near its inner edge the disk develops both a strong $m=1$ asymmetry and growing orbital eccentricity. Because the MHD stresses introduce more matter into the gap, the total torque per unit disk mass is $\simeq 14$ times larger than found previously. The inner boundary accretion rate per unit disk mass is $\simeq 40$ times greater than found from previous hydrodynamical calculations. The implied binary shrinkage rate is determined by a balance between the rate at which the binary gains angular momentum by accretion and loses it by gravitational torque. The large accretion rate brings these two rates nearly into balance, but in net, we find that $\dot a/a < 0$, and its magnitude is about 2.7 times larger than predicted by the earlier hydrodynamic simulations. If the binary comprises two massive black holes, the accretion rate may be great enough for one or both to be AGN, with consequences for the physical state of the gas both in the disk body and in its inner gap.

Recent general relativistic magneto-hydrodynamic (MHD) simulations of accretion onto black holes have shown that, contrary to the basic assumptions of the Novikov-Thorne model, there can be substantial magnetic stress throughout the plunging region. Additional dissipation and radiation can therefore be expected. We use data from a particularly well-resolved simulation of accretion onto a non-spinning black hole to compute both the radiative efficiency of such a flow and its spectrum if all emitted light is radiated with a thermal spectrum whose temperature matches the local effective temperature. This disk is geometrically thin enough (H/r ~= 0.06) that little heat is retained in the flow. In terms of light reaching infinity (i.e., after allowance for all relativistic effects and for photon capture by the black hole), we find that the radiative efficiency is at least ~=6-10% greater than predicted by the Novikov-Thorne model (complete radiation of all heat might yield another ~6%). We also find that the spectrum more closely resembles the Novikov-Thorne prediction for a/M ~= 0.2--0.3 than for the correct value, a/M=0. As a result, if the spin of a non-spinning black hole is inferred by model-fitting to a Novikov-Thorne model with known black hole mass, distance, and inclination, the inferred a/M is too large by ~= 0.2--0.3.

Discretized numerical simulations are a powerful tool for investigation of nonlinear MHD turbulence in accretion disks. However, confidence in their quantitative predictions requires a demonstration that further refinement of the spatial gridscale would not result in any significant change. This has yet to be accomplished, particularly for global disk simulations. In this paper, we combine data from previously published stratified shearing box simulations and new global disk simulations to calibrate several quantitative diagnostics by which one can estimate progress toward numerical convergence. Using these diagnostics, we find that the established criterion for an adequate numerical description of linear growth of the magneto-rotational instability (the number of cells across a wavelength of the fastest-growing vertical wavenumber mode) can be extended to a criterion for adequate description of nonlinear MHD disk turbulence, but the standard required is more stringent. We also find that azimuthal resolution, which has not often been extensively examined in previous studies, can significantly affect the evolution of the poloidal magnetic field. We further analyze the comparative resolution requirements of a small sample of initial magnetic field geometries; not surprisingly, more complicated initial field geometries require higher spatial resolution. Otherwise, they tend to evolve to qualitatively similar states if evolved for sufficient time. Applying our quantitative resolution criteria to a sample of previously published global simulations, we find that, with perhaps a single exception, they are significantly under-resolved, and therefore underestimate the magnetic turbulence and resulting stress levels throughout the accretion flow.

We investigate the effect of shear viscosity and Ohmic resistivity on the magnetorotational instability (MRI) in vertically stratified accretion disks through a series of local simulations with the Athena code. First, we use a series of unstratified simulations to calibrate physical dissipation as a function of resolution and background field strength; the effect of the magnetic Prandtl number, Pm = viscosity/resistivity, on the turbulence is captured by ~32 grid zones per disk scale height, H. In agreement with previous results, our stratified disk calculations are characterized by a subthermal, predominately toroidal magnetic field that produces MRI-driven turbulence for |z| < 2 H. Above |z| = 2 H, magnetic pressure dominates and the field is buoyantly unstable. Large scale radial and toroidal fields are also generated near the mid-plane and subsequently rise through the disk. The polarity of this mean field switches on a roughly 10 orbit period in a process that is well-modeled by an alpha-omega dynamo. Turbulent stress increases with Pm but with a shallower dependence compared to unstratified simulations. For sufficiently large resistivity, on the order of cs H/1000, where cs is the sound speed, MRI turbulence within 2 H of the mid-plane undergoes periods of resistive decay followed by regrowth. This regrowth is caused by amplification of toroidal field via the dynamo. This process results in large amplitude variability in the stress on 10 to 100 orbital timescales, which may have relevance for partially ionized disks that are observed to have high and low accretion states.

We explore the parameter dependence of inner disk stress in black hole accretion by contrasting the results of a number of simulations, all employing 3-d general relativistic MHD in a Schwarzschild spacetime. Five of these simulations were performed with the intrinsically conservative code HARM3D, which allows careful regulation of the disk aspect ratio, H/R; our simulations span a range in H/R from 0.06 to 0.17. We contrast these simulations with two previously reported simulations in a Schwarzschild spacetime in order to investigate possible dependence of the inner disk stress on magnetic topology. In all cases, much care was devoted to technical issues: ensuring adequate resolution and azimuthal extent, and averaging only over those time-periods when the accretion flow is in approximate inflow equilibrium. We find that the time-averaged radial-dependence of fluid-frame electromagnetic stress is almost completely independent of both disk thickness and poloidal magnetic topology. It rises smoothly inward at all radii (exhibiting no feature associated with the ISCO) until just outside the event horizon, where the stress plummets to zero. Reynolds stress can also be significant near the ISCO and in the plunging region; the magnitude of this stress, however, depends on both disk thickness and magnetic topology. The two stresses combine to make the net angular momentum accreted per unit rest-mass 7-15% less than the angular momentum of the ISCO.

Resistivity and viscosity have a significant role in establishing the energy levels in turbulence driven by the magnetorotational instability (MRI) in local astrophysical disk models. This study uses the Athena code to characterize the effects of a constant shear viscosity \nu and Ohmic resistivity \eta in unstratified shearing box simulations with a net toroidal magnetic flux. A previous study of shearing boxes with zero net magnetic field performed with the ZEUS code found that turbulence dies out for values of the magnetic Prandtl number, P_m = \nu/\eta, below P_m \sim 1; for P_m \gtrsim 1, time- and volume-averaged stress levels increase with P_m. We repeat these experiments with Athena and obtain consistent results. Next, the influence of viscosity and resistivity on the toroidal field MRI is investigated both for linear growth and for fully-developed turbulence. In the linear regime, a sufficiently large \nu or \eta can prevent MRI growth; P_m itself has little direct influence on growth from linear perturbations. By applying a range of values for \nu and \eta to an initial state consisting of fully developed turbulence in the presence of a background toroidal field, we investigate their effects in the fully nonlinear system. Here, increased viscosity enhances the turbulence, and the turbulence decays only if the resistivity is above a critical value; turbulence can be sustained even when P_m < 1, in contrast to the zero net field model. While we find preliminary evidence that the stress converges to a small range of values when \nu and \eta become small enough, the influence of dissipation terms on MRI-driven turbulence for relatively large \eta and \nu is significant, independent of field geometry.

We report on a global, three-dimensional GRMHD simulation of an accretion torus embedded in a large scale vertical magnetic field orbiting a Schwarzschild black hole. This simulation investigates how a large scale vertical field evolves within a turbulent accretion disk and whether global magnetic field configurations suitable for launching jets and winds can develop. We find that a "coronal mechanism" of magnetic flux motion, which operates largely outside the disk body, dominates global flux evolution. In this mechanism, magnetic stresses driven by orbital shear create large-scale half-loops of magnetic field that stretch radially inward and then reconnect, leading to discontinuous jumps in the location of magnetic flux. In contrast, little or no flux is brought in directly by accretion within the disk itself. The coronal mechanism establishes a dipole magnetic field in the evacuated funnel around the orbital axis with a field intensity regulated by a combination of the magnetic and gas pressures in the inner disk. These results prompt a reevaluation of previous descriptions of magnetic flux motion associated with accretion. Local pictures are undercut by the intrinsically global character of magnetic flux. Formulations in terms of an "effective viscosity" competing with an "effective resistivity" are undermined by the nonlinearity of of the magnetic dynamics and the fact that the same turbulence driving mass motion (traditionally identified as "viscosity") can alter magnetic topology.

Magnetic fields connecting the immediate environs of rotating black holes to large distances appear to be the most promising mechanism for launching relativistic jets, an idea first developed by Blandford and Znajek in the mid-1970s. To enable an understanding of this process, we provide a brief introduction to dynamics and electromagnetism in the spacetime near black holes. We then present a brief summary of the classical Blandford-Znajek mechanism and its conceptual foundations. Recently, it has become possible to study these effects in much greater detail using numerical simulation. After discussing which aspects of the problem can be handled well by numerical means and which aspects remain beyond the grasp of such methods, we summarize their results so far. Simulations have confirmed that processes akin to the classical Blandford-Znajek mechanism can launch powerful electromagnetically-dominated jets, and have shown how the jet luminosity can be related to black hole spin and concurrent accretion rate. However, they have also shown that the luminosity and variability of jets can depend strongly on magnetic field geometry. We close with a discussion of several important open questions.

(abridged) We apply a second-order Godunov code, Athena, to studies of the magnetorotational instability using unstratified shearing box simulations with a uniform net vertical field and a sinusoidally varying zero net vertical field. The Athena results agree well with similar studies that used different numerical algorithms. We conduct analyses to study the flow of energy from differential rotation to turbulent fluctuations to thermalization. A study of the temporal correlation between the time derivatives of volume-averaged energy components shows that energy injected into turbulent fluctuations dissipates on a timescale of $\Omega^{-1}$, where $\Omega$ is the orbital frequency of the local domain. Magnetic dissipation dominates over kinetic dissipation, although not by as great a factor as the ratio of magnetic to kinetic energy. We Fourier-transform the magnetic and kinetic energy evolution equations and, using the assumption that the time-averaged energies are constant, determine the level of numerical dissipation as a function of length scale and resolution. By modeling numerical dissipation as if it were physical in origin, we characterize numerical resistivity and viscosity in terms of effective Reynolds and Prandtl numbers. The resulting effective magnetic Prandtl number is $\sim 2$, independent of resolution or initial field geometry. MRI simulations with effective Reynolds and Prandtl numbers determined by numerical dissipation are not equivalent to those where these numbers are set by physical resistivity and viscosity. These results serve, then, as a baseline for future shearing box studies where dissipation is controlled by the inclusion of explicit viscosity and resistivity.

Numerical simulation of magnetohydrodynamic (MHD) turbulence makes it possible to study accretion dynamics in detail. However, special effort is required to connect inflow dynamics (dependent largely on angular momentum transport) to radiation (dependent largely on thermodynamics and photon diffusion). To this end we extend the flux-conservative, general relativistic MHD code HARM from axisymmetry to full 3D. The use of an energy conserving algorithm allows the energy dissipated in the course of relativistic accretion to be captured as heat. The inclusion of a simple optically thin cooling function permits explicit control of the simulated disk's geometric thickness as well as a direct calculation of both the amplitude and location of the radiative cooling associated with the accretion stresses. Fully relativistic ray-tracing is used to compute the luminosity received by distant observers. For a disk with aspect ratio H/r ~ 0.1 accreting onto a black hole with spin parameter a/M = 0.9, we find that there is significant dissipation beyond that predicted by the classical Novikov-Thorne model. However, much of it occurs deep in the potential, where photon capture and gravitational redshifting can strongly limit the net photon energy escaping to infinity. In addition, with these parameters and this radiation model, significant thermal and magnetic energy remains with the gas and is accreted by the black hole. In our model, the net luminosi ty reaching infinity is 6% greater than the Novikov-Thorne prediction. If the accreted thermal energy were wholly radiated, the total luminosity of the accretion flow would be ~20% greater than the Novikov-Thorne value.

General Relativistic (GR) Magnetohydrodynamic (MHD) simulations of black hole accretion find significant magnetic stresses near and inside the innermost stable circular orbit (ISCO), suggesting that such flows could radiate in a manner noticeably different from the prediction of the standard model, which assumes that there are no stresses in that region. We provide estimates of how phenomenologically interesting parameters like the ``radiation edge", the innermost ring of the disc from which substantial thermal radiation escapes to infinity, may be altered by stresses near the ISCO. These estimates are based on data from a large number of three-dimensional GRMHD simulations combined with GR ray-tracing. For slowly spinning black holes ($a/M<0.9$), the radiation edge lies well inside where the standard model predicts, particularly when the system is viewed at high inclination. For more rapidly spinning black holes, the contrast is smaller. At fixed total luminosity, the characteristic temperature of the accretion flow increases between a factor of $1.2-2.4$ over that predicted by the standard model, whilst at fixed mass accretion rate, there is a corresponding enhancement of the accretion luminosity which may be anywhere from tens of percent to order unity. When all these considerations are combined, we find that, for fixed black hole mass, luminosity, and inclination angle, our uncertainty in the characteristic temperature of the radiation reaching distant observers due to uncertainty in dissipation profile (around a factor of 3) is {\it greater} than the uncertainty due to a complete lack of knowledge of the black hole's spin (around a factor of 2) and furthermore that spin estimates based on the stress-free inner boundary condition provide an upper limit to $a/M$.

A new code for astrophysical magnetohydrodynamics (MHD) is described. The code has been designed to be easily extensible for use with static and adaptive mesh refinement. It combines higher-order Godunov methods with the constrained transport (CT) technique to enforce the divergence-free constraint on the magnetic field. Discretization is based on cell-centered volume-averages for mass, momentum, and energy, and face-centered area-averages for the magnetic field. Novel features of the algorithm include (1) a consistent framework for computing the time- and edge-averaged electric fields used by CT to evolve the magnetic field from the time- and area-averaged Godunov fluxes, (2) the extension to MHD of spatial reconstruction schemes that involve a dimensionally-split time advance, and (3) the extension to MHD of two different dimensionally-unsplit integration methods. Implementation of the algorithm in both C and Fortran95 is detailed, including strategies for parallelization using domain decomposition. Results from a test suite which includes problems in one-, two-, and three-dimensions for both hydrodynamics and MHD are given, not only to demonstrate the fidelity of the algorithms, but also to enable comparisons to other methods. The source code is freely available for download on the web.

Because the magneto-rotational instability is capable of exponentially amplifying weak preexisting magnetic fields, it might be hoped that the character of the magnetic field in accretion disks is independent of the nature of the seed field. However, the divergence-free nature of magnetic fields in highly conducting fluids ensures that their large-scale topology is preserved, no matter how greatly the field intensity is changed. By performing global two-dimensional and three-dimensional general relativistic magnetohydrodynamic disk simulations with several different topologies for the initial magnetic field, we explore the degree to which the character of the flows around black holes depends on the initial topology. We find that while the qualitative properties of the accretion flow are nearly independent of field topology, jet-launching is very sensitive to it: a sense of vertical field consistent for at least an inner disk inflow time is essential to the support of strong jets.

The inward flow of matter through accretion disks is driven by MHD turbulence. Global general relativistic MHD simulations shed quantitative light on this process, revealing a number of aspects of accretion previously unrecognized. Among them are strong stresses in the marginally stable and plunging regions near the black hole and electromagnetically-dominated conical relativistic jets that can form spontaneously from the accretion flow. The energy release associated with both of these effects can significantly augment the classical energy-release estimates based on purely hydrodynamic models.

We present an exact three-dimensional wave solution to the shearing sheet equations of motion. The existence of this solution argues against transient amplification as a route to turbulence in unmagnetized disks. Moreover, because the solution covers an extensive dynamical range in wavenumber space, it is an excellent test of the dissipative properties of numerical codes.

We use a relativistic ray-tracing code to calculate the light curves observed from a global general relativistic magneto-hydrodynamic simulation of an accretion flow onto a Schwarzschild black hole. We apply three basic emission models to sample different properties of the time-dependent accretion disk. With one of these models, which assumes thermal blackbody emission and free-free absorption, we can predict qualitative features of the high-frequency power spectrum from stellar-mass black holes in the "Thermal Dominant" state. The simulated power spectrum is characterized by a power law of index Gamma ~ 3 and total rms fractional variance of <~ 2% above 10 Hz. For each emission model, we find that the variability amplitude should increase with increasing inclination angle. On the basis of a newly-developed formalism for quantifying the significance of quasi-periodic oscillations (QPOs) in simulation data, we find that these simulations are able to identify any such features with (rms/mean) amplitudes >~ 1 % near the orbital frequency at the inner-most stable orbit. Initial results indicate the existence of transient QPO peaks with frequency ratios of nearly 2:3 at a 99.9% confidence limit, but they are not generic features because at any given time they are seen only from certain observer directions. Additionally, we present detailed analysis of the azimuthal structure of the accretion disk and the evolution of density perturbations in the inner disk. These "hot spot" structures appear to be roughly self-similar over a range of disk radii, with a single characteristic size \delta\phi=25 deg and \delta r/r=0.3, and typical lifetimes T_l ~ 0.3 T_orb.

Simulations of black hole accretion have shown that magnetic stresses are present near and inside the innermost stable circular orbit (ISCO). This finding suggests that such flows may be more luminous than predicted by the standard accretion disk model. Here we apply a prescription for heat dissipation within the simulated accretion flows to estimate their implied radiative efficiency. We assume that dissipation is proportional to the current density squared, and find that the resulting azimuthally-averaged and shell-integrated radial profile is well-matched to the radial heat dissipation profile of the standard disk model for the region outside the ISCO, particularly when it is adjusted to account for additional stress at the ISCO. In contrast to the standard model, however, the dissipation profile derived from the current density continues past the ISCO and through the plunging region. The total predicted dissipation rate is between $\simeq 30%$ and $\simeq 100%$ greater than that predicted by the standard model, depending on the black hole spin. Most of the additional dissipation takes place just outside the ISCO. To predict luminosities, we assume instantaneous radiation and zero optical depth, but allow for photon capture. The net radiative efficiency seen by a distant observer is increased relative to the standard model by $\simeq 25%$--80%, with the largest fractional increase for intermediate black hole spins because the increase in dissipation from enhanced stress that occurs for rapid spin is partially offset by the increased likelihood that the additional photons will be captured by the hole.

We compute a series of three-dimensional general relativistic magnetohydrodynamic simulations of accretion flows in the Kerr metric to investigate the properties of the unbound outflows that result. The overall strength of these outflows increases sharply with increasing black hole rotation rate, but a number of generic features are found in all cases. The mass in the outflow is concentrated in a hollow cone whose opening angle is largely determined by the effective potential for matter orbiting with angular momentum comparable to that of the innermost stable circular orbit. The dominant force accelerating the matter outward comes from the pressure of the accretion disk's corona. The principal element that shapes the outflow is therefore the centrifugal barrier preventing accreting matter from coming close to the rotation axis. Inside the centrifugal barrier, the cone contains very little matter and is dominated by electromagnetic fields that rotate at a rate tied closely to the rotation of the black hole. These fields carry an outward-going Poynting flux whose immediate energy source is the rotating spacetime of the Kerr black hole. When the spin parameter a/M of the black hole exceeds ~0.9, the energy carried to infinity by these outflows can be comparable to the nominal radiative efficiency predicted in the Novikov-Thorne model. Similarly, the expelled angular momentum can be comparable to that accreted by the black hole. Both the inner electromagnetic part and the outer matter part can contribute in significant fashion to the energy and angular momentum of the outflow.

This paper continues the analysis of a set of general relativistic 3D MHD simulations of accreting tori in the Kerr metric with different black hole spins. We focus on bound matter inside the initial pressure maximum, where the time-averaged motion of gas is inward and an accretion disk forms. We use the flows of mass, angular momentum, and energy in order to understand dynamics in this region. The sharp reduction in accretion rate with increasing black hole spin reported in Paper I of this series is explained by a strongly spin-dependent outward flux of angular momentum conveyed electromagnetically; when a/M > 0.9, this flux can be comparable to the inward angular momentum flux carried by the matter. In all cases, there is outward electromagnetic angular momentum flux throughout the flow; in other words, contrary to the assertions of traditional accretion disk theory, there is in general no "stress edge", no surface within which the stress is zero. The retardation of accretion in the inner disk by electromagnetic torques also alters the radial distribution of surface density, an effect that may have consequences for observable properties such as Compton reflection. The net accreted angular momentum is sufficiently depressed by electromagnetic effects that in the most rapidly-spinning black holes mass growth can lead to spindown. Spinning black holes also lose energy by Poynting flux; this rate is also a strongly increasing function of black hole spin, rising to 10% or more of the rest-mass accretion rate at very high spin. As the black hole spins faster, the path of the Poynting flux changes from being predominantly within the accretion disk to predominantly within the funnel outflow.

We have carried out fully relativistic numerical simulations of accretion disks in the Kerr metric. In this paper we focus on the unbound outflows that emerge self-consistently from the accretion flow. These outflows are found in the axial funnel region and consist of two components: a hot, fast, tenuous outflow in the axial funnel proper, and a colder, slower, denser jet along the funnel wall. Although a rotating black hole is not required to produce these unbound outflows, their strength is enhanced by black hole spin. The funnel-wall jet is excluded from the axial funnel due to elevated angular momentum, and is also pressure-confined by a magnetized corona. The tenuous funnel outflow accounts for a significant fraction of the energy transported to large distances in the higher-spin simulations. We compare the outflows observed in our simulations with those seen in other simulations.

Observations are providing increasingly detailed quantitative information about the accretion flows that power such high energy systems as X-ray binaries and active galactic nuclei. Analytic models of such systems must rely on assumptions such as regular flow geometry and a simple, parameterized stress. Global numerical simulations offer a way to investigate the basic physical dynamics of accretion flows without these assumptions. For black hole accretion studies one solves the equations of general relativistic magnetohydrodynamics. Magnetic fields are of fundamental importance to the structure and evolution of accretion disks because magnetic turbulence is the source of the anomalous stress that drives accretion. We have developed a three-dimensional general relativistic magnetohydrodynamic simulation code to evolve time-dependent accretion systems self-consistently. Recent global simulations of black hole accretion disks suggest that the generic structure of the accretion flow is usefully divided into five regimes: the main disk, the inner disk, the corona, the evacuated funnel, and the funnel wall jet. The properties of each of these regions are summarized.

We present a detailed analysis of the magnetic field structure found in general relativistic 3D MHD simulations of accreting tori in the Kerr metric with different black hole spins. Among the properties analyzed are the field strength as a function of position and black hole spin, the shapes of field lines, the degree to which they connect different regions, and their degree of tangling. We investigate prior speculations about the structure of the magnetic fields and discuss how frequently certain configurations are seen in the simulations. We also analyze the distribution of current density, with a view toward identifying possible locations for magnetic energy dissipation.

This is the first in a series of papers that investigate the properties of accretion flows in the Kerr metric through three-dimensional, general relativistic magnetohydrodynamic simulations of tori with a near-Keplerian initial angular velocity profile. We study four models with increasing black hole spin, from a/M=0 to 0.998, for which the structural parameters of the initial tori are maintained nearly constant. The subsequent accretion flows arise self-consistently from stresses and turbulence created by the magnetorotational instability. We investigate the overall evolution and the late-time global structure in the resulting non-radiative accretion flows, including the magnetic fields within the disks, the properties of the flow in the plunging region, and the flux of conserved quantities into the black hole. Independent of black hole spin, the global structure is described in terms of five regions: the main disk body, the coronal envelope, the inner disk, consisting of an inner torus and plunging region, an evacuated axial funnel, and a bi-conical outflow confined to the corona-funnel boundary. We find evidence for lower accretion rates, stronger funnel-wall outflows, and increased stress in the near hole region with increasing black hole spin.

This paper presents an initial survey of the properties of accretion flows in the Kerr metric from three-dimensional, general relativistic magnetohydrodynamic simulations of accretion tori. We consider three fiducial models of tori around rotating, both prograde and retrograde, and nonrotating black holes; these three fiducial models are also contrasted with axisymmetric simulations and a pseudo-Newtonian simulation with equivalent initial conditions to delineate the limitations of these approximations.

The processes of planet formation and migration depend intimately on the interaction between planetesimals and the gaseous disks in which they form. The formation of gaps in the disk can severely limit the mass of the planet and its migration toward the protostar. We investigate the process of gap formation through magnetohydrodynamic simulations in which internal stress arises self-consistently from turbulence generated by the magnetorotational instability. The simulations investigate three different planetary masses and two disk temperatures to bracket the tidal (thermal) and viscous gap opening conditions. The results are in general qualitative agreement with previous simulations of gap formation, but show significant differences. In the presence of MHD turbulence, the gaps produced are shallower and asymmetrically wider than those produced with pure hydrodynamics. The rate of gap formation is also slowed, with accretion occurring across the developing gap. Viscous hydrodynamics does not adequately describe the evolution, however, because planets capable of producing gaps also may be capable of affecting the level MHD turbulence in different regions of the disk.

Chaotic flow is studied in a series of numerical magnetohydrodynamical simulations that use the shearing box formalism. This mimics important features of local accretion disk dynamics. The magnetorotational instability gives rise to flow turbulence, and quantitative chaos parameters, such as the largest Lyapunov exponent, can be measured. Linear growth rates appear in these exponents even when the flow is fully turbulent. The extreme sensitivity to initial conditions associated with chaotic flows has practical implications, the most important of which is that hundreds of orbital times are needed to extract a meaningful average for the stress. If the evolution time in a disk is less than this, the classical $\alpha$ formalism will break down.

This paper describes the development and testing of a general relativistic magnetohydrodynamic (GRMHD) code to study ideal MHD in the fixed background of a Kerr black hole. The code is a direct extension of the hydrodynamic code of Hawley, Smarr, and Wilson, and uses Evans and Hawley constrained transport (CT) to evolve the magnetic fields. Two categories of test cases were undertaken. A one dimensional version of the code (Minkowski metric) was used to verify code performance in the special relativistic limit. The tests include Alfv\'en wave propagation, fast and slow magnetosonic shocks, rarefaction waves, and both relativistic and non-relativistic shock tubes. A series of one- and two-dimensional tests were also carried out in the Kerr metric: magnetized Bondi inflow, a magnetized inflow test due to Gammie, and two-dimensional magnetized constant-$l$ tori that are subject to the magnetorotational instability.

This paper presents results of three-dimensional simulations of global hydrodynamic instabilities in black hole tori, extending earlier work by Hawley to Kerr spacetimes. This study probes a three-dimensional parameter space of torus angular momentum, torus size, and black hole angular momentum. We have observed the growth of the Papaloizou-Pringle instability for a range of torus configurations and the resultant formation of m=1 planets. We have also observed the quenching of this instability in the presence of early accretion flows; however, in one simulation both early accretion and planet formation occurred. Though most of the conclusions reached in Hawley's earlier work on Schwarzschild black holes carry over to Kerr spacetime, the presence of frame dragging in the Kerr geometry adds an element of complexity to the simulations; we have seen especially clear examples of this phenomenon in the accretion flows that arise from retrograde tori.

We review numerical simulations of MHD turbulence. The last decade has witnessed fundamental advances both in the technical capabilities of direct numerical simulation, and in our understanding of key physical processes. Magnetic fields tap directly into the free energy sources in a sufficiently ionized gas. The result is that adverse angular velocity and adverse temperature gradients, not the classical angular momentum and entropy gradients, destabilize laminar and stratified flow. This has profound consequences for astrophysical accretion flows, and has opened the door to a new era of numerical simulation experiments.}

The principles underlying a proposed class of black hole accretion models are examined. The flows are generally referred to as ``convection-dominated,'' and are characterized by inward transport of angular momentum by thermal convection and outward viscous transport, vanishing mass accretion, and vanishing local energy dissipation. We examine the viability of these ideas by explicitly calculating the leading order angular momentum transport of axisymmetric modes in magnetized, differentially rotating, stratified flows. The modes are destabilized by the generalized magnetorotational instability. There are no inward transporting modes at all unless the magnitude of the (imaginary) Brunt frequency is comparable to the epicyclic frequency, a condition requiring substantial levels of dissipation. When inward transporting modes do exist, they appear at long wavelengths, unencumbered by magnetic tension. Moreover, very general thermodynamic principles prohibit the complete recovery of irreversible dissipative energy losses, a central feature of convection-dominated models. Dissipation of the free energy of differential rotation inevitably requires outward angular momentum transport. Our results are in good agreement with global MHD simulations, which find significant levels of outward transport and energy dissipation, whether or not destabilizing entropy gradients are present.

We analyze three-dimensional magnetohydrodynamic (MHD) simulations of a nonradiative accretion flow around a black hole using a pseudo-Newtonian potential. The flow originates from a torus initially centered at 100 gravitational (Schwarzschild) radii. Accretion is driven by turbulent stresses generated self-consistently by the magnetorotational instability. The resulting flow has three well-defined dynamical components: a hot, thick, rotationally-dominated Keplerian disk; a surrounding magnetized corona with vigorous circulation and outflow; and a magnetically-confined jet along the centrifugal funnel wall. Inside of 10 gravitational radii, the disk becomes very hot, more toroidal, and highly intermittent. These results contrast sharply with quasi-spherical, self-similar viscous models. There are no significant dynamical differences between simulations that include resistive heating and those that do not. We conclude by deducing some simple radiative properties of our solutions, and apply the results to the accretion-powered Galactic center source Sgr A*.

What is meant by the "inner edge" of an accretion disk around a black hole depends on the property that defines the edge. We discuss four such definitions using data from recent high-resolution numerical simulations. These are: the "turbulence edge", where flux-freezing becomes more important than turbulence in determining the magnetic field structure; the "stress edge", where plunging matter loses dynamical contact with the outer accretion flow; the "reflection edge", the smallest radius capable of producing significant X-ray reflection features; and the "radiation edge", the innermost place from which significant luminosity emerges. All these edges are dependent on the accretion rate and are non-axisymmetric and time-variable. Although all are generally located in the vicinity of the marginally stable orbit, significant displacements can occur, and data interpretations placing the disk edge precisely at this point can be misleading. If observations are to be used successfully as diagnostics of accretion in strong gravity, the models used to interpret them must take careful account of these distinctions.

New three dimensional magnetohydrodynamic simulations of accretion disk dynamics in a pseudo-Newtonian Paczynski-Wiita potential are presented. These have finer resolution in the inner disk than any previously reported. Finer resolution leads to increased magnetic field strength, greater accretion rate, and greater fluctuations in the accretion rate. One simulation begins with a purely poloidal magnetic field, the other with a purely toroidal field. Compared to the poloidal initial field simulation, a purely toroidal initial field takes longer to reach saturation of the magnetorotational instability and produces less turbulence and weaker magnetic field energies. For both initial field configurations, magnetic stresses continue across the marginally stable orbit; measured in units corresponding to the Shakura-Sunyaev alpha parameter, the stress grows from ~0.1 in the disk body to as much as ~10 deep in the plunging region. Matter passing the inner boundary of the simulation has ~10% greater binding energy and ~10% smaller angular momentum than it did at the marginally stable orbit. Both the mass accretion rate and the integrated stress fluctuate widely on a broad range of timescales.

We present a global magnetohydrodynamic (MHD) three dimensional simulation of a nonradiative accretion flow originating in a pressure supported torus. The evolution is controlled by the magnetorotational instability which produces turbulence. The flow forms a nearly Keplerian disk. The total pressure scale height in this disk is comparable to the vertical size of the initial torus. Gas pressure dominates only near the equator; magnetic pressure is more important in the surrounding atmosphere. A magnetically dominated bound outflow is driven from the disk. The accretion rate through the disk exceeds the final rate into the hole, and a hot torus forms inside 10 r_g. Hot gas, pushed up against the centrifugal barrier and confined by magnetic pressure, is ejected in a narrow, unbound, conical outflow. The dynamics are controlled by magnetic turbulence, not thermal convection, and a hydrodynamic alpha model is inadequate to describe the flow. The limitations of two dimensional MHD simulations are also discussed.

This paper presents a series of global three dimensional accretion disk simulations carried out in the cylindrical limit in which the vertical component of the gravitational field is neglected. The simulations use a cylindrical pseudo-Newtonian potential to model the main dynamical properties of the Schwarzschild metric. The disks are initially constant density with a Keplerian angular momentum distribution and contain a weak toroidal or vertical field. These simulations reaffirm many of the conclusions of previous local simulations. The magnetorotational instability grows rapidly and produces MHD turbulence with a significant Maxwell stress which drives accretion. Tightly-wrapped low-$m$ spiral waves are prominent. In some simulations radial variations in Maxwell stress concentrate gas into rings, creating substantial spatial inhomogeneities. There is a nonzero stress at the marginally stable orbit which produces a small decline in specific angular momentum inside the last stable orbit. Detailed comparisons between simulations are used to examine the effects of computational domain and equation of state. Simulations that begin with vertical fields have greater field amplification and higher ratios of stress to magnetic pressure compared with those beginning with toroidal fields. In contrast to MHD, hydrodynamics alone neither creates nor sustains turbulence.

We present a detailed three dimensional magnetohydrodynamic (MHD) simulation describing the inner region of a disk accreting onto a black hole. To avoid the technical complications of general relativity, the dynamics are treated in Newtonian fashion using the pseudo-Newtonian Pacz\'ynski-Wiita potential. The disk evolves due to angular momentum transport which is produced naturally from MHD turbulence generated by the magnetorotational instability. We find that the resulting stress is continuous across the marginally stable orbit, in contradiction with the widely-held assumption that the stress should go to zero there. As a consequence, the specific angular momentum of the matter accreted into the hole is smaller than the specific angular momentum at the marginally stable orbit. The disk exhibits large fluctuations in almost every quantity, both spatially and temporally. In particular, the ratio of stress to pressure (the local analog of the Shakura-Sunyaev $\alpha$ parameter) exhibits both systematic gradients and large fluctuations; from $\sim 10^{-2}$ in the disk midplane at large radius, it rises to $\sim 10$ both at a few gas density scaleheights above the plane at large radius, and near the midplane well inside the plunging region. Driven in part by large-amplitude waves excited near the marginally stable orbit, both the mass accretion rate and the integrated stress exhibit large fluctuations whose Fourier power spectra are smooth "red" power-laws stretching over several orders of magnitude in timescale.

We present three-dimensional magnetohydrodynamic simulations of the nonlinear evolution of the magnetorotational instability (MRI) with a non-zero Ohmic resistivity. The properties of the saturated state depend on the initial magnetic field configuration. In simulations with an initial uniform vertical field, the MRI is able to support angular momentum transport even for large resistivities through the quasi-periodic generation of axisymmetric radial channel solutions rather than through the maintenance of anisotropic turbulence. Simulations with zero net flux show that the angular momentum transport and the amplitude of magnetic energy after saturation are significantly reduced by finite resistivity, even at levels where the linear modes are only slightly affected. This occurs at magnetic Reynolds numbers expected in low, cool states of dwarf novae, these results suggest that finite resistivity may account for the low and high angular momentum transport rates inferred for these systems.

Global time-dependent simulations provide a means to investigate time-dependent dynamic evolution in accretion disks. This paper seeks to extend previous local simulations by beginning a systematic effort to develop fully global three-dimensional simulations. The nonlinear development of the magnetorotational instability is investigated using a time-explicit finite difference code written in cylindrical coordinates. The equations of ideal magnetohydrodynamics are solved with the assumption of an adiabatic equation of state. Both a Newtonian, and a pseudo-Newtonian potential are used. Two simplifications are also explored: a cylindrical gravitational potential (the ``cylindrical disk''), and axisymmetry. The results from those simulations are compared with fully three dimensional global simulations. The global simulations investigate initial equilibrium pressure supported accretion tori with two different initial field geometries: poloidal fields that are constant along initial equidensity surfaces, and toroidal fields with a constant ratio of gas to magnetic pressure. In both cases the magnetorotational instability rapidly develops, and the torus becomes turbulent. The resulting turbulence transports angular momentum, and the torus develops an angular momentum distribution that is near Keplerian.

The dynamical state of the solar nebula depends critically upon whether or not the gas is magnetically coupled. The presence of a subthermal field will cause laminar flow to break down into turbulence. Magnetic coupling, in turn, depends upon the ionization fraction of the gas. The inner most region of the nebula ($\lesssim 0.1$ AU) is magnetically well-coupled, as is the outermost region ($\gtrsim 10$ AU). The magnetic status of intermediate scales ($\sim 1$ AU) is less certain. It is plausible that there is a zone adjacent to the inner disk in which turbulent heating self-consistently maintains the requisite ionization levels. But the region adjacent to the active outer disk is likely to be magnetically ``dead.'' Hall currents play a significant role in nebular magnetohydrodynamics. Though still occasionally argued in the literature, there is simply no evidence to support the once standard claim that differential rotation in a Keplerian disk is prone to break down into shear turbulence by nonlinear instabilities. There is abundant evidence---numerical, experimental, and analytic---in support of the stabilizing role of Coriolis forces. Hydrodynamical turbulence is almost certainly not a source of enhanced turbulence in the solar nebula, or in any other astrophysical accretion disk.

We employ a variety of numerical simulations in the local shearing box system to investigate in greater depth the local hydrodynamic stability of Keplerian differential rotation. In particular we explore the relationship of Keplerian shear to the nonlinear instabilities known to exist in simple Cartesian shear. The Coriolis force is the source of linear stabilization in differential rotation. We exploit the formal equivalence of constant angular momentum flows and simple Cartesian shear to examine the transition from stability to nonlinear instability. The manifestation of nonlinear instability in shear flows is known to be sensitive to initial perturbation and to the amount of viscosity; marginally (linearly) stable differentially rotating flows exhibit this same sensitivity. Keplerian systems, however, are completely stable; the strength of the stabilizing Coriolis force easily overwhelms any destabilizing nonlinear effects. In fact, nonlinear effects speed the decay of applied turbulence by producing a rapid cascade of energy to high wavenumbers where dissipation occurs. Our conclusions are tested with grid resolution experiments and by comparison with results from a code that employs an alternative numerical algorithm. The properties of hydrodynamic differential rotation are contrasted with magnetohydrodynamic differential rotation. The kinetic stress couples to the vorticity which limits turbulence, while the magnetic stress couples to the shear which promotes turbulence. Thus magnetohydrodynamic turbulence is uniquely capable of acting as a turbulent angular momentum transport mechanism in disks.

We carry out three-dimensional magnetohydrodynamical simulations of the magnetorotational (Balbus-Hawley) instability in weakly-ionized plasmas. We adopt a formulation in which the ions and neutrals each are treated as separate fluids coupled only through a collisional drag term. Ionization and recombination processes are not considered. The linear stability of the ion-neutral system has been previously considered by Blaes and Balbus (1994). Here we extend their results into the nonlinear regime by computing the evolution of Keplerian angular momentum distribution in the local shearing box approximation. We find significant turbulence and angular momentum transport when the collisional frequency is on order 100 times the orbital frequency. At higher collision rates, the two-fluid system studied here behaves much like the fully ionized systems studied previously. At lower collision rates the evolution of the instability is determined primarily by the properties of the ions, with the neutrals acting as a sink for the turbulence. Since in this regime saturation occurs when the magnetic field is superthermal with respect to the ion pressure, we find the amplitude of the magnetic energy and the corresponding angular momentum transport rate is proportional to the ion density. Our calculations show the ions and neutrals are essentially decoupled when the collision frequency is less than 0.01 times the orbital frequency; in this case the ion fluid behaves as in the single fluid simulations and the neutrals remain quiescent. We find that purely toroidal initial magnetic field configurations are unstable to the magnetorotational instability across the range of coupling frequencies.