We introduce a theory for the development of a transitional column density $\Sigma_{\rm TP}$ between the lognormal and the power-law forms of the probability distribution function (PDF) in a molecular cloud. Our turbulent magnetohydrodynamic simulations show that the value of $\Sigma_{\rm TP}$ increases as the strength of both the initial magnetic field and turbulence increases. We develop an analytic expression for $\Sigma_{\rm TP}$ based on the interplay of turbulence, a (strong) magnetic field, and gravity. The transition value $\Sigma_{\rm TP}$ scales with $\mathcal{M}^2_{\rm 0}$, the square of the initial sonic Mach number, and $\beta_{0}$, the initial ratio of gas pressure to magnetic pressure. We fit the variation of $\Sigma_{\rm TP}$ among different model clouds as a function of $\mathcal{M}^2_{\rm 0} \beta_{0}$, or equivalently the square of the initial Alfv\'enic Mach number $\mathcal{M}^2_{\rm A0}$. This implies that the transition value $\Sigma_{\rm TP}$ is an imprint of cloud initial conditions and is set by turbulent compression of a magnetic cloud. Physically, the value of $\Sigma_{\rm TP}$ denotes the boundary above which the mass-to-flux ratio becomes supercritical and gravity drives the evolution.

We examine the linear stability of a filamentary cloud permeated by a perpendicular magnetic field. The initial magnetic field is assumed to be uniform and perpendicular to the cloud axis. The model cloud is assumed to have a Plummer-like density profile and to be supported against the self-gravity by turbulence. The effects of turbulence are taken into account by enhancing the effective pressure of a low density gas. We derive the effective pressure as a function of the density from the condition of the hydrostatic balance. It is shown that the model cloud is more unstable against radial collapse, when the radial density slope is shallower. When the magnetic field is mildly strong, the radial collapse is suppressed. If the displacement vanishes in the region very far from the cloud axis, the model cloud is stabilized completely by a mildly strong magnetic field. If rearrangement of the magnetic flux tubes is permitted, the model cloud is unstable even when the magnetic field is extremely strong. The stability depends on the outer boundary condition as in case of the isothermal cloud. The growth rate of the rearrangement mode is smaller when the radial density slope is shallower.

Simulations generally show that non-self-gravitating clouds have a lognormal column density ($\Sigma$) probability distribution function (PDF), while self-gravitating clouds with active star formation develop a distinct power-law tail at high column density. Although the growth of the power law can be attributed to gravitational contraction leading to the formation of condensed cores, it is often debated if an observed lognormal shape is a direct consequence of supersonic turbulence alone, or even if it is really observed in molecular clouds. In this paper we run three-dimensional magnetohydrodynamic simulations including ambipolar diffusion with different initial conditions to see the effect of strong magnetic fields and nonlinear initial velocity perturbations on the evolution of the column density PDFs. Our simulations show that column density PDFs of clouds with supercritical mass-to-flux ratio, with either linear perturbations or nonlinear turbulence, quickly develop a power-law tail such that $dN/d \log \Sigma \propto \Sigma^{-\alpha}$ with index $\alpha \simeq 2$. Interestingly, clouds with subcritical mass-to-flux ratio also proceed directly to a power-law PDF, but with a much steeper index $\alpha \simeq 4$. This is a result of gravitationally-driven ambipolar diffusion. However, for nonlinear perturbations with a turbulent spectrum ($v_{k}^{2} \propto k^{-4}$), the column density PDFs of subcritical clouds do retain a lognormal shape for a major part of the cloud evolution, and only develop a distinct power-law tail with index $\alpha \simeq 2$ at greater column density when supercritical pockets are formed.

We examine the linear stability of an isothermal filamentary cloud permeated by a perpendicular magnetic field. Our model cloud is assumed to be supported by gas pressure against the self-gravity in the unperturbed state. For simplicity, the density distribution is assumed to be symmetric around the axis. Also for simplicity, the initial magnetic field is assumed to be uniform and turbulence is not taken into account. The perturbation equation is formulated to be an eigenvalue problem. The growth rate is obtained as a function of the wavenumber for fragmentation along the axis and the magnetic field strength. The growth rate depends critically on the outer boundary. If the displacement vanishes in the region very far from the cloud axis (fixed boundary), cloud fragmentation is suppressed by a moderate magnetic field, which means the plasma beta is below 1.67 on the cloud axis. If the displacement is constant along the magnetic field in the region very far from the cloud, the cloud is unstable even when the magnetic field is infinitely strong. The cloud is deformed by circulation in the plane perpendicular to the magnetic field. The unstable mode is not likely to induce dynamical collapse, since it is excited even when the whole cloud is magnetically subcritical. For both the boundary conditions the magnetic field increases the wavelength of the most unstable mode. We find that the magnetic force suppresses compression perpendicular to the magnetic field especially in the region of low density.

Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on the achievements that have been provided by centimeter and meter wavelength radio observations. The article also introduces various methods to analyze linear polarization data, including synchrotron radiation, Faraday rotation, depolarization, and Faraday tomography.

We develop a magnetic ribbon model for molecular cloud filaments. These result from turbulent compression in a molecular cloud in which the background magnetic field sets a preferred direction. We use our model to calculate a synthetic observed relation between apparent width in projection versus observed column density. The relationship is relatively flat, in rough agreement with the observations, and unlike the simple expectation based on a Jeans length argument.

We develop a magnetic ribbon model for molecular cloud filaments. These result from turbulent compression in a molecular cloud in which the background magnetic field sets a preferred direction. We argue that this is a natural model for filaments and is based on the interplay between turbulence, strong magnetic fields, and gravitationally-driven ambipolar diffusion, rather than pure gravity and thermal pressure. An analytic model for the formation of magnetic ribbons that is based on numerical simulations is used to derive a lateral width of a magnetic ribbon. This differs from the thickness along the magnetic field direction, which is essentially the Jeans scale. We use our model to calculate a synthetic observed relation between apparent width in projection versus observed column density. The relationship is relatively flat, similar to observations, and unlike the simple expectation based on a Jeans length argument.

Magnetic fields play essential roles in various astronomical objects. Radio astronomy has revealed that magnetic fields are ubiquitous in our Universe. However, the real origin and evolution of magnetic fields is poorly proven. In order to advance our knowledge of cosmic magnetism in coming decades, the Square Kilometre Array (SKA) should have supreme sensitivity than ever before, which provides numerous observation points in the cosmic space. Furthermore, the SKA should be designed to facilitate wideband polarimetry so as to allow us to examine sightline structures of magnetic fields by means of depolarization and Faraday Tomography. The SKA will be able to drive cosmic magnetism of the interstellar medium, the Milky Way, galaxies, AGN, galaxy clusters, and potentially the cosmic web which may preserve information of the primeval Universe. The Japan SKA Consortium (SKA-JP) Magnetism Science Working Group (SWG) proposes the project "Resolving 4-D Nature of Magnetism with Depolarization and Faraday Tomography", which contains ten scientific use cases.

We clarify the mechanism of accelerated core formation by large-scale nonlinear flows in subcritical magnetic clouds by finding a semi-analytical formula for the core formation time and describing the physical processes that lead to them. Recent numerical simulations show that nonlinear flows induce rapid ambipolar diffusion that leads to localized supercritical regions that can collapse. Here, we employ non-ideal magnetohydrodynamic simulations including ambipolar diffusion for gravitationally stratified sheets threaded by vertical magnetic fields. One of the horizontal dimensions is eliminated, resulting in a simpler two-dimensional simulation that can clarify the basic process of accelerated core formation. A parameter study of simulations shows that the core formation time is inversely proportional to the square of the flow speed when the flow speed is greater than the Alfv\'en speed. We find a semi-analytical formula that explains this numerical result. The formula also predicts that the core formation time is about three times shorter than that with no turbulence, when the turbulent speed is comparable to the Alfv\'en speed.

We investigate effects of cosmic-rays on the linear growth of the Kelvin-Helmholtz instability. Cosmic-rays are treated as an adiabatic gas and allowed to diffuse along magnetic field lines. We calculated the dispersion relation of the instability for various sets of two free parameters, the ratio of the cosmic-ray pressure to the thermal gas pressure and the diffusion coefficient. Including cosmic-ray effects, a shear layer is more destabilized and the growth rates can be enhanced in comparison with the ideal magnetohydrodynamical case. Whether the growth rate is effectively enhanced or not depends on the diffusion coefficient of cosmic-rays. We obtain the criterion for effective enhancement by comparing the growing time scale of the instability with the diffusion time scale of cosmic-rays. These results can be applied to various astrophysical phenomena where a velocity shear is present, such as outflows from star-forming galaxies, AGN jet, channel flows resulting from the nonlinear development of the magnetorotational instability, and galactic disks.

We carried out global three-dimensional magneto-hydrodynamic 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 magneto-rotational 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 $\beta$ ($ = p_{\rm gas}/p_{\rm mag}$) becomes less than 5 near the disk surface, magnetic flux escapes from the disk by 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 MRI drive quasi-periodic reversal of azimuthal magnetic fields in timescale of 10 rotation period. 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.

Relativistic Sweet-Parker type magnetic reconnection is investigated by relativistic resistive magnetohydrodynamic (RRMHD) simulations. As an initial setting, we assume anti-parallel magnetic fields and a spatially uniform resistivity. A perturbation imposed on the magnetic fields triggers magnetic reconnection around a current sheet, and the plasma inflows into the reconnection region. The inflows are then heated due to ohmic dissipation in the diffusion region, and finally become relativistically hot outflows. The outflows are not accelerated to ultra-relativistic speeds (i.e., Lorentz factor ~ 1), even when the magnetic energy dominates the thermal and rest mass energies in the inflow region. Most of the magnetic energy in the inflow region is converted into the thermal energy of the outflow during the reconnection process. The energy conversion from magnetic to thermal energy in the diffusion region results in an increase in the plasma inertia. This prevents the outflows from being accelerated to ultra-relativistic speeds. We find that the reconnection rate R obeys the scaling relation R S^{-0.5}, where S is the Lundquist number. This feature is the same as that of non-relativistic reconnection. Our results are consistent with the theoretical predictions of Lyubarsky (2005) for Sweet-Parker type magnetic reconnection.

We employ the three-dimensional magnetohydrodynamic simulation including ambipolar diffusion to study the gravitationally-driven fragmentation of subcritical molecular clouds, in which the gravitational fragmentation is stabilized as long as magnetic flux-freezing applies. The simulations show that the cores in an initially subcritical cloud generally develop gradually over an ambipolar diffusion time, which is about a few $\times 10^7$ years in a typical molecular cloud. On the other hand, the formation of collapsing cores in subcritical clouds is accelerated by supersonic nonlinear flows. Our parameter study demonstrates that core formation occurs faster as the strength of the initial flow speed in the cloud increases. We found that the core formation time is roughly proportional to the inverse of the square root of the enhanced density created by the supersonic nonlinear flows. The density dependence is similar to that derived in quasistatically contracting magnetically supported clouds, although the core formation conditions are created by the nonlinear flows in our simulations. We have also found that the accelerated formation time is not strongly dependent on the initial strength of the magnetic field if the cloud is highly subcritical. Our simulation shows that the core formation time in our model subcritical clouds is several $\times 10^6$ years, due to the presence of large-scale supersonic flows ($\sim 3$ times sound speed). Once a collapsing core forms, the density, velocity, and magnetic field structure of the core does not strongly depend on the initial strength of the velocity fluctuation. The infall velocities of the cores are subsonic and the magnetic field lines show weak hourglass shapes.

We propose a primordial-origin model for the composite configurations of global magnetic fields in spiral galaxies. We show that uniform tilted magnetic field wound up into a rotating disk galaxy can evolve into composite magnetic configurations comprising bisymmetric spiral (S=BSS), axisymmetric spiral (A=ASS), plane-reversed spiral (PR), and/or ring (R) fields in the disk, and vertical (V) fields in the center. By MHD simulations we show that these composite galactic fields are indeed created from weak primordial uniform field, and that the different configurations can co-exist in the same galaxy. We show that spiral fields trigger the growth of two-armed gaseous arms. The centrally accumulated vertical fields are twisted and produce jet toward the halo. We find that the more vertical was the initial uniform field, the stronger is the formed magnetic field in the galactic disk.

Alfven waves created by sub-photospheric motions or by magnetic reconnection in the low solar atmosphere seem good candidates for coronal heating. However, the corona is also likely to be heated more directly by magnetic reconnection, with dissipation taking place in current sheets. Distinguishing observationally between these two heating mechanisms is an extremely difficult task. We perform 1.5-dimensional MHD simulations of a coronal loop subject to each type of heating and derive observational quantities that may allow these to be differentiated.

The unusual mushroom-shaped HI cloud, GW 123.4--1.5, is hundreds of parsecs in size but does not show any correlations to HI shells or chimney structures. To investigate the origin and velocity structure of GW 123.4--1.5, we perform three-dimensional hydrodynamical simulations of the collision of a high-velocity cloud with the Galactic disk. We also perform a parameter study of the density, radius, and incident angle of the impact cloud. The numerical experiments indicate that we reproduce the mushroom-shaped structure which resembles GW 123.4--1.5 in shape, size, position-velocity across the cap of the mushroom, and the density ratio between the mushroom and surrounding gas. GW 123.4--1.5 is expected to be formed by the almost head-on collision of a HVC with velocity $\sim 100 \kms$ and mass $\sim 10^5 \Msun$ about $5 \times 10^7 \yr$ ago. A mushroom-shaped structure like GW 123.4--1.5 must be infrequent on the Galactic plane, because the head-on collision which explains the mushroom structure seems rare for observed HVCs. HVC-disk collision explains not only the origin of the mushroom-shaped structure but also the formation of a variety of structures like shells, loops, and vertical structures in our Galaxy.

We demonstrate that the formation of collapsing cores in subcritical clouds is accelerated by nonlinear flows, by performing three-dimensional non-ideal MHD simulations. An initial random supersonic (and trans-Alfvenic) turbulent-like flow is input into a self-gravitating gas layer that is threaded by a uniform magnetic field (perpendicular to the layer) such that the initial mass-to-flux ratio is subcritical. Magnetic ambipolar diffusion occurs very rapidly initially due to the sharp gradients introduced by the turbulent flow. It subsequently occurs more slowly in the traditional near-quasistatic manner, but in regions of greater mean density than present in the initial state. The overall timescale for runaway growth of the first core(s) is several times, 10^6 yr, even though previous studies have found a timescale of several times, 10^7 yr when starting with linear perturbations and similar physical parameters. Large-scale supersonic flows exist in the cloud and provide an observationally testable distinguishing characteristic from core formation due to linear initial perturbations. However, the nonlinear flows have decayed sufficiently that the relative infall motions onto the first core are subsonic, as in the case of starting from linear initial perturbations. The ion infall motions are very similar to those of neutrals; however, they lag the neutral infall in directions perpendicular to the mean magnetic field direction and lead the neutral infall in the direction parallel to the mean magnetic field.

We employ the first fully three-dimensional simulation to study the role of magnetic fields and ion-neutral friction in regulating gravitationally-driven fragmentation of molecular clouds. The cores in an initially subcritical cloud develop gradually over an ambipolar diffusion time while the cores in an initially supercritical cloud develop in a dynamical time. The infall speeds on to cores are subsonic in the case of an initially subcritical cloud, while an extended (\ga 0.1 pc) region of supersonic infall exists in the case of an initially supercritical cloud. These results are consistent with previous two-dimensional simulations. We also found that a snapshot of the relation between density (rho) and the strength of the magnetic field (B) at different spatial points of the cloud coincides with the evolutionary track of an individual core. When the density becomes large, both relations tend to B \propto \rho^{0.5}.

We consider MHD waves as a heating source of cool cores of galaxy clusters. In particular, we focus on transverse waves (Alfven waves), because they can propagate a longer distance than longitudinal waves (sound waves). Using MHD simulations, we found that the transverse waves can stably heat a cool core if the wave period is large enough (>~ 10^8 yr). Moreover, the longitudinal waves that are created as a by-product of the nonlinear evolution of the transverse waves could be observed as the 'ripples' found in cool cores.

We use numerical simulations to study the effect of nonlinear MHD waves in a stratified, self-gravitating molecular cloud that is bounded by a hot and tenuous external medium. In a previous paper, we had shown the details of a standard model and studied the effect of varying the dimensionless amplitude. In this paper, we present the results of varying two other important free parameters: beta_0, the initial ratio of gas to magnetic pressure at the cloud midplane, and the dimensionless frequency of driving. Furthermore, we present the case of a temporally random driving force. Our results demonstrate that a very important consideration for the actual level of turbulent support against gravity is the ratio of driving wavelength lambda_0 to the the size of the initial non-turbulent cloud; maximum cloud expansion is achieved when this ratio is close to unity. The best consistency with the observational correlation of magnetic field strength, turbulent line width, and density is achieved by cloud models with beta_0 approx 1. We also calculate the spatial power spectra of the turbulent clouds, and show that significant power is developed on scales larger than the scale length H_0 of the initial cloud, even if the input wavelength of turbulence lambda_0 approx H_0. The cloud stratification and resulting increase of Alfven speed toward the cloud edge allows for a transfer of energy to wavelengths significantly larger than lambda_0. This explains why the relevant time scale for turbulent dissipation is the crossing time over the cloud scale rather than the crossing time over the driving scale.

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 analyzed the results of non-linear resistive magnetohydrodynamical (MHD) simulations of jet formation to study the acceleration mechanism of axisymmetric, resistive MHD jets. The initial state is a constant angular momentum, polytropic torus threaded by weak uniform vertical magnetic fields. The time evolution of the torus is simulated by applying the CIP-MOCCT scheme extended for resistive MHD equations. We carried out simulations up to 50 rotation period at the innermost radius of the disk created by accretion from the torus. The acceleration forces and the characteristics of resistive jets were studied by computing forces acting on Lagrangian test particles. Since the angle between the rotation axis of the disk and magnetic field lines is smaller in resistive models than in ideal MHD models, magnetocentrifugal acceleration is smaller. The effective potential along a magnetic field line has maximum around $z \sim 0.5r_0$ in resistive models, where $r_0$ is the radius where the density of the initial torus is maximum. Jets are launched after the disk material is lifted to this height by pressure gradient force. Even in this case, the main acceleration force around the slow magnetosonic point is the magnetocentrifugal force. The power of the resistive MHD jet is comparable to the mechanical energy liberated in the disk by mass accretion. Joule heating is not essential for the formation of jets.

We propose that the mushroom-shaped structure of the Galactic worm GW 123.4--1.5 is created by a cloud collision with the Galactic gas disk. A hydrodynamic simulation shows that a mushroom-shaped structure is created after the cloud crosses the Galactic midplane. The lifetime of the mushroom-shaped structure is of order the dynamical time scale of the disk, \sim 10^7 years. We find that the velocities across the cap of the mushroom-shaped structure in the simulation are consistent with the observed values. The simulation also predicts a structure on the opposite side of the Galactic plane which is created by the Kelvin-Helmholtz instability after the cloud passes through the disk.

We perform numerical simulations of nonlinear MHD waves in a gravitationally stratified molecular cloud that is bounded by a hot and tenuous external medium. We study the relation between the strength of the turbulence and various global properties of a molecular cloud, within a 1.5-dimensional approximation. Under the influence of a driving source of Alfvenic disturbances, the cloud is lifted up by the pressure of MHD waves and reaches a steady-state characterized by oscillations about a new time-averaged equilibrium state. The nonlinear effect results in the generation of longitudinal motions and many shock waves; however, the wave kinetic energy remains predominantly in transverse, rather than longitudinal, motions. There is an approximate equipartition of energy between the transverse velocity and fluctuating magnetic field (aspredicted by small-amplitude theory) in the region of the stratified cloud which contains most of the mass; however, this relation breaks down in the outer regions, particularly near the cloud surface, where the motions have a standing-wave character. This means that the Chandrasekhar-Fermi formula applied to molecular clouds must be significantly modified in such regions. Models of an ensemble of clouds show that, for various strengths of the input energy, the velocity dispersion in the cloud $\sigma \propto Z^{0.5}$, where $Z$ is a characteristic size of the cloud.Furthermore, $\sigma$ is always comparable to the mean Alfven velocity of the cloud, consistent with observational results.

We propose the Galactic flare model for the origin of the X-ray gas in the Galactic halo. For this purpose, we examine the magnetic reconnection triggered by Parker instability (magnetic buoyancy instability), by performing the two-dimensional resistive numerical magnetohydrodynamic simulations. As a result of numerical simulations, the system evolves as following phases: Parker instability occurs in the Galactic disk. In the nonlinear phase of Parker instability, the magnetic loop inflates from the Galactic disk into the Galactic halo, and collides with the anti-parallel magnetic field, so that the current sheets are created in the Galactic halo. The tearing instability occurs, and creates the plasmoids (magnetic islands). Just after the plasmoid ejection, further current-sheet thinning occurs in the sheet, and the anomalous resistivity sets in. Petschek reconnection starts, and heats the gas quickly in the Galactic halo. It also creates the slow and fast shock regions in the Galactic halo. The magnetic field ($B\sim 3 \mu$G), for example, can heat the gas ($n\sim 10^{-3}$ cm$^{-3}$) to temperature of $\sim 10^6$ K via the reconnection in the Galactic halo. The gas is accelerated to Alfv\'en velocity ($\sim 300$ km s$^{-1}$). Such high velocity jets are the evidence of the Galactic flare model we present in this paper, if the Doppler shift of the bipolar jet is detected in the Galactic halo. Full size figures are available at http://www.kwasan.kyoto-u.ac.jp/~tanuma/study/ApJ2002/ApJ2002.html

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.

We examine the magnetic reconnection triggered by a supernova (or a point explosion) in interstellar medium, by performing two-dimensional resistive magnetohydrodynamic (MHD) numerical simulations with high spatial resolution. We found that the magnetic reconnection starts long after a supernova shock (fast-mode MHD shock) passes a current sheet. The current sheet evolves as follows: (i) Tearing-mode instability is excited by the supernova shock, and the current sheet becomes thin in its nonlinear stage. (ii) The current-sheet thinning is saturated when the current-sheet thickness becomes comparable to that of Sweet-Parker current sheet. After that, Sweet-Parker type reconnection starts, and the current-sheet length increases. (iii) ``Secondary tearing-mode instability'' occurs in the thin Sweet-Parker current sheet. (iv) As a result, further current-sheet thinning occurs and anomalous resistivity sets in, because gas density decreases in the current sheet. Petschek type reconnection starts and heats interstellar gas. Magnetic energy is released quickly while magnetic islands are moving in the current sheet during Petschek type reconnection. The released magnetic energy is determined by the interstellar magnetic field strength, not energy of initial explosion nor distance to explosion. We suggest that magnetic reconnection is a possible mechanism to generate X-ray gas in Galaxy.

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.

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.