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.

Using two-dimensional magnetohydrodynamic simulations, we investigated the gas dynamics around a black hole plunging into a molecular cloud. In these calculations, we assumed a parallel-magnetic-field layer in the cloud. The size of the accelerated region is far larger than the Bondi-Hoyle-Lyttleton radius, being approximately inversely proportional to the Alfv\'en Mach number for the plunging black hole. Our results successfully reproduce the "Y" shape in position velocity maps of the "Bullet" in the W44 molecular cloud. The size of the Bullet is also reproduced within an order of magnitude using a reasonable parameter set. This consistency supports the shooting model of the Bullet, according to which an isolated black hole plunged into a molecular cloud to form a compact broad-velocity-width feature.

We perform global three-dimensional (3D) radiation-hydrodynamic (RHD) simulations of out- flow from supercritical accretion flow around a 10 Msun black hole. We only solve the outflow part, starting from the axisymmetric 2D simulation data in a nearly steady state but with small perturbations in a sinusoidal form being added in the azimuthal direction. The mass accretion rate onto the black hole is ~10^2 L_E/c^2 in the underlying 2D simulation data and the outflow rate is ~10 L_E/c^2 (with LE and c being the Eddington luminosity and speed of light, respectively). We first confirm the emergence of clumpy outflow, which was discovered by the 2D RHD simulations, above the photosphere located at a few hundreds of Schwarzschild radii (r_S) from the central black hole. As prominent 3D features we find that the clumps have the shape of a torn sheet, rather than a cut string, and that they are rotating around the central black hole with a sub-Keplerian velocity at a distance of ~10^3 r_S from the center. The typical clump size is ~30 r_S or less in the radial direction, and is more elongated in the angular directions, ~hundreds of r_S at most. The sheet separation ranges from 50 to 150 r_S. We expect stochastic time variations when clumps pass across the line of the sight of a distant observer. Variation timescales are estimated to be several seconds for a black hole with mass of ten to several tens of Msun, in rough agreement with the observations of some ultra-luminous X-ray sources.

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.

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.