The role of radiation in general relativistic magnetohydrodynamic (GRMHD) accretion simulations is discussed through axisymmetric simulations of the evolution of an initial torus seeded with a weak magnetic field. The paper compares and contrasts the rate of accretion onto a Kerr black hole and mass flux out out of the initial torus at large radii in the GRMHD code of De Villiers and Hawley and a newly developed radiative GRMHD code. This rGRMHD code currently operates in the diffusion approximation, restricting the study of radiative effects to the bound portion of the accretion disk/jet system. However, these preliminary findings suggest that radiative effects do play a potentially significant role in regulating the accretion flow.

This paper provides an alternate derivation of the equations used in the GRMHD code of De Villiers and Hawley using Stokes Theorem. This derivation places the published form of the discretized equations of GRMHD in a broader context, and suggests possible future directions for numerical work.

This paper presents a survey of a set of simulations of accretion disks orbiting a rapidly rotating Kerr black hole and embedded in a large-scale initial magnetic field. Each simulation uses a common state for the initial torus, including an MRI seed-field consisting of poloidal loops along isodensity contours, and differ only in the strength of the initial large-scale magnetic field in relation to these poloidal loops. Simulations with a weak large-scale initial field differ little from simulations that evolve from only the poloidal field. Simulations where the large-scale initial field distorts or completely overwhelms the poloidal loops show more extensive regions of turbulence, due to the action of the MRI on the large-scale field. However, the overall structure of the late-time state of the simulations remains qualitatively unchanged, that is each simulation sees the emergence of a turbulent accretion disk, axial jets, a funnel-wall outflow and an extensive corona.

We have carried out 2D and 3D general relativistic magnetohydrodynamic simulations of jets launched self-consistently from accretion disks orbiting Kerr black holes and applied the results to the inner engine of the collapsar model of gamma-ray bursts. The accretion flow launches energetic jets in the axial funnel region of the disk/jet system, as well as a substantial coronal wind. The jets feature knot-like structures of extremely hot, ultra-relativistic gas; the gas in these knots begins at moderate velocities near the inner engine, and is accelerated to ultra-relativistic velocities (Lorentz factors of 50, and higher) by the Lorentz force in the axial funnel. The increase in jet velocity takes place in an acceleration zone extending to at least a few hundred gravitational radii from the inner engine. The overall energetics of the jets are strongly spin-dependent, with high-spin black holes producing the highest energy and mass fluxes. In addition, with high-spin black holes, the ultra-relativistic outflow is cylindrically collimated within a few hundred gravitational radii of the black hole, whereas in the zero- spin case the jet retains a constant opening angle of approximately 16 degrees. The simulations also show that the coronal wind, though considerably slower and colder than the jets, also carries a significant amount of mass and energy. When simulation data is scaled to the physical dimensions of a collapsar the jets operate for a period ranging from 0.1 to 1.4 seconds, until the accretion disk is depleted, delivering $10^{48}$ to $10^{49}$ erg.

In this paper and a companion work, we report on the first global numerical simulations of self-gravitating magnetized tori, subject in particular to the influence of the magnetorotational instability (MRI). In this work, paper I, we restrict our calculations to the study of the axisymmetric evolution of such tori. Our goals are twofold: (1) to investigate how self-gravity influences the global structure and evolution of the disks; and (2) to determine whether turbulent density inhomogeneities can be enhanced by self-gravity in this regime. As in non self-gravitating models, the linear growth of the MRI is followed by a turbulent phase during which angular momentum is transported outward. As a result, self-gravitating tori quickly develop a dual structure composed of an inner thin Keplerian disk fed by a thicker self-gravitating disk, whose rotation profile is close to a Mestel disk. Our results show that the effects of self-gravity enhance density fluctuations much less than they smooth the disk, and giving it more coherence. We discuss the expected changes that will occur in 3D simulations, the results of which are presented in a companion paper.

We present 3D magnetohydrodynamic (MHD) numerical simulations of the evolution of self--gravitating and weakly magnetized disks with an adiabatic equation of state. Such disks are subject to the development of both the magnetorotational and gravitational instabilities, which transport angular momentum outward. As in previous studies, our hydrodynamical simulations show the growth of strong m=2 spiral structure. This spiral disturbance drives matter toward the central object and disappears when the Toomre parameter Q has increased well above unity. When a weak magnetic field is present as well, the magnetorotational instability grows and leads to turbulence. In that case, the strength of the gravitational stress tensor is lowered by a factor of about~2 compared to the hydrodynamical run and oscillates periodically, reaching very small values at its minimum. We attribute this behavior to the presence of a second spiral mode with higher pattern speed than the one which dominates in the hydrodynamical simulations. It is apparently excited by the high frequency motions associated with MHD turbulence. The nonlinear coupling between these two spiral modes gives rise to a stress tensor that oscillates with a frequency which is a combination of the frequencies of each of the modes. This interaction between MHD turbulence and gravitational instabilities therefore results in a smaller mass accretion rate onto the central object.

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 global 2D and 3D simulations of self-gravitating magnetized tori. We used the 2D calculations to demonstrate that the properties of the MRI are not affected by the presence of self-gravity: MHD turbulence and enhanced angular momentum transport follow the linear growth of the instability. In 3D, we have studied the interaction between an $m=2$ gravitational instability and MHD turbulence. We found its strength to be significantly decreased by the presence of the latter, showing that both instabilities strongly interact in their non-linear phases. We discuss the consequences of these results.

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.

We present the first global simulations of self-gravitating magnetized tori. The simulations are performed with Zeus-2D and GLOBAL. We find the magnetorotational instability (MRI) to behave similarly in a self-gravitating environment as in previous simulations of non self-gravitating systems: enhancement of turbulent angular momentum transport follows the linear phase. The torus quickly develops a two component structure composed of an inner thick disk in Keplerian rotation and an outer massive disk. We compare this result with zero mass global simulations in 2D, and also present preliminary results of 3D simulations.

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.

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.