ABSTRACT
ATLASGAL is an 870-µm dust survey of 420 deg2 the inner Galactic plane and has been used to identify ∼10 000 dense molecular clumps. Dedicated follow-up observations and complementary ...surveys are used to characterize the physical properties of these clumps, map their Galactic distribution, and investigate the evolutionary sequence for high-mass star formation. The analysis of the ATLASGAL data is ongoing: We present an up-to-date version of the catalogue. We have classified 5007 clumps into four evolutionary stages (quiescent, protostellar, young stellar objects and H ii regions) and find similar numbers of clumps in each stage, suggesting a similar lifetime. The luminosity-to-mass (Lbol/Mfwhm) ratio curve shows a smooth distribution with no significant kinks or discontinuities when compared to the mean values for evolutionary stages indicating that the star formation process is continuous and that the observational stages do not represent fundamentally different stages or changes in the physical mechanisms involved. We compare the evolutionary sample with other star formation tracers (methanol and water masers, extended green objects and molecular outflows) and find that the association rates with these increases as a function of evolutionary stage, confirming that our classification is reliable. This also reveals a high association rate between quiescent sources and molecular outflows, revealing that outflows are the earliest indication that star formation has begun and that star formation is already ongoing in many of the clumps that are dark even at 70 µm.
Context.
Understanding the details of the formation process of massive (i.e.
M
≳ 8–10
M
⊙
) stars is a long-standing problem in astrophysics. They form and evolve very quickly, and almost their ...entire formation process takes place deeply embedded in their parental clumps. Together with the fact that these objects are rare and at a relatively large distance, this makes observing them very challenging.
Aims.
We present a method for deriving accurate timescales of the evolutionary phases of the high-mass star formation process.
Methods.
We modelled a representative number of massive clumps of the ATLASGAL-TOP100 sample that cover all the evolutionary stages. The models describe an isothermal collapse and the subsequent warm-up phase, for which we followed the chemical evolution. The timescale of each phase was derived by comparing the results of the models with the properties of the sources of the ATLASGAL-TOP100 sample, taking into account the mass and luminosity of the clumps, and the column densities of methyl acetylene (CH
3
CCH), acetonitrile (CH
3
CN), formaldehyde (H
2
CO), and methanol (CH
3
OH).
Results.
We find that the molecular tracers we chose are affected by the thermal evolution of the clumps, showing steep ice evaporation gradients from 10
3
to 10
5
AU during the warm-up phase. We succeed in reproducing the observed column densities of CH
3
CCH and CH
3
CN, but H
2
CO and CH
3
OH agree less with the observed values. The total (massive) star formation time is found to be ~5.2 × 10
5
yr, which is defined by the timescales of the individual evolutionary phases of the ATLASGAL-TOP100 sample: ~5 × 10
4
yr for 70-μm weak, ~1.2 × 10
5
yr for mid-IR weak, ~2.4 × 10
5
yr for mid-IR bright, and ~1.1 × 10
5
yr for HII-region phases.
Conclusions.
With an appropriate selection of molecular tracers that can act as chemical clocks, our model allows obtaining robust estimates of the duration of the individual phases of the high-mass star formation process. It also has the advantage of being capable of including additional tracers aimed at increasing the accuracy of the estimated timescales.
ABSTRACT
An estimate of the degree of CO-depletion (fD) provides information on the physical conditions occurring in the innermost and densest regions of molecular clouds. A key parameter in these ...studies is the size of the depletion radius, i.e. the radius within which the C-bearing species, and in particular CO, are largely frozen on to dust grains. A strong depletion state (i.e. fD > 10, as assumed in our models) is highly favoured in the innermost regions of dark clouds, where the temperature is <20 K and the number density of molecular hydrogen exceeds a few × 104 cm−3. In this work, we estimate the size of the depleted region by studying the Infrared Dark Cloud (IRDC) G351.77−0.51. Continuum observations performed with the Herschel Space Observatory and the LArge APEX BOlometer CAmera, together with APEX C18O and C17O J = 2→1 line observations, allowed us to recover the large-scale beam- and line-of-sight-averaged depletion map of the cloud. We built a simple model to investigate the depletion in the inner regions of the clumps in the filament and the filament itself. The model suggests that the depletion radius ranges from 0.02 to 0.15 pc, comparable with the typical filament width (i.e. ∼0.1 pc). At these radii, the number density of H2 reaches values between 0.2 and 5.5 × 105 cm−3. These results provide information on the approximate spatial scales on which different chemical processes operate in high-mass star-forming regions and also suggest caution when using CO for kinematical studies in IRDCs.
Disentangling the different stages of the star formation process, in particular in the high-mass regime, is a challenge in astrophysics. Chemical clocks could help alleviate this problem, but their ...evolution strongly depends on many parameters, leading to degeneracy in the interpretation of the observational data. One of these uncertainties is the degree of CO depletion. We present here the first self-consistent magneto-hydrodynamic simulations of high-mass, star-forming regions at different scales, fully coupled with a nonequilibrium chemical network, which includes C-N-O bearing molecules. Depletion and desorption processes are treated time dependently. The results show that full CO depletion (i.e., all gas-phase CO frozen-out on the surface of dust grains) can be reached very quickly, in one-third or even smaller fractions of the freefall time, whether the collapse proceeds on slow or fast timescales. This leads to a high level of deuteration in a short time, both for typical tracers like N2H+, as well as for the main ion H 3 + , the latter being in general larger and more extended. N2 depletion is slightly less efficient, and no direct effects on N-bearing molecules and deuterium fractionation are observed. We show that CO depletion is not the only driver of deuteration, and that there is a strong impact on Dfrac when changing the grain size. We finally apply a two-dimensional Gaussian point-spread function to our results to mimic observations with single-dish and interferometers. Our findings suggest that the low-values observed in high-mass star-forming clumps are in reality masking a full-depletion stage in the inner 0.1 pc region.
Context. Deuteration has been suggested to be a reliable chemical clock of star-forming regions due to its strong dependence on density and temperature changes during cloud contraction. In ...particular, the H3+ isotopologues (e.g. ortho-H2D+) seem to act as good proxies of the evolutionary stages of the star formation process. While this has been widely explored in low-mass star-forming regions, in the high-mass counterparts only a few studies have been pursued, and the reliability of deuteration as a chemical clock remains inconclusive. Aims. We present a large sample of o-H2D+ observations in high-mass star-forming regions and discuss possible empirical correlations with relevant physical quantities to assess its role as a chronometer of star-forming regions through different evolutionary stages. Methods. APEX observations of the ground-state transition of o-H2D+ were analysed in a large sample of high-mass clumps selected from the ATLASGAL survey at different evolutionary stages. Column densities and beam-averaged abundances of o-H2D+ with respect to H2, X(o-H2D+), were obtained by modelling the spectra under the assumption of local thermodynamic equilibrium. Results. We detect 16 sources in o-H2D+ and find clear correlations between X(o-H2D+) and the clump bolometric luminosity and the dust temperature, while only a mild correlation is found with the CO-depletion factor. In addition, we see a clear correlation with the luminosity-to-mass ratio, which is known to trace the evolution of the star formation process. This would indicate that the deuterated forms of H3+ are more abundant in the very early stages of the star formation process and that deuteration is influenced by the time evolution of the clumps. In this respect, our findings would suggest that the X(o-H2D+) abundance is mainly affected by the thermal changes rather than density changes in the gas. We have employed these findings together with observations of H13CO+, DCO+, and C17O to provide an estimate of the cosmic-ray ionisation rate in a sub-sample of eight clumps based on recent analytical work. Conclusions. Our study presents the largest sample of o-H2D+ in star-forming regions to date. The results confirm that the deuteration process is strongly affected by temperature and suggests that o-H2D+ can be considered a reliable chemical clock during the star formation processes, as proved by its strong temporal dependence.
Analyzing the kinematics of filamentary molecular clouds is a crucial step toward understanding their role in the star formation process. Therefore, we study the kinematics of 283 filament candidates ...in the inner Galaxy, that were previously identified in the ATLASGAL dust continuum data. The 13CO(2 – 1) and C18O(2 – 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of the Inner Galactic Inter Stellar Medium) allows us to analyze the kinematics of these targets and to determine their physical properties at a resolution of 30′′ and 0.25 km s−1. To do so, we developed an automated algorithm to identify all velocity components along the line-of-sight correlated with the ATLASGAL dust emission, and derive size, mass, and kinematic properties for all velocity components. We find two-third of the filament candidates are coherent structures in position-position-velocity space. The remaining candidates appear to be the result of a superposition of two or three filamentary structures along the line-of-sight. At the resolution of the data, on average the filaments are in agreement with Plummer-like radial density profiles with a power-law exponent of p ≈ 1.5 ± 0.5, indicating that they are typically embedded in a molecular cloud and do not have a well-defined outer radius. Also, we find a correlation between the observed mass per unit length and the velocity dispersion of the filament of m ∝ σv2 $m \propto \sigma_{\rm{v}}^2$ m∝ σ v 2 . We show that this relation can be explained by a virial balance between self-gravity and pressure. Another possible explanation could be radial collapse of the filament, where we can exclude infall motions close to the free-fall velocity.
Context. The assumption of a gas-to-dust mass ratio γ is a common approach to estimate the basic properties of molecular clouds, such as total mass and column density of molecular hydrogen, from ...(sub)mm continuum observations of the dust. In the Milky Way a single value is used at all galactocentric radii, independently of the observed metallicity gradients. Both models and extragalactic observations suggest that this quantity increases for decreasing metallicity Z, typical of the outer regions in disks, where fewer heavy elements are available to form dust grains. Aims. We aim to investigate the variation of the gas-to-dust ratio as a function of galactocentric radius and metallicity, to allow a more accurate characterisation of the quantity of molecular gas across the galactic disk, as derived from observations of the dust. Methods. Observations of the optically thin C18O (2–1) transition were obtained with the APEX telescope for a sample of 23 massive and dense star-forming regions in the far outer Galaxy (galactocentric distance greater than 14 kpc). From the modelling of this line and of the spectral energy distribution of the selected clumps we computed the gas-to-dust ratio and compared it to that of well-studied sources from the ATLASGAL TOP100 sample in the inner galactic disk. Results. The gradient in γ is found to be 0.087+0.047-0.025 dex kpc-1 (or equivalently γ ∝ Z-1.4+0.3-1.0). The dust-to-metal ratio, decreases with galactocentric radius, which is the most common situation also for external late-type galaxies. This suggests that grain growth dominates over destruction. The predicted γ is in excellent agreement with the estimates in Magellanic clouds, for the appropriate value of Z.
Context. Massive-star formation and the processes involved are still poorly understood. The ATLASGAL survey provides an ideal basis for detailed studies of large numbers of massive-star forming ...clumps covering the whole range of evolutionary stages. The ATLASGAL Top100 is a sample of clumps selected by their infrared and radio properties to be representative for the whole range of evolutionary stages. Aims. The ATLASGAL Top100 sources are the focus of a number of detailed follow-up studies that will be presented in a series of papers. In the present work we use the dust continuum emission to constrain the physical properties of this sample and identify trends as a function of source evolution. Methods. We determine flux densities from mid-infrared to submillimeter wavelength (8-870 mu m) images and use these values to fit their spectral energy distributions and determine their dust temperature and flux. Combining these with recent distances from the literature including maser parallax measurements we determine clump masses, luminosities and column densities. Results. We define four distinct source classes from the available continuum data and arrange these into an evolutionary sequence. This begins with sources found to be dark at 70 mu m, followed by 24 mu m weak sources with an embedded 70 mu m source, continues through mid-infrared bright sources and ends with infrared bright sources associated with radio emission (i.e., Hii regions). We find trends for increasing temperature, luminosity, and column density with the proposed evolution sequence, confirming that this sample is representative of different evolutionary stages of massive star formation. Our sources span temperatures from approximately 11 to 41 K, with bolometric luminosities in the range 57 L sub(middot in circle)-3.8 x 10 super(6)L sub(middot in circle). The highest masses reach 4.3 x 10 super(4)M sub(middot in circle) and peak column densities up to 1.1 x 10 super(24) cm super(-1), and therefore have the potential to form the most massive O-type stars. We show that at least 93 sources (85%) of this sample have the ability to form massive stars and that most are gravitationally unstable and hence likely to be collapsing. Conclusions. The highest column density ATLASGAL sources cover the whole range of evolutionary stages from the youngest to the most evolved high-mass-star forming clumps. Study of these clumps provides a unique starting point for more in-depth research on massive-star formation in four distinct evolutionary stages whose well defined physical parameters afford more detailed studies. As most of the sample is closer than 5 kpc, these sources are also ideal for follow-up observations with high spatial resolution.