We present the results of three-dimensional global resistive magnetohydrodynamic (MHD) simulations of black hole accretion flows. General relativistic effects are simulated by using the pseudo-Newtonian potential. Initial state is an equilibrium model of a torus threaded by weak toroidal magnetic fields. As the magnetorotational instability (MRI) grows in the torus, mass accretes to the black hole by losing the angular momentum. We found that in the innermost plunging region, non-axisymmetric accretion flow creates bisymmetric spiral magnetic fields and current sheets. Mass accretion along the spiral channel creates one armed spiral density distribution. Since the accreting matter carries in magnetic fields which subsequently are stretched and amplified due to differential rotation, current density increases inside the channel. Magnetic reconnection taking place in the current sheet produces slow mode shock waves which propagate away from the reconnection site. Magnetic energy release in the innermost plunging region can be the origin of X-ray shots observed in black hole candidates. Numerical simulations reproduced soft X-ray excess preceding the peak of the shots, X-ray hardening at the peak of the shot, and hard X-ray time lags.Comment: 31 pages, 17 figures (reduced resolution). Accepted for publication in ApJ. High resolution figures can be downloaded from http://www.astro.phys.s.chiba-u.ac.jp/~machida/apj2002.htm

We carried out three-dimensional global resistive magnetohydrodynamic simulations of the cooling instability in optically thin hot black hole accretion flows by assuming bremsstrahlung cooling. General relativistic effects are simulated using the pseudo-Newtonian potential. The cooling instability grows when the density of the accretion disk becomes sufficiently large. We found that as the instability grows the accretion flow changes from an optically thin, hot, gas pressure-supported state (low/hard state) to a cooler, magnetically supported, quasi-steady state. During this transition, the magnetic pressure exceeds the gas pressure because the disk shrinks in the vertical direction while almost conserving the toroidal magnetic flux. Since further vertical contraction of the disk is suppressed by magnetic pressure, the cool disk stays in an optically thin, spectrally hard state. The magnetically supported disk exists for a time scale much longer than the thermal time scale, and comparable to the accretion time scale. We examined the stability of the magnetically supported disk analytically, assuming that the toroidal magnetic flux is conserved, and found it to be thermally and secularly stable. Our findings may explain why black hole candidates stay in luminous, hard state even when their luminosity exceeds the threshold for the onset of the cooling instability.

The central few hundred parsecs of the Milky Way host a massive black hole and exhibit very violent gas motion and high temperatures in molecular gas. The origin of these properties has been a mystery for the past four decades. Wide-field imaging of the ¹² CO (rotational quantum number J = 1 to 0) 2.6-millimeter spectrum has revealed huge loops of dense molecular gas with strong velocity dispersions in the galactic center. We present a magnetic flotation model to explain that the formation of the loops is due to magnetic buoyancy caused by the Parker instability. The model has the potential to offer a coherent explanation for the origin of the violent motion and extensive heating of the molecular gas in the galactic center.

Two-dimensional (2D) and three-dimensional (3D) hydrodynamical simulations have been intensively performed recently to clarify the accretion-flow structure in the low-radiation-efficiency limit. However, the results depend critically on the parameterized magnitude of the viscosity, which, in principle, should be determined self-consistently by MHD simulations. We analyzed the structure of 3D MHD accretion flows initially threaded by weak toroidal magnetic fields, and found for the first time large-scale convective motions dominating near to the black hole. Radial profiles of each physical quantity include: the density, $ \rho \propto r^{-0.5}$; radial velocity, $ v_r\propto r^{-1.5}$; temperature, $ T \propto r^{-1.0}$; and field strength, $ B^2 \propto r^{-1.5}$. Although the flow structure, itself, appears to be similar to those obtained by hydrodynamic simulations, the observational appearance is distinct. Unlike non-magnetic models, in which radiation is dominant at the outermost convective zones because of outward energy flow by convection, substantial accretion energy can be released in the vicinity of a black hole in MHD flow via magnetic reconnection. Such reconnection leads to sporadic flare events, thus producing variability in out-going radiation, as is commonly observed in objects with black-hole accretion.

We carried out global three-dimensional magnetohydrodynamic simulations of dynamo activities in galactic gaseous disks without assuming equatorial symmetry. Numerical results indicate the growth of azimuthal magnetic fields non-symmetric to the equatorial plane. As the magnetorotational instability (MRI) grows, the mean strength of magnetic fields is amplified until the magnetic pressure becomes as large as 10% of the gas pressure. When the local plasma β (=p gas/p mag) becomes less than 5 near the disk surface, magnetic flux escapes from the disk by the Parker instability within one rotation period of the disk. The buoyant escape of coherent magnetic fields drives dynamo activities by generating disk magnetic fields with opposite polarity to satisfy the magnetic flux conservation. The flotation of the azimuthal magnetic flux from the disk and the subsequent amplification of disk magnetic field by the MRI drive quasi-periodic reversal of azimuthal magnetic fields on a timescale of 10 rotation periods. Since the rotation speed decreases with radius, the interval between the reversal of azimuthal magnetic fields increases with radius. The rotation measure computed from the numerical results shows symmetry corresponding to a dipole field.