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.

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.

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.

Filament S242 is 25 pc long with massive clumps and YSO clusters concentrated in its end regions; it is considered a good example of edge collapse. We mapped this filament in the $^{12}$CO(1-0) and $^{13}$CO(1-0) lines. A large-scale velocity gradient along filament S242 has been detected; the relative velocity between the two end-clumps is $\sim$ 3 km s$^{-1}$, indicating an approaching motion between them. These signatures are consistent with the filament S242 being formed through the collapse of a single elongated entity, where an effect known as "gravitational focusing" drives the ends of the filament to collapse (edge collapse). Based on this picture, we estimate a collapse timescale of $\sim$ 4.2 Myr, which is the time needed for a finite and elongated entity evolving to the observed filament S242. For the whole filament, we find that increases in surface densities lead to increases in velocity dispersion, which can be consistently explained as the result of self-gravity. We also calculated the contribution of longitudinal collapse to the observed velocity dispersion and found it to be the dominant effect in driving the gas motion near the end-clumps. We propose that our filament S242 is formed through a two-stage collapse model, where the edge collapse of a truncated filament is followed by a stage of longitudinal accretion toward the dense end-clumps.

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

Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $\mu$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.

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

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

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

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

(Abridged) ASHES has been designed to systematically characterize the earliest stages and to constrain theories of high-mass star formation. A total of 12 massive (>500 $M_{\odot}$), cold (<15 K), 3.6-70 $\mu$m dark prestellar clump candidates, embedded in IRDCs, were carefully selected in the pilot survey to be observed with ALMA. We mosaiced each clump (~1 arcmin^2) in dust and line emission with the 12m/7m/TP arrays at 224 GHz, resulting in ~1.2" resolution (~4800 AU). As the first paper of the series, we concentrate on the dust emission to reveal the clump fragmentation. We detect 294 cores, from which 84 (29%) are categorized as protostellar based on outflow activity or 'warm core' line emission. The remaining 210 (71%) are considered prestellar core candidates. The number of detected cores is independent of the mass sensitivity range of the observations. On average, more massive clumps tend to form more cores. We find a large population of low-mass (<1 M) cores and no high-mass (>30 $M_{\odot}$) prestellar cores. The most massive prestellar core has a mass of 11 $M_{\odot}$. From the prestellar CMF, we derive a power law index of 1.17+-0.1, slightly shallower than Salpeter (1.35). We use the MST technique to characterize the separation between cores and their spatial distribution, and derive mass segregation ratios. While there is a range of core masses and separations detected in the sample, the mean separation and mass of cores are well explained by thermal fragmentation and are inconsistent with turbulent Jeans fragmentation. The core spatial distribution is well described by hierarchical subclustering rather than centrally peaked clustering. There is no conclusive evidence of mass segregation. We test several theoretical conditions, and conclude that overall, competitive accretion and global hierarchical collapse scenarios are favored over the turbulent core accretion scenario.

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 a multi-wavelength study of the Planck cold clump G181.84+0.31, which is located at the northern end of the extended filamentary structure S242. We have extracted 9 compact dense cores from the SCUBA-2 850 um map, and we have identified 18 young stellar objects (YSOs, 4 Class I and 14 Class II) based on their Spitzer, Wide-field Infrared Survey Explorer (WISE) and Two-Micron All-Sky Survey (2MASS) near- and mid-infrared colours. The dense cores and YSOs are mainly distributed along the filamentary structures of G181.84 and are well traced by HCO$^{+}$(1-0) and N$_{2}$H$^{+}$(1-0) spectral-line emission. We find signatures of sequential star formation activities in G181.84: dense cores and YSOs located in the northern and southern sub-structures are younger than those in the central region. We also detect global velocity gradients of about 0.8$\pm$0.05 km s$^{-1}$pc$^{-1}$ and 1.0$\pm$0.05 km s$^{-1}$pc$^{-1}$ along the northern and southern sub-structures, respectively, and local velocity gradients of 1.2$\pm$0.1 km s$^{-1}$pc$^{-1}$ in the central substructure. These results may be due to the fact that the global collapse of the extended filamentary structure S242 is driven by an edge effect, for which the filament edges collapse first and then further trigger star formation activities inward. We identify three substructures in G181.84 and estimate their critical masses per unit length, which are $\sim$ 101$\pm$15 M$_{\odot}$ pc$^{-1}$, 56$\pm$8 M$_{\odot}$ pc$^{-1}$ and 28$\pm$4 M$_{\odot}$ pc$^{-1}$, respectively. These values are all lower than the observed values ($\sim$ 200 M$_{\odot}$ pc$^{-1}$), suggesting that these sub-structures are gravitationally unstable.

The formation scenario of brown dwarfs is still unclear because observational studies to investigate its initial condition are quite limited. Our systematic survey of nearby low-mass star-forming regions using the Atacama Compact Array (aka Morita array) and the IRAM 30 m telescope in 1.2 mm continuum has identified a centrally concentrated starless condensation with a central H$_2$ volume density of $\sim$10$^6$ cm$^{-3}$, MC5-N, connected to a narrow (width $\sim$0.03 pc) filamentary cloud in the Taurus L1495 region. The mass of the core is $\sim$0.2-0.4 $M_{\odot}$, which is an order of magnitude smaller than typical low-mass prestellar cores. Taking into account a typical core to star formation efficiency for prestellar cores ($\sim$20%-40%) in nearby molecular clouds, brown dwarf(s) or very low-mass star(s) may be going to be formed in this core. We have found possible substructures at the high-density portion of the core, although much higher angular resolution observation is needed to clearly confirm them. The subsequent N$_2$H$^+$ and N$_2$D$^+$ observations using the Nobeyama 45 m telescope have confirmed the high-deuterium fractionation ($\sim$30%). These dynamically and chemically evolved features indicate that this core is on the verge of proto-brown dwarf or very low-mass star formation and is an ideal source to investigate the initial conditions of such low-mass objects via gravitational collapse and/or fragmentation of the filamentary cloud complex.

Filaments are ubiquitous structures in molecular clouds and play an important role in the mass assembly of stars. We present results of dynamical stability analyses for filaments in the infrared dark cloud G14.225$-$0.506, where a delayed onset of massive star formation was reported in the two hubs at the convergence of multiple filaments of parsec length. Full-synthesis imaging is performed with the Atacama Large Millimeter/submillimeter Array (ALMA) to map the $\mathrm{N_2H^+} \; (1-0)$ emission in two hub-filament systems with a spatial resolution of $\sim 0.034 \; \mathrm{pc}$. Kinematics are derived from sophisticated spectral fitting algorithm that accounts for line blending, large optical depth, and multiple velocity components. We identify five velocity coherent filaments and derive their velocity gradients with principal component analysis. The mass accretion rates along the filaments are up to $10^{-4} \; \mathrm{M_\odot \, \mathrm{yr^{-1}}}$ and are significant enough to affect the hub dynamics within one free-fall time ($\sim 10^5 \; \mathrm{yr}$). The $\mathrm{N_2H^+}$ filaments are in equilibrium with virial parameter $\alpha_\mathrm{vir} \sim 1.2$. We compare $\alpha_\mathrm{vir}$ measured in the $\mathrm{N_2H^+}$ filaments, $\mathrm{NH_3}$ filaments, $870 \; \mu\mathrm{m}$ dense clumps, and $3 \; \mathrm{mm}$ dense cores. The decreasing trend in $\alpha_\mathrm{vir}$ with decreasing spatial scales persists, suggesting an increasingly important role of gravity at small scales. Meanwhile, $\alpha_\mathrm{vir}$ also decreases with decreasing non-thermal motions. In combination with the absence of high-mass protostars and massive cores, our results are consistent with the global hierarchical collapse scenario.

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.

How stellar feedback from high-mass stars (e.g., H{\sc ii} regions) influences the surrounding interstellar medium and regulates new star formation is still unclear. To address this question, we observed the G9.62+0.19 complex in 850 $\mu$m continuum with the JCMT/POL-2 polarimeter. An ordered magnetic field has been discovered in its youngest clump, the G9.62 clump. The magnetic field strength is determined to be $\sim$1 mG. Magnetic field plays a larger role than turbulence in supporting the clump. However, the G9.62 clump is still unstable against gravitational collapse even if thermal, turbulent, and magnetic field support are taken into account all together. The magnetic field segments in the outskirts of the G9.62 clump seem to point toward the clump center, resembling a dragged-in morphology, indicating that the clump is likely undergoing magnetically-regulated global collapse. However, The magnetic field in its central region is aligned with the shells of the photodissociation regions (PDRs) and is approximately parallel to the ionization (or shock) front, indicating that the magnetic field therein is likely compressed by the expanding H{\sc ii} regions that formed in the same complex.

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

We present the results of CS J=2-1 mapping observations towards 39 massive star-forming regions selected from the previous CO line survey of cold IRAS sources with high-velocity CO flows along the Galactic plane (Yang et al. 2002). All sources are detected in CS J=2-1 showing the existence of CS clumps around the IRAS sources. However, one-third of the sources are not deeply embedded in the dense clumps by comparison of the central powering IRAS sources and the morphologies of CS clumps. Physical parameters of the dense molecular clumps are presented. We have identified 12 high-mass protostellar object (HMPO) candidates by checking the association between the dense cores and the IRAS sources, the detection of water maser, and the radio properties towards the IRAS sources. We find that the HMPO sources are characterized by low FIR luminosity to virial mass ratios since they are in very early evolutionary stages when the massive protostars have not reached their full luminosities, which are typical for zero-age main sequence stars. Large turbulent motion in the HMPO sources may be largely due to the large kinetic energy ejected by the central protostars formed in the dense clumps. However, alternative means or undetected outflows may also be responsible for the turbulence in the clumps.

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

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.

In order to understand the initial conditions and early evolution of star formation in a wide range of Galactic environments, we carried out an investigation of 64 \textit{Planck} Galactic Cold Clumps (PGCCs) in the second quadrant of the Milky Way. Using the $^{13}$CO and C$^{18}$O $J = 1 - 0$ lines, and 850\,$\mu$m continuum observations, we investigated cloud fragmentation and evolution associated with star formation. We extracted 468 clumps and 117 cores from the $^{13}$CO line and 850\,$\mu$m continuum maps, respectively. We make use of the Bayesian Distance Calculator and derived the distances of all 64 PGCCs. We found that in general, the mass-size plane follows a relation of $m\sim r^{1.67}$. At a given scale, the masses of our objects are around 1/10 of that of typical Galactic massive star-forming regions. Analysis of the clump and core masses, virial parameters, densities, and mass-size relation suggests that the PGCCs in our sample have a low core formation efficiency ($\sim$3.0\%), and most PGCCs are likely low-mass star-forming candidates. Statistical study indicates that the 850\,$\mu$m cores are more turbulent, more optically thick, and denser than the $^{13}$CO clumps for star formation candidates, suggesting that the 850\,$\mu$m cores are likely more appropriate future star-formation candidates than the $^{13}$CO clumps.

Using Galactic Plane surveys, we have selected a massive (1200 M$_\odot$), cold (14 K) 3.6-70 $\mu$m dark IRDC G331.372-00.116. This IRDC has the potential to form high-mass stars and, given the absence of current star formation signatures, it seems to represent the earliest stages of high-mass star formation. We have mapped the whole IRDC with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.1 and 1.3 mm in dust continuum and line emission. The dust continuum reveals 22 cores distributed across the IRDC. In this work, we analyze the physical properties of the most massive core, ALMA1, which has no molecular outflows detected in the CO (2-1), SiO (5-4), and H$_2$CO (3-2) lines. This core is relatively massive ($M$ = 17.6 M$_\odot$), subvirialized (virial parameter $\alpha_{vir}=M_{vir}/M=0.14$), and is barely affected by turbulence (transonic Mach number of 1.2). Using the HCO$^+$ (3-2) line, we find the first detection of infall signatures in a relatively massive, prestellar core (ALMA1) with the potential to form a high-mass star. We estimate an infall speed of 1.54 km s$^{-1}$ and a high accretion rate of 1.96 $\times$ 10$^{-3}$ M$_\odot$ yr$^{-1}$. ALMA1 is rapidly collapsing, out of virial equilibrium, more consistent with competitive accretion scenarios rather than the turbulent core accretion model. On the other hand, ALMA1 has a mass $\sim$6 times larger than the clumps Jeans mass, being in an intermediate mass regime ($M_{J}=2.7<M\lesssim$ 30 M$_\odot$), contrary to what both the competitive accretion and turbulent core accretion theories predict.

We have investigated the TUKH122 prestellar core in the Orion A cloud using ALMA 3 mm dust continuum, N$_2$H$^+$ ($J=1-0$), and CH$_3$OH ($J_K=2_K-1_K$) molecular line observations. Previous studies showed that TUKH122 is likely on the verge of star formation because the turbulence is almost dissipated and chemically evolved among other starless cores in the Orion A cloud. By combining ALMA 12-m and ACA data, we recover extended emission with a resolution of $\sim5"$ corresponding to 0.01 pc and identify 6 condensations with a mass range of $0.1-0.4$ $M_\odot$ and a radius of $\lesssim0.01$ pc. These condensations are gravitationally bound following a virial analysis and are embedded in the filament including the elongated core with a mass of $\sim29$ $M_\odot$ and a radial density profile of $r^{-1.6}$ derived by {\it Herschel}. The separation of these condensations is $\sim0.035$ pc, consistent with the thermal jeans length at a density of $4.4\times10^5$ cm$^{-3}$. This density is similar to the central part of the core. We also find a tendency that the N$_2$H$^+$ molecule seems to deplete at the dust peak condensation. This condensation may be beginning to collapse because the linewidth becomes broader. Therefore, the fragmentation still occurs in the prestellar core by thermal Jeans instability and multiple stars are formed within the TUKH122 prestellar core. The CH$_3$OH emission shows a large shell-like distribution and surrounds these condensations, suggesting that the CH$_3$OH molecule formed on dust grains is released into gas phase by non-thermal desorption such as photoevaporation caused by cosmic-ray induced UV radiation.

Planck Galactic Cold Clumps (PGCCs) possibly represent the early stages of star formation. To understand better the properties of PGCCs, we studied 16 PGCCs in the L1495 cloud with molecular lines and continuum data from Herschel, JCMT/SCUBA-2 and the PMO 13.7 m telescope. Thirty dense cores were identified in 16 PGCCs from 2-D Gaussian fitting. The dense cores have dust temperatures of $T_{\rm d}$ = 11-14 K, and H$_{2}$ column densities of $N_{\rm H_{2}}$ = 0.36-2.5$\times10^{22}$ cm$^{-2}$. We found that not all PGCCs contain prestellar objects. In general, the dense cores in PGCCs are usually at their earliest evolutionary stages. All the dense cores have non-thermal velocity dispersions larger than the thermal velocity dispersions from molecular line data, suggesting that the dense cores may be turbulence-dominated. We have calculated the virial parameter $\alpha$ and found that 14 of the dense cores have $\alpha$ $<$ 2, while 16 of the dense cores have $\alpha$ $>$ 2. This suggests that some of the dense cores are not bound in the absence of external pressure and magnetic fields. The column density profiles of dense cores were fitted. The sizes of the flat regions and core radii decrease with the evolution of dense cores. CO depletion was found to occur in all the dense cores, but is more significant in prestellar core candidates than in protostellar or starless cores. The protostellar cores inside the PGCCs are still at a very early evolutionary stage, sharing similar physical and chemical properties with the prestellar core candidates.

Filamentary structures are ubiquitous in high-mass star-forming molecular clouds. Their relation with high-mass star formation is still to be understood. Here we report interferometric observations toward 8 filamentary high-mass star-forming clouds. A total of 50 dense cores are identified in these clouds, most of which present signatures of high-mass star formation. Five of them are not associated with any star formation indicators, hence are prestellar core candidates. Evolutionary phases of these cores and their linewidths, temperatures, NH$_3$ abundances, and virial parameters are found to be correlated. In a sub-sample of 4 morphologically well-defined filaments, we find that their fragmentation can not be solely explained by thermal or turbulence pressure support. We also investigate distributions of gas temperatures and non-thermal motions along the filaments, and find a spatial correlation between non-thermal linewidths and star formation activities. We find evidence of gas flows along these filaments, and derive an accretion rate along filaments of $\sim$10$^{-4}$ M$_\odot$ yr$^{-1}$. These results suggest a strong relationship between massive filaments and high-mass star formation, through i) filamentary fragmentation in very early evolutionary phases to form dense cores, ii) accretion flows along filaments that are important for the growth of dense cores and protostars, and iii) enhancement of non-thermal motion in the filaments by the feedback or accretion during star formation.

We make the first data release (DR1) of the ongoing ESO Public Survey SAMPLING (http://dx.doi.org/10.7910/DVN/0L8NHX). DR1 comprises of 124 fields distributed in $70^\circ < l < 216^\circ$, $-35^\circ < b < 25^\circ$. The 12CO and 13CO (2-1) cubes are gridded in $8"$ pixels, with an effective resolution of $36"$. The channel width is 0.33 km/s and the RMS noise is $T_{\rm mb}<0.2$ K. Once completed, SAMPLING and complementary surveys will initiate the first major step forward to characterize molecular clouds and star formation on truly Galactic scales.

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.

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 ~2000 PGCCs in J=1-0 transitions of CO isotopologues and ~1000 PGCCs in 850 micron continuum emisison. 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, the length and the mean line-mass (M/L) of the G26 filament are ~6200 Msun, ~12 pc and ~500 Msun/pc, respectively. Ten massive clumps including eight starless ones are found along the filament. The most massive Clump as a whole may be still in global collapse while its denser part seems to be undergoing expansion due to outflow feedback. The fragmentation in G26 filament from cloud scale to clump scale is in agreement with gravitational fragmentation of an isothermal, non-magnetized, 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 due to large-scale compression flows evidenced by the temperature and velocity gradients across its natal cloud.

At the Astrochemistry Focus Group, we discussed what is still missing in our understanding even with new knowledge given at this conference, and what can be done for that within 10 years from now. Still missing in understanding are UV-photons and cosmic-rays interactions with icy dust grains, Sulphur and Phosphorus chemistry, Metallicity effect, Duration (time) effect, COM formation and destruction, phase transition, dust-gas interface, dust evolution, etc. What we should do are multi-scale high spectral resolution molecular observations, laboratory work, theory, radiative transfer, etc. We need careful modeling without simplifying things. Next, I introduce our research on Planck cold clumps. 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\%. 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.

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

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.

Stellar feedback from high-mass stars (e.g., H{\sc ii} regions) can strongly influence the surrounding interstellar medium and regulate star formation. Our new ALMA observations reveal sequential high-mass star formation taking place within one sub-virial filamentary clump (the G9.62 clump) in the G9.62+0.19 complex. The 12 dense cores (MM 1-12) detected by ALMA are at very different evolutionary stages, from starless core phase to UC H{\sc ii} region phase. Three dense cores (MM6, MM7/G, MM8/F) are associated with outflows. The mass-velocity diagrams of outflows associated with MM7/G and MM8/F can be well fitted with broken power laws. The mass-velocity diagram of SiO outflow associated with MM8/F breaks much earlier than other outflow tracers (e.g., CO, SO, CS, HCN), suggesting that SiO traces newly shocked gas, while the other molecular lines (e.g., CO, SO, CS, HCN) mainly trace the ambient gas continuously entrained by outflow jets. Five cores (MM1, MM3, MM5, MM9, MM10) are massive starless core candidates whose masses are estimated to be larger than 25 M$_{\sun}$, assuming a dust temperature of $\leq$ 20 K. The shocks from the expanding H{\sc ii} regions ("B" \& "C") to the west may have great impact on the G9.62 clump through compressing it into a filament and inducing core collapse successively, leading to sequential star formation. Our findings suggest that stellar feedback from H{\sc ii} regions may enhance the star formation efficiency and suppress the low-mass star formation in adjacent pre-existing massive clumps.

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

Understanding the formation of wide binary systems of very low mass stars (M $\le$ 0.1 Msun) is challenging. The most obvious route is via widely separated low-mass collapsing fragments produced through turbulent fragmentation of a molecular core. However, close binaries/multiples from disk fragmentation can also evolve to wide binaries over a few initial crossing times of the stellar cluster through tidal evolution. Finding an isolated low mass wide binary system in the earliest stage of formation, before tidal evolution could occur, would prove that turbulent fragmentation is a viable mechanism for (very) low mass wide binaries. Here we report high resolution ALMA observations of a known wide-separation protostellar binary, showing that each component has a circumstellar disk. The system is too young to have evolved from a close binary and the disk axes are misaligned, providing strong support for the turbulent fragmentation model. Masses of both stars are derived from the Keplerian rotation of the disks; both are very low mass stars.

We present the results from a series of ground-based radio observations toward a Planck Galactic Cold Clump (PGCC), PGCC G108.84-00.81, which is located in one curved filamentary cloud in the vicinity of an extended HII region Sh2-152 and SNR G109.1-1.0. PGCC G108.84-00.81 is mainly composed of two clumps, "G108-N" and "G108-S". In the 850 micron dust continuum emission map, G108-N is shown as one component while G108-S is fragmented into four components. There is no infrared source associated with G108-N while there are two infrared sources (IRS 1 and IRS 2) associated with G108-S. The total mass of G108-N is larger than the jeans mass, suggesting that G108-N is gravitationally unstable and a potential place for a future star formation. The clump properties of G108-N and G108-S such as the gas temperature and the column density, are not distinctly different. However, G108-S is slightly more evolved than G108-N, in the consideration of the CO depletion factor, molecular abundances, and association with infrared sources. G108-S seems to be affected by the compression from Sh2-152, while G108-N is relatively protected from the external effect

The Infrared Dark Cloud (IRDC) G028.23-00.19 hosts a massive (1,500 Msun), cold (12 K), and 3.6-70 um IR dark clump (MM1) that has the potential to form high-mass stars. We observed this prestellar clump candidate with the SMA (~3.5" resolution) and JVLA (~2.1" resolution) in order to characterize the early stages of high-mass star formation and to constrain theoretical models. Dust emission at 1.3 mm wavelength reveals 5 cores with masses <15 Msun. None of the cores currently have the mass reservoir to form a high-mass star in the prestellar phase. If the MM1 clump will ultimately form high-mass stars, its embedded cores must gather a significant amount of additional mass over time. No molecular outflows are detected in the CO (2-1) and SiO (5-4) transitions, suggesting that the SMA cores are starless. By using the NH3 (1,1) line, the velocity dispersion of the gas is determined to be transonic or mildly supersonic (DeltaV_nt}/DeltaV_th ~1.1-1.8). The cores are not highly supersonic as some theories of high-mass star formation predict. The embedded cores are 4 to 7 times more massive than the clump thermal Jeans mass and the most massive core (SMA1) is 9 times less massive than the clump turbulent Jeans mass. These values indicate that neither thermal pressure nor turbulent pressure dominates the fragmentation of MM1. The low virial parameters of the cores (0.1-0.5) suggest that they are not in virial equilibrium, unless strong magnetic fields of ~1-2 mG are present. We discuss high-mass star formation scenarios in a context based on IRDC G028.23-00.19, a study case believed to represent the initial fragmentation of molecular clouds that will form high-mass stars.

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 have performed a dense core survey toward the Infrared Dark Cloud G14.225-0.506 at 3 mm continuum emission with the Atacama Large Millimeter/Submillimeter Array (ALMA). This survey covers the two hub-filament systems with an angular resolution of $\sim3$\arcsec ($\sim0.03$ pc). We identified 48 dense cores. Twenty out of the 48 cores are protostellar due to their association with young stellar objects (YSOs) and/or X-ray point-sources, while the other 28 cores are likely prestellar and unrelated with known IR or X-ray emission. Using APEX 870 $\mu$m continuum emission, we also identified the 18 clumps hosting these cores. Through virial analysis using the ALMA N$_2$H$^+$ and VLA/Effelsberg NH$_3$ molecular line data, we found a decreasing trend in the virial parameter with decreasing scales from filaments to clumps, and then to cores. The virial parameters of $0.1-1.3$ in cores, indicate that cores are likely undergoing dynamical collapse. The cumulative Core Mass Function (CMF) for the prestellar cores candidates has a power law index of $\alpha=1.6$, with masses ranging from 1.5 to 22 $M_\odot$. We find no massive prestellar or protostellar cores. Previous studies suggest that massive O-tpye stars have not been produced yet in this region. Therefore, high-mass stars should be formed in the prestellar cores by accreting a significant amount of gas from the surrounding medium. Another possibility is that low-mass YSOs become massive by accreting from their parent cores that are fed by filaments. These two possibilities might be consistent with the scenario of global hierarchical collapse.

We observed 146 Galactic clumps in HCN (4-3) and CS (7-6) with the Atacama Submillimeter Telescope Experiment (ASTE) 10-m telescope. A tight linear relationship between star formation rate and gas mass traced by dust continuum emission was found for both Galactic clumps and the high redshift (z>1) star forming galaxies (SFGs), indicating a constant gas depletion time of ~100 Myr for molecular gas in both Galactic clumps and high z SFGs. However, low z galaxies do not follow this relation and seem to have a longer global gas depletion time. The correlations between total infrared luminosities (L_TIR) and molecular line luminosities (L'_mol) of HCN (4-3) and CS (7-6) are tight and sublinear extending down to clumps with LTIR~10^{3} L_sun. These correlations become linear when extended to external galaxies. A bimodal behavior in the LTIR--L'mol correlations was found for clumps with different dust temperature, luminosity-to-mass ratio, and sigmaline/sigmavir. Such bimodal behavior may be due to evolutionary effects. The slopes of LTIR--L'mol correlations become more shallow as clumps evolve. We compared our results with lower J transition lines in wu et al. (2010). The correlations between clump masses and line luminosities are close to linear for low effective excitation density tracers but become sublinear for high effective excitation density tracers for clumps with LTIR larger than LTIR~10^4.5 Lsun. High effective excitation density tracers cannot linearly trace the total clump masses, leading to a sublinear correlations for both Mclump-L'mol and LTIR-L'mol relations.

Most bipolar outflows are associated with individual young stellar objects and have small opening angles. Here we report the discovery of an extremely wide-angle ($\sim$180$\arcdeg$) bipolar outflow ("EWBO") in a cluster forming region AFGL 5142 from low-velocity emission of the HCN (3-2) and HCO$^{+}$ (3-2) lines. This bipolar outflow is along a north-west to south-east direction with a line-of-sight flow velocity of about 3 km~s$^{-1}$ and is spatially connected to the high-velocity jet-like outflows. It seems to be a collection of low-velocity material entrained by the high-velocity outflows due to momentum feedback. The total ejected mass and mass loss rate due to both high velocity jet-like outflows and the "EWBO" are $\sim$24.5 M$_{\sun}$ and $\sim1.7\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, respectively. Global collapse of the clump is revealed by the "blue profile" in the HCO$^{+}$ (1-0) line. A hierarchical network of filaments was identified in NH$_{3}$ (1,1) emission. Clear velocity gradients of the order of 10 km~s$^{-1}$~pc$^{-1}$ are found along filaments, indicating gas inflow along the filaments. The sum of the accretion rate along filaments and mass infall rate along the line of sight is $\sim$3.1$\times10^{-3}$ M$_{\sun}$~yr$^{-1}$, which exceeds the total mass loss rate, indicating that the central cluster is probably still gaining mass. The central cluster is highly fragmented and 22 condensations are identified in 1.1 mm continuum emission. The fragmentation process seems to be determined by thermal pressure and turbulence. The magnetic field may not play an important role in fragmentation.

We report the detection of a wing component in NH$_3$ emission toward the starless core TUKH122 with subthermal turbulence in the Orion A cloud. This NH$_3$ core is suggested to be on the verge of star formation because the turbulence inside the NH$_3$ core is almost completely dissipated, and also because it is surrounded by CCS, which resembles the prestellar core L1544 in Taurus showing infall motions. Observations were carried out with the Nobeyama 45 m telescope at 0.05 km s$^{-1}$ velocity resolution. We find that the NH$_3$ line profile consists of two components. The quiescent main component has a small linewidth of 0.3 km s$^{-1}$ dominated by thermal motions, and the red-shifted wing component has a large linewidth of 1.36 km s$^{-1}$ representing turbulent motions. These components show kinetic temperatures of 11 K and $<$ 30 K, respectively. Furthermore, there is a clear velocity offset between the NH$_3$ quiescent gas ($VLSR=3.7$ km s$^{-1}$) and the turbulent gas ($VLSR=4.4$ km s$^{-1}$). The centroid velocity of the turbulent gas corresponds to that of the surrounding gas traced by the $^{13}$CO ($J=1-0$) and CS ($J=2-1$) lines. LVG model calculations for CS and CO show that the turbulent gas has a temperature of $8-13$ K and an H$_2$ density of ~ $10^4$ cm$^{-3}$, suggesting that the temperature of the turbulent component is also ~ 10 K. The detections of both NH$_3$ quiescent and wing components may indicate a sharp transition from the turbulent parent cloud to the quiescent dense core.

The scientific impact of a facility is the most important measure of its success. Monitoring and analysing the scientific return can help to modify and optimise operations and adapt to the changing needs of scientific research. The methodology that we have developed to monitor the scientific productivity of the ALMA Observatory, as well as the first results, are described. We focus on the outcome of the first cycle (Cycle 0) of ALMA Early Science operations. Despite the fact that only two years have passed since the completion of Cycle 0 and operations have already changed substantially, this analysisconfirms the effectiveness of the underlying concepts. We find that ALMA is fulfilling its promise as a transformational facility for the observation of the Universe in the submillimetre.

We have analyzed the angular momentum of the molecular cloud cores in the Orion A giant molecular cloud observed in the N2H+ J = 1-0 line with the Nobeyama 45 m radio telescope. We have measured the velocity gradient using position velocity diagrams passing through core centers, and made sinusoidal fitting against the position angle. 27 out of 34 N2H+ cores allowed us to measure the velocity gradient without serious confusion. The derived velocity gradient ranges from 0.5 to 7.8 km/s/pc. We marginally found that the specific angular momentum J/M (against the core radius R) of the Orion N2H+ cores tends to be systematically larger than that of molecular cloud cores in cold dark clouds obtained by Goodman et al., in the J/M-R relation. The ratio beta of rotational to gravitational energy is derived to be beta = 10^{-2.3+/-0.7}, and is similar to that obtained for cold dark cloud cores in a consistent definition. The large-scale rotation of the integral-shaped filament of the Orion A giant molecular cloud does not likely govern the core rotation at smaller scales.

We are performing a series of observations with ground-based telescopes toward Planck Galactic cold clumps (PGCCs) in the $\lambda$ Orionis complex in order to systematically investigate the effects of stellar feedback. In the particular case of PGCC G192.32-11.88, we discovered an extremely young Class 0 protostellar object (G192N) and a proto-brown dwarf candidate (G192S). G192N and G192S are located in a gravitationally bound bright-rimmed clump. The velocity and temperature gradients seen in line emission of CO isotopologues indicate that PGCC G192.32-11.88 is externally heated and compressed. G192N probably has the lowest bolometric luminosity ($\sim0.8$ L$_{\sun}$) and accretion rate (6.3$\times10^{-7}$ M$_{\sun}$~yr$^{-1}$) when compared with other young Class 0 sources (e.g. PACS Bright Red sources (PBRs)) in the Orion complex. It has slightly larger internal luminosity ($0.21\pm0.01$ L$_{\sun}$) and outflow velocity ($\sim$14 km~s$^{-1}$) than the predictions of first hydrostatic cores (FHSCs). G192N might be among the youngest Class 0 sources, which are slightly more evolved than a FHSC. Considering its low internal luminosity ($0.08\pm0.01$ L$_{\odot}$) and accretion rate (2.8$\times10^{-8}$ M$_{\sun}$~yr$^{-1}$), G192S is an ideal proto-brown dwarf candidate. The star formation efficiency ($\sim$0.3\%-0.4\%) and core formation efficiency ($\sim$1\%) in PGCC G192.32-11.88 are significantly smaller than in other giant molecular clouds or filaments, indicating that the star formation therein is greatly suppressed due to stellar feedback.

We present observations of HDCO and H2CO emission toward a sample of 15 Class 0 protostars in the Orion A and B clouds. Of these, eleven protostars are Herschel-identified PACS Bright Red Sources (PBRS) and four are previously identified protostars. Our observations revealed the chemical properties of the PBRS envelope for the first time. The column densities of HDCO and H2CO are derived from single dish observations at an angular resolution of ~20 arcsec (~8400 AU). The degree of deuteration in H2CO ([HDCO]/[H2CO]) was estimated to range from 0.03 to 0.31. The deuterium fractionation of most PBRS (70%) is similar to that of the non-PBRS sources. Three PBRS (30%) exhibit high deuterium fractionation, larger than 0.15. The large variation of the deuterium fractionation of H2CO in the whole PBRS sample may reflect the diversity in the initial conditions of star forming cores. There is no clear correlation between the [HDCO]/[H2CO] ratio and the evolutionary sequence of protostars.

We have observed the HC3N (J=10-9) and N2H+ (J=1-0) lines toward the Vela C molecular clouds with the Mopra 22 m telescope to study chemical characteristics of dense cores. The intensity distributions of these molecules are similar to each other at an angular resolution of 53", corresponding to 0.19 pc suggesting that these molecules trace the same dense cores. We identified 25 local peaks in the velocity-integrated intensity maps of the HC3N and/or N2H+ emission. Assuming LTE conditions, we calculated the column densities of these molecules and found a tendency that N2H+/HC3N abundance ratio seems to be low in starless regions while it seems to be high in star-forming regions, similar to the tendencies in the NH3/CCS, NH3/HC3N, and N2H+/CCS abundance ratios found in previous studies of dark clouds and the Orion A GMC. We suggest that carbon chain molecules, including HC3N, may trace chemically young molecular gas and N-bearing molecules, such as N2H+, may trace later stages of chemical evolution in the Vela C molecular clouds. It may be possible that the N2H+/HC3N abundance ratio of ~ 1.4 divides the star-forming and starless peaks in the Vela C, although it is not as clear as those in NH3/CCS, NH3/HC3N, and N2H+/CCS for the Orion A GMC. This less clear separation may be caused by our lower spatial resolution or the misclassification of star-forming and starless peaks due to the larger distance of the Vela C. It might be also possible that the HC3N (J=10-9) transition is not a good chemical evolution tracer compared with CCS (J=4-3 and 7-6) transitions.

W49N is a mini-starburst in the Milky Way and thus an ideal laboratory for high-mass star formation studies. Due to its large distance (11.1$_{-0.7}^{+0.9}$ kpc), the kinematics inside and between the dense molecular clumps in W49N are far from well understood. The SMA observations resolved the continuum emission into two clumps. The molecular line observation of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) suggests that the two clumps have a velocity difference of $\sim$7 km~s$^{-1}$. The eastern clump is very close to two radio sources "G1" and "G2", and the western clump coincides with a radio source "B". The HCN (3-2) line reveals an extremely energetic outflow, which is among the most energetic molecular outflows in the Milky Way. This is the first report of high-velocity molecular outflow detection in W49N. The outflow jet might be in precession, which could account for the distribution, velocity and rotation of water maser spots. Three absorption systems are identified in HCO$^{+}$ (3-2) spectra. The absorption features are blueshifted with respect to the emission of SO$_{2}$ (28$_{4,24}$-28$_{3,25}$) lines, indicating that a cold layer is expanding in front of the warm gas. Further analysis indicates that the expansion is decelerated from the geometric expansion centers.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

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

We present Atacama Large Millimeter/submillimeter Array (ALMA) 1.3 mm continuum images of the asteroid 3 Juno obtained with an angular resolution of 0.042 arcseconds (60 km at 1.97 AU). The data were obtained over a single 4.4 hr interval, which covers 60% of the 7.2 hr rotation period, approximately centered on local transit. A sequence of ten consecutive images reveals continuous changes in the asteroid's profile and apparent shape, in good agreement with the sky projection of the three-dimensional model of the Database of Asteroid Models from Inversion Techniques. We measure a geometric mean diameter of 259pm4 km, in good agreement with past estimates from a variety of techniques and wavelengths. Due to the viewing angle and inclination of the rotational pole, the southern hemisphere dominates all of the images. The median peak brightness temperature is 215pm13 K, while the median over the whole surface is 197pm15 K. With the unprecedented resolution of ALMA, we find that the brightness temperature varies across the surface with higher values correlated to the subsolar point and afternoon areas, and lower values beyond the evening terminator. The dominance of the subsolar point is accentuated in the final four images, suggesting a reduction in the thermal inertia of the regolith at the corresponding longitudes, which are possibly correlated to the location of the putative large impact crater. These results demonstrate ALMA's potential to resolve thermal emission from the surface of main belt asteroids, and to measure accurately their position, geometric shape, rotational period, and soil characteristics.

We present initial results of very high resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of the $z$=3.042 gravitationally lensed galaxy HATLAS J090311.6+003906 (SDP.81). These observations were carried out using a very extended configuration as part of Science Verification for the 2014 ALMA Long Baseline Campaign, with baselines of up to 15 km. We present continuum imaging at 151, 236 and 290 GHz, at unprecedented angular resolutions as fine as 23 milliarcseconds (mas), corresponding to an un-magnified spatial scale of ~180 pc at z=3.042. The ALMA images clearly show two main gravitational arc components of an Einstein ring, with emission tracing a radius of ~1.5". We also present imaging of CO(10-9), CO(8-7), CO(5-4) and H2O line emission. The CO emission, at an angular resolution of ~170 mas, is found to broadly trace the gravitational arc structures but with differing morphologies between the CO transitions and compared to the dust continuum. Our detection of H2O line emission, using only the shortest baselines, provides the most resolved detection to date of thermal H2O emission in an extragalactic source. The ALMA continuum and spectral line fluxes are consistent with previous Plateau de Bure Interferometer and Submillimeter Array observations despite the impressive increase in angular resolution. Finally, we detect weak unresolved continuum emission from a position that is spatially coincident with the center of the lens, with a spectral index that is consistent with emission from the core of the foreground lensing galaxy.

We investigate the effects of Supergiant Shells (SGSs) and their interaction on dense molecular clumps by observing the Large Magellanic Cloud (LMC) star forming regions N48 and N49, which are located between two SGSs, LMC 4 and LMC 5. $^{12}$CO ($J$=3-2, 1-0) and $^{13}$CO ($J$=1-0) observations with the ASTE and Mopra telescopes have been carried out towards these regions. A clumpy distribution of dense molecular clumps is revealed with 7 pc spatial resolution. Large velocity gradient analysis shows that the molecular hydrogen densities ($n({\rm H}_2)$) of the clumps are distributed from low to high density ($10^3$-$10^5$ cm$^{-3}$) and their kinetic temperatures ($T_{\rm kin}$) are typically high (greater than $50$ K). These clumps seem to be in the early stages of star formation, as also indicated from the distribution of H$\alpha$, young stellar object candidates, and IR emission. We found that the N48 region is located in the high column density HI envelope at the interface of the two SGSs and the star formation is relatively evolved, whereas the N49 region is associated with LMC 5 alone and the star formation is quiet. The clumps in the N48 region typically show high $n({\rm H}_2)$ and $T_{\rm kin}$, which are as dense and warm as the clumps in LMC massive cluster-forming areas (30 Dor, N159). These results suggest that the large-scale structure of the SGSs, especially the interaction of two SGSs, works efficiently on the formation of dense molecular clumps and stars.

The physical and chemical properties of prestellar cores, especially massive ones, are still far from being well understood due to the lack of a large sample. The low dust temperature ($<$14 K) of Planck cold clumps makes them promising candidates for prestellar objects or for sources at the very initial stages of protostellar collapse. We have been conducting a series of observations toward Planck cold clumps (PCCs) with ground-based radio telescopes. In general, when compared with other star forming samples (e.g. infrared dark clouds), PCCs are more quiescent, suggesting that most of them may be in the earliest phase of star formation. However, some PCCs are associated with protostars and molecular outflows, indicating that not all PCCs are in a prestellar phase. We have identified hundreds of starless dense clumps from the mapping survey with the Purple Mountain Observatory (PMO) 13.7-m telescope. Follow-up observations suggest that these dense clumps are ideal targets to search for prestellar objects.

Telescope bibliographies are integral parts of observing facilities. They are used to associate the published literature with archived observational data, to measure an observatory's scientific output through publication and citation statistics, and to define guidelines for future observing strategies. The ESO and NRAO librarians as well as NAOJ jointly maintain the ALMA (Atacama Large Millimeter/submillimeter Array) bibliography, a database of refereed papers that use ALMA data. In this paper, we illustrate how relevant articles are identified, which procedures are used to tag entries in the database and link them to the correct observations, and how results are communicated to ALMA stakeholders and the wider community. Efforts made to streamline the process will be explained and evaluated, and a first analysis of ALMA papers published after two years of observations will be given.

We imaged two starless molecular cloud cores, TUKH083 and TUKH122, in the Orion A giant molecular cloud in the CCS and ammonia (NH$_3$) emission with the Very Large Array. TUKH122 contains one NH$_3$ core "TUKH122-n," which is elongated and has a smooth oval boundary. Where observed, the CCS emission surrounds the NH$_3$ core. This configuration resembles that of the N$_2$H$^+$ and CCS distribution in the Taurus starless core L1544, a well-studied example of a dense prestellar core exhibiting infall motions. The linewidth of TUKH122-n is narrow (0.20 km s$^{-1}$) in the NH$_3$ emission line and therefore dominated by thermal motions. The smooth oval shape of the core boundary and narrow linewidth in NH$_3$ seem to imply that TUKH122-n is dynamically relaxed and quiescent. TUKH122-n is similar to L1544 in the kinetic temperature (10 K), linear size (0.03 pc), and virial mass ($\sim$ 2 $M_{\odot}$). Our results strongly suggest that TUKH122-n is on the verge of star formation. TUKH122-n is embedded in the 0.2 pc massive (virial mass $\sim$ 30 $M_{\odot}$) turbulent parent core, while the L1544 NH$_3$ core is embedded in the 0.2 pc less-massive (virial mass $\sim$ 10 $M_{\odot}$) thermal parent core. TUKH083 shows complicated distribution in NH$_3$, but was not detected in CCS. The CCS emission toward TUKH083 appears to be extended, and is resolved out in our interferometric observations.

We report on spatially resolved X-ray spectroscopy of the north-eastern part of the mixed morphology supernova remnant (SNR) W28 with {\it XMM-Newton}. The observed field of view includes a prominent and twisted shell emission forming the edge of this SNR as well as part of the center-filled X-ray emission brightening toward the south-west edge of the field of view. The shell region spectra are in general represented by an optically thin thermal plasma emission in collisional ionization equilibrium with a temperature of $\sim$0.3~keV and a density of $\sim$10~cm$^{-3}$, which is much higher than the density obtained for inner parts. In contrast, we detected no significant X-ray flux from one of the TeV $\gamma$-ray peaks with an upper-limit flux of 2.1$\times$10$^{-14}$ erg cm$^{-2}$ s$^{-1}$ in the 2--10~keV band. The large flux ratio of TeV to X-ray, larger than 16, and the spatial coincidence of the molecular cloud and the TeV $\gamma$-ray emission site indicate that the TeV $\gamma$-ray of W28 is $\pi^{0}$-decay emission originating from collisions between accelerated protons and molecular cloud protons. Comparing the spectrum in the TeV band and the X-ray upper limit, we obtained a weak upper limit on the magnetic field strength $B\lesssim$ 1500~$\mu$G.

We have mapped six molecular cloud cores in the Orion A giant molecular cloud (GMC), whose kinetic temperatures range from 10 to 30 K, in CCS and N2H+ with Nobeyama 45 m radio telescope to study their chemical characteristics. We identified 31 intensity peaks in the CCS and N2H+ emission in these molecular cloud cores. It is found for cores with temperatures lower than ~ 25 K that the column density ratio of N(N2H+)/N(CCS) is low toward starless core regions while it is high toward star-forming core regions, in case that we detected both of the CCS and N2H+ emission. This is very similar to the tendency found in dark clouds (kinetic temperature ~ 10 K). The criterion found in the Orion A GMC is N(N2H+)/N(CCS) ~ 2-3. In some cases, the CCS emission is detected toward protostars as well as the N2H+ emission. Secondary late-stage CCS peak in the chemical evolution caused by CO depletion may be a possible explanation for this. We found that the chemical variation of CCS and N2H+ can also be used as a tracer of evolution in warm (10-25 K) GMC cores. On the other hand, some protostars do not accompany N2H+ intensity peaks but are associated with dust continuum emitting regions, suggesting that the N2H+ abundance might be decreased due to CO evaporation in warmer star-forming sites.

The properties of tidally induced arms provide a means to study molecular cloud formation and the subsequent star formation under environmental conditions which in principle are different from quasi stationary spiral arms. We report the properties of a newly discovered molecular gas arm of likely tidal origin at the south of NGC 4039 and the overlap region in the Antennae galaxies, with a resolution of 1"68 x 0"85, using the Atacama Large Millimeter/submillimeter Array science verification CO(2-1) data. The arm extends 3.4 kpc (34") and is characterized by widths of ~ 200 pc (2") and velocity widths of typically \DeltaV ~ 10-20 km/s . About 10 clumps are strung out along this structure, most of them unresolved, with average surface densities of \Sigma_gas ~ 10-100 Msun pc^{-2}, and masses of (1-8) x 10^6 Msun. These structures resemble the morphology of beads on a string, with an almost equidistant separation between the beads of about 350 pc, which may represent a characteristic separation scale for giant molecular associations. We find that the star formation efficiency at a resolution of 6" (600 pc) is in general a factor of 10 higher than in disk galaxies and other tidal arms and bridges. This arm is linked, based on the distribution and kinematics, to the base of the western spiral arm of NGC 4039, but its morphology is different to that predicted by high-resolution simulations of the Antennae galaxies.

We analyze the HCO+ 3-2 and H13CO+ 3-2 line profiles of 27 high-mass star-forming regions to identify asymmetries that are suggestive of mass inflow. Three quantitative measures of line asymmetry are used to indicate whether a line profile is blue, red or neither - the ratio of the temperature of the blue and red peaks, the line skew and the dimensionless parameter delta_v. We find nine HCO+ 3-2 line profiles with a significant blue asymmetry and four with significant red asymmetric profiles. Comparing our HCO+ 3-2 results to HCN 3-2 observations from Wu et al. (2003, 2010), we find that eight of the blue and three of red have profiles with the same asymmetry in HCN. The eight sources with blue asymmetries in both tracers are considered strong candidates for inflow. Quantitative measures of the asymmetry (e.g. delta_v) tend to be larger for HCN. This, combined with possible HCO+ abundance enhancements in outflows, suggests that HCN may be a better tracer of inflow. Understanding the behavior of common molecular tracers like HCO+ in clumps of different masses is important for properly analyzing the line profiles seen in a sample of sources representing a broad range of clump masses. Such studies will soon be possible with the large number of sources with possible self-absorption seen in spectroscopic follow-up observations of clumps identified in the Bolocam Galactic Plane Survey.

The NGC 1333 IRAS 4A protobinary was observed in the 1.3 cm and 6.9 mm continuum and the ammonia and SiO lines, with an angular resolution of about 0.4 arcseconds. The continuum maps show the circumstellar structures of the two protostars, A1 and A2. The A1 system is brighter and more massive than the A2 system. The ratio of mass, including dense gas and protostar, is about 6. The properties of the circumstellar disks and outflows suggest that A1 may be younger than A2. The deflected part of the northeastern jet of A2 is bright in the SiO line, and the distance between the brightest peak and deflection point suggests that the enhancement of SiO takes about 100 yr after the collision with a dense core. The ammonia maps show a small structure that seems to be a part of the obstructing core. The outflow properties were studied by comparing interferometric maps of SiO, ammonia, formaldehyde, and HCN lines. Their overall structures agree well, suggesting that these species are excited by the same mechanism. However, the intensity distributions show that SiO is chemically unique. SiO may be directly linked to the primary jet while the other species may be tracing the entrained ambient gas.

We present observations of HCO+ and H^13CO+, N2H+, HCS+, HCN and HN^13C, SO and ^34SO, CCH, SO_2, and CH_3OH-E towards a sample of 27 high-mass clumps coincident with water maser emission. All transitions are observed with or convolved to nearly identical resolution (30"), allowing for inter-comparison of the clump properties derived from the mapped transitions. We find N2H+ emission is spatially differentiated compared to the dust and the other molecules towards a few very luminous cores (10 of 27) and the N2H+ integrated intensity does not correlate well with dust continuum flux. We calculate the effective excitation density, n_eff, the density required to excite a 1 K line in T_kin=20 K gas for each molecular tracer. The intensity of molecular tracers with larger effective excitation densities (n_eff > 10^5 cm^-3) appear to correlate more strongly with the submillimeter dust continuum intensity. The median sizes of the clumps are anti-correlated with the n_eff of the tracers (which span more than three orders of magnitude). Virial mass is not correlated with n_eff, especially where the lines are optically thick as the linewidths may be broadened significantly by non-virial motions. The median mass surface density and median volume density of the clumps is correlated with n_eff indicating the importance of understanding the excitation conditions of the molecular tracer when deriving the average properties of an ensemble of cores.

We report new mid-infrared observations of the remarkable object IRAS 15099-5856 using the space telescopes AKARI and Spitzer, which demonstrate the presence of prominent crystalline silicate emission in this bright source. IRAS 15099-5856 has a complex morphology with a bright central compact source (IRS1) surrounded by knots, spurs, and several extended (~4') arc-like filaments. The source is seen only at >= 10 um. The Spitzer MIR spectrum of IRS1 shows prominent emission features from Mg-rich crystalline silicates, strong [Ne II] 12.81 um and several other faint ionic lines. We model the MIR spectrum as thermal emission from dust and compare with the Herbig Be star HD 100546 and the luminous blue variable R71, which show very similar MIR spectra. Molecular line observations reveal two molecular clouds around the source, but no associated dense molecular cores. We suggest that IRS1 is heated by UV radiation from the adjacent O star Muzzio 10 and that its crystalline silicates most likely originated in a mass outflow from the progenitor of the supernova remnant (SNR) MSH 15-52. IRS1, which is embedded in the SNR, could have been shielded from the SN blast wave if the progenitor was in a close binary system with Muzzio 10. If MSH15-52 is a remnant of Type Ib/c supernova (SN Ib/c), as has been previously proposed, this would confirm the binary model for SN Ib/c. IRS1 and the associated structures may be the relics of massive star death, as shaped by the supernova explosion, the pulsar wind and the intense ionizing radiation of the embedded O star.

The bipolar jet of the NGC 1333 IRAS 4A2 protostar shows a velocity gradient in the direction perpendicular to the jet axis. This lateral velocity gradient can be seen throughout the jet imaged in a silicon monoxide line, 2500-8700 AU from the driving source, and is consistent with the rotation of the accretion disk. If this gradient is caused by the rotation of the jet around its axis, the average specific angular momentum is about 1.5 x 10^21 cm^2 s^-1. Comparison of the kinematics between the jet and the disk suggests that the jet-launching region on the disk has a radius of about 2 AU, which supports the disk-wind models. The angular momentum transported away by the jet seems to be large enough for the protostar to accrete matter from the disk, confirming the crucial role of jets in the early phase of star formation process.

In order to precisely determine temperature and density of molecular gas in the Large Magellanic Cloud, we made observations of optically thin $^{13}$CO($J=3-2$) transition by using the ASTE 10m telescope toward 9 peaks where $^{12}$CO($J=3-2$) clumps were previously detected with the same telescope. The molecular clumps include those in giant molecular cloud (GMC) Types I (with no signs of massive star formation), II (with HII regions only), and III (with HII regions and young star clusters). We detected $^{13}$CO($J=3-2$) emission toward all the peaks and found that their intensities are 3 -- 12 times lower than those of $^{12}$CO($J=3-2$). We determined the intensity ratios of $^{12}$CO($J=3-2$) to $^{13}$CO($J=3-2$), $R^{12/13}_{3-2}$, and $^{13}$CO($J=3-2$) to $^{13}$CO($J=1-0$), $R^{13}_{3-2/1-0}$, at 45$\arcsec$ resolution. These ratios were used for radiative transfer calculations in order to estimate temperature and density of the clumps. The parameters of these clumps range kinetic temperature $T\mathrm{_{kin}}$ = 15 -- 200 K, and molecular hydrogen gas density $n(\mathrm{H_2})$ = 8$\times 10^2$ -- 7$\times 10^3$ cm$^{-3}$. We confirmed that the higher density clumps show higher kinetic temperature and that the lower density clumps lower kinetic temperature at a better accuracy than in the previous work. The kinetic temperature and density increase generally from a Type I GMC to a Type III GMC. We interpret that this difference reflects an evolutionary trend of star formation in molecular clumps. The $R^{13}_{3-2/1-0}$ and kinetic temperature of the clumps are well correlated with H$\alpha$ flux, suggesting that the heating of molecular gas $n(\mathrm{H_2})$ = $10^3$ -- $10^4$ cm$^{-3}$ can be explained by stellar FUV photons.

We have observed molecular cloud cores in the Orion A giant molecular cloud (GMC) in CCS, HC3N, DNC, and HN13C to study their chemical characteristics. We have detected CCS in the Orion A GMC for the first time. CCS was detected in about a third of the observed cores. The cores detected in CCS are not localized but are widely distributed over the Orion A GMC. The CCS peak intensity of the core tends to be high in the southern region of the Orion A GMC. The HC3N peak intensity of the core also tends to be high in the southern region, while there are HC3N intense cores near Orion KL, which is not seen in CCS. The core associated with Orion KL shows broad HC3N line profile, and star formation activity near Orion KL seems to enhance the HC3N emission. The column density ratio of NH3 to CCS is lower near the middle of the filament, and is higher toward the northern and southern regions along the Orion A GMC filament. This ratio is known to trace the chemical evolution in nearby dark cloud cores, but seems to be affected by core gas temperature in the Orion A GMC: cores with low NH3 to CCS column density ratios tend to have warmer gas temperature. The value of the column density ratio of DNC to HN13C is generally similar to that in dark cloud cores, but becomes lower around Orion KL due to higher gas temperature.

The NGC 1333 IRAS 4A protobinary was observed in the ammonia (2, 2) and (3, 3) lines with an angular resolution of 0.3 arcsec. The ammonia emission source of IRAS 4A2 is elongated in the direction perpendicular to the bipolar jet and has a size of 0.55 arcsec or 130 AU. This emission structure was interpreted as a circumstellar disk around the IRAS 4A2 protostar, and the rotation kinematics of the disk was investigated by making a position-velocity diagram along the major axis. Assuming a Keplerian rotation, the disk has a rotation velocity of 1.8 km/s at a radius of 20 AU, which implies a central object of about 0.08 solar masses. The collapse age of the protostar is about 50,000 yr. The mass, accretion rate, and age are consistent with what are expected from the standard theory of low-mass star formation. If IRAS 4A2 grows at this rate, it may become a star similar to the Sun.

New 12CO J=4-3 and 13CO J=3-2 observations of the N159 region in the Large Magellanic Cloud have been made. The 12CO J=4-3 distribution is separated into three clumps. These new measurements toward the three clumps are used in coupled calculations of molecular rotational excitation and line radiation transfer, along with other transitions of the 12CO as well as the isotope transitions of 13CO. The temperatures and densities are determined to be ~70-80K and ~3x10^3 cm-3 in N159W and N159E and ~30K and ~1.6x10^3 cm-3 in N159S. These results are compared with the star formation activity. The N159E clump is associated with embedded cluster(s) as observed at 24 micron and the derived high temperature is explained as due to the heating by these sources. The N159E clump is likely responsible for a dark lane in a large HII region by the dust extinction. The N159W clump is associated with embedded clusters mainly toward the eastern edge of the clump only. These clusters show offsets of 20"-40" from the 12CO J=4-3 peak and are probably responsible for heating indicated by the derived high temperature. The N159W clump exhibits no sign of star formation toward the 12CO J=4-3 peak position and its western region. We suggest that the N159W peak represents a pre-star-cluster core of ~105M_sol which deserves further detailed studies. Note that recent star formation took place between N159W and N159E as indicated by several star clusters and HII regions, while the natal molecular gas toward the stars have already been dissipated by the ionization and stellar winds of the OB stars. The N159S clump shows little sign of star formation as is consistent with the lower temperature and somewhat lower density. The N159S clump is also a candidate for future star formation.

Circularly polarized 3.5 cm continuum emission was detected toward three radio sources in the R CrA region using the Very Large Array. The Class I protostar IRS 5b persistently showed polarized radio emission with a constant helicity over 8 yr, which suggests that its magnetosphere has a stable configuration. There is a good correlation between the Stokes I and Stokes V fluxes, and the fractional polarization is about 0.17. During active phases the fractional polarization is a weakly decreasing function of Stokes I flux, which suggests that IRS 5b is phenomenologically similar to other types of flare stars such as RS CVn binaries. The variability timescale of the polarized flux is about a month, and the magnetosphere of IRS 5b must be very large in size. The Class I protostar IRS 7A was detected once in circularly polarized radio emission, even though IRS 7A drives a thermal radio jet. This detection implies that the radio emission from the magnetosphere of a young protostar can escape the absorption by the partially ionized wind at least once in a while. The properties of IRS 7A and IRS 5b suggests that Class I protostars have organized peristellar magnetic fields of a few kilogauss and that the detectability of magnetospheric emission may depend on the evolutionary status of protostar. Also reported is the detection of circularly polarized radio emission toward the variable radio source B5.

The R CrA region was observed in the 3.5 and 6.2 cm continuum with high angular resolutions (0.6--1.7 arcseconds) using the Very Large Array. Archival data sets were also analyzed for comparison, which provided angular resolutions up to 0.3 arcseconds. A cluster of young stellar objects was detected, and a rich array of star forming activities was revealed. IRS 7A showed an enhanced outflow activity recently. The main peak of IRS 7A positionally coincides with an X-ray source, which suggests that the X-ray emission is directly related to the central protostar. The Class 0 source SMA 2 is associated with a double radio source, B 9a and 9b, and seems to be driving two outflows. The B 9 complex is probably a multiple-protostar system. Both B 9a and 9b are nonthermal radio sources with negative spectral indices. IRS 7B is a compact radio source surrounded by an extended structure. The compact source corresponds to the Class 0/I source SMA 1, and it is also closely associated with an X-ray source, suggesting that magnetic activities start early in the protostellar stage of evolution. The extended structure of IRS 7B may be a bipolar outflow. IRS 5 was resolved into two sources with a separation of 0.9 arcseconds. Both IRS 5a and 5b display radio flares and X-ray emission, suggesting that energetic magnetic processes are active in both members. The month-scale active phase of IRS 5b implies that the flare activity must involve large-scale magnetic fields. During the strong flare event of IRS 5b in 1998, IRS 5a also showed an enhanced level of radio emission. This concurrent activity suggests that IRS 5 may be an interacting young binary system, but the interaction mechanism is unknown. Alternatively, what was seen in the radio images could be a circumbinary halo.

Using a sub-pixel event repositioning technique, we spatially resolved X-ray emission from the infrared double system IRS 5 in the R Corona Australis molecular cloud with ~0.8" separation. As far as we know, this result - obtained from 8 Chandra archival observations between 2000 and 2005 - is the first X-ray study of individual sources in a Class I protostar binary system with a projected separation of less than 200 AU. We extracted light curves and spectra of the individual sources using a two-dimensional image fitting method. IRS 5a at the south, the source which was brighter in the near-infrared, showed three X-ray flares lasting >20 ksec, reminiscent of X-ray flares from pre-main sequence stars, while the northern source (IRS 5b) was quiescent in X-rays in all the observations except for a 2005 August 9 observation with a factor of ~2 flux enhancement. In quiescence, these sources showed almost identical X-ray spectra, with NH ~4e22 cm-2, kT ~2 keV, and log Lx ~30.2-3 ergs s-1. IRS 5a showed plasma at temperatures up to kT ~5-6 keV during flares, while the column density of IRS 5b increased by a factor of 2 during an observation on 2005 August 9. We discuss the evolutionary stages and variation of the X-ray activity of these sources.

The ``integral-shaped filament'' of the Orion A giant molecular cloud was mapped in N2H+, and its northern end, the OMC-2/3 region was observed also in HC3N and CCS. The N2H+ distribution is similar to the dust continuum distribution, except for the central part of the Orion Nebula. The distribution of H13CO+ holds resemblance to that of dust continuum, but the N2H+ distribution looks more similar to dust continuum distribution. The N-bearing molecules, N2H+ and NH3 seem to be more intense in OMC-2, compared with the H13CO+ and CS distribution. We identified 34 cloud cores from N2H+ data. Over the Orion Nebula region, the N2H+ linewidth is large (1.1-2.1 km/s). In the OMC-2/3 region, it becomes moderate (0.5-1.3 km/s), and it is smaller (0.3-1.1 km/s) in the south of the Orion Nebula. On the other hand, the gas kinetic temperature of the quiescent cores observed in N2H+ is rather constant (~ 20 K) over the $\int$-shaped filament. We detected no CCS emission in the OMC-2/3 region. In general, N2H+ and HC3N distribution is quite similar in the OMC-2/3 region, but we observed displacement between N2H+ and HC3N over 2' scale in OMC-3, which has a chain of Class 0-I protostars (candidates).

We have carried out sub-mm 12CO(J=3-2) observations of 6 giant molecular clouds (GMCs) in the Large Magellanic Cloud (LMC) with the ASTE 10m sub-mm telescope at a spatial resolution of 5 pc and very high sensitivity. We have identified 32 molecular clumps in the GMCs and revealed significant details of the warm and dense molecular gas with n(H2) $\sim$ 10$^{3-5}$ cm$^{-3}$ and Tkin $\sim$ 60 K. These data are combined with 12CO(J=1-0) and 13CO(J=1-0) results and compared with LVG calculations. We found that the ratio of 12CO(J=3-2) to 12CO(J=1-0) emission is sensitive to and is well correlated with the local Halpha flux. We interpret that differences of clump propeties represent an evolutionary sequence of GMCs in terms of density increase leading to star formation.Type I and II GMCs (starless GMCs and GMCs with HII regions only, respectively) are at the young phase of star formation where density does not yet become high enough to show active star formation and Type III GMCs (GMCs with HII regions and young star clusters) represents the later phase where the average density is increased and the GMCs are forming massive stars. The high kinetic temperature correlated with \Halpha flux suggests that FUV heating is dominant in the molecular gas of the LMC.

The NGC 1333 IRAS 4A protobinary was observed in the ammonia (2, 2) and (3, 3) lines and in the 1.3 cm continuum with a high resolution (about 1.0 arcsec). The ammonia maps show two compact sources, one for each protostar, and they are probably protostellar accretion disks. The disk associated with IRAS 4A2 is seen nearly edge-on and shows an indication of rotation. The A2 disk is brighter in the ammonia lines but dimmer in the dust continuum than its sibling disk, with the ammonia-to-dust flux ratios different by about an order of magnitude. This difference suggests that the twin disks have surprisingly dissimilar characters, one gas-rich and the other dusty. The A2 disk may be unusually active or hot, as indicated by its association with water vapor masers. The existence of two very dissimilar disks in a binary system suggests that the formation process of multiple systems has a controlling agent lacking in the isolated star formation process and that stars belonging to a multiple system do not necessarily evolve in phase with each other.

We present results of submillimeter observations of a low-mass young stellar objects (YSOs) condensation in the Cederblad 110 region of the Chamaeleon I dark cloud with Atacama Submillimeter Telescope Experiment. Our HCO+(J=4-3) map reveals a dense molecular gas with an extent of ~ 0.1 pc, which is a complex of two envelopes associated with class I sources Ced110 IRS4 and IRS11 and a very young object Cha-MMS1. The other two class I sources in this region, IRS6 and NIR89, are located outside the clump and have no extended HCO+ emission. HCO+ abundance is calculated to be 2.6 x 10^{-10} for MMS1 and 3.4 x 10^{-9} for IRS4, which are comparable to the reported value for other young sources. Bipolar outflows from IRS4 and IRS6 are detected in our 12CO(J=3-2) map. The outflow from IRS4 seems to collide with Cha-MMS1. The outflow has enough momentum to affect gas motion in MMS1, although no sign has been detected to indicate that a triggered star formation has occurred.

On the basis of near-infrared imaging observations, we derived visual extinction (Av) distribution toward ten Bok globules through measurements of both the color excess (E_{H-K}) and the stellar density at J, H, and Ks (star count). Radial column density profiles for each globule were analyzed with the Bonnor-Ebert sphere model. Using the data of our ten globules and four globules in the literature, we investigated the stability of globules on the basis of xi_max, which characterizes the Bonnor-Ebert sphere as well as the stability of the equilibrium state against the gravitational collapse. We found that more than half of starless globules are located near the critical state (xi_max = 6.5 +/- 2). Thus, we suggest that a nearly critical Bonnor-Ebert sphere characterizes the typical density structure of starless globules. Remaining starless globules show clearly unstable states (xi_max > 10). Since unstable equilibrium states are not long maintained, we expect that these globules are on the way to gravitational collapse or that they are stabilized by non-thermal support. It was also found that all the star-forming globules show unstable solutions of xi_max >10, which is consistent with the fact that they have started gravitational collapse. We investigated the evolution of a collapsing gas sphere whose initial condition is a nearly critical Bonnor-Ebert sphere. We found that the column density profiles of the collapsing sphere mimic those of the static Bonnor-Ebert spheres in unstable equilibrium. The collapsing gas sphere resembles marginally unstable Bonnor-Ebert spheres for a long time. We found that the frequency distribution of xi_max for the observed starless globules is consistent with that from model calculations of the collapsing sphere.

Lynds 1498 is a pre-protostellar core (PPC) and was one of the initial objects toward which molecular depletion and differentiation was detected. Despite the considerable scrutiny of L1498, there has not been a extensive study of the density and temperature structure as derived from radiative transfer modeling of dust continuum observations. We present deep SCUBA observations of L1498 at 850 and 450 micron, high resolution BEARS maps of the N2H+ 1-0 transition, CSO observations of the N2H+ 3-2 transition, and GBT observations of the C3S 4-3 transition. We also present a comparison of derived properties between L1498 and nearby PPCs that have been observed at far-infrared and submillimeter wavelengths. We present a more realistic treatment of PPC heating which varies the strength of the ISRF, Sisrf, and includes attenuation of the ISRF due to dust grains at the outer radius of the core, Av. The best-fitted model consists of a Bonner-Ebert sphere with a central density of 1 - 3 x 10^4 cm-3, R_o ~ 0.29 pc, 0.5 <= Sisrf <= 1, Av ~ 1 mag, and a nearly isothermal temperature profile of ~ 10.5 K for OH8 opacities. C3S emission shows a central depletion hole while N2H+ emission is centrally peaked. The observed depletions of C3S and H2CO, the modest N2H+ abundance, and a central density that is an order of magnitude lower than other modeled PPCs suggests that L1498 may be a forming PPC. Our derived temperature and density profile will improve modeling of molecular line observations that will explicate the core's kinematical and chemical state. (abridged)

N2H+ observations of molecular cloud cores in Taurus with the Nobeyama 45 m radio telescope are reported. We compare ``cores with young stars'' with ``cores without young stars''. The differences in core radius, linewidth, and core mass are small. Linewidth is dominated by thermal motions in both cases. N2H+ maps show that the intensity distribution does not differ much between cores without stars and those with stars. This is in contrast to the result previously obtained in H13CO+ toward Taurus molecular cloud cores. Larger degree of depletion of H13CO+ in starless cores will be one possible explanation for this difference. We studied the physical state of molecular cloud cores in terms of ``critical pressure'' for the surface (external) pressure. There is no systematic difference between starless cores and cores with stars in this analysis. Both are not far from the critical state for pressure equilibrium. We suggest that molecular cloud cores in which thermal support is dominated evolve toward star formation by keeping close to the critical state. This result is in contrast with that obtained in the intermediate-mass star forming region OMC-2/3, where molecular cloud cores evolve by decreasing the critical pressure appreciably. We investigate the radial distribution of the integrated intensity. Cores with stars are found to have shallow (-1.8 to -1.6) power-law density profiles.

We have compared the intensity distribution of molecular line emission with that of dust continuum emission, and modeled molecular line profiles in three different preprotostellar cores in order to test how dynamical evolution is related to chemical evolution, and whether we can use different chemical tracers to identify specific dynamical evolutionary stages. We used dust continuum emission to obtain the input density and temperature structures by calculating radiative transfer of dust emission. Our results show that chemical evolution is dependent on dynamical processes, which can give different evolutionary timescales, as well as the density structure of the core.

We have compared molecular line emission to dust continuum emission and modeled molecular lines using Monte Carlo simulations in order to study the depletion of molecules and the ionization fraction in three preprotostellar cores, L1512, L1544, and L1689B. L1512 is much less dense than L1544 and L1689B, which have similar density structures. L1689B has a different environment from those of L1512 and L1544. We used density and temperature profiles, calculated by modeling dust continuum emission in the submillimeter, for modeling molecular line profiles. In addition, we have used molecular line profiles and maps observed in several different molecules toward the three cores. We find a considerable diversity in chemical state among the three cores. The molecules include those sensitive to different timescales of chemical evolution such as CCS, the isotopes of CO and HCO+, DCO+, and N2H+. The CO molecule is significantly depleted in L1512 and L1544, but not in L1689B. CCS may be in the second enhancement of its abundance in L1512 and L1544 because of the significant depletion of CO molecules. N2H+ might already start to be depleted in L1512, but it traces very well the distribution of dust emission in L1544. On the other hand, L1689B may be so young that N2H+ has not reached its maximum yet. The ionization fraction has been calculated using H13CO+ and DCO+. This study suggests that chemical evolution depends on the absolute timescale during which a core stays in a given environment as well as its density structure.

We present observations of the 3P1-3P0 fine structure transition of atomic carbon [CI], the J=3-2 transition of CO, as well as of the J=1-0 transitions of 13CO and C18O toward DR15, an HII region associated with two mid-infrared dark clouds (IRDCs). The 13CO and C18O J=1-0 emissions closely follow the dark patches seen in optical wavelength, showing two self-gravitating molecular cores with masses of 2000 Msun and 900 Msun, respectively, at the positions of the catalogued IRDCs. Our data show a rough spatial correlation between [CI] and 13CO J=1-0. Bright [CI] emission occurs in relatively cold gas behind the molecular cores, neither in highly excited gas traced by CO J=3-2 emission nor in HII region/molecular cloud interface. These results are inconsistent with those predicted by standard photodissociation region (PDR) models, suggesting an origin for interstellar atomic carbon unrelated to photodissociation processes.

We have observed the Orion Molecular Clouds 2 and 3 (OMC-2 and OMC-3) with the Chandra X-ray Observatory (CXO). The northern part of OMC-3 is found to be particularly rich in new X-ray features; four hard X-ray sources are located in and along the filament of cloud cores. Two sources coincide positionally with the sub$mm$-$mm$ dust condensations of MMS 2 and 3 or an outflow radio source VLA 1, which are in a very early phase of star formation. The X-ray spectra of these sources show an absorption column of (1-3) x 10^23 H cm-2. Assuming a moderate temperature plasma, the X-ray luminosity in the 0.5-10 keV band is estimated to be ~10^30 erg s^-1 at a distance of 450 pc. From the large absorption, positional coincidence and moderate luminosity, we infer that the hard X-rays are coming from very young stellar objects embedded in the molecular cloud cores. We found another hard X-ray source near the edge of the dust filament. The extremely high absorption of 3 x 10^23 H cm^-2 indicates that the source must be surrounded by dense gas, suggesting that it is either a YSO in an early accretion phase or a Type II AGN (e.g. a Seyfert 2), although no counterpart is found at any other wavelength. In contrast to the hard X-ray sources, soft X-ray sources are found spread around the dust filaments, most of which are identified with IR sources in the T Tauri phase.

Large scale mapping observations of the 3P1-3P0 fine structure transition of atomic carbon (CI, 492 GHz) and the J=3-2 transition of CO (346 GHz) toward the Orion A molecular cloud have been carried out with the Mt. Fuji submillimeter-wave telescope. The observations cover 9 square degrees, and include the Orion nebula M42 and the L1641 dark cloud complex. The CI emission extends over almost the entire region of the Orion A cloud and is surprisingly similar to that of 13CO(J=1-0).The CO(J=3-2) emission shows a more featureless and extended distribution than CI.The CI/CO(J=3-2) integrated intensity ratio shows a spatial gradient running from the north (0.10) to the south (1.2) of the Orion A cloud, which we interpret as a consequence of the temperature gradient. On the other hand, the CI/13CO(J=1-0) intensity ratio shows no systematic gradient. We have found a good correlation between the CI and 13CO(J=1-0) intensities over the Orion A cloud. This result is discussed on the basis of photodissociation region models.

We have mapped the C I 3P1-3P0 line at 492 GHzin three molecular clouds immersed in weak ultraviolet radiation fields, TMC-1, L134N, and IC 5146. In all three clouds, the CI peak TA* ~ 1 K, with very small dispersion. The spatial C I distribution is extended and rather smooth. The J = 2-1 transitions of CO isotopomers were observed at the same angular resolution as C I. The C I peak TA* is typically a third of the peak TA* of 13CO J = 2-1, and the C I emission is usually more extended than emission in 13CO or C18O J=2-1. The C I linewidth is close to the 13CO J = 2-1 linewidth, larger than the C18O J = 2-1 line width, and smaller than the 12CO J = 2-1 linewidth. The uniformity of the C I peak TA* is remarkable for a line in the Wien portion of the Planck function and indicates a very uniform excitation temperature. This uniformity is best explained if the line if opaque and thermalized. If so, the CI line probes kinetic temperature in clouds exposed to low ultraviolet fluxes. This conclusion has significant implications for the thermal balance in such clouds. At Av ~ 2, these clouds have a remarkably constant temperature from place to place and from cloud to cloud (7.9+/-0.8 K). Photodissociation region models of clouds immersed in the mean interstellar radiation field tend to predict stronger lines than we see, but this may be an artifact of assumptions about the temperature.