Numerical simulations of the force-free magnetodynamics (FFMD) of the electromagnetic field around a spinning black hole are useful to investigate the dynamic electromagnetic processes around a spinning black hole, such as the emergence of the Blandford--Znajek mechanism. To reveal the basic physics of magnetic fields around a black hole through the dynamic process, we use one-dimensional (1D) FFMD along the axisymmetric magnetic surface, which provides a relatively simple, sufficiently precise, and powerful tool to analyze the dynamic process around a spinning black hole. We review the analytic and numerical aspects of 1D FFMD for an arbitrary magnetic surface around a black hole. In addition, we also show some numerical simulation test results for three types of magnetic surfaces at the equatorial plane of the black hole.

We performed numerical simulations of general relativistic magnetohydrodynamics with uniform resistivity to investigatethe occurrence of magnetic reconnection in a split-monopole magnetic field around a Schwarzschild black hole. We found that magnetic reconnection happens near the black hole at its equatorial plane. The magnetic reconnection has a point-like reconnection region and slow shock waves, as in the Petschek reconnection model. The magnetic reconnection rate decreases as the resistivity becomes smaller. When the global magnetic Reynolds number is $10^4$ or larger, the magnetic reconnection rate increases linearly with time from $2 \tau_{\rm S}$ to $\sim 10 \tau_{\rm S}$ ($\tau_{\rm S}=r_{\rm S}/c, r_{\rm S}$ is the Schwarzschild radius and $c$ is the speed of light). The linear increase of the reconnection rate agrees with the magnetic reconnection in the Rutherford regime of the tearing mode instability.

Recent general relativistic MHD simulations have suggested that relativistic jets from active galactic nuclei (AGNs) have been powered by rotational energy of central black holes. Some mechanisms of extraction of black hole rotational energy have been proposed, like the Penrose process, Blandford-Znajek mechanism, MHD Penrose process, and superradiance. The Blandford-Znajek mechanism is the most promising mechanism for the engines of the relativistic jets from AGNs. However, an intuitive interpretation of this mechanism with causality is not clarified yet, while the Penrose process has a clear interpretation for the causal energy extraction from the black hole with negative energy. In this paper, we present a formula to build physical intuition so that in the Blandford-Znajek mechanism as well as other electromagnetic processes, negative electromagnetic energy plays an important role to extract the rotational energy of the black holes causally.

Group velocity of electromagnetic waves in plasmas derived by standard relativistic resistive MHD (resistive RMHD) equations is superluminal. If we assume that the group velocity represents the propagation velocity of a signal, we have to worry about the causality problem. That is, some acausal phenomena may be induced, such that information transportation to the absolute past and spontaneous decrease in the entropy. Here, we tried to find the acausal phenomena using standard resistive RMHD numerical simulations in the suggested situation of the acausal phenomena. The calculation results showed that even in such situations no acausal effect happens. The numerical result with respect to the velocity limit of the information transportation is consistent with a linear theory of wave train propagation. Our results assure that we can use these equations without problems of acausal phenomena.

In almost all of plasma theories for astrophysical objects, we have assumed the charge quasi-neutrality of unmagnetized plasmas in global scales. This assumption has been justified because if there is a charged plasma, it induces electric field which attracts the opposite charge, and this opposite charge reduces the charge separation. Here, we report a newly discovered instability which causes a charge separation in a rotating plasma inside of an innermost stable circular orbit (ISCO) around a black hole. The growth rate of the instability is smaller than that of the disk instability even in the unstable disk region and is forbidden in the stable disk region outside of the ISCO. However, this growth rate becomes comparable to that of the disk instability when the plasma density is much lower than a critical density inside of the ISCO. In such case, the charge separation instability would become apparent and cause the charged accretion into the black hole, thus charge the hole up.

To study phenomena of plasmas around rotating black holes, we have derived a set of 3+1 formalism of generalized general relativistic magnetohydrodynamic (GRMHD) equations. Especially, we investigated general relativistic phenomena with respect to the Ohm's law. We confirmed the electromotive force due to the gravitation, centrifugal force, and frame-dragging effect in plasmas near the black holes. These effects are significant only in the local small-scale phenomena compared to the scale of astrophysical objects. We discuss the possibility of magnetic reconnection, which is triggered by one of these effects in a small-scale region and influences the plasmas globally. We clarify the conditions of applicability of the generalized GRMHD, standard resistive GRMHD, and ideal GRMHD for plasmas in black hole magnetospheres.

We derived one-fluid equations based on a relativistic two-fluid approximation of e$^\pm$ pair plasma and electron-ion plasma to reveal the specific relativistic nature of their behavior. Assuming simple condition on the relativistic one-fluid equations, we propose generalized relativistic magnetohydrodynamic (RMHD) equations which satisfy causality. We show the linear analyses of these equations regarding various plasma waves to show the validity of the generalized RMHD equations derived here and to reveal the distinct properties of the pair plasma and electron-ion plasma. The distinct properties relate to (i) the inertia effect of electric charge, (ii) the momentum of electric current, (iii) the relativistic Hall effect, (iv) the thermal electromotive force, and (v) the thermalized energy exchange between the two fluids. Using the generalized RMHD equations, we also clarify the condition that we can use standard RMHD equations and that we need the distinct RMHD equations of pair and electron-ion plasmas. The standard RMHD is available only when the relative velocity of the two fluids is nonrelativistic, a difference of the enthalpy densities of the two fluids is much smaller than the total enthalpy density, and the above distinct properties of the pair/electron-ion plasma are negligible. We discuss a general relativistic version of the equations applicable to the pair and electron-ion plasmas in black hole magnetospheres. We find the effective resistivity due to shear of frame ragging around a rotating black hole.

We investigate mechanisms of energy extraction from a rotating black hole in terms of negative energy-at-infinity. In addition to the Penrose process through particle fission, the Blandford-Znajek mechanism by magnetic tension, and the magnetohydrodynamic Penrose process, we examine energy extraction from a black hole caused by magnetic reconnection in the ergosphere. The reconnection redistributes the angular momentum efficiently to yield the negative energy-atinfinity. We derive a condition for the process to operate in a simple situation, where the plasma is incompressible and the magnetic energy is converted completely to the plasma kinetic energy locally. Astrophysical situations of magnetic reconnection around the black holes are also discussed.

We have performed several simulations of black hole systems (non-rotating, black hole spin parameter a=0.0 and rapidly rotating, a=0.95) with a geometrically thin Keplerian disk using the newly developed RAISHIN code. The simulation results show the formation of jets driven by the Lorentz force and the gas pressure gradient. The jets have mildly relativistic speed (\sim 0.4c). The matter is continuously supplied from the accretion disk and the jet propagates outward until each applicable terminal simulation time (non-rotating: t/\tau_S = 275 and rotating: t/tau_S = 200, \tau_S \equiv r_S/c). It appears that a rotating black hole creates an additional, faster, and more collimated matter-dominated inner outflow (\sim 0.5c) formed and accelerated by the twisted magnetic field resulting from frame-dragging in the black hole ergosphere. This is the first known simulation confirming the formation of an inner magnetically-driven, matter-dominated jet by the frame-dragging effect from a black hole co-rotating with a thin Keplerian disk threaded by a vertical magnetic field. This result indicates that jet kinematic structure depends on black hole rotation and on the initial magnetic field configuration and strength.

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity nearly 0.3c) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around r = 3r_S. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfven waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

We have performed a fully three-dimensional general relativistic magnetohydrodynamic (GRMHD) simulation of jet formation from a thin accretion disk around a Schwarzschild black hole with a free-falling corona. The initial simulation results show that a bipolar jet (velocity $\sim 0.3c$) is created as shown by previous two-dimensional axisymmetric simulations with mirror symmetry at the equator. The 3-D simulation ran over one hundred light-crossing time units ($\tau_{\rm S} = r_{\rm S}/c$ where $r_{\rm S} \equiv 2GM/c^2$) which is considerably longer than the previous simulations. We show that the jet is initially formed as predicted due in part to magnetic pressure from the twisting the initially uniform magnetic field and from gas pressure associated with shock formation in the region around $r = 3 r_{\rm S}$. At later times, the accretion disk becomes thick and the jet fades resulting in a wind that is ejected from the surface of the thickened (torus-like) disk. It should be noted that no streaming matter from a donor is included at the outer boundary in the simulation (an isolated black hole not binary black hole). The wind flows outwards with a wider angle than the initial jet. The widening of the jet is consistent with the outward moving torsional Alfv\'{e}n waves (TAWs). This evolution of disk-jet coupling suggests that the jet fades with a thickened accretion disk due to the lack of streaming material from an accompanying star.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of collapsars including a rotating black hole. This paper is an extension of our previous paper (Mizuno et al. 2004). The current calculation focuses on the effect of black hole rotation using general relativistic MHD with simplified microphysics. Initially, we assume that the core collapse is failed in this star. A few solar mass rotating black hole is inserted by hand into the calculation. We consider two cases, the co-rotating case and counter-rotating case with respect to the black hole rotation. Although the counter-rotating case may be unrealistic for collapsar, we perform as a maximally dragging case of magnetic field. The simulation results show the formation of a disk-like structure and the generation of a jet-like outflow near the central black hole. The jet-like outflow propagates outwardly with the twisted magnetic field and becomes collimated. We have found that the jets are generated and accelerated mainly by the magnetic field. The total jet velocity in the rotating black hole case is comparable to that of the non-rotating black hole case (Mizuno et al. 2004), 0.3c. When the rotation of the black hole is faster, the magnetic field is twisted strongly owing to the frame-dragging effect. The magnetic energy stored by the twisting magnetic field is converted to kinetic energy of the jet directly rather than propagating as an Alfven wave. Thus, as the rotation of the black hole becomes faster, the poloidal velocity of the jet becomes faster. In the rapidly rotating black hole case the jet-like outflow can be produced by the frame dragging effect only through the twisting of magnetic field even if there is no stellar rotation.

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). The current calculation focuses on general relativistic MHD with simplified microphysics (we ignore neutrino cooling, physical equation of state, and photodisintegration). Initially, we assume that the core collapse has failed in this star. A few $M_{\odot}$ black hole is inserted by hand into the calculation. The simulations presented in the paper follow the accretion of gas into a black hole that is assumed to have formed before the calculation begins.The simulation results show the formation of a disk-like structure and the generation of a jetlike outflow inside the shock wave launched at the core bounce. We have found that the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic (0.3c). The velocity of the jet becomes larger as the initial rotational velocity of stellar matter gets faster. On the other hand, the dependence on the initial magnetic field strength is a bit more complicated: the velocity of the jet increases with the initial field strength in the weak field regime, then is saturated at some intermediate field strength, and decreases beyond the critical field strength. These results are related to the stored magnetic energy determined by the balance between the propagation time of the Alfven wave and the rotation time of the disk (or twisting time).

We have performed 2.5-dimensional general relativistic magnetohydrodynamic (MHD) simulations of the gravitational collapse of a magnetized rotating massive star as a model of gamma ray bursts (GRBs). This simulation showed the formation of a disk-like structure and the generation of a jet-like outflow inside the shock wave launched at the core bounce. We have found the jet is accelerated by the magnetic pressure and the centrifugal force and is collimated by the pinching force of the toroidal magnetic field amplified by the rotation and the effect of geometry of the poloidal magnetic field. The maximum velocity of the jet is mildly relativistic ($\sim$ 0.3 c).

We have performed a first fully 3-D GRMHD simulation with Schwarzschild black hole with a free falling corona. The initial simulation results show that a jet is created as in the previous simulations using the axisymmetric geometry with the mirror symmetry at the equator. However, the time to generate the jet is longer than in the 2-D axisymmetric simulations. We expect that due to the additional freedom in the azimuthal dimension without axisymmetry with respect to the z axis and reflection symmetry with respect to the equatorial plane, the dynamics of jet formation can be modified. Further simulations are required for the study of instabilities along the azimuthal direction such as accretion-eject instability

We have performed a first fully 3-D GRMHD simulation with Schwarzschild black hole with a free falling corona. The initial simulation results show that a jet is created as in previous axisymmetric simulations. However, the time to generate the jet is longer than in the 2-D simulations. We expect that due to the additional azimuthal dimension the dynamics of jet formation can be modified.

To investigate the formation mechanism of relativistic jets in active galactic nuclei and micro-quasars, we have developed a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we report on the first numerical simulation of jet formation in a rapidly-rotating (a=0.95) Kerr black hole magnetosphere. We study cases in which the Keplerian accretion disk is both co-rotating and counter-rotating with respect to the black hole rotation. In the co-rotating disk case, our results are almost the same as those in Schwarzschild black hole cases: a gas pressure-driven jet is formed by a shock in the disk, and a weaker magnetically-driven jet is also generated outside the gas pressure-driven jet. On the other hand, in the counter-rotating disk case, a new powerful magnetically-driven jet is formed inside the gas pressure-driven jet. The newly found magnetically-driven jet in the latter case is accelerated by a strong magnetic field created by frame dragging in the ergosphere. Through this process, the magnetic field extracts the energy of the black hole rotation.

Recent observations of Galactic Black Hole Candidates (BHCs) suggest that those that are superluminal jet sources have more rapid black hole spin rates than otherwise normal BHCs. This provides observational support for models of relativistic jet formation that extract rotational energy of the central black hole. To investigate this mechanism, we have developed a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we report on the first numerical simulation of the formation of a relativistic jet in a rapidly-rotating (a=0.95) Kerr black hole magnetosphere. We assume that the initial velocity of the disk is zero. We find that the maximum velocity of the jet reaches 0.93c (Lorentz factor, 2.7) and the terminal velocity of the jet is 0.85c (Lorentz factor, 1.9). On the other hand, for a non-rotating (a=0) Schwarzschild black hole, the maximum outflow velocity is less than 0.6c for initial magnetospheric conditions similar to to those of the Kerr black hole case. These numerical results show the importance of the rapidly rotating black hole for the relativistic jet formation.