We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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 magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $\theta_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of BHR 71 was found to be $B_{\rm pos} = 8.8 - 15.0$ $\mu$G, indicating that the BHR 71 core is magnetically supercritical with $\lambda = 1.44 - 2.43$. Taking into account the effect of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of BHR 71 was $M_{\rm cr} = 14.5-18.7$ M$_{\odot}$, which is consistent with the observed core mass of $M_{\rm core} \approx 14.7$ M$_{\odot}$ (Yang et al. 2017). We conclude that BHR 71 is in a condition close to a kinematically critical state, and the magnetic field direction lies close to the plane of sky. Since BHR 71 is a star-forming core, a significantly subcritical condition (i.e., the magnetic field direction deviating from the plane of sky) is unlikely, and collapsed from a condition close to a kinematically critical state. There are two possible scenarios to explain the curved magnetic fields of BHR 71, one is an hourglass-like field structure due to mass accumulation and the other is the Inoue \& Fukui (2013) mechanism, which proposes the interaction of the core with a shock wave to create curved magnetic fields wrapping around the core.

In this study, the detailed magnetic field structure of the dense protostellar core Barnard 335 (B335) was revealed based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. Magnetic fields pervading B335 were mapped using 24 stars after subtracting unrelated ambient polarization components, for the first time revealing that they have an axisymmetrically distorted hourglass-shaped structure toward the protostellar core. On the basis of simple two- and three-dimensional magnetic field modeling, magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $90^{\circ} \pm 7^{\circ}$ and $50^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B335 was determined to be $30.2 \pm 17.7$ $\mu {\rm G}$. The critical mass of B335, evaluated using both magnetic and thermal/turbulent support against collapse, was determined to be $M_{\rm cr} = 3.37 \pm 0.94$ ${\rm M}_{\odot}$, which is identical to the observed core mass of $M_{\rm core}=3.67$ M$_{\odot}$. We thus concluded that B335 started its contraction from a condition near equilibrium. We found a linear relationship in the polarization versus extinction diagram, up to $A_V \sim 15$ mag toward the stars with the greatest obscuration, which verified that our observations and analysis provide an accurate depiction of the core.

The detailed magnetic field structure of the starless dense core CB81 (L1774, Pipe 42) in the Pipe Nebula was determined based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. The magnetic fields pervading CB81 were mapped using 147 stars and axisymmetrically distorted hourglass-like fields were identified. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $4^{\circ} \pm 8^{\circ}$ and $20^{\circ} \pm 20^{\circ}$, respectively. The total magnetic field strength of CB81 was found to be $7.2 \pm 2.3$ $\mu{\rm G}$. Taking into account the effects of thermal/turbulent pressure and magnetic fields, the critical mass of CB81 was calculated to be $M_{\rm cr}=4.03 \pm 0.40$ M$_{\odot}$, which is close to the observed core mass of $M_{\rm core}=3.37 \pm 0.51$ M$_{\odot}$. We thus conclude that CB81 is in a condition close to the critical state. In addition, a spatial offset of $92''$ was found between the center of magnetic field geometry and the dust extinction distribution; this offset structure could not have been produced by self-gravity. The data also indicate a linear relationship between polarization and extinction up to $A_V \sim 30$ mag going toward the core center. This result confirms that near-infrared polarization can accurately trace the overall magnetic field structure of the core.

Observational data for the hourglass-like magnetic field toward the starless dense core FeSt 1-457 were compared with a flux freezing magnetic field model (Myers et al. 2018). Fitting of the observed plane-of-sky magnetic field using the flux freezing model gave a residual angle dispersion comparable with the results based on a simple three-dimensional parabolic model. The best-fit parameters for the flux freezing model were a line-of-sight magnetic inclination angle of $\gamma_{\rm mag} = 35^{\circ} \pm 15^{\circ}$ and a core center to ambient (background) density contrast of $\rho_{\rm c} / \rho_{\rm bkg} = 75$. The initial density for core formation ($\rho_0$) was estimated to be $\rho_{\rm c} / 75 = 4670$ cm$^{-3}$, which is about one order of magnitude higher than the expected density ($\sim 300$ cm$^{-3}$) for the inter-clump medium of the Pipe Nebula. FeSt 1-457 is likely to have been formed from the accumulation of relatively dense gas, and the relatively dense background column density of $A_V \simeq 5$ mag supports this scenario. The initial radius (core formation radius) $R_0$ and the initial magnetic field strength $B_0$ were obtained to be 0.15 pc ($1.64 R$) and $10.8-14.6$ $\mu$G, respectively. We found that the initial density $\rho_0$ is consistent with the mean density of the nearly critical magnetized filament with magnetic field strength $B_0$ and radius $R_0$. The relatively dense initial condition for core formation can be naturally understood if the origin of the core is the fragmentation of magnetized filaments.

The magnetic field structure, kinematical stability, and evolutionary status of the starless dense core Barnard 68 (B68) are revealed based on the near-infrared polarimetric observations of background stars, measuring the dichroically polarized light produced by aligned dust grains in the core. After subtracting unrelated ambient polarization components, the magnetic fields pervading B68 are mapped using 38 stars and axisymmetrically distorted hourglass-like magnetic fields are obtained, although the evidence for the hourglass field is not very strong. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles on the plane-of-sky and in the line-of-sight direction are determined to be $47^{\circ} \pm 5^{\circ}$ and $20^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B68 is obtained to be $26.1 \pm 8.7$ $\mu {\rm G}$. The critical mass of B68, evaluated using both magnetic and thermal/turbulent support, is $M_{\rm cr} = 2.30 \pm 0.20$ ${\rm M}_{\odot}$, which is consistent with the observed core mass of $M_{\rm core}=2.1$ M$_{\odot}$, suggesting nearly critical state. We found a relatively linear relationship between polarization and extinction up to $A_V \sim 30$ mag toward the stars with deepest obscuration. Further theoretical and observational studies are required to explain the dust alignment in cold and dense regions in the core.

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

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

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 relationship between submillimeter (submm) dust emission polarization and near-infrared (NIR) $H$-band polarization produced by dust dichroic extinction was studied for the cold starless dense core FeSt 1-457. The distribution of polarization angles ($90^{\circ}$-rotated for submm) and degrees were found to be very different between at submm and NIR wavelengths. The mean polarization angles for FeSt 1-457 at submm and NIR wavelengths are $132.1^{\circ} \pm 22.0^{\circ}$ and $2.7^{\circ} \pm 16.2^{\circ}$, respectively. The correlation between $P_H$ and $A_V$ was found to be linear from outermost regions to relatively dense line of sight of $A_V \approx 25$ mag, indicating that NIR polarization reflects overall polarization (magnetic field) structure of the core at least in this density range. The flat $P_H/A_V$ versus $A_V$ correlations were confirmed, and the polarization efficiency was found to be comparable to the observational upper limit (Jones 1989). On the other hand, as reported by Alves et al., submm polarization degrees show clear linearly decreasing trend against $A_V$ from $A_V \approx 20$ mag to the densest center ($A_V \approx 41$ mag), appearing as "polarization hole" structure. The power law index for the $P_{\rm submm}$ versus $A_V$ relationship was obtained to be $\approx -1$, indicating that the alignment for the submm sensitive dust is lost. These very different polarization distributions at submm and NIR wavelengths suggest that (1) there is different radiation environment at these wavelengths or (2) submm-sensitive dust is localized or the combination of them.

In the present study, the magnetic field scaling on density, $|B| \propto \rho^{\kappa}$, was revealed in a single starless core for the first time. The $\kappa$ index of $0.78 \pm 0.10$ was obtained toward the starless dense core FeSt 1-457 based on the analysis of the radial distribution of the polarization angle dispersion of background stars measured at the near-infrared wavelengths. The result prefers $\kappa = 2/3$ for the case of isotropic contraction, and the difference of the observed value from $\kappa = 1/2$ is 2.8 sigma. The distribution of the ratio of mass to magnetic flux was evaluated. FeSt 1-457 was found to be magnetically supercritical near the center ($\lambda \approx 2$), whereas nearly critical or slightly subcritical at the core boundary ($\lambda \approx 0.98$). Ambipolar-diffusion-regulated star formation models for the case of moderate magnetic field strength may explain the physical status of FeSt 1-457. The mass-to-flux ratio distribution for typical dense cores (critical Bonnor--Ebert sphere with central $\lambda=2$ and $\kappa=1/2$--$2/3$) was calculated and found to be magnetically critical/subcritical at the core edge, which indicates that typical dense cores are embedded in and evolve from magnetically critical/subcritical diffuse surrounding medium.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

The relationship between dust polarization and extinction was determined for the cold dense starless molecular cloud core FeSt 1-457 based on the background star polarimetry of dichroic extinction at near-infrared wavelengths. Owing to the known (three-dimensional) magnetic field structure, the observed polarizations from the core were corrected by considering (a) the subtraction of the ambient polarization component, (b) the depolarization effect of inclined distorted magnetic fields, and (c) the magnetic inclination angle of the core. After these corrections, a linear relationship between polarization and extinction was obtained for the core in the range up to $A_V \approx 20$ mag. The initial polarization vs. extinction diagram changed dramatically after the corrections of (a) to (c), with the correlation coefficient being refined from 0.71 to 0.79. These corrections should affect the theoretical interpretation of the observational data. The slope of the finally obtained polarization--extinction relationship is $P_H / E_{H-K_s} = 11.00 \pm 0.72$ $\%$ ${\rm mag}^{-1}$, which is close to the statistically estimated upper limit of the interstellar polarization efficiency (Jones 1989). This consistency suggests that the upper limit of interstellar polarization efficiency might be determined by the observational viewing angle toward polarized astronomical objects.

In this study, we investigate the line emissions from cold molecular gas based on our previous "radiation-driven fountain model" (Wada 2016), which reliably explains the spectral energy distribution of the nearest type 2 Seyfert galaxy, the Circinus galaxy. Using a snapshot of the best-fit radiation-hydrodynamic model for the central r < 16 pc, in which non-equilibrium X-ray-dominated region chemistry is solved, we conduct post-processed, non-local thermodynamic equilibrium radiation transfer simulations for the CO lines. We obtain a spectral line energy distribution with a peak around J=6, and its distribution suggests that the lines are not thermalized. However, for a given line-of-sight, the optical depth distribution is highly non-uniform between $\tau \ll 1$ and $\tau \gg 1$. The CO-to-H2 conversion factor (X_CO), which can be directly obtained from the results, is not a constant and depends strongly on the integrated intensity, and it differs from the fiducial value for local objects. X_CO exhibits a large dispersion of more than one order of magnitude, reflecting the non-uniform internal structure of a "torus." We also found that the physical conditions differ between grid cells on a scale of a few parsecs along the observed lines of sight; therefore, a specific observed line ratio does not necessarily represent a single physical state of the ISM.

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.

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.

Three dimensional (3D) magnetic field information on molecular clouds and cores is important for revealing their kinematical stability (magnetic support) against gravity which is fundamental for studying the initial conditions of star formation. In the present study, the 3D magnetic field structure of the dense starless core FeSt 1-457 is determined based on the near-infrared polarimetric observations of the dichroic polarization of background stars and simple 3D modeling. With an obtained angle of line-of-sight magnetic inclination axis $\theta_{\rm inc}$ of $45^{\circ}\pm10^{\circ}$ and previously determined plane-of-sky magnetic field strength $B_{\rm pol}$ of $23.8\pm12.1$ $\mu{\rm G}$, the total magnetic field strength for FeSt 1-457 is derived to be $33.7\pm18.0$ $\mu{\rm G}$. The critical mass of FeSt 1-457, evaluated using both magnetic and thermal/turbulent support is ${M}_{\rm cr} = 3.70\pm0.92$ ${\rm M}_{\odot}$, which is identical to the observed core mass, $M_{\rm core}=3.55\pm0.75$ ${\rm M}_{\odot}$. We thus conclude that the stability of FeSt 1-457 is in a condition close to the critical state. Without infalling gas motion and no associated young stars, the core is regarded to be in the earliest stage of star formation, i.e., the stage just before the onset of dynamical collapse following the attainment of a supercritical condition. These properties would make FeSt 1-457 one of the best starless cores for future studies of the initial conditions of star formation.

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

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 mechanisms causing millimeter-wave polarization in protoplanetary disks are under debate. To disentangle the polarization mechanisms, we observe the protoplanetary disk around HL Tau at 3.1 mm with the Atacama Large Millimeter/submillimeter Array (ALMA), which had polarization detected with CARMA at 1.3 mm. We successfully detect the ring-like azimuthal polarized emission at 3.1 mm. This indicates that dust grains are aligned with the major axis being in the azimuthal direction, which is consistent with the theory of radiative alignment of elongated dust grains, where the major axis of dust grains is perpendicular to the radiation flux. Furthermore, the morphology of the polarization vectors at 3.1 mm is completely different from those at 1.3 mm. We interpret that the polarization at 3.1 mm to be dominated by the grain alignment with the radiative flux producing azimuthal polarization vectors, while the self-scattering dominates at 1.3 mm and produces the polarization vectors parallel to the minor axis of the disk. By modeling the total polarization fraction with a single grain population model, the maximum grain size is constrained to be $100{\rm~\mu m}$, which is smaller than the previous predictions based on the spectral index between ALMA at 3 mm and VLA at 7 mm.

Magnetic fields are believed to play an important role in controlling the stability and contraction of dense condensations of gas and dust leading to the formation of stars and planetary systems. In the present study, the magnetic field of FeSt 1-457, a cold starless molecular cloud core, was mapped on the basis of the polarized near-infrared light from 185 background stars after being dichroically absorbed by dust aligned with the magnetic field in the core. A distinct "hourglass-shaped" magnetic field was identified in the region of the core, which was interpreted as the first evidence of a magnetic field structure distorted by mass condensation in a starless core. The steep curvature of the magnetic field lines obtained in the present study indicates that the distortion was mainly created during the formation phase of the dense core. The derived mass-to-magnetic flux ratio indicates that the core is in a magnetically supercritical state. However, the stability of the core can be considered to be in a nearly critical state if the additional contributions from the thermal and turbulent support are included. Further diffusion of the magnetic field and/or turbulent dissipation would cause the onset of dynamical collapse of the core. The geometrical relationship between the direction of the magnetic field lines and the elongation of the core was found to be in good agreement with the theoretical predictions for the formation of Sun-like stars under the influence of a magnetic field.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

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.

Filaments are considered to be basic structures and molecular clouds consist of filaments. Filaments are often observed as extending in the direction perpendicular to the interstellar magnetic field. The structure of filaments has been studied based on a magnetohydrostatic equilibrium model (Tomisaka 2014). Here, we simulate the expected polarization pattern for isothermal magnetohydrostatic filaments. The filament exhibits a polarization pattern in which the magnetic field is apparently perpendicular to the filament when observed from the direction perpendicular to the magnetic field. When the line-of-sight is parallel to the global magnetic field, the observed polarization pattern is dependent on the center-to-surface density ratio for the filament and the concentration of the gas mass toward the central magnetic flux tube. Filaments with low center-to-surface density ratios have an insignificant degree of polarization when observed from the direction parallel to the global magnetic field. However, models with a large center-to-surface density ratio have polarization patterns that indicate the filament is perpendicularly threaded by the magnetic field. When mass is heavily concentrated at the central magnetic flux tube, which can be realized by the ambipolar diffusion process, the polarization pattern is similar to that expected for a low center-to-surface density contrast.

We have constructed two types of analytical models for an isothermal filamentary cloud supported mainly by magnetic tension. The first one describes an isolated cloud while the second considers filamentary clouds spaced periodically. Both the models assume that the filamentary clouds are highly flattened. The former is proved to be the asymptotic limit of the latter in which each filamentary cloud is much thinner than the distance to the neighboring filaments. We show that these models reproduce main features of the 2D equilibrium model of Tomisaka (2014) for filamentary cloud threaded by perpendicular magnetic field. It is also shown that the critical mass to flux ratio is $ M /\Phi = (2 \pi \sqrt{G}) ^{-1} $, where $ M $, $ \Phi $ and $ G $ denote the cloud mass, the total magnetic flux of the cloud, and the gravitational constant, respectively. This upper bound coincides with that for an axisymmetric cloud supported by poloidal magnetic fields. We applied the variational principle for studying the Jeans instability of the first model. Our model cloud is unstable against fragmentation as well as the filamentary clouds threaded by longitudinal magnetic field. The fastest growing mode has a wavelength several times longer than the cloud diameter. The second model describes quasi-static evolution of filamentary molecular cloud by ambipolar diffusion.

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.

Herschel observation has recently revealed that interstellar molecular clouds consist of many filaments. Polarization observations in optical and infrared wavelengths indicate that the magnetic field often runs perpendicular to the filament. In this paper, the magnetohydrostatic configuration of isothermal gas is studied, in which the thermal pressure and the Lorentz force are balanced against the self-gravity and the magnetic field is globally perpendicular to the axis of the filament. The model is controlled by three parameters: center-to-surface density ratio (\rho_c/\rho_s), plasma \beta of surrounding interstellar gas (\beta_0) and the radius of the hypothetical parent cloud normalized by the scale-height (R'_0), although there remains a freedom how the mass is distributed against the magnetic flux (mass loading). In the case that R'_0 is small enough, the magnetic field plays a role in confining the gas. However, the magnetic field generally has an effect in supporting the cloud. There is a maximum line-mass (mass per unit length) above which the cloud is not supported against the gravity. Compared with the maximum line-mass of non-magnetized cloud (2c_s^2/G, where c_s and G represent respectively the isothermal sound speed and the gravitational constant), that of the magnetized filament is larger than the non-magnetized one. The maximum line-mass is numerically obtained as \lambda_max \simeq 0.24 \Phi_cl/G^1/2 + 1.66 c_s^2/G, where \Phi_cl represents one half of the magnetic flux threading the filament per unit length. The maximum mass of the filamentary cloud is shown to be significantly affected by the magnetic field when the magnetic flux per unit length exceeds \Phi_cl \gtrsim 3 pc \mu G (c_s/190 m s^-1)^2.

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 molecular evolution from a molecular cloud core to a first hydrostatic core in three spatial dimensions. We perform a radiation hydrodynamic simulation in order to trace fluid parcels, in which molecular evolution is investigated, using a gas-phase and grain-surface chemical reaction network. We derive spatial distributions of molecular abundances and column densities in the core harboring the first core. We find that the total of gas and ice abundances of many species in a cold era (10 K) remain unaltered until the temperature reaches ~500 K. The gas abundances in the warm envelope and the outer layer of the first core (T < 500 K) are mainly determined via the sublimation of ice-mantle species. Above 500 K, the abundant molecules, such as H2CO, start to be destroyed, and simple molecules, such as CO, H2O and N2 are reformed. On the other hand, some molecules are effectively formed at high temperature; carbon-chains, such as C2H2 and cyanopolyynes, are formed at the temperature of >700 K. We also find that large organic molecules, such as CH3OH and HCOOCH3, are associated with the first core (r < 10 AU). Although the abundances of these molecules in the first core stage are comparable or less than in the protostellar stage (hot corino), reflecting the lower luminosity of the central object, their column densities in our model are comparable to the observed values toward the prototypical hot corino, IRAS 16293-2422. We propose that these large organic molecules can be good tracers of the first cores.

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.

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

A first core is a first hydrostatic object formed in the course of dynamical contraction of a molecular cloud core. Since the inflow pattern changes drastically before and after the first core formation, it is regarded as a milestone in the star formation process. In order to identify the first core from a mapping observation, the features expected for the first core are studied for CS rotation transitions at radio wavelengths. The non-LTE radiation transfer is calculated for the results of radiation magnetohydrodynamical simulations of the contraction of the magnetized molecular cloud core in rotation (Tomida et al. 2010a). We use the Monte-Carlo method to solve the non-LTE radiation transfer in a nested grid hierarchy. In the first core phase, an outflow arises from the vicinity of the first core due to the twisted magnetic field amplified by the rotation motion of the contracting gas disk. The disk and outflow system has several characteristic observational features: (i) in the pole-on view, relatively opaque lines indicate asymmetry in the emission lines in which the blue side is stronger than the red side (an infall signature of the envelope); (ii) in the edge-on view, the disk has a signature of simultaneous rotation and infall, i.e., the integrated intensity of the reaching side is brighter than that of the receding side and the gradient in the intensity-weighted velocity is larger in the reaching side; (iii) the observed outflow indicates rotation around the rotation axis. The size of the outflow gives the approximate age after the first core is formed, since the outflow is not expected for the earlier runaway isothermal collapse phase.

The observational expectation of polarization measurements of thermal dust radiation is investigated to find information on molecular outflows based on magnetohydrodynamical (MHD) and radiation transfer simulations. There are two major proposed models for the driving of molecular outflows: (1) molecular gas is accelerated by a magnetic pressure gradient or magnetocentrifugal wind mechanism before the magnetic field and molecular gas are decoupled, (2) the linear momentum of a highly collimated jet is transferred to the ambient molecular gas. In order to distinguish between these two models, it is crucial to observe the configuration of the magnetic field. An observation of a toroidal magnetic field is strong evidence that the first of the models is appropriate. In this paper, we calculated the polarization distribution of thermal dust radiation due to the alignment of dust grains along the magnetic field using molecular outflow data calculated by two-dimensional axisymmetric MHD simulations. An asymmetric distribution around the z-axis is characteristic for magnetic fields composed of both poloidal and toroidal components. We determined that the outflow has a low polarization degree compared with the envelope and that the envelope and outflow have different polarization directions (B-vector), namely, the magnetic field within the envelope is parallel to the global magnetic field lines while the magnetic field of the outflow is perpendicular to it. Thus we have demonstrated that the point-symmetric (rather than axisymmetric) distributions of low polarization regions indicate that molecular outflows are likely to be magnetically driven. Observations of this polarization distribution with tools such as ALMA would confirm the origin of the molecular outflow.

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.

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

Origin of the magnetic field ubiquitous in the Universe is studied based on the Biermann mechanism, which is expected to work in the non-barotropic region. We perform a series of two-dimensional MHD simulations of the first generation supernova remnant (SNR) expanding in the inhomogeneous interstellar matter (ISM) and study the Biermann mechanism working in the interior of the SNR. Especially, we pay attention to the relaxation process of electron and ion temperatures via the Coulomb interaction. In the early SNR in which the electron temperature is much lower than the ion temperature, the Biermann mechanism is ineffective, since the gradient of electron pressure is small. Magnetic fields begin to be generated just behind the shock front when the electron temperature is sufficiently relaxed. Assuming the explosion energy of 10^52 erg, the total magnetic energy generated reaches about 10^26 erg and does not depend strongly on the parameters of either SNR or ISM. Analytic expression to estimate the magnetic total energy is presented and it is shown this agrees well with our numerical results. Finally we evaluate the expected amplitude of magnetic fields in protogalaxies as ~10^-19 G, which is sufficient for seed fields of the observed galactic magnetic fields.

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

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.

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 investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

Recent observations suggest molecular line ratios in millimeter and submillimeter bands may be a good tool to reveal the long-standing question on the origin of energy sources in obscured active galaxies -- AGN and/or starburst. Observations of actual molecular medium show in general inhomogeneous structures as well as high-resolution hydrodynamic simulations do. In order for precise interpretation of emergent line emission from the inhomogeneous molecular gas to probe the dominant energy source of active galaxies, we study characteristic features of emergent intensities via three-dimensional non-LTE (non-local thermodynamic equilibrium) line transfer simulations. Our results succeeded in making clear 1) the necessary conditions for HCN/HCO$^{+}$-dichotomy, and 2) importance of clumpiness on intensity ratio and its interpretation. These results are obtained for the first time by our realistic three-dimensional simulations, and line transfer simulations will be a powerful tool to comprehensive studies of extragalactic interstellar medium (ISM) in forthcoming ALMA (Atacama Large Millimeter/submillimeter Array) era.

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

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

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.

We explore millimeter line diagnostics of an obscuring molecular torus modeled by a hydrodynamic simulation with three-dimensional nonLTE radiative transfer calculations. Based on the results of high-resolution hydrodynamic simulation of the molecular torus around an AGN, we calculate intensities of HCN and HCO^{+} rotational lines as two representative high density tracers. The three-dimensional radiative transfer calculations shed light on a complicated excitation state in the inhomogeneous torus, even though a spatially uniform chemical structure is assumed. Our results suggest that HCN must be much more abundant than HCO^{+} in order to obtain a high ratio ($R_{HCN/HCO+}\sim 2$) observed in some of the nearby galaxies. There is a remarkable dispersion in the relation between integrated intensity and column density, indicative of possible shortcomings of HCN(1-0) and HCO^{+}(1-0) lines as high density tracers. The internal structures of the inhomogeneous molecular torus down to subparsec scale in external galaxies will be revealed by the forthcoming Atacama Large Millimeter/submillimeter Array (ALMA). The three-dimensional radiative transfer calculations of molecular lines with high-resolution hydrodynamic simulation prove to be a powerful tool to provide a physical basis for molecular line diagnostics of the central regions of external galaxies.

Developing a theory of low-mass star formation ($\sim 0.1$ to 3~M$_{\odot}$) remains one of the most elusive and important goals of theoretical astrophysics. The star-formation process is the outcome of the complex dynamics of interstellar gas involving non-linear interactions of turbulence, gravity, magnetic field and radiation. The evolution of protostellar condensations, from the moment they are assembled by turbulent flows to the time they reach stellar densities, spans an enormous range of scales, resulting in a major computational challenge for simulations. Since the previous Protostars and Planets conference, dramatic advances in the development of new numerical algorithmic techniques have been successfully implemented on large scale parallel supercomputers. Among such techniques, Adaptive Mesh Refinement and Smooth Particle Hydrodynamics have provided frameworks to simulate the process of low-mass star formation with a very large dynamic range. It is now feasible to explore the turbulent fragmentation of molecular clouds and the gravitational collapse of cores into stars self-consistently within the same calculation. The increased sophistication of these powerful methods comes with substantial caveats associated with the use of the techniques and the interpretation of the numerical results. In this review, we examine what has been accomplished in the field and present a critique of both numerical methods and scientific results. We stress that computational simulations should obey the available observational constraints and demonstrate numerical convergence. Failing this, results of numerical simulations do not advance our understanding of low-mass star formation.

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 estimate the polarized thermal dust emission from MHD simulations of protostellar collapse and outflow formation in order to investigate alignment of outflows with magnetic fields. The polarization maps indicate that alignment of an outflow with the magnetic field depends on the field strength inside the cloud core; the direction of the outflow, projected on the plane of the sky, is aligned preferentially with the mean polarization vector for a cloud core with a magnetic field strength of 80 microgauss, while it does not tend to be aligned for 50 microgauss as long as the 1000 AU scale is considered. The direction of the magnetic field at the cloud center is probed by the direction of the outflow. In addition, the magnetic field at the cloud center can be revealed by ALMA even when the source is embedded deeply in the envelope. The Chandrasekhar-Fermi formula is examined using the polarization maps, indicating that the field strength predicted by the formula should be corrected by a factor of 0.24 - 0.44. The correction factor has a tendency to be lower for a cloud core with a weaker magnetic field.

We present the evolution path of the first core in rotating molecular cloud core and the evolution during the second collapse phase using axisymmetric numerical calculations. The structure and evolution of the rotating first core is characterized by the angular momentum $J_{\rm core}$ and mass $M_{\rm core}$, both of which increase with time by accretion from the infalling envelope. We find the evolution path of it can be considered as a sequence of equilibrium solutions with a constant $J_{\rm core}/M_{\rm core}^2$. Such evolution paths are in good agreement with the results of three-dimensional numerical simulations. The evolution paths is completely different from spherical case. Except for extremely small rotation rate of $J_{\rm core}/M_{\rm core}^2 < 0.015 G/(\sqrt{2} c_{\rm iso})$, the first core does not increase in the central density beyond a certain maximum value and can not begins the second collapse direct and grows in mass by accretion. Such a first core has a very long time-scale and seems to suffer from a non-axisymmetric evolution. In the usual case, we expect that the time-scale of a rotating first core becomes several thousand years and several percent of the prestellar cores contain a first core at least. Such a rotating first core has a small average luminosity of $L_{\rm core} = 0.003-0.03$($\dot{M}_{\rm core}/\dot{M}_{\rm ref}$) L$_{\odot}$.

The evolution of supernova remnants (SNRs) is studied, with particular attention to the effect of magnetic fields with axisymmetric two-dimensional magnetohydrodynamical simulations. The evolution of magnetic SNRs is the same as non-magnetic ones in the adiabatic Sedov stage. After a thin shell is formed, the shell is driven by the pressure of the hot interior gas (bubble). Evolution in the pressure-driven snow-plow phase is much affected by the magnetic field. The shell sweeping the magnetic field lines thickens owing to the magnetic pressure force. After $5\times 10^5{\rm yr}$ - $2\times 10^6{\rm yr}$, the inner boundary of the thick shell begins to contract. This compresses the hot bubble radially and maintains its thermal pressure. Thus, the bubble forms a prolate spheroidal shape and becomes thinner and thinner, since it expands in a direction parallel to the magnetic field for $B_0 \gtrsim 3\mu {\rm G}$. Finally, the bubble contracts. The porosity of the hot low-density gas in ISM is reduced, taking the effect of the magnetic field into account.

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.

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.

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 have performed three-dimensional, non-LTE radiative transfer calculations for 12CO and 13CO lines, applying them to our high-resolution hydrodynamic models of the `torus' around a SMBH in an active galactic nucleus. The hydrodynamic simulations reveal inhomogeneous and turbulent gas structure on a sub-pc scale in a circum-nuclear starburst region. Thick disks interlaced with filaments, clumps and holes are naturally formed due to the interplay among energy feedback from supernovae, self-gravity of the gas, galactic rotation, and radiative cooling. The intensity maps of the molecular lines for the circum-nuclear disks show a clumpy structure reflecting the intrinsic inhomogeneity and turbulent motion of the gas disk. The fine structure of the `torus' could be resolved in the nearby active galaxies using ALMA. We also found that the X-factor is not uniformly distributed in the central 100 pc region. The X-factor derived for 12CO(J=1-0) intensity depends strongly on the intensity, whereas the X for 12CO(3-2) is nearly constant over two orders of magnitude of the intensity. The suggested conversion factor for the molecular gas mass is X_CO (3-2) ~ 0.27 x 10^{20} cm^{-2} (K km s^{-1})^{-1}.

The collapse of slowly rotating molecular cloud cores threaded by magnetic fields is investigated by high-resolution numerical simulation. Outflow formation in the collapsing cloud cores is also followed. In the models examined, the cloud core and parent cloud rotate rigidly and are initially threaded by a uniform magnetic field. The simulations show that the cloud core collapses along the magnetic field lines. The magnetic field in the dense region of the cloud core rotates faster than that of the parent cloud as a consequence of spin-up of the central region during the collapse. The cloud core exhibits significant precession of the rotation axis, magnetic field, and disk orientation, with precession highest in the models with low initial field strength ($\lesssim 20 \mu {\rm G}$). Precession in models with initial fields of $\sim 40 \mu {\rm G}$ is suppressed by strong magnetic braking. Magnetic braking transfers angular momentum form the central region and acts more strongly on the component of angular momentum oriented perpendicular to the magnetic field. After the formation of an adiabatic core, outflow is ejected along the local magnetic field lines. Strong magnetic braking associated with the outflow causes the direction of angular momentum to converge with that of the local magnetic field, resulting in the convergence of the local magnetic field, angular momentum, outflow, and disk orientation by the outflow formation phase. The magnetic field of a young star is inclined at an angle of no more than $30^\circ$ from that of the parent cloud at initial field strengths of $\sim 20 \mu {\rm G}$, while at an initial field strength of $\sim 40 \mu {\rm G}$, the magnetic field of the young star is well aligned with that of the parent cloud.

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.

Collapse of the rotating magnetized molecular cloud core is studied with the axisymmetric magnetohydrodynamical (MHD) simulations. Due to the change of the equation of state of the interstellar gas, the molecular cloud cores experience several different phases as collapse proce eds. In the isothermal run-away collapse ($n \lesssim 10^{10}{\rm H_2 cm}^{-3}$), a pseudo-disk is formed and it continues to contract till the opaque core is fo rmed at the center. In this disk, a number of MHD fast and slow shock pairs appear running parallelly to the disk. After the equation of state becomes hard, an adiabatic core is formed, which is separated from the isothermal contracting pseudo-disk by the accretion shock front facing radially outwards. By the effect of the magnetic tension, the angular momentum is transferred from the disk mid-plane to the surface. The gas with excess angular momentum near the surface is finally ejected, which explains the molecular bipolar outflow. Two types of outflows are observed. When the poloidal magnetic field is strong (magnetic energy is comparable to the thermal one), a U-shaped outflow is formed in which fast moving gas is confined to the wall whose shape looks like a capit al letter U. The other is the turbulent outflow in which magnetic field lines and velocity fi elds are randomly oriented. In this case, turbulent gas moves out almost perpendicularly from the disk. The continuous mass accretion leads to the quasistatic contraction of the first core. A second collapse due to dissociation of H$_2$ in the first core follows. Finally another quasistatic core is again formed by atomic hydrogen (the second core). It is found that another outflow is ejected around the second atomic core, which seems to correspond to the optical jets or the fast neutral winds.

Dynamical contraction of a slowly-rotating magnetized cloud is studied using the magnetohydrodynamical (MHD) simulations. In the isothermal stage ($n \la n_{\rm A} \sim 10^{10}{\rm cm}^{-3}$), the cloud evolves similarly to that expected from Larson-Penston self-similar solution and experiences the run-away collapse. However, after the central density exceeds $\sim n_{\rm A}$, an accretion disk is formed around the adiabatic core. Just outside the core, outflow is ejected by the effect of magnetic torque (magnetocentrifugal wind). Since $\sim$ 10% of the mass is ejected with almost all the angular momentum, the specific angular momentum of the protostellar core reduces to that observed in main-sequence stars.

If the angular momentum of the molecular cloud core were conserved during the star formation process, a new-born star would rotate much faster than its fission speed. This constitutes the angular momentum problem of new-born stars. In this paper, the angular momentum transfer in the contraction of a rotating magnetized cloud is studied with axisymmetric MHD simulations. Owing to the large dynamic range covered by the nested-grid method, the structure of the cloud in the range from 10 AU to 0.1 pc is explored. First, the cloud experiences a run-away collapse, and a disk forms perpendicularly to the magnetic field, in which the central density increases greatly in a finite time-scale. In this phase, the specific angular momentum j of the disk decreases to $\simeq 1/3$ of the initial cloud. After the central density of the disk exceeds $\sim 10^{10}{\rm cm}^{-3}$, the infall on to the central object develops. In this accretion stage, the rotation motion and thus the toroidal magnetic field drive the outflow. The angular momentum of the central object is transferred efficiently by the outflow as well as the effect of the magnetic stress. In 7000 yr from the core formation, the specific angular momentum of the central $0.17M_\odot$ decreases a factor of 10^{-4} from the initial value (i.e. from $10^{20}{\rm cm^2 s^{-1}}$ to $10^{16}{\rm cm^2 s^{-1}}$).

Dynamical collapses of magnetized molecular cloud cores are studied with magnetohydrodynamical simulations from the run-away collapse phase to the accretion phase. In the run-away collapse phase, a disk threaded by magnetic field lines is contracting due to its self-gravity and its evolution is well expressed by a self-similar solution. The central density increases greatly in a finite time scale and reaches a density at which an opaque core is formed at the center. After that, matter accretes to the newly formed core (accretion phase). In this stage, a rotationally supported disk is formed in a cloud core without magnetic fields. In contrast, the disk continues to contract in the magnetized cloud core, since the magnetic fields transfer angular momentum from the disk. Its rotation motion winds up the threading magnetic field lines. Eventually, strong toroidal magnetic fields are formed and begin to drive the outflow, even if there is no toroidal field component initially. Bipolar molecular outflows observed in protostar candidates are naturally explained by this model.

The importance of the interstellar magnetic field is studied in relation to the evolutions of superbubbles with a three-dimensional (3D) numerical magnetohydrodynamical (MHD) simulation. A superbubble is a large supernova remnant driven by sequential supernova explosions in an OB association. Its evolution is affected by the density stratification in the galactic disk. After the size reaches 2--3 times the density scale-height, the superbubble expands preferentially in the $z$-direction. Finally it can punch out the gas disk (blow-out). On the other hand, the magnetic field running parallel to the galactic disk has an effect to prevent from expanding in the direction perpendicular to the field. The density stratification and the magnetic fields have completely opposite effects on the evolution of the superbubble. We present results of 3D MHD simulation in which both effects are included. As a result, it is concluded that when the magnetic field has a much larger scale-height than the density, even for a model that the bubble would blow out from the disk if the magnetic field were absent, the magnetic field with the strength of 5 $\mu$G can confine the bubble in $|z| \alt 300$ pc for $\simeq$20 Myr (confinement). In a model that the field strength decreases in the halo in proportion to $B \propto \rho^{1/2}$, the superbubble eventually blows out like a model of B=0 even if the magnetic field in the mid-plane is as strong as $B=5\mu$G.

We present ASCA and ROSAT observations of SN 1986J covering the period 1991 August to 1996 January. From observations with the ROSAT HRI and PSPC, we find that the 0.5-2.5 keV flux decreased proportional to $t^{-2}$ during this period; the ASCA data are consistent with this result and extend it to the 2-10 keV band. ASCA spectra from 1994 January and 1996 January are consistent with thermal emission from a solar metallicity plasma at an equilibrium temperature kT = 5-7.5 keV, somewhat hotter than that observed from other X-ray supernovae. These spectra also show a clear Fe K emission line at 6.7 keV with FWHM < 20,000 km/s (90% confidence). This limit on the line width is consistent with the reverse shock model of Chevalier & Fransson (1994), but does not rule out the clumpy wind model of Chugai (1993).

Accretion flow in a contracting magnetized isothermal cloud was studied using magnetohydrodynamical simulations and a nested grid technique. First, the interstellar magnetized cloud experiences a ``runaway collapse'' phase, in which the central density increases drastically within a finite time scale. Finally, it enters an accretion phase, in which inflowing matter accretes onto a central high-density disk or a new-born star. We found that the accretion rate reaches (4 -- 40) $\times c_s^3/G$, where $c_s$ and $G$ represent the isothermal sound speed and the gravitational constant, respectively. This is much larger than the standard accretion rate of $0.975c_s^3/G$ for a hydrostatic isothermal spherical cloud (Shu 1977, AAA19.065.044). Due to the effect of an extra infall velocity achieved in the runaway phase ($\sim 2 c_s$), the accretion rate is boosted. This rate declines with time in contrast to Shu's solution, but keeps $\gtsim 2.5 c_s^3/G$. The observed gas infall rate around proto-stars such as L1551 IRS 5 and HL Tau is also discussed.

Fragmentation process in a cylindrical magnetized cloud is studied with the nested grid method. The nested grid scheme use 15 levels of grids with different spatial resolution overlaid subsequently, which enables us to trace the evolution from the molecular cloud density $\sim 100 {\rm cm}^{-3}$ to that of the protostellar disk $\sim 10^{10} {\rm cm} ^{-3}$ or more. Fluctuation with small amplitude grows by the gravita- tional instability. It forms a disk perpendicular to the magnetic fields which runs in the direction parallel to the major axis of the cloud. Matter accrets on to the disk mainly flowing along the magnetic fields and this makes the column density increase. The radial inflow, whose velocity is slower than that perpendicular to the disk, is driven by the increase of the gravity. While the equation of state is isothermal and magnetic fields are perfectly coupled with the matter, which is realized in the density range of $\rho \lesssim 10^{10} {\rm cm}^{-3}$, never stops the contraction. The structure of the contracting disk reaches that of a singular solution as the density and the column density obey $\rho(r)\propto r^{-2}$ and $\sigma(r) \propto r^{-1}$, respectively. The magnetic field strength on the mid-plane is proportional to $\rho(r)^{1/2}$ and further that of the center ($B_c$) evolves as proportional to the square root of the gas density ($\propto \rho_c^{1/2}$). It is shown that isothermal clouds experience ``run-away'' collapses. The evolution after the equation of state becomes hard is also discussed.

Gravitational collapse of the cylindrical elongated cloud is studied by numerical magnetohydrodynamical simulations. In the infinitely long cloud in hydrostatic configuration, small perturbations grow by the gravitational instability. The most unstable mode indicated by a linear perturbation theory grows selectively even from a white noise. The growth rate agrees with that calculated by the linear theory. First, the density-enhanced region has an elongated shape, i.e., prolate spheroidal shape. As the collapse proceeds, the high-density fragment begins to contract mainly along the symmetry axis. Finally, a spherical core is formed in the non-magnetized cloud. In contrast, an oblate spheroidal dense disk is formed in a cloud in which the magnetic pressure is nearly equal to the thermal one. The radial size of the disk becomes proportional to the initial characteristic density scale-height in the r-direction. As the collapse proceeds, a slowly contracting dense part is formed (approximately < 10% in mass) inside of the fast contracting disk. And this is separated from other part of the disk whose inflow velocity is accelerated as reaching the center of the core. From arguments on the Jeans mass and the magnetic critical mass, it is concluded that the fragments formed in a cylindrical elongated cloud can not be supported against the self- gravity and it will eventually collapse.

Reanalysis of Einstein IPC data and new observations from the GINGA LAC indicate the presence of extended X-ray emission (10-50 kpc) around the starburst galaxy M82. Here we model this emission by calculating numerical hydrodynamic simulations of the starburst event to much later times and larger scales than previously considered. For our models, we adopt a supernova rate of 0.1 ${\rm yr}^{-1}$, and an extended low-density static halo that is bound to the galaxy. There are three stages to the evolution of the wind-blown bubble and the propagation of the shock front: the bubble expands in an almost uniform density disk gas, with a deceleration of the shock front ($t \alt $ 3.6 Myr); breakout from the disk and the upward acceleration of the shock front (3.6 Myr $\alt t \alt$ 18 Myr); propagation into the halo, leading to a more spherical system and shock deceleration (18 Myr $\alt t$). For a halo density of $10^{-3} {\rm cm}^{-3}$, the outflow reaches a distance of 40-50 kpc from the center of the starburst galaxy in 50 Myr. We calculate the time evolution of the X-ray luminosity and find that the extended starburst emits $3\times 10^{39}\lcgs$ to $10^{40}\lcgs$ in the GINGA LAC band and $\sim 10^{41}\lcgs$ in the Einstein or ROSAT HRI band. The degree of the ionization equilibrium in the outflow and its effect on the iron K$\alpha$ line emission are discussed.