Using local three dimensional radiation hydrodynamics simulations, the nonlinear outcome of gravitational instability in an irradiated protoplanetary disc is investigated in a parameter space of the surface density $\Sigma$ and the radius $r$. Starting from laminar flow, axisymmetric self-gravitating density waves grow first. Their self-gravitating degree becomes larger when $\Sigma$ is larger or the cooling time is shorter at larger radii. The density waves eventually collapse owing to non-axisymmetric instability, which results in either fragmentation or gravito-turbulence after a transient phase. The boundaries between the two are found at $r \sim 75$ AU as well as at the $\Sigma$ that corresponds to the initial Toomre's parameter of $\sim 0.2$. The former boundary corresponds to the radius where the cooling time becomes short, approximating unity. Even when gravito-turbulence is established around the boundary radius, such a short cooling time inevitably makes the fluctuation of $\Sigma$ large enough to trigger fragmentation. On the other hand, when $\Sigma$ is beyond the latter boundary (i.e. the initial Toomre's parameter is less than $\sim 0.2$), the initial laminar flow is so unstable against self-gravity that it evolves into fragmentation regardless of the radius or, equivalently, the cooling time. Runaway collapse follows fragmentation when the mass concentration at the centre of a bound object is high enough that the temperature exceeds the H$_2$ dissociation temperature.

Warm Comptonization models for the soft X-ray excess in AGN do not self-consistently explain the relationship between the Comptonizing medium and the underlying accretion disc. Because of this, they cannot directly connect the fitted Comptonization temperatures and optical depths to accretion disc parameters. Since bulk velocities exceed thermal velocities in highly radiation pressure dominated discs, in these systems bulk Comptonization by turbulence may provide a physical basis in the disc itself for warm Comptonization models. We model the dependence of bulk Comptonization on fundamental accretion disc parameters, such as mass, luminosity, radius, spin, inner boundary condition, and $\alpha$. In addition to constraining warm Comptonization models, our model can help distinguish contributions from bulk Comptonization to the soft X-ray excess from those due to other physical mechanisms, such as absorption and reflection. By linking the time variability of bulk Comptonization to fluctuations in the disc vertical structure due to magnetorotational instability (MRI) turbulence, our results show that observations of the soft X-ray excess can be used to study disc turbulence in the radiation pressure dominated regime. Because our model connects bulk Comptonization to one dimensional vertical structure temperature profiles in a physically intuitive way, it will be useful for understanding this effect in future simulations run in new regimes.

We present the results of local, vertically stratified, radiation magnetohydrodynamic shearing box simulations of magnetorotational instability (MRI) turbulence for a (hydrogen poor) composition applicable to accretion disks in AM CVn type systems. Many of these accreting white dwarf systems are helium analogues of dwarf novae (DNe). We utilize frequency-integrated opacity and equation of state tables appropriate for this regime to accurately portray the relevant thermodynamics. We find bistability of thermal equilibria in the effective temperature, surface mass density plane typically associated with disk instabilities. Along this equilibrium curve (i.e. the S-curve) we find that the stress to thermal pressure ratio $\alpha$ varied with peak values of $\sim 0.15$ near the tip of the upper branch. Similar to DNe, we found enhancement of $\alpha$ near the tip of the upper branch caused by convection; this increase in $\alpha$ occurred despite our choice of zero net vertical magnetic flux. Two notable differences we find between DN and AM CVn accretion disk simulations are that AM CVn disks are capable of exhibiting persistent convection in outburst, and ideal MHD is valid throughout quiescence for AM CVns. In contrast, DNe simulations only show intermittent convection, and non-ideal MHD effects are likely important in quiescence. By combining our previous work with these new results, we also find that convective enhancement of the MRI is anticorrelated with mean molecular weight.

Using radiation hydrodynamics simulations in a local stratified shearing box with realistic equations of state and opacities, we explored the outcome of self-gravity at 50 AU in a protoplanetary disc irradiated by the central star. We found that gravito-turbulence is sustained for a finite range of the surface density, from $\sim 80$ to $\sim$ 250 gcm$^{-2}$. The disk is laminar below the range while fragments above it. In the range of gravito-turbulence, the Toomre parameter decreases monotonically from $\sim 1$ to $\sim 0.7$ as the surface density increases while an effective cooling time is almost constant at $\sim 4$ in terms of the inverse of the orbital frequency. The turbulent motions are supersonic at all heights, which dissipates through both shock waves and compressional heating. The compressional motions, occurring near the midplane, create upward flows, which not only contribute to supporting the disc but also to transporting the dissipated energy to the disc surfaces. The irradiation does not affect much the gravito-turbulence near the midplane unless the grazing angle is larger than 0.32. We also show that a simple cooling function with a constant cooling time does not approximate the realistic cooling.

Vertically stratified shearing box simulations of magnetorotational turbulence commonly exhibit a so-called butterfly diagram of quasi-periodic azimuthal field reversals. However, in the presence of hydrodynamic convection, field reversals no longer occur. Instead, the azimuthal field strength fluctuates quasi-periodically while maintaining the same polarity, which can either be symmetric or antisymmetric about the disc midplane. Using data from the simulations of Hirose et al. (2014), we demonstrate that the lack of field reversals in the presence of convection is due to hydrodynamic mixing of magnetic field from the more strongly magnetized upper layers into the midplane, which then annihilate field reversals that are starting there. Our convective simulations differ in several respects from those reported in previous work by others, in which stronger magnetization likely plays a more important role than convection.

Bulk velocities exceed thermal velocities for sufficiently radiation pressure dominated accretion flows. We model the contribution of bulk Comptonization to the soft X-ray excess in AGN. Bulk Comptonization is due to both turbulence and the background shear. We calculate spectra both taking into account and not taking into account bulk velocities using scaled data from radiation magnetohydrodynamic (MHD) shearing box simulations. We characterize our results with temperatures and optical depths to make contact with other warm Comptonization models of the soft excess. We chose our fiducial mass, $M = 2 \times 10^6 M_{\odot}$, and accretion rate, $L/L_{\rm Edd} = 2.5$, to correspond to those fit to the super-Eddington narrow line Seyfert 1 (NLS1) RE1034+396. The temperatures, optical depths, and Compton $y$ parameters we find broadly agree with those fit to RE1034+396. The effect of bulk Comptonization is to shift the Wien tail to higher energy and lower the gas temperature, broadening the spectrum. Observations of the soft excess in NLS1s can constrain the properties of disc turbulence if the bulk Comptonization contribution can be separated out from contributions from other physical effects, such as reflection and absorption.

The phenomenological Disc Instability Model has been successful in reproducing the observed light curves of dwarf nova outbursts by invoking an enhanced Shakura-Sunyaev $\alpha$ parameter $\sim0.1-0.2$ in outburst compared to a low value $\sim0.01$ in quiescence. Recent thermodynamically consistent simulations of magnetorotational (MRI) turbulence with appropriate opacities and equation of state for dwarf nova accretion discs have found that thermal convection enhances $\alpha$ in discs in outburst, but only near the hydrogen ionization transition. At higher temperatures, convection no longer exists and $\alpha$ returns to the low value comparable to that in quiescence. In order to check whether this enhancement near the hydrogen ionization transition is sufficient to reproduce observed light curves, we incorporate this MRI-based variation in $\alpha$ into the Disc Instability Model, as well as simulation-based models of turbulent dissipation and convective transport. These MRI-based models can successfully reproduce observed outburst and quiescence durations, as well as outburst amplitudes, albeit with different parameters from the standard Disc Instability Models. The MRI-based model lightcurves exhibit reflares in the decay from outburst, which are not generally observed in dwarf novae. However, we highlight the problematic aspects of the quiescence physics in the Disc Instability Model and MRI simulations that are responsible for this behavior.

A large set of numerical simulations of magnetohydrodynamic (MHD) turbulence induced by the magnetorotational instability (MRI) is presented. Revisiting the previous survey conducted by Sano et al. (2004), we investigate the gas pressure dependence of the saturation level. In ideal MHD simulations, the gas pressure dependence is found to be very sensitive to the choice of a numerical scheme. This is because the numerical magnetic Prandtl number varies according to the scheme as well as the pressure, which considerably affects the results. The saturation level is more sensitive to the numerical magnetic Prandtl number than the pressure. In MHD simulations with explicit viscosity and resistivity, the saturation level increases with the physical magnetic Prandtl number, and it is almost independent of the gas pressure when the magnetic Prandtl number is constant. This is indicative of the incompressible turbulence saturated by the secondary tearing instability.

Using radiation magnetohydrodynamics simulations with realistic opacities and equation of state, and zero net magnetic flux, we have explored thermodynamics in the inner part of protoplanetary discs where magnetic turbulence is expected. The thermal equilibrium curve consists of the upper, lower, and middle branches. The upper (lower) branch corresponds to hot (cool) and optically very (moderately) thick discs, respectively, while the middle branch is characterized by convective energy transport near the midplane. Convection is also the major energy transport process near the low surface density end of the upper branch. There, convective motion is fast with Mach numbers reaching $\gtrsim 0.01$, and enhances both magnetic turbulence and cooling, raising the ratio of vertically-integrated shear stress to vertically-integrated pressure by a factor of several. This convectively enhanced ratio seems a robust feature in accretion discs having an ionization transition. We have also examined causes of the S-shaped thermal equilibrium curve, as well as the thermal stability of the equilibrium solutions. Finally, we compared our results with the disc instability models used to explain FU Ori outbursts. Although the thermal equilibrium curve in our results also exhibits bistability, the surface density contrast across the bistability is an order of magnitude smaller, and the stress-to-pressure ratios in both upper and lower branches are two orders of magnitude greater, than those favored in the disc instability models. It therefore appears likely that FU Ori outbursts are not due solely to a thermal-viscous limit cycle resulting from accretion driven by local magnetic turbulence.

We present the results of local, vertically stratified, radiation MHD shearing box simulations of MRI turbulence appropriate for the hydrogen ionizing regime of dwarf nova and soft X-ray transient outbursts. We incorporate the frequency-integrated opacities and equation of state for this regime, but neglect non-ideal MHD effects and surface irradiation, and do not impose net vertical magnetic flux. We find two stable thermal equilibrium tracks in the effective temperature versus surface mass density plane, in qualitative agreement with the S-curve picture of the standard disk instability model. We find that the large opacity at temperatures near $10^4$K, a corollary of the hydrogen ionization transition, triggers strong, intermittent thermal convection on the upper stable branch. This convection strengthens the magnetic turbulent dynamo and greatly enhances the time-averaged value of the stress to thermal pressure ratio $\alpha$, possibly by generating vertical magnetic field that may seed the axisymmetric magnetorotational instability, and by increasing cooling so that the pressure does not rise in proportion to the turbulent dissipation. These enhanced stress to pressure ratios may alleviate the order of magnitude discrepancy between the $\alpha$-values observationally inferred in the outburst state and those that have been measured from previous local numerical simulations of magnetorotational turbulence that lack net vertical magnetic flux.

Young stars with masses 2-8 Suns, called the Herbig Ae and Be stars, often show a near-infrared excess too large to explain with a hydrostatically-supported circumstellar disk of gas and dust. At the same time the accretion flow carrying the circumstellar gas to the star is thought to be driven by magneto-rotational turbulence, which according to numerical MHD modeling yields an extended low-density atmosphere supported by the magnetic fields. We demonstrate that the base of the atmosphere can be optically-thick to the starlight and that the parts lying near 1 AU are tall enough to double the fraction of the stellar luminosity reprocessed into the near-infrared. We generate synthetic spectral energy distributions (SEDs) using Monte Carlo radiative transfer calculations with opacities for sub-micron silicate and carbonaceous grains. The synthetic SEDs closely follow the median Herbig SED constructed recently by Mulders and Dominik, and in particular match the large near-infrared flux, provided the grains have a mass fraction close to interstellar near the disk's inner rim.

Turbulence driven by magnetorotational instability (MRI) affects planetesimal formation by inducing diffusion and collisional fragmentation of dust particles. We examine conditions preferred for planetesimal formation in MRI-inactive "dead zones" using an analytic dead-zone model based on our recent resistive MHD simulations. We argue that successful planetesimal formation requires not only a sufficiently large dead zone (which can be produced by tiny dust grains) but also a sufficiently small net vertical magnetic flux (NVF). Although often ignored, the latter condition is indeed important since the NVF strength determines the saturation level of turbulence in MRI-active layers. We show that direct collisional formation of icy planetesimal across the fragmentation barrier is possible when the NVF strength is lower than 10 mG (for the minimum-mass solar nebula model). Formation of rocky planetesimals via the secular gravitational instability is also possible within a similar range of the NVF strength. Our results indicate that the fate of planet formation largely depends on how the NVF is radially transported in the initial disk formation and subsequent disk accretion processes.

Turbulence driven by magnetorotational instability (MRI) crucially affects the evolution of solid bodies in protoplanetary disks. On the other hand, small dust particles stabilize MRI by capturing ionized gas particles needed for the coupling of the gas and magnetic fields. To provide an empirical basis for modeling the coevolution of dust and MRI, we perform three-dimensional, ohmic-resistive MHD simulations of a vertically stratified shearing box with an MRI-inactive "dead zone" of various sizes and with a net vertical magnetic flux of various strengths. We find that the vertical structure of turbulence is well characterized by the vertical magnetic flux and three critical heights derived from the linear analysis of MRI in a stratified disk. In particular, the turbulent structure depends on the resistivity profile only through the critical heights and is insensitive to the details of the resistivity profile. We discover scaling relations between the amplitudes of various turbulent quantities (velocity dispersion, density fluctuation, vertical diffusion coefficient, and outflow mass flux) and vertically integrated accretion stresses. We also obtain empirical formulae for the integrated accretion stresses as a function of the vertical magnetic flux and the critical heights. These empirical relations allow to predict the vertical turbulent structure of a protoplanetary disk for a given strength of the magnetic flux and a given resistivity profile.

We examine heating and cooling in protostellar disks using 3-D radiation-MHD calculations of a patch of the Solar nebula at 1 AU, employing the shearing-box and flux-limited radiation diffusion approximations. The disk atmosphere is ionized by stellar X-rays, well-coupled to magnetic fields, and sustains a turbulent accretion flow driven by magneto-rotational instability, while the interior is resistive and magnetically dead. The turbulent layers heat by absorbing the light from the central star and by dissipating the magnetic fields. They are optically-thin to their own radiation and cool inefficiently. The optically-thick interior in contrast is heated only weakly, by re-emission from the atmosphere. The interior is colder than a classical viscous model, and isothermal. The magnetic fields support an extended atmosphere that absorbs the starlight 1.5 times higher than the hydrostatic viscous model. The disk thickness thus measures not the internal temperature, but the magnetic field strength. Fluctuations in the fields move the starlight-absorbing surface up and down. The height ranges between 13% and 24% of the radius over timescales of several orbits, with implications for infrared variability. The fields are buoyant, so the accretion heating occurs higher in the atmosphere than the stresses. The heating is localized around current sheets, caused by magneto-rotational instability at lower elevations and by Parker instability at higher elevations. Gas in the sheets is heated above the stellar irradiation temperature, even though accretion is much less than irradiation power when volume-averaged. The hot optically-thin current sheets might be detectable through their line emission.

Standard models of radiation supported accretion disks generally assume that diffusive radiation flux is solely responsible for vertical heat transport. This requires that heat must be generated at a critical rate per unit volume if the disk is to be in hydrostatic and thermal equilibrium. This raises the question of how heat is generated and how energy is transported in MHD turbulence. By analysis of a number of radiation/MHD stratified shearing-box simulations, we show that the divergence of the diffusive radiation flux is indeed capped at the critical rate, but deep inside the disk, substantial vertical energy flux is also carried by advection of radiation. Work done by radiation pressure is a significant part of the energy budget, and much of this work is dissipated later through damping by radiative diffusion. We show how this damping can be measured in the simulations, and identify its physical origins. Radiative damping accounts for as much as tens of percent of the total dissipation, and is the only realistic physical mechanism for dissipation of turbulence that can actually be resolved in numerical simulations of accretion disks. Buoyancy associated with dynamo-driven, highly magnetized, nearly-isobaric nonlinear slow magnetosonic fluctuations is responsible for the radiation advection flux, and also explains the persistent periodic magnetic upwelling seen at all values of the radiation to gas pressure ratio. The intimate connection between radiation advection and magnetic buoyancy is the first example we know of in astrophysics in which a dynamo has direct impact on the global energetics of a system.

We study the properties of the turbulence driven by the magnetorotational instability (MRI) in a stratified shearing box with outflow boundary conditions and an equation of state determined by self-consistent dissipation and radiation losses. A series of simulations with increasing resolution are performed within a fixed computational box. We achieve numerical convergence with respect to radial and azimuthal resolution. As vertical resolution is improved, the ratio of stress to pressure increases only slowly, but the absolute levels of both the stress and the pressure increase noticeably. The vertical correlation length of the magnetic field within the core of the disk is highly time-dependent, but averaged over time, it is $\simeq 3$ times larger than found in previous unstratifed simulations. This correlation length decreases slowly as vertical resolution increases. We suggest that this greater degree of correlation in stratified than unstratified disks, even near the midplane, is due to field buoyancy. We further show that the undulatory Parker instability drives the magnetic field upwelling at several scaleheights from the midplane that is characteristic of stratified MHD-turbulent disks.

This article reports on three-dimensional (3-D) MHD simulations of non-rotating and rapidly rotating black holes and the adjacent black hole accretion disk magnetospheres. A particular emphasis is placed on the vertical magnetic flux that is advected inward from large radii and threads the equatorial plane near the event horizon. In both cases of non-rotating and rotating black holes, the existence of a significant vertical magnetic field in this region is like a switch that creates powerful jets. There are many similarities in the vertical flux dynamics in these two cases in spite of the tremendous enhancement of azimuthal twisting of the field lines and enhancement of the jet power because of an "ergospheric disk" in the Kerr metric. A 3-D approach is essential because two-dimensional axisymmetric flows are incapable of revealing the nature of vertical flux near a black hole. Poloidal field lines from the ergospheric accretion region have been visualized in 3-D and much of the article is devoted to a formal classification of the different manifestations of vertical flux in the Kerr case.

We present the results of a series of radiation-MHD simulations of a local patch of an accretion disk, with fixed vertical gravity profile but with different surface mass densities and a broad range of radiation to gas pressure ratios. Each simulation achieves a thermal equilibrium that lasts for many cooling times. After averaging over times long compared to a cooling time, we find that the vertically integrated stress is approximately proportional to the vertically-averaged total thermal (gas plus radiation) pressure. We map out--for the first time on the basis of explicit physics--the thermal equilibrium relation between stress and surface density: the stress decreases (increases) with increasing surface mass density when the simulation is radiation (gas) pressure dominated. The dependence of stress on surface mass density in the radiation pressure dominated regime suggests the possibility of a Lightman-Eardley inflow instability, but global simulations or shearing box simulations with much wider radial boxes will be necessary to confirm this and determine its nonlinear behavior.

We present the results of one and three-dimensional radiative transfer calculations of polarized spectra emerging from snapshots of radiation magnetohydrodynamical simulations of the local vertical structure of black hole accretion disks. The simulations cover a wide range of physical regimes relevant for the high/soft state of black hole X-ray binaries. We constrain the uncertainties in theoretical spectral color correction factors due to the presence of magnetic support of the disk surface layers and strong density inhomogeneities. For the radiation dominated simulation, magnetic support increases the color correction factor by about ten percent, but this is largely compensated by a ten percent softening due to inhomogeneities. We also compute the effects of inhomogeneities and Faraday rotation on the resulting polarization. Magnetic fields in the simulations are just strong enough to produce significant Faraday depolarization near the spectral peak of the radiation field. X-ray polarimetry may therefore be a valuable diagnostic of accretion disk magnetic fields, being able to directly test simulations of magnetorotational turbulence.

When the accretion rate is more than a small fraction of Eddington, the inner regions of accretion disks around black holes are expected to be radiation-dominated. However, in the alpha-model, these regions are also expected to be thermally unstable. In this paper, we report two 3-d radiation MHD simulations of a vertically-stratified shearing box in which the ratio of radiation to gas pressure is ~ 10, and yet no thermal runaway occurs over a timespan ~ 40 cooling times. Where the time-averaged dissipation rate is greater than the critical dissipation rate that creates hydrostatic equilibrium by diffusive radiation flux, the time-averaged radiation flux is held to the critical value, with the excess dissipated energy transported by radiative advection. Although the stress and total pressure are well-correlated as predicted by the alpha-model, we show that stress fluctuations precede pressure fluctuations, contrary to the usual supposition that the pressure controls the saturation level of the magnetic energy. This fact explains the thermal stability. Using a simple toy-model, we show that independently-generated magnetic fluctuations can drive radiation pressure fluctuations, creating a correlation between the two while maintaining thermal stability.

We have employed a 3-d energy-conserving radiation MHD code to simulate the vertical structure and thermodynamics of a shearing box whose parameters were chosen so that the radiation and gas pressures would be comparable. The upper layers of this disk segment are magnetically-dominated, creating conditions appropriate for both photon bubble and Parker instabilities. We find little evidence for photon bubbles, even though the simulation has enough spatial resolution to see them and their predicted growth rates are high. On the other hand, there is strong evidence for Parker instabilities, and they appear to dominate the evolution of the magnetically supported surface layers. The disk photosphere is complex, with large density inhomogeneities at both the scattering and effective (thermalization) photospheres of the evolving horizontally-averaged structure. Both the dominant magnetic support and the inhomogeneities are likely to have strong effects on the spectrum and polarization of thermal photons emerging from the disk atmosphere. The inhomogeneities are also large enough to affect models of reflection spectra from the atmospheres of accretion disks.

We explore the thermodynamic and global structural properties of a local patch of an accretion disk whose parameters were chosen so that radiation pressure and gas pressure would be comparable in magnitude. Heating, radiative transport, and cooling are computed self-consistently with the structure by solving the equations of radiation MHD in the shearing-box approximation. Using a fully 3-d and energy-conserving code, we compute the structure and energy balance of this disk segment over a span of more than forty cooling times. As is also true when gas pressure dominates, the disk's upper atmosphere is magnetically-supported. However, unlike the gas-dominated case, no steady-state is reached; instead, the total (i.e., radiation plus gas) energy content fluctuates by factors of 3--4 over timescales of several tens of orbits, with no secular trend. Because the radiation pressure varies much more than the gas pressure, the ratio of radiation pressure to gas pressure varies over the approximate range 0.5--2. The volume-integrated dissipation rate generally increases with increasing total energy, but the mean trend is somewhat slower than linear, and the instantaneous dissipation rate is often a factor of two larger or smaller than the mean for that total energy level. Locally, the dissipation rate per unit volume scales approximately in proportion to the current density; the time-average dissipation rate per unit mass is proportional to m^{-1/2}, where m is the horizontally-averaged mass column density to the nearer of the top or bottom surface. As in our earlier study of a gas-dominated shearing-box, we find that energy transport is completely dominated by radiative diffusion, with Poynting flux carrying less than 1% of the energy lost from the box.

Stellar atmosphere models of ionized accretion disks have generally neglected the contribution of magnetic fields to the vertical hydrostatic support, although magnetic fields are widely believed to play a critical role in the transport of angular momentum. Simulations of magnetorotational turbulence in a vertically stratified shearing box geometry show that magnetic pressure support can be dominant in the upper layers of the disk. We present calculations of accretion disk spectra that include this magnetic pressure support, as well as a vertical dissipation profile based on simulation. Magnetic pressure support generically produces a more vertically extended disk atmosphere with a larger density scale height. This acts to harden the spectrum compared to models that neglect magnetic pressure support. We estimate the significance of this effect on disk-integrated spectra by calculating an illustrative disk model for a stellar mass black hole, assuming that similar magnetic pressure support exists at all radii.

(shortened) We calculate the vertical structure of a local patch of an accretion disk in which heating by dissipation of MRI-driven MHD turbulence is balanced by radiative cooling. Heating, radiative transport, and cooling are computed self-consistently with the structure by solving the equations of radiation MHD in the shearing-box approximation. Using a fully 3-d and energy-conserving code, we compute the structure of this disk segment over a span of more than five cooling times. After a brief relaxation period, a statistically steady-state develops. Measuring height above the midplane in units of the scale-height H predicted by a Shakura-Sunyaev model, we find that the disk atmosphere stretches upward, with the photosphere rising to about 7H, in contrast to the approximately 3H predicted by conventional analytic models. This more extended structure, as well as fluctuations in the height of the photosphere, may lead to departures from Planckian form in the emergent spectra. Dissipation is distributed across the region within roughly 3H of the midplane, but is very weak at greater altitudes. Because fluctuations in the dissipation are particularly strong away from the midplane, the emergent radiation flux can track dissipation fluctuations with a lag that is only 0.1--0.2 times the mean cooling time of the disk. Long timescale asymmetries in the dissipation distribution can also cause significant asymmetry in the flux emerging from the top and bottom surfaces of the disk. Radiative diffusion dominates Poynting flux in the vertical energy flow throughout the disk.

WIDGET is a robotic telescope for monitoring the HETE-2 field-of-view to detect Gamma-ray Burst optical flashes or possible optical precursors. The system has 62degx62deg wide field-of-view which covers about 80% of HETE-2 one with a 2kx2k Apogee U10 CCD camera and a Canon EF 24mm f/1.4 wide-angle lens without a bandpass filter. WIDGET has been in operation since June 2004 at Akeno observing site where is about 200 km apart from Tokyo. Typical limiting magnitude with S/N=3 at the site is V=10mag for 5 seconds exposure and V=11mag for 30 seconds exposure. We had already six coincident observations with HETE-2 position alerts. It was, however, cloudy for all cases due to rainy season in Japan. Expected number of coincident observations under clear sky is about 5 events per year. We will extend the system in early 2005 for Swift era to monitor optical transients in wider field-of-view, multi-color or polarization modes.

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.

The Kelvin--Helmholtz (KH) and tearing instabilities are likely to be important for the process of fast magnetic reconnection that is believed to explain the observed explosive energy release in solar flares. Theoretical studies of the instabilities, however, typically invoke simplified initial magnetic and velocity fields that are not solutions of the governing magnetohydrodynamic (MHD) equations. In the present study, the stability of a reconnecting current sheet is examined using a class of exact global MHD solutions for steady state incompressible magnetic reconnection, discovered by Craig & Henton. Numerical simulation indicates that the outflow solutions where the current sheet is formed by strong shearing flows are subject to the KH instability. The inflow solutions where the current sheet is formed by a fast and weakly sheared inflow are shown to be tearing unstable. Although the observed instability of the solutions can be interpreted qualitatively by applying standard linear results for the KH and tearing instabilities, the magnetic field and plasma flow, specified by the Craig--Henton solution, lead to the stabilization of the current sheet in some cases. The sensitivity of the instability growth rate to the global geometry of magnetic reconnection may help in solving the trigger problem in solar flare research.

A selection of results from the general relativistic MHD accretion simulations described in the previous talk are presented. We find that the magnetic field strength increases sharply with decreasing radius and is also enhanced near rapidly-spinning black holes. The greater magnetic field strength associated with rapid black hole rotation leads to a large outward electromagnetic angular momentum flux that substantially reduces both the mean accretion rate and the net accreted angular momentum per unit rest-mass. This electromagnetic stress strongly violates the traditional guess that the accretion stress vanishes at and inside the marginally stable orbit. Possible observational consequences include a constraint on the maximum spin of black holes, enhancement to the radiative efficiency, and concentration of fluorescent Fe Kalpha to the innermost part of the accretion disk.

Distributions of Faraday rotation measure (FRM) and the projected magnetic field derived by a 3-dimensional simulation of MHD jets are investigated based on our "sweeping magnetic twist model". FRM and Stokes parameters were calculated to be compared with radio observations of large scale wiggled AGN jets on kpc scales. We propose that the FRM distribution can be used to discuss the 3-dimensional structure of magnetic field around jets and the validity of existing theoretical models, together with the projected magnetic field derived from Stokes parameters. In the previous paper, we investigated the basic straight part of AGN jets by using the result of a 2-dimensional axisymmetric simulation. The derived FRM distribution has a general tendency to have a gradient across the jet axis, which is due to the toroidal component of the magnetic field generated by the rotation of the accretion disk. In this paper, we consider the wiggled structure of the AGN jets by using the result of a 3-dimensional simulation. Our numerical results show that the distributions of FRM and the projected magnetic field have a clear correlation with the large scale structure of the jet itself, namely, 3-dimensional helix. Distributions, seeing the jet from a certain direction, show a good matching with those in a part of 3C449 jet. This suggests that the jet has a helical structure and that the magnetic field (especially the toroidal component) plays an important role in the dynamics of the wiggle formation because it is due to a current-driven helical kink instability in our model.

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.

Using the numerical data of MHD simulation for AGN jets based on our ``sweeping magnetic twist model'', we calculated the Faraday rotation measure (FRM) and the Stokes parameters to compare with observations. We propose that the FRM distribution can be used to discuss the 3-dimensional structure of magnetic field around jets, together with the projected magnetic field derived from the Stokes parameters. In the present paper, we supposed the basic straight part of AGN jet, and used the data of axisymmetric simulation. The FRM distribution we derived has a general tendency to have gradient across the jet axis, which is due to the toroidal component of the helical magnetic field generated by the rotation of the accretion disk. This kind of gradient in the FRM distribution is actually observed in some AGN jets (e.g. Asada et al. 2002), which suggests helical magnetic field around the jets and thus supports our MHD model. Following this success, we are now extending our numerical observation to the wiggled part of the jets using the data of 3-dimensional simulation based on our model in the following paper.