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 present measurements of the thermal Sunyaev-Zel'dovich effect (SZE) toward SPT-CL J2334-4243 (the Phoenix galaxy cluster) at z=0.597 by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. The SZE is imaged at 5" resolution (corresponding to the physical scale of 23kpc/h) within 200kpc/h from the central galaxy with the peak signal-to-noise ratio exceeding 11. Combined with the Chandra X-ray image, the ALMA SZE data further allow for non-parametric deprojection of electron temperature, density, and entropy. Our method can minimize contamination by the central AGN and the X-ray absorbing gas within the cluster, both of which largely affect the X-ray spectrum. We find no significant asymmetry or disturbance in the SZE image within the current measurement errors. The detected SZE signal agrees well with the average pressure profile of local cool-core clusters. Unlike any other known clusters, however, gas temperature drops by at least a factor of 5 toward the center. We identify ~6x10^{11} M_sun cool gas with temperature ~3keV in the inner 20kpc/h. Taken together, our results imply that the gas is indeed cooling efficiently and nearly isobarically down to this radius in the Phoenix cluster.

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 report ALMA observations of NGC 1333 IRAS 4A, a young low-mass protostellar binary, referred as 4A1 and 4A2. With multiple H$_2$CO transitions and HNC (4$-$3) observed at a resolution of 0.25" ($\sim$70 au), we investigate the gas kinematics of 4A1 and 4A2. Our results show that, on the large angular scale ($\sim$10"), 4A1 and 4A2 each display a well-collimated outflow along the N-S direction, and an S-shaped morphology is discerned in the outflow powered by 4A2. On the small scale ($\sim$0.3"), 4A1 and 4A2 exhibit distinct spectral features toward the continuum centroid, with 4A1 showing simple symmetric profiles predominantly in absorption and 4A2 demonstrating rather complicated profiles in emission as well as in absorption. Based on radiative transfer modeling exercises, we find that the physical parameters inferred from earlier low-resolution observations cannot be directly extrapolated down to the 4A1 inner region. Possible reasons for the discrepancies between the observed and modelled profiles are discussed. We constrain the mass infall rate in 4A1 to be at most around 3$\times$10$^{-5}$ M$_{\odot}$ year$^{-1}$ at the layer of 75 au. For the kinematics of the 4A2 inner envelope, the absorbing dips in the H$_2$CO spectra are skewed toward the redshifted side and likely signatures of inward motion. These absorbing dips are relatively narrow. This is, like the case for 4A1, significantly slower than the anticipated inflow speed. We estimate a mass infall rate of 3.1$-$6.2 $\times$ 10$^{-5}$ M$_\odot$ year$^{-1}$ at the layer of 100 au in 4A2.

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

We analyze $^{12}$CO($J$=2-1), $^{13}$CO($J$=2-1), C$^{18}$O ($J$=2-1), and 1.3 mm continuum maps of the $\rho$ Ophiuchi A photo-dissociation region (PDR) obtained with ALMA. Layered structures of the three CO isotopologues with an angular separation of 10 arcsec = 6.6$\times$10$^{-3}$ pc = 1400 au are clearly detected around the Be star, S1 (i.e., each front of emission shifts from the near to far side in order of $^{12}$CO, $^{13}$CO, and C$^{18}$O). We estimate the spatial variations of $X$($^{13}$CO)/$X$(C$^{18}$O) abundance ratios, and find that the abundance ratio is as high as 40 near the emission front, and decreases to the typical value in the solar system of 5.5 in a small angular scale of 4 arcsec = 2.6$\times$10$^{-3}$ pc = 560 au. We also find that the $I$($^{12}$CO(2-1))/$I$($^{13}$CO(2-1)) intensity ratio is very high ($>$21) in the flat-spectrum young stellar object, GY-51, located in the PDR. The enhancement of the ratios indicates that the UV radiation significantly affects the CO isotopologues via selective dissociation in the overall $\rho$ Ophiuchi A PDR, and that the $\rho$ Ophiuchi A PDR has a plane-parallel structure.

We report high angular resolution observations of NGC1333 IRAS4A, a protostellar binary including A1 and A2, at 0.84 mm with the Atacama Large Millimeter/submillimeter Array. From the continuum observations, we suggest that the dust emission from the A1 core is optically thick, and A2 is predominantly optically thin. The A2 core, exhibiting a forest of spectral lines including complex molecules, is a well known hot corino as suggested by previous works. More importantly, we report, for the first time, the solid evidence of complex organic molecules (COMs), including CH3OH, ^13CH3OH, CH2DOH, CH3CHO associated with the A1 core seen in absorption. The absorption features mostly arise from a compact region around the continuum peak position of the A1 core. Rather than originating from a larger common envelope surrounding the protobinary, the COM features are associated with individual cores A1 and A2. Considering the signatures observed in both continuum and spectral lines, we propose two possible scenarios for IRAS 4A1 - the COM absorption lines may arise from a hot-corino-like atmosphere at the surface of an optically-thick circumstellar disk around A1, or the absorption may arise from different layers of a temperature-stratified dense envelope.

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 investigate the gas structures around young binary stars by using three-dimensional numerical simulations. Each model exhibits circumstellar disks, spiral arms, and a circumbinary disk with an inner gap or cavity. The circumbinary disk has an asymmetric pattern rotating at an angular velocity of approximately one-fourth of the binary orbit of the moderate-temperature models. Because of this asymmetry, the circumbinary disk has a density bump and a vortex, both of which continue to exist until the end of our calculation. The density bump and vortex are attributed to enhanced angular momentum, which is promoted by the gravitational torque of the stars. In a hot model ($c \ge 2.0$), the asymmetry rotates considerably more slowly than in the moderate-temperature models. The cold models ($c \le 0.02$) exhibit eccentric circumbinary disks, the precession of which is approximated by a secular motion of the ballistic particles. The asymmetry in the circumbinary disk does not depend on the mass ratio, but it becomes less clear as the specific angular momentum of the infalling envelope increases. The relative accretion rate onto the stars is sensitive to the angular momentum of the infalling envelope. For envelopes with constant angular momentum, the secondary tends to have a higher accretion rate than the primary, except in very low angular momentum cases. For envelopes with a constant angular velocity, the primary has a higher accretion rate than the secondary because gas with low specific angular momentum falls along the polar directions.

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.

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

RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern characteristic of gas sloshing, whereas we find no significant residual in the Sunyaev-Zel'dovich effect (SZE) image from the Atacama Large Millimeter/submillimeter Array (ALMA). We estimate the equation of state of perturbations in the gas from the X-ray and SZE residual images. The inferred velocity is 420 +310 -420 km s-1, which is much lower than the adiabatic sound speed of the intracluster medium in the core. We thus conclude that the perturbation is nearly isobaric, and gas sloshing motion is consistent with being in pressure equilibrium. Next, we report evidence for gas stripping of an infalling subcluster, which likely shock-heats gas to high temperature well in excess of 20 keV. Using mass distribution inferred from strong lensing images of the Hubble Space Telescope (HST), we find that the mass peak is located away from the peak position of stripped gas with statistical significance of > 5{\sigma}. Unlike for the gas sloshing, the velocity inferred from the equation of state of the excess hot gas is comparable to the adiabatic sound speed expected for the 20 keV intracluster medium. All of the results support that the southeast substructure is created by a merger. On the other hand, the positional offset between the mass and the gas limits the self-interaction cross section of dark matter to be less than 3.7 h-1 cm2 g-1 (95% CL).

We present results from our SMA observations and data analyses of the SMA archival data of the Class I protostar IRAS 04169+2702. The high-resolution (~0.5") $^{13}$CO (3-2) image cube shows a compact ($r$ ~< 100 au) structure with a northwest (blue) to southeast (red) velocity gradient, centered on the 0.9-mm dust-continuum emission. The direction of the velocity gradient is orthogonal to the axis of the molecular outflow as seen in the SMA $^{12}$CO (2-1) data. A similar gas component is seen in the SO (6$_5$-5$_4$) line. On the other hand, the C$^{18}$O (2-1) emission traces a more extended ($r$ ~400 au) component with the opposite, northwest (red) to southeast (blue) velocity gradient. Such opposite velocity gradients in the different molecular lines are also confirmed from direct fitting to the visibility data. We have constructed models of a forward-rotating and counter-rotating Keplerian disk and a protostellar envelope, including the SMA imaging simulations. The counter-rotating model could better reproduce the observed velocity channel maps, although we could not obtain statistically significant fitting results. The derived model parameters are; Keplerian radius of 200 au, central stellar mass of 0.1 $M_{solar}$, and envelope rotational and infalling velocities of 0.20 km s$^{-1}$ and 0.16 km s$^{-1}$, respectively. One possible interpretation for these results is the effect of the magnetic field in the process of disk formation around protostars, $i.e.$, Hall effect.

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.

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.

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.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

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 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 present observations at 7 mm that fully resolve the two circumstellar disks, and a reanalyses of archival observations at 3.5 cm that resolve along their major axes the two ionized jets, of the class I binary protostellar system L1551 NE. We show that the two circumstellar disks are better fit by a shallow inner and steep outer power-law than a truncated power-law. The two disks have very different transition radii between their inner and outer regions of $\sim$18.6 AU and $\sim$8.9 AU respectively. Assuming that they are intrinsically circular and geometrically thin, we find that the two circumstellar disks are parallel with each other and orthogonal in projection to their respective ionized jets. Furthermore, the two disks are closely aligned if not parallel with their circumbinary disk. Over an interval of $\sim$10 yr, source B (possessing the circumsecondary disk) has moved northwards with respect to and likely away from source A, indicating an orbital motion in the same direction as the rotational motion of their circumbinary disk. All the aforementioned elements therefore share the same axis for their angular momentum, indicating that L1551 NE is a product of rotationally-driven fragmentation of its parental core. Assuming a circular orbit, the relative disk sizes are compatible with theoretical predictions for tidal truncation by a binary system having a mass ratio of $\sim$0.2, in agreement with the reported relative separations of the two protostars from the center of their circumbinary disk. The transition radii of both disks, however, are a factor of $\gtrsim$1.5 smaller than their predicted tidally-truncated radii.

We present the first image of the thermal Sunyaev-Zel'dovich effect (SZE) obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Combining 7-m and 12-m arrays in Band 3, we create an SZE map toward a galaxy cluster RXJ1347.5-1145 with 5 arc-second resolution (corresponding to the physical size of 20 kpc/h), the highest angular and physical spatial resolutions achieved to date for imaging the SZE, while retaining extended signals out to 40 arc-seconds. The 1-sigma statistical sensitivity of the image is 0.017 mJy/beam or 0.12 mK_CMB at the 5 arc-second full width at half maximum. The SZE image shows a good agreement with an electron pressure map reconstructed independently from the X-ray data and offers a new probe of the small-scale structure of the intracluster medium. Our results demonstrate that ALMA is a powerful instrument for imaging the SZE in compact galaxy clusters with unprecedented angular resolution and sensitivity. As the first report on the detection of the SZE by ALMA, we present detailed analysis procedures including corrections for the missing flux, to provide guiding methods for analyzing and interpreting future SZE images by ALMA.

Either bulk rotation or local turbulence is widely invoked to drive fragmentation in collapsing cores so as to produce multiple star systems. Even when the two mechanisms predict different manners in which the stellar spins and orbits are aligned, subsequent internal or external interactions can drive multiple systems towards or away from alignment thus masking their formation process. Here, we demonstrate that the geometrical and dynamical relationship between the binary system and its surrounding bulk envelope provide the crucial distinction between fragmentation models. We find that the circumstellar disks of the binary protostellar system L1551 IRS 5 are closely parallel not just with each other but also with their surrounding flattened envelope. Measurements of the relative proper motion of the binary components spanning nearly 30 yr indicate an orbital motion in the same sense as the envelope rotation. Eliminating orbital solutions whereby the circumstellar disks would be tidally truncated to sizes smaller than are observed, the remaining solutions favor a circular or low-eccentricity orbit tilted by up to $\sim$25$^\circ$ from the circumstellar disks. Turbulence-driven fragmentation can generate local angular momentum to produce a coplanar binary system, but which bears no particular relationship with its surrounding envelope. Instead, the observed properties conform with predictions for rotationally-driven fragmentation. If the fragments were produced at different heights or on opposite sides of the midplane in the flattened central region of a rotating core, the resulting protostars would then exhibit circumstellar disks parallel with the surrounding envelope but tilted from the orbital plane as is observed.

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.

Ground-based observations at terahertz (THz) frequencies are a newly explorable area of astronomy for the next ten years. We discuss science cases for a first-generation 10-m class THz telescope, focusing on the Greenland Telescope as an example of such a facility. We propose science cases and provide quantitative estimates for each case. The largest advantage of ground-based THz telescopes is their higher angular resolution (~ 4 arcsec for a 10-m dish), as compared to space or airborne THz telescopes. Thus, high-resolution mapping is an important scientific argument. In particular, we can isolate zones of interest for Galactic and extragalactic star-forming regions. The THz windows are suitable for observations of high-excitation CO lines and [N II] 205 um lines, which are scientifically relevant tracers of star formation and stellar feedback. Those lines are the brightest lines in the THz windows, so that they are suitable for the initiation of ground-based THz observations. THz polarization of star-forming regions can also be explored since it traces the dust population contributing to the THz spectral peak. For survey-type observations, we focus on ``sub-THz'' extragalactic surveys, whose uniqueness is to detect galaxies at redshifts z ~ 1--2, where the dust emission per comoving volume is the largest in the history of the Universe. Finally we explore possibilities of flexible time scheduling, which enables us to monitor active galactic nuclei, and to target gamma-ray burst afterglows. For these objects, THz and submillimeter wavelength ranges have not yet been explored.

We performed the first spectral-line survey at 82--106 GHz and 335--355 GHz toward the outflow-shocked region, OMC 2-FIR 4, the outflow driving source, FIR 3, and the northern outflow lobe, FIR 3N. We detected 120 lines of 20 molecular species. The line profiles are found to be classifiable into two types: one is a single Gaussian component with a narrow ($<$ 3 km s$^{-1}$) width and another is two Gaussian components with narrow and wide ($>$ 3km s$^{-1}$) widths. The narrow components for the most of the lines are detected at all positions, suggesting that they trace the ambient dense gas. For CO, CS, HCN, and HCO$^{+}$, the wide components are detected at all positions, suggesting the outflow origin. The wide components of C$^{34}$S, SO, SiO, H$^{13}$CN, HC$^{15}$N, H$_2^{13}$CO, H$_2$CS, HC$_3$N, and CH$_3$OH are detected only at FIR 4, suggesting the outflow-shocked gas origin. The rotation diagram analysis revealed that the narrow components of C$_2$H and H$^{13}$CO$^+$ show low temperatures of 12.5$\pm$1.4 K, while the wide components show high temperatures of 20--70 K. This supports our interpretation that the wide components trace the outflow and/or outflow-shocked gas. We compared observed molecular abundances relative to H$^{13}$CO$^+$ with those of the outflow-shocked region, L1157 B1, and the hot corino, IRAS 16293-2422. Although we cannot exclude a possibility that the chemical enrichment in FIR 4 is caused by the hot core chemistry, the chemical compositions in FIR 4 are more similar to those in L1157 B1 than those in IRAS 16293-2422.

We performed mapping observations of the Class I protostellar binary system L1551 NE in the C$^{18}$O ($J$=3-2), $^{13}$CO ($J$=3-2), CS ($J$=7-6), and SO ($J_N$=7$_8$-6$_7$) lines with Atacama Submillimeter Telescope Experiment (ASTE). The ASTE C$^{18}$O data are combined with our previous SMA C$^{18}$O data, which show a $r \sim$300-AU scale Keplerian disk around the protostellar binary system. The C$^{18}$O maps show a $\sim$20000-AU scale protostellar envelope surrounding the central Keplerian circumbinary disk. The envelope exhibits a northeast (blue) - southwest (red) velocity gradient along the minor axis, which can be interpreted as a dispersing gas motion with an outward velocity of 0.3 km s$^{-1}$, while no rotational motion in the envelope is seen. In addition to the envelope, two $\lesssim$4000 AU scale, high-velocity ($\gtrsim$1.3 km s$^{-1}$) redshifted $^{13}$CO and CS emission components are found to $\sim$40$^{\prime\prime}$ southwest and $\sim$20$^{\prime\prime}$ west of the protostellar binary. These redshifted components are most likely outflow components driven from the neighboring protostellar source L1551 IRS 5, and are colliding with the envelope in L1551 NE. The net momentum, kinetic and internal energies of the L1551 IRS 5 outflow components are comparable to those of the L1551 NE envelope, and the interactions between the outflows and the envelope are likely to cause the dissipation of the envelope and thus suppression of the further growth of the mass and mass ratio of the central protostellar binary in L1551 NE.

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 present Atacama Large Millimeter/submillimeter Array (ALMA) observations from the 2014 Long Baseline Campaign in dust continuum and spectral line emission from the HL Tau region. The continuum images at wavelengths of 2.9, 1.3, and 0.87 mm have unprecedented angular resolutions of 0.075 arcseconds (10 AU) to 0.025 arcseconds (3.5 AU), revealing an astonishing level of detail in the circumstellar disk surrounding the young solar analogue HL Tau, with a pattern of bright and dark rings observed at all wavelengths. By fitting ellipses to the most distinct rings, we measure precise values for the disk inclination (46.72pm0.05 degrees) and position angle (+138.02pm0.07 degrees). We obtain a high-fidelity image of the 1.0 mm spectral index ($\alpha$), which ranges from $\alpha\sim2.0$ in the optically-thick central peak and two brightest rings, increasing to 2.3-3.0 in the dark rings. The dark rings are not devoid of emission, we estimate a grain emissivity index of 0.8 for the innermost dark ring and lower for subsequent dark rings, consistent with some degree of grain growth and evolution. Additional clues that the rings arise from planet formation include an increase in their central offsets with radius and the presence of numerous orbital resonances. At a resolution of 35 AU, we resolve the molecular component of the disk in HCO+ (1-0) which exhibits a pattern over LSR velocities from 2-12 km/s consistent with Keplerian motion around a ~1.3 solar mass star, although complicated by absorption at low blue-shifted velocities. We also serendipitously detect and resolve the nearby protostars XZ Tau (A/B) and LkHa358 at 2.9 mm.

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.

We present the results of continuum and 12CO(3-2) and CH3OH(7-6) line observations of IRAS16547-4247 made with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of ~0.5". The 12CO(3-2) emission shows two high-velocity outflows whose driving sources are located within the dust continuum peak. The alignment of these outflows do not coincide with that of the wide-angle, large scale, bipolar outflow detected with APEX in previous studies. The CH3OH(7-6) line emission traces an hourglass structure associated with the cavity walls created by the outflow lobes. Taking into account our results together with the position of the H2O and class I CH3OH maser clusters, we discuss two possible scenarios that can explain the hourglass structure observed in IRAS16547-4247: (1) precession of a biconical jet, (2) multiple, or at least two, driving sources powering intersecting outflows. Combining the available evidence, namely, the presence of two cross-aligned bipolar outflows and two different H2O maser groups, we suggest that IRAS16547-4247 represents an early formation phase of a protocluster.

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 report the ALMA observation of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), and 13CO (3-2) lines at a ~1.6 times higher resolution and a ~6 times higher sensitivity than those of our previous SMA observations, which revealed a r ~300 AU-scale circumbinary disk in Keplerian rotation. The 0.9-mm continuum shows two opposing U-shaped brightenings in the circumbinary disk, and exhibits a depression between the circumbinary disk and the circumstellar disk of the primary protostar. The molecular lines trace non-axisymmetric deviations from Keplerian rotation in the circumbinary disk at higher velocities relative to the systemic velocity, where our previous SMA observations could not detect the lines. In addition, we detect inward motion along the minor axis of the circumbinary disk. To explain the newly-observed features, we performed a numerical simulation of gas orbits in a Roche potential tailored to the inferred properties of L1551 NE. The observed U-shaped dust features coincide with locations where gravitational torques from the central binary system are predicted to impart angular momentum to the circumbinary disk, producing shocks and hence density enhancements seen as a pair of spiral arms. The observed inward gas motion coincides with locations where angular momentum is predicted to be lowered by the gravitational torques. The good agreement between our observation and model indicates that gravitational torques from the binary stars constitute the primary driver for exchanging angular momentum so as to permit infall through the circumbinary disk of L1551 NE.

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.

After more than 30 years of investigations, the nature of gas-grain interactions at low temperatures remains an unresolved issue in astrochemistry. Water ice is the dominant ice found in cold molecular clouds, however, there is only one region where cold (~10 K) water vapor has been detected - L1544. This study aims to shed light on ice desorption mechanisms under cold cloud conditions by expanding the sample. The clumpy distribution of methanol in dark clouds testifies to transient desorption processes at work -- likely to also disrupt water ice mantles. Therefore, the Herschel HIFI instrument was used to search for cold water in a small sample of prominent methanol emission peaks. We report detections of the ground-state transition of o-H2O (J = 1_10 - 1_01) at 556.9360 GHz toward two positions in the cold molecular cloud Barnard 5. The relative abundances of methanol and water gas support a desorption mechanism which disrupts the outer ice mantle layers, rather than causing complete mantle removal.

We present SMA interferometric observations of the 12CO J=2-1, 13CO (J=2-1), and C18O (J=2-1) lines and 225 GHz continuum emission and SMT single-dish observations of C18O (J=2-1) toward L43, a protostellar object in transition from Class I to II. The 225 GHz continuum emission shows a weak (~23.6 mJy), compact (<1000 AU) component associated with the central protostar. Our simulated observations show that it can be explained by dust thermal emission arising from an envelope which has a hole or a constant intensity region within a few hundred AU of the protostar. This suggests the disappearance or a lower concentration distribution of the envelope on a small scale. The 12CO and 13CO emission exhibit molecular outflows to the south and north. The C18O emission shows two molecular blobs, which correspond to the reflection nebulosity seen in near-infrared images, while there is no C18O emission associated with the protostar. The near-infrared features are likely due to the scattering at the positions of the blobs. The visible scattering features should result from the optical thinness of the envelope material, which is consistent with the less-concentrated distribution in the continuum emission. From single-dish observations, we found that the mass of the envelope (~1.5 Msun) + protostar (~0.5 Msun) is comparable with the virial mass of M_vir=1.0 Msun within 40 arcsec. This suggests that the envelope is likely gravitationally bound. We suggest that the protostellar envelope of L43 has been disappearing by consumption through accretion, at least in the close vicinity of the protostar.

We present the various science cases for building Band 1 receivers as part of ALMA's ongoing Development Program. We describe the new frequency range for Band 1 of 35-52 GHz, a range chosen to maximize the receiver suite's scientific impact. We first describe two key science drivers: 1) the evolution of grains in protoplanetary disks and debris disks, and 2) molecular gas in galaxies during the era of re-ionization. Studies of these topics with Band 1 receivers will significantly expand ALMA's Level 1 Science Goals. In addition, we describe a host of other exciting continuum and line science cases that require ALMA's high sensitivity and angular resolution. For example, ALMA Band 1 continuum data will probe the Sunyaev-Zel'dovich Effect in galaxy clusters, Very Small Grains and spinning dust, ionized jets from young stars, spatial and flaring studies of Sgr A*, the acceleration sites of solar flares, pulsar wind nebulae, radio supernovae, and X-ray binaries. Furthermore, ALMA Band 1 line data will probe chemical differentiation in cloud cores, complex carbon chain molecules, extragalactic radio recombination lines, masers, magnetic fields through Zeeman effect measurements, molecular outflows from young stars, the co-evolution of star formation and active galactic nuclei, and the molecular content of galaxies at z ~ 3. ALMA provides similar to better sensitivities than the JVLA over 35-50 GHz, with differences increasing with frequency. ALMA's smaller antennas and shorter baselines, greater number of baselines, and single-dish capabilities, however, give it a significant edge for observing extended emission, making wide-field maps (mosaics), or attaining high image fidelity, as required by the described science cases.

We report follow-up observations of the Class I binary protostellar system L1551 NE in the C18O (3--2) line with the SMA in its compact and subcompact configurations. Our previous observations at a higher angular resolution in the extended configuration revealed a circumbinary disk exhibiting Keplerian motion. The combined data having more extensive spatial coverage (~140 - 2000 AU) verify the presence of a Keplerian circumbinary disk, and reveals for the first time a distinct low-velocity (~< +-0.5 km s-1 from the systemic velocity) component that displays a velocity gradient along the minor axis of the circumbinary disk. Our simple model that reproduces the main features seen in the Position-Velocity diagrams comprises a circumbinary disk exhibiting Keplerian motion out to a radius of ~300 AU, beyond which the gas exhibits pure infall at a constant velocity of ~0.6 km s-1. The latter is significantly smaller than the expected free-fall velocity of ~2.2 km s-1 onto the L1551 NE protostellar mass of ~0.8 Msun at ~300 AU, suggesting that the infalling gas is decelerated as it moves into regions of high gas pressure in the circumbinary disk. The discontinuity in angular momenta between the outer infalling gas and inner Keplerian circumbinary disk implies an abrupt transition in the effectiveness at which magnetic braking is able to transfer angular momentum outwards, a result perhaps of the different plasma beta and ionization fractions between the outer and inner regions of the circumbinary disk.

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 present a set of compelling science cases for the ALMA Band 1 receiver suite. For these cases, we assume in tandem the updated nominal Band 1 frequency range of 35-50 GHz with a likely extension up to 52 GHz; together these frequencies optimize the Band 1 science return. The scope of the science cases ranges from nearby stars to the re-ionization edge of the Universe. Two cases provide additional leverage on the present ALMA Level One Science Goals and are seen as particularly powerful motivations for building the Band 1 Receiver suite: (1) detailing the evolution of grains in protoplanetary disks, as a complement to the gas kinematics, requires continuum observations out to ~35 GHz (~9mm); and (2) detecting CO 3-2 line emission from galaxies like the Milky Way during the epoch of re-ionization, i.e., 6 < z < 10, also requires Band 1 receiver coverage. The range of Band 1 science is wide, however, and includes studies of very small dust grains in the ISM, pulsar wind nebulae, radio supernovae, X-ray binaries, the Galactic Center (i.e., Sgr A*), dense cloud cores, complex carbon-chain molecules, masers, magnetic fields in the dense ISM, jets and outflows from young stars, distant galaxies, and galaxy clusters (i.e., the Sunyaev-Zel'dovich Effect). A comparison of ALMA and the Jansky VLA (JVLA) at the same frequencies of Band 1 finds similar sensitivity performance at 40-50 GHz, with a slight edge for ALMA at higher frequencies (e.g., within a factor of 2 for continuum observations). With its larger number of instantaneous baselines, however, ALMA Band 1data will have greater fidelity than those from the JVLA at similar frequencies.

We present combined Submillimeter-Array (SMA) + single-dish images of the (sub)millimeter dust continuum emission toward two prestellar cores SM1 and B2-N5 in the nearest star cluster forming region, $\rho$ Ophiuchus. Our combined images indicate that SM1 and B2-N5 consist of three and four condensations, respectively, with masses of $10^{-2}-10^{-1}M_\odot$ and sizes of a few hundred AU. The individual condensations have mean densities of $10^8-10^9$ cm$^{-3}$ and the masses are comparable to or larger than the critical Bonner-Ebert mass, indicating that the self-gravity plays an important role in the dynamical evolution of the condensations. The coalescence timescale of these condensations is estimated to be about $10^4$ yr, which is comparable to the local gravitational collapse timescale, suggesting that merging of the condensations, instead of accretion, plays an essential role in the star formation process. These results challenge the standard theory of star formation, where a single, rather featureless prestellar core collapses to form at most a couple of condensations, each of which potentially evolves into a protostar that is surrounded by a rotating disk where planets are created.

We present imaging simulations of the Sunyaev-Zel'dovich effect of galaxy clusters for the Atacama Large Millimeter/submillimeter Array (ALMA) including the Atacama Compact Array (ACA). In its most compact configuration at 90GHz, ALMA will resolve the intracluster medium with an effective angular resolution of 5 arcsec. It will provide a unique probe of shock fronts and relativistic electrons produced during cluster mergers at high redshifts, that are hard to spatially resolve by current and near-future X-ray detectors. Quality of image reconstruction is poor with the 12m array alone but improved significantly by adding ACA; expected sensitivity of the 12m array based on the thermal noise is not valid for the Sunyaev-Zel'dovich effect mapping unless accompanied by an ACA observation of at least equal duration. The observations above 100 GHz will become excessively time-consuming owing to the narrower beam size and the higher system temperature. On the other hand, significant improvement of the observing efficiency is expected once Band 1 is implemented in the future.

We present SubMillimeter-Array observations of a Keplerian disk around the Class I protobinary system L1551 NE in 335 GHz continuum emission and submillimeter line emission in 13CO (J=3-2) and C18O (J=3-2) at a resolution of ~120 x 80 AU. The 335-GHz dust-continuum image shows a strong central peak closely coincident with the binary protostars and likely corresponding to circumstellar disks, surrounded by a ~600 x 300 AU feature elongated approximately perpendicular to the [Fe II] jet from the southern protostellar component suggestive of a circumbinary disk. The 13CO and C18O images confirm that the circumbinary continuum feature is indeed a rotating disk; furthermore, the C18O channel maps can be well modeled by a geometrically-thin disk exhibiting Keplerian rotation. We estimate a mass for the circumbinary disk of ~0.03-0.12 Msun, compared with an enclosed mass of ~0.8 Msun that is dominated by the protobinary system. Compared with several other Class I protostars known to exhibit Keplerian disks, L1551 NE has the lowest bolometric temperature (~91 K), highest envelope mass (~0.39 Msun), and the lowest ratio in stellar mass to envelope + disk + stellar mass (~0.65). L1551 NE may therefore be the youngest protostellar object so far found to exhibit a Keplerian disk. Our observations present firm evidence that Keplerian disks around binary protostellar systems, ``Keplerian circumbinary disks', can exist. We speculate that tidal effects from binary companions could transport angular momenta toward the inner edge of the circumbinary disk and create the Keplerian circumbinary disk.

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.

It remains a key challenge to establish the molecular content of different components of low-mass protostars, like their envelopes and disks, and how this depends on the evolutionary stage and/or environment of the young stars. Observations at submillimeter wavelengths provide a direct possibility to study the chemical composition of low-mass protostars through transitions probing temperatures up to a few hundred K in the gas surrounding these sources. This paper presents a large molecular line survey of the deeply embedded protostellar binary IRAS 16293-2422 from the Submillimeter Array (SMA) - including images of individual lines down to approximately 1.5-3" (190-380 AU) resolution. More than 500 individual transitions are identified related to 54 molecular species (including isotopologues) probing temperatures up to about 550 K. Strong chemical differences are found between the two components in the protostellar system with a separation between, in particular, the sulfur- and nitrogen-bearing species and oxygen-bearing complex organics. The action of protostellar outflow on the ambient envelope material is seen in images of CO and SiO and appear to influence a number of other species, including (deuterated) water, HDO. The effects of cold gas-phase chemistry is directly imaged through maps of CO, N2D+ and DCO+, showing enhancements of first DCO+ and subsequently N2D+ in the outer envelope where CO freezes-out on dust grains.

We have made mapping observations of L1551 IRS 5, L1551NE, L723, and L43 and single-point observations of IRAS 16293-2422 in the submillimeter CS (J = 7-6) and HCN (J = 4-3) lines with ASTE. Including our previous ASTE observations of L483 and B335, we found a clear linear correlation between the source bolometric luminosities and the total integrated intensities of the submillimeter lines (I_CS ~L_bol^0.92). The combined ASTE + SMA CS (7-6) image of L1551 IRS 5 exhibits an extended (~2000 AU) component tracing the associated reflection nebula at the west and southwest, as well as a compact (< 500 AU) component centered on the protostellar position. The emission peaks of the CS and HCN emissions in L1551 NE are not located at the protostellar position but offset (~1400 AU) toward the associated reflection nebula at the west. With the statistical analyses, we confirmed the opposite velocity gradients of the CS (7-6) emission to those of the millimeter lines along the outflow direction, which we reported in our early paper. The magnitudes of the submillimeter velocity gradients are estimated to be (9.7\pm1.7) \times 10-3 km s-1 arcsec-1 in L1551 IRS 5 and (7.6\pm2.4) \times 10-3 km s-1 arcsec-1 in L483. We suggest that the "skewed" submillimeter molecular emissions toward the associated reflection nebulae at a few thousands AU scale trace the warm (> 40 K) walls of the envelope cavities, excavated by the associated outflows and irradiated by the central protostars directly. The opposite velocity gradients along the outflow direction likely reflect the dispersing gas motion at the wall of the cavity in the envelopes perpendicular to the outflow.

We have observed dense gas around the Very Low-Luminosity Ob jects (VeLLOs) L1521F-IRS and IRAM 04191+1522 in carbon-chain and organic molecular lines with the Nobeyama 45 m telescope. Towards L1521F-IRS, carbon-chain lines of CH3CCH (50-40), C4H (17/2-15/2), and C3H2 (212-101) are 1.5 - 3.5 times stronger than those towards IRAM 04191+1522, and the abundances of the carbon-chain molecules towards L1521F-IRS are 2 to 5 times higher than those towards IRAM 04191+1522. Mapping observations of these carbon-chain molecular lines show that in L1521F the peak positions of these carbon-chain molecular lines are different from each other and there is no emission peak towards the VeLLO position, while in IRAM 04191+1522 these carbon-chain lines are as weak as the detection limits except for the C3H2 line. The observed chemical differentiation between L1521F and IRAM 04191+1522 suggests that the evolutionary stage of L1521F-IRS is younger than that of IRAM 04191+1522, consistent with the extent of the associated outflows seen in the 13CO (1-0) line. The non-detection of the organic molecular lines of CH3OH (6-2-7-1 E) and CH3CN (60-50) implies that the warm (~ 100 K) molecular-desorbing region heated by the central protostar is smaller than ~ 100 AU towards L1521F-IRS and IRAM 04191+1522, suggesting the young age of these VeLLOs. We propose that the chemical status of surrounding dense gas can be used to trace the evolutionary stages of VeLLOs.

We present new, wide and deep images in the 1.1 mm continuum and the $^{12}$CO ($J$=1-0) emission toward the northern part of the Orion A Giant Molecular Cloud (Orion-A GMC). The 1.1 mm data were taken with the AzTEC camera mounted on the Atacama Submillimeter Telescope Experiment (ASTE) 10 m telescope in Chile, and the $^{12}$CO ($J$=1-0) data were with the 25 beam receiver (BEARS) on the NRO 45 m telescope in the On-The-Fly (OTF) mode. The present AzTEC observations are the widest $(\timeform{1.D7}$ $\times$ $\timeform{2.D3}$, corresponding to 12 pc $\times$ 17 pc) and the highest-sensitivity ($\sim$9 mJy beam$^{-1}$) 1.1 mm dust-continuum imaging of the Orion-A GMC with an effective spatial resolution of $\sim$ 40$\arcsec$. The $^{12}$CO ($J$=1-0) image was taken over the northern $\timeform{1D.2} \times\timeform{1D.2}$ (corresponding 9 pc $\times$ 9 pc) area with a sensitivity of 0.93 K in $T_{\rm MB}$, a velocity resolution of 1.0 km s$^{-1}$, and an effective spatial resolution of 21$\arcsec$. With these data, together with the MSX 8 $\mu$m, Spitzer 24 $\mu$m and the 2MASS data, we have investigated the detailed structure and kinematics of molecular gas associated with the Orion-A GMC and have found evidence for interactions between molecular clouds and the external forces that may trigger star formation. Two types of possible triggers were revealed; 1) Collision of the diffuse gas on the cloud surface, particularly at the eastern side of the OMC-2/3 region, and 2) Irradiation of UV on the pre-existing filaments and dense molecular cloud cores. Our wide-field and high-sensitivity imaging have provided the first comprehensive view of the potential sites of triggered star formation in the Orion-A GMC.

In order to study how outflows from protostars influence the physical and chemical conditions of the parent molecular cloud, we have observed Barnard 1 (B1) main core, which harbors four Class 0 and three Class I sources, in the CO (J=1-0), CH3OH (J_K=2_K-1_K), and the SiO (J=1-0) lines using the Nobeyama 45 m telescope. We have identified three CO outflows in this region; one is an elongated (~ 0.3 pc) bipolar outflow from a Class 0 protostar B1-c in the submillimeter clump SMM 2, another is a rather compact (~ 0.1 pc) outflow from a Class I protostar B1 IRS in the clump SMM 6, and the other is an extended outflow from a Class I protostar in SMM 11. In the western lobe of the SMM 2 outflow, both the SiO and CH3OH lines show broad redshifted wings with the terminal velocities of 25 km/s and 13 km/s, respectively. It is likely that the shocks caused by the interaction between the outflow and ambient gas enhance the abundance of SiO and CH3OH in the gas phase. The total energy input rate by the outflows (1.1x10^{-3} Lsun) is smaller than the energy loss rate (8.5x10^{-3} Lsun$) through the turbulence decay in B1 main core, which suggests that the outflows can not sustain the turbulence in this region. Since the outflows are energetic enough to compensate the dissipating turbulence energy in the neighboring, more evolved star forming region NGC 1333, we suggest that the turbulence energy balance depends on the evolutionary state of the star formation in molecular clouds.

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.

We present a Nobeyama 45 m Radio Telescope map and Australia Telescope Compact Array pointed observations of N2H+ 1-0 emission towards the clustered, low mass star forming Oph B Core within the Ophiuchus molecular cloud. We compare these data with previously published results of high resolution NH3 (1,1) and (2,2) observations in Oph B. We use 3D Clumpfind to identify emission features in the single-dish N2H+ map, and find that the N2H+ `clumps' match well similar features previously identified in NH3 (1,1) emission, but are frequently offset to clumps identified at similar resolution in 850 micron continuum emission. Wide line widths in the Oph B2 sub-Core indicate non-thermal motions dominate the Core kinematics, and remain transonic at densities n ~ 3 x 10^5 cm^-3 with large scatter and no trend with N(H2). Non-thermal motions in Oph B1 and B3 are subsonic with little variation, but also show no trend with H2 column density. Over all Oph B, non-thermal N2H+ line widths are substantially narrower than those traced by NH3, making it unlikely NH3 and N2H+ trace the same material, but the v_LSR of both species agree well. We find evidence for accretion in Oph B1 from the surrounding ambient gas. The NH3/N2H+ abundance ratio is larger towards starless Oph B1 than towards protostellar Oph B2, similar to recent observational results in other star-forming regions. Small-scale structure is found in the ATCA N2H+ 1-0 emission, where emission peaks are again offset from continuum emission. In particular, the ~1 M_Sun B2-MM8 clump is associated with a N2H+ emission minimum and surrounded by a broken ring-like N2H+ emission structure, suggestive of N2H+ depletion. We find a strong general trend of decreasing N2H+ abundance with increasing N(H2) in Oph B which matches that found for NH3.

We have conducted millimeter interferometric observations of the Orion Molecular Cloud-2 (OMC-2) FIR 6 region at an angular resolution of $\sim$ 4$\arcsec$ - 7$\arcsec$ with the Nobeyama Millimeter Array (NMA). In the 3.3 mm continuum emission we detected dusty core counterparts of the previously identified FIR sources (FIR 6a, 6b, 6c, and 6d), and moreover, resolved FIR 6a into three dusty cores. The size and mass of these cores are estimated to be 1100-5900 AU and 0.19-5.5 M$_{\odot}$, respectively. We found that in the $^{12}$CO ($J$=1--0) emission FIR 6b, 6c, and 6d eject the molecular outflow and that the FIR 6c outflow also exhibits at least two collimated jet-like components in the SiO ($J$=2--1) emission. At the tip of one of the SiO components there appears abrupt increase of the SiO line width ($\sim$ 15 km s$^{-1}$), where the three resolved cores in FIR 6a seem to delineate the tip. These results imply the presence of the interaction and the bowshock front between the FIR 6c molecular outflow and FIR 6a.

The multiple protostellar system L1551 IRS5 exhibits a large-scale bipolar molecular outflow. We have studied this outflow within ~4000 AU of its driving source(s) with the SubMillimeter Array. Our CO(2-1) image at ~4" (~560 AU) resolution reveals three distinct components: 1) an X-shaped structure spanning ~20" from center with a similar symmetry axis and velocity pattern as the large-scale outflow; 2) an S-shaped structure spanning ~10" from center also with an opposite velocity pattern to the large-scale outflow; and 3) a compact central component spanning ~1.4" from center again with a similar symmetry axis and velocity pattern as the large-scale outflow. The X-shaped component likely comprises the limb-brightened walls of a cone-shaped cavity excavated by the outflows from the two main protostellar components. The compact central component likely comprises material newly entrained by one or both outflows from the two main protostellar components. The S-shaped component mostly likely comprises a precessing outflow with its symmetry axis inclined in the opposite sense to the plane of the sky than the other two components. This outflow may be driven by a recently reported candidate third protostellar component in L1551 IRS5, whose circumstellar disk is misaligned relative to the two main protostellar components. Gravitational interactions between this protostellar component and its more massive northern neighbor may be causing the circumstellar disk and hence outflow of this component to precess.

Using the ASTE 10 m submillimeter telescope and the 1.4 m Infrared Survey Facility (IRSF), we performed an extensive outflow survey in the Orion Molecular Cloud -2 and -3 region. Our survey, which includes 41 potential star-forming sites, has been newly compiled using multi-wavelength data based on millimeter- and submillimeter-continuum observations as well as radio continuum observations. From the CO (3-2) observations performed with the ASTE 10 m telescope, we detected 14 CO molecular outflows, seven of which were newly identified. This higher detection rate, as compared to previous CO (1-0) results in the same region, suggests that CO (3-2) may be a better outflow tracer. Physical properties of these outflows and their possible driving sources were derived. Derived parameters were compared with those of CO outflows in low- and high-mass starforming regions. We show that the CO outflow momentum correlates with the bolometric luminosity of the driving source and with the envelope mass, regardless of the mass of the driving sources. In addition to these CO outflows, seven sources having NIR features suggestive of outflows were also identified.

We have carried out the observations of the OMC-2 FIR 3/4 region with the NMA and ASTE in the H$^{13}$CO$^{+}$ (1--0), $^{12}$CO (3--2, 1--0), SiO ($v$=0, $J$=2--1), CS (2--1), and CH$_3$OH ($J_K$=7$_K$--6$_K$) lines and in the 3.3 mm continuum emission. Our NMA observations in the H$^{13}$CO$^{+}$ emission have revealed 0.07 pc-scale dense gas associated with FIR 4. The $^{12}$CO (3--2,1--0) emission shows high-velocity blue and red shifted components at the both north-east and south-west of FIR 3, suggesting a molecular outflow nearly along the plane of the sky driven by FIR 3. The SiO and the CH$_{3}$OH emission are detected around the interface between the outflow and the dense gas. Furthermore, the $^{12}$CO (1--0) emission shows an L-shaped structure in the P-V diagram. These results imply presence of the shock due to the interaction between the molecular outflow driven by FIR 3 and the dense gas associated with FIR 4. Moreover, our high angular-resolution observations of FIR 4 in the 3.3 mm continuum emission have first found that FIR 4 consists of eleven dusty cores. The separation among these cores is on the same order of the Jeans length, suggesting that the fragmentation into these cores has been caused by the gravitational instability. The time scale of the fragmentation is similar to the time scale of the interaction between the molecular outflow and the dense gas. We suggest that the interaction between the molecular outflow from FIR 3 and the dense gas associated with FIR 4 triggered the fragmentation into these dusty cores, and hence the next generation the cluster formation.

We have mapped the proto-binary source IRAS 16293-2422 in CO 2-1, 13CO 2-1, and CO 3-2 with the Submillimeter Array (SMA). The maps with resolution of 1".5-5" reveal a single small scale (~3000 AU) bipolar molecular outflow along the east-west direction. We found that the blueshifted emission of this small scale outflow mainly extends to the east and the redshifted emission to the west from the position of IRAS 16293A. A comparison with the morphology of the large scale outflows previously observed by single-dish telescopes at millimeter wavelengths suggests that the small scale outflow may be the inner part of the large scale (~15000 AU) E-W outflow. On the other hand, there is no clear counterpart of the large scale NE-SW outflow in our SMA maps. Comparing analytical models to the data suggests that the morphology and kinematics of the small scale outflow can be explained by a wide-angle wind with an inclination angle of ~30-40 degrees with respect to the plane of the sky. The high resolution CO maps show that there are two compact, bright spots in the blueshifted velocity range. An LVG analysis shows that the one located 1" to the east of source A is extremely dense, n(H_2)~10^7 cm^-3, and warm, T_kin >55 K. The other one located 1" southeast of source B has a higher temperature of T_kin >65 K but slightly lower density of n(H_2)~10^6 cm^-3. It is likely that these bright spots are associated with the hot core-like emission observed toward IRAS 16293. Since both two bright spots are blueshifted from the systemic velocity and are offset from the protostellar positions, they are likely formed by shocks.

We have extensively mapped a sample of dense molecular clouds (L1512, TMC-1C, L1262, Per 7, L1389, L1251E) in lines of HC3N, CH3OH, SO and C^{18}O. We demonstrate that a high degree of chemical differentiation is present in all of the observed clouds. We analyse the molecular maps for each cloud, demonstrating a systematic chemical differentiation across the sample, which we relate to the evolutionary state of the cloud. We relate our observations to the cloud physical, kinematical and evolutionary properties, and also compare them to the predictions of simple chemical models. The implications of this work for understanding the origin of the clumpy structures and chemical differentiation observed in dense clouds are discussed.

(Abridged) With the SMA we have made high angular-resolution (~1" = 160 AU) observations of the protobinary system IRAS 16293-2422 in the J = 4-3 lines of HCN and HC^15N, and in the continuum at 354.5 GHz. The HCN (4-3) line was also observed using the JCMT to supply missing short spacing information. Submillimeter continuum emission is detected from the individual binary components with a separation of ~5". The HC^15N (4-3) emission has revealed a compact (~500 AU) flattened structure (P.A. = -16 degree) associated with Source A, and shows a velocity gradient along the projected minor axis, which can be interpreted as an infalling gas motion. Our HCN image including the short-spacing information shows an extended (~3000 AU) circumbinary envelope as well as the compact structure associated with Source A. A toy model consisting of a flattened structure with radial infall towards a 1 Msun central star reproduces the HCN/HC^15N position-velocity diagram along the minor axis of the HC^15N emission. In the extended envelope there is also a North-East (Blue) to South-West (Red) velocity gradient across the binary alignment, which is likely to reflect gas motion in the swept-up dense gas associated with the molecular outflow from Source A. At Source B, there is only a weak compact structure with much narrower line widths (~2 km/s) seen in the optically-thin HC^15N emission than that at Source A (>10 km/s), and there is no clearly defined bipolar molecular outflow associated with Source B. These results imply the different evolutionary stages between Source A and B in the common circumbinary envelope.

We present ~2'' resolution CO(3-2), HCO+(4-3) and 880micron continuum images of the luminous infrared galaxy NGC6240 obtained at the Submillimeter Array. We find that the spatially resolved CO(3-2), HCO+(4-3) and the 880micron emission peaks between the two nuclear components that are both known to harbor AGNs. Our Large Velocity Gradient (LVG) analysis performed on each velocity channel suggests that the peak of the molecular gas emission traced in our observations is warm (T = 20 - 100K), dense (nH2 = 10^(5.0 - 5.4) cm^-3) and moderately optically thin (tau = 0.2 - 2) in the central 1 kpc. We also find large column densities of ~10^(23) cm^(-2). Such extreme conditions are observed over ~300 km s^(-1) centered around the CO derived systemic velocity. The derived molecular gas mass from the CO(3-2) emission and a CO-to-H2 conversion factor commonly used for ULIRGs is (6.9 +/- 1.7) x 10^9 Msun, and this is consistent with the mass derived from previous CO(2-1) observations. The gas is highly turbulent in the central kpc (Delta v_(FWZI) ~ 1175 km s^(-1)). Furthermore, possible inflow or outflow activity is suggested from the CO(3-2) velocity distribution. We tentatively state that 3.5 x 10^8 Msun of isolated CO(3-2) emission seen west of the northern disk may be associated with outflows from starburst superwinds, but the gas outflow scenario from one of the central AGN is not completely ruled out. Piecing all of the information together, the central region of NGC 6240 harbors 2 AGNs, ~10^(10) Msun of molecular gas mass, 5 x 10^7 Msun of dust mass, and has possible evidence of inflow and outflow activity.

This paper presents a large spectral line and continuum survey of 8 deeply embedded, low-mass protostellar cores using the Submillimeter Array. Each source was observed in high excitation lines of some of the most common molecular species, CO, HCO+, CS, SO, H2CO, CH3OH and SiO. Line emission from 11 species originating from warm and dense gas have been imaged at high angular resolution (1-3"; typically 200-600 AU) together with continuum emission at 230 GHz (1.3 mm) and 345 GHz (0.8 mm). Compact continuum emission is observed for all sources which likely originates in marginally optically thick circumstellar disks, with typical lower limits to their masses of 0.1 M_sun (1-10% of the masses of their envelopes) and having a dust opacity law with beta approximately 1. Prominent outflows are present in CO 2-1 observations in all sources: the most diffuse outflows are found in the sources with the lowest ratios of disk-to-envelope mass, and it is suggested that these sources are in a phase where accretion of matter from the envelope has almost finished and the remainder of the envelope material is being dispersed by the outflows. Other characteristic dynamical signatures are found with inverse P Cygni profiles indicative of infalling motions seen in the 13CO 2-1 lines toward NGC1333-IRAS4A and -IRAS4B. Outflow-induced shocks are present on all scales in the protostellar environments and are most clearly traced by the emission of CH3OH in NGC1333-IRAS4A and -IRAS4B. These observations suggest that the emission of CH3OH and H2CO from these proposed "hot corinos" are related to the shocks caused by the protostellar outflows. Only one source, NGC1333-IRAS2A, has evidence for hot, compact CH3OH emission coincident with the embedded protostar.

With the ASTE telescope, we have made observations of three low-mass protostellar envelopes around L483, B335, and L723 in the submillimeter CS ($J$=7--6) and HCN ($J$=4--3) lines. We detected both the CS and HCN lines toward all the targets, and the typical CS intensity ($\sim$ 1.0 K in T$_{B}$) is twice higher than that of the HCN line. Mapping observations of L483 in these lines have shown that the submillimeter emissions in the low-mass protostellar envelope are resolved, exhibit a western extension from the central protostar, and that the deconvolved size is $\sim$ 5500 AU $\times$ 3700 AU (P.A. = 78$^{\circ}$) in the HCN emission. The extent of the submillimeter emissions in L483 implies the presence of higher-temperature ($\gtrsim$ 40 K) gas at 4000 AU away from the central protostar, which suggests that we need to take 2-dimensional radiative transfer models with a flattened disklike envelope and bipolar cavity into account to explain the temperature structure inside the low-mass protostellar envelope. The position-velocity diagrams of these submillimeter lines in L483 and B335 exhibit different velocity gradients from those found in the previous millimeter observations. In particular, along the axis of the associated molecular outflow the sense of the velocity gradient traced by the submillimeter lines is opposite to that of the millimeter observations or the associated molecular outflow, both in L483 and B335. We suggest that expanding gas motions at the surface of the flattened disklike envelope around the protostar, which is irradiated from the central star directly, are the origin of the observed submillimeter velocity structure.

We report the results of H13CO+(1-0), CO(1-0), and 3.3 mm dust continuum observations toward one of the strongest mm-wave sources in OMC-3, MMS 7, with the Nobeyama Millimeter Array (NMA) and the Nobeyama 45 m telescope. With the NMA, we detected centrally-condensed 3.3 mm dust-continuum emission which coincides with the MIR source and the free-free jet. Our combined H13CO+ observations have revealed a disk-like envelope. The size and the mass of the disk-like envelope are 0.15 times 0.11 pc and 5.1 - 9.1 M_sun, respectively. The combined map also shows that the outer portion of the disk-like envelope has a fan-shaped structure which delineates the rim of the CO(1-0) outflow observed with the NMA. The position-velocity (P-V) diagrams in the H13CO+ (1-0) emission show that the velocity field in the disk-like envelope is composed of a dispersing gas motion and a possible rigid-like rotation. The mass dispersing rate is estimated to be (3.4 - 6.0) times 10^-5 M_sun/yr, which implies that MMS 7 has an ability to disperse ~10 M_sun during the protostellar evolutional time of a few times 10^5 yr. The specific angular momentum of the possible rotation in the disk-like envelope is nearly two orders of magnitude larger than that in low-mass cores. The turn-over point of the power law of the angular momentum distribution in the disk-like envelope (< 0.007 pc), which is likely to be related to the outer radius of the central mass accretion, is similar to the size of the 3.3 mm dust condensation. The intermediate-mass protostar MMS 7 is in the last stage of the main accretion phase and that the substantial portion of the outer gas has already been dispersed, while the mass accretion may still be on-going at the innermost region traced by the dusty condensation.

Using the recently-commissioned Submillimeter Array (SMA), we have detected several complex organic molecules, including (CH3)2O, C2H5OH, C2H5CN, and tentatively CH2CDCN, toward the protostellar hot cores of IRAS 16293-2422 at arcsecond-resolution (~< 400 AU in radius). Vibrationally excited transitions of SO, SO2 and HCN with energy levels up to 1800 K were also observed. In addition to the other organic molecules (HC3N, CH2CO, CH3OH, CH2CHCN and HCOOCH3) previously reported by us (Kuan et al. 2004), these results clearly indicate the existence of a rich organic chemistry in low-mass `hot corinos'. From the observation of optically thin HC15N emission, we conclude I16293A is a rotating circumstellar disk lying along the north-south direction ~10 deg to the east and with an inclination ~30 deg to the sky. We suggest that the observed vibrational SO and SO2 emission may originate from shock waves near or in the circumstellar disks. Between the two cores, we find a strong anticorrelation in emission from C2H5OH and C2H5CN. The relative contribution of gas phase and grain-surface chemistries to the production of the observed complex molecules is discussed. We point out the shortcomings underlying recent claims that all the O-bearing organics are formed on grains. The presence of so many well-known interstellar molecules in solar-type hot corinos strengthens the link between molecular cloud chemistry, the starting materials of protoplanetary disks such as the protosolar nebula, and the composition of comets. Establishing the fine details of this connection is crucial in answering fundamental questions concerning the importance of galactic astrochemistry for astrobiology.

Arcsecond-resolution spectral observations toward the protobinary system IRAS 16293-2422 at 344 and 354 GHz were conducted using the Submillimeter Array. Complex organic molecules such as CH3OH and HCOOCH3 were detected. Together with the rich organic inventory revealed, it clearly indicates the existence of two, rather than one, compact hot molecular cores (smaller than or equal to 400 AU in radius) associated with each of the protobinary components identified by their dust continuum emission in the inner star-forming core.

We have imaged the circumstellar envelope around the binary protostar L1551 IRS 5 in CS (J=7-6) and 343 GHz continuum emission at ~ 3 arcsec resolution using the Submillimeter Array. The continuum emission shows an elongated structure (~ 220 x 100 AU) around the binary perpendicular to the axis of the associated radio jet. The CS emission extends over ~ 400 AU, appears approximately circularly symmetric, and shows a velocity gradient from southeast (blueshifted) to northwest (redshifted). The direction of the velocity gradient is different from that observed in C18O (J=1-0). This may be because rotation is more dominant in the CS envelope than the C18O envelope, in which both infall and rotation exist. The CS emission may be divided into two velocity components: (1) a "high" velocity disk-like structure surrounding the protostar, +-1.0 - 1.5 km s-1 from the systemic velocity, and (2) a "low" velocity structure, located southwest of the protostar, < 1.0 km s -1 from the systemic velocity. The high-velocity component traces warm and dense gas with kinematics consistent with rotation around the protostar. The low-velocity component may arise from dense gas entrained in the outflow. Alternatively, this component may trace infalling and rotating gas in an envelope with a vertical structure.

We present ~2"-4" aperture synthesis observations of the circumstellar disk surrounding the nearby young star TW Hya in the CO J=2--1 and J=3--2 lines and associated dust continuum obtained with the partially completed Submillimeter Array. The extent and peak flux of the 230 and 345 GHz dust emission follow closely the predictions of the irradiated accretion disk model of Calvet et al. (2002). The resolved molecular line emission extends to a radius of at least 200 AU, the full extent of the disk visible in scattered light, and shows a clear pattern of Keplerian rotation. Comparison of the images with 2D Monte Carlo models constrains the disk inclination angle to 7+/-1 degrees. The CO emission is optically thick in both lines, and the kinetic temperature in the line formation region is ~20K. Substantial CO depletion, by an order of magnitude or more from canonical dark cloud values, is required to explain the characteristics of the line emission.

We have carried out mapping observations of the molecular core associated with the young Class 0 protostar, IRAM 04191+1522, in the CH3OH (JK=2K-1K) and C34S (J=2-1) lines using the 45 m telescope at Nobeyama Radio Observatory. As well as an elongated envelope associated with the protostellar formation (size \~0.07 pc x 0.04 pc, mass ~ 2.3 Mo), two compact (~ 0.03 pc) condensations were found in the CH3OH line at the southern edge of the elongated envelope, where the blueshifted CO outflow emerging from the protostar is located. In contrast to the elongated envelope, those compact CH3OH condensations show much larger line width (up to 2.0 km s-1) with centroid velocities blueshifted by ~ 0.8 km s-1. The compact condensations have momenta (~ 0.06 Mo km s-1) comparable to that of the blueshifted molecular outflow. These results suggest that the compact condensations are probably formed in the course of interaction between the outflow and the ambient gas surrounding the protostar, and that such interaction may cause dissipation of a part of the ambient gas. No drastic, localized enhancement of the CH3OH abundance is, however, observed toward the compact condensations, implying that there seems to be no significant shock heating at the compact condensations. This may be because the CO outflow velocity (< 10 km s-1) is too low to cause effective heating to release CH3OH on dust grains into gas phase.

We have investigated the evolution and distribution of molecules in collapsing prestellar cores via numerical chemical models, adopting the Larson-Penston solution and its delayed analogues to study collapse. Molecular abundances and distributions in a collapsing core are determined by the balance among the dynamical, chemical and adsorption time scales. When the central density n_H of a prestellar core with the Larson-Penston flow rises to 3 10^6 cm^{-3}, the CCS and CO column densities are calculated to show central holes of radius 7000 AU and 4000 AU, respectively, while the column density of N2H+ is centrally peaked. These predictions are consistent with observations of L1544. If the dynamical time scale of the core is larger than that of the Larson-Penston solution owing to magnetic fields, rotation, or turbulence, the column densities of CO and CCS are smaller, and their holes are larger than in the Larson-Penston core with the same central gas density. On the other hand, N2H+ and NH3 are more abundant in the more slowly collapsing core. Therefore, molecular distributions can probe the collapse time scale of prestellar cores. Deuterium fractionation has also been studied via numerical calculations. The deuterium fraction in molecules increases as a core evolves and molecular depletion onto grains proceeds. When the central density of the core is n_H=3 10^6 cm^{-3}, the ratio DCO+/HCO+ at the center is in the range 0.06-0.27, depending on the collapse time scale and adsorption energy; this range is in reasonable agreement with the observed value in L1544.