Based mainly on X-ray observations, studies are made on interactions between the intra-cluster medium (ICM) in clusters of galaxies and their member galaxies. Through (magneto)hydrodynamic and gravitational channels, the moving galaxies are expected to drag the ICM around them, and transfer to the ICM some fraction of their dynamical energies on cosmological time scales. This hypothesis is in line with several observations, including the possible cosmological infall of galaxies towards the cluster center, found over redshifts of z~1 to z~0. Further assuming that the energy lost by the galaxies is first converted into ICM turbulence and then dissipated, this picture can explain the subsonic and uniform ICM turbulence, measured with Hitomi in the core region of the Perseus cluster. The scenario may also explain several other unanswered problems regarding clusters of galaxies, including what prevents the ICM from the expected radiative cooling, how the various mass components in nearby clusters have attained different radial distributions, and how a thermal stability is realized between hot and cool ICM components that co-exist around cD galaxies. This view is also considered to pertain to the general scenario of galaxy evolution, including their environmental effects.

We present the results of two-temperature magnetohydrodynamic simulations of the propagation of sub-relativistic jets of active galactic nuclei. The dependence of the electron and ion temperature distributions on the fraction of electron heating fe at the shock front is studied for fe=0, 0.05, and 0.2. Numerical results indicate that in sub-relativistic, rarefied jets, the jet plasma crossing the terminal shock forms a hot, two-temperature plasma in which the ion temperature is higher than the electron temperature. The two-temperature plasma expands and forms a backflow referred to as a cocoon, in which the ion temperature remains higher than the electron temperature for longer than 100 Myr. Electrons in the cocoon are continuously heated by ions through Coulomb collisions, and the electron temperature thus remains at Te > 10^9 K in the cocoon. X-ray emissions from the cocoon are weak because the electron number density is low. Meanwhile, soft X-rays are emitted from the shocked intracluster medium surrounding the cocoon. Mixing of the jet plasma and the shocked intracluster medium through the Kelvin--Helmholtz instability at the interface enhances X-ray emissions around the contact discontinuity between the cocoon and shocked intracluster medium.

We report a search for astronomical neutrinos in the energy region from several GeV to TeV in the direction of the blazar TXS0506+056 using the Super-Kamiokande detector following the detection of a 100 TeV neutrino from the same location by the IceCube collaboration. Using Super-Kamiokande neutrino data across several data samples observed from April 1996 to February 2018 we have searched for both a total excess above known backgrounds across the entire period as well as localized excesses on smaller time scales in that interval. No significant excess nor significant variation in the observed event rate are found in the blazar direction. Upper limits are placed on the electron and muon neutrino fluxes at 90\% confidence level as $6.03 \times 10^{-7}$ and $4.52 \times 10^{-7}$ to $9.26 \times 10^{-10}$ [${\rm erg}/{\rm cm}^2/{\rm s}$], respectively.

Supernova detection is a major objective of the Super-Kamiokande (SK) experiment. In the next stage of SK (SK-Gd), gadolinium (Gd) sulfate will be added to the detector, which will improve the ability of the detector to identify neutrons. A core-collapse supernova will be preceded by an increasing flux of neutrinos and anti-neutrinos, from thermal and weak nuclear processes in the star, over a timescale of hours; some of which may be detected at SK-Gd. This could provide an early warning of an imminent core-collapse supernova, hours earlier than the detection of the neutrinos from core collapse. Electron anti-neutrino detection will rely on inverse beta decay events below the usual analysis energy threshold of SK, so Gd loading is vital to reduce backgrounds while maximising detection efficiency. Assuming normal neutrino mass ordering, more than 200 events could be detected in the final 12 hours before core collapse for a 15-25 solar mass star at around 200 pc, which is representative of the nearest red supergiant to Earth, $\mathrm{\alpha}$Ori (Betelgeuse). At a statistical false alarm rate of 1 per century, detection could be up to 10 hours before core collapse, and a pre-supernova star could be detected by SK-Gd up to 600 pc away. A pre-supernova alert could be provided to the astrophysics community following gadolinium loading.

We present a new magnetohydrodynamic (MHD) simulation code with the aim of providing accurate numerical solutions to astrophysical phenomena where discontinuities, shock waves, and turbulence are inherently important. The code implements the HLLD approximate Riemann solver, the fifth-order-monotonicity-preserving interpolation (MP5) scheme, and the hyperbolic divergence cleaning method for a magnetic field. This choice of schemes significantly improved numerical accuracy and stability, and saved computational costs in multidimensional problems. Numerical tests of one- and two-dimensional problems showed the advantages of using the high-order scheme by comparing with results from a standard second-order TVD MUSCL scheme. The present code enabled us to explore long-term evolution of a three-dimensional accretion disk around a black hole, in which compressible MHD turbulence caused continuous mass accretion via nonlinear growth of the magneto-rotational instability (MRI). Numerical tests with various computational cell sizes exhibited a convergent picture of the early nonlinear growth of the MRI in a global model, and indicated that the MP5 scheme has more than twice the resolution of the MUSCL scheme in practical applications.

We studied roles of the magnetic field on the gas dynamics in the Galactic bulge by a three-dimensional global magnetohydrodynamical simulation data, particularly focusing on vertical flows that are ubiquitously excited by magnetic activity. In local regions where the magnetic filed is stronger, it is frequently seen that fast down-flows slide along inclined magnetic field lines that are associated with buoyantly rising magnetic loops. The vertical velocity of these down-flows reaches ~ 100 km s$^{-1}$ near the foot-point of the loops by the gravitational acceleration toward the Galactic plane. The two footpoints of rising magnetic loops are generally located at different radial locations and the field lines are deformed by the differential rotation. The angular momentum is transported along the field lines, and the radial force balance breaks down. As a result, a fast downflow is often observed only at the one footpoint located at the inner radial position. The fast downflow compresses the gas to form a dense region near the footpoint, which will be important in star formation afterward. Furthermore, the horizontal components of the velocity are also fast near the foot-point because the down-flow is accelerated along the magnetic sliding slope. As a result, the high-velocity flow creates various characteristic features in a simulated position-velocity diagram, depending on the viewing angle.

A possibility of time-delayed radio brightenings of Sgr A* triggered by the pericenter passage of the G2 cloud is studied by carrying out global three-dimensional magnetohydrodynamic simulations taking into account the radiative cooling of the tidal debris of the G2 cloud. Magnetic fields in the accretion flow are strongly perturbed and re-organized after the passage of G2. We have found that the magnetic energy in the accretion flow increases by a factor 3-4 in 5-10 years after the pericenter passage of G2 by a dynamo mechanism driven by the magneto-rotational instability. Since this B-field amplification enhances the synchrotron emission from the disk and the outflow, the radio and the infrared luminosity of Sgr A* is expected to increase around A.D. 2020. The time-delay of the radio brightening enables us to determine the rotation axis of the preexisting disk.

The formation mechanism of CO clouds observed with NANTEN2 and Mopra telescope toward the stellar cluster Westerlund 2 is studied by three-dimensional magnetohydrodynamic (MHD) simulations taking into account the interstellar cooling. These molecular clouds show a peculiar shape composing of an arc-shaped cloud in one side of a TeV-ray source HESS J1023-575 and a linear distribution of clouds (jet clouds) in another side. We propose that these clouds are formed by the interaction of a jet with interstellar neutral hydrogen (HI) clumps. By studying the dependence of the shape of dense cold clouds formed by shock compression and cooling on the filling factor of HI clumps, we found that the density distribution of HI clumps determines the shape of molecular clouds formed by the jet-cloud interaction; arc-clouds are formed when the filling factor is large. On the other hand, when the filling factor is small, molecular clouds align with the jet. The jet propagates faster in models with small filling factors.

We carried out 2.5-dimensional resistive MHD simulations to study the formation mechanism of molecular loops observed by Fukui et al. (2006) at Galactic central region. Since it is hard to form molecular loops by uplifting dense molecular gas, we study the formation mechanism of molecular gas in rising magnetic arcades. This model is based on the in-situ formation model of solar prominences, in which prominences are formed by cooling instability in helical magnetic flux ropes formed by imposing converging and shearing motion at footpoints of the magnetic arch anchored to the solar surface. We extended this model to Galactic center scale (a few hundreds pc). Numerical results indicate that magnetic reconnection taking place in the current sheet formed inside the rising magnetic arcade creates dense blobs confined by the rising helical magnetic flux ropes. Thermal instability taking place in the flux ropes forms dense molecular filaments floating at high Galactic latitude. The mass of the filament increases with time, and can exceed 10^5 solar mass.

The vertical flow structure in the galactic center region remains poorly understood. We analyzed the MHD simulation data by Suzuki et al. (2015) for better understanding. As a result, we found the fast downflows with a speed of ~100 km/s near the foot-points of magnetic loops. These downflows are flowing along a magnetic field line and accelerated by the gravity. The direction of the fast flows is changed by the magnetic loop geometry, as it moves. As a result, not only vertical motions but also azimuthal and radial motions are excited. This feature could be relevant to the observed high velocity dispersion in the position-velocity diagram.

By performing a global magnetohydrodynamical (MHD) simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic centre (GC) region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternative to the bar-potential-driven model for the parallelogram shape of the central molecular zone. In addition, Parker instability (magnetic buoyancy) creates vertical magnetic structure, which would correspond to observed molecular loops, and frequently excited vertical flows. Furthermore, the time-averaged net gas flow is directed outward, whereas the flows are highly time dependent, which would contribute to the outflow from the bulge.

By performing a global magneto-hydrodynamical simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic center region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. In addition, the magnetic pressure-gradient force also drives radial flows in a similar manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternative to the bar-potential-driven model for the parallelogram-shape of the central molecular zone. This is a natural extension into the central few 100 pc of the magnetic activity, which is observed as molecular loops at radii from a few 100 pc to 1 kpc. Furthermore, the time-averaged net gas flow is directed outward, whereas the flows are highly time-dependent, which we discuss from a viewpoint of the outflow from the bulge.

The formation mechanism of the jet-aligned CO clouds found by NANTEN CO observations is studied by magnetohydrodynamical (MHD) simulations taking into account the cooling of the interstellar medium. Motivated by the association of the CO clouds with the enhancement of HI gas density, we carried out MHD simulations of the propagation of a supersonic jet injected into the dense HI gas. We found that the HI gas compressed by the bow shock ahead of the jet is cooled down by growth of the cooling instability triggered by the density enhancement. As a result, cold dense sheath is formed around the interface between the jet and the HI gas. The radial speed of the cold, dense gas in the sheath is a few km/s almost independent of the jet speed. Molecular clouds can be formed in this region. Since the dense sheath wrapping the jet reflects waves generated in the cocoon, the jet is strongly perturbed by the vortices of the warm gas in the cocoon, which breaks up the jet and forms a secondary shock in the HI-cavity drilled by the jet. The particle acceleration at the shock can be the origin of radio and X-ray filaments observed near the eastern edge of W50 nebula surrounding the galactic jet source SS433.

We have made new CO observations of two molecular clouds, which we call "jet" and "arc" clouds, toward the stellar cluster Westerlund 2 and the TeV gamma-ray source HESS J1023-575. The jet cloud shows a linear structure from the position of Westerlund 2 on the east. In addition, we have found a new counter jet cloud on the west. The arc cloud shows a crescent shape in the west of HESS J1023-575. A sign of star formation is found at the edge of the jet cloud and gives a constraint on the age of the jet cloud to be ~Myrs. An analysis with the multi CO transitions gives temperature as high as 20 K in a few places of the jet cloud, suggesting that some additional heating may be operating locally. The new TeV gamma-ray images by H.E.S.S. correspond to the jet and arc clouds spatially better than the giant molecular clouds associated with Westerlund 2. We suggest that the jet and arc clouds are not physically linked with Westerlund 2 but are located at a greater distance around 7.5 kpc. A microquasar with long-term activity may be able to offer a possible engine to form the jet and arc clouds and to produce the TeV gamma-rays, although none of the known microquasars have a Myr age or steady TeV gamma-rays. Alternatively, an anisotropic supernova explosion which occurred ~Myr ago may be able to form the jet and arc clouds, whereas the TeV gamma-ray emission requires a microquasar formed after the explosion.

We have discovered two molecular features at radial velocities of -35 km/s and 0 km/s toward the infrared Double Helix Nebula (DHN) in the Galactic center with NANTEN2. The two features show good spatial correspondence with the DHN. We have also found two elongated molecular ridges at these two velocities distributed vertically to the Galactic plane over 0.8 degree. The two ridges are linked by broad features in velocity and are likely connected physically with each other. The ratio between the 12CO J=2-1 and J=1-0 transitions is 0.8 in the ridges which is larger than the average value 0.5 in the foreground gas, suggesting the two ridges are in the Galactic center. An examination of the K band extinction reveals a good coincidence with the CO 0 km/s ridge and is consistent with a distance of 8 +/-2 kpc. We discuss the possibility that the DHN was created by a magnetic phenomenon incorporating torsional Alfv\'en waves launched from the circumnuclear disk (Morris, Uchida & Do 2006) and present a first estimate of the mass and energy involved in the DHN.

The tidal disruption event by a supermassive black hole in Swift J1644+57 can trigger limit-cycle oscillations between a supercritically accreting X-ray bright state and a subcritically accreting X-ray dim state. Time evolution of the debris gas around a black hole with mass $M=10^{6} {\MO}$ is studied by performing axisymmetric, two-dimensional radiation hydrodynamic simulations. We assumed the $\alpha$-prescription of viscosity, in which the viscous stress is proportional to the total pressure. The mass supply rate from the outer boundary is assumed to be ${\dot M}_{\rm supply}=100L_{\rm Edd}/c^2$, where $L_{\rm Edd}$ is the Eddington luminosity, and $c$ is the light speed. Since the mass accretion rate decreases inward by outflows driven by radiation pressure, the state transition from a supercritically accreting slim disk state to a subcritically accreting Shakura-Sunyaev disk starts from the inner disk and propagates outward in a timescale of a day. The sudden drop of the X-ray flux observed in Swift J1644+57 in August 2012 can be explained by this transition. As long as ${\dot M}_{\rm supply}$ exceeds the threshold for the existence of a radiation pressure dominant disk, accumulation of the accreting gas in the subcritically accreting region triggers the transition from a gas pressure dominant Shakura-Sunyaev disk to a slim disk. This transition takes place at $t {\sim}~50/({\alpha}/0.1)$ days after the X-ray darkening. We expect that if $\alpha > 0.01$, X-ray emission with luminosity $\gtrsim 10^{44}$ ${\rm erg}{\cdot}{\rm s}^{-1}$ and jet ejection will revive in Swift J1644+57 in 2013--2014.

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.

We present global solutions of optically thin, two-temperature black hole accretion disks incorporating magnetic fields. We assume that the {\pi}{\phi}-component of the Maxwell stress is proportional to the total pressure, and prescribe the radial dependence of the magnetic flux advection rate in order to complete the set of basic equations. We obtained magnetically supported (low-{\beta}) disk solutions, whose luminosity exceeds the maximum luminosity for an advection-dominated accretion flow (ADAF), L > 0.4 {\alpha}^2 L_Edd, where L_Edd is the Eddington luminosity. The accretion flow is composed of the outer ADAF, a luminous hot accretion flow (LHAF) inside the transition layer from the outer ADAF to the low-{\beta} disk, the low-{\beta} disk, and the inner ADAF. The low-{\beta} disk region becomes wider as the mass-accretion rate increases further. In the low-{\beta} disk, the magnetic heating balances the radiative cooling, and the electron temperature decreases from ~ 10^9.5 K to ~ 10^8 K as the luminosity increases. These results are consistent with the anti-correlation between the energy cutoff in X-ray spectra (hence the electron temperature) and the luminosity when L > 0.1 L_Edd, observed in the bright/hard state during the bright hard-to-soft transitions of transient outbursts in galactic black hole candidates.

We obtained self-similar solutions of relativistically expanding magnetic loops by assuming axisymmetry and a purely radial flow. The stellar rotation and the magnetic fields in the ambient plasma are neglected. We include the Newtonian gravity of the central star. These solutions are extended from those in our previous work (Takahashi, Asano, & Matsumoto 2009) by taking into account discontinuities such as the contact discontinuity and the shock. The global plasma flow consists of three regions, the outflowing region, the post shocked region, and the ambient plasma. They are divided by two discontinuities. The solutions are characterized by the radial velocity, which plays a role of the self-similar parameter in our solutions. The shock Lorentz factor gradually increases with radius. It can be approximately represented by the power of radius with the power law index of 0.25. We also carried out magnetohydrodynamic simulations of the evolution of magnetic loops to study the stability and the generality of our analytical solutions. We used the analytical solutions as the initial condition and the inner boundary conditions. We confirmed that our solutions are stable over the simulation time and that numerical results nicely recover the analytical solutions. We then carried out numerical simulations to study the generality of our solutions by changing the power law index \delta of the ambient plasma density \rho_0 \propto r^{-\delta}. We alter the power law index \delta from 3.5 in the analytical solutions. The analytical solutions are used as the initial conditions inside the shock in all simulations. We observed that the shock Lorentz factor increases with time when \delta is larger than 3, while it decreases with time when \delta is smaller than 3. The shock Lorentz factor is proportional to t^{(\delta-3)/2}. These results are consistent with the analytical studies by Shapiro (1979).

We have carried out 12CO(J =2-1) and 12CO(J =3-2) observations at spatial resolutions of 1.0-3.8 pc toward the entirety of loops 1 and 2 and part of loop 3 in the Galactic center with NANTEN2 and ASTE. These new results revealed detailed distributions of the molecular gas and the line intensity ratio of the two transitions, R3-2/2-1. In the three loops, R3-2/2-1 is in a range from 0.1 to 2.5 with a peak at ~ 0.7 while that in the disk molecular gas is in a range from 0.1 to 1.2 with a peak at 0.4. This supports that the loops are more highly excited than the disk molecular gas. An LVG analysis of three transitions, 12CO J =3-2 and 2-1 and 13CO J =2-1, toward six positions in loops 1 and 2 shows density and temperature are in a range 102.2 - 104.7 cm-3 and 15-100 K or higher, respectively. Three regions extended by 50-100 pc in the loops tend to have higher excitation conditions as characterized by R3-2/2-1 greater than 1.2. The highest ratio of 2.5 is found in the most developed foot points between loops 1 and 2. This is interpreted that the foot points indicate strongly shocked conditions as inferred from their large linewidths of 50-100 km s-1, confirming the suggestion by Torii et al. (2010b). The other two regions outside the foot points suggest that the molecular gas is heated up by some additional heating mechanisms possibly including magnetic reconnection. A detailed analysis of four foot points have shown a U shape, an L shape or a mirrored-L shape in the b-v distribution. It is shown that a simple kinematical model which incorporates global rotation and expansion of the loops is able to explain these characteristic shapes.

Fukui et al. (2006) discovered two huge molecular loops in the Galactic center located in (l, b) ~ (355 deg-359 deg, 0 deg-2 deg) in a large velocity range of -180-40 km s^-1. Following the discovery, we present detailed observational properties of the two loops based on NANTEN 12CO(J=1-0) and 13CO(J=1-0) datasets at 10 pc resolution including a complete set of velocity channel distributions and comparisons with HI and dust emissions as well as with the other broad molecular features. We find new features on smaller scales in the loops including helical distributions in the loop tops and vertical spurs. The loops have counterparts of the HI gas indicating that the loops include atomic gas. The IRAS far infrared emission is also associated with the loops and was used to derive an X-factor of 0.7(+/-0.1){\times}10^20 cm^-2 (K km s^-1)^-1 to convert the 12CO intensity into the total molecular hydrogen column density. From the 12CO, 13CO, H I and dust datasets we estimated the total mass of loops 1 and 2 to be ~1.4 {\times} 106 Msun and ~1.9 {\times} 10^6 Msun, respectively, where the H I mass corresponds to ~10-20% of the total mass and the total kinetic energy of the two loops to be ~10^52 ergs. An analysis of the kinematics of the loops yields that the loops are rotating at ~47 km s-1 and expanding at ~141 km s^-1 at a radius of 670 pc from the center. Fukui et al. (2006) presented a model that the loops are created by the magnetic flotation due to the Parker instability with an estimated magnetic field strength of ~150 {\mu}G. We present comparisons with the recent numerical simulations of the magnetized nuclear disk by Machida et al. (2009) and Takahashi et al. (2009) and show that the theoretical results are in good agreements with the observations. The helical distributions also suggest that some magnetic instability plays a role similarly to the solar helical features.

Fukui et al. (2006) discovered two molecular loops in the Galactic center and argued that the foot points of the molecular loops, two bright spots at both loops ends, represent the gas accumulated by the falling motion along the loops, subsequent to magnetic flotation by the Parker instability. We have carried out sensitive CO observations of the foot points toward l=356 deg at a few pc resolution in the six rotational transitions of CO; 12CO(J=1-0, 3-2, 4-3, 7-6), 13CO(J=1-0) and C18O(J=1-0). The high resolution image of 12CO (J=3-2) has revealed the detailed distribution of the high excitation gas including U shapes, the outer boundary of which shows sharp intensity jumps accompanying strong velocity gradients. An analysis of the multi-J CO transitions shows that the temperature is in a range from 30-100 K and density is around 10^3-10^4 cm^-3, confirming that the foot points have high temperature and density although there is no prominent radiative heating source such as high mass stars in or around the loops. We argue that the high temperature is likely due to the shock heating under C-shock condition caused by the magnetic flotation. We made a comparison of the gas distribution with theoretical numerical simulations and note that the U shape is consistent with numerical simulations. We also find that the region of highest temperature of ~100 K or higher inside the U shape corresponds to the spur having an upward flow, additionally heated up either by magnetic reconnection or bouncing in the interaction with the narrow neck at the bottom of the U shape. We note these new findings further reinforce the magnetic floatation interpretation.

We obtained thermal equilibrium solutions for optically thin, two-temperature black hole accretion disks incorporating magnetic fields. The main objective of this study is to explain the bright/hard state observed during the bright/slow transition of galactic black hole candidates. We assume that the energy transfer from ions to electrons occurs via Coulomb collisions. Bremsstrahlung, synchrotron, and inverse Compton scattering are considered as the radiative cooling processes. In order to complete the set of basic equations, we specify the magnetic flux advection rate. We find magnetically supported (low-beta), thermally stable solutions. In these solutions, the total amount of the heating via the dissipation of turbulent magnetic fields goes into electrons and balances the radiative cooling. The low-$\beta$ solutions extend to high mass accretion rates and the electron temperature is moderately cool. High luminosities and moderately high energy cutoffs in the X-ray spectrum observed in the bright/hard state can be explained by the low-beta solutions.

We determine the spin of a supermassive black hole in the context of discseismology by comparing newly detected quasi-periodic oscillations (QPOs) of radio emission in the Galactic centre, Sagittarius A* (Sgr A*), as well as infrared and X-ray emissions with those of the Galactic black holes. We find that the spin parameters of black holes in Sgr A* and in Galactic X-ray sources have a unique value of $\approx 0.44$ which is smaller than the generally accepted value for supermassive black holes, suggesting evidence for the angular momentum extraction of black holes during the growth of supermassive black holes. Our results demonstrate that the spin parameter approaches the equilibrium value where spin-up via accretion is balanced by spin-down via the Blandford-Znajek mechanism regardless of its initial spin. We anticipate that measuring the spin of black holes by using QPOs will open a new window for exploring the evolution of black holes in the Universe.

We have discovered a molecular dome-like feature towards $355^{\circ} \leq l \leq 359^{\circ}$ and $0^{\circ} \leq b \leq 2^{\circ}$. The large velocity dispersions of 50--100 km s$^{-1}$ of this feature are much larger than those in the Galactic disk and indicate that the feature is located in the Galactic center, probably within $\sim1$ kpc of Sgr A$^{*}$. The distribution has a projected length of $\sim600$ pc and height of $\sim300$ pc from the Galactic disk and shows a large-scale monotonic velocity gradient of $\sim130$ km s $^{-1}$ per $\sim600$ pc. The feature is also associated with HI gas having a more continuous spatial and velocity distribution than that of $^{12}$CO. We interpret the feature as a magnetically floated loop similar to loops 1 and 2 and name it "loop 3". Loop 3 is similar to loops 1 and 2 in its height and length but is different from loops 1 and 2 in that the inner part of loop 3 is filled with molecular emission. We have identified two foot points at the both ends of loop 3. HI, $^{12}$CO and $^{13}$CO datasets were used to estimate the total mass and kinetic energy of loop 3 to be $\sim3.0 \times 10^{6} \Mo$ and $\sim1.7 \times 10^{52}$ ergs. The huge size, velocity dispersions and energy are consistent with the magnetic origin the Parker instability as in case of loops 1 and 2 but is difficult to be explained by multiple stellar explosions. We argue that loop 3 is in an earlier evolutionary phase than loops 1 and 2 based on the inner-filled morphology and the relative weakness of the foot points. This discovery indicates that the western part of the nuclear gas disk of $\sim1$ kpc radius is dominated by the three well-developed magnetically floated loops and suggests that the dynamics of the nuclear gas disk is strongly affected by the magnetic instabilities.

We carried out two-dimensional magnetohydrodynamic simulations of the Galactic gas disk to show that the dense loop-like structures discovered by the Galactic center molecular cloud survey by NANTEN 4 m telescope can be formed by the buoyant rise of magnetic loops due to the Parker instability. At the initial state, we assumed a gravitationally stratified disk consisting of the cool layer ($T \sim 10^3$ K), warm layer ($T \sim 10^4$ K), and hot layer ($T \sim 10^5$ K). Simulation box is a local part of the disk containing the equatorial plane. The gravitational field is approximated by that of a point mass at the galactic center. The self-gravity, and the effects of the galactic rotation are ignored. Numerical results indicate that the length of the magnetic loops emerging from the disk is determined by the scale height of the hot layer ($\sim$ 100 pc at 1 kpc from the Galactic center). The loop length, velocity gradient along the loops and large velocity dispersions at their foot points are consistent with the NANTEN observations. We also show that the loops become top-heavy when the curvature of the loop is sufficiently small, so that the rising loop accumulates the overlying gas faster than sliding it down along the loop. This mechanism is similar to that of the formation of solar chromospheric arch filaments. The molecular loops emerge from the low temperature layer just like the dark filaments observed in the H$\alpha$ image of the emerging flux region of the sun.

We present the results of global three dimensional magneto-hydrodynamic simulations of black hole accretion flows. We focus on the dependence of numerical results on the gas temperature Tout supplied from the outer region. General relativistic effects are taken into account using the pseudo-Newtonian potential. We ignore the radiative cooling of the accreting gas. The initial state is a torus whose density maximum is at 35rs or 50rs from the gravitating center, where rs is the Schwarzschild radius. The torus is initially threaded by a weak azimuthal magnetic field. We found that mass accretion rate and the mass outflow rate strongly depend on the temperature of the initial torus. The ratio of the average Maxwell stress generated by the magneto-rotational instability (MRI) to gas pressure is alpha ~0.05 in the hot torus and alpha ~ 0.01 in the cool torus. In the cool model, a constant angular momentum inner torus is formed around 4-8rs. This inner torus deforms itself from a circle to a crescent quasi-periodically. During this deformation, the mass accretion rate, the magnetic energy and the Maxwell stress increase. As the magnetic energy is released, the inner torus returns to a circular shape and starts the next cycle. Power spectral density (PSD) of the time variation of the mass accretion rate in the cool model has a low frequency peak around 10Hz when we assumed a 10 solar mass black hole. The PSD of the hot model is flat in 1-30Hz. The slope of the PSD in the cool model is steeper than that in the hot model in 30-100Hz. The mass outflow rate in the low temperature model also shows quasi-periodic oscillation. Intermittent outflows are created in both models. The outflow speed is $0.01c - 0.05c$.

Supercritical accretion flows inevitably produce radiation-pressure driven outflows, which will Compton up-scatter soft photons from the underlying accretion flow, thereby making hard emission. We perform two dimensional radiation hydrodynamic simulations of supercritical accretion flows and outflows, incorporating such Compton scattering effects, and demonstrate that there appears a new hard spectral state at higher photon luminosities than that of the slim-disk state. In this state, as the photon luminosity increases, the photon index decreases and the fraction of the hard emission increases. The Compton $y$-parameter is of the order of unity (and thus the photon index will be $\sim 2$) when the apparent photon luminosity is ${\sim}30L_{\rm E}$ (with $L_{\rm E}$ being the Eddington luminosity) for nearly face-on sources. This explains the observed spectral hardening of the ULX NGC1313 X-2 in its brightening phase and thus supports the model of supercritical accretion onto stellar mass black holes in this ULX.

A survey for the molecular clouds in the Galaxy with NANTEN mm telescope has discovered molecular loops in the Galactic center region. The loops show monotonic gradients of the line of sight velocity along the loops and the large velocity dispersions towards their foot points. It is suggested that these loops are explained in terms of the buoyant rise of magnetic loops due to the Parker instability. We have carried out global three-dimensional magneto-hydrodynamic simulations of the gas disk in the Galactic center. The gravitational potential is approximated by the axisymmetric potential proposed by Miyamoto & Nagai (1975). At the initial state, we assume a warm (~ 10^4 K) gas torus threaded by azimuthal magnetic fields. Self-gravity and radiative cooling of the gas are ignored. We found that buoyantly rising magnetic loops are formed above the differentially rotating, magnetically turbulent disk. By analyzing the results of global MHD simulations, we have identified individual loops, about 180 in the upper half of the disk, and studied their statistical properties such as their length, width, height, and velocity distributions along the loops. Typical length and height of a loop are 1kpc and 200pc, respectively. The line of sight velocity changes linearly along a loop and shows large dispersions around the foot-points. Numerical results indicate that loops emerge preferentially from the region where magnetic pressure is large. We argue that these properties are consistent with those of the molecular loops discovered by NANTEN.

We present new thermal equilibrium solutions for optically thin and thick disks incorporating magnetic fields. The purpose of this paper is to explain the bright hard state and the bright/slow transition observed in the rising phases of outbursts in BHCs. On the basis of the results of 3D MHD simulations, we assume that magnetic fields inside the disk are turbulent and dominated by the azimuthal component and that the azimuthally averaged Maxwell stress is proportional to the total pressure. We prescribe the magnetic flux advection rate to determine the azimuthal magnetic flux at a given radius. We find magnetically supported, thermally stable solutions for both optically thin and thick disks, in which the heating enhanced by the strong magnetic field balances the radiative cooling. The temperature in a low-$\beta$ disk is lower than that in an ADAF/RIAF but higher than that in a standard disk. We also study the radial dependence of the thermal equilibrium solutions. The optically thin, low-$\beta$ branch extends to $ \dot M \gtrsim 0.1 {\dot M}_{\rm Edd}$, in which the temperature anti-correlates with the mass accretion rate. Thus optically thin low-$\beta$ disks can explain the bright hard state. Optically thick, low-$\beta$ disks have the radial dependence of the effective temperature $T_{\rm eff} \propto \varpi^{-3/4}$. Such disks will be observed as staying in a high/soft state. Furthermore, limit cycle oscillations between an optically thick low-$\beta$ disk and a slim disk will occur because the optically thick low-$\beta$ branch intersects with the radiation pressure dominated standard disk branch. These limit cycle oscillations will show a smaller luminosity variation than that between a standard disk and a slim disk.

We obtained self-similar solutions of relativistically expanding magnetic loops taking into account the azimuthal magnetic fields. We neglect stellar rotation and assume axisymmetry and a purely radial flow. As the magnetic loops expand, the initial dipole magnetic field is stretched into the radial direction. When the expansion speed approaches the light speed, the displacement current reduces the toroidal current and modifies the distribution of the plasma lifted up from the central star. Since these self-similar solutions describe the free expansion of the magnetic loops, i.e., $Dv/Dt=0$, the equations of motion are similar to those of the static relativistic magnetohydrodynamics. This allows us to estimate the total energy stored in the magnetic loops by applying the virial theorem. This energy is comparable to that of the giant flares observed in magnetars.

Steep gradients of temperature and density, called cold fronts, are observed by Chandra in a leading edge of subclusters moving through the intracluster medium (ICM). The presence of cold fronts indicates that thermal conduction across the front is suppressed by magnetic fields. We carried out three-dimensional magnetohydrodynamic (MHD) simulations including anisotropic thermal conduction of a subcluster moving through a magnetically turbulent ICM. We found that turbulent magnetic fields are stretched and amplified by shear flows along the interface between the subcluster and the ambient ICM. Since magnetic fields reduce the efficiency of thermal conduction across the front, the cold front survives at least 1 Gyr. We also found that a moving subcluster works as an amplifier of magnetic fields. Numerical results indicate that stretched turbulent magnetic fields accumulate behind the subcluster and are further amplified by vortex motions. The moving subcluster creates a long tail of ordered magnetic fields, in which the magnetic field strength attains plasma beta < 10.

We obtained steady solutions of optically thin, single temperature, magnetized black hole accretion disks assuming thermal bremsstrahlung cooling. Based on the results of 3D MHD simulations of accretion disks, we assumed that the magnetic fields inside the disk are turbulent and dominated by azimuthal component. We decomposed magnetic fields into an azimuthally averaged mean field and fluctuating fields. We also assumed that the azimuthally averaged Maxwell stress is proportional to the total pressure. The radial advection rate of the azimuthal magnetic flux $\dot \Phi$ is prescribed as being proportional to $\varpi^{- \zeta}$, where $\varpi$ is the radial coordinate and $\zeta$ is a parameter which parameterizes the radial variation of $\dot \Phi$. We found that when accretion rate $\dot M$ exceeds the threshold for the onset of the thermal instability, a magnetic pressure dominated new branch appears. Thus the thermal equilibrium curve of optically thin disk has a 'Z'-shape in the plane of surface density and temperature. This indicates that as the mass accretion rate increases, a gas pressure dominated optically thin hot accretion disk undergoes a transition to a magnetic pressure dominated, optically thin cool disk. This disk corresponds to the X-ray hard, luminous disk in black hole candidates observed during the transition from a low/hard state to a high/soft state. We also obtained global steady transonic solutions containing such a transition layer.

We carried out global three-dimensional resistive magnetohydrodynamic simulations of galactic gaseous disks to investigate how the galactic magnetic fields are amplified and maintained. We adopt a steady axisymmetric gravitational potential given by \cite*{Miy80}. As the initial condition, we assume a warm ($T\sim10^{5} \mathrm{K}$) rotating gas torus centered at $\varpi=10 \mathrm{kpc}$ threaded by weak azimuthal magnetic fields. Numerical results indicate that in differentially rotating galactic gaseous disks, magnetic fields are amplified due to magneto-rotational instability and magnetic turbulence develops. After the amplification of magnetic energy saturates, the disk stays in a quasi-steady state. The mean azimuthal magnetic field increases with time and shows reversals with period of 1Gyr (2Gyr for a full cycle). The amplitude of $B_{\phi}$ near the equatorial plane is $B_{\phi}\sim1.5\mu \mathrm{G}$ at $\varpi=5 \mathrm{kpc}$. The magnetic fields show large fluctuations whose standard deviation is comparable to the mean field. The mean azimuthal magnetic field in the disk corona has direction opposite to the mean magnetic field inside the disk. The mass accretion rate driven by the Maxwell stress is $\sim10^{-3}M_{\odot}/ \mathrm{yr}$ at $\varpi=2.5 \mathrm{kpc}$ when the mass of the initial torus is $\sim5\times10^{8}M_{\odot}$.

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

Recent high resolution observations by the {\it Chandra} X-ray satellite revealed various substructures in hot X-ray emitting plasmas in cluster of galaxies. For example, {\it Chandra} revealed the existence of sharp discontinuities in the surface brightness at the leading edge of subclusters in merging clusters (e.g., Abell 3667), where the temperature drops sharply across the fronts. These sharp edges are called cold fronts. We present results of three-dimensional (3D) magnetohydrodynamic simulations of the interaction between a dense subcluster plasma and ambient magnetized intracluster medium. Anisotropic heat conduction along magnetic field lines is included. At the initial state, magnetic fields are assumed to be uniform and transverse to the motion of the dense subcluster. Since magnetic fields ahead of the subcluster slip toward the third direction in the 3D case, the strength of magnetic fields in this region can be reduced compared to that in the 2D case. Nevertheless, a cold front can be maintained because the magnetic field lines wrapping around the forehead of the subcluster suppress the heat conduction across them. On the other hand, when the magnetic field is absent, a cold front cannot be maintained because isotropic heat conduction from the hot ambient plasma rapidly heats the cold subcluster plasma.

We study the magnetic interaction between a neutron star and its disk by solving the time-dependent relativistic force-free equations. At the initial state, we assume that the dipole magnetic field of the neutron star connects the neutron star and its equatorial disk, which deeply enters into the magnetosphere of the neutron star. Magnetic fields are assumed to be frozen to the star and the disk. The rotation of the neutron star and the disk is imposed as boundary conditions. We apply Harten-Lax-van Leer (HLL) method to simulate the evolution of the star-disk system. We carry out simulations for (1) a disk inside the corotation radius, in which the disk rotates faster than the star, and (2) a disk outside the corotation radius, in which the neutron star rotates faster than the disk. Numerical results indicate that for both models, the magnetic field lines connecting the disk and the star inflate as they are twisted by the differential rotation between the disk and the star. When the twist angle exceeds pi radian, the magnetic field lines expand with speed close to the light speed. This mechanism can be the origin of relativistic outflows observed in binaries containing a neutron star.

High resolution observations of cluster of galaxies by Chandra have revealed the existence of an X-ray emitting comet-like galaxy C153 in the core of cluster of galaxies A2125. The galaxy C153 moving fast in the cluster core has a distinct X-ray tail on one side, obviously due to ram pressure stripping, since the galaxy C153 crossed the central region of A2125. The X-ray emitting plasma in the tail is substantially cooler than the ambient plasma. We present results of two-dimensional magnetohydrodynamic simulations of the time evolution of a subclump like C153 moving in magnetized intergalactic matter. Anisotropic heat conduction is included. We found that the magnetic fields are essential for the existence of the cool X-ray tail, because in non-magnetized plasma the cooler subclump tail is heated up by isotropic heat conduction from the hot ambient plasma and does not form such a comet-like tail.

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.

Recent Chandra observations of clusters of galaxies revealed the existence of a sharp ridge in the X-ray surface brightness where the temperature drops across the front. This front is called the cold front. We present the results of two-dimensional magnetohydrodynamic simulations of the time evolution of a dense subcluster plasma moving in a cluster of galaxies. Anisotropic heat conduction along the magnetic field lines is included. In the models without magnetic fields, the numerical results indicate that the heat conduction from the hot ambient plasma heats the cold dense plasma of the subcluster and diffuses out the cold front. When magnetic fields exist in a cluster of galaxies, however, cold fronts can be maintained because the heat conduction across the magnetic field lines is suppressed. We found that, even when the magnetic fields in a cluster of galaxies are disordered, heat conduction across the front is restricted because the magnetic field lines are stretched along the front. Numerical results reproduced the X-ray intensity distribution observed in the A3667 cluster of galaxies.

We present the results of 2.5-dimensional resistive magnetohydrodynamic (MHD) simulations of the magnetic interaction between a weakly magnetized neutron star and its accretion disk. General relativistic effects are simulated by using the pseudo-Newtonian potential. We find that well-collimated jets traveling along the rotation axis of the disk are formed by the following mechanism: (1) The magnetic loops connecting the neutron star and the disk are twisted due to the differential rotation between the neutron star and the disk. (2) Twist injection from the disk initiates expansion of the loop. (3) The expanding magnetic loops create a magnetic tower in which accelerated disk material travel as collimated bipolar jets. The propagation speed of the working surface of the jet is the order of 10% of the speed of light ($\sim 0.1c$). (4) Magnetic reconnection taking place inside the expanding magnetic loops injects hot bubbles intermittently into the magnetic tower. The ejection speed of the bubble is the order of the local Alfv\'{e}n speed of the launching point and $\sim 0.2c$ in our simulations. (5) The hot bubbles moving inside the tower catch up with the working surface of the jet. High energy electrons created by the magnetic reconnection are a plausible source of radio emission. Our model can explain the formation process of a narrow jet from a weakly magnetized ($|{\boldmath$B_{*}$}|\le 10^{9}$ gauss) neutron star and the correlation between radio flares of the core and of the lobe observed in Sco X-1.

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

The detectability of gravitational wave (GW) radiation from accretion disks is discussed based on various astrophysical contexts. In order to emit GW radiation, the disk shape should lose axial symmetry. We point out that a significant deformation is plausible in non-radiative hot accretion disks because of enhanced magnetic activity, whereas it is unlikely for standard-type cool disks. We have analyzed the 3D magnetohydrodynamical (MHD) simulation data of magnetized accretion flow, finding non-axisymmetric density patterns. The corresponding ellipticity is $\epsilon \sim 0.01$. The expected time variations of GW radiation are overall chaotic, but there is a hint of quasi-periodicity. GW radiation has no interesting consequence, however, in the case of close binaries, because of very tiny disk masses. GW radiation is not significant, either, for AGN because of very slow rotation velocities. The most promising case can be found in gamma-ray bursts or supernovae, in which a massive torus (or disk) with a solar mass or so may be formed around a stellar-mass compact object as the result of a merger of compact objects, or by the fallback of exploded material towards the center in a supernova. Although much more intense GW radiation is expected before the formation of the torus, the detection of GW radiation in the subsequent accretion phase is of great importance, since it will provide a good probe to investigating their central engines.

We carried out two dimensional high-resolution magnetohydrodynamic (MHD) simulations of an accretion disk around a weakly magnetized neutron star. General relativistic effects are taken into account by using pseudo-Newtonian potential of Pacz\'{y}nski, B., P. J. Witta (1980). When magnetic loops connect the neutron star and the accretion disk, the twist injection from the disk or from the rotating neutron star triggers expansion of the loops. Since the expanding magnetic loops prevent inflow toward the magnetic poles of the neutron star, disk matter accumulates on the boundary between the magnetosphere and the disk. Magnetic reconnection taking place in the loops creates a channel along which the disk matter can accrete and unloads the magnetosphere. This process produces quasi-periodic variation of the accretion flow in the innermost region of the disk. We found two kinds of oscillations. One is the magnetospheric oscillation regulated by magnetic reconnection. The other is the radial disk oscillation. The typical frequency of the oscillations is 100 Hz to 2 kHz. Furthermore, we predict that QPO sources inevitably accompany X-ray flares by magnetic reconnection and bipolar outflows of hot X-ray emitting plasma similar to the optical jets in protostars.

We carried out 2.5-dimensional resistive magnetohydrodynamic simulations to study the effects of magnetic diffusivity on magnetically driven mass accretion and jet formation. We found that (1) when the normalized magnetic diffusivity, is small, mass accretion and jet formation take place intermittently; (2) when diffusivity is middle, the system evolves toward a quasi-steady state; and the system evolves toward a quasi-steady state; and (3) when diffusivity is large, the accretion/mass outflow rate decreases with diffusivity and approaches 0. The results of these simulations indicate magnetic braking provide a mass accretion rate which is sufficient to explain the activity of AGNs.

Outflows and convective motions in accretion flows have been intensively discussed recently in the context of advection-dominated accretion flow (ADAF) based on two-dimensional (2D) and three-dimensional (3D) hydrodynamical simulations. We, however, point that without proper treatments of the disk magnetic fields, a major source of viscosity, one can never derive general, firm conclusions concerning the occurrence of outflows and convection. We analyzed the 3D MHD numerical simulation data of magnetized accretion flows initially threaded by weak toroidal magnetic fields, finding large-scale convective motions dominating near the black hole. In contrast, outflows occur only temporarily and are not very significant in our simulations. If there grow strong vertical fields somehow, however, formation of bi-polar jets is inevitable. It is claimed that radiation could be dominant at the outermost zones of the convective disks because of outward energy flow by convection, however, this is no longer the case in convective MHD flows, since accretion energy can be released via magnetic reconnection in the inner parts. Such reconnection leads to sporadic flare events, thus producing substantial variability in out-going radiation.

The stability of nonaxisymmetric perturbations in differentially rotating astrophysical accretion disks is analyzed by fully incorporating the properties of shear flows. We verify the presence of discrete unstable eigenmodes with complex and pure imaginary eigenvalues, without any artificial disk edge boundaries, unlike Ogilvie & Pringle(1996)'s claim. By developing the mathematical theory of a non-self-adjoint system, we investigate the nonlocal behavior of eigenmodes in the vicinity of Alfven singularities at omega_D=omega_A, where omega_D is the Doppler-shifted wave frequency and omega_A=k_// v_A is the Alfven frequency. The structure of the spectrum of discrete eigenmodes is discussed and the magnetic field and wavenumber dependence of the growth rate are obtained. Exponentially growing modes are present even in a region where the local dispersion relation theory claims to have stable eigenvalues. The velocity field created by an eigenmode is obtained, which explains the anomalous angular momentum transport in the nonlinear stage of this stability.

Rapid fluctuation with a frequency dependence of $1/f^{\alpha}$ (with $\alpha \simeq 1 - 2$) is characteristic of radiation from black-hole objects. Its origin remains poorly understood. We examine the three-dimensional magnetohydrodynamical (MHD) simulation data, finding that a magnetized accretion disk exhibits both $1/f^\alpha$ fluctuation (with $\alpha \simeq 2$) and a fractal magnetic structure (with the fractal dimension of $D \sim 1.9$). The fractal field configuration leads reconnection events with a variety of released energy and of duration, thereby producing $1/f^\alpha$ fluctuations.

We present the results of three-dimensional global magnetohydrodynamic (MHD) simulations of the Parker-shearing instability in a differentially rotating torus initially threaded by toroidal magnetic fields. An equilibrium model of magnetized torus is adopted as an initial condition. When $\beta_0 = P_{\rm gas}/P_{\rm mag} \sim 1$ at the initial state, magnetic flux buoyantly escapes from the disk and creates loop-like structures similar to those in the solar corona. Inside the torus, growth of non-axisymmetric magneto-rotational (or Balbus & Hawley) instability generates magnetic turbulence. Magnetic field lines are tangled in small scale but in large scale they show low azimuthal wave number spiral structure. After several rotation period, the system oscillates around a state with $\beta \sim 5$. We found that magnetic pressure dominated ($\beta < 1$) filaments are created in the torus. The volume filling factor of the region where $\beta \leq 0.3$ is 2-10%. Magnetic energy release in such low-$\beta$ regions may lead to violent flaring activities in accretion disks and in galactic gas disks.

We propose a model of hard X-ray flares in protostars observed by ASCA satellite. Assuming that the dipole magnetic field of the protostar threads the protostellar disk, we carried out 2.5-dimensional magnetohydrodynamic (MHD) simulations of the disk-star interaction. The closed magnetic loops connecting the central star and the disk are twisted by the rotation of the disk. As the twist accumulates, magnetic loops expand and finally approach to the open field configuration. A current sheet is formed inside the expanding loops. In the presence of resistivity, magnetic reconnection takes place in the current sheet. Outgoing magnetic island and post flare loops are formed as a result of the reconnection. The time scale of this `flare' is the order of the rotation period of the disk. The released magnetic energy partly goes into the thermal energy and heats up the flaring plasma up to $10^8$ K. The length of the flaring loop is several times of the radius of the central star, consistent with observations. The speed of the hot plasmoid ejected by the reconnection is $200-400 $ km s$^{-1}$ when the footpoint of the loop is at 0.03 AU from 1 M$_\odot$ protostar. The hot plasma outflow can explain the speed and mass flow rate of optical jets. Dense, cold, magnetically accelerated wind ($v \sim 150-250$ km s$^{-1}$) emanates from the surface of the disk along the partially open magnetic field lines threading the disk. This dense, cold wind may correspond to high velocity neutral winds.

The response of advection-dominated accretion disks to local disturbances is examined by one-dimensional numerical simulations. It is generally believed that advection-dominated disks are thermally stable. We, however, find that any disurbance added onto accretion flow at large radii does not decay so rapidly that it can move inward with roughly the free-fall velocity. Although disturbances continue to be present, the global disk structure will not be modified largely. This can account for persistent hard X-ray emission with substantial variations observed in active galactic nuclei and stellar black hole candidates during the hard state. Moreover, when the disturbance reaches the innermost parts, an acoustic wave emerges, propagating outward as a shock wave. The resultant light variation is roughly (time) symmetric and is quite reminiscent of the observed X-ray shots of Cygnus X-1.

The effect of the magnetic skew on the Parker instability is investigated by means of the linear stability analysis for a gravitationally stratified gas layer permeated by a horizontal magnetic field. When the magnetic field is skewed (i.e., the field line direction is a function of the height), the wavelength of the most unstable mode is $ \lambda \; \sim \; 10 H $ where $ H $ is the pressure scale height. The growth rate of the short wavelength modes is greatly reduced when the gradient in the magnetic field direction exceeds 0.5 radian per scale height. Our results indicate that the Parker instability in a skewed magnetic field preferentially forms large scale structures like giant molecular clouds.