We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks.

We have studied the structure and kinematics of the dense molecular gas in the Orion Molecular Cloud 1 (OMC1) region with the N2H+ 3-2 line. The 6'x9' (~0.7x1.1 pc) region surrounding the Orion KL core has been mapped with the Submillimeter Array (SMA) and the Submillimeter Telescope (SMT). The combined SMA and SMT image having a resolution of ~5.4" (~2300 au) reveals multiple filaments with a typical width of 0.02-0.03 pc. On the basis of the non-LTE analysis using the N2H+ 3-2 and 1-0 data, the density and temperature of the filaments are estimated to be ~10^7 cm-3 and ~15-20 K, respectively. The core fragmentation is observed in three massive filaments, one of which shows the oscillations in the velocity and intensity that could be the signature of core-forming gas motions. The gas kinetic temperature is significantly enhanced in the eastern part of OMC1, likely due to the external heating from the high mass stars in M42 and M43. In addition, the filaments are colder than their surrounding regions, suggesting the shielding from the external heating due to the dense gas in the filaments. The OMC1 region consists of three sub-regions, i.e. north, west, and south of Orion KL, having different radial velocities with sharp velocity transitions. There is a north-to-south velocity gradient from the western to the southern regions. The observed velocity pattern suggests that dense gas in OMC1 is collapsing globally toward the high-mass star-forming region, Orion Nebula 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 ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature.

Episodic accretion has been used to explain the wide range of protostellar luminosities, but its origin and influence on the star forming process are not yet fully understood. We present an ALMA survey of N$_2$H$^+$ ($1-0$) and HCO$^+$ ($3-2$) toward 39 Class 0 and Class I sources in the Perseus molecular cloud. N$_2$H$^+$ and HCO$^+$ are destroyed via gas-phase reactions with CO and H$_2$O, respectively, thus tracing the CO and H$_2$O snowline locations. A snowline location at a much larger radius than that expected from the current luminosity suggests that an accretion burst has occurred in the past which has shifted the snowline outward. We identified 18/18 Class 0 and 9/10 Class I post-burst sources from N$_2$H$^+$, and 7/17 Class 0 and 1/8 Class I post-burst sources from HCO$^+$.The accretion luminosities during the past bursts are found to be $\sim10-100~L_\odot$. This result can be interpreted as either evolution of burst frequency or disk evolution. In the former case, assuming that refreeze-out timescales are 1000 yr for \ce{H2O} and 10,000 yr for CO, we found that the intervals between bursts increases from 2400 yr in the Class 0 to 8000 yr in the Class I stage. This decrease in the burst frequency may reflect that fragmentation is more likely to occur at an earlier evolutionary stage when the young stellar object is more prone to instability.

We present molecular line and dust continuum observations of a Planck-detected cold cloud, G074.11+00.11. The cloud consists of a system of curved filaments and a central star-forming clump. The clump is associated with several infrared sources and H2O maser emission. We aim to determine the mass distribution and gas dynamics within the clump, to investigate if the filamentary structure seen around the clump repeats itself on a smaller scale, and to estimate the fractions of mass contained in dense cores and filaments. The velocity distribution of pristine dense gas can be used to investigate the global dynamical state of the clump, the role of filamentary inflows, filament fragmentation and core accretion. We use molecular line and continuum observations from single dish observatories and interferometric facilities to study the kinematics of the region. The molecular line observations show that the central clump may have formed as a result of a large-scale filament collision. The central clump contains three compact cores. Assuming a distance of 2.3 kpc, based on Gaia observations and a three-dimensional extinction method of background stars, the mass of the central clump exceeds 700 solar masses, which is roughly 25% of the total mass of the cloud. Our virial analysis suggests that the central clump and all identified substructures are collapsing. We find no evidence for small-scale filaments associated with the cores. Our observations indicate that the clump is fragmented into three cores with masses in the range [10,50] solar masses and that all three are collapsing. The presence of an H2O maser emission suggests active star formation. However, the CO lines show only weak signs of outflows. We suggest that the region is young and any processes leading to star formation have just recently begun.

This white paper gives a brief summary of Galactic continuum surveys with the next generation 850 micron camera at the James Clerk Maxwell Telescope (JCMT) in the next decade. This new camera will have mapping speeds at least 10-20 times faster than the present SCUBA-2 camera, and will enable deep (<10 mJy/beam) and extensive continuum surveys such as a wider (-5 deg.<l<240 deg., |b|<2-5 deg.) Galactic Plane survey, a larger Gould Belt survey, and a follow-up survey of all Planck compact objects visible from the northern hemisphere. These surveys will provide a complete census of molecular clouds, filaments, and dense cores across the Galaxy, vital for studying star formation in various environments.

We present (sub)millimeter imaging at 0.5'' resolution of the massive star-forming region G358.93-0.03 acquired in multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz methanol maser emission. Using SMA and ALMA, we have discovered 14 new Class II methanol maser lines ranging in frequency from 199 GHz to 361 GHz, which originate mostly from vt=1 torsionally-excited transitions and include one vt=2 transition. The latter detection provides the first observational evidence that Class II maser pumping involves levels in the vt=2 state. The masers are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure 1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt=2 line showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1 has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core (MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of the sources during the one month interval, and the total protocluster emission remains comparable to prior single dish measurements.

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.

The HH 211 protostellar system is currently the youngest Class 0 system found with a rotating disk. We have mapped it at ~ 50 au (0.16") resolution, studying its magnetic field morphology with dust polarization in continuum at 232 and 358 GHz and its kinematics in C18O J=2-1 line. A flattened envelope extending out to ~ 400 au from the disk is detected in the continuum and C18O, slightly misaligned with the disk by 8 degree. It is spiraling inwards and expected to transform into a rotating disk at ~ 20 au, consistent with the disk radius estimated before. It appears to have a constant specific angular momentum and itself can result from an inside-out collapse of an extended envelope detected before in NH$_3$. In the flattened envelope, the polarization is mainly due to the magnetically aligned dust grains, inferring a highly pinched poloidal field morphology there. Thus, both the kinematics and field morphology support that the flattened envelope is a pseudodisk formed as the infalling gas is guided by the field lines to the equatorial plane. Interestingly, a point symmetric polarization distribution is also seen in the flattened envelope, implying that the pinched field lines also have a significant toroidal component generated by the rotation. No significant loss of angular momentum and thus no clear magnetic braking are detected in the flattened envelope around the disk probably because of the large misalignment between the axis of the rotation and the axis of the magnetic field in the cloud core.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

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 present ALMA long-baseline observations toward the Class 0 protostar IRAS 16253-2429 (hereafter IRAS 16253) with a resolution down to 0.12" (~15 au). The 1.3 mm dust continuum emission has a deconvolved Gaussian size of 0.16" x 0. 07" (20 au x 8.8 au), likely tracing an inclined dusty disk. Interestingly, the position of the 1.38 mm emission is offset from that of the 0.87 mm emission along the disk minor axis. Such an offset may come from a torus-like disk with very different optical depths between these two wavelengths. Furthermore, through CO (2 - 1) and C18O (2 - 1) observations, we study rotation and infall motions in this disk-envelope system and infer the presence of a Keplerian disk with a radius of 8 - 32 au. This result suggests that the disk could have formed by directly evolving from a first core, because IRAS16253 is too young to gradually grow a disk to such a size considering the low rotation rate of its envelope. In addition, we find a quadruple pattern in the CO emission at low velocity, which may originate from CO freeze out at the disk/envelope midplane. This suggests that the "cold disk" may appear in the early stage, implying a chemical evolution for the disk around this proto-brown dwarf (or very low-mass protostar) different from that of low-mass stars.

The filamentary infrared dark cloud L1287 is actively forming a dense cluster of low-mass YSOs at its inner $\sim$0.1 pc region. To help understand the origin of this low-mass YSO cluster, the present work aims at resolving the gas structures and kinematics. We have performed $\sim$1$"$ angular resolution ($\sim$930 AU) SMA observations at $\sim$1.3 mm. From a $\sim$2$"$ resolution 1.3 mm continuum image we identified six dense cores, namely SMA1-6 with masses in the range of $\sim0.4-4$ M$_\odot$. From a $\sim$1$"$ resolution 1.3 mm continuum image, we find a high fragmentation level, with 14 compact millimeter sources within 0.1 pc (two of them associated with the known accretion outburst YSOs RNO 1C and RNO 1B). The dense gas tracer DCN (3--2) traces well the dust continuum emission and shows a clear velocity gradient along the NW-SE direction centered at SMA3. There is another velocity gradient with opposite direction around the most luminous YSO IRAS 00338+6312. The fragmentation within 0.1 pc in L1287 is very high compared to other regions at the same spatial scales. The incoherent motions of dense gas flows are sometimes interpreted by being influenced by (proto)stellar feedback (e.g., outflows), which is not yet ruled out in this particular target source. The directions of the velocity gradients traced by DCN are approximately perpendicular to those of the dominant CO outflow(s). Therefore, we alternatively hypothesize that the velocity gradients revealed by DCN trace the convergence from the $\gtrsim$0.1 pc scales infalling motion towards the rotational motions around the more compact ($\sim0.02$ pc) sources. This global molecular gas converging flow may feed the formation of the dense low-mass YSO cluster. IRAS 00338+6312 is the most likely powering source of the dominant CO outflow. A compact blue-shifted outflow from RNO 1C is also identified.

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.

The physical and chemical properties of the $\rho$ Ophiuchi A core were studied using 1.3 mm continuum and molecular lines such as C$^{18}$O, C$^{17}$O, CH$_3$OH and H$_2$CO observed with the Submillimeter Array (SMA). The continuum and C$^{18}$O data were combined with the single-dish data obtained with the IRAM 30m telescope and the James Clerk Maxwell Telescope (JCMT), respectively. The combined 1.3 mm continuum map reveals three major sources, SM1, SM1N and VLA1623 embedded in the extended emission running along the north-south direction, and two additional compact condensations in the continuum ridge connecting SM1 and VLA1623.The spatial distribution of the C$^{18}$O emission is significantly different from that of the continuum emission; the C$^{18}$O emission is enhanced at the eastern and western edges of the continuum ridge, with its peak brightness temperature of 40--50 K. This supports the picture that the $\rho$-Oph A core is heated externally from the nearby stars Oph S1 and HD147889. In contrast, the C$^{18}$O intensity is lower than 15--20 K at the center of the ridge where the continuum emission is bright. The C$^{18}$O abundance decreases inside the ridge, and shows anti-correlation with the N$_2$H$^+$ abundance. However, both C$^{18}$O and N$_2$H$^+$ show strong depletion at the Class 0 protostar VLA1623, implying that the dense gas surrounding VLA1623 is colder than the freeze-out temperature of N$_2$. The blue- and red-shifted components of CH$_3$OH and H$_2$CO lines are seen at SM1, suggesting outflow activity of embedded source in SM1, although the spatial distributions do not show clear bipolarity.

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

HH 211-mms is one of the youngest Class 0 protostellar systems in Perseus at ~ 235 pc away. We have mapped its central region at up to ~ 7 AU (0.03") resolution. A dusty disk is seen deeply embedded in a flattened envelope, with an intensity jump in dust continuum at ~ 350 GHz. It is nearly edge-on and is almost exactly perpendicular to the jet axis. It has a size of ~ 30 au along the major axis. It is geometrically thick, indicating that the (sub)millimeter light emitting grains have yet to settle to the midplane. Its inner part is expected to have transformed into a Keplerian rotating disk with a radius of ~ 10 au. A rotating disk atmosphere and a compact rotating bipolar outflow are detected in SO. The outflow fans out from the inner disk surfaces and is rotating in the same direction as the flattened envelope, and hence could trace a disk wind carrying away angular momentum from the inner disk. From the rotation of the disk atmosphere, the protostellar mass is estimated to be <~ 50 M_Jup. Together with results from the literature, our result favors a model where the disk radius grows linearly with the protostellar mass, as predicted by models of pre-stellar dense core evolution that asymptotes to an $r^{-1}$ radial profile for both the column density and angular velocity.

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.

Magnetic field is one of the key agents that play a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained due to the limitations in observations. We study the magnetic field in the massive infrared dark cloud G035.39-00.33 from dust continuum polarization observations at 850 $\micron$ with SCUBA-2/POL-2 at JCMT. The magnetic field tends to be perpendicular to the densest part of the main filament (F$_{M}$), whereas it has a less defined relative orientation in the rest of the structure, where it tends to be parallel to some diffuse regions. A mean plane-of-the-sky magnetic field strength of $\sim$50 $\mu$G for F$_{M}$ is obtained using Davis-Chandrasekhar-Fermi method. Based on $^{13}$CO (1-0) line observations, we suggest a formation scenario of F$_{M}$ due to large-scale ($\sim$10 pc) cloud-cloud collision. Using additional NH$_3$ line data, we estimate that F$_{M}$ will be gravitationally unstable if it is only supported by thermal pressure and turbulence. The northern part of F$_{M}$, however, can be stabilized by a modest additional support from the local magnetic field. The middle and southern parts of F$_{M}$ are likely unstable even if the magnetic field support is taken into account. We claim that the clumps in F$_{M}$ may be supported by turbulence and magnetic fields against gravitational collapse. Finally, we identified for the first time a massive ($\sim$200 M$_{\sun}$), collapsing starless clump candidate, "c8", in G035.39-00.33. The magnetic field surrounding "c8" is likely pinched, hinting at an accretion flow along the filament.

Based on the 850 $\mu$m dust continuum data from SCUBA-2 at James Clerk Maxwell Telescope (JCMT), we compare overall properties of Planck Galactic Cold Clumps (PGCCs) in the $\lambda$ Orionis cloud to those of PGCCs in the Orion A and B clouds. The Orion A and B clouds are well known active star-forming regions, while the $\lambda$ Orionis cloud has a different environment as a consequence of the interaction with a prominent OB association and a giant Hii region. PGCCs in the $\lambda$ Orionis cloud have higher dust temperatures ($Td=16.13\pm0.15$ K) and lower values of dust emissivity spectral index ($ \beta=1.65\pm0.02$) than PGCCs in the Orion A (Td=13.79$\pm 0.21$K, $\beta=2.07\pm0.03$) and Orion B ($Td=13.82\pm0.19$K, $\beta=1.96\pm0.02$) clouds. We find 119 sub-structures within the 40 detected PGCCs and identify them as cores. Of total 119 cores, 15 cores are discovered in the $\lambda$ Orionis cloud, while 74 and 30 cores are found in the Orion A and B clouds, respectively. The cores in the $\lambda$ Orionis cloud show much lower mean values of size R=0.08 pc, column density N(H2)=$(9.5\pm1.2) \times 10^{22}$ cm$^{-2}$, number density n(H2)=$(2.9 \pm 0.4)\times10^{5}$ cm$^{-3}$, and mass $M_{core}$=$1.0\pm0.3$ M$_{\odot}$ compared to the cores in the Orion A (R=0.11pc, $N(H2)=(2.3\pm0.3) \times 10^{23}$ cm$^{-2}$, n(H2)=$(3.8\pm0.5) \times 10^{5}$cm$^{-3}$, and $M_{core}$=$2.4 \pm 0.3$ M$_{\odot}$) and Orion B (R=0.16pc, N(H2)=$(3.8 \pm 0.4) \times 10^{23}$cm$^{-2}$, n(H2)=$(15.6\pm1.8)\times10^{5}$ cm$^{-3}$, and $M_{core}$= $2.7\pm0.3$ M$_{\odot}$) clouds. These core properties in the $\lambda$ Orionis cloud can be attributed to the photodissociation and external heating by the nearby Hii region, which may prevent the PGCCs from forming gravitationally bound structures and eventually disperse them. These results support the idea of negative stellar feedback on core formation.

We performed the new JVLA full polarization observations at 40-48 GHz (6.3-7.5 mm) towards the nearby ($d$ $=$147$\pm$3.4 pc) Class 0 YSO IRAS 16293-2422, and compare with the previous SMA observations reported by Rao et al. (2009; 2014). We observed the quasar J1407+2827 which is weakly polarized and can be used as a leakage term calibrator for $<$9 GHz observations, to gauge the potential residual polarization leakage after calibration. We did not detect Stokes Q, U, and V intensities from the observations of J1407+2827, and constrain (3-$\sigma$) the residual polarization leakage after calibration to be $\lesssim$0.3\%. We detect linear polarization from one of the two binary components of our target source, IRAS\,16293-2422\,B. The derived polarization position angles from our observations are in excellent agreement with those detected from the previous observations of the SMA, implying that on the spatial scale we are probing ($\sim$50-1000 au), the physical mechanisms for polarizing the continuum emission do not vary significantly over the wavelength range of $\sim$0.88-7.5 mm. We hypothesize that the observed polarization position angles trace the magnetic field which converges from large scale to an approximately face-on rotating accretion flow. In this scenario, magnetic field is predominantly poloidal on $>$100 au scales, and becomes toroidal on smaller scales. However, this interpretation remains uncertain due to the high dust optical depths at the central region of IRAS\,16293-2422\,B and the uncertain temperature profile. We suggest that dust polarization at wavelengths comparable or longer than 7\,mm may still trace interstellar magnetic field. Future sensitive observations of dust polarization in the fully optically thin regime will have paramount importance for unambiguously resolving the magnetic field configuration.

HH 212 is a Class 0 protostellar system found to host a "hamburger"-shaped dusty disk with a rotating disk atmosphere and a collimated SiO jet at a distance of ~ 400 pc. Recently, a compact rotating outflow has been detected in SO and SO2 toward the center along the jet axis at ~ 52 au (0.13") resolution. Here we resolve the compact outflow into a small-scale wide-opening rotating outflow shell and a collimated jet, with the observations in the same S-bearing molecules at ~ 16 au (0.04") resolution. The collimated jet is aligned with the SiO jet, tracing the shock interactions in the jet. The wide-opening outflow shell is seen extending out from the inner disk around the SiO jet and has a width of ~ 100 au. It is not only expanding away from the center, but also rotating around the jet axis. The specific angular momentum of the outflow shell is ~ 40 au km/s. Simple modeling of the observed kinematics suggests that the rotating outflow shell can trace either a disk wind or disk material pushed away by an unseen wind from the inner disk or protostar. We also resolve the disk atmosphere in the same S-bearing molecules, confirming the Keplerian rotation there.

Episodic accretion has been proposed as a solution to the long-standing luminosity problem in star formation; however, the process remains poorly understood. We present observations of line emission from N2H+ and CO isotopologues using the Atacama Large Millimeter/submillimeter Array (ALMA) in the envelopes of eight Very Low Luminosity Objects (VeLLOs). In five of the sources the spatial distribution of emission from N2H+ and CO isotopologues shows a clear anti-correlation. It is proposed that this is tracing the CO snow line in the envelopes: N2H+ emission is depleted toward the center of these sources in contrast to the CO isotopologue emission which exhibits a peak. The positions of the CO snow lines traced by the N2H+ emission are located at much larger radii than those calculated using the current luminosities of the central sources. This implies that these five sources have experienced a recent accretion burst because the CO snow line would have been pushed outwards during the burst due to the increased luminosity of the central star. The N2H+ and CO isotopologue emission from DCE161, one of the other three sources, is most likely tracing a transition disk at a later evolutionary stage. Excluding DCE161, five out of seven sources (i.e., ~70%) show signatures of a recent accretion burst. This fraction is larger than that of the Class 0/I sources studied by J{\o}rgensen et al. (2015) and Frimann et al. (2016) suggesting that the interval between accretion episodes in VeLLOs is shorter than that in Class 0/I sources.

Analysis of all-sky Planck submillimetre observations and the IRAS 100um data has led to the detection of a population of Galactic cold clumps. The clumps can be used to study star formation and dust properties in a wide range of Galactic environments. Our aim is to measure dust spectral energy distribution (SED) variations as a function of the spatial scale and the wavelength. We examine the SEDs at large scales using IRAS, Planck, and Herschel data. At smaller scales, we compare with JCMT/SCUBA-2 850um maps with Herschel data that are filtered using the SCUBA-2 pipeline. Clumps are extracted using the Fellwalker method and their spectra are modelled as modified blackbody functions. According to IRAS and Planck data, most fields have dust colour temperatures T_C ~ 14-18K and opacity spectral index values of beta=1.5-1.9. The clumps/cores identified in SCUBA-2 maps have T~ 13K and similar beta values. There are some indications of the dust emission spectrum becoming flatter at wavelengths longer than 500um. In fits involving Planck data, the significance is limited by the uncertainty of the corrections for CO line contamination. The fits to the SPIRE data give a median beta value slightly above 1.8. In the joint SPIRE and SCUBA-2 850um fits the value decreases to beta ~1.6. Most of the observed T-beta anticorrelation can be explained by noise. The typical submillimetre opacity spectral index beta of cold clumps is found to be ~1.7. This is above the values of diffuse clouds but lower than in some previous studies of dense clumps. There is only tentative evidence of T-beta anticorrelation and beta decreasing at millimetre wavelengths.

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.

Young stellar objects (YSOs) may undergo periods of active accretion (outbursts), during which the protostellar accretion rate is temporarily enhanced by a few orders of magnitude. Whether or not these accretion outburst YSOs possess similar dust/gas reservoirs to each other, and whether or not their dust/gas reservoirs are similar as quiescent YSOs, are issues not yet clarified. The aim of this work is to characterize the millimeter thermal dust emission properties of a statistically significant sample of long and short duration accretion outburst YSOs (i.e., FUors and EXors) and the spectroscopically identified candidates of accretion outbursting YSOs (i.e., FUor-like objects). We have carried out extensive Submillimeter Array (SMA) observations mostly at $\sim$225 GHz (1.33 mm) and $\sim$272 GHz (1.10 mm), from 2008 to 2017. We covered accretion outburst YSOs located at $<$1 kpc distances from the solar system. We analyze all the existing SMA data of such objects, both published and unpublished, in a coherent way to present a millimeter interferometric database of 29 objects. We obtained 21 detections at $>$3-$\sigma$ significance. Detected sources except for the two cases of V883 Ori and NGC 2071 MM3 were observed with $\sim$1$"$ angular resolution. Overall our observed targets show a systematically higher millimeter luminosity distribution than those of the $M_{*}>$0.3 $M_{\odot}$ Class II YSOs in the nearby ($\lesssim$400 pc) low-mass star-forming molecular clouds (e.g., Taurus, Lupus, Upp Scorpio, and Chameleon I). In addition, at 1 mm our observed confirmed binaries or triple-system sources are systematically fainter than the rest of the sources even though their 1 mm fluxes are broadly distributed. We may have detected $\sim$30-60\% millimeter flux variability from V2494 Cyg and V2495 Cyg, from the observations separated by $\sim$1 year.

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.

The central problem in forming a star is the angular momentum in the circumstellar disk which prevents material from falling into the central stellar core. An attractive solution to the "angular momentum problem" appears to be the ubiquitous (low-velocity and poorly-collimated) molecular outflows and (high-velocity and highly-collimated) protostellar jets accompanying the earliest phase of star formation that remove angular momentum at a range of disk radii. Previous observations suggested that outflowing material carries away the excess angular momentum via magneto-centrifugally driven winds from the surfaces of circumstellar disks down to ~ 10 AU scales, allowing the material in the outer disk to transport to the inner disk. Here we show that highly collimated protostellar jets remove the residual angular momenta at the ~ 0.05 AU scale, enabling the material in the innermost region of the disk to accrete toward the central protostar. This is supported by the rotation of the jet measured down to ~ 10 AU from the protostar in the HH 212 protostellar system. The measurement implies a jet launching radius of ~ 0.05_{-0.02}^{+0.05} AU on the disk, based on the magneto-centrifugal theory of jet production, which connects the properties of the jet measured at large distances to those at its base through energy and angular momentum conservation.

HH 212 is a nearby (400 pc) Class 0 protostellar system recently found to host a "hamburger"-shaped dusty disk with a radius of ~ 60 AU, deeply embedded in an infalling-rotating flattened envelope. We have spatially resolved this envelope-disk system with the Atacama Large Millimeter/submillimeter Array at up to ~ 16 AU (0.04") resolution. The envelope is detected in HCO+ J=4-3 down to the dusty disk. Complex organic molecules (COMs) and doubly deuterated formaldehyde (D2CO) are detected above and below the dusty disk within ~ 40 AU of the central protostar. The COMs are methanol (CH3OH), deuterated methanol (CH2DOH), methyl mercaptan (CH3SH), and formamide (NH2CHO, a prebiotic precursor). We have modeled the gas kinematics in HCO+ and COMs, and found a centrifugal barrier at a radius of ~ 44 AU, within which a Keplerian rotating disk is formed. This indicates that HCO+ traces the infalling-rotating envelope down to centrifugal barrier and COMs trace the atmosphere of a Keplerian rotating disk within the centrifugal barrier. The COMs are spatially resolved for the first time, both radially and vertically, in the atmosphere of a disk in the earliest, Class 0 phase of star formation. Our spatially resolved observations of COMs favor their formation in the disk rather than a rapidly infalling (warm) inner envelope. The abundances and spatial distributions of the COMs provide strong constraints on models of their formation and transport in low-mass star formation.

In the earliest (so-called "Class 0") phase of sunlike (low-mass) star formation, circumstellar disks are expected to form, feeding the protostars. However, such disks are difficult to resolve spatially because of their small sizes. Moreover, there are theoretical difficulties in producing such disks in the earliest phase, due to the retarding effects of magnetic fields on the rotating, collapsing material (so-called "magnetic braking"). With the Atacama Large Millimeter/submillimeter Array (ALMA), it becomes possible to uncover such disks and study them in detail. HH 212 is a very young protostellar system. With ALMA, we not only detect but also spatially resolve its disk in dust emission at submillimeter wavelength. The disk is nearly edge-on and has a radius of ~ 60 AU. Interestingly, it shows a prominent equatorial dark lane sandwiched between two brighter features, due to relatively low temperature and high optical depth near the disk midplane. For the first time, this dark lane is seen at submillimeter wavelength, producing a "hamburger"-shaped appearance that is reminiscent of the scattered-light image of an edge-on disk in optical and near infrared. Our observations open up an exciting possibility of directly detecting and characterizing small disks around the youngest protostars through high-resolution imaging with ALMA, which provides strong constraints on theories of disk formation.

We are motivated by the recent measurements of dust opacity indices beta around young stellar objects (YSOs), which suggest that efficient grain growth may have occurred earlier than the Class I stage. The present work makes use of abundant archival interferometric observations at submillimeter,millimeter, and centimeter wavelength bands to examine grain growth signatures in the dense inner regions (<1000 AU) of nine Class 0 YSOs. A systematic data analysis is performed to derive dust temperatures, optical depths, and dust opacity indices based on single-component modified black body fittings to the spectral energy distributions (SEDs). The fitted dust opacity indices (beta) are in a wide range of 0.3 to 2.0 when single-component SED fitting is adopted. Four out of the nine observed sources show beta lower than 1.7, the typical value of the interstellar dust. Low dust opacity index (or spectral index) values may be explained by the effect of dust grain growth, which makes beta<1.7. Alternatively, the very small observed values of beta may be interpreted by the presence of deeply embedded hot inner disks, which only significantly contribute to the observed fluxes at long wavelength bands. This possibility can be tested by the higher angular resolution imaging observations of ALMA, or more detailed sampling of SEDs in the millimeter and centimeter bands. The beta values of the remaining five sources are close to or consistent with 1.7, indicating that grain growth would start to significantly reduce the values of beta no earlier than the late-Class 0 stage for these YSOs.

Preplanetary nebulae and planetary nebulae are evolved, mass-losing stellar objects that show a wide variety of morphologies. Many of these nebulae consist of outer structures that are nearly spherical (spiral/shell/arc/halo) and inner structures that are highly asymmetric (bipolar/multipolar). The coexistence of such geometrically distinct structures is enigmatic because it hints at the simultaneous presence of both wide and close binary interactions, a phenomenon that has been attributed to stellar binary systems with eccentric orbits. Here we report new high-resolution molecular-line observations of the circumstellar spiral-shell pattern of AFGL 3068, an asymptotic giant branch star transitioning to the preplanetary nebula phase. The observations clearly reveal that the dynamics of the mass loss is influenced by the presence of an eccentric-orbit binary. This quintessential object opens a new window on the nature of deeply embedded binary stars through the circumstellar spiral-shell patterns that reside at distances of several thousand Astronomical Units from the stars.

We resolved FU Ori at 29-37 GHz using the JVLA with $\sim$0$''$.07 resolution, and performed the complementary JVLA 8-10 GHz observations, the SMA 224 GHz and 272 GHz observations, and compared with archival ALMA 346 GHz observations to obtain the SEDs. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at $>$29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency). The more sophisticated models for SEDs suggest that the $>$29 GHz emission is contributed by a combination of free-free emission from ionized gas, and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks ($\lesssim$0.1 au) has been sublimated, and thus the disks are no more well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal and the MRI instabilities when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate.

We observed thirteen Planck cold clumps with the James Clerk Maxwell Telescope/SCUBA-2 and with the Nobeyama 45 m radio telescope. The N$_2$H$^+$ distribution obtained with the Nobeyama telescope is quite similar to SCUBA-2 dust distribution. The 82 GHz HC$_3$N, 82 GHz CCS, and 94 GHz CCS emission are often distributed differently with respect to the N$_2$H$^+$ emission. The CCS emission, which is known to be abundant in starless molecular cloud cores, is often very clumpy in the observed targets. We made deep single-pointing observations in DNC, HN$^{13}$C, N$_2$D$^+$, cyclic-C$_3$H$_2$ toward nine clumps. The detection rate of N$_2$D$^+$ is 50\%. Furthermore, we observed the NH$_3$ emission toward 15 Planck cold clumps to estimate the kinetic temperature, and confirmed that most of targets are cold ($\lesssim$ 20 K). In two of the starless clumps observe, the CCS emission is distributed as it surrounds the N$_2$H$^+$ core (chemically evolved gas), which resembles the case of L1544, a prestellar core showing collapse. In addition, we detected both DNC and N$_2$D$^+$. These two clumps are most likely on the verge of star formation. We introduce the Chemical Evolution Factor (CEF) for starless cores to describe the chemical evolutionary stage, and analyze the observed Planck cold clumps.

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.

$Aims.$ We study the relation between the jet and the outflow in the IRAS 04166+2706 protostar. This Taurus protostar drives a molecular jet that contains multiple emission peaks symmetrically located from the central source. The protostar also drives a wide-angle outflow consisting of two conical shells. $Methods.$ We have used the Atacama Large Millimeter/submillimeter Array (ALMA) interferometer to observe two fields along the IRAS 04166+2706 jet. The fields were centered on a pair of emission peaks that correspond to the same ejection event, and were observed in CO(2-1), SiO(5-4), and SO(65-54). $ Results.$ Both ALMA fields present spatial distributions that are approximately elliptical and have their minor axes aligned with the jet direction. As the velocity increases, the emission in each field moves gradually across the elliptical region. This systematic pattern indicates that the emitting gas in each field lies in a disk-like structure that is perpendicular to the jet axis and is expanding away from the jet. A small degree of curvature in the first-moment maps indicates that the disks are slightly curved in the manner expected for bow shocks moving away from the IRAS source. A simple geometrical model confirms that this scenario fits the main emission features. $Conclusions.$ The emission peaks in the IRAS 04166+2706 jet likely represent internal bow shocks where material is being ejected laterally away from the jet axis. While the linear momentum of the ejected gas is dominated by the component in the jet direction, the sideways component is not negligible, and can potentially affect the distribution of gas in the surrounding outflow and core.

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

We report new Karl G. Jansky Very Large Array (JVLA), 0$"$.5 angular resolution observations of linearly polarized continuum emission at 6.9 mm, towards the Class 0 young stellar object (YSO) NGC1333 IRAS4A. This target source is a collapsing dense molecular core, which was resolved at short wavelengths to have hourglass shaped B-field configuration. We compare these 6.9 mm observations with previous polarization Submillimeter Array (SMA) observations at 0.88 mm, which have comparable angular resolution ($\sim$0$"$7). We found that at the same resolution, the observed polarization position angles at 6.9 mm are slightly deviated from those observed at 0.88 mm. Due to the lower optical depth of the emission at 6.9 mm, and the potential effect of dust grain growth, the new JVLA observations are likely probing B-field alignments in regions interior to those sampled by the previous polarization observations at higher frequencies. Our understanding can be improved by more sensitive observations, and observations for the more extended spatial scales.

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

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

HH 212 is a nearby (400 pc) highly collimated protostellar jet powered by a Class 0 source in Orion. We have mapped the inner 80" (~ 0.16 pc) of the jet in SiO (J=8-7) and CO (J=3-2) simultaneously at ~ 0.5 resolution with the Atacama Millimeter/Submillimeter Array at unprecedented sensitivity. The jet consists of a chain of knots, bow shocks, and sinuous structures in between. As compared to that seen in our previous observations with the Submillimeter Array, it appears to be more continuous, especially in the northern part. Some of the knots are now seen associated with small bow shocks, with their bow wings curving back to the jet axis, as seen in pulsed jet simulations. Two of them are reasonably resolved, showing kinematics consistent with sideways ejection, possibly tracing the internal working surfaces formed by a temporal variation in the jet velocity. In addition, nested shells are seen in CO around the jet axis connecting to the knots and bow shocks, driven by them. The proper motion of the jet is estimated to be ~ 115+-50 km/s, comparing to our previous observations. The jet has a small semi-periodical wiggle, with a period of ~ 93 yrs. The amplitude of the wiggle first increases with the distance from the central source and then stays roughly constant. One possible origin of the wiggle could be the kink instability in a magnetized jet.

HH 211 is a young Class 0 protostellar system, with a flattened envelope, a possible rotating disk, and a collimated jet. We have mapped it with the Submillimeter Array in 341.6 GHz continuum and SiO J=8-7 at ~ 0.6 resolution. The continuum traces the thermal dust emission in the flattened envelope and the possible disk. Linear polarization is detected in the continuum in the flattened envelope. The field lines implied from the polarization have different orientations, but they are not incompatible with current gravitational collapse models, which predict different orientation depending on the region/distance. Also, we might have detected for the first time polarized SiO line emission in the jet due to the Goldreich-Kylafis effect. Observations at higher sensitivity are needed to determine the field morphology in the jet.

Two submm/mm sources in the Barnard 1b (B1-b) core, B1-bN and B1-bS, have been studied in dust continuum, H13CO+ J=1-0, CO J=2-1, 13CO J=2-1, and C18O J=2-1. The spectral energy distributions of these sources from the mid-IR to 7 mm are characterized by very cold temperatures of T_dust < 20 K and low bolometric luminosities of 0.15-0.31 L_sun. The internal luminosities of B1-bN and B1-bS are estimated to be <0.01-0.03 L_sun and ~0.1-0.2 L_sun, respectively. Millimeter interferometric observations have shown that these sources have already formed central compact objects of ~100 AU sizes. Both B1-bN and B1-bS are driving the CO outflows with low characteristic velocities of ~2-4 km/s. The fractional abundance of H13CO+ at the positions of B1-bN and B1-bS is lower than the canonical value by a factor of 4-8. This implies that significant fraction of CO is depleted onto dust grains in dense gas surrounding these sources. The observed physical and chemical properties suggest that B1-bN and B1-bS are in the earlier evolutionary stage than most of the known Class 0 protostars. Especially, the properties of B1-bN agree with those of the first hydrostatic core predicted by the MHD simulations. The CO outflow was also detected in the mid-IR source located at ~15" from B1-bS. Since the dust continuum emission was not detected in this source, the circumstellar material surrounding this source is less than 0.01 M_sun. It is likely that the envelope of this source was dissipated by the outflow from the protostar that is located to the southwest of B1-b.

HH 212 is a nearby (400 pc) Class 0 protostellar system showing several components that can be compared with theoretical models of core collapse. We have mapped it in 350 GHz continuum and HCO+ J=4-3 emission with ALMA at up to ~ 0.4" resolution. A flattened envelope and a compact disk are seen in continuum around the central source, as seen before. The HCO+ kinematics shows that the flattened envelope is infalling with small rotation (i.e., spiraling) into the central source, and thus can be identified as a pseudodisk in the models of magnetized core collapse. Also, the HCO+ kinematics shows that the disk is rotating and can be rotationally supported. In addition, to account for the missing HCO+ emission at low-redshifted velocity, an extended infalling envelope is required, with its material flowing roughly parallel to the jet axis toward the pseudodisk. This is expected if it is magnetized with an hourglass B-field morphology. We have modeled the continuum and HCO+ emission of the flattened envelope and disk simultaneously. We find that a jump in density is required across the interface between the pseudodisk and the disk. A jet is seen in HCO+ extending out to ~ 500 AU away from the central source, with the peaks upstream of those seen before in SiO. The broad velocity range and high HCO+ abundance indicate that the HCO+ emission traces internal shocks in the jet.

The molecular outflow from IRAS 04166+2706 was mapped with the Submillimeter Array (SMA) at 350 GHz continuum and CO J = 3$-$2 at an angular resolution of ~1 arcsec. The field of view covers the central arc-minute, which contains the inner four pairs of knots of the molecular jet. On the channel map, conical structures are clearly present in the low velocity range (|V$-$V$_0$|$<$10 km s$^{-1}$), and the highly collimated knots appear in the Extremely High Velocity range (EHV, 50$>$|V$-$V$_0$|$>$30 km s$^{-1}$). The higher angular resolution of ~1 arcsec reveals the first blue-shifted knot (B1) that was missing in previous PdBI observation of Sant\'iago-Garc\'ia et al. (2009) at an offset of ~6 arcsec to the North-East of the central source. This identification completes the symmetric sequence of knots in both the blue- and red-shifted lobes of the outflow. The innermost knots R1 and B1 have the highest velocities within the sequence. Although the general features appear to be similar to previous CO J = 2$-$1 images in Sant\'iago-Garc\'ia et al. (2009), the emission in CO J = 3$-$2 almost always peaks further away from the central source than that of CO J = 2$-$1 in the red-shifted lobe of the channel maps. This gives rise to a gradient in the line-ratio map of CO J = 3$-$2/J = 2$-$1 from head to tail within a knot. A large velocity gradient (LVG) analysis suggests that the differences may reflect a higher gas kinetic temperature at the head. We also explore possible constraints imposed by the non-detection of SiO J = 8$-$7.

We use the Submillimeter Array to observe, at 1.4 mm, the blue-lobe of the L1157 outflow at high spatial resolution (~ 3"). We detected SiO, H_2CO, and CH_3OH lines from several molecular clumps that constitute the outflow. All three molecules were detected along the wall of the inner cavity that is supposedly related with the later ejection event. On the other hand, no emission was detected towards positions related to an old ejection episode, likely due to space filtering from the interferometer. The H_2CO and CH_3OH emission is detected only at velocities close to the systemic velocity. The spatial distributions of the H_2CO and CH_3OH are similar. These emission lines trace the U-shaped structure seen in the mid-infrared image. In contrast, the SiO emission is detected in wider velocity range with a peak at ~14 km s/s blue-shifted from the systemic velocity. The SiO emission is brightest at the B1 position, which corresponds to the apex of the U-shaped structure. There are two compact SiO clumps along the faint arc-like feature to the east of the U-shaped structure. At the B1 position, there are two velocity components; one is a compact clump with a size of ~1500 AU seen in the high-velocity and the other is an extended component with lower velocities. The kinematic structure at the B1 position is different from that expected in a single bow shock. It is likely that the high-velocity SiO clump at the B1 position is kinetically independent from the low-velocity gas. The line ratio between SiO (5--4) and SiO (2--1) suggests that the high velocity SiO clumps consist of high density gas of n ~ 10^5 - 10^6 cm^-3, which is comparable to the density of the bullets in the extremely high velocity (EHV) jets. It is likely that the high-velocity SiO clumps in L1157 have the same origin as the EHV bullets.

Two submm/mm sources in the Barnard 1b (B1-b) core, B1-bN and B1-bS, have been observed with the Submillimeter Array (SMA) and the Submillimeter Telescope (SMT). The 1.1 mm continuum map obtained with the SMA reveals that the two sources contain spatially compact components, suggesting that they harbor protostars. The N2D+ and N2H+ J=3-2 maps were obtained by combining the SMA and SMT data. The N2D+ map clearly shows two peaks at the continuum positions. The N2H+ map also peaks at the continuum positions, but is more dominated by the spatially extended component. The N2D+/N2H+ ratio was estimated to be ~ 0.2 at the positions of both B1-bN and B1-bS. The derived N2D+/N2H+ ratio is comparable to those of the prestellar cores in the late evolutionary stage and the class 0 protostars in the early evolutionary stage. Although B1-bN is bright in N2H+ and N2D+, this source was barely seen in H13CO+. This implies that the depletion of carbon-bearing molecules is significant in B1-bN. The chemical property suggests that B1-bN is in the earlier evolutionary stage as compared to B1-bS with the H13CO+ counterpart. The N2H+ and N2D+ lines show that the radial velocities of the two sources are different by ~ 0.9 km s-1. However, the velocity pattern along the line through B1-bN and B1-bS suggests that these two sources were not formed out of a single rotating cloud. It is likely that the B1-b core consists of two velocity components, each of which harbors a very young source.

The protostellar jet driven by L1448C was observed in the SiO J=8-7 and CO J=3-2 lines and 350 GHz dust continuum at ~1" resolution with the Submillimeter Array (SMA). A narrow jet from the northern source L1448C(N) was observed in the SiO and the high-velocity CO. The jet consists of a chain of emission knots with an inter-knot spacing of ~2" (500 AU) and a semi-periodic velocity variation. The innermost pair of knots, which are significant in the SiO map but barely seen in the CO, are located at ~1" (250 AU) from the central source, L1448C(N). Since the dynamical time scale for the innermost pair is only ~10 yr, SiO may have been formed in the protostellar wind through the gas-phase reaction, or been formed on the dust grain and directly released into the gas phase by means of shocks. It is found that the jet is extremely active with a mechanical luminosity of ~7 L_sun, which is comparable to the bolometric luminosity of the central source (7.5 L_sun). The mass accretion rate onto the protostar derived from the mass-loss rate is ~10^{-5} M_sun/yr. Such a high mass accretion rate suggests that the mass and the age of the central star are 0.03-0.09 M_sun and (4-12)x10^3 yr, respectively, implying that the central star is in the very early stage of protostellar evolution. The low-velocity CO emission delineates two V-shaped shells with a common apex at L1448C(N). The kinematics of these shells are reproduced by the model of a wide opening angle wind. The co-existence of the highly-collimated jets and the wide-opening angle shells can be explained by the unified X-wind model" in which highly-collimated jet components correspond to the on-axis density enhancement of the wide-opening angle wind. The CO $J$=3--2 map also revealed the second outflow driven by the southern source L1448C(S) located at ~8.3" (2000 AU) from L1448C(N).

HH 211 is a highly collimated jet originating from a nearby young Class 0 protostar. Here is a follow-up study of the jet with our previous observations at unprecedented resolution up to ~ 0.3" in SiO (J=8-7), CO (J=3-2), and SO (N_J=8_9-7_8). SiO, CO, and SO can all be a good tracer of the HH 211 jet, tracing the internal shocks in the jet. Although the emissions of these molecules show roughly the same morphology of the jet, there are detailed differences. In particular, the CO emission traces the jet closer to the source than the SiO and SO emissions. In addition, in the better resolved internal shocks, both the CO and SO emission are seen slightly ahead of the SiO emission. The jet is clearly seen on both sides of the source with more than one cycle of wiggle. The wiggle is reflection-symmetric about the source and can be reasonably fitted by an orbiting source jet model. The best-fit parameters suggest that the source itself could be a very low-mass protobinary with a total mass of ~ 60 M_Jup and a binary separation of ~ 4.6 AU. The abundances of SiO and SO in the gas phase are found to be highly enhanced in the jet as compared to the quiescent molecular clouds, even close to within 300 AU from the source where the dynamical time scale is <10 yrs. The abundance enhancements of these molecules are closely related to the internal shocks. The detected SiO is either the consequence of the release of Si-bearing material from dust grains or of its formation via gas chemistry in the shocks. The SO, on the other hand, seems to form via gas chemistry in the shocks.

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.

HH 211 is a nearby young protostellar system with a highly collimated jet. We have mapped it in 352 GHz continuum, SiO (J=8-7), and HCO+ (J=4-3) emission at up to ~ 0.2" resolution with the Submillimeter Array (SMA). The continuum source is now resolved into two sources, SMM1 and SMM2, with a separation of ~ 84 AU. SMM1 is seen at the center of the jet, probably tracing a (inner) dusty disk around the protostar driving the jet. SMM2 is seen to the southwest of SMM1 and may trace an envelope-disk around a small binary companion. A flattened envelope-disk is seen in HCO+ around SMM1 with a radius of ~ 80 AU perpendicular to the jet axis. Its velocity structure is consistent with a rotation motion and can be fitted with a Keplerian law that yields a mass of ~ 50+-15 Jupiter mass (a mass of a brown dwarf) for the protostar. Thus, the protostar could be the lowest mass source known to have a collimated jet and a rotating flattened envelope-disk. A small-scale (~ 200 AU) low-speed (~ 2 km/s) outflow is seen in HCOP+ around the jet axis extending from the envelope-disk. It seems to rotate in the same direction as the envelope-disk and may carry away part of the angular momentum from the envelope-disk. The jet is seen in SiO close to ~ 100 AU from SMM1. It is seen with a "C-shaped" bending. It has a transverse width of <~ 40 AU and a velocity of ~ 170+-60 km/s. A possible velocity gradient is seen consistently across its innermost pair of knots, with ~ 0.5 km/s at ~ 10 AU, consistent with the sense of rotation of the envelope-disk. If this gradient is an upper limit of the true rotational gradient of the jet, then the jet carries away a very small amount of angular momentum of ~ 5 AU km/s and thus must be launched from the very inner edge of the disk near the corotation radius.

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 mapped the protostellar jet HH 211 in 342 GHz continuum, SiO ($J=8-7$), and CO ($J=3-2$) emission at $\sim$ \arcs{1} resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum at the center of the jet, tracing an envelope and a possible optically thick compact disk (with a size $<$ 130 AU) around the protostar. A knotty jet is seen in CO and SiO as in \H2{}, but extending closer to the protostar. It consists of a chain of knots on each side of the protostar, with an interknot spacing of $\sim$ \arcs{2}$-$\arcs{3} or 600$-$900 AU and the innermost pair of knots at only $\sim$ \arcsa{1}{7} or 535 AU from the protostar. These knots likely trace unresolved internal (bow) shocks (i.e., working surfaces) in the jet, with a velocity range up to $\sim$ 25 \vkm{}. The two-sided mass-loss rate of the jet is estimated to be $\sim (0.7-2.8)\times 10^{-6}$ \solarmass{} yr$^{-1}$. The jet is episodic, precessing, and bending. A velocity gradient is seen consistently across two bright SiO knots (BK3 and RK2) perpendicular to the jet axis, with $\sim$ 1.5$\pm$0.8 \vkm{} at $\sim$ 30$\pm$15 AU, suggesting a presence of a jet rotation. The launching radius of the jet, derived from the potential jet rotation, is $\sim$ 0.15$-$0.06 AU in the inner disk.

We report on the results of a Submillimeter Array interferometric observation of the proto-planetary nebula CRL 618 in the 12CO J=6-5 line. With the new capability of SMA enabling us to use two receivers at a time, we also observed simultaneously in the 12CO J=2-1 and 13CO J=2-1 lines. The 12CO J=6-5 and 13CO J=2-1 lines were first interferometrically observed toward CRL 618. The flux of the high velocity component of the 12CO J=6-5 line is almost fully recovered, while roughly 80% of the flux of the low velocity component is resolved out. The low recovery rate suggests that the emission region of the low velocity component of the 12CO J=6-5 line is largely extended. Continuum emission is detected both at 230 and 690 GHz. The flux of the 690 GHz continuum emission seems to be partially resolved out, suggesting dust emission partly contaminates the 690 GHz continuum flux. The cavity structure, which has been confirmed in a previous observation in the 12CO J=2-1 line, is not clearly detected in the 12CO J=6-5 line, and only the south wall of the cavity is detected. This result suggests that the physical condition of the molecular envelope of CRL 618 is not exactly axial symmetric.

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

HH 212 is a nearby (460 pc) protostellar jet discovered in H$_2$ powered by a Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. It is highly collimated and symmetric with matched pairs of bow shocks on either side of the source. We have mapped it in 850 $\mu$m continuum, SiO ($J=8-7$), CO ($J=3-2$), SO ($N_J=8_9-7_8$), HCO$^+$ ($J = 4-3$), and H$^{13}$CO$^+$ ($J = 4-3$) emission simultaneously at $\sim$ 1$''$ resolution with the Submillimeter Array (SMA). Thermal dust emission is seen in continuum around the source, mainly arising from an inner envelope (i.e., the inner part of a previously seen flattened envelope) and a possible disk. The inner envelope is also seen with rotation in CO, HCO$^+$, and probably SO. Like H$_2$ emission, CO and SiO emission are seen along the jet axis but extending closer to the source, tracing the bow shocks with a broad range of velocities and the continuous structures in between. SO emission is seen only around the source, forming a jetlike structure extending along the jet axis from the source, likely tracing the jet near the launching region. The jet is episodic and bending. It may also be slightly precessing as the jetlike SO structure shows a slight S-shaped symmetry about the source. A hint of jet rotation is also seen across the jet axis. Internal outflow shells are seen in CO and HCO$^+$, associated with the bow shocks in the inner part of the jet. The bases of the HCO$^+$ shells are seen with a hint of rotation similar to that seen in the inner envelope, probably consisted mainly of the material extended from the inner envelope and even the possible disk.

We present millimeter interferometric observations of the molecular envelope around the red supergiant VY CMa with the SubMillimeter Array (SMA). The high angular resolution (< 2'') allows us to derive the structure of the envelope as observed in the 1.3 mm continuum, 12CO(2-1), 13CO(2-1) and SO(6,5-5,4) lines emission. The circumstellar envelope is resolved into three components: a dense, compact and dusty central component, embedded in a more diffuse and extended envelope plus a high velocity component. We construct a simple model, consisting of a spherically symmetric slowly expanding envelope and bipolar outflows with a wide opening angle (~ 120 deg.) viewed close to the line of sight (i = 15 deg.). Our model can explain the main features of the SMA data and previous single-dish CO multi-line observations. An episode of enhanced mass loss along the bipolar direction is inferred from our modelling. The SMA data provide a better understanding of the complicated morphology seen in the optical/IR high resolution observations.

We report a dual-band observation at 223 and 654 GHz (460 micron) toward an ultracompact (UC) HII region, G240.31+0.07, with the Submillimeter Array. With a beam size of 1"5 X 0"8, the dust continuum emission is resolved into two clumps, with clump A coincident well with an H2O maser and the UC HII region. The newly discovered clump, B, about 1"3 (~8.3 X 10^3 AU) to the southwest of clump A, is also associated with H2O masers and may be a more recent star-forming site. The continuum flux densities imply an opacity spectral index of \beta = 1.5 +- 0.3, suggestive of a value lower than the canonical 2.0 found in the interstellar medium and in cold, massive cores. The presence of hot (~100 K) molecular gas is derived by the brightness ratio of two H2CO lines in the 223 GHz band. A radial velocity difference of 2.5 +- 0.4 km/s is found between the two clumps in C18O (6-5) emission. The total (nebular and stellar) mass of roughly 58 Msun in the central region is close to, but not by far larger than, the minimum mass required for the two clumps to be gravitationally bound for binary rotation. Our continuum data do not suggest a large amount of matter associated with the H2 knots that were previously proposed to arise from a massive disk or envelope.

We observed the HH 211 jet in the submillimeter continuum and the CO(3-2) and SiO(8-7) transitions with the Submillimeter Array. The continuum source detected at the center of the outflow shows an elongated morphology, perpendicular to the direction of the outflow axis. The high-velocity emission of both molecules shows a knotty and highly collimated structure. The SiO(8-7) emission at the base of the outflow, close to the driving source, spans a wide range of velocities, from -20 up to 40 km s^{-1}. This suggests that a wide-angle wind may be the driving mechanism of the HH 211 outflow. For distances greater than 5" (1500 AU) from the driving source, emission from both transitions follows a Hubble-law behavior, with SiO(8-7) reaching higher velocities than CO(3-2), and being located upstream of the CO(3-2) knots. This indicates that the SiO(8-7) emission is likely tracing entrained gas very close to the primary jet, while the CO(3-2) is tracing less dense entrained gas. From the SiO(5-4) data of Hirano et al. we find that the SiO(8-7)/SiO(5-4) brightness temperature ratio along the jet decreases for knots far from the driving source. This is consistent with the density decreasing along the jet, from (3-10)x10^6 cm^{-3} at 500 AU to (0.8-4)x10^6 cm^{-3} at 5000 AU from the driving source.

We have mapped the SiO J=5-4 line at 217GHz from the HH211 molecular outflow with the Submillimeter Array (SMA). The high resolution map (1.6''x0.9'') shows that the SiO J=5-4 emission comes from the central narrow jet along the outflow axis with a width of ~0.8'' (~250 AU) FWHM. The SiO jet consists of a chain of knots separated by 3-4'' (~1000 AU) and most of the SiO knots have counterparts in shocked H_2 emission seen in a new, deep VLT near-infrared image of the outflow. A new, innermost pair of knots are discovered at just +/-2'' from the central star. The line ratio between the SiO J=5-4 data and upper limits from the SiO J=1-0 data of Chandler & Richer (2001) suggests that these knots have a temperature in excess of 300-500 K and a density of (0.5-1) x10^7 cm^{-3}. The radial velocity measured for these knots is ~30 km/s, comparable to the maximum velocity seen in the entire jet. The high temperature, high density, and velocity structure observed in this pair of SiO knots suggest that they are closely related to the primary jet launched close to the protostar.

HH 212 is a highly collimated jet discovered in H2 powered by a young Class 0 source, IRAS 05413-0104, in the L1630 cloud of Orion. We have mapped around it in 1.33 mm continuum, 12CO ($J=2-1$), 13CO ($J=2-1$), C18O ($J=2-1$), and SO ($J_K = 6_5-5_4$) emission at $\sim$ \arcs{2.5} resolution with the Submillimeter Array. A dust core is seen in the continuum around the source. A flattened envelope is seen in C18O around the source in the equator perpendicular to the jet axis, with its inner part seen in 13CO. The structure and kinematics of the envelope can be roughly reproduced by a simple edge-on disk model with both infall and rotation. In this model, the density of the disk is assumed to have a power-law index of $p=-1.5$ or -2, as found in other low-mass envelopes. The envelope seems dynamically infalling toward the source with slow rotation because the kinematics is found to be roughly consistent with a free fall toward the source plus a rotation of a constant specific angular momentum. A 12CO outflow is seen surrounding the H2 jet, with a narrow waist around the source. Jetlike structures are also seen in 12CO near the source aligned with the H2 jet at high velocities. The morphological relationship between the H2 jet and the 12CO outflow, and the kinematics of the 12CO outflow along the jet axis are both consistent with those seen in a jet-driven bow shock model. SO emission is seen around the source and the H2 knotty shocks in the south, tracing shocked emission around them.

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.

Using the partially completed Submillimeter Array with five antennas, we have observed the CO J=2-1 and 3-2 emission from the envelope surrounding the carbon star V Hya. The high-angular resolution (2"-4") maps show that V Hya is powering a bipolar molecular jet having an extreme velocity of 70-185 km/s. The axis of this high velocity jet is perpendicular to the major axis of the flattened disk-like envelope, which is expanding with a velocity of ~16 km/s.There is a third kinematic component, a medium-velocity wind having a deprojected velocity of 40-120 km/s moving along the disk plane. Both the high velocity jet and the medium velocity wind have a dynamical time scale of a few hundred years. The flattened structure and the collimated jet observed in V Hya suggests that the formation of asymmetrical structure proceeds while the central star is still in the AGB phase

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 present the CO J=6-5, 4-3, and 3-2 spectra from the blueshifted gas of the outflow driven by the low-mass class 0 protostar in the L1157 dark cloud. Strong submillimeter CO emission lines with T_mb > 30 K have been detected at 63" (~0.13 pc) south from the protostar. It is remarkable that the blue wings in the submillimeter lines are stronger by a factor of 3-4 than that of the CO J=1-0 emission line. The CO line ratios suggest that the blueshifted lobe of this outflow consists of moderately dense gas of n(H_2) = (1-3)x10^4 cm^-3 heated to T_kin = 50-170 K.It is also suggested that the kinetic temperature of the outflowing gas increases from ~80 K near the protostar to ~170 K at the shocked region in the lobe center, toward which the largest velocity dispersion of the CO emission is observed. A remarkable correlation between the kinetic temperature and velocity dispersion of the CO emission along the lobe provides us with direct evidence that the molecular gas at the head of the jet-driven bow shock is indeed heated kinematically. The lower temperature of ~80 K measured at the other shocked region near the end of the lobe is explained if this shock is in a later evolutionary stage, in which the gas has been cooled mainly through radiation of the CO rotational lines.

We have mapped the thermal emission line of SiO (v = 0; J = 2-1) associated with the quadrupolar molecular outflow driven by the very cold far-infrared source IRAS 16293-2422. The SiO emission is significantly enhanced in the northeastern red lobe and at the position ~50" east of the IRAS source. Strong SiO emission observed at ~50" east of the IRAS source presents evidence for a dynamical interaction between a part of the eastern blue lobe and the dense ambient gas condensation, however, such an interaction is unlikely to be responsible for producing the quadrupolar morphology. The SiO emission in the northeastern red lobe shows the spatial and velocity structure similar to those of the CO outflow, suggesting that the SiO emission comes from the molecular outflow in the northeastern red lobe itself. The observed velocity structure is reproduced by a simple spatio-kinematic model of bow shock with a shock velocity of 19-24 km/s inclined by 30-45 deg from the plane of the sky. This implies that the northeastern red lobe is independent of the eastern blue lobe and that the quadrupolar structure is due to two separate bipolar outflows. The SiO emission observed in the western red lobe has a broad pedestal shape with low intensity. Unlike the SiO emission in the northeastern red lobe, the spatial extent of the SiO emission in the western red lobe is restricted to its central region. The spatial and velocity structures and the line profiles suggest that three different types of the SiO emission are observed in this outflow; the SiO emission arises from the interface between the outflowing gas and the dense ambient gas clump, the SiO emission coming from the outflow lobe itself, and the broad SiO emission with low intensity observed at the central region of the outflow lobe.

We have measured the ortho-to-para ratio of ammonia in the blueshifted gas of the L1157 outflow by observing the six metastable inversion lines from (J, K) = (1, 1) to (6, 6). The highly excited (5, 5) and (6, 6) lines were first detected in the low-mass star forming regions. The rotational temperature derived from the ratio of four transition lines from (3, 3) to (6, 6) is 130-140 K, suggesting that the blueshifted gas is heated by a factor of ~10 as compared to the quiescent gas. The ortho-to-para ratio of the NH3 molecules in the blueshifted gas is estimated to be 1.3--1.7, which is higher than the statistical equilibrium value. This ratio provides us with evidence that the NH3 molecules have been evaporated from dust grains with the formation temperature between 18 and 25 K. It is most likely that the NH3 molecules on dust grains have been released into the gas phase through the passage of strong shock waves produced by the outflow. Such a scenario is supported by the fact that the ammonia abundance in the blueshifted gas is enhanced by a factor of ~5 with respect to the dense quiescent gas.