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 report ALMA Cycle 4 observations of the Class I binary protostellar system L1551 IRS 5 in the 0.9-mm continuum emission, C18O (J=3-2), OCS (J=28-27), and four other Band 7 lines. At ~0.07" (= 10 au) resolution in the 0.9 mm emission, two circumstellar disks (CSDs) associated with the binary protostars are separated from the circumbinary disk (CBD). The CBD is resolved into two spiral arms, one connecting to the CSD around the northern binary source, Source N, and the other to Source S. As compared to the CBD in the neighboring protobinary system L1551 NE, the CBD in L1551 IRS 5 is more compact (r ~150 au) and the m=1 mode of the spirals found in L1551 NE is less obvious in L1551 IRS 5. Furthermore, the dust and molecular-line brightness temperatures of CSDs and CBD reach >260 K and >100 K, respectively, in L1551 IRS 5, much hotter than those in L1551 NE. The gas motions in the spiral arms are characterized by rotation and expansion. Furthermore, the transitions from the CBD to the CSD rotations at around the L2 and L3 Lagrangian points and gas motions around the L1 point are identified. Our numerical simulations reproduce the observed two spiral arms and expanding gas motion as a result of gravitational torques from the binary, transitions from the CBD to the CSD rotations, and the gas motion around the L1 point. The higher temperature in L1551 IRS 5 likely reflects the inferred FU-Ori event.

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.

We present our observational results of the 0.87 mm polarized dust emission in the Class 0 protostar B335 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA) at a 0.2" (20 au) resolution. We compared our data at 0.87 mm with those at 1.3 mm from the ALMA archive. The observed polarization orientations at the two wavelengths are consistent within the uncertainty, and the polarization percentages are systematically higher at 1.3 mm than 0.87 mm by a factor of ~1.7, suggesting that the polarized emission originates from magnetically aligned dust grains. We inferred the magnetic field orientations from the observed polarization orientations. We found that the magnetic field changes from ordered and highly pinched to more complicated and asymmetric structures within the inner 100 au scale of B335, and the magnetic field connects to the center along the equatorial plane as well as along the directions which are ~40-60 degrees from the equatorial plane. We performed non-ideal MHD simulations of collapsing dense cores. We found that similar magnetic field structures appear in our simulations of dense cores with the magnetic field and rotational axis slightly misaligned by 15 degrees but not in those with the aligned magnetic field and rotational axis. Our results suggest that the midplane of the inner envelope within the inner 100 au scale of B335 could be warped because of the misaligned magnetic field and rotational axis, and the magnetic field could be dragged by the warped accretion flows.

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

We present CS(7-6) line maps toward the central parsec of the Galactic Center (GC), conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The primary goal is to find and characterize the gas structure in the inner cavity of the circumnuclear disk (CND) in high resolution (1.3"=0.05 pc). Our large field-of-view mosaic maps -- combining interferometric and single-dish data that recover extended emission - provide a first homogeneous look to resolve and link the molecular streamers in the CND with the neutral nuclear filaments newly detected within the central cavity of the CND. We find that the nuclear filaments are rotating with Keplerian velocities in a nearly face-on orbit with an inclination angle of ~10-20 degree (radius <= 0.5 pc). This is in contrast to the CND which is highly inclined at ~65-80 degree (radius ~2-5 pc). Our analysis suggests a highly warped structure from the CND to the nuclear filaments. This result may hint that the nuclear filaments and the CND were created by different external clouds passing by Sgr A*.

We investigate the silicate feature of the two Class I components of HH250-IRS, a resolved binary system with a separation of $0''53$ driving a Herbig-Haro flow. Each component has its own circumstellar envelope, and the system is surrounded by a circumbinary disk. We have carried out low resolution spectroscopy in the 8-13$\mu$m range using VISIR at ESO's Very Large Telescope. The silicate features of both sources are clearly different. The NW component has a broad, smooth absorption profile lacking structure. The SE component shows the silicate feature in emission, with structure longwards of 9.5$\mu$m indicating the presence of crystalline dust in the dominant form of forsterite. The apparent lack of an absorption feature caused by foreground dust is probably due to the filling of the band with emission by amorphous silicates in the envelope of the object. Despite their virtually certain coevality, the differences in the components of the HH250-IRS binary are most likely due to markedly different circumstellar environments. The NW component displays an unevolved envelope, whereas dust growth and crystallization has taken place in the SE component. The weak or absent signatures of enstatite in the latter are fairly unusual among envelopes with crystalline dust, and we tentatively relate it to a possible wide gap or an inner truncation of the disk already hinted in previous observations by a drop in the $L'$-band flux, which might indicate that the SE component could actually be a very close binary. We speculate that the clear differences between the silicate feature spectra of both components of HH250-IRS may be due either to disk evolution sped up by multiplicity, or by accretion variability leading to episodes of crystal formation.

We present our observational results of the 1.1 mm continuum and the HCO+ (3-2) line in HL Tau at angular resolutions of 0.1" obtained with ALMA and our data analysis of the 2.9 mm and 1.1 mm continuum and the HCO+ (3-2) and (1-0) lines of the HL Tau disk. The Keplerian rotation of the HL Tau disk is well resolved in the HCO+ (3-2) emission, and the stellar mass is estimated to be 2.1+/-0.2 Msun with a disk inclination angle of 47 deg. The radial profiles of the HCO+ column density and excitation temperature are measured with the LTE analysis of the two transitions of the HCO+ emission. An HCO+ gas gap at a radius of 30 au, where the column density drops by a factor of 4-8, is found in the HCO+ column density profile, coincident with the dust gap traced by the continuum emission. No other clear HCO+ gas gap is seen. This HCO+ gas gap can be opened by a planet with mass of 0.5-0.8 Mj, which is comparable to the planet mass adopted in numerical simulations to form the dust gap at the same radius in the HL Tau disk. In addition to the disk component, a one-arm spiral with a length of ~3" (520 au) stretching out from the inner disk is observed in the HCO+ (3-2) emission. The observed velocity structures along the spiral suggest an infalling and rotating gas stream toward the inner disk.

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.

(Abbreviated) We investigate the interplay between magnetic (B) field, gravity, and turbulence in the fragmentation process of cores within the filamentary infrared dark cloud G34.43+00.24. We observe the magnetic field (B) morphology across G34.43 and compare with the kinematics obtained from N2H+ across the filament. We derive local velocity gradients from N2H+, tracing motion in the plane of sky, and compare with the observed local B field orientations in the plane of sky. Besides a large-scale east-west velocity gradient, we find a close alignment between local B field orientations and local velocity gradients toward the MM1/MM2 ridge. This local correlation in alignment suggests that gas motions are influenced by the magnetic field morphology or vice versa. Additionally, this alignment seems to be getting even closer with increasing integrated emission in N2H+, possibly indicating that a growing gravitational pull is more and more aligning B field and gas motion. We analyze and quantify B field, gravity, turbulence, and their relative importance toward the MM1, MM2 and MM3 regions with various techniques over two scales, a larger clump area at 2 pc scale and the smaller core area at 0.6 pc scale. While gravitational energy, B field, and turbulent pressure all grow systematically from large to small scale, the ratios among the three constituents develop clearly differently over scale. We propose that this varying relative importance between B field, gravity, and turbulence over scale drives and explains the different fragmentation types seen at sub-pc scale (no fragmentation in MM1; aligned fragmentation in MM2; clustered fragmentation in MM3). We discuss uncertainties, subtleties, and the robustness of our conclusion, and we stress the need of a multi-scale joint analysis to understand the dynamics in these systems.

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.

HL Tau is a Class I-II protostar embedded in an infalling and rotating envelope and possibly associated with a planet forming disk, and it is co-located in a 0.1 pc molecular cloud with two nearby young stellar objects. Our ALMA observations revealed two arc-like structures on a 1000 au scale connected to the disk, and their kinematics could not be explained with any conventional model of infalling and rotational motions. In this work, we investigate the nature of these arc-like structures connected to the HL Tau disk. We conducted new observations in the 13CO and C18O (3-2; 2-1) lines with JCMT and IRAM 30m, and obtained the ACA data with the 7-m array. With the single-dish, ACA, and ALMA data, we analyzed the gas motions on both 0.1 pc and 1000 au scales in the HL Tau region. We constructed new kinematical models of an infalling and rotating envelope with the consideration of relative motion between HL Tau and the envelope. By including the relative motion between HL Tau and its protostellar envelope, our kinematical model can explain the observed velocity features in the arc-like structures. The morphologies of the arc-like structures can also be explained with an asymmetric initial density distribution in our model envelope. In addition, our single-dish results support that HL Tau is located at the edge of a large-scale (0.1 pc) expanding shell driven by the wind or outflow from XZ Tau, as suggested in the literature. The estimated expanding velocity of the shell is comparable to the relative velocity between HL Tau and its envelope in our kinematical model. These results hints that the large-scale expanding motion likely impacts the protostellar envelope around HL Tau and affects its gas kinematics. We found that the mass infalling rate from the envelope onto the HL Tau disk can be decreased by a factor of two due to this impact by the large-scale expanding shell.

We present our analysis of the magnetic field structures from 6000 au to 100 au scales in the Class 0 protostar B335 inferred from our JCMT POL-2 observations and the ALMA archival polarimetric data. To interpret the observational results, we perform a series of (non-)ideal MHD simulations of the collapse of a rotating non-turbulent dense core, whose initial conditions are adopted to be the same as observed in B335, and generate synthetic polarization maps. The comparison of our JCMT and simulation results suggests that the magnetic field on a 6000 au scale in B335 is pinched and well aligned with the bipolar outflow along the east-west direction. Among all our simulations, the ALMA polarimetric results are best explained with weak magnetic field models having an initial mass-to-flux ratio of 9.6. However, we find that with the weak magnetic field, the rotational velocity on a 100 au scale and the disk size in our simulations are larger than the observational estimates by a factor of several. An independent comparison of our simulations and the gas kinematics in B335 observed with the SMA and ALMA favors strong magnetic field models with an initial mass-to-flux ratio smaller than 4.8. We discuss two possibilities resulting in the different magnetic field strengths inferred from the polarimetric and molecular-line observations, (1) overestimated rotational-to-gravitational energy in B335 and (2) additional contributions in the polarized intensity due to scattering on a 100 au scale.

We report our SMA observations of the Class \RN{1}-\RN{2} protostar HL Tau in the $^{13}$CO (2--1), C$^{18}$O (2--1), SO(5$_6$--4$_5$), and the 1.3 mm dust-continuum emission and our analyses of the ALMA long baseline data of the HCO$^{+}$ (1--0) emission. The 1.3 mm continuum emission observed with the SMA shows compact ($\sim$0$\farcs$8 $\times$ 0$\farcs$5) and extended ($\sim$6$\farcs$5 $\times$ 4$\farcs$3) components, tracing the protoplanetary disk and the protostellar envelope, respectively. The $^{13}$CO, C$^{18}$O, and HCO$^+$ show a compact ($\sim$ 200 AU) component at velocities higher than 3 km s$^{-1}$ from the systemic velocity and an extended ($\sim$ 1000 AU) component at the lower velocities. The high-velocity component traces the Keplerian rotating disk, and the low-velocity component traces the infalling envelope. The HCO$^+$ high-velocity component is fitted with a Keplerian disk model with a central stellar mass of 1.4 $M_{\odot}$. The radial intensity profiles of the $^{13}$CO and C$^{18}$O along the disk major axis are fitted with a disk+envelope model, and the gas masses of the disk and envelope are estimated to be $2\mbox{--}40\times10^{-4}$ $M_{\odot}$ and $2.9\times10^{-3}$ $M_{\odot}$, respectively. The disk dust mass has been estimated to be $1\mbox{--}3 \times 10^{-3}$ $M_{\odot}$ in the literature. Thus, our estimated disk gas mass suggests that the gas-to-dust mass ratio in the disk is $<$10, a factor of ten lower than the estimated ratio in the envelope. We discuss the possible gas depletion or CO depletion in the planet-forming candidate HL Tau in the context of disk and envelope evolution.

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

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.

In recent ALMA surveys, the gas distributions and velocity structures of most of the protoplanetary disks can still not be imaged at high S/N due to the short integration time. In this work, we re-analyzed the ALMA 13CO (3-2) and C18O (3-2) data of 88 young stellar objects in Lupus with the velocity-aligned stacking method to enhance S/N and to study the kinematics and disk properties traced by molecular lines. This method aligns spectra at different positions in a disk based on the projected Keplerian velocities at their positions and then stacks them. This method enhances the S/N ratios of molecular-line data and allows us to obtain better detections and to constrain dynamical stellar masses and disk orientations. We obtain 13CO detections in 41 disks and C18O detections in 18 disks with 11 new detections in 13CO and 9 new detections in C18O after applying the method. We estimate the disk orientations and the dynamical stellar masses from the 13CO data. Our estimated dynamical stellar masses correlate with the spectroscopic stellar masses, and in a subsample of 16 sources, where the inclination angles are better constrained, the two masses are in a good agreement within the uncertainties and with a mean difference of 0.15 Msun. With more detections of fainter disks, our results show that high gas masses derived from the 13CO and C18O lines tend to be associated with high dust masses estimated from the continuum emission. Nevertheless, the scatter is large (0.9 dex), implying large uncertainties in deriving the disk gas mass from the line fluxes. We find that with such large uncertainties it is expected that there is no correlation between the disk gas mass and the mass accretion rate with the current data. Deeper observations to detect disks with gas masses <1E-5 Msun in molecular lines are needed to investigate the correlation between the disk gas mass and the mass accretion rate.

Ambipolar diffusion can cause a velocity drift between ions and neutrals. This is one of the non-ideal MHD effects proposed to enable the formation of Keplerian disks with sizes of tens of au. To observationally study ambipolar diffusion in collapsing protostellar envelopes, we compare gas kinematics traced by the H13CO+ (3-2) and C18O (2-1) lines in the Class 0 protostar B335 obtained with our ALMA observations. A central compact (~1"-2") component that is elongated perpendicular to the outflow direction and exhibits a clear velocity gradient along the outflow direction is observed in both lines and most likely traces the infalling flattened envelope. We constructed kinematical models to fit the observed velocity structures and to measure the infalling velocities of the ionized and neutral gas on a 100 au scale in B335. The infalling velocities in the H13CO+ and C18O emission are both measured to be 0.85+/-0.2 km/s at a radius of 100 au, suggesting that the velocity drift between the ionized and neutral gas is at most 0.3 km/s at a radius of 100 au. The Hall parameter for H13CO+ is estimated to be >>1 on a 100 au scale in B335, so that H13CO+ is expected to be attached to the magnetic field. Our non-detection or upper limit of the velocity drift between ions and neutrals could suggest that the magnetic field remains rather well coupled to the bulk neutral material on a 100 au scale in this source, and that any significant field-matter decoupling, if present, likely occurs only on a smaller scale, leading to an accumulation of magnetic flux and thus efficient magnetic braking in the inner envelope. However, because of our limited angular resolution, we cannot rule out a significant ambipolar drift only in the midplane of the infalling envelope. Future observations with higher angular resolutions (~0.1") are needed to examine this possibility and ambipolar diffusion on a smaller scale.

We present the first ALMA dust polarization observations towards the high-mass star-forming regions W51 e2, e8, and W51 North in Band 6 (230 GHz) with a resolution around 0.26" ($\sim5$mpc). Polarized emission in all three sources is clearly detected and resolved. Measured relative polarization levels are between 0.1\% and 10\%. While the absolute polarization shows complicated structures, the relative polarization displays the typical anti-correlation with Stokes $I$, though with a large scatter. Inferred magnetic (B) field morphologies are organized and connected. Detailed substructures are resolved, revealing new features such as cometary-shaped B-field morphologies in satellite cores, symmetrically converging B-field zones, and possibly streamlined morphologies. The local B-field dispersion shows some anti-correlation with the relative polarization. Moreover, lowest polarization percentages together with largest dispersions coincide with B-field convergence zones. We put forward $\sin\omega$, where $\omega$ is the measurable angle between a local B-field orientation and local gravity, as a measure of how effectively the B-field can oppose gravity. Maps of $\sin\omega$ for all three sources show organized structures that suggest a locally varying role of the B-field, with some regions where gravity can largely act unaffectedly, possibly in a network of narrow magnetic channels, and other regions where the B-field can work maximally against gravity.

How important is the magnetic (B-) field when compared to gravity and turbulence in the star-formation process? Does its importance depend on scale and location? We summarize submm dust polarization observations towards the large filamentary infrared dark cloud G34 and towards a dense core in the high-mass star-forming region W51. We detect B-field orientations that are either perpendicular or parallel to the G34 filament axis. These B-field orientations further correlate with local velocity gradients. Towards three cores in G34 we find a varying importance between B-field, gravity, and turbulence that seems to dictate varying types of fragmentation. At highest resolution towards the gravity-dominated collapsing core W51 e2 we resolve new B-field features, such as converging B-field lines and possibly magnetic channels.

Herbig-Haro flows are signposts of recent major accretion and outflow episodes. We aim to determine the nature and properties of the little-known outflow source HH 250-IRS, which is embedded in the Aquila clouds. We have obtained adaptive optics-assisted L-band images with the NACO instrument on the Very Large Telescope (VLT), together with N- and Q-band imaging with VISIR also on the VLT. Using the SINFONI instrument on the VLT we carried out H- and K-band integral field spectroscopy of HH 250-IRS, complemented with spectra obtained with the SpeX instrument at the InfraRed Telescope Facility (IRTF) in the JHKL bands. Finally, the SubMillimeter Array (SMA) interferometer was used to study the circumstellar environment of HH 250-IRS at 225 and 351 GHz with CO (2-1) and CO (3-2) maps and 0.9 mm and 1.3 mm continuum images. The HH 250-IRS source is resolved into a binary with 0''53 separation, corresponding to 120 AU at the adopted distance of 225 pc. The individual components show heavily veiled spectra with weak CO absorption indicative of late-type stars. Both are Class I sources, but their spectral energy distributions between 1.5 $\mu$m and 19 $\mu$m differ markedly and suggest the existence of a large cavity around one of the components. The millimeter interferometric observations indicate that the gas mainly traces a circumbinary envelope or disk, while the dust emission is dominated by one of the circumstellar envelopes. HH 250-IRS is a new addition to the handful of multiple systems where the individual stellar components, the circumstellar disks and a circumbinary disk can be studied in detail, and a rare case among those systems in which a Herbig-Haro flow is present.

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.

The interaction between a supermassive black hole (SMBH) and the surrounding material is of primary importance in modern astrophysics. The detection of the molecular 2-pc circumnuclear disk (CND) immediately around the Milky Way SMBH, SgrA*, provides an unique opportunity to study SMBH accretion at sub-parsec scales. Our new wide-field CS(2-1) map toward the Galactic center (GC) reveals multiple dense molecular streamers originated from the ambient clouds 20-pc further out, and connecting to the central 2 parsecs of the CND. These dense gas streamers appear to carry gas directly toward the nuclear region and might be captured by the central potential. Our phase-plot analysis indicates that these streamers show a signature of rotation and inward radial motion with progressively higher velocities as the gas approaches the CND and finally ends up co-rotating with the CND. Our results might suggest a possible mechanism of gas feeding the CND from 20 pc around 2 pc in the GC. In this paper, we discuss the morphology and the kinematics of these streamers. As the nearest observable Galactic nucleus, this feeding process may have implications for understanding the processes in extragalactic nuclei.

The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8" (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2" (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_\odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_\odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau.

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.

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 report the ALMA Cycle 2 observations of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), 13CO (3-2), SO (7_8-6_7), and the CS (7-6) emission. At 0.18" (= 25 AU) resolution, ~4-times higher than that of our Cycle 0 observations, the circumbinary disk as seen in the 0.9-mm emission is shown to be comprised of a northern and a southern spiral arm, with the southern arm connecting to the circumstellar disk around Source B. The western parts of the spiral arms are brighter than the eastern parts, suggesting the presence of an m=1 spiral mode. In the C18O emission, the infall gas motions in the inter-arm regions and the outward gas motions in the arms are identified. These observed features are well reproduced with our numerical simulations, where gravitational torques from the binary system impart angular momenta to the spiral-arm regions and extract angular momenta from the inter-arm regions. Chemical differentiation of the circumbinary disk is seen in the four molecular species. Our Cycle 2 observations have also resolved the circumstellar disks around the individual protostars, and the beam-deconvolved sizes are 0.29" X 0.19" (= 40 X 26 AU) (P.A. = 144 deg) and 0.26" X 0.20" (= 36 X 27 AU) (P.A. = 147 deg) for Sources A and B, respectively. The position and inclination angles of these circumstellar disks are misaligned with that of the circumbinary disk. The C18O emission traces the Keplerian rotation of the misaligned disk around Source A.

We present ALMA 1.3 mm continuum, 12CO, C18O, and SO data for the Class 0 protostars, Lupus 3 MMS, IRAS 15398-3559, and IRAS 16253-2429 at resolutions of ~100 AU. By measuring a rotational profile in C18O, a 100 AU Keplerian disk around a 0.3 Msun protostar is observed in Lupus 3 MMS. No 100 AU Keplerian disks are observed in IRAS 15398-3559 and IRAS 16253-2429. Nevertheless, embedded compact (<30 AU) continuum components are detected. The C18O emission in IRAS 15398-3559 shows signatures of infall with a constant angular momentum. IRAS 16253-2429 exhibits signatures of infall and rotation, but its rotational profile is unresolved. By fitting the C18O data with our kinematic models, the protostellar masses and the disk radii are inferred to be 0.01 Msun and 20 AU in IRAS 15398-3559, and 0.03 Msun and 6 AU in IRAS 16253-2429. By comparing the specific angular momentum profiles from 10,000 to 100 AU in 8 Class 0 and I protostars, we find that the evolution of envelope rotation can be described with conventional inside-out collapse models. In comparison with a sample of 18 protostars with known disk radii, our results reveal signs of disk growth, with the disk radius increasing as M*^{0.8+/-0.14} or t^{1.09+\-0.37} in the Class 0 stage, where M* is the protostellar mass and t is the age. The disk growth rate slows down in the Class I stage. Besides, we find a hint that the mass accretion rate declines as t^{-0.26+\-0.04} from the Class 0 to I stages.

We introduce a new stacking method in Keplerian disks that (1) enhances signal-to-noise ratios (S/N) of detected molecular lines and (2) that makes visible otherwise undetectable weak lines. Our technique takes advantage of the Keplerian rotational velocity pattern. It aligns spectra according to their different centroid velocities at their different positions in a disk and stacks them. After aligning, the signals are accumulated in a narrower velocity range as compared to the original line width without alignment. Moreover, originally correlated noise becomes de-correlated. Stacked and aligned spectra, thus, have a higher S/N. We apply our method to ALMA archival data of DCN (3-2), DCO+ (3-2), N2D+ (3-2), and H2CO (3_0,3-2_0,2), (3_2,2-2_2,1), and (3_2,1-2_2,0) in the protoplanetary disk around HD 163296. As a result, (1) the S/N of the originally detected DCN (3-2), DCO+ (3-2), and H2CO (3_0,3-2_0,2) and N2D+ (3-2) lines are boosted by a factor of >4-5 at their spectral peaks, implying one order of magnitude shorter integration times to reach the original S/N; and (2) the previously undetectable spectra of the H2CO (3_2,2-2_2,1) and (3_2,1-2_2,0) lines are materialized at more than 3 sigma. These dramatically enhanced S/N allow us to measure intensity distributions in all lines with high significance. The principle of our method can not only be applied to Keplerian disks but also to any systems with ordered kinematic patterns.

We conducted IRAM-30m C18O (2-1) and SMA 1.3mm continuum, 12CO (2-1), and C18O (2-1) observations toward the Class 0/I protostar L1455 IRS1 in Perseus. The IRAM results show L1455 IRS1 located in a dense core of 0.05 pc in size with a mass of 0.54M_sun, connecting to a filamentary structure. Inside the dense core, compact components of 350 AU and 1500 AU are detected in the SMA 1.3 mm continuum and C18O, with a velocity gradient in the latter one perpendicular to a bipolar outflow in 12CO, likely tracing a rotational motion. We measure a rotational velocity profile ~ r^-0.75 that becomes shallower at a turning radius of ~200 AU which is approximately the radius of the 1.3 mm continuum component. These results hint the presence of a Keplerian disk with a radius <200 AU around L1455 IRS1 with a protostellar mass of about 0.28 M_sun. We derive a core rotation that is about one order of magnitude faster than expected. A significant velocity gradient along a filament towards IRS1 indicates that this filament is dynamically important, providing a gas reservoir and possibly responsible for the faster-than-average core rotation. Previous polarimetric observations show a magnetic field aligned with the outflow axis and perpendicular to the associated filament on a 0.1~pc scale, while on the inner 1000 AU scale, the field becomes perpendicular to the outflow axis.This change in magnetic field orientations is consistent with our estimated increase in rotational energy from large to small scales that overcomes the magnetic field energy, wrapping the field lines and aligning them with the disk velocity gradient. These results are discussed in the context of the interplay between filament, magnetic field, and gas kinematics from large to small scale. Possible emerging trends are explored with a sample of 8 Class 0/I protostars.

We have analyzed the HCO+ (1-0) data of the Class I-II protostar, HL Tau, obtained from the Atacama Large Millimeter/Submillimeter Array long baseline campaign. We generated the HCO+ image cube at an angular resolution of ~0.07 (~10 AU), and performed azimuthal averaging on the image cube to enhance the signal-to-noise ratio and measure the radial profile of the HCO+ integrated intensity. Two gaps at radii of ~28 AU and ~69 AU and a central cavity are identified in the radial intensity profile. The inner HCO+ gap is coincident with the millimeter continuum gap at a radius of 32 AU. The outer HCO+ gap is located at the millimeter continuum bright ring at a radius of 69 AU and overlaps with the two millimeter continuum gaps at radii of 64 AU and 74 AU. On the contrary, the presence of the central cavity is likely due to the high optical depth of the 3 mm continuum emission and not the depletion of the HCO+ gas. We derived the HCO+ column density profile from its intensity profile. From the column density profile, the full-width-half-maximum widths of the inner and outer HCO+ gaps are both estimated to be ~14 AU, and their depths are estimated to be ~2.4 and ~5.0. These results are consistent with the expectation from the gaps opened by forming (sub-)Jovian mass planets, while placing tight constraints on the theoretical models solely incorporating the variation of dust properties and grain sizes.

We highlight distinct and systematic observational features of magnetic field morphologies in polarized submm dust continuum. We illustrate this with specific examples and show statistical trends from a sample of 50 star-forming regions.

CIT 6 is a carbon star in the transitional phase from the asymptotic giant branch (AGB) to the protoplanetary nebulae (pPN). Observational evidences of two point sources in the optical, circumstellar arc segments in an HC$_3$N line emission, and a bipolar nebula in near-infrared provide strong support for the presence of a binary companion. Hence, CIT 6 is very attractive for studying the role of companions in the AGB-pPN transition. We have carried out high resolution $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ observations of CIT 6 with the Submillimeter Array combined with the Submillimeter Telescope (single-dish) data. The $^{12}$CO channel maps reveal a spiral-shell pattern connecting the HC$_3$N segments in a continuous form, and an asymmetric outflow corresponding to the near-infrared bipolar nebula. Rotation of the $^{12}$CO channel peak position may be related to the inner spiral winding and/or the bipolar outflow. An eccentric orbit binary is suggested for the presences of an anisotropic mass loss to the west and a double spiral pattern. The lack of interarm emission to the west may indicate a feature corresponding to the periastron passage of a highly eccentric orbit of the binary. Spatially-averaged radial and spectral profiles of $^{12}$CO $J=2-1$ and $^{13}$CO $J=2-1$ are compared with simple spherical radiative transfer models, suggesting a change of $^{12}$CO/$^{13}$CO abundance ratio from $\sim30$ to $\sim50$ inward in the CSE of CIT 6. The millimeter continuum emission is decomposed into extended dust thermal emission (spectral index $\sim-2.4$) and compact emission from radio photosphere (spectral index $\sim-2.0$).

We have conducted ALMA cycle 2 observations in the 1.3 mm continuum and in the C18O (2-1) and SO (5_6-4_5) lines at a resolution of ~0.3" toward the Class 0 protostar B335. The 1.3 mm continuum, C18O, and SO emission all show central compact components with sizes of ~40-180 AU within more extended components. The C18O component shows signs of infalling and rotational motion. By fitting simple kinematic models to the C18O data, the protostellar mass is estimated to be 0.05 Msun. The specific angular momentum, on a 100 AU scale, is ~4.3E-5 km/s*pc. A similar specific angular momentum, ~3E-5 to 5E-5 km/s*pc, is measured on a 10 AU scale from the velocity gradient observed in the central SO component, and there is no clear sign of an infalling motion in the SO emission. By comparing the infalling and rotational motion, our ALMA results suggest that the observed rotational motion has not yet reached Keplerian velocity neither on a 100 AU nor even on a 10 AU scale. Consequently, the radius of the Keplerian disk in B335 (if present) is expected to be 1-3 AU. The expected disk radius in B335 is one to two orders of magnitude smaller than those of observed Keplerian disks around other Class 0 protostars. Based on the observed infalling and rotational motion from 0.1 pc to inner 100 AU scales, there are two possible scenarios to explain the presence of such a small Keplerian disk in B335: magnetic braking and young age. If our finding is the consequence of magnetic braking, ~50% of the angular momentum of the infalling material within a 1000 AU scale might have been removed, and the magnetic field strength on a 1000 AU scale is estimated to be ~200 uG. If it is young age, the infalling radius in B335 is estimated to be ~2700 AU, corresponding to a collapsing time scale of ~5E4 yr.

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.

We perform imaging and analyses of SMA 1.3 mm continuum, C18O (2-1) and 12CO (2-1) line data of 17 Class 0 and 0/I protostars to study their gas kinematics on a 1,000-AU scale. Continuum and C18O (2-1) emission are detected toward all the sample sources and show central primary components with sizes of ~600-1,500 AU associated with protostars. The velocity gradients in C18O (2-1) have wide ranges of orientations from parallel to perpendicular to the outflows, with magnitudes from ~1 to ~530 km/s/pc. We construct a simple kinematic model to reproduce the observed velocity gradients, estimate the infalling and rotational velocities, and infer the disk radii and the protostellar masses. The inferred disk radii range from <5 AU to >500 AU with estimated protostellar masses from <0.1 Msun to >1 Msun. Our results hint that both large and small disks are possibly present around Class 0 protostars, which could be a sign of disk growth at the Class 0 stage. In addition, the directions of the overall velocity gradients in 7 out of the 17 sources are close to perpendicular to their outflow axes, which is a signature of significant rotational motions. From our model fitting, the specific angular momenta in these sources are estimated to be >2E-4 km/s*pc, suggesting that magnetic braking is unlikely efficient on a 1,000-AU scale in these Class 0 and 0/I sources. In a sub-sample with observed magnetic field orientations, we find no source with large specific angular momenta together with closely aligned magnetic field and outflow axes. This possibly hints that the magnetic field, if originally aligned with the rotational axis, can play a role in removing angular momentum from infalling material at the Class 0 stage. We discuss our results in comparison with theoretical models of collapsing dense cores with and without magnetic fields in the context of disk formation.

Submillimeter dust polarization measurements of a sample of 50 star-forming regions, observed with the SMA and the CSO covering pc-scale clouds to mpc-scale cores, are analyzed in order to quantify the magnetic field importance. The magnetic field misalignment $\delta$ -- the local angle between magnetic field and dust emission gradient -- is found to be a prime observable, revealing distinct distributions for sources where the magnetic field is preferentially aligned with or perpendicular to the source minor axis. Source-averaged misalignment angles $\langle|\delta|\rangle$ fall into systematically different ranges, reflecting the different source-magnetic field configurations. Possible bimodal $\langle|\delta|\rangle$-distributions are found for the separate SMA and CSO samples. Combining both samples broadens the distribution with a wide maximum peak at small $\langle|\delta|\rangle$-values. Assuming the 50 sources to be representative, the prevailing source-magnetic field configuration is one that statistically prefers small magnetic field misalignments $|\delta|$. When interpreting $|\delta|$ together with an MHD force equation, as developed in the framework of the polarization-intensity gradient method, a sample-based scaling fits the magnetic field tension-to-gravity force ratio $\langle\Sigma_B\rangle$ versus $\langle|\delta|\rangle$ with $\langle\Sigma_B\rangle = 0.116 \cdot \exp(0.047\cdot \langle|\delta|\rangle)\pm 0.20$ (mean error), providing a way to estimate the relative importance of the magnetic field, only based on measurable field misalignments $|\delta|$. The force ratio $\Sigma_B$ discriminates systems that are collapsible on average ($\langle \Sigma_B\rangle <1$) from other molecular clouds where the magnetic field still provides enough resistance against gravitational collapse ($\langle \Sigma_B\rangle >1$) (abridged).

We present combined SubMillimeter Array (SMA) + Atacama Submillimeter Telescope Experiment (ASTE) images of the Class I protobinary L1551 IRS 5 in the CS ($J$ = 7--6) line, the submillimeter images of L1551 IRS 5 with the most complete spatial sampling ever achieved ($0''.9$ -- $36''$). The SMA image of L1551 IRS 5 in the 343 GHz dust-continuum emission is also presented, which shows an elongated feature along the northwest to southeast direction ($\sim$160 AU $\times$ 80 AU), perpendicular to the associated radio jets. The combined SMA+ASTE images show that the high-velocity ($\gtrsim$1.5 km s$^{-1}$) CS emission traces the structure of the dust component and shows a velocity gradient along the major axis, which is reproduced by a geometrically-thin Keplerian-disk model with a central stellar mass of $\sim$0.5 $M_{\odot}$. The low-velocity ($\lesssim$1.3 km s$^{-1}$) CS emission shows an extended ($\sim$1000 AU) feature that exhibits slight south (blueshifted) to north (redshifted) emission offsets, which is modeled with a rotating and infalling envelope with a conserved angular momentum. The rotational motion of the envelope connects smoothly to the inner Keplerian rotation at a radius of $\sim$64 AU. The infalling velocity of the envelope is $\sim$three times lower than the free-fall velocity toward the central stellar mass of 0.5 $M_{\odot}$. These results demonstrate transition from the infalling envelope to the Keplerian disk, consistent with the latest theoretical studies of disk formation. We suggest that sizable ($r\sim$50--200 AU) Keplerian disks are already formed when the protostars are still deeply embedded in the envelopes.

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.

We performed SMA observations in the C18O (2-1) emission line toward six Class 0 and I protostars, to study rotational motions of their surrounding envelopes and circumstellar material on 100 to 1000 AU scales. C18O (2-1) emission with intensity peaks located at the protostellar positions is detected toward all the six sources. The rotational velocities of the protostellar envelopes as a function of radius were measured from the Position-Velocity diagrams perpendicular to the outflow directions passing through the protostellar positions. Two Class 0 sources, B335 and NGC 1333 IRAS 4B, show no detectable rotational motion, while L1527 IRS (Class 0/I) and L1448-mm (Class 0) exhibit rotational motions with radial profiles of Vrot ~ r^{-1.0+/-0.2} and ~ r^{-1.0+/-0.1}, respectively. The other Class I sources, TMC-1A and L1489 IRS, exhibit the fastest rotational motions among the sample, and their rotational motions have flatter radial profiles of Vrot ~ r^{-0.6+/-0.1} and ~ r^{-0.5+/-0.1}, respectively. The rotational motions with the radial dependence of ~ r^{-1} can be interpreted as rotation with a conserved angular momentum in a dynamically infalling envelope, while those with the radial dependence of ~ r^{-0.5} can be interpreted as Keplerian rotation. These observational results demonstrate categorization of rotational motions from infalling envelopes to Keplerian-disk formation. Models of the inside-out collapse where the angular momentum is conserved are discussed and compared with our observational results.

We made C18O (2-1) and CS (7-6) images of the protostellar envelope around B335 with a high spatial dynamic range from ~10000 to ~400 AU, by combining the Submillimeter Array and single-dish data. The C18O emission shows an extended (~10000 AU) structure as well as a compact (~1500 AU) component concentrated at the protostellar position. The CS emission shows a compact (~900 AU) component surrounding the protostar, plus a halo-like (~3000 AU) structure elongated along the east-west direction. At higher velocities (|dV| >~0.3 km s^-1), the CS emission is stronger and more extended than the C18O emission. Physical conditions of the envelope were examined through an LVG model. At |dV| >~0.3 km s^-1, the gas temperature is higher (>40 K) than that at |dV| <~0.3 km s^-1, whereas the gas density is lower (<10^6 cm^-3). We consider that the higher-temperature and lower-density gas at |dV| >~0.3 km s^-1 is related to the associated outflow, while the lower-temperature and higher-density gas at |dV| <~0.3 km s^-1 is the envelope component. From the inspection of the positional offsets in the velocity channel maps, the radial profile of the specific angular momentum of the envelope rotation in B335 was revealed at radii from ~10^4 down to ~10^2 AU. The specific angular momentum decreases down to the radius of ~370 AU, and then appears to be conserved within that radius. A possible scenario of the evolution of envelope rotation is discussed.

We have performed detailed imaging and analyses of SMA observations in 230 GHz continuum, 12CO (2-1), 13CO (2-1), and C18O (2-1) emission toward B335, a Bok globule with an embedded Class 0 source (Lbol ~ 1.5 Lsun). We report the first discover of high-velocity 12CO jets with a size of ~ 900 AU x 1500 AU along the E-W direction in B335. The estimated mass-loss rate (~ 2.3 x 10^-7 Msun yr^-1) and the momentum flux (~ 3.7 x 10^-5 Msun yr^-1 km s^-1) of the 12CO jets in B335 are one order of magnitude lower than those of other 12CO jets in more luminous sources such as HH 211 (Lbol ~ 3.6 Lsun) and HH 212 (Lbol ~ 14 Lsun). The weaker jet activity in B335 could be due to the lower active accretion onto the central protostar. The C18O emission shows a compact (~ 1500 AU) condensation associated with the central protostar, and it likely traces the protostellar envelope around B335. The envelope exhibits a velocity gradient from the east (blueshifted) to west (redshifted) that can be interpreted as an infalling motion. The estimated central stellar mass, the mass infalling rate, and the accretion luminosity are 0.04 Msun, 6.9 x 10^-6 Msun yr^-1, and 2.1 Lsun, respectively. On the other hand, there is no clear velocity gradient perpendicular to the outflow axis in the C18O envelope, suggesting little envelope rotation on a hundred-AU scale. The upper limits of the rotational velocity and specific angular momentum were estimated to be 0.04 km s^-1 and 7.0 x 10^-5 km s^-1 pc at a radius of 370 AU, respectively. The specific angular momentum and the inferred Keplerian radius (~ 6 AU) in B335 are 1 - 2 orders of magnitude smaller than those in other more-evolved sources. Possible scenarios to explain the lower specific angular momentum in B335 are discussed.