The formation and early evolution of low mass young stellar objects (YSOs) are investigated using three-dimensional non-ideal magneto-hydrodynamics simulations. We investigate the evolution of YSOs up to ~ 10^4 yr after protostar formation, at which protostellar mass reaches ~ 0.1 M_\odot . We particularly focus on the impact of the dust model on the evolution. We found that a circumstellar disk is formed in all simulations regardless of the dust model. Disk size is approximately 10 AU at the protostar formation epoch, and it increases to several tens of AU at ~ 10^4 yr after protostar formation. Disk mass is comparable to central protostellar mass and gravitational instability develops. In the simulations with small dust size, the warp of the pseudodisk develops ~ 10^4 yr after protostar formation. The warp strengthens magnetic braking in the disk and decreases disk size. Ion-neutral drift can occur in the infalling envelope under the conditions that the typical dust size is a \gtrsim 0.2{\mu}m and the protostar (plus disk) mass is M \gtrsim 0.1 M_\odot. The outflow activity is anti-correlated to the dust size and the strong outflow appears with small dust grains.

The formation of circumstellar disks is investigated using three-dimensional resistive magnetohydrodynamic simulations, in which the initial prestellar cloud has a misaligned rotation axis with respect to the magnetic field. We examine the effects of (i) the initial angle difference between the global magnetic field and the cloud rotation axis ($\theta_0$) and (ii) the ratio of the thermal to gravitational energy ($\alpha_0$). We study $16$ models in total and calculate the cloud evolution until $\sim \! 5000$ yr after protostar formation. Our simulation results indicate that an initial non-zero $\theta_0$ ($> 0$) promotes the disk formation but tends to suppress the outflow driving, for models that are moderately gravitationally unstable, $\alpha_0 \lesssim 1$. In these models, a large-sized rotationally-supported disk forms and a weak outflow appears, in contrast to a smaller disk and strong outflow in the aligned case ($\theta_0 = 0$). Furthermore, we find that when the initial cloud is highly unstable with small $\alpha_0$, the initial angle difference $\theta_0$ does not significantly affect the disk formation and outflow driving.

The formation of a close binary system is investigated using a three-dimensional resistive magnetohydrodynamics simulation. Starting from a prestellar cloud, the cloud evolution is calculated until about 400 yr after protostar formation. Fragmentation occurs in the gravitationally collapsing cloud and two fragments evolve into protostars. The protostars orbit each other and a protobinary system appears. A wide-angle low-velocity outflow emerges from the circumbinary streams that encloses two protostars, while each protostar episodically drives high-velocity jets. Thus, the two high-velocity jets are surrounded by the low-velocity circumbinary outflow. The speed of the jets exceeds $\gtrsim 100 km s^{-1}$. Although the jets have a collimated structure, they are swung back on the small scale and are tangled at the large scale due to the binary orbital motion. A circumstellar disk also appears around each protostar. In the early main accretion phase, the binary orbit is complicated, while the binary separation is within $<30$ au. For the first time, all the characteristics of protobinary systems recently observed with large telescopes are reproduced in a numerical simulation.

We have performed survey-type observations in 1 mm continuum and molecular lines toward dense cores (32 prestellar + 7 protostellar) with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in the Taurus molecular clouds using the ALMA-ACA (Atacama Compact Array) stand-alone mode with an angular resolution of 6.$''$5 ($\sim$900 au). The primary purpose of this study is to investigate the innermost part of dense cores toward understanding the initial condition of star formation. In the protostellar cores, contributions from protostellar disks dominate the observed continuum flux with a range of 35-90\% except for the very low-luminosity object. For the prestellar cores, we have successfully confirmed continuum emission from dense gas with a density of $\gtrsim$3 $\times$10$^5$ cm$^{-3}$ toward approximately one-third of the targets. Thanks to the lower spatial frequency coverage with the ACA-7 m array, the detection rate is significantly higher than that of the previous surveys, which have 0 or 1 continuum detected sources among large number of starless samples using the ALMA Main array. The statistical counting method tells us that the lifetime of the prestellar cores until protostar formation therein approaches the free-fall time as the density increases. Among the prestellar cores, at least two targets have possible internal substructures, which are detected in continuum emission with the size scale of $\sim$1000 au if we consider the molecular line (C$^{18}$O and N$_2$D$^{+}$) distributions. These results suggest that small scale fragmentation/coalescence processes within a $\lesssim$0.1 pc scale region occur to determine the final core mass leading to protostars before the dynamical collapse whose central density range is $\sim$(0.3-1) $\times$10$^6$ cm$^{-3}$.

The first hydrostatic core, the first quasi-hydrostatic object formed during the star formation process, is still the observational missing link between the prestellar and protostellar phases, mainly due to its short lifetime. Although we have not established a clear method to identify this rare object, recent theoretical studies predict that the first core has millimeter continuum emission and low-velocity outflow with a wide opening angle. An extensive continuum/outflow survey toward a large number of $"$starless$"$ cores in nearby star-forming regions works as a pathfinder. We observed 32 prestellar cores in Taurus with an average density of $\gtrsim$10$^5$ cm$^{-3}$ in 1.3 mm continuum and molecular lines using the ALMA-ACA stand-alone mode. Among the targets, MC35-mm centered at one of the densest $"$starless$"$ cores in Taurus has blueshifted/redshifted wings in the $^{12}$CO (2-1) line, indicating that there is deeply embedded object driving molecular outflow. The observed velocities and sizes of the possible outflow lobes are 2-4 km s$^{-1}$, and $\sim$2 $\times$10$^3$ au, respectively, and the dynamical time is calculated to be $\sim$10$^3$ yr. In addition to this, the core is one of the strongest N$_2$D$^{+}$ (3-2) emitter in our sample. All the observed signatures do not conflict with any of the theoretical predictions about the first hydrostatic core so far, and thus MC35-mm is unique as the only first-core candidate in the Taurus molecular cloud.

We observed polarization of the SiO rotational transitions from Orion Source I (SrcI) to probe the magnetic field in bipolar outflows from this high mass protostar. Both 43 GHz $J$=1-0 and 86 GHz $J$=2-1 lines were mapped with $\sim$20 AU resolution, using the Very Large Array (VLA) and Atacama Large Millimeter/Submillimeter Array (ALMA), respectively. The $^{28}$SiO transitions in the ground vibrational state are a mixture of thermal and maser emission. Comparison of the polarization position angles in the $J$=1-0 and $J$=2-1 transitions allows us to set an upper limit on possible Faraday rotation of $10^{4}$ radians m$^{-2}$, which would twist the $J$=2-1 position angles typically by less than 10 degrees. The smooth, systematic polarization structure in the outflow lobes suggests a well ordered magnetic field on scales of a few hundred AU. The uniformity of the polarization suggests a field strength of $\sim$30 milli-Gauss. It is strong enough to shape the bipolar outflow and possibly lead to sub-Keplerian rotation of gas at the base of the outflow. The strikingly high fractional linear polarizations of 80-90% in the $^{28}$SiO $v$=0 masers require anisotropic pumping. We measured circular polarizations of 60% toward the strongest maser feature in the $v$=0 $J$=1-0 peak. Anisotropic resonant scattering (ARS) is likely to be responsible for this circular polarization. We also present maps of the $^{29}$SiO $v$=0 $J$=2-1 maser and several other SiO transitions at higher vibrational levels and isotopologues.

We have observed the Class I protostar L1489 IRS with the Atacama Millimeter/submillimeter Array (ALMA) in Band 6. The C$^{18}$O $J=$2-1 line emission shows flattened and non-axisymmetric structures in the same direction as its velocity gradient due to rotation. We discovered that the C$^{18}$O emission shows dips at a radius of ~200-300 au while the 1.3 mm continuum emission extends smoothly up to r~400 au. At the radius of the C$^{18}$O dips, the rotational axis of the outer portion appears to be tilted by ~15 degrees from that of the inner component. Both the inner and outer components with respect to the C$^{18}$O dips exhibit the $r^{-0.5}$ Keplerian rotation profiles until r~600 au. These results not only indicate that a Keplerian disk extends up to ~600 au but also that the disk is warped. We constructed a three dimensional warped disk model rotating at the Keplerian velocity, and demonstrated that the warped disk model reproduces main observed features in the velocity channel maps and the PV diagrams. Such a warped disk system can form by mass accretion from a misaligned envelope. We also discuss a possible disk evolution scenario based on comparisons of disk radii and masses between Class I and Class II sources.

The accretion phase of star formation is investigated in magnetically-dominated clouds that have an initial subcritical mass-to-flux ratio. We employ nonideal magnetohydrodynamic simulations that include ambipolar diffusion and ohmic dissipation. During the early prestellar phase the mass-to-flux ratio rises toward the critical value for collapse, and during this time the angular momentum of the cloud core is reduced significantly by magnetic braking. Once a protostar is formed in the core, the accretion phase is characterized by the presence of a small amount of angular momentum but a large amount of magnetic flux in the near-protostellar environment. The low angular momentum leads to a very small (or even nonexistent) disk and weak outflow, while the large magnetic flux can lead to an interchange instability that rapidly removes flux from the central region. The effective magnetic braking in the early collapse phase can even lead to a counter-rotating disk and outflow, in which the rotation direction of the disk and outflow is opposite to that of the infalling envelope. The solutions with a counter-rotating disk, tiny disk, or nonexistent disk (direct collapse) are unique outcomes that are realized in collapse from magnetically-dominated clouds with an initial subcritical mass-to-flux ratio.

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.

A possible origin of kHz QPOs in low-mass X-ray binaries is proposed. Recent numerical MHD simulations of accretion disks with turbulent magnetic fields of MRI definitely show the presence of two-armed spiral structure in quasi-steady state of accretion disks. In such deformed disks, two-armed ($m=2$) c-mode ($n=1$) oscillations are excited by wave-wave resonant instability. Among these excited oscillations, the fundamental in the radial direction ($n_r=0$) will be the higher kHz QPO of a twin QPOs, and the first overtone ($n_r=1$) in the radial direction will be the lower kHz QPO of the twin. A possible cause of the twin high-frequency QPOs (HFQPOs) in BH X-ray binaries is also discussed.

Using resistive magnetohydrodynamics simulations, the propagation of protostellar jets, the formation of circumstellar discs and the configuration of magnetic fields are investigated from the prestellar cloud phase until $\sim$500 yr after protostar formation. As the initial state, we prepare magnetized rotating clouds, in which the rotation axis is misaligned with the global magnetic field by an angle $\theta_0$. We calculate the cloud evolution for nine models with different $\theta_0 (=$ 0, 5, 10, 30, 45, 60, 80, 85, 90$^\circ$). Our simulations show that there is no significant difference in the physical quantities of the protostellar jet, such as the mass and momentum, among the models except for the model with $\theta_0=90^\circ$. On the other hand, the directions of the jet, disc normal and magnetic field are never aligned with each other during the early phase of star formation except for the model with $\theta_0=0^\circ$. Even when the rotation axis of the prestellar cloud is slightly inclined to the global magnetic field, the directions of the jet, disc normal and local magnetic field differ considerably, and they randomly change over time. Our results indicate that it is very difficult to extract any information from the observations of the directions of the outflow, disc and magnetic field at the scale of $\lesssim 1000$ au. Thus, we cannot use such observations to derive any restrictions on the star formation process.

We have observed the submillimeter continuum condensations SMM2, SMM4, SMM9, and SMM11 in the star forming cluster Serpens Main using the Atacama Large Millimeter/submillimeter Array during Cycle 3 in the 1.3 mm continuum, 12CO J=2-1, SO J_N=6_5-5_4, and C18O J=2-1 lines at an angular resolution of ~0.55 (240 au). Sixteen sources have been detected in the 1.3 mm continuum, which can be classified into three groups. Group 1 consists of six sources showing extended continuum emission and bipolar/monopolar 12CO outflows. Although all the Group 1 members are classified as Class 0 protostars, our observations suggest evolutionary trends among them in terms of 12CO outflow dynamical time, SO emission distribution, C18O fractional abundance, and continuum morphology. Group 2 consists of four sources associated with a continuum filamentary structure and no 12CO outflows. Central densities estimated from the 1.3 mm continuum intensity suggest that they are prestellar sources in a marginally Jeans unstable state. Group 3 consists of six Spitzer sources showing point-like 1.3 mm continuum emission and clumpy 12CO outflows. These features of Group 3 suggest envelope dissipation, preventing disk growth from the present size, r <~ 60 au. The Group 3 members are protostars that may be precursors to the T Tauri stars associated with small disks at tens-au radii identified in recent surveys.

This article summarizes the work presented at the workshop 'The Power of Faraday Tomography: towards 3D mapping of cosmic magnetic fields', held in Miyazaki, Japan, in Spring 2018. We place the various oral and poster presentations given at the workshop in a broader perspective and present some highlight results from every presenter.

We report on the first bird's-eye view of the innermost accretion disk around the high-mass protostellar object G353.273+0.641, taken by Atacama Large Millimter/submillimeter Array long-baselines. The disk traced by dust continuum emission has a radius of 250 au, surrounded by the infalling rotating envelope traced by thermal CH$_3$OH lines. This disk radius is consistent with the centrifugal radius estimated from the specific angular momentum in the envelope. The lower-limit envelope mass is $\sim$5-7 M$_{\odot}$ and accretion rate onto the stellar surface is 3 $\times$ 10$^{-3}$ M$_{\odot}$ yr$^{-1}$ or higher. The expected stellar age is well younger than 10$^{4}$ yr, indicating that the host object is one of the youngest high-mass objects at present. The disk mass is 2-7 M$_{\odot}$, depending on the dust opacity index. The estimated Toomre's $Q$ parameter is typically 1-2 and can reach 0.4 at the minimum. These $Q$ values clearly satisfy the classical criteria for the gravitational instability, and are consistent with the recent numerical studies. Observed asymmetric and clumpy structures could trace a spiral arm and/or disk fragmentation. We found that 70$\%$ of the angular momentum in the accretion flow could be removed via the gravitational torque in the disk. Our study has indicated that the dynamical nature of a self-gravitating disk could dominate the early phase of high-mass star formation. This is remarkably consistent with the early evolutionary scenario of a low-mass protostar.

The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka Morita array) and the IRAM 30 m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H$_2$ volume density of $\sim$10$^6$ cm$^{-3}$, MC5-N, connected to a narrow (width $\sim$0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim$0.2-0.4 $M_{\odot}$, which is an order of magnitude smaller than typical low-mass prestellar cores. Taking into account a typical core to star formation efficiency for prestellar cores ($\sim$20%-40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N$_2$H$^+$ and N$_2$D$^+$ observations using the Nobeyama 45 m telescope have confirmed the high-deuterium fractionation ($\sim$30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.

The evolution of collapsing clouds embedded in different star-forming environments is investigated using three-dimensional non-ideal magnetohydrodynamics simulations considering different cloud metallicities ($Z/\thinspace Z_\odot$ = 0, 10$^{-5}$, 10$^{-4}$, 10$^{-3}$, 10$^{-2}$, 10$^{-1}$ and 1) and ionisation strengths ($C_\zeta$=0, 0.01, 1 and 10, where $C_\zeta$ is a coefficient controlling the ionisation intensity and $C_\zeta=1$ corresponds to the ionisation strength of nearby star-forming regions). With all combinations of these considered values of $Z/\thinspace Z_\odot$ and $C_\zeta$, 28 different star-forming environments are prepared and simulated. The cloud evolution in each environment is calculated until the central density reaches $n\approx10^{16}\,{\rm cm}^{-3}$ just before protostar formation, and the outflow driving conditions are derived. An outflow appears when the (first) adiabatic core forms in a magnetically active region where the magnetic field is well coupled with the neutral gas. In cases where outflows are driven, their momentum fluxes are always comparable to the observations of nearby star-forming regions. Thus, these outflows should control the mass growth of the protostars as in the local universe. Roughly, an outflow appears when $Z/\thinspace Z_\odot>10^{-4}$ and $C_\zeta \ge 0.01$. It is expected that the transition of the star formation mode from massive stars to normal solar-type stars occurs when the cloud metallicity is enhanced to the range of $Z/\thinspace Z_\odot\approx 10^{-4}$--$10^{-3}$, above which relatively low-mass stars would preferentially appear as a result of strong mass ejection.

We present the observational evidence of a pseudo-disc around the proto-brown dwarf Mayrit 1701117, the driving source of the large-scale HH~1165 jet. Our analysis is based on ALMA $^{12}$CO (2-1) line and 1.37 mm continuum observations at an angular resolution of $\sim$0.4$^{\prime\prime}$. The pseudo-disc is a bright feature in the CO position-velocity diagram (PVD), elongated in a direction perpendicular to the jet axis, with a total (gas+dust) mass of $\sim$0.02 M$_{\odot}$, size of 165-192 AU, and a velocity spread of $\pm$2 km s$^{-1}$. The large velocity gradient is a combination of infalling and rotational motions, indicating a contribution from a pseudo-disc and an unresolved inner Keplerian disc. There is weak emission detected in the H$_{2}$CO (3-2) and N$_{2}$D$^{+}$ (3-2) lines. H$_{2}$CO emission likely probes the inner Keplerian disc where CO is expected to be frozen, while N$_{2}$D$^{+}$ possibly originates from an enhanced clump at the outer edge of the pseudo-disc. We have considered various models (core collapse, disc fragmentation, circum-binary disc) that can fit both the observed CO spectrum and the position-velocity offsets. The observed morphology, velocity structure, and the physical dimensions of the pseudo-disc are consistent with the predictions from the core collapse simulations for brown dwarf formation. From the best model fit, we can constrain the age of the proto-brown dwarf system to be $\sim$30,000-40,000 yr. A comparison of the H$_{2}$ column density derived from the CO line and 1.37 mm continuum emission indicates that only about 2% of the CO is depleted from the gas phase.

Starting from a prestellar core with a size of $1.2\times10^4$ AU, we calculate the evolution of a gravitationally collapsing core until $\sim2000$ yr after protostar formation using a three-dimensional resistive magnetohydrodynamic simulation, in which the protostar is resolved with a spatial resolution of $5.6\times10^{-3}$ AU. Following protostar formation, a rotationally supported disk is formed. Although the disk size is as small as $\sim2-4$ AU, it remains present until the end of the simulation. Since the magnetic field dissipates and the angular momentum is then not effectively transferred by magnetic effects, the disk surface density gradually increases and spiral arms develop due to gravitational instability. The disk angular momentum is then transferred mainly by gravitational torques, which induce an episodic mass accretion onto the central protostar. The episodic accretion causes a highly time-variable mass ejection (the high-velocity jet) near the disk inner edge, where the magnetic field is well coupled with the neutral gas. As the mass of the central protostar increases, the jet velocity gradually increases and exceeds $\sim100$ km s$^{-1}$. The jet opening angle widens with time at its base, while the jet keeps a very good collimation on the large scale. In addition, a low-velocity outflow is driven from the disk outer edge. A cavity-like structure, a bow shock and several knots, all of which are usually observed in star-forming regions, are produced in the outflowing region.

Recent observations uncover various phenomena around the protostar such as misalignment between the outflow and magnetic field, precession of the jet, and time variability of the ejected clumps, whose origins are under debate. We perform a three-dimensional resistive magnetohydrodynamics simulation of the protostar formation in a star-forming core whose rotation axis is tilted at an angle $45^\circ$ with respect to the initial magnetic field, in which the protostar is resolved with a spatial resolution of $0.01\,$au. In low-dense outer region, the prestellar core contracts along the magnetic field lines due to the flux freezing. In high-dense inner region, on the other hand, the magnetic dissipation becomes efficient and weakens the magnetic effects when the gas number density exceeds about $10^{11}\,{\rm cm^{-3}}$. Then, the normal direction of the flattened disc is aligned with the angular momentum vector. The outflow, jet, and protostellar ejection are driven from different scales of the circumstellar disc and spout in different directions normal to the warped disc. These axes do not coincide with the global magnetic field direction and vary with time. This study demonstrates that a couple of misalignment natures reported by observations can be simultaneously reproduced only by assuming the star-forming core rotating around a different direction from the magnetic field.

We present high-resolution images of the submillimeter SiO line emissions of a massive young stellar object Orion Source I using the Atacama Large Millimeter/ Submillimeter Array (ALMA) at band 8. We detected the 464 GHz SiO v=4 J=11-10 line in Source I, which is the first detection of the SiO v=4 line in star-forming regions, together with the 465 GHz 29SiO v=2 J=11-10 and the 428 GHz SiO v=2 J=10-9 lines with a resolution of 50 AU. The 29SiO v=2 J=11-10 and SiO v=4 J=11-10 lines have compact structures with the diameter of <80 AU. The spatial and velocity distribution suggest that the line emissions are associated with the base of the outflow and the surface of the edge-on disk. In contrast, SiO v=2 J=10-9 emission shows a bipolar structure in the direction of northeast-southwest low-velocity outflow with ~200 AU scale. The emission line exhibits a velocity gradient along the direction of the disk elongation. With the assumption of the ring structure with Keplerian rotation, we estimated the lower limit of the central mass to be 7 solar mass and the radius of 12 AU< r <26 AU.

The Hall effect plays a significant role in star formation because it induces rotation in the infalling envelope, which in turn affects the formation and evolution of the circumstellar disk. The importance of the Hall effect varies with the Hall coefficient, and this coefficient is determined by the fractional abundances of charged species. These abundance values are primarily based on the size and quantity of dust grains as well as the cosmic ray intensity, which respectively absorb and create charged species. Thus, the Hall coefficient varies with both the properties of dust grains and the cosmic ray rate (or ionization source). In this study, we explore the dependence of the Hall coefficient on the grain size and cosmic ray ionization rate using a simplified chemical network model. Following this, using an analytic model, we estimate the typical size of a circumstellar disk induced solely by the Hall effect. The results show that the disk grows during the main accretion phase to a size of ${\sim}$ 3 - 100 au, with the actual size depending on the parameters. These findings suggest that the Hall effect greatly affects circumstellar disk formation, especially in the case that the dust grains have a typical size of ${\sim}$ 0.025 ${\mu}$m - 0.075 ${\mu}$m.

Using the $\approx$15km ALMA long baselines, we imaged the Stokes $I$ emission and linearly polarized intensity ($PI$) in the 1.1-mm continuum band of a very young intermediate-mass protostellar source, MMS 6, in the Orion Molecular Cloud-3. The achieved angular resolution, $0''.02{\times}0''.03$ ($\approx$10 AU), shows for the first time a wealth of data on the dust emission polarization in the central 200 AU of a protostar. The $PI$ peak is offset to the south-west (SW) by $\approx$20 AU with respect to the Stokes $I$ peak. Its polarization degree is 11 % with its $E$-vector orientation of P.A.${\approx}135^{\circ}$. A partial ring-like structure with a radius of $\approx$80 AU is detected in $PI$ but not in the Stokes $I$. NW (north-west) and SE (south-east) parts of the ring are bright with a high polarization degree of $\gtrsim$10 %, and their $E$-vector orientations are roughly orthogonal to those observed near the center. We also detected arm-like polarized structures, extending to 1000 AU scale to the north, with the $E$-vectors aligned along the minor axis of the structures. We explored possible origins of the polarized emission comparing with magnetohydrodynamical (MHD) simulations of the toroidal wrapping of the magnetic field. The simulations are consistent with the $PI$ emission in the ring-like and the extended arm-like structures observed with ALMA. However, the current simulations do not completely reproduce observed polarization characteristics in the central 50 AU. Although the self-scattering model can explain the polarization pattern and positional offset between the Stokes $I$ and $PI$, this model is not able to reproduce the observed high degree of polarization.

Both high- and low-velocity outflows are occasionally observed around a protostar by molecular line emission. The high-velocity component is called `Extremely High-Velocity (EHV) flow,' while the low-velocity component is simply referred as `(molecular) outflow.' This study reports a newly found EHV flow and outflow around MMS $5$ in the Orion Molecular Cloud 3 observed with ALMA. In the observation, CO $J$=2--1 emission traces both the EHV flow ($|v_{\rm{LSR}} - v_{\rm{sys}}|$ $\simeq$ 50--100 $\rm{km\ s^{-1}}$) and outflow ($|v_{\rm{LSR}} - v_{\rm{sys}}|$ $\simeq$ 10--50 $\rm{km\ s^{-1}}$). On the other hand, SiO $J$=5--4 emission only traces the EHV flow. The EHV flow is collimated and located at the root of the V-shaped outflow. The CO outflow extends up to $\sim$ 14,000\,AU with a position angle (P.A.) of $\sim79^\circ$ and the CO redshifted EHV flow extends to $\sim$11,000 AU with P.A. $\sim96^\circ$. The EHV flow is smaller than the outflow, and the dynamical timescale of the EHV flow is shorter than that of the outflow by a factor of $\sim 3$. The flow driving mechanism is discussed based on the size, time scale, axis difference between the EHV flow and outflow, and the periodicity of the knots. Our results are consistent with the nested wind scenario, although the jet entrainment scenario could not completely be ruled out.

The effect of misalignment between the magnetic field $\magB$ and the angular momentum $\Jang$ of molecular cloud cores on the angular momentum evolution during the gravitational collapse is investigated by ideal and non-ideal MHD simulations. For the non-ideal effect, we consider the ohmic and ambipolar diffusion. Previous studies that considered the misalignment reported qualitatively contradicting results. Magnetic braking was reported as being either strengthened or weakened by misalignment in different studies. We conducted simulations of cloud-core collapse by varying the stability parameter $\alpha$ (the ratio of the thermal to gravitational energy of the core) with and without including magnetic diffusion The non-ideal MHD simulations show the central angular momentum of the core with $\theta=0^\circ$ ($\Jang \parallel \magB$) being always greater than that with $\theta=90^\circ$ ($\Jang \perp \magB$), independently of $\alpha$, meaning that circumstellar disks form more easily form in a core with $\theta=0^\circ$. The ideal MHD simulations, in contrast, show the the central angular momentum of the core with $\theta=90^\circ$ being greater than with $\theta=0^\circ$ for small $\alpha$, and is smaller for large $\alpha$. Inspection of the angular momentum evolution of the fluid elements reveals three mechanisms contributing to the evolution of the angular momentum: (i) magnetic braking in the isothermal collapse phase, (ii) selective accretion of the rapidly (for $\theta=90^\circ$ ) or slowly (for $\theta=0^\circ$) rotating fluid elements to the central region, and (iii) magnetic braking in the first-core and the disk. The difference between the ideal and non-ideal simulations arises from the different efficiencies of (iii).

Observational study of galactic magnetic fields is limited by projected observables. Comparison with numerical simulations is helpful to understand the real structures, and observational visualization of numerical data is an important task. Machida et al. (2018) have reported Faraday depth maps obtained from numerical simulations. They showed that the relation between azimuthal angle and Faraday depth depends on the inclination angle. In this paper, we investigate 100MHz to 10GHz radio synchrotron emission from spiral galaxies, using the data of global three-dimensional magneto-hydrodynamic simulations. We model internal and external Faraday depolarization at small scales and assume a frequency independent depolarization. It is found that the internal and external Faraday depolarization becomes comparable inside the disk and the dispersion of Faraday depth becomes about 4rad/m^{2} for face-on view and 40rad/m2 for edge-on view, respectively. The internal depolarization becomes ineffective in the halo. Because of the magnetic turbulence inside the disk, frequency independent depolarization works well and the polarization degree becomes 0.3 at high frequency. When the observed frequency is in the 100 MHz band, polarized intensity vanishes in the disk, while that from the halo can be observed. Because the remaining component of polarized intensity is weak in the halo and the polarization degree is about a few %, it may be difficult to observe that component. These results indicate that the structures of global magnetic fields in spiral galaxies could be elucidated, if broadband polarimetry such as that with the Square Kilometre Array is achieved.

Observational study of galactic magnetic fields is hampered by the fact that the observables only probe various projections of the magnetic fields. Comparison with numerical simulations is helpful to understand the real structures, and observational visualization of numerical data is an important task. In this paper, we investigate 8~GHz radio synchrotron emission from spiral galaxies, using the data of global three-dimensional magneto-hydrodynamic simulations. We assume a frequency independent depolarization in our observational visualization. We find that the appearance of the global magnetic field depends on the viewing angle: a face-on view seemingly has hybrid magnetic field types combining axisymmetric modes with higer order modes; at a viewing angle of $\sim 70\degr$, the galaxy seems to contain a ring-like magnetic field structure; while in edge-on view, only field structure parallel to the disk can be seen. The magnetic vector seen at 8~GHz traces the global magnetic field inside the disk. These results indicate that the topology of global magnetic field obtained from the relation between azimuthal angle and Faraday depth strongly depends on the viewing anglue of the galaxy. As one of the examples, we compare our results at a viewing angle of $25\degr$ with the results of IC342. The relation between azimuthal angle and Faraday depth of the numerical result shows a tendency similar to IC342, such as the peak numbers of the Faraday depth.

On the basis of various data such as ALMA, JVLA, Chandra, {\it Herschel}, and {\it Spitzer}, we confirmed that two protostellar candidates in Oph-A are bona fide protostars or proto-brown dwarfs (proto-BDs) in extremely early evolutionary stages. Both objects are barely visible across infrared (IR, i.e., near-IR to far-IR) bands. The physical nature of the cores is very similar to that expected in first hydrostatic cores (FHSCs), objects theoretically predicted in the evolutionary phase prior to stellar core formation with gas densities of $\sim$ 10$^{11-12}$ cm$^{-3}$. This suggests that the evolutionary stage is close to the FHSC formation phase. The two objects are associated with faint X-ray sources, suggesting that they are in very early phase of stellar core formation with magnetic activity. In addition, we found the CO outflow components around both sources which may originate from the young outflows driven by these sources. The masses of these objects are calculated to be $\sim 0.01-0.03$ $M_\odot$ from the dust continuum emission. Their physical properties are consistent with that expected from the numerical model of forming brown dwarfs. These facts (the X-ray detection, CO outflow association, and FHSC-like spectral energy distributions) strongly indicate that the two objects are proto-BDs or will be in the very early phase of protostars which will evolve more massive protostars if they gain enough mass from the surroundings. The ages of these two objects are likely to be within $\sim 10^3$ years after the protostellar core (or second core) formation, taking into account the outflow dynamical times ($\lesssim$ 500 yrs).

We report ALMA Cycle 3 observations in CO isotopes toward a dense core, MC27/L1521F in Taurus, which is considered to be at an early stage of multiple star formation in a turbulent environment. Although most of the high-density parts of this core are considered to be as cold as $\sim$10 K, high-angular resolution ($\sim$20 au) observations in $^{12}$CO ($J$ = 3--2) revealed complex warm ($>$15--60 K) filamentary/clumpy structures with the sizes from a few tens of au to $\sim$1,000 au. The interferometric observations of $^{13}$CO and C$^{18}$O show that the densest part with arc-like morphologies associated with the previously identified protostar and condensations are slightly redshifted from the systemic velocity of the core. We suggest that the warm CO clouds may be consequences of shock heating induced by interactions among the different density/velocity components that originated from the turbulent motions in the core. However, such a small-scale and fast turbulent motion does not correspond to a simple extension of the line-width-size relation (i.e., Larson'{}s law), and thus the actual origin remains to be studied. The high-angular resolution CO observations are expected to be essential in detecting small-scale turbulent motions in dense cores and to investigate protostar formation therein.

We have observed the submillimeter continuum condensation SMM4 in Serpens Main using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 3 in 1.3 mm continuum, 12CO J=2-1, SO J_N=6_5-5_4, and C18O J=2-1 lines at angular resolutions of ~0.55" (240 au). The 1.3 mm continuum emission shows that SMM4 is spatially resolved into two protostars embedded in the same core: SMM4A showing a high brightness temperature, 18 K, with little extended structure and SMM4B showing a low brightness temperature, 2 K, with compact and extended structures. Their separation is ~2100 au. Analysis of the continuum visibilities reveals a disk-like structure with a sharp edge at r~240 au in SMM4A, and a compact component with a radius of 56 au in SMM4B. The 12CO emission traces fan-shaped and collimated outflows associated with SMM4A and SMM4B, respectively. The blue and red lobes of the SMM4B outflow have different position angles by ~30 deg. Their inclination and bending angles in the 3D space are estimated at i_b~ 36 deg, i_r~70 deg, and alpha~40 deg, respectively. The SO emission traces shocked regions, such as cavity walls of outflows and the vicinity of SMM4B. The C18O emission mainly traces an infalling and rotating envelope around SMM4B. The C18O fractional abundance in SMM4B is ~50 times smaller than that of the interstellar medium. These results suggest that SMM4A is more evolved than SMM4B. Our studies in Serpens Main demonstrate that continuum and line observations at millimeter wavelengths allow us to differentiate evolutionary phases of protostars within the Class 0 phase.

The Telescope Array observatory utilizes fluorescence detectors and surface detectors to observe air showers produced by ultra high energy cosmic rays in the Earth's atmosphere. Cosmic ray events observed in this way are termed hybrid data. The depth of air shower maximum is related to the mass of the primary particle that generates the shower. This paper reports on shower maxima data collected over 8.5 years using the Black Rock Mesa and Long Ridge fluorescence detectors in conjunction with the array of surface detectors. We compare the means and standard deviations of the observed $X_{\mathrm{max}}$ distributions with Monte Carlo $X_{\mathrm{max}}$ distributions of unmixed protons, helium, nitrogen, and iron, all generated using the QGSJet~II-04 hadronic model. We also perform an unbinned maximum likelihood test of the observed data, which is subjected to variable systematic shifting of the data $X_{\mathrm{max}}$ distributions to allow us to test the full distributions, and compare them to the Monte Carlo to see which elements are not compatible with the observed data. For all energy bins, QGSJet~II-04 protons are found to be compatible with Telescope Array hybrid data at the 95% confidence level after some systematic $X_{\mathrm{max}}$ shifting of the data. Three other QGSJet~II-04 elements are found to be compatible using the same test procedure in an energy range limited to the highest energies where data statistics are sparse.

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.

We report a polarization analysis of the eastern region of W50, observed with the Australia Telescope Compact Array (ATCA) at 1.4 - 3.0 GHz. In order to study the physical structures in the region where the SS433 jet and W50 interact, we obtain an intrinsic magnetic field vector map of that region. We find that the orientation of the intrinsic magnetic field vectors are aligned along the total intensity structures, and that there are characteristic, separate structures related to the jet, the bow shock, and the terminal shock. The Faraday rotation measures (RMs), and the results of Faraday Tomography suggest that a high intensity, filamentary structure in the north-south direction of the eastern-edge region can be separated into at least two parts to the north and south. The results of Faraday Tomography also show that there are multiple components along the line of sight and/or within the beam area. In addition, we also analyze the X-ray ring-like structure observed with XMM-Newton. While the possibility still remains that this X-ray ring is real, it seems that the structure is not ring-like at radio wavelengths. Finally, we suggest that the structure is a part of the helical structure that coils the eastern ear of W50.

In nearby star-forming clouds, amplification and dissipation of the magnetic field are known to play crucial roles in the star-formation process. The star-forming environment varies from place to place and era to era in galaxies. In the present study, amplification and dissipation of magnetic fields in star-forming clouds are investigated under different environments using magnetohydrodynamics (MHD) simulations. We consider various star-forming environments in combination with the metallicity and the ionization strength, and prepare prestellar clouds having two different mass-to-flux ratios. We calculate the cloud collapse until protostar formation using ideal and non-ideal (inclusion and exclusion of Ohmic dissipation and ambipolar diffusion) MHD calculations to investigate the evolution of the magnetic field. We perform 288 runs in total and show the diversity of the density range within which the magnetic field effectively dissipates, depending on the environment. In addition, the dominant dissipation process (Ohmic dissipation or ambipolar diffusion) is shown to strongly depend on the star-forming environment. Especially, for the primordial case, magnetic field rarely dissipates without ionization source, while it efficiently dissipates when very weak ionization sources exist in the surrounding environment. The results of the present study help to clarify star formation in various environments.

One of the outstanding problems in star-formation theory concerns the transfer of angular momentum such that mass can accrete onto a newly born young stellar object (YSO). From a theoretical standpoint, outflows and jets are predicted to play an essential role in angular momentum transfer and their rotation motions have been reported for both low- and high-mass YSOs. However, little quantitative discussion on outflow launching mechanisms have been presented for high-mass YSOs due to a lack of observational data. Here we present a clear signature of rotation in the bipolar outflow driven by Orion Source I, a high-mass YSO candidate, using the Atacama Large Millimeter/Submillimeter Array (ALMA). A rotational transition of silicon monoxide (Si18O) reveals a velocity gradient perpendicular to the outflow axis which is consistent with that of the circumstellar disk traced by a high-excitation water (H2O) line. The launching radii and outward velocity of the outflow are estimated to be >10 au and 10 km s-1, respectively. These parameters rule out a possibility that the observed outflow is produced by entrainment of a high-velocity jet, and that contribution from stellar-wind or X-wind which have smaller launching radii are significant in the case of Source I. Thus, present results provide a convincing evidence of a rotating outflow directly driven by the magneto-centrifugal disk wind launched by a high-mass YSO candidate.

We present observational results of the submillimeter H2O and SiO lines toward a candidate high-mass young stellar object Orion Source I using ALMA. The spatial structures of the high excitation lines at lower-state energies of >2500 K show compact structures consistent with the circumstellar disk and/or base of the northeast-southwest bipolar outflow with a 100 au scale. The highest excitation transition, the SiO (v=4) line at band 8, has the most compact structure. In contrast, lower-excitation transitions are more extended than 200 au tracing the outflow. Almost all the line show velocity gradients perpendicular to the outflow axis suggesting rotation motions of the circumstellar disk and outflow. While some of the detected lines show broad line profiles and spatially extended emission components indicative of thermal excitation, the strong H2O lines at 321 GHz, 474 GHz, and 658 GHz with brightness temperatures of >1000 K show clear signatures of maser action.

The driving mechanism of massive outflows observed in high-mass star-forming regions is investigated using three-dimensional magnetohydrodynamics (MHD) and protostellar evolution calculations. In our previous paper, we showed that the mass outflow rate depends strongly on the mass accretion rate onto the circumstellar disk around a high-mass protostar, and massive outflows may be driven by the magnetic effect in high-mass star-forming cores. In the present study, in order to verify that the MHD disk wind is the primary driving mechanism of massive outflows, we quantitatively compare outflow properties obtained through simulations and observations. Since the outflows obtained through simulations are slightly younger than those obtained through observations, the time-integrated quantities of outflow mass, momentum, and kinetic energy are slightly smaller than those obtained through observations. On the other hand, time-derivative quantities of mass ejection rate, outflow momentum flux, and kinetic luminosity obtained through simulations are in very good agreement with those obtained through observations. This indicates that the MHD disk wind greatly contributes to the massive outflow driving from high-mass protostars, and the magnetic field might significantly control the high-mass star formation process.

We report ALMA observations in 0.87 mm continuum and $^{12}$CO ($J$ = 3--2) toward a very low-luminosity ($<$0.1 $L_{\odot}$) protostar, which is deeply embedded in one of the densest core MC27/L1521F, in Taurus with an indication of multiple star formation in a highly dynamical environment. The beam size corresponds to $\sim$20 AU, and we have clearly detected blueshifted/redshifted gas in $^{12}$CO associated with the protostar. The spatial/velocity distributions of the gas show there is a rotating disk with a size scale of $\sim$10 AU, a disk mass of $\sim$10$^{-4}$ $M_{\odot}$ and a central stellar mass of $\sim$0.2 $M_{\odot}$. The observed disk seems to be detached from the surrounding dense gas, although it is still embedded at the center of the core whose density is $\sim$10$^{6}$ cm$^{-3}$. The current low-outflow activity and the very low luminosity indicate that the mass accretion rate onto the protostar is extremely low in spite of a very early stage of star formation. We may be witnessing the final stage of the formation of $\sim$0.2 $M_{\odot}$ protostar. However, we cannot explain the observed low luminosity with the standard pre-main-sequence evolutionary track unless we assume cold accretion with an extremely small initial radius of the protostar ($\sim$0.65 $R_\odot$). These facts may challenge our current understanding of the low mass star formation, in particular the mass accretion process onto the protostar and the circumstellar disk.

We report our current SMA and ALMA studies of disk and planet formation around protostars. We have revealed that $r \gtrsim$100 AU scale disks in Keplerian rotation are ubiquitous around Class I sources. These Class I Keplerian disks are often embedded in rotating and infalling protostellar envelopes. The infalling speeds of the protostellar envelopes are typically $\sim$ 3-times smaller than the free-fall velocities, and the rotational profiles follow the $r^{-1}$ profile, that is, rotation with the conserved specific angular momentum. Our latest high-resolution ($\sim$0$\farcs$5) ALMA studies, as well as the other studies in the literature, have unveiled that $r \sim$100-AU scale Keplerian disks are also present in several Class 0 protostars, while in the other Class 0 sources the inferred upper limits of the Keplerian disks are very small ($r \lessim$20 AU). Our recent data analyses of the ALMA long baseline data of the Class I-II source HL Tau have revealed gaps in molecular gas as well as in dust in the surrounding disk, suggesting the presence of sub-Jovian planets in the disk. These results imply that disk and planet formation should be completed in the protostellar stage.

The driving mechanism of protostellar outflows and jets and their effects on the star formation process obtained from recent theoretical and numerical studies are described. Low-velocity outflows are driven by an outer region of the circumstellar disk, while high-velocity jets are driven near an inner edge of the disk. The disk angular momentum is effectively transferred by magnetic effects in the outflow and jet driving regions where the magnetic field is well coupled with neutral gas. On the other hand, in a high density gas region of the disk (or intermediate region), the magnetic field dissipates and is decoupled from neutral gas. Thus, in such a magnetically inactive region, no outward flow appears and the disk angular momentum is not effectively transferred by magnetic effects. Therefore, in the disk intermediate region, the disk surface density continues to increase and gravitational instability occurs and produce a non-axisymmetric (or spiral) structure. After spiral arms sufficiently develop, the disk angular momentum is transferred by gravitational torque and a large amount of the disk mass accretes onto the central protostar from the circumstellar disk. The episodic accretion induces time-variable high-velocity jets. The jets do not significantly contribute to a dynamical evolution of the protostar and circumstellar disk, while the low-velocity outflows can eject a large fraction of the infalling gas and determine the final stellar mass.

We report the discovery of an extremely young protostar, SMM11, located in the associated submillimeter condensation SMM11 in the Serpens Main cluster using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 3 at 1.3 mm and an angular resolution of ~0.5"~210 AU. SMM11 is a Class 0 protostar without any counterpart at 70 um or shorter wavelengths. The ALMA observations show 1.3 mm continuum emission associated with a collimated 12CO bipolar outflow. Spitzer and Herschel data show that SMM11 is extremely cold (T_bol=26 K) and faint (L_bol<~0.9 Lsun). We estimate the inclination angle of the outflow to be ~80 deg, almost parallel to the plane of the sky, from simple fitting using wind-driven-shell model. The continuum visibilities consist of Gaussian and power-law components, suggesting a spherical envelope with a radius of ~600 AU around the protostar. The estimated low C18O abundance, X(C18O)=1.5-3x10^-10, is also consistent with its youth. The high outflow velocity, a few 10 kms^-1 at a few 1000 AU, is much higher than theoretical simulations of first hydrostatic cores and we suggest that SMM11 is a transitional object right after the second collapse of the first core.

Star formation in magnetically subcritical clouds is investigated using a three-dimensional non-ideal magneto-hydrodynamics simulation. Since rapid cloud collapse is suppressed until the magnetic flux is sufficiently removed from the initially magnetically subcritical cloud by ambipolar diffusion, it takes > 5-10t_ff to form a protostar, where t_ff is the freefall timescale of the initial cloud. The angular momentum of the star forming cloud is efficiently transferred to the interstellar medium before the rapid collapse begins, and the collapsing cloud has a very low angular momentum. Unlike the magnetically supercritical case, no large-scale low-velocity outflow appears in such a collapsing cloud due to the short lifetime of the first core. Following protostar formation, a very weak high-velocity jet, which has a small momentum and might disappear at a later time, is driven near the protostar, while the circumstellar disc does not grow during the early mass accretion phase. The results show that the star formation process in magnetically subcritical clouds is qualitatively different from that in magnetically supercritical clouds.

We perform three-dimensional radiation non-ideal magnetohydrodynamics simulations and investigate the impact of the Hall effect on the angular momentum evolution in the collapsing cloud cores in which the magnetic field $\magB$ and angular momentum $\angJ$ are misaligned with each other. We find that the Hall effect notably changes the magnetic torques in the pseudo-disk, and strengthens and weakens the magnetic braking in cores with an acute and obtuse relative angles between $\magB$ and $\angJ$, respectively. This suggests that the bimodal evolution of the disk size may occur in early disk evolutionary phase even if $\magB$ and $\angJ$ are randomly distributed. We show that a counter-rotating envelope form in the upper envelope of the pseudo-disk in cloud cores with obtuse relative angles. We also find that a counter-rotating region forms at the midplane of the pseudo-disk in cloud cores with acute relative angles. The former and latter types of counter-rotating envelopes may be associated with the YSOs with a large ($r\sim100$ AU) and small ($r\lesssim10$ AU) disks, respectively.

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 have newly observed the Class 0/I protostar L1527 IRS using the Atacama Large Millimeter/submillimeter Array (ALMA) during its Cycle 1 in 220 GHz dust continuum and C18O (J=2-1) line emissions with a ~2 times higher angular resolution (~0.5") and ~4 times better sensitivity than our ALMA Cycle 0 observations. Continuum emission shows elongation perpendicular to the associated outflow, with a deconvolved size of 0.53"x0.15". C18O emission shows similar elongation, indicating that both emissions trace the disk and the flattened envelope surrounding the protostar. The velocity gradient of the C18O emission along the elongation due to rotation of the disk/envelope system is re-analyzed, identifying Keplerian rotation proportional to r^-0.5 more clearly than the Cycle 0 observations. The Keplerian-disk radius and the dynamical stellar mass are kinematically estimated to be ~74 AU and ~0.45 Mo, respectively. The continuum visibility is fitted by models without any annulus averaging, revealing that the disk is in hydrostatic equilibrium. The best-fit model also suggests a density jump by a factor of ~5 between the disk and the envelope, suggesting that disks around protostars can be geometrically distinguishable from the envelope from a viewpoint of density contrast. Importantly, the disk radius geometrically identified with the density jump is consistent with the radius kinematically estimated. Possible origin of the density jump due to the mass accretion from the envelope to the disk is discussed. C18O observations can be reproduced by the same geometrical structures derived from the dust observations, with possible C18O freeze-out and localized C18O desorption.

The relation between the mass accretion rate onto the circumstellar disc and the rate of mass ejection by magnetically driven winds is investigated using three-dimensional magnetohydrodynamics simulations. Using a spherical cloud core with a varying ratio of thermal to gravitational energy, which determines the mass accretion rate onto the disc, to define the initial conditions, the outflow propagation for approximately 10^4 yr after protostar formation is then calculated for several cloud cores. The mass ejection rate and accretion rate are comparable only when the magnetic energy of the initial cloud core is comparable to the gravitational energy. Consequently, in strongly magnetised clouds a higher mass accretion rate naturally produces both massive protostars and massive outflows. The simulated outflow mass, momentum, kinetic energy and momentum flux agree well with observations, indicating that massive stars form through the same mechanism as low-mass stars but require a significantly strong magnetic field to launch massive outflows.

We investigate the formation of circumstellar disks and outflows subsequent to the collapse of molecular cloud cores with the magnetic field and turbulence. Numerical simulations are performed by using an adaptive mesh refinement to follow the evolution up to $\sim 1000$~yr after the formation of a protostar. In the simulations, circumstellar disks are formed around the protostars; those in magnetized models are considerably smaller than those in nonmagnetized models, but their size increases with time. The models with stronger magnetic field tends to produce smaller disks. During evolution in the magnetized models, the mass ratios of a disk to a protostar is approximately constant at $\sim 1-10$\%. The circumstellar disks are aligned according to their angular momentum, and the outflows accelerate along the magnetic field on the $10-100$~au scale; this produces a disk that is misaligned with the outflow. The outflows are classified into two types: a magneto-centrifugal wind and a spiral flow. In the latter, because of the geometry, the axis of rotation is misaligned with the magnetic field. The magnetic field has an internal structure in the cloud cores, which also causes misalignment between the outflows and the magnetic field on the scale of the cloud core. The distribution of the angular momentum vectors in a core also has a non-monotonic internal structure. This should create a time-dependent accretion of angular momenta onto the circumstellar disk. Therefore, the circumstellar disks are expected to change their orientation as well as their sizes in the long-term evolutions.

We study formation and long-term evolution of a circumstellar disk in a collapsing molecular cloud core using a resistive magnetohydrodynamic simulation. While the formed circumstellar disk is initially small, it grows as accretion continues and its radius becomes as large as 200 AUs toward the end of the Class-I phase. A pair of grand-design spiral arms form due to gravitational instability in the disk, and they transfer angular momentum in the highly resistive disk. Although the spiral arms disappear in a few rotations as expected in a classical theory, new spiral arms form recurrently as the disk soon becomes unstable again by gas accretion. Such recurrent spiral arms persist throughout the Class-0 and I phase. We then perform synthetic observations and compare our model with a recent high-resolution observation of a young stellar object Elias 2-27, whose circumstellar disk has grand design spiral arms. We find good agreement between our theoretical model and the observation. Our model suggests that the grand design spiral arms around Elias 2-27 are consistent with material arms formed by gravitational instability. If such spiral arms commonly exist in young circumstellar disks, it implies that young circumstellar disks are considerably massive and gravitational instability is the key process of angular momentum transport.

We have measured continuum flux densities of a high-mass protostar candidate, a radio source I in the Orion KL region (Orion Source I) using the Atacama Large Millimeter/Submillimeter Array (ALMA) at band 8 with an angular resolution of 0.1". The continuum emission at 430, 460, and 490 GHz associated with Source I shows an elongated structure along the northwest-southeast direction perpendicular to the so-called low-velocity bipolar outflow. The deconvolved size of the continuum source, 90 au times 20 au, is consistent with those reported previously at other millimeter/submillimeter wavelength. The flux density can be well fitted to the optically thick black-body spectral energy distribution (SED), and the brightness temperature is evaluated to be 700-800 K. It is much lower than that in the case of proton-electron or H- free-free radiations. Our data are consistent with the latest ALMA results by Plambeck & Wright (2016), in which the continuum emission have been proposed to arise from the edge-on circumstellar disk via thermal dust emission, unless the continuum source consists of an unresolved structure with the smaller beam filling factor.

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.

Angular momentum distribution and its redistribution are of crucial importance in formation and evolution of circumstellar disks. Many molecular line observations toward young stellar objects indicate that radial distributions of the specific angular momentum $j$ are more or less universal. In small scales, typically $R\lesssim$ 100 AU, the specific angular momenta distribute like $j\propto r^{1/2}$, indicating existence of rotationally supported disk. In outer regions, $R\gtrsim$ 5000 AU, $j$ increases as the radius increases and the slope is steeper than unity, which is supposed to reflect the original angular momentum distributions in the maternal molecular clouds. And lastly there is a connecting region, 100 AU $\lesssim R \lesssim$ 5000 AU, in which $j$-distribution looks almost flat. While this is often interpreted as a consequence of conservation of the specific angular momentum, it actually is insufficient and requires a stronger condition that the initial distribution of $j$ must be spatially uniform. However, this requirement is unrealistic and inconsistent with observations. In this work, we propose a simple alternative explanation; the flat $j$ profile is produced by strong prolongation due to the large velocity gradient in the accreting flow no matter what the initial $j$ distribution is. We provide a simple analytic model for gravitational collapse of molecular clouds. This model can be used to estimate ages of protostars based solely on the observed rotation profile. We demonstrate its validity in comparison with hydrodynamic simulations, and apply the model to young stellar objects such as L1527 IRS, TMC-1A and B335.

We present the results of ALMA observations of dust continuum emission and molecular rotational lines toward a dense core MC27 (aka L1521F) in Taurus, which is considered to be at a very early stage of star formation. The detailed column density distribution on size scales from a few tens AU to ~10,000 AU scale are revealed by combining the ALMA (12 m array + 7 m array) data with the published/unpublished single-dish data. The high angular resolution observations at 0.87 mm with a synthesized beam size of ~0."74 x 0."32 reveal that a protostellar source, MMS-1, is not spatially resolved and lacks associated gas emission, while a starless high-density core, MMS-2, has substructures both in dust and molecular emission. The averaged radial column density distribution of the inner part of MC27/L1521F (r $\lesssim$ 3000 AU) is N(H2) ~r$^{-0.4}$, clearly flatter than that of the outer part, ~r$^{-1.0}$. The complex velocity/spatial structure obtained with previous ALMA observations is located inside the inner flatter region, which may reflect the dynamical status of the dense core.

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.

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.

We have observed the Class I protostar TMC-1A with Atacama Millimeter/submillimeter Array (ALMA) in 12CO and C18O (J=2-1), and 1.3-mm dust continuum emission. Continuum emission with a deconvolved size of 0.50"x0.37", perpendicular to the 12CO outflow, is detected. It most likely traces a circumstellar disk around TMC-1A, as previously reported. In contrast, a more extended structure is detected in C18O although it is still elongated with a deconvolved size of 3.3"x2.2", indicating that C18O traces mainly a flattened envelope surrounding the disk and the central protostar. C18O shows a clear velocity gradient perpendicular to the outflow at higher velocities, indicative of rotation, while an additional velocity gradient along the outflow is found at lower velocities. The radial profile of the rotational velocity is analyzed in detail, finding that it is given as a power-law \propto r^{-a} with an index of ~0.5 at higher velocities. This suggests that the rotation at higher velocities can be explained as Keplerian rotation orbiting a protostar with a dynamical mass of 0.68 Mo (inclination-corrected). The additional velocity gradient of C18O along the outflow is considered to be mainly infall motions in the envelope. Position-Velocity diagrams made from models consisting of an infalling envelope and a Keplerian disk are compared with the observations, revealing that the observed infall velocity is ~0.3 times smaller than the free fall velocity yielded by the dynamical mass of the protostar. Magnetic fields could be responsible for the slow infall velocity. A possible scenario of Keplerian disk formation is discussed.

The formation process of circumstellar disks is still controversial because of the interplay of complex physical processes that occurs during the gravitational collapse of prestellar cores. In this study, we investigate the effect of the Hall current term on the formation of the circumstellar disk using three- dimensional simulations. In our simulations, all non-ideal effects as well as the radiation transfer are considered. The size of the disk is significantly affected by a simple difference in the inherent properties of the prestellar core, namely whether the rotation vector and the magnetic field are parallel or anti- parallel. In the former case, only a very small disk (< 1 AU) is formed. On the other hand, in the latter case, a massive and large (> 20 AU) disk is formed in the early phase of protostar formation. Since the parallel and anti-parallel properties do not readily change, we expect that the parallel and anti-parallel properties are also important in the subsequent disk evolution and the difference between the two cases is maintained or enhanced. This result suggests that the disk size distribution of the Class 0 young stellar objects is bimodal. Thus, the disk evolution can be categorized into two cases and we may call the parallel and anti-parallel systems as Ortho-disk and Para-disk, respectively. We also show that the anti-rotating envelopes against the disk-rotation appear with a size of 200 AU. We predict that the anti-rotating envelope will be found in the future observations.

We investigate the formation and evolution of a first core, protostar, and circumstellar disc with a three-dimensional non-ideal (including both Ohmic and ambipolar diffusion) radiation magnetohydrodynamics simulation. We found that the magnetic flux is largely removed by magnetic diffusion in the first core phase and that the plasma $\beta$ of the centre of the first core becomes large, $\beta>10^4$. Thus, proper treatment of first core phase is crucial in investigating the formation of protostar and disc. On the other hand, in an ideal simulation, $\beta\sim 10$ at the centre of the first core. The simulations with magnetic diffusion show that the circumstellar disc forms at almost the same time of protostar formation even with a relatively strong initial magnetic field (the value for the initial mass-to-flux ratio of the cloud core relative to the critical value is $\mu=4$). The disc has a radius of $r \sim 1$ AU at the protostar formation epoch. We confirm that the disc is rotationally supported. We also show that the disc is massive ($Q\sim 1$) and that gravitational instability may play an important role in the subsequent disc evolution.

The main accretion phase of star formation is investigated in clouds with different metallicities in the range of 0 \le Z \le Z_\odot, resolving the protostellar radius. Starting from a near-equilibrium prestellar cloud, we calculate the cloud evolution up to \sim100 yr after the first protostar formation. The star formation process considerably differs between clouds with lower (Z \le 10^-4 Z_\odot) and higher (Z > 10^-4 Z_\odot) metallicities. Fragmentation frequently occurs and many protostars appear without forming a stable circumstellar disc in lower-metallicity clouds. In these clouds, although protostars mutually interact and some are ejected from the cloud centre, many remain as a small stellar cluster. In contrast, higher-metallicity clouds produce a single protostar surrounded by a nearly stable rotation-supported disc. In these clouds, although fragmentation occasionally occurs in the disc, the fragments migrate inwards and finally fall onto the central protostar. The difference in cloud evolution is due to different thermal evolutions and mass accretion rates. The thermal evolution of the cloud determines the emergence and lifetime of the first core. The first core develops prior to the protostar formation in higher-metallicity clouds, whereas no (obvious) first core appears in lower-metallicity clouds. The first core evolves into a circumstellar disc with a spiral pattern, which effectively transfers the angular momentum outwards and suppresses frequent fragmentation. In lower-metallicity clouds, the higher mass accretion rate increases the disc surface density within a very short time, rendering the disc unstable to self-gravity and inducing vigorous fragmentation.

The transport of angular momentum by magnetic fields is a crucial physical process in formation and evolution of stars and disks. Because the ionization degree in star forming clouds is extremely low, non-ideal magnetohydrodynamic (MHD) effects such as ambipolar diffusion and Ohmic dissipation work strongly during protostellar collapse. These effects have significant impacts in the early phase of star formation as they redistribute magnetic flux and suppress angular momentum transport by magnetic fields. We perform three-dimensional nested-grid radiation magnetohydrodynamic (RMHD) simulations including Ohmic dissipation and ambipolar diffusion. Without these effects, magnetic fields transport angular momentum so efficiently that no rotationally supported disk is formed even after the second collapse. Ohmic dissipation works only in a relatively high density region within the first core and suppresses angular momentum transport, enabling formation of a very small rotationally supported disk after the second collapse. With both Ohmic dissipation and ambipolar diffusion, these effects work effectively in almost the entire region within the first core and significant magnetic flux loss occurs. As a result, a rotationally supported disk is formed even before a protostellar core forms. The size of the disk is still small, about 5 AU at the end of the first core phase, but this disk will grow later as gas accretion continues. Thus the non-ideal MHD effects can resolve the so-called magnetic braking catastrophe while maintaining the disk size small in the early phase, which is implied from recent interferometric observations.

A protostellar jet and outflow are calculated for \sim 270 yr following the protostar formation using a three dimensional magnetohydrodynamics simulation, in which both the protostar and its parent cloud are spatially resolved. A high-velocity (\sim100km/s) jet with good collimation is driven near the disk's inner edge, while a low-velocity (<10km/s) outflow with a wide opening angle appears in the outer-disk region. The high-velocity jet propagates into the low-velocity outflow, forming a nested velocity structure in which a narrow high-velocity flow is enclosed by a wide low-velocity flow. The low-velocity outflow is in a nearly steady state, while the high-velocity jet appears intermittently. The time-variability of the jet is related to the episodic accretion from the disk onto the protostar, which is caused by gravitational instability and magnetic effects such as magnetic braking and magnetorotational instability. Although the high-velocity jet has a large kinetic energy, the mass and momentum of the jet are much smaller than those of the low-velocity outflow. A large fraction of the infalling gas is ejected by the low-velocity outflow. Thus, the low-velocity outflow actually has a more significant effect than the high-velocity jet in the very early phase of the star formation.

We investigate structure of self-gravitating disks, their fragmentation and evolution of the fragments (the clumps) using both analytic approach and three-dimensional radiation hydrodynamics simulations starting from molecular cores. The simulations show that non-local radiative transfer determines disk temperature. We find the disk structure is well described by an analytical model of quasi-steady self-gravitating disk with radial radiative transfer. Because the radiative process is not local and radiation from the interstellar medium cannot be ignored, the local radiative cooling would not be balanced with the viscous heating in a massive disk around a low mass star. In our simulations, there are cases in which the disk does not fragment even though it satisfies the fragmentation criterion based on disk cooling time ($Q \sim 1$ and $\Omega t_{\rm cool}\sim 1$). This indicates that at least the criterion is not sufficient condition for fragmentation. We determine the parameter range for the host cloud core in which disk fragmentation occurs. In addition, we show that the temperature evolution of the center of the clump is close to that of typical first cores and the minimum initial mass of clumps to be about a few Jupiter mass.

We report Atacama Large Millimeter/submillimeter Array (ALMA) cycle 0 observations of C$^{18}$O ($J=2-1$), SO ($J_N= 6_5-5_4$) and 1.3mm dust continuum toward L1527 IRS, a class 0 solar-type protostar surrounded by an infalling and rotating envelope. C$^{18}$O emission shows strong redshifted absorption against the bright continuum emission associated with L1527 IRS, strongly suggesting infall motions in the C$^{18}$O envelope. The C$^{18}$O envelope also rotates with a velocity mostly proportional to $r^{-1}$, where $r$ is the radius, while the rotation profile at the innermost radius (54 AU) may be shallower than $r^{-1}$, suggestive of formation of a Keplerian disk around the central protostar of 0.3 Mo in dynamical mass. SO emission arising from the inner part of the C$^{18}$O envelope also shows rotation in the same direction as the C$^{18}$O envelope. The rotation is, however, rigid-body like which is very different from the differential rotation shown by C$^{18}$O. In order to explain the line profiles and the position-velocity (PV) diagrams of C$^{18}$O and SO observed, simple models composed of an infalling envelope surrounding a Keplerian disk of 54 AU in radius orbiting a star of 0.3 Mo are examined. It is found that in order to reproduce characteristic features of the observed line profiles and PV diagrams, the infall velocity in the model has to be smaller than the free-fall velocity yielded by a star of 0.3 Mo. Possible reasons for the reduced infall velocities are discussed.

We have conducted ALMA observations in the 1.3 mm continuum and 12CO (2-1), C18O (2-1) and SO (5_6-4_5) lines toward L1489 IRS, a Class I protostar surrounded by a Keplerian disk and an infalling envelope. The Keplerian disk is clearly identified in the 12CO and C18O emission, and its outer radius (~700 AU) and mass (~0.005 Msun) are comparable to those of disks around T Tauri stars. The protostellar mass is estimated to be 1.6 Msun with the inclination angle of 66 deg. In addition to the Keplerian disk, there are blueshifted and redshifted off-axis protrusions seen in the C18O emission pointing toward the north and the south, respectively, adjunct to the middle part of the Keplerian disk. The shape and kinematics of these protrusions can be interpreted as streams of infalling flows with a conserved angular momentum following parabolic trajectories toward the Keplerian disk, and the mass infalling rate is estimated to be ~5E-7 Msun/yr. The specific angular momentum of the infalling flows (~2.5E-3 km/s*pc) is comparable to that at the outer radius of the Keplerian disk (~4.8E-3 km/s*pc). The SO emission is elongated along the disk major axis and exhibits a linear velocity gradient along the axis, which is interpreted as that the SO emission primarily traces a ring region in the flared Keplerian disk at radii of ~250-390 AU. The local enhancement of the SO abundance in the ring region can be due to the accretion shocks at the centrifugal radius where the infalling flows fall onto the disk. Our ALMA observations unveiled both the Keplerian disk and the infalling gas onto the disk, and the disk can further grow by accreting material and angular momenta from the infalling gas.

Starless dense cores eventually collapse dynamically, which forms protostars inside them, and the physical properties of the cores determine the nature of the forming protostars. We report ALMA observations of dust continuum emission and molecular rotational lines toward MC27 or L1521F, which is considered to be very close to the first protostellar core phase. We found a few starless high-density cores, one of which has a very high density of ~10^7 cm^-3, within a region of several hundred AU around a very low-luminosity protostar detected by Spitzer. A very compact bipolar outflow with a dynamical timescale of a few hundred years was found toward the protostar. The molecular line observation shows several cores with an arc-like structure, possibly due to the dynamical gas interaction. These complex structures revealed in the present observations suggest that the initial condition of star formation is highly dynamical in nature, which is considered to be a key factor in understanding fundamental issues of star formation such as the formation of multiple stars and the origin of the initial mass function of stars.

We investigate accretion of solid materials onto circumplanetary disks from heliocentric orbits rotating in protoplanetary disks, which is a key process for the formation of regular satellite systems. In the late stage of gas-capturing phase of giant planet formation, the accreting gas from protoplanetary disks forms circumplanetary disks. Since the accretion flow toward the circumplanetary disks affects the particle motion through gas drag force, we use hydrodynamic simulation data for the gas drag term to calculate the motion of solid materials. We consider wide range of size for the solid particles ($10^{-2}$-$10^6$m), and find that the accretion efficiency of the solid particles peaks around 10m-sized particles because energy dissipation of drag with circum-planetary disk gas in this size regime is most effective. The efficiency for particles larger than 10m size becomes lower because gas drag becomes less effective. For particles smaller than 10m, the efficiency is lower because the particles are strongly coupled with the back-ground gas flow, which prevent particles from accretion. We also find that the distance from the planet where the particles are captured by the circumplanetary disks is in a narrow range and well described as a function of the particle size.

We report the results of our three-dimensional radiation hydrodynamics simulation of collapsing unmagnetized molecular cloud cores. We investigate the formation and evolution of the circumstellar disk and the clumps formed by disk fragmentation. Our simulation shows that disk fragmentation occurs in the early phase of circumstellar disk evolution and many clumps form. The clump can be represented by a polytrope sphere of index $n \sim 3$ and $n \gtrsim 4$ at central temperature $T_c \lesssim100$ K and $T_c \gtrsim 100$ K, respectively. We demonstrate, numerically and theoretically, that the maximum mass of the clump, beyond which it inevitably collapses, is $\sim 0.03 M_{\odot}$. The entropy of the clump increases during its evolution, implying that evolution is chiefly determined by mass accretion from the disk rather than by radiative cooling. Although most of the clumps rapidly migrate inward and finally fall onto the protostar, a few clumps remain in the disk. The central density and temperature of the surviving clump rapidly increase and the clump undergoes a second collapse within 1000 - 2000 years after its formation. In our simulation, three second cores of masses $0.2\msun$, $0.15\msun$, and $0.06\msun$ formed. These are protostars or brown dwarfs rather than protoplanets. For the clumps to survive as planetary-mass objects, the rapid mass accretion should be prevented by some mechanisms.

The formation of Population III stars is investigated using resistive magnetohydrodynamic simulations. Starting from a magnetized primordial prestellar cloud, we calculate the cloud evolution several hundreds of years after first protostar formation, resolving the protostellar radius. When the natal minihalo field strength is weaker than B \lesssim 10^-13 (n/1 cm^-3)^-2/3 G (n is the hydrogen number density), magnetic effects can be ignored. In this case, fragmentation occurs frequently and a stellar cluster forms, in which stellar mergers and mass exchange between protostars contribute to the mass growth of these protostars. During the early gas accretion phase, the most massive protostar remains near the cloud centre, whereas some of the less massive protostars are ejected. The magnetic field significantly affects Population III star formation when B_amb \gtrsim 10^-12 (n/1 cm^-3)^-2/3 G. In this case, because the angular momentum around the protostar is effectively transferred by both magnetic braking and protostellar jets, the gas falls directly onto the protostar without forming a disk, and only a single massive star forms. In addition, a massive binary stellar system appears when B_amb \sim 10^-12 (n/1cm^-3)^-2/3 G. Therefore, the magnetic field determines the end result of the formation process (cluster, binary or single star) for Population III stars. Moreover, no persistent circumstellar disk appears around the protostar regardless of the magnetic field strength, which may influence the further evolution of Population III stars.

The formation of a circumstellar disk in collapsing cloud cores is investigated with three-dimensional magnetohydrodynamic simulations. We prepare four types of initial cloud having different density profiles and calculate their evolution with or without a sink. To investigate the effect of magnetic dissipation on disk formation, the Ohmic dissipation is considered in some models. Calculations show that disk formation is very sensitive to both the initial cloud configuration and the sink treatment. The disk size considerably differs in clouds with different density profiles even when the initial clouds have almost the same mass-to-flux ratio. Only a very small disk (\sim 10 AU in size) appears in clouds with a uniform density profile, whereas a large disk (\sim 100 AU in size) forms in clouds with a Bonnor-Ebert density profile. In addition, a large sink accretion radius numerically impedes disk formation during the main accretion phase and tends to foster the misleading notion that disk formation is completely suppressed by magnetic braking. The protostellar outflow is also greatly affected by the sink properties. A sink accretion radius of \lesssim 1 AU and sink threshold density of \gtrsim 10^ 13 cm^-3 are necessary for investigating disk formation during the main accretion phase.

We investigate the formation process of self-gravitating protoplanetary disks in unmagnetized molecular clouds. The angular momentum is redistributed by the action of gravitational torques in the massive disk during its early formation. We develop a simplified one-dimensional accretion disk model that takes into account the infall of gas from the envelope onto the disk and the transfer of angular momentum in the disk with an effective viscosity. First we evaluate the gas accretion rate from the cloud core onto the disk by approximately estimating the effects of gas pressure and gravity acting on the cloud core. We formulate the effective viscosity as a function of the Toomre Q parameter that measures the local gravitational stability of the rotating thin disk. We use a function for viscosity that changes sensitively with Q when the disk is gravitationally unstable. We find a strong self-regulation mechanism in the disk evolution. During the formation stage of protoplanetary disks, the evolution of the surface density does not depend on the other details of the modeling of effective viscosity, such as the prefactor of the viscosity coefficient. Next, to verify our model, we compare the time evolution of the disk calculated with our formulation with that of three-dimensional hydrodynamical simulations. The structures of the resultant disks from the one-dimensional accretion disk model agree well with those of the three-dimensional simulations. Our model is a useful tool for the further modeling of chemistry, radiative transfer, and planet formation in protoplanetary disks.

The evolution of protostellar outflow is investigated with resistive magneto-hydrodynamic nested-grid simulations that cover a wide range of spatial scales (\sim 1AU - 1pc). We follow cloud evolution from the pre-stellar core stage until the infalling envelope dissipates long after the protostar formation. We also calculate protostellar evolution to derive protostellar luminosity with time-dependent mass accretion through a circumstellar disk. The protostellar outflow is driven by the first core prior to protostar formation and is directly driven by the circumstellar disk after protostar formation. The opening angle of the outflow is large in the Class 0 stage. A large fraction of the cloud mass is ejected in this stage, which reduces the star formation efficiency to \sim 50%. After the outflow breaks out from the natal cloud, the outflow collimation is gradually improved in the Class I stage. The head of the outflow travels more than \sim 10^5AU in \sim 10^5 yr. The outflow momentum, energy and mass derived in our calculations agree well with observations. In addition, our simulations show the same correlations among outflow momentum flux, protostellar luminosity and envelope mass as those in observations. These correlations differ between Class 0 and I stages, which is explained by different evolutionary stages of the outflow; in the Class 0 stage, the outflow is powered by the accreting mass and acquires its momentum from the infalling envelope; in the Class I stage, the outflow enters the momentum-driven snow-plough phase. Our results suggest that protostellar outflow should determine the final stellar mass and significantly affect the early evolution of low-mass protostars.

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 report the first three-dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without Ohmic dissipation.We take into account many physical processes required to study star formation processes, including a realistic equation of state. We follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail.We can calculate only about one year after the formation of the protostellar cores with our direct three-dimensional RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully describe the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one-dimensional simulations using full radiation transfer, but differ considerably from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models because of efficient angular momentum transport by magnetic fields, but Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. The formed disks are still small (less than 0.35 AU) because we simulate only the earliest evolution. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive MHD models.

The configuration and evolution of the magnetic field in star-forming cores are investigated in order to directly compare simulations and observations. We prepare four different initial clouds having different magnetic field strengths and rotation rates, in which magnetic field lines are aligned/misaligned with the rotation axis. First, we calculate the evolution of such clouds from the prestellar stage until long after protostar formation. Then, we calculate the polarization of thermal dust emission expected from the simulation data. We create polarization maps with arbitrary viewing angles and compare them with observations. Using this procedure, we confirmed that the polarization distribution projected on the celestial plane strongly depends on the viewing angle of the cloud. Thus, by comparing the observations with the polarization map predicted by the simulations, we can roughly determine the angle between the direction of the global magnetic field and the line of sight. The configuration of the polarization vectors also depends on the viewing angle. We find that an hourglass configuration of magnetic field lines is not always realized in a collapsing cloud when the global magnetic field is misaligned with the cloud rotation axis. Depending on the viewing angle, an S-shaped configuration of the magnetic field (or the polarization vectors) appears early in the protostellar accretion phase. This indicates that not only the magnetic field but also the cloud rotation affects the dynamical evolution of such a cloud. In addition, by comparing the simulated polarization with actual observations, we can estimate properties of the host cloud such as the evolutionary stage, magnetic field strength, and rotation rate.

We investigate the formation and evolution of circumstellar disks in turbulent cloud cores until several 104 years after protostar formation using smoothed particle hydrodynamics (SPH) calculations. The formation and evolution process of circumstellar disk in turbulent cloud cores differs substantially from that in rigidly rotating cloud cores. In turbulent cloud cores, a filamentary structure appears before the protostar formation and the protostar forms in the filament. If the turbulence is initially sufficiently strong, the remaining filament twists around the protostar and directly becomes a rotation-supported disk. Upon formation, the disk orientation is generally misaligned with the angular momentum of its host cloud core and it dynamically varies during the main accretion phase, even though the turbulence is weak. This is because the angular momentum of the entire cloud core is mainly determined by the large scale velocity field whose wavelength is comparable to the cloud scale, whereas the angular momentum of the disk is determined by the local velocity field where the protostar forms and these two velocity fields do not correlate with each other. In the case of disk evolution in a binary or multiple stars, the disks are misaligned with each other at least during the main accretion phase, because there is no correlation between the velocity fields around the position where each protostar forms. In addition, each disk is also misaligned with the binary orbital plane. Such misalignment can explain the recent observations of misaligned disks and misaligned protostellar outflows.

We measured polarized dust emission at 350um towards the high-mass star forming massive dense clump IRAS 20126+4104 using the SHARC II Polarimeter, SHARP, at the Caltech Submillimeter Observatory. Most of the observed magnetic field vectors agree well with magnetic field vectors obtained from a numerical simulation for the case when the global magnetic field lines are inclined with respect to the rotation axis of the dense clump. The results of the numerical simulation show that rotation plays an important role on the evolution of the massive dense clump and its magnetic field. The direction of the cold CO 1-0 bipolar outflow is parallel to the observed magnetic field within the dense clump as well as the global magnetic field, as inferred from optical polarimetry data, indicating that the magnetic field also plays a critical role in an early stage of massive star formation. The large-scale Keplerian disk of the massive (proto)star rotates in almost opposite sense to the clump's envelope. The observed magnetic field morphology and the counter-rotating feature of the massive dense clump system provide hints to constrain the role of magnetic fields in the process of high mass star formation.

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 investigate gas accretion flow onto a circumplanetary disk from a protoplanetary disk in detail by using high-resolution three-dimensional nested-grid hydrodynamic simulations, in order to provide a basis of formation processes of satellites around giant planets. Based on detailed analyses of gas accretion flow, we find that most of gas accretion onto circumplanetary disks occurs nearly vertically toward the disk surface from high altitude, which generates a shock surface at several scale heights of the circumplanetary disk. The gas that has passed through the shock surface moves inward because its specific angular momentum is smaller than that of the local Keplerian rotation, while gas near the midplane in the protoplanetary disk cannot accrete to the circumplanetary disk. Gas near the midplane within the planet's Hill sphere spirals outward and escapes from the Hill sphere through the two Lagrangian points L$_1$ and L$_2$. We also analyze fluxes of accreting mass and angular momentum in detail and find that the distributions of the fluxes onto the disk surface are well described by power-law functions and that a large fraction of gas accretion occurs at the outer region of the disk, i.e., at about 0.1 times the Hill radius. The nature of power-law functions indicates that, other than the outer edge, there is no specific radius where gas accretion is concentrated. These source functions of mass and angular momentum in the circumplanetary disk would provide us with useful constraints on the structure and evolution of the circumplanetary disk, which is important for satellite formation.

Star formation efficiency controlled by the protostellar outflow in a single cloud core is investigated by three-dimensional resistive MHD simulations. Starting from the prestellar cloud core, the star formation process is calculated until the end of the main accretion phase. In the calculations, the mass of the prestellar cloud is parameterized. During the star formation, the protostellar outflow is driven by the circumstellar disk. The outflow extends also in the transverse direction until its width becomes comparable to the initial cloud scale, and thus, the outflow has a wide opening angle of >40 degrees. As a result, the protostellar outflow sweeps up a large fraction of the infalling material and ejects it into the interstellar space. The outflow can eject at most over half of the host cloud mass, significantly decreasing star formation efficiency. The outflow power is stronger in clouds with a greater initial mass. Thus, the protostellar outflow effectively suppresses star formation efficiency in a massive cloud. The outflow weakens significantly and disappears in several free-fall timescales of the initial cloud after the cloud begins to collapse. The natal prestellar core influences the lifetime and size of the outflow. At the end of the main accretion phase, a massive circumstellar disk comparable in mass to the protostar remains. Calculations show that typically, ~30% of the initial cloud mass is converted into the protostar and ~20% remains in the circumstellar disk, while ~40% is ejected into the interstellar space by the protostellar outflow. Therefore, a single cloud core typically has a star formation efficiency of 30-50%.

Using resistive magnetohydrodynamics simulation, we investigate circumstellar disk formation in a strongly magnetized cloud. As the initial state, an isolated cloud core embedded in a low-density interstellar medium with a uniform magnetic field is adopted. The cloud evolution is calculated until almost all gas inside the initial cloud falls onto either the circumstellar disk or a protostar, and a part of the gas is ejected into the interstellar medium by the protostellar outflow driven by the circumstellar disk. In the early main accretion phase, the disk size is limited to \sim 10 AU because the angular momentum of the circumstellar disk is effectively transferred by both magnetic braking and the protostellar outflow. In the later main accretion phase, however, the circumstellar disk grows rapidly and exceeds 100 AU by the end of the main accretion phase. This rapid growth of the circumstellar disk is caused by the depletion of the infalling envelope, while magnetic braking is effective when the infalling envelope is more massive than the circumstellar disk. The infalling envelope cannot brake the circumstellar disk when the latter is more massive than the former. In addition, the protostellar outflow weakens and disappears in the later main accretion phase, because the outflow is powered by gas accretion onto the circumstellar disk. Although the circumstellar disk formed in a magnetized cloud is considerably smaller than that in an unmagnetized cloud, a circumstellar disk exceeding 100 AU can form even in a strongly magnetized cloud.

We have carried out hydro-dynamical simulations to investigate the formation and evolution of protostar and circumstellar disks from the prestellar cloud. As the initial state, we adopt the molecular cloud core with two non-dimensional parameters representing the thermal and rotational energies. With these parameters, we make 17 models and calculate the cloud evolution ~ 10^4 years after the protostar formation. We find that early evolution of the star-disk system can be qualitatively classified into four modes: the massive disk, early fragmentation, late fragmentation, and protostar dominant modes. In the 'massive disk mode' to which the majority of models belong, the disk mass is greater than the protostellar mass for over 10^4 years and no fragmentation occurs in the circumstellar disk. The collapsing cloud shows fragmentation before the protostar formation in the 'early fragmentation mode'. The circumstellar disk shows fragmentation after the protostar formation in the 'late fragmentation mode', in which the secondary star substantially gains its mass from the circumstellar disk after fragmentation and it has a mass comparable to that of the primary star. The protostellar mass rapidly increases and exceeds the circumstellar disk mass in the 'protostar dominant mode'. This mode appears only when the initial molecular cloud core has a considerably small rotational energy. Comparison of our results with observations indicates that a majority of protostar has a fairly massive disk during the main accretion phase: the circumstellar disk mass is comparable to or more massive than the protostar. It is expected that such a massive disk promotes gas-giant formation by gravitational instability in a subsequent evolution stage.

A first adiabatic core is a transient object formed in the early phase of star formation. The observation of a first core is believed to be difficult because of its short lifetime and low luminosity. On the basis of radiation hydrodynamic simulations, we propose a novel theoretical model of first cores, Exposed Long-lifetime First core (ELF). In the very low-mass molecular core, the first core evolves slowly and lives longer than 10,000 years because the accretion rate is considerably low. The evolution of ELFs is different from that of ordinary first cores because radiation cooling has a significant effect there. We also carry out radiation transfer calculation of dust-continuum emission from ELFs to predict their observational properties. ELFs have slightly fainter but similar SEDs to ordinary first cores in radio wavelengths, therefore they can be observed. Although the probabilities that such low mass cores become gravitationally unstable and start to collapse are low, we still can expect that a considerable number of ELFs can be formed because there are many low-mass molecular cloud cores in star-forming regions that can be progenitors of ELFs.

The formation and evolution of a circumstellar disk in magnetized cloud cores is investigated from prestellar core stage until sim 10^4 yr after protostar formation. In the circumstellar disk, fragmentation first occurs due to gravitational instability in a magnetically inactive region, and substellar-mass objects appear. The substellar-mass objects lose their orbital angular momenta by gravitational interaction with the massive circumstellar disk and finally fall onto the protostar. After this fall, the circumstellar disk increases its mass by mass accretion and again induces fragmentation. The formation and falling of substellar-mass objects are repeated in the circumstellar disk until the end of the main accretion phase. In this process, the mass of fragments remain small, because the circumstellar disk loses its mass by fragmentation and subsequent falling of fragments before it becomes very massive. In addition, when fragments orbit near the protostar, they disturb the inner disk region and promote mass accretion onto the protostar. The orbital motion of substellar-mass objects clearly synchronizes with the time variation of the accretion luminosity of the protostar. Moreover, as the objects fall, the protostar shows a strong brightening for a short duration. The intermittent protostellar outflows are also driven by the circumstellar disk whose magnetic field lines are highly tangled owing to the orbital motion of fragments. The time-variable protostellar luminosity and intermittent outflows may be a clue for detecting planetary-mass objects in the circumstellar disk.

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.

The formation and evolution of the circumstellar disk in the collapsing molecular cloud is investigated from the prestellar stage resolving both the molecular cloud core and the protostar itself. In the collapsing cloud, the first adiabatic core appears prior to the protostar formation. Reflecting the thermodynamics of the collapsing gas, the first core is much more massive than the protostar. When the molecular cloud has no angular momentum, the first core falls onto the protostar and disappears a few years after the protostar formation. On the other hand, when the molecular cloud has an angular momentum, the first core does not disappear even after the protostar formation, and directly evolves into the circumstellar disk with a Keplerian rotation. There are two paths for the formation of the circumstellar disk. When the initial cloud has a considerably small rotational energy, two nested disks appear just after the protostar formation. During the early main accretion phase, the inner disk increases its size and merges with the outer disk (i.e. first core) to form a single circumstellar disk with a Keplerian rotation. On the other hand, when the molecular cloud has a rotational energy comparable to observations, a single Keplerian disk that corresponds to the first core already exists prior to the protostar formation. In such a cloud, the first core density gradually increases, maintaining the Keplerian rotation and forms the protostar inside it. Thus, the protostar is born in the Keplerian disk. In other words, a massive disk already exists before the protostar formation. In each case, the protostar at its formation is already surrounded by a massive circumstellar disk. The circumstellar disk is about 10-100 times more massive than the protostar in the main accretion disk. Such disks are favourable sites for the formation of binary companions and gas-giant planets.

In star forming regions, we can observe different evolutionary stages of various objects and phenomena such as molecular clouds, protostellar jets and outflows, circumstellar disks, and protostars. However, it is difficult to directly observe the star formation process itself, because it is veiled by the dense infalling envelope. Numerical simulations can unveil the star formation process in the collapsing gas cloud. Recently, some studies showed protostar formation from the prestellar core stage, in which both molecular clouds and protostars are resolved with sufficient spatial resolution. These simulations showed fragmentation and binary formation, outflow and jet driving, and circumstellar disk formation in the collapsing gas clouds. In addition, the angular momentum transfer and dissipation process of the magnetic field in the star formation process were investigated. In this paper, I briefly review recent developments in numerical simulations of low-mass star formation.

We use resistive magnetohydrodynamical simulations with the nested grid technique to study the formation of protoplanetary disks around protostars from molecular cloud cores that provide the realistic environments for planet formation. We find that gaseous planetary-mass objects are formed in the early evolutionary phase by gravitational instability in regions that are decoupled from the magnetic field and surrounded by the injection points of the magnetohydrodynamical outflows during the formation phase of protoplanetary disks. Magnetic decoupling enables massive disks to form and these are subject to gravitational instability, even at ~ 10 AU. The frequent formation of planetary-mass objects in the disk suggests the possibility of constructing a hybrid planet formation scenario, where the rocky planets form later under the influence of the giant planets in the protoplanetary disk.

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.

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.

The evolution of primordial collapsing clouds and formation of proto-Population III stars are investigated using three-dimensional ideal MHD simulation. We calculated the evolution of magnetized primordial clouds from the prestellar stage until the epoch after the proto-Population III star formation, spatially resolving both parsec-scale clouds and sub-AU scale protostars. The formation process of proto-population III star is characterized by the ratio of rotational to magnetic energy of the parent cloud. When the rotational energy is larger than the magnetic energy, fragmentation occurs in the collapsing primordial cloud before the proto-Population III star formation and binary or multiple system appears. Instead, when the magnetic energy is larger than the rotational energy, strong jet with >100km s^-1 is driven by circumstellar disk around the proto-population III star without fragmentation. Thus, even in the early universe, the magnetic field plays an important role in the star formation process.

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 investigate gas accretion onto a protoplanet, by considering the thermal effect of gas in three-dimensional hydrodynamical simulations, in which the wide region from a protoplanetary gas disk to a Jovian radius planet is resolved using the nested-grid method. We estimate the mass accretion rate and growth timescale of gas giant planets. The mass accretion rate increases with protoplanet mass for M_p < M_cri, while it becomes saturated or decreases for M_p > M_cri, where M_cri = 0.036 M_Jup (a_p/1AU)^0.75, and M_Jup and a_p are the Jovian mass and the orbital radius, respectively. The growth timescale of a gas giant planet or the timescale of the gas accretion onto the protoplanet is about 10^5 yr, that is two orders of magnitude shorter than the growth timescale of the solid core. The thermal effects barely affect the mass accretion rate because the gravitational energy dominates the thermal energy around the protoplanet. The mass accretion rate obtained in our local simulations agrees quantitatively well with those obtained in global simulations with coarser spatial resolution. The mass accretion rate is mainly determined by the protoplanet mass and the property of the protoplanetary disk. We find that the mass accretion rate is correctly calculated when the Hill or Bondi radius is sufficiently resolved. Using the oligarchic growth of protoplanets, we discuss the formation timescale of gas giant planets.

We perform a three-dimensional nested-grid radiation magneto-hydrodynamics (RMHD) simulation with self-gravity to study the early phase of the low-mass star formation process from a rotating molecular cloud core to a first adiabatic core just before the second collapse begins. Radiation transfer is handled with the flux-limited diffusion approximation, operator-splitting and implicit time-integrator. In the RMHD simulation, the outer region of the first core attains a higher entropy and the size of first core is larger than that in the magnetohydrodynamics simulations with the barotropic approximation. Bipolar molecular outflow consisting of two components is driven by magnetic Lorentz force via different mechanisms, and shock heating by the outflow is observed. Using the RMHD simulation we can predict and interpret the observed properties of star-forming clouds, first cores and outflows with millimeter/submillimeter radio interferometers, especially the Atacama Large Millimeter/submillimeter Array (ALMA).

The formation and evolution of the circumstellar disk in unmagnetized molecular clouds is investigated using three-dimensional hydrodynamic simulations from the prestellar core until the end of the main accretion phase. In collapsing clouds, the first (adiabatic) core with a size of ~10AU forms prior to the formation of the protostar. At its formation, the first core has a thick disk-like structure, and is mainly supported by the thermal pressure. After the protostar formation, it decreases the thickness gradually, and becomes supported by the centrifugal force. We found that the first core is a precursor of the circumstellar disk. This indicates that the circumstellar disk is formed before the protostar formation with a size of ~10AU, which means that no protoplanetary disk smaller than <10AU exists. Reflecting the thermodynamics of the collapsing gas, at the protostar formation epoch, the circumstellar disk has a mass of ~0.01-0.1 solar mass, while the protostar has a mass of ~10^-3 solar mass. Thus, just after the protostar formation, the circumstellar disk is about 10-100 times more massive than the protostar. Even in the main accretion phase that lasts for ~10^5yr, the circumstellar disk mass dominates the protostellar mass. Such a massive disk is unstable to gravitational instability, and tends to show fragmentation. Our calculations indicate that the planet or brown-dwarf mass object may form in the circumstellar disk in the main accretion phase. In addition, the mass accretion rate onto the protostar shows strong time variability that is caused by the perturbation of proto-planets and/or the spiral arms in the circumstellar disk. Such variability provides a useful signature for detecting the planet-sized companion in the circumstellar disk around very young protostars.

Star formation in relic HII regions of the first stars is investigated using magneto-hydrodynamical simulations with a nested grid method that covers 10 orders of magnitude in spatial scale and 20 orders of magnitude in density contrast. Due to larger fraction of H_2 and HD molecules, its prestellar thermal evolution is considerably different from that in the first star formation. Reflecting the difference, two hydrostatic cores appear in a nested manner: a protostar is enclosed by a transient hydrostatic core, which appears during the prestellar collapse. If the initial natal core rotates fast at a rate with rotational to gravitational energy ratio \beta_0 > 0.01-0.1, the transient hydrostatic core fragments to \sim 10 M_\odot sub-cores at density \sim 10^9 cm^-3. With smaller rotation energy, fragmentation occurs at higher density while a single protostar forms without fragmentation if rotation is extremely slow with \beta_0 < 10^-6 -10^-5. If magnetic field is present, these threshold values of \beta_0 is boosted owing to angular momentum transport by the magnetic breaking. Magnetic field also drives the protostellar outflows. With strong magnetic field, two distinct outflows are observed: The slower one emanates from the transient hydrostatic core, while the faster one from the protostar. These flows may affect the final stellar mass by ejecting some of masses in the initial core, and also may play some role in driving and maintenance of interstellar turbulence in young galaxies.

We examine emission from a young protostellar object (YPO) with three-dimensional ideal MHD simulations and three-dimensional non-local thermodynamic equilibrium (non-LTE) line transfer calculations, and show the first results. To calculate the emission field, we employed a snapshot result of an MHD simulation having young bipolar outflows and a dense protostellar disk (a young circumstellar disk) embedded in an infalling envelope. Synthesized line emission of two molecular species (CO and SiO) show that subthermally excited SiO lines as a high density tracer can provide a better probe of the complex velocity field of a YPO, compared to fully thermalized CO lines. In a YPO at the earliest stage when the outflows are still embedded in the collapsing envelope, infall, rotation and outflow motions have similar speeds. We find that the combined velocity field of these components introduces a great complexity in the line emissions through varying optical thickness and emissivity, such as asymmetric double-horn profiles. We show that the rotation of the outflows, one of the features that characterizes an outflow driven by magneto-centrifugal forces, appears clearly in velocity channel maps and intensity-weighted mean velocity (first moment of velocity) maps. The somewhat irregular morphology of the line emission at this youngest stage is dissimilar to a more evolved object such as young Class 0. High angular resolution observation by e.g., the Atacama Large Millimeter/submillimeter Array (ALMA) telescope can reveal these features.

We investigate the driving mechanism of outflows and jets in star formation process using resistive MHD nested grid simulations. We found two distinct flows in the collapsing cloud core: Low-velocity outflows (sim 5 km/s) with a wide opening angle, driven from the first adiabatic core, and high-velocity jets (sim 50 km/s) with good collimation, driven from the protostar. High-velocity jets are enclosed by low-velocity outflow. The difference in the degree of collimation between the two flows is caused by the strength of the magnetic field and configuration of the magnetic field lines. The magnetic field around an adiabatic core is strong and has an hourglass configuration. Therefore, the low-velocity outflow from the adiabatic core are driven mainly by the magnetocentrifugal mechanism and guided by the hourglass-like field lines. In contrast, the magnetic field around the protostar is weak and has a straight configuration owing to Ohmic dissipation in the high-density gas region. Therefore, high-velocity jet from the protostar are driven mainly by the magnetic pressure gradient force and guided by straight field lines. Differing depth of the gravitational potential between the adiabatic core and the protostar cause the difference of the flow speed. Low-velocity outflows correspond to the observed molecular outflows, while high-velocity jets correspond to the observed optical jets. We suggest that the protostellar outflow and the jet are driven by different cores (the first adiabatic core and protostar), rather than that the outflow being entrained by the jet.

The protostellar outflows have indispensable role in the formation of single stars, because they carry off the excess angular momentum from the centre of the shrinking gas cloud, and permits further collapse to form a star. On the other hand, a significant fraction of stars is supposed to be born as binaries with circumbinary disk that are frequently observed. Here, we investigate the evolution of a magnetized rotating cloud using three-dimensional resistive MHD nested-grid code, and show that the outflow is driven by the circumbinary disk and has an important role even in the binary formation. After the adiabatic core formation in the collapsing cloud core, the magnetic flux is significantly removed from the centre of the cloud by the Ohmic dissipation. Since this removal makes the magnetic braking ineffective, the adiabatic core continuously acquires the angular momentum to induce fragmentation and subsequent binary formation. The magnetic field accumulates in the circumbinary disk where the removal and accretion of magnetic field are balanced, and finally drives circumbinary outflow. This result explains the spectacular morphology of some specific young stellar objects such as L1551 IRS5. We can infer that most of the bipolar molecular outflows observed by low density tracers (i.e., CO) would correspond to circumbinary or circum-multiple outflows found in this report, since most of the young stellar objects are supposed to be binaries or multiples.

The fragmentation process in collapsing clouds with various metallicities is studied using three-dimensional nested-grid hydrodynamics. Initial clouds are specified by three parameters: cloud metallicity, initial rotation energy and initial cloud shape. For different combinations of these parameters, we calculate 480 models in total and study cloud evolution, fragmentation conditions, orbital separation and binary frequency. For the cloud to fragment during collapse, the initial angular momentum must be higher than a threshold value, which decreases with decreasing metallicity. Although the exact fragmentation conditions depend also on the initial cloud shape, this dependence is only modest. Our results indicate a higher binary frequency in lower-metallicity gas. In particular, with the same median rotation parameter as in the solar neighbourhood, a majority of stars are born as members of binary/multiple systems for < 10^-4 Z_sun. With initial mass < 0.1 M_sun, if fragments are ejected in embryo from the host clouds by multi-body interaction, they evolve to substellar-mass objects. This provides a formation channel for low-mass stars in zero- or low-metallicity environments.

Brown dwarf formation and star formation efficiency are studied using a nested grid simulation that covers five orders of magnitude in spatial scale (10^4 - 0.1AU). Starting with a rotating magnetized compact cloud with a mass of 0.22 M_sun, we follow the cloud evolution until the end of main accretion phase. Outflow of about 5 km/s emerges about 100 yr before the protostar formation and does not disappear until the end of the calculation. The mass accretion rate declines from 10^-6 M_sun/yr to 10^-8 - 10^-12 M_sun/yr in a short time (about 10^4 yr) after the protostar formation. This is because (1) a large fraction of mass is ejected from the host cloud by the protostellar outflow and (2) the gas escapes from the host cloud by the thermal pressure. At the end of the calculation, 74% (167 M_Jup) of the total mass (225 M_Jup) is outflowing from the protostar, in which 34% (77 M_Jup) of the total mass is ejected by the protostellar outflow with supersonic velocity and 40% (90 M_Jup) escapes with subsonic velocity. On the other hand, 20% (45 M_Jup) is converted into the protostar and 6% (13 M_Jup) remains as the circumstellar disk. Thus, the star formation efficiency is epsilon = 0.2. The resultant protostellar mass is in the mass range of brown dwarfs. Our results indicate that brown dwarfs can be formed in compact cores in the same manner as hydrogen-burning stars, and the magnetic field and protostellar outflow are essential in determining the star formation efficiency and stellar mass.

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$.

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 construct a new family of analytic models of black hole accretion disks in dynamical equilibria. Our construction is based on assuming distributions of angular momentum and entropy. For a particular choice of the distribution of angular momentum, we calculate the shapes of equipressure surfaces. The equipressure surfaces we find are similar to those in thick, slim and thin disks, and to those in ADAFs.

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.

The formation of a circumplanetary disk and accretion of angular momentum onto a protoplanetary system are investigated using 3D hydrodynamical simulations. The local region around a protoplanet in a protoplanetary disk is considered with sufficient spatial resolution: the region from outside the Hill sphere to the Jovian radius is covered by the nested-grid method. To investigate the thermal effects of the circumplanetary disk, various equations of state are adopted. Large thermal energy around the protoplanet slightly changes the structure of the circumplanetary disk. Compared with a model adopting an isothermal equation of state, in a model with an adiabatic equation of state, the protoplanet's gas envelope extends farther, and a slightly thick disk appears near the protoplanet. However, different equations of state do not affect the acquisition process of angular momentum for the protoplanetary system. Thus, the specific angular momentum acquired by the system is fitted as a function only of the protoplanet's mass. A large fraction of the total angular momentum contributes to the formation of the circumplanetary disk. The disk forms only in a compact region in very close proximity to the protoplanet. Adapting the results to the solar system, the proto-Jupiter and Saturn have compact disks in the region of r < 21 r_J (r < 0.028 r_HJ) and r < 66 r_S (r < 0.061 r_HS), respectively, where r_J (r_HJ) and r_S (r_HS) are the Jovian and Saturnian (Hill) radius, respectively. The surface density has a peak in these regions due to the balance between centrifugal force and gravity of the protoplanet. The size of these disks corresponds well to the outermost orbit of regular satellites around Jupiter and Saturn. Regular satellites may form in such compact disks around proto-gas giant planets.

Jet driving and fragmentation process in collapsing primordial cloud are studied using three-dimensional MHD nested grid simulations. Starting from a rotating magnetized spherical cloud with the number density of n=10^3 cm^-3, we follow the evolution of the cloud up to the stellar density n=10^22 cm^-3. We calculate 36 models parameterizing the initial magnetic and rotational energies (\gamma_0, \beta_0). In the collapsing primordial clouds, the cloud evolutions are characterized by the ratio of the initial rotational to magnetic energy, \gamma_0/\beta_0. The Lorentz force significantly affects the cloud evolution when \gamma_0 > \beta_0, while the centrifugal force is more dominant than the Lorentz force when \beta_0 > \gamma_0. When the cloud rotates rapidly with angular velocity of \Omega_0 > 10^-17 (n/10^3 cm^-3)^2/3 s^-1 and \beta_0 > \gamma_0, fragmentation occurs before the protostar is formed, but no jet appears after the protostar formation. On the other hand, a strong jet appears after the protostar formation without fragmentation when the initial cloud has the magnetic field of B_0 > 10^-9 (n/10^3 cm^-3)^2/3 G and \gamma_0 > \beta_0. Our results indicate that proto-Population III stars frequently show fragmentation and protostellar jet. Population III stars are therefore born as binary or multiple stellar systems, and they can drive strong jets, which disturb the interstellar medium significantly, as well as in the present-day star formation, and thus they may induce the formation of next generation stars.

Cloud evolution for various metallicities is investigated by three-dimensional nested grid simulations, in which the initial ratio of rotational to gravitational energy of the host cloud \beta_0 (=10^-1 - 10^-6) and cloud metallicity Z (=0 - Z_\odot) are parameters. Starting from a central number density of n = 10^4 cm^-3, cloud evolution for 48 models is calculated until the protostar is formed (n \simeq 10^23 cm^-3) or fragmentation occurs. The fragmentation condition depends both on the initial rotational energy and cloud metallicity. Cloud rotation promotes fragmentation, while fragmentation tends to be suppressed in clouds with higher metallicity. Fragmentation occurs when \beta_0 > 10^-3 in clouds with solar metallicity, while fragmentation occurs when \beta_0 > 10^-5 in the primordial gas cloud. Clouds with lower metallicity have larger probability of fragmentation, which indicates that the binary frequency is a decreasing function of cloud metallicity. Thus, the binary frequency at the early universe (or lower metallicity environment) is higher than at present day (or higher metallicity environment). In addition, binary stars born from low-metallicity clouds have shorter orbital periods than those from high-metallicity clouds. These trends are explained in terms of the thermal history of the collapsing cloud.

We study the effect of poloidal magnetic field on type I planetary migration by linear perturbation analysis in the shearing-sheet approximation and the analytic results are compared with numerical calculations. In contrast to the unmagnetized case, the basic equations that describe the wake due to the planet in the disk allow magnetic resonances at which density perturbation diverges. In order to simplify the problem, we consider the case without magneto-rotational instability. We perform two sets of analyses: two-dimensional and three-dimensional. In two-dimensional analysis, we find the generalization of the torque formula previously known in unmagnetized case. In three-dimensional calculations, we focus on the disk with very strong magnetic field and derive a new analytic formula for the torque exerted on the planet. We find that when Alfven velocity is much larger than sound speed, two-dimensional torque is suppressed and three-dimensional modes dominate, in contrast to the unmagnetized case.

We investigate the accretion of angular momentum onto a protoplanet system using three-dimensional hydrodynamical simulations. We consider a local region around a protoplanet in a protoplanetary disk with sufficient spatial resolution. We describe the structure of the gas flow onto and around the protoplanet in detail. We find that the gas flows onto the protoplanet system in the vertical direction crossing the shock front near the Hill radius of the protoplanet, which is qualitatively different from the picture established by two-dimensional simulations. The specific angular momentum of the gas accreted by the protoplanet system increases with the protoplanet mass. At Jovian orbit, when the protoplanet mass M_p is M_p < 1 M_J, where M_J is Jovian mass, the specific angular momentum increases as j \propto M_p. On the other hand, it increases as j \propto M_p^2/3 when the protoplanet mass is M_p > 1 M_J. The stronger dependence of the specific angular momentum on the protoplanet mass for M_p < 1 M_J is due to thermal pressure of the gas. The estimated total angular momentum of a system of a gas giant planet and a circumplanetary disk is two-orders of magnitude larger than those of the present gas giant planets in the solar system. A large fraction of the total angular momentum contributes to the formation of the circumplanetary disk. We also discuss the satellite formation from the circumplanetary disk.

The fragmentation process of primordial-gas cores during prestellar collapse is studied using three-dimensional nested-grid hydrodynamics. Starting from the initial central number density of n \sim10^3 cm^-3, we follow the evolution of rotating spherical cores up to the stellar density n \simeq 10^{22} cm^-3. An initial condition of the cores is specified by three parameters: the ratios of the rotation and thermal energies to the gravitational energy (\beta_0, and \alpha_0, respectively), and the amplitude of the bar-mode density perturbation (A_\phi). Cores with rotation \beta_0 > 10^{-6} are found to fragment during the collapse. The fragmentation condition hardly depends on either the initial thermal energy \alpha_0 or amplitude of bar-mode perturbation A_\phi. Since the critical rotation parameter for fragmentation is lower than that expected in the first star formation, binaries or multiples are also common for the first stars.

Fragmentation and binary formation processes are studied using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n=10^4 cm^-3) to the stellar core (n \simeq 10^22 cm^-3). We calculated 147 models with different initial magnetic, rotational, and thermal energies, and the amplitudes of the non-axisymmetric perturbation. In a collapsing cloud, fragmentation is mainly controlled by the initial ratio of the rotational to the magnetic energy, regardless of the initial thermal energy and amplitude of the non-axisymmetric perturbation. When the clouds have large rotational energies in relation to magnetic energies, fragmentation occurs in the low-density evolution phase (10^12 cm^-3 < n < 10^15 cm^-3) with separations of 3-300 AU. Fragments that appeared in this phase are expected to evolve into wide binary systems. On the other hand, fragmentation does not occur in the low-density evolution phase, when initial clouds have large magnetic energies in relation to the rotational energies. In these clouds, fragmentation only occurs in the high-density evolution phase (n > 10^17 cm^-3) after the clouds experience significant reduction of the magnetic field owing to Ohmic dissipation in the period of 10^12 cm^-3 < n < 10^15 cm^-3. Fragments appearing in this phase have separations of < 0.3 AU, and are expected to evolve into close binary systems. As a result, we found two typical fragmentation epochs, which cause different stellar separations. Although these typical separations are disturbed in the subsequent gas accretion phase, we might be able to observe two peaks of binary separations in extremely young stellar groups.

The evolution of the magnetic field and angular momentum in the collapsing cloud core is studied using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n=10^4 cm^-3) to the stellar core (n \simeq 10^22 cm^-3). The magnetic field strengths at the center of the clouds converge to a certain value as the clouds collapse, when the clouds have the same angular momenta but different strengths of the magnetic fields at the initial state. For 10^12 cm^-3 < n < 10^16 cm^-3, Ohmic dissipation considerably removes the magnetic field from the collapsing cloud core, and the magnetic field lines, which are strongly twisted for n <10^12 cm^-3, is de-collimated. The magnetic field lines are twisted and amplified again for nc > 10^16 cm^-3, because the magnetic field is recoupled with the warm gas. Finally, protostars at their formation epoch have 0.1-1kG of the magnetic fields, which are comparable to observations. The magnetic field strength of protostar slightly depends on the angular momentum of the host cloud. The protostar formed from the slowly rotating cloud core has a stronger magnetic field. The evolution of the angular momentum is closely related to the evolution of the magnetic field. The angular momentum in the collapsing cloud is removed by the magnetic effect. The formed protostars have 0.1-2 days of the rotation period at their formation epoch, which are slightly shorter than the observation. This indicates that the further removal mechanism of the angular momentum such as interaction between the protostar and disk, wind gas or jet is important in further evolution of the protostar.

The driving mechanism of jets and outflows in star formation process is studied using resistive MHD nested grid simulations. We calculated cloud evolution from the molecular cloud core to the stellar core. In the collapsing cloud core, we found two distinct flows: Low-velocity flows (sim 5 km/s) with a wide opening angle, driven from the adiabatic core, and high-velocity flows (sim 30 km/s) with good collimation, driven from the protostar. High-velocity flows are enclosed by low-velocity flows after protostar formation. The difference in the degree of collimation between the two flows is caused by the strength of the magnetic field and configuration of the magnetic field lines. The magnetic field around an adiabatic core is strong and has an hourglass configuration; therefore, flows from the adiabatic core are driven mainly by the magnetocentrifugal mechanism and guided by the hourglass-like field lines. In contrast, the magnetic field around the protostar is weak and has a straight configuration owing to Ohmic dissipation in the high-density gas region. Therefore, flows from the protostar are driven mainly by the magnetic pressure gradient force and guided by straight field lines. Differing depth of the gravitational potential between the adiabatic core and the protostar cause the difference of the flow speed. Low-velocity flows correspond to the observed molecular outflows, while high-velocity flows correspond to the observed optical jets. We suggest that the outflow and the jet are driven by different cores, rather than that the outflow being entrained by the jet.

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 investigate outflows driven by a giant protoplanet using three-dimensional MHD nested grid simulations. We consider a local region around the protoplanet in the protoplanetary disk, and calculate three models: (a) unmagnetized disk model, (b) magnetized disk model having magnetic field azimuthally parallel to the disk, and (c) magnetic field perpendicular to the disk. Outflows with velocities, at least, 10 km/s are driven by the protoplanets in both magnetized disk models, while outflow does not appear in unmagnetized disk model. Tube-like outflows along the azimuthal direction of the protoplanetary disk appear in model with magnetic field being parallel to the disk. In this model, the magnetically dominated regions (i.e., density gap) are clearly contrasted from other regions and spiral waves appear near the protoplanet. On the other hand, in model with magnetic field being perpendicular to the disk, outflows are driven by a protoplanet with cone-like structure just as seen in the outflow driven by a protostar. Magnetic field lines are strongly twisted near the protoplanet and the outflows have well-collimated structures in this model.These outflows can be landmarks for searching exo-protoplanets in their formation stages. Our results indicate that the accretion rate onto the protoplanet tend to have a larger value than that expected from previous hydrodynamical calculations, since a fraction of the angular momentum of circum-planetary disk is removed by outflows, enhanced non-axisymmetric patterns caused by magnetic field, and magnetic braking. Possible implications for observation are also briefly discussed.

The protostellar jets driven by the formation of the first stars are studied by using MHD nested grid simulations. Starting from a slowly rotating spherical cloud of 5.1 times 10^4 Msun permeated by a uniform magnetic field, we follow the evolution from the central number density n = 10^3 cm^-3 to n simeq 10^23 cm^-3. We calculate four models that differ in initial magnetic field strengths and angular velocities. In all models, protostars of simeq 10^-3 Msun are formed at n simeq 10^22 cm^-3 in accordance with one-dimensional calculations. By this epoch, the magnetic flux density is amplified by 10 orders of magnitude from the initial value. Consequently, the formed protostar possesses the magnetic field of \sim 10^6 G that is much larger than the flux density of the present counterparts, reflecting the fact that the dissipation of a magnetic field is ineffective in primordial gas clouds. If the initial magnetic field B > 10^-9 (n/10^3 cm^-3)^2/3 G, the protostellar jet is launched and its velocities reaches \sim 70 km s^-1 by the time the protostellar mass becomes (4-6) times 10^-3 Msun, and a fraction (3-10%) of the accreting matter is blown off from the central region. Owing to the interaction of these ejecta with surrounding matter, expanding bow shocks are created at both heads of the jet. If this jet continues to sweep out the surrounding gas that otherwise accretes onto the central star or circumstellar disk, the final mass of the first star can be substantially reduced. In addition, dense post-shock regions behind the bow shocks are expected to promote the chemical reactions (formation of H_2 and HD), and this provides possible environments for subsequent low-mass star formation in the early universe.

The stellar core formation and high speed jets driven by the formed core are studied by using three-dimensional resistive MHD nested grid simulations. Starting with a Bonnor-Ebert isothermal cloud rotating in a uniform magnetic field, we calculate the cloud evolution from the molecular cloud core (n = 10^6 cm^-3, r_c = 4.6 times 10^4 AU) to the stellar core (n \simeq 10^23 cm^-3, r_c \simeq 1 solar radius). We resolve cloud structure over 7 orders of magnitude in spatial extent and over 17 orders of magnitude in density contrast. For comparison, we calculate two models: resistive and ideal MHD models. Both models have the same initial condition, but the former includes dissipation process of magnetic field while the latter does not. The magnetic fluxes in resistive MHD model are extracted from the first core during 10^12 cm^-3 < n < 10^16 cm^-3 by Ohmic dissipation. Magnetic flux density of the formed stellar core (n \simeq 10^20 cm^-3) in resistive MHD model is two orders of magnitude smaller than that in ideal MHD model. Since magnetic braking is less effective in resistive MHD model, rapidly rotating stellar core (the second core) is formed. After stellar core formation, the magnetic field of the core is largely amplified both by magneto-rotational instability and the shearing motion between the stellar core and ambient medium. As a consequence, high speed (simeq 45 km,s^-1) jets are driven by the second core, which results in strong mass ejection. A cocoon-like structure around the second core also forms with clear bow shocks.

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 studied the collapse of rotating molecular cloud cores with inclined magnetic fields, based on three-dimensional numerical simulations.The numerical simulations start from a rotating Bonnor-Ebert isothermal cloud in a uniform magnetic field. The magnetic field is initially taken to be inclined from the rotation axis. As the cloud collapses, the magnetic field and rotation axis change their directions. When the rotation is slow and the magnetic field is relatively strong, the direction of the rotation axis changes to align with the magnetic field, as shown earlier by Matsumoto & Tomisaka. When the magnetic field is weak and the rotation is relatively fast, the magnetic field inclines to become perpendicular to the rotation axis. In other words, the evolution of the magnetic field and rotation axis depends on the relative strength of the rotation and magnetic field. Magnetic braking acts to align the rotation axis and magnetic field, while the rotation causes the magnetic field to incline through dynamo action. The latter effect dominates the former when the ratio of the angular velocity to the magnetic field is larger than a critical value \Omega_0/ B_0 > 0.39 G^1/2 c_s^-1, where B_0, \Omega_0, G, and c_s^-1 denote the initial magnetic field, initial angular velocity, gravitational constant, and sound speed, respectively. When the rotation is relatively strong, the collapsing cloud forms a disk perpendicular to the rotation axis and the magnetic field becomes nearly parallel to the disk surface in the high density region. A spiral structure appears due to the rotation and the wound-up magnetic field in the disk.

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.

Subsequent to Paper I, the evolution and fragmentation of a rotating magnetized cloud are studied with use of three-dimensional MHD nested-grid simulations. After the isothermal runaway collapse, an adiabatic gas forms a protostellar first core at the center of the cloud. When the isothermal gas is stable for fragmentation in a contracting disk, the adiabatic core often breaks into several fragments. Conditions for fragmentation and binary formation are studied. All the cores which show fragmentations are geometrically thin, as the diameter-to-thickness ratio is larger than 3. Two patterns of fragmentation are found. (1) When a thin disk is supported by centrifugal force, the disk fragments through a ring configuration (ring fragmentation). This is realized in a fast rotating adiabatic core as Omega >0.2 tau_ff^-1, where Omega and tau_ff represent the angular rotation speed and the free-fall time of the core, respectively. (2) On the other hand, the disk is deformed to an elongated bar in the isothermal stage for a strongly magnetized or rapidly rotating cloud. The bar breaks into 2 - 4 fragments (bar fragmentation). Even if a disk is thin, the disk dominated by the magnetic force or thermal pressure is stable and forms a single compact body. In either ring or bar fragmentation mode, the fragments contract and a pair of outflows are ejected from the vicinities of the compact cores. The orbital angular momentum is larger than the spin angular momentum in the ring fragmentation. On the other hand, fragments often quickly merge in the bar fragmentation, since the orbital angular momentum is smaller than the spin angular momentum in this case. Comparison with observations is also shown.

We discuss evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of > 5 times 10 ^10 cm^-3 . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches 1.5 (n/5 times 10^10 cm^-3)^1/2 mG, while if the core is rotationally dominated the angular velocity approaches 2.57 times 10^-3, (n/5 times 10^10 cm^-3)^1/2 yr^-1, where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds 5 times 10^10 cm^-3. Fragmentation of the very dense core and emergence of outflows from fragments are shown in the subsequent paper.

The evolution of a gas shell, swept by the supernova remnant of a massive first generation star, is studied with H_2 and HD chemistry taken into account. When a first-generation star explodes as a supernova, H_2 and HD molecules are formed in the swept gas shell and effectively cool the gas shell to temperatures of 32 K - 154 K. If the supernova remnant can sweep to gather the ambient gas, the gas shell comes to be dominated by its self-gravity, and hence, is expected to fragment. Our result shows that for a reasonable range of temperatures (200 K - 1000 K) of interstellar gas, the formation of second-generation stars can be triggered by a single supernova or hypernova.

We discuss the origin of HE0107-5240 which is the most metal poor star yet observed ([Fe/H] = -5.3). Its discovery has an important bearing on the question of the observability of "first generation" stars. In common with other metal-poor stars (-4 < [Fe/H] < -2.5), HE0107-5240 shows a peculiar abundance pattern (CNO rich, moderate Na rich). We elaborate the binary scenario on the basis of the evolution and nucleosynthesis of extremely metal-poor, low-mass model stars and discuss the possibility of discriminating this scenario from others. In our picture, iron peak elements arise in surface layers of the component stars by accretion of gas from the polluted primordial cloud. To explain the observed C, N, O, and Na enhancements as well as 12C/13C ratio, we suppose that the currently observed star accreted matter from a AGB companion. To estimate the abundances in the matter transferred in the binary, we rely on the results of computations of model stars constructed with up-to-date input physics. Nucleosynthesis in the helium flash convection with hydrogen injected is followed, allowing us to discuss the abundances of s-process elements, in addition to explaining the origin of the observed O and Na enrichments. From the observed abundances, we conclude that HE0107-5240 has evolved from a wide binary with a primary of initial mass, 1.2 ~ 3 Msun. We estimated the present binary separation of ~ 34 AU and period of ~ 150 years. We also conclude that the abundance distribution of the heavy s-process elements, may hold the key to understand the origin of HE0107-5240. An enhancement of [Pb/Fe] = 1 ~ 2 should be observed. If the enhancement of s-process elements is not detected, HE0107-5240 may be a first generation secondary in a binary system with a primary of mass less than 2.5 Msun.

This is the first paper about the fragmentation and mass outflow in the molecular cloud by using three-dimensional MHD nested-grid simulations. The binary star formation process is studied paying particular attention to the fragmentation of a rotating magnetized molecular cloud. We assume an isothermal rotating and magnetized cylindrical cloud in hydrostatic balance. Non-axisymmetric as well as axisymmetric perturbations are added to the initial state and the subsequent evolutions are studied. The evolution is characterized by three parameters: the amplitude of the non-axisymmetric perturbations, the rotation speed, and the magnetic field strength. As a result, it is found that non-axisymmetry hardly evolves in the early phase, but begins to grow after the gas contracts and forms a thin disk. Disk formation is strongly promoted by the rotation speed and the magnetic field strength. There are two types of fragmentation: fragmentation from a ring and that from a bar. Thin adiabatic cores fragment if a thickness is smaller than 1/4 of the radius. For the fragments to survive, they should be formed in a heavily elongated barred core or a flat round disk. In the models showing fragmentation, outflows from respective fragments are found as well as those driven by the rotating bar or the disk.

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.

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.

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.