Search

Leveraging the high resolution, high sensitivity, and wide frequency coverage of the Atacama Large Millimeter/submillimeter Array (ALMA), the QUARKS survey, standing for "Querying Underlying mechanisms of massive star formation with ALMA-Resolved gas Kinematics and Structures", is observing 139 massive star-forming clumps at ALMA Band 6 ($\lambda\sim$ 1.3 mm). This paper introduces the Atacama Compact Array (ACA) 7-m data. Combining multi-wavelength data, we provide the first edition of QUARKS atlas, offering insights into the multiscale and multiphase interstellar medium in high-mass star formation. The ACA 1.3 mm catalog includes 207 continuum sources that are called ACA sources. Their gas kinetic temperatures are estimated using three formaldehyde (H$_2$CO) transitions with a non-LTE radiation transfer model, and the mass and density are derived from a dust emission model. The ACA sources are massive (16-84 percentile values of 6-160 $M_{\odot}$), gravity-dominated ($M\propto R^{1.1}$) fragments within massive clumps, with supersonic turbulence ($\mathcal{M}>1$) and embedded star-forming protoclusters. We find a linear correlation between the masses of the fragments and the massive clumps, with a ratio of 6% between the two. When considering the fragments as representative of dense gas, the ratio indicates a dense gas fraction (DGF) of 6%, although with a wide scatter ranging from 1% to 10%. If we consider the QUARKS massive clumps to be what is observed at various scales, then the size-independent DGF indicates a self-similar fragmentation or collapsing mode in protocluster formation. With the ACA data over four orders of magnitude of luminosity-to-mass ratio ($L/M$), we find that the DGF increases significantly with $L/M$, which indicates clump evolutionary stage. We observed a limited fragmentation at the subclump scale, which can be explained by dynamic global collapse process., Comment: 24 pages, 7 figures. Accepted for publication in Research in Astronomy and Astrophysics. QUARKS atlas link: https://drive.google.com/file/d/1KTqXxCDduYepvLd9kIvZVSSytK48OmfL/view?usp=sharing

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet., Comment: 16 pages, 7 figures. accepted for publication in ApJ

We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, respectively. Several tens of dendrogram structures have been extracted in the position-position-velocity space of H13CO+, which include 21 small-scale leaves and 20 larger-scale branches. Overall, their gas motions are supersonic but they exhibit the interesting behavior where leaves tend to be less dynamically supersonic than the branches. For the larger-scale, branch structures, the observed velocity-size relation (i.e., velocity variation/dispersion versus size) are seen to follow the Larson scaling exponent while the smaller-scale, leaf structures show a systematic deviation and display a steeper slope. We argue that the origin of the observed kinematics of the branch structures is likely to be a combination of turbulence and gravity-driven ordered gas flows. In comparison, gravity-driven chaotic gas motion is likely at the level of small-scale leaf structures. The results presented in our previous paper and this current follow-up study suggest that the main driving mechanism for mass accretion/inflow observed in G34 varies at different spatial scales. We therefore conclude that a scale-dependent combined effect of turbulence and gravity is essential to explain the star-formation processes in G34., Comment: 11 pages, 6 figures, and 1 table. To appear in MNRAS

Protostellar outflows and jets play a vital role in star formation as they carry away excess angular momentum from the inner disk surface, allowing the material to be transferred toward the central protostar. Theoretically, low velocity and poorly collimated outflows appear from the beginning of the collapse, at the first hydrostatic core (FHSC) stage. With growing protostellar core mass, high-density jets are launched which entrain an outflow from the infalling envelope. Until now, molecular jets have been observed at high velocity ($\gtrsim$ 100 km/s) in early Class\,0 protostars. We, for the first time, detect a dense molecular jet in SiO emission with small-velocity ($\sim$ 4.2 km\,s$^{-1}$, deprojected $\sim$ 24 km\,s$^{-1}$) from source G208.89-20.04Walma (hereafter, G208Walma) using ALMA Band\,6 observations. This object has some characteristics of FHSCs, such as a small outflow/jet velocity, extended 1.3\,mm continuum emission, and N$_2$D$^+$ line emission. Additional characteristics, however, are typical of early protostars: collimated outflow and SiO jet. The full extent of the outflow corresponds to a dynamical time scale of $\sim$ 930$^{+200}_{-100}$ years. The spectral energy distribution also suggests a very young source having an upper limit of T$_{bol}$ $\sim$ 31 K and L$_{bol}$ $\sim$ 0.8 L$_\sun$. We conclude that G208Walma is likely in the transition phase from FHSC to protostar, and the molecular jet has been launched within a few hundred years of initial collapse. Therefore, G208Walma may be the earliest object discovered in the protostellar phase with a molecular jet., Comment: 16 pages, 7 figures. accepted for publication in ApJ

We present a comprehensive study of the gas kinematics associated with density structures at different spatial scales in the filamentary infrared dark cloud, G034.43+00.24 (G34). This study makes use of the H13CO+ (1-0) molecular line data from the ALMA Three-millimeter Observations of Massive Star-forming regions (ATOMS) survey, which has spatial and velocity resolution of 0.04 pc and 0.2 km/s, respectively. Several tens of dendrogram structures have been extracted in the position-position-velocity space of H13CO+, which include 21 small-scale leaves and 20 larger-scale branches. Overall, their gas motions are supersonic but they exhibit the interesting behavior where leaves tend to be less dynamically supersonic than the branches. For the larger-scale, branch structures, the observed velocity-size relation (i.e., velocity variation/dispersion versus size) are seen to follow the Larson scaling exponent while the smaller-scale, leaf structures show a systematic deviation and display a steeper slope. We argue that the origin of the observed kinematics of the branch structures is likely to be a combination of turbulence and gravity-driven ordered gas flows. In comparison, gravity-driven chaotic gas motion is likely at the level of small-scale leaf structures. The results presented in our previous paper and this current follow-up study suggest that the main driving mechanism for mass accretion/inflow observed in G34 varies at different spatial scales. We therefore conclude that a scale-dependent combined effect of turbulence and gravity is essential to explain the star-formation processes in G34., Comment: 11 pages, 6 figures, and 1 table. To appear in MNRAS

We present velocity-resolved Stratospheric Observatory for Infrared Astronomy (SOFIA)/upgrade German REceiver for Astronomy at Terahertz Frequencies observations of [O I] and [C II] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [O I] emission toward only the protostar and [C II] emission toward the protostar and the redshifted outflow. The [O I] emission has a width of similar to 100 km s(-1) only in the blueshifted velocity, suggesting an origin in shocked gas. The [C ii] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing redshifted [O I] emission. Fitting the [O I] line profile with two Gaussian components, we find one component at the source velocity with a width of similar to 20 km s(-1) and another extremely broad component at -30 km s(-1) with a width of 87.5 km s(-1), the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [O I], high-J CO, and H2O lines in the cavity shocks indicate that [O I] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of similar to 1.5 x 10(-4). Attributing the extremely broad [O I] component to atomic winds, we estimate an intrinsic mass-loss rate of (1.3 +/- 0.8) x 10(-6) M-circle dot yr(-1). The intrinsic mass-loss rates derived from low-J CO, [O i], and H i are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows.

We present velocity-resolved Stratospheric Observatory for Infrared Astronomy (SOFIA)/upgrade German REceiver for Astronomy at Terahertz Frequencies observations of [O I] and [C II] lines toward a Class I protostar, L1551 IRS 5, and its outflows. The SOFIA observations detect [O I] emission toward only the protostar and [C II] emission toward the protostar and the redshifted outflow. The [O I] emission has a width of similar to 100 km s(-1) only in the blueshifted velocity, suggesting an origin in shocked gas. The [C ii] lines are narrow, consistent with an origin in a photodissociation region. Differential dust extinction from the envelope due to the inclination of the outflows is the most likely cause of the missing redshifted [O I] emission. Fitting the [O I] line profile with two Gaussian components, we find one component at the source velocity with a width of similar to 20 km s(-1) and another extremely broad component at -30 km s(-1) with a width of 87.5 km s(-1), the latter of which has not been seen in L1551 IRS 5. The kinematics of these two components resemble cavity shocks in molecular outflows and spot shocks in jets. Radiative transfer calculations of the [O I], high-J CO, and H2O lines in the cavity shocks indicate that [O I] dominates the oxygen budget, making up more than 70% of the total gaseous oxygen abundance and suggesting [O]/[H] of similar to 1.5 x 10(-4). Attributing the extremely broad [O I] component to atomic winds, we estimate an intrinsic mass-loss rate of (1.3 +/- 0.8) x 10(-6) M-circle dot yr(-1). The intrinsic mass-loss rates derived from low-J CO, [O i], and H i are similar, supporting the model of momentum-conserving outflows, where the atomic wind carries most momentum and drives the molecular outflows.

Gas at high Galactic latitude is a relatively little-noticed component of the interstellar medium. In an effort to address this, forty-one Planck Galactic Cold Clumps at high Galactic latitude (HGal; $|b|>25^{\circ}$) were observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=1-0 lines, using the Purple Mountain Observatory 13.7-m telescope. $^{12}$CO (1-0) and $^{13}$CO (1-0) emission was detected in all clumps while C$^{18}$O (1-0) emission was only seen in sixteen clumps. The highest and average latitudes are $71.4^{\circ}$ and $37.8^{\circ}$, respectively. Fifty-one velocity components were obtained and then each was identified as a single clump. Thirty-three clumps were further mapped at 1$^\prime$ resolution and 54 dense cores were extracted. Among dense cores, the average excitation temperature $T_{\mathrm{ex}}$ of $^{12}$CO is 10.3 K. The average line widths of thermal and non-thermal velocity dispersions are $0.19$ km s$^{-1}$ and $0.46$ km s$^{-1}$ respectively, suggesting that these cores are dominated by turbulence. Distances of the HGal clumps given by Gaia dust reddening are about $120-360$ pc. The ratio of $X_{13}$/$X_{18}$ is significantly higher than that in the solar neighbourhood, implying that HGal gas has a different star formation history compared to the gas in the Galactic disk. HGal cores with sizes from $0.01-0.1$ pc show no notable Larson's relation and the turbulence remains supersonic down to a scale of slightly below $0.1$ pc. None of the HGal cores which bear masses from 0.01-1 $M_{\odot}$ are gravitationally bound and all appear to be confined by outer pressure., Comment: 35 pages, 13 figures

We have identified 453 compact dense cores in 3 mm continuum emission maps in the ATOMS (ALMA Three-millimeter Observations of Massive Star-forming regions) survey, and compiled three catalogues of high-mass star forming cores. One catalogue, referred to as H/UC-HII catalogue, includes 89 cores that enshroud hyper/ultra compact (H/UC) HII regions as characterized by associated compact H40alpha emission. A second catalogue, referred to as pure s-cHMC, includes 32 candidate Hot Molecular Cores (HMCs) showing rich spectra (N>20lines) of complex organic molecules (COMs) but not associated with H/UC-HII regions. The third catalogue, referred to as pure w-cHMC, includes 58 candidate HMCs with relatively low levels of COM richness and not associated with H/UC-HII regions. These three catalogues of dense cores provide an important foundation for future studies of the early stages of high-mass star formation across the Milky Way. We also find that nearly half of H/UC-HII cores are candidate HMCs. From the number counts of COM-containing and H/UC-HII cores, we suggest that the duration of high-mass protostellar cores showing chemically rich features is at least comparable to the lifetime of H/UC-HII regions. For cores in the H/UC-HII catalogue, the width of the H40alpha line increases as the core size decreases, suggesting that the non-thermal dynamical and/or pressure line-broadening mechanisms dominate on the smaller scales of the H/UC-HII cores., Comment: 17 pages, five tables, and 11 figures. Accepted for publication at MNRAS

Gas at high Galactic latitude is a relatively little noticed component of the interstellar medium. In an effort to address this, 41 Planck Galactic Cold Clumps at high Galactic latitude (HGal; divide b divide > 25 degrees) were observed in (CO)-C-12, (CO)-C-13, and (CO)-O-18 J = 1-0 lines, using the Purple Mountain Observatory 13.7 m telescope. (CO)-C-12 (1-0) and (CO)-C-13 (1-0) emission was detected in all clumps, while (CO)-O-18 (1-0) emission was only seen in 16 clumps. The highest and average latitudes are 71.degrees 4 and 37.degrees 8, respectively. Fifty-one velocity components were obtained, and then each was identified as a single clump. Thirty-three clumps were further mapped at 1 ' resolution, and 54 dense cores were extracted. Among dense cores, the average excitation temperature T (ex) of (CO)-C-12 is 10.3 K. The average line widths of thermal and nonthermal velocity dispersions are 0.19 and 0.46 km s(-1), respectively, suggesting that these cores are dominated by turbulence. Distances of the HGal clumps given by Gaia dust reddening are about 120-360 pc. The ratio of X (13)/X (18) is significantly higher than that in the solar neighborhood, implying that HGal gas has a different star formation history compared to the gas in the Galactic disk. HGal cores with sizes from 0.01 to 0.1 pc show no notable Larson's relation, and the turbulence remains supersonic down to a scale of slightly below 0.1 pc. None of the HGal cores that bear masses from 0.01 to 1 M (circle dot) are gravitationally bound, and all appear to be confined by outer pressure.

We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the J = 1-0 lines of (CO)-C-13, (CO)-O-18, HCN, HCO+, and N2H+ and the J= 2-1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission maps, are presented and compared to previous work. The CS line traces dust column density over more than one order of magnitude, and the N2H+ line best traces the highest column density regime (log(NH2) > 22.8). Line luminosities, integrated over the cloud, are compared to those seen in other galaxies. The HCO+-to-HCN luminosity ratio in the Orion A cloud is similar to that of starburst galaxies, while that in the Ophiuchus cloud is in between those of active galactic nuclei and starburst galaxies.

We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in (CO)-C-13, (CO)-O-18, HCO+ J = 1-0, and CS J = 2-1. First, the scaling relations derived from PCA of the (CO)-C-13 maps show that the velocity difference (delta v) for a given spatial scale (L) is the highest in the integral-shaped filament (ISF) and L1688, where the most active star formation occurs in the two clouds. The delta v increases with the number density and total bolometric luminosity of the protostars in the subregions. Second, in the ISF and L1688 regions, the delta v of (CO)-O-18, HCO+, and CS are generally higher than that of (CO)-C-13, which implies that the dense gas is more turbulent than the diffuse gas in the star-forming regions; stars form in dense gas, and dynamical activities associated with star formation, such as jets and outflows, can provide energy into the surrounding gas to enhance turbulent motions.

Gas at high Galactic latitude is a relatively little noticed component of the interstellar medium. In an effort to address this, 41 Planck Galactic Cold Clumps at high Galactic latitude (HGal; divide b divide > 25 degrees) were observed in (CO)-C-12, (CO)-C-13, and (CO)-O-18 J = 1-0 lines, using the Purple Mountain Observatory 13.7 m telescope. (CO)-C-12 (1-0) and (CO)-C-13 (1-0) emission was detected in all clumps, while (CO)-O-18 (1-0) emission was only seen in 16 clumps. The highest and average latitudes are 71.degrees 4 and 37.degrees 8, respectively. Fifty-one velocity components were obtained, and then each was identified as a single clump. Thirty-three clumps were further mapped at 1 ' resolution, and 54 dense cores were extracted. Among dense cores, the average excitation temperature T (ex) of (CO)-C-12 is 10.3 K. The average line widths of thermal and nonthermal velocity dispersions are 0.19 and 0.46 km s(-1), respectively, suggesting that these cores are dominated by turbulence. Distances of the HGal clumps given by Gaia dust reddening are about 120-360 pc. The ratio of X (13)/X (18) is significantly higher than that in the solar neighborhood, implying that HGal gas has a different star formation history compared to the gas in the Galactic disk. HGal cores with sizes from 0.01 to 0.1 pc show no notable Larson's relation, and the turbulence remains supersonic down to a scale of slightly below 0.1 pc. None of the HGal cores that bear masses from 0.01 to 1 M (circle dot) are gravitationally bound, and all appear to be confined by outer pressure.

We investigate the effect of star formation on turbulence in the Orion A and Ophiuchus clouds using principal component analysis (PCA). We measure the properties of turbulence by applying PCA on the spectral maps in (CO)-C-13, (CO)-O-18, HCO+ J = 1-0, and CS J = 2-1. First, the scaling relations derived from PCA of the (CO)-C-13 maps show that the velocity difference (delta v) for a given spatial scale (L) is the highest in the integral-shaped filament (ISF) and L1688, where the most active star formation occurs in the two clouds. The delta v increases with the number density and total bolometric luminosity of the protostars in the subregions. Second, in the ISF and L1688 regions, the delta v of (CO)-O-18, HCO+, and CS are generally higher than that of (CO)-C-13, which implies that the dense gas is more turbulent than the diffuse gas in the star-forming regions; stars form in dense gas, and dynamical activities associated with star formation, such as jets and outflows, can provide energy into the surrounding gas to enhance turbulent motions.

We have used the Taeduk Radio Astronomy Observatory to observe the Orion A and Ophiuchus clouds in the J = 1-0 lines of (CO)-C-13, (CO)-O-18, HCN, HCO+, and N2H+ and the J= 2-1 line of CS. The fully sampled maps with uniform noise levels are used to create moment maps. The variations of the line intensity and velocity dispersion with total column density, derived from dust emission maps, are presented and compared to previous work. The CS line traces dust column density over more than one order of magnitude, and the N2H+ line best traces the highest column density regime (log(NH2) > 22.8). Line luminosities, integrated over the cloud, are compared to those seen in other galaxies. The HCO+-to-HCN luminosity ratio in the Orion A cloud is similar to that of starburst galaxies, while that in the Ophiuchus cloud is in between those of active galactic nuclei and starburst galaxies.

Gas at high Galactic latitude is a relatively little-noticed component of the interstellar medium. In an effort to address this, forty-one Planck Galactic Cold Clumps at high Galactic latitude (HGal; $|b|>25^{\circ}$) were observed in $^{12}$CO, $^{13}$CO and C$^{18}$O J=1-0 lines, using the Purple Mountain Observatory 13.7-m telescope. $^{12}$CO (1-0) and $^{13}$CO (1-0) emission was detected in all clumps while C$^{18}$O (1-0) emission was only seen in sixteen clumps. The highest and average latitudes are $71.4^{\circ}$ and $37.8^{\circ}$, respectively. Fifty-one velocity components were obtained and then each was identified as a single clump. Thirty-three clumps were further mapped at 1$^\prime$ resolution and 54 dense cores were extracted. Among dense cores, the average excitation temperature $T_{\mathrm{ex}}$ of $^{12}$CO is 10.3 K. The average line widths of thermal and non-thermal velocity dispersions are $0.19$ km s$^{-1}$ and $0.46$ km s$^{-1}$ respectively, suggesting that these cores are dominated by turbulence. Distances of the HGal clumps given by Gaia dust reddening are about $120-360$ pc. The ratio of $X_{13}$/$X_{18}$ is significantly higher than that in the solar neighbourhood, implying that HGal gas has a different star formation history compared to the gas in the Galactic disk. HGal cores with sizes from $0.01-0.1$ pc show no notable Larson's relation and the turbulence remains supersonic down to a scale of slightly below $0.1$ pc. None of the HGal cores which bear masses from 0.01-1 $M_{\odot}$ are gravitationally bound and all appear to be confined by outer pressure., Comment: 35 pages, 13 figures

We have identified 453 compact dense cores in 3 mm continuum emission maps in the ATOMS (ALMA Three-millimeter Observations of Massive Star-forming regions) survey, and compiled three catalogues of high-mass star forming cores. One catalogue, referred to as H/UC-HII catalogue, includes 89 cores that enshroud hyper/ultra compact (H/UC) HII regions as characterized by associated compact H40alpha emission. A second catalogue, referred to as pure s-cHMC, includes 32 candidate Hot Molecular Cores (HMCs) showing rich spectra (N>20lines) of complex organic molecules (COMs) but not associated with H/UC-HII regions. The third catalogue, referred to as pure w-cHMC, includes 58 candidate HMCs with relatively low levels of COM richness and not associated with H/UC-HII regions. These three catalogues of dense cores provide an important foundation for future studies of the early stages of high-mass star formation across the Milky Way. We also find that nearly half of H/UC-HII cores are candidate HMCs. From the number counts of COM-containing and H/UC-HII cores, we suggest that the duration of high-mass protostellar cores showing chemically rich features is at least comparable to the lifetime of H/UC-HII regions. For cores in the H/UC-HII catalogue, the width of the H40alpha line increases as the core size decreases, suggesting that the non-thermal dynamical and/or pressure line-broadening mechanisms dominate on the smaller scales of the H/UC-HII cores., Comment: 17 pages, five tables, and 11 figures. Accepted for publication at MNRAS

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation., Comment: published on MNRAS

Using wide-field $^{13}$CO ($J=1-0$) data taken with the Nobeyama 45-m telescope, we investigate cloud structures of the infrared dark cloud complex in M17 with SCIMES. In total, we identified 118 clouds that contain 11 large clouds with radii larger than 1 pc. The clouds are mainly distributed in the two representative velocity ranges of 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$. By comparing with the ATLASGAL catalog, we found that the majority of the $^{13}$CO clouds with 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$ are likely located at distances of 2 kpc (Sagittarius arm) and 3 kpc (Scutum arm), respectively. Analyzing the spatial configuration of the identified clouds and their velocity structures, we attempt to reveal the origin of the cloud structure in this region. Here we discuss three possibilities: (1) overlapping with different velocities, (2) cloud oscillation, and (3) cloud-cloud collision. From the position-velocity diagrams, we found spatially-extended faint emission between $\sim$ 20 km s$^{-1}$ and $\sim$ 35 km s$^{-1}$, which is mainly distributed in the spatially-overlapped areas of the clouds. We also found that in some areas where clouds with different velocities overlapped, the magnetic field orientation changes abruptly. The distribution of the diffuse emission in the position-position-velocity space and the bending magnetic fields appear to favor the cloud-cloud collision scenario compared to other scenarios. In the cloud-cloud collision scenario, we propose that two $\sim$35 km s$^{-1}$ foreground clouds are colliding with clouds at $\sim$20 km s$^{-1}$ with a relative velocity of 15 km s$^{-1}$. These clouds may be substructures of two larger clouds having velocities of $\sim$ 35 km s$^{-1}$ ($\gtrsim 10^3 $ M$_{\odot}$) and $\sim$ 20 km s$^{-1}$ ($\gtrsim 10^4 $ M$_{\odot}$), respectively., Comment: 34 pages, 13 figures, Accepted for publication in PASJ

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales., Comment: Published on MNRAS. The full tables are included in Tables.pdf or Tables.tex files, which can be downloaded from source files

Planck Galactic Cold Clumps (PGCCs) are contemplated to be the ideal targets to probe the early phases of star formation. We have conducted a survey of 72 young dense cores inside PGCCs in the Orion complex with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3\,mm (band 6) using three different configurations (resolutions $\sim$ 0$\farcs$35, 1$\farcs$0, and 7$\farcs$0) to statistically investigate their evolutionary stages and sub-structures. We have obtained images of the 1.3\,mm continuum and molecular line emission ($^{12}$CO, and SiO) at an angular resolution of $\sim$ 0$\farcs$35 ($\sim$ 140\,au) with the combined arrays. We find 70 substructures within 48 detected dense cores with median dust-mass $\sim$ 0.093\,M$_{\sun}$ and deconvolved size $\sim$ 0$\farcs$27. Dense substructures are clearly detected within the central 1000\,au of four candidate prestellar cores. The sizes and masses of the substructures in continuum emission are found to be significantly reduced with protostellar evolution from Class\,0 to Class\,I. We also study the evolutionary change in the outflow characteristics through the course of protostellar mass accretion. A total of 37 sources exhibit CO outflows, and 20 ($>$50\%) show high-velocity jets in SiO. The CO velocity-extents ($\Delta$Vs) span from 4 to 110 km/s with outflow cavity opening angle width at 400\,au ranging from $[\Theta_{obs}]_{400}$ $\sim$ 0$\farcs$6 to 3$\farcs$9, which corresponds to 33$\fdg$4$-$125$\fdg$7. For the majority of the outflow sources, the $\Delta$Vs show a positive correlation with $[\Theta_{obs}]_{400}$, suggesting that as protostars undergo gravitational collapse, the cavity opening of a protostellar outflow widens and the protostars possibly generate more energetic outflows., Comment: 43 pages, 25 figures (14 in the main text, 9 in the appendix), accepted for publication in the Astrophysical Journal Supplement Series (ApJS)

We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($\lambda$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $\lambda$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc., Comment: 55 pages, 12 Figures, 8 Tables, accepted for publication in ApJS

The "ATOMS," standing for {\it ALMA Three-millimeter Observations of Massive Star-forming regions}, survey has observed 146 active star forming regions with ALMA Band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the "ATOMS" survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS $J = 2-1$, HCO$^+$ $J = 1-0$ and HCN $J = 1-0$, are found to show extended gas emission in low density regions within the clump; less than 25\% of their emission is from dense cores. SO, CH$_3$OH, H$^{13}$CN and HC$_3$N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over $\sim$1 pc), which may be caused by slow shocks from large--scale colliding flows or H{\sc ii} regions. Stellar feedback from an expanding H{\sc ii} region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analyzed with other survey data, e.g., "MALT90", "Orion B", "EMPIRE", "ALMA\_IMF", and "ALMAGAL", to deepen our understandings of "dense gas" star formation scaling relations and massive proto-cluster formation., Comment: published on MNRAS

We report studies of the relationships between the total bolometric luminosity ($L_{\rm bol}$ or $L_{\rm TIR}$) and the molecular line luminosities of $J=1-0$ transitions of H$^{13}$CN, H$^{13}$CO$^+$, HCN, and HCO$^+$ with data obtained from ACA observations in the "ATOMS" survey of 146 active Galactic star forming regions. The correlations between $L_{\rm bol}$ and molecular line luminosities $L'_{\rm mol}$ of the four transitions all appear to be approximately linear. Line emission of isotopologues shows as large scatters in $L_{\rm bol}$-$L'_{\rm mol}$ relations as their main line emission. The log($L_{\rm bol}$/$L'_{\rm mol}$) for different molecular line tracers have similar distributions. The $L_{\rm bol}$-to-$L'_{\rm mol}$ ratios do not change with galactocentric distances ($R_{\rm GC}$) and clump masses ($M_{\rm clump}$). The molecular line luminosity ratios (HCN-to-HCO$^+$, H$^{13}$CN-to-H$^{13}$CO$^+$, HCN-to-H$^{13}$CN and HCO$^+$-to-H$^{13}$CO$^+$) all appear constant against $L_{\rm bol}$, dust temperature ($T_{\rm d}$), $M_{\rm clump}$ and $R_{\rm GC}$. Our studies suggest that both the main lines and isotopologue lines are good tracers of the total masses of dense gas in Galactic molecular clumps. The large optical depths of main lines do not affect the interpretation of the slopes in star formation relations. We find that the mean star formation efficiency (SFE) of massive Galactic clumps in the "ATOMS" survey is reasonably consistent with other measures of the SFE for dense gas, even those using very different tracers or examining very different spatial scales., Comment: Published on MNRAS. The full tables are included in Tables.pdf or Tables.tex files, which can be downloaded from source files

Using wide-field $^{13}$CO ($J=1-0$) data taken with the Nobeyama 45-m telescope, we investigate cloud structures of the infrared dark cloud complex in M17 with SCIMES. In total, we identified 118 clouds that contain 11 large clouds with radii larger than 1 pc. The clouds are mainly distributed in the two representative velocity ranges of 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$. By comparing with the ATLASGAL catalog, we found that the majority of the $^{13}$CO clouds with 10 $-$ 20 km s$^{-1}$ and 30 $-$ 40 km s$^{-1}$ are likely located at distances of 2 kpc (Sagittarius arm) and 3 kpc (Scutum arm), respectively. Analyzing the spatial configuration of the identified clouds and their velocity structures, we attempt to reveal the origin of the cloud structure in this region. Here we discuss three possibilities: (1) overlapping with different velocities, (2) cloud oscillation, and (3) cloud-cloud collision. From the position-velocity diagrams, we found spatially-extended faint emission between $\sim$ 20 km s$^{-1}$ and $\sim$ 35 km s$^{-1}$, which is mainly distributed in the spatially-overlapped areas of the clouds. We also found that in some areas where clouds with different velocities overlapped, the magnetic field orientation changes abruptly. The distribution of the diffuse emission in the position-position-velocity space and the bending magnetic fields appear to favor the cloud-cloud collision scenario compared to other scenarios. In the cloud-cloud collision scenario, we propose that two $\sim$35 km s$^{-1}$ foreground clouds are colliding with clouds at $\sim$20 km s$^{-1}$ with a relative velocity of 15 km s$^{-1}$. These clouds may be substructures of two larger clouds having velocities of $\sim$ 35 km s$^{-1}$ ($\gtrsim 10^3 $ M$_{\odot}$) and $\sim$ 20 km s$^{-1}$ ($\gtrsim 10^4 $ M$_{\odot}$), respectively., Comment: 34 pages, 13 figures, Accepted for publication in PASJ

The collapse of a protostellar envelope results in the growth of a protostar and the development of a protoplanetary disk, playing a critical role during the early stages of star formation. Characterizing the gas infall in the envelope constrains the dynamical models of star formation. We present unambiguous signatures of infall, probed by optically thick molecular lines, toward an isolated embedded protostar, BHR 71 IRS1. The three-dimensional radiative transfer calculations indicate that a slowly rotating infalling envelope model following the "inside-out" collapse reproduces the observations of both lines, as well as the low-velocity emission of the HCN line. The envelope has a model-derived ag lines, where outflows or a Keplerian disk may contribute. The ALMA observations serendipitously discover the emission of complex organic molecules (COMs) concentrated within a radius of 100 au, indicating that BHR 71 IRS1 harbors a hot corino. Eight species of COMs are identified, including CH3OH and CH3OCHO, along with H2CS, SO2 and HCN v(2) = 1. The emission of methyl formate and C-13-methanol shows a clear velocity gradient within a radius of 50 au, hinting at an unresolved Keplerian rotating disk.

The collapse of a protostellar envelope results in the growth of a protostar and the development of a protoplanetary disk, playing a critical role during the early stages of star formation. Characterizing the gas infall in the envelope constrains the dynamical models of star formation. We present unambiguous signatures of infall, probed by optically thick molecular lines, toward an isolated embedded protostar, BHR 71 IRS1. The three-dimensional radiative transfer calculations indicate that a slowly rotating infalling envelope model following the "inside-out" collapse reproduces the observations of both lines, as well as the low-velocity emission of the HCN line. The envelope has a model-derived ag lines, where outflows or a Keplerian disk may contribute. The ALMA observations serendipitously discover the emission of complex organic molecules (COMs) concentrated within a radius of 100 au, indicating that BHR 71 IRS1 harbors a hot corino. Eight species of COMs are identified, including CH3OH and CH3OCHO, along with H2CS, SO2 and HCN v(2) = 1. The emission of methyl formate and C-13-methanol shows a clear velocity gradient within a radius of 50 au, hinting at an unresolved Keplerian rotating disk.

The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns., Comment: Science white paper submitted to the Astro2020 Decadal Survey

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes., Comment: Submitted to the US Astronomy and Astrophysics Decadal Survey (Astro2020)

The gaseous mass of protoplanetary disks is a fundamental quantity in planet formation. The presence of gas is necessary to assemble planetesimals, it determines timescales of giant planet birth, and it is an unknown factor for a wide range of properties of planet formation, from chemical abundances (X/H) to the mass efficiency of planet formation. The gas mass obtained from traditional tracers, such as dust thermal continuum and CO isotopologues, are now known to have significant (1 - 2 orders of magnitude) discrepancies. Emission from the isotopologue of H2, hydrogen deuteride (HD), offers an alternative measurement of the disk gas mass. Of all of the regions of the spectrum, the far-infrared stands out in that orders of magnitude gains in sensitivity can be gleaned by cooling a large aperture telescope to 8 K. Such a facility can open up a vast new area of the spectrum to exploration. One of the primary benefits of this far-infrared revolution would be the ability to survey hundreds of planet-forming disks in HD emission to derive their gaseous masses. For the first time, we will have statistics on the gas mass as a function of evolution, tracing birth to dispersal as a function of stellar spectral type. These measurements have broad implications for our understanding of the time scale during which gas is available to form giant planets, the dynamical evolution of the seeds of terrestrial worlds, and the resulting chemical composition of pre-planetary embryos carrying the elements needed for life. Measurements of the ground-state line of HD requires a space-based observatory operating in the far-infrared at 112 microns., Comment: Science white paper submitted to the Astro2020 Decadal Survey

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes., Comment: Submitted to the US Astronomy and Astrophysics Decadal Survey (Astro2020)

Our ALMA observations of HCO+ and HCN show such redshifted absorption toward an isolated core, BHR 71. Both lines show a similar redshifted absorption profile. We also found emissions of complex organic molecules (COMs) around 345 GHz from a compact region centered on the continuum source, which is barely resolved with a beam of 0″27, corresponding to ∼50 AU.

Magnetic field plays a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained owing 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 μm with SCUBA-2/POL-2 at JCMT for the first time. 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 ~50 μG for F M is obtained using the Davis–Chandrasekhar–Fermi method. Based on 13CO (1–0) line observations, we suggest a formation scenario of F M due to large-scale (~10 pc) cloud–cloud collision. Using additional NH3 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 (~200 M ⊙), 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.

Magnetic field plays a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained owing 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 μm with SCUBA-2/POL-2 at JCMT for the first time. 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 ~50 μG for F M is obtained using the Davis–Chandrasekhar–Fermi method. Based on 13CO (1–0) line observations, we suggest a formation scenario of F M due to large-scale (~10 pc) cloud–cloud collision. Using additional NH3 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 (~200 M ⊙), 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.

Magnetic field plays a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained owing 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 μm with SCUBA-2/POL-2 at JCMT for the first time. 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 ~50 μG for F M is obtained using the Davis–Chandrasekhar–Fermi method. Based on 13CO (1–0) line observations, we suggest a formation scenario of F M due to large-scale (~10 pc) cloud–cloud collision. Using additional NH3 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 (~200 M ⊙), 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.

The low dust temperatures (< 14 K) of Planck Galactic cold clumps (PGCCs) make them ideal targets to probe the initial conditions and very early phase of star formation. TOP-SCOPE is a joint survey program targeting similar to 2000 PGCCs in J = 1-0 transitions of CO isotopologues and similar to 1000 PGCCs in 850 mu m continuum emission. The objective of the TOP-SCOPE survey and the joint surveys (SMT 10 m, KVN 21 m, and NRO 45 m) is to statistically study the initial conditions occurring during star formation and the evolution of molecular clouds, across a wide range of environments. The observations, data analysis, and example science cases for these surveys are introduced with an exemplar source, PGCC G26.53+0.17 (G26), which is a filamentary infrared dark cloud (IRDC). The total mass, length, and mean line mass (M/L) of the G26 filament are similar to 6200 M-circle dot, similar to 12 pc, and similar to 500 M-circle dot pc(-1), respectively. Ten massive clumps, including eight starless ones, are found along the filament. The most massive clump as a whole may still be in global collapse, while its denser part seems to be undergoing expansion owing to outflow feedback. The fragmentation in the G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, nonmagnetized, and turbulent supported cylinder. A bimodal behavior in dust emissivity spectral index (beta) distribution is found in G26, suggesting grain growth along the filament. The G26 filament may be formed owing to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.

Magnetic field plays a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained owing 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 mu m with SCUBA-2/POL-2 at JCMT for the first time. 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 similar to 50 mu G for F-M is obtained using the Davis-Chandrasekhar-Fermi method. Based on (CO)-C-13 (1-0) line observations, we suggest a formation scenario of F-M due to large-scale (similar to 10 pc) cloud-cloud collision. Using additional NH3 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 (similar to 200 M-circle dot, 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.

The low dust temperatures (< 14 K) of Planck Galactic cold clumps (PGCCs) make them ideal targets to probe the initial conditions and very early phase of star formation. TOP-SCOPE is a joint survey program targeting similar to 2000 PGCCs in J = 1-0 transitions of CO isotopologues and similar to 1000 PGCCs in 850 mu m continuum emission. The objective of the TOP-SCOPE survey and the joint surveys (SMT 10 m, KVN 21 m, and NRO 45 m) is to statistically study the initial conditions occurring during star formation and the evolution of molecular clouds, across a wide range of environments. The observations, data analysis, and example science cases for these surveys are introduced with an exemplar source, PGCC G26.53+0.17 (G26), which is a filamentary infrared dark cloud (IRDC). The total mass, length, and mean line mass (M/L) of the G26 filament are similar to 6200 M-circle dot, similar to 12 pc, and similar to 500 M-circle dot pc(-1), respectively. Ten massive clumps, including eight starless ones, are found along the filament. The most massive clump as a whole may still be in global collapse, while its denser part seems to be undergoing expansion owing to outflow feedback. The fragmentation in the G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, nonmagnetized, and turbulent supported cylinder. A bimodal behavior in dust emissivity spectral index (beta) distribution is found in G26, suggesting grain growth along the filament. The G26 filament may be formed owing to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.

Magnetic field plays a crucial role in shaping molecular clouds and regulating star formation, yet the complete information on the magnetic field is not well constrained owing 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 mu m with SCUBA-2/POL-2 at JCMT for the first time. 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 similar to 50 mu G for F-M is obtained using the Davis-Chandrasekhar-Fermi method. Based on (CO)-C-13 (1-0) line observations, we suggest a formation scenario of F-M due to large-scale (similar to 10 pc) cloud-cloud collision. Using additional NH3 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 (similar to 200 M-circle dot, 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.

Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s., Comment: 12 pages, including 10 figures. Accepted for publication in Astrophysical Journal. Revised following comments from community

The V380 Ori NE bipolar outflow was imaged in the SiO and CO J = 1 - 0 lines, and dense cores in L1641 were observed in the 2.0-0.89 mm continuum. The highly collimated SiO jet shows point-symmetric oscillation patterns in both position and velocity, which suggests that the jet axis is precessing and the driving source may belong to a non-coplanar binary system. By considering the position and velocity variabilities together, accurate jet parameters were derived. The protostellar system is viewed nearly edge-on, and the jet has a flow speed of 35 km/s and a precession period of 1600 years. The CO outflow length gives a dynamical timescale of 6300 years, and the protostar must be extremely young. The inferred binary separation of 6-70 au implies that this protobinary system may have been formed through the disk instability process. The continuum spectra of L1641 dense cores indicate that the emission comes from dust, and the fits with modified blackbody functions give emissivity power indices of beta = 0.3-2.2. The emissivity index shows a positive correlation with the molecular line width, but no strong correlation with bolometric luminosity or temperature. V380 Ori NE has a particularly low value of beta = 0.3, which tentatively suggests the presence of millimeter-sized dust grains. Because the dust growth takes millions of years, much longer than the protostellar age, this core may have produced large grains in the starless core stage. HH 34 MMS and HH 147 MMS also have low emissivity indices., Comment: To appear in the Astrophysical Journal Supplement Series

Gas in protostellar disks provides the raw material for giant planet formation and controls the dynamics of the planetesimal-building dust grains. Accurate gas mass measurements help map the observed properties of planet-forming disks onto the formation environments of known exoplanets. Rare isotopologues of carbon monoxide (CO) have been used as gas mass tracers for disks in the Lupus star-forming region, with an assumed interstellar CO/H$_2$ abundance ratio. Unfortunately, observations of T-Tauri disks show that CO abundance is not interstellar---a finding reproduced by models that show CO abundance decreasing both with distance from the star and as a function of time. Here we present radiative transfer simulations that assess the accuracy of CO-based disk mass measurements. We find that the combination of CO chemical depletion in the outer disk and optically thick emission from the inner disk leads observers to underestimate gas mass by more than an order of magnitude if they use the standard assumptions of interstellar CO/H$_2$ ratio and optically thin emission. Furthermore, CO abundance changes on million-year timescales, introducing an age/mass degeneracy into observations. To reach factor of a few accuracy for CO-based disk mass measurements, we suggest that observers and modelers adopt the following strategies: (1) select the low-$J$ transitions; (2) observe multiple CO isotopologues and use either intensity ratios or normalized line profiles to diagnose CO chemical depletion; and (3) use spatially resolved observations to measure the CO abundance distribution., Comment: 22 pages, 16 figures, Accepted for publication in ApJ

Turbulence is the leading candidate for angular momentum transport in protoplanetary disks and therefore influences disk lifetimes and planet formation timescales. However, the turbulent properties of protoplanetary disks are poorly constrained observationally. Recent studies have found turbulent speeds smaller than what fully-developed MRI would produce (Flaherty et al. 2015, 2017). However, existing studies assumed a constant CO/H2 ratio of 0.0001 in locations where CO is not frozen-out or photo-dissociated. Our previous studies of evolving disk chemistry indicate that CO is depleted by incorporation into complex organic molecules well inside the freeze-out radius of CO. We consider the effects of this chemical depletion on measurements of turbulence. Simon et al. (2015) suggested that the ratio of the peak line flux to the flux at line center of the CO J=3-2 transition is a reasonable diagnostic of turbulence, so we focus on that metric, while adding some analysis of the more complex effects on spatial distribution. We simulate the emission lines of CO based on chemical evolution models presented in Yu et al. (2016), and find that the peak-to-trough ratio changes as a function of time as CO is destroyed. Specifically, a CO-depleted disk with high turbulent velocity mimics the peak-to-trough ratios of a non-CO-depleted disk with lower turbulent velocity. We suggest that disk observers and modelers take into account the possibility of CO depletion when using line peak-to-trough ratios to constrain the degree of turbulence in disks. Assuming that CO/H2 = 0.0001 at all disk radii can lead to underestimates of turbulent speeds in the disk by at least 0.2 km/s., Comment: 12 pages, including 10 figures. Accepted for publication in Astrophysical Journal. Revised following comments from community