We present a multi-wavelength study (from X-ray to millimetre) of the nearby low-luminosity active galactic nucleus (LLAGN) NGC 7213. We combine the information from the different bands to characterize the source in terms of contribution from the AGN and the host-galaxy interstellar medium (ISM). This approach allows us to provide a coherent picture of the role of the AGN and its impact, if any, on the star formation and molecular gas properties of the host galaxy. We focused our study on archival ALMA Cycle 1 observations, where the CO(2-1) emission line has been used as a tracer of the molecular gas. Using the 3DBarolo code on ALMA data, we performed the modelling of the molecular gas kinematics traced by the CO(2-1) emission, finding a rotationally dominated pattern. The host-galaxy molecular gas mass was estimated from the integrated CO(2-1) emission line obtained with APEX data, assuming an $\alpha_{CO}$ conversion factor. By using the ALMA data, we would have underestimated the gas masses by a factor $\sim$3, given the filtering out of the large scale emission in interferometric observations. We also performed a complete X-ray spectral analysis on archival observations, revealing a relatively faint and unobscured AGN. The AGN results to be too faint to significantly affect the properties of the host-galaxy, such as star formation activity and molecular gas kinematics and distribution.

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 alleviating 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 non-equilibrium 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 free-fall 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 N$_2$H$^+$, as well as for the main ion H$_3^+$, the latter being in general larger and more extended. N$_2$ 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 $D_{frac}$ 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.

An estimate of the degree of CO-depletion ($f_D$) 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 onto dust grains. A strong depletion state (i.e. $f_D>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 $\times$10$^{4}$ 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 C$^{18}$O and C$^{17}$O J=2$\rightarrow$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.$\sim$0.1 pc). At these radii, the number density of H$_2$ reaches values between 0.2 and 5.5$\times$10$^{5}$ 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.

Young massive stars are usually found embedded in dense and massive molecular clumps and are known for being highly obscured and distant. During their formation process, deuteration is regarded as a potentially good indicator of the formation stage. Therefore, proper observations of such deuterated molecules are crucial, but still, hard to perform. In this work, we test the observability of the transition o-H$_2$D$^+(1_{10}$-$1_{11})$, using a synthetic source, to understand how the physical characteristics are reflected in observations through interferometers and single-dish telescopes. In order to perform such tests, we post-processed a magneto-hydrodynamic simulation of a collapsing magnetized core using the radiative transfer code POLARIS. Using the resulting intensity distributions as input, we performed single-dish (APEX) and interferometric (ALMA) synthetic observations at different evolutionary times, always mimicking realistic configurations. Finally, column densities were derived to compare our simulations with real observations previously performed. Our derivations for o-H$_2$D$^+$ are in agreement with values reported in the literature, in the range of 10$^{\!10-11}$cm$^{\!-2}$ and 10$^{\!12-13}$cm$^{\!-2}$ for single-dish and interferometric measurements, respectively.

ALMA products are stored in the Science Archive in the form of FITS images. It is a common idea that the FITS image headers should collect in their keywords all the information that an archive User might want to search for in order to quickly select, compare, or discard datasets. With this perspective in mind, we first present a short description of the current status of the ALMA FITS archive and images. We realized that at the moment most of the parameters that could be useful for a general User are still missing in the archived data. We then provide a CASA task generating the image header keywords that we suggest to be relevant for the scientific exploitation of the ALMA archival data. The proposed tool could be also applied to several types of interferometer data and part of it is implemented in a web interface. An example of the scientific application of the keywords is also discussed.

We have conducted a 3-mm molecular-line survey towards 570 high-mass star-forming clumps, using the Mopra telescope. The sample is selected from the 10,000 clumps identified by the ATLASGAL survey and includes all of the most important embedded evolutionary stages associated with massive star formation, classified into five distinct categories (quiescent, protostellar, young stellar objects, \hii\ regions and photo-dominated regions). The observations were performed in broadband mode with frequency coverage of 85.2 to 93.4\,GHz and a velocity resolution of $\sim$0.9\,\kms, detecting emission from 26 different transitions. We find significant evolutionary trends in the detection rates, integrated line intensities, and abundances of many of the transitions and also identify a couple of molecules that appear to be invariant to changes in the dust temperature and evolutionary stage (N$_2$H$^+$\,(1-0) and HN$^{13}$C\,(1-0)). We use the K-ladders for CH$_3$C$_2$H\,(5-4) and CH$_3$CH\,(5-4) to calculate the rotation temperatures and find $\sim$1/3 of the quiescent clumps have rotation temperatures that suggest the presence of an internal heating source. These sources may constitute a population of very young protostellar objects that are still dark at 70\,\mum\ and suggest that the fraction of truly quiescent clumps may only be a few per cent. We also identify a number of line ratios that show a strong correlation with the evolutionary stage of the embedded objects and discuss their utility as diagnostic probes of evolution.

The diffuse and translucent molecular clouds traced in absorption along the line of sight to strong background sources have so far been investigated mainly in the spectral domain because of limited angular resolution or small sizes of the background sources. We aim to resolve and investigate the spatial structure of molecular clouds traced by several molecules detected in absorption along the line of sight to SgrB2(N). We have used spectral line data from the EMoCA survey performed with ALMA, taking advantage of the high sensitivity and angular resolution. We identify, on the basis of the spectral analysis of c-C3H2 across the field of view, 15 main velocity components along the line of sight to SgrB2(N) and several components in the envelope of SgrB2. The c-C3H2 column densities reveal two categories of clouds. Clouds in Category I (3 kpc arm, 4 kpc arm, and some GC clouds) have smaller c-C3H2 column densities, smaller linewidths, and smaller widths of their column density PDFs than clouds in Category II (Scutum arm, Sgr arm, and other GC clouds). To investigate the spatial structure we derive opacity maps for the following molecules: c-C3H2, H13CO+, 13CO, HNC, HN13C, HC15N, CS, C34S, 13CS, SiO, SO, and CH3OH. These maps reveal that most molecules trace relatively homogeneous structures that are more extended than the field of view defined by the background continuum emission (about 15", that is 0.08-0.6pc depending on the distance). SO and SiO show more complex structures with smaller clumps of size ~5-8". Our analysis suggests that the driving of the turbulence is mainly solenoidal in the investigated clouds. On the basis of HCO+, we conclude that most line-of-sight clouds towards SgrB2 are translucent, including all clouds where complex organic molecules were recently detected. We also conclude that CCH and CH are good probes of H2 in both diffuse and translucent clouds.

Context: In cold and dense gas prior to the formation of young stellar objects, heavy molecular species (including CO) are accreted onto dust grains. Under these conditions H$_3^+$ and its deuterated isotopologues become more abundant, enhancing the deuterium fraction of molecules such as N$_2$H$^+$ that are formed via ion-neutral reactions. Because this process is extremely temperature sensitive, the abundance of these species is likely linked to the evolutionary stage of the source. Aims: We investigate how the abundances of o-H$_2$D$^+$ and N$_2$D$^+$ vary with evolution in high-mass clumps. Methods: We observed with APEX the ground-state transitions of o-H$_2$D$^+$ near 372 GHz, and N$_2$D$^+$(3-2) near 231 GHz for three massive clumps in different evolutionary stages. The sources were selected within the G351.77-0.51 complex to minimise the variation of initial chemical conditions, and to remove distance effects. We modelled their dust continuum emission to estimate their physical properties, and also modelled their spectra under the assumption of local thermodynamic equilibrium to calculate beam-averaged abundances. Results: We find an anticorrelation between the abundance of o-H$_2$D$^+$ and that of N$_2$D$^+$, with the former decreasing and the latter increasing with evolution. With the new observations we are also able to provide a qualitative upper limit to the age of the youngest clump of about 10$^5$ yr, comparable to its current free-fall time. Conclusions: We can explain the evolution of the two tracers with simple considerations on the chemical formation paths, depletion of heavy elements, and evaporation from the grains. We therefore propose that the joint observation and the relative abundance of o-H$_2$D$^+$ and N$_2$D$^+$ can act as an efficient tracer of the evolutionary stages of the star-formation process.

High-mass stars are formed within massive molecular clumps, where a large number of stars form close together. The evolution of the clumps with different masses and luminosities is mainly regulated by its high-mass stellar content and the formation of such objects is still not well understood. In this work, we characterise the mid-J CO emission in a statistical sample of 99 clumps (Top100) selected from the ATLASGAL survey that are representative of the Galactic proto-cluster population. High-spatial resolution APEX-CHAMP+ maps of the CO(6-5) and CO(7-6) transitions were obtained and combined with additional single-pointing APEX-FLASH+ spectra of the CO(4-3) line. We study the correlations of the CO line luminosities and profiles for the three CO transitions with the clump properties and investigate if and how they change as a function of the evolution. All sources were detected above 3-$\sigma$ in all three CO transitions and most of the sources exhibit broad CO emission likely associated with molecular outflows. We found that the extension of the mid-J CO emission is correlated with the size of the dust emission traced by the Herschel-PACS 70 $\mu$m maps. The CO line luminosity is correlated with the luminosity and mass of the clumps. However, it does not correlate with the L/M ratio. The dependency of the CO luminosity with the properties of the clumps is steeper for higher-J transitions. Our data seem to exclude that this trend is biased by self-absorption features in the CO emission, but rather suggest that different J transitions arise from different regions of the inner envelope. Moreover, high-mass clumps show similar trends in CO luminosity as lower mass clumps, but are systematically offset towards larger values, suggesting that higher column density and/or temperature (of unresolved) CO emitters are found inside high-mass clumps.

G351.776-0.527 is among the most massive, closest, and youngest filaments in the inner Galactic plane and therefore it is an ideal laboratory to study the kinematics of dense gas and mass replenishment on a large scale. In this paper, we present far-infrared (FIR) and submillimetre wavelength continuum observations combined with spectroscopic C$^{18}$O (2-1) data of the entire region to study its temperature, mass distribution, and kinematics. The structure is composed of a main elongated region with an aspect ratio of $\sim 23$, which is associated with a network of filamentary structures. The main filament has a remarkably constant width of 0.2 pc. The total mass of the network (including the main filament) is $\geq 2600$ M$_\odot$, while we estimate a mass of $\sim 2000$ M$_\odot$ for the main structure. Therefore, the network harbours a large reservoir of gas and dust that could still be accreted onto the main structure. From the analysis of the gas kinematics, we detect two velocity components in the northern part of the main filament. The data also reveal velocity oscillations in C$^{18}$O along the spine in the main filament and in at least one of the branches. Considering the region as a single structure, we find that it is globally close to virial equilibrium indicating that the entire structure is approximately in a stable state.

Analysing the kinematics of filamentary molecular clouds is a crucial step towards 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 $^{13}$CO(2 - 1) and C$^{18}$O(2 - 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of the Inner Galactic Inter Stellar Medium) allows us to analyse the kinematics of these targets and to determine their physical properties at a resolution of 30 arcsec and 0.25 km/s. 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 \sim \sigma_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.

We present a catalogue of ALMA flux density measurements of 754 calibrators observed between August 2012 and September 2017, for a total of 16,263 observations in different bands and epochs. The flux densities were measured reprocessing the ALMA images generated in the framework of the ALMACAL project, with a new code developed by the Italian node of the European ALMA Regional Centre. A search in the online databases yielded redshift measurements for 589 sources ($\sim$78 per cent of the total). Almost all sources are flat-spectrum, based on their low-frequency spectral index, and have properties consistent with being blazars of different types. To illustrate the properties of the sample we show the redshift and flux density distributions as well as the distributions of the number of observations of individual sources and of time spans in the source frame for sources observed in bands 3 (84$-$116 GHz) and 6 (211$-$275 GHz). As examples of the scientific investigations allowed by the catalogue we briefly discuss the variability properties of our sources in ALMA bands 3 and 6 and the frequency spectra between the effective frequencies of these bands. We find that the median variability index steadily increases with the source-frame time lag increasing from 100 to 800 days, and that the frequency spectra of BL Lacs are significantly flatter than those of flat-spectrum radio quasars. We also show the global spectral energy distributions of our sources over 17 orders of magnitude in frequency.

Fragmentation of massive dense molecular clouds is the starting point in the formation of rich clusters and massive stars. Theory and numerical simulations indicate that the population of the fragments (number, mass, diameter, separation) resulting from the gravitational collapse of such clumps is probably regulated by the balance between the magnetic field and the other competitors of self-gravity, in particular turbulence and protostellar feedback. We have observed 11 massive, dense and young star-forming clumps with the Atacama Large Millimeter Array (ALMA) in the thermal dust continuum emission at $\sim 1$~mm with an angular resolution of 0.25 arcseconds with the aim of determining their population of fragments. We find fragments on sub-arcsecond scales in 8 out of the 11 sources. The ALMA images indicate two different fragmentation modes: a dominant fragment surrounded by companions with much smaller mass and size, and many ($\geq 8$) fragments with a gradual change in masses and sizes. On average, the largest number of fragments is found towards the warmer and more massive clumps. Also, the warmer clumps tend to form fragments with larger mass and size. To understand the role of the different physical parameters to regulate the final population of the fragments, we have simulated the collapse of a massive clump of $100$ and $300$ M$_{\odot}$ having different magnetic support. The simulations indicate that: (1) fragmentation is inhibited when the initial turbulence is low, independent of the other physical parameters. (2) a filamentary distribution of the fragments is favoured in a highly magnetised clump. We conclude that the clumps that show many fragments distributed in a filamentary-like structure are likely characterised by a strong magnetic field, while the others are possible also in a weaker magnetic field.

The local Seyfert 2 galaxy NGC5135, thanks to its almost face-on appearance, a bulge overdensity of stars, the presence of a large-scale bar, an AGN and a Supernova Remnant, is an excellent target to investigate the dynamics of inflows, outflows, star formation and AGN feedback. Here we present a reconstruction of the gas morphology and kinematics in the inner regions of this galaxy, based on the analysis of Atacama Large Millimeter Array (ALMA) archival data. To our purpose, we combine the available $\sim$100 pc resolution ALMA 1.3 and 0.45 mm observations of dust continuum emission, the spectroscopic maps of two transitions of the CO molecule (tracer of molecular mass in star forming and nuclear regions), and of the CS molecule (tracer of the dense star forming regions) with the outcome of the SED decomposition. By applying the $^{\rm 3D}$BAROLO software (3D-Based Analysis of Rotating Object via Line Observations), we have been able to fit the galaxy rotation curves reconstructing a 3D tilted-ring model of the disk. Most of the observed emitting features are described by our kinematic model. We also attempt an interpretation for the emission in few regions that the axisymmetric model fails to reproduce. The most relevant of these is a region at the northern edge of the inner bar, where multiple velocity components overlap, as a possible consequence of the expansion of a super-bubble.

We aim to directly determine the kinetic temperature and spatial density with formaldehyde for the $\sim$100 brightest ATLASGAL-selected clumps at 870 $\mu$m representing various evolutionary stages of high-mass star formation. Ten transitions ($J$ = 3-2 and 4-3) of ortho- and para-H$_2$CO near 211, 218, 225, and 291 GHz were observed with the APEX 12 m telescope. Using non-LTE models with RADEX, we derive the gas kinetic temperature and spatial density using the measured p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$, 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$, and 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ ratios. The gas kinetic temperatures derived from the p-H$_2$CO 3$_{21}$-2$_{20}$/3$_{03}$-2$_{02}$ and 4$_{22}$-3$_{21}$/4$_{04}$-3$_{03}$ line ratios are high, ranging from 43 to $>$300 K with an unweighted average of 91 $\pm$ 4 K. Deduced $T_{\rm kin}$ values from the $J$ = 3-2 and 4-3 transitions are similar. Spatial densities of the gas derived from the p-H$_2$CO 4$_{04}$-3$_{03}$/3$_{03}$-2$_{02}$ line ratios yield 0.6-8.3 $\times$ 10$^6$ cm$^{-3}$ with an unweighted average of 1.5 ($\pm$0.1) $\times$ 10$^6$ cm$^{-3}$. A comparison of kinetic temperatures derived from p-H$_2$CO, NH$_3$, and the dust emission indicates that p-H$_2$CO traces a distinctly higher temperature than the NH$_3$ (2,2)/(1,1) transitions and the dust, tracing heated gas more directly associated with the star formation process. The H$_2$CO linewidths are found to be correlated with bolometric luminosity and increase with the evolutionary stage of the clumps, which suggests that higher luminosities tend to be associated with a more turbulent molecular medium. It seems that the spatial densities measured with H$_2$CO do not vary significantly with the evolutionary stage of the clumps. However, averaged gas kinetic temperatures derived from H$_2$CO increase with time through the evolution of the clumps.

(Abridged) Context: The assumption of a gas-to-dust mass ratio (\gamma) 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 C\$^{18}\$O (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 gas-to-dust ratio is found to be 0.087 dex/kpc (or equivalently \gamma\ proportional to Z\$^{-1.4}\$). 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 gas-to-dust ratio is in excellent agreement with the estimates in Magellanic clouds, for the appropriate value of Z.

Deuterium fractionation is considered as an important process to infer the chemical ages of prestellar cores in filaments. We present here the first magneto-hydrodynamical simulations including a chemical network to study deuterium fractionation in magnetized and turbulent filaments and their substructures. The filaments typically show widespread deuterium fractionation with average values $\gtrsim0.01$. For individual cores of similar age, we observe the deuteration fraction to increase with time, but also to be independent of their average properties such as density, virial or mass-to-magnetic flux ratio. We further find a correlation of the deuteration fraction with core mass, average H$_2$ density and virial parameter only at late evolutionary stages of the filament and attribute this to the lifetime of the individual cores. Specifically, chemically old cores reveal higher deuteration fractions. Within the radial profiles of selected cores, we notice differences in the structure of the deuteration fraction or surface density, which we can attribute to their different turbulent properties. High deuteration fractions of the order $0.01-0.1$ may be reached within approximately $200$~kyrs, corresponding to two free-fall times, as defined for cylindrical systems, of the filaments

Abridged: ATLASGAL is an unbiased 870 micron submillimetre survey of the inner Galactic plane. It provides a large and systematic inventory of all massive, dense clumps in the Galaxy (>1000 Msun) and includes representative samples of all embedded stages of high-mass star formation. Here we present the first detailed census of the properties (velocities, distances, luminosities and masses) and spatial distribution of a complete sample of ~8000 dense clumps located in the Galactic disk. We derive highly reliable velocities and distances to ~97% of the sample and use mid- and far-infrared survey data to develop an evolutionary classification scheme that we apply to the whole sample. Comparing the evolutionary subsamples reveals trends for increasing dust temperatures, luminosities and line-widths as a function of evolution indicating that the feedback from the embedded proto-clusters is having a significant impact on the structure and dynamics of their natal clumps. We find 88\,per\,cent are already associated with star formation at some level. We also find the clump mass to be independent of evolution suggesting that the clumps form with the majority of their mass in-situ. We estimate the statistical lifetime of the quiescent stage to be ~5 x 10^4 yr for clump masses ~1000 Msun decreasing to ~1 x 10^4 yr for clump masses >10000 Msun. We find a strong correlation between the fraction of clumps associated with massive stars and peak column density. The fraction is initially small at low column densities but reaching 100\,per\,cent for column densities above 10^{23} cm^{-2}; there are no clumps with column density clumps above this value that are not already associated with massive star formation. All of the evidence is consistent with a dynamic view of star formation wherein the clumps form rapidly and are initially very unstable so that star formation quickly ensues.

(Abridged) Aims: We aim to use the progressive heating of the gas caused by the feedback of high-mass young stellar objects (YSOs) to prove the statistical validity of the most common schemes used to define an evolutionary sequence for high-mass clumps, and characterise the sensitivity of different tracers to this process. Methods: From the spectroscopic follow-ups of the ATLASGAL TOP100 sample, we selected several multiplets of CH3CN, CH3CCH, and CH3OH emission lines to derive and compare the physical properties of the gas in the clumps along the evolutionary sequence. Our findings are compared with results obtained from CO isotopologues, dust, and NH3 from previous studies on the same sample. Results: The chemical properties of each species have a major role on the measured physical properties. Low temperatures are traced by NH3, CH3OH, and CO (in the early phases), the warm and dense envelope can be probed with CH3CN, CH3CCH, and, in evolved sources via CO isotopologues. CH3OH and CH3CN are also abundant in the hot cores, and their high-excitation transitions may be good tools to study the kinematics in the hot gas surrounding the YSOs that these clumps are hosting. All tracers show, to different degrees, progressive warming with evolution. The relation between gas temperature and L/M is reproduced by a toy model of a spherical, internally heated clump. Conclusions: The evolutionary sequence defined for the clumps is statistically valid and we could identify the processes dominating in different intervals of L/M. For L/M<2Lsun/Msun a large quantity of gas is still being accumulated and compressed at the bottom of the potential well. Between 2Lsun/Msun<L/M<40Lsun/Msun the YSOs gain mass and increase in L; the first hot cores appear around L/M=10Lsun/Msun. Finally, for L/M>40Lsun/Msun HII regions become common, showing that dissipation of the parental clump dominates.

High-mass stars are expected to form from dense prestellar cores. Their precise formation conditions are widely discussed, including their virial condition, which results in slow collapse for super-virial cores with strong support by turbulence or magnetic fields, or fast collapse for sub-virial sources. To disentangle their formation processes, measurements of the deuterium fractions are frequently employed to approximately estimate the ages of these cores and to obtain constraints on their dynamical evolution. We here present 3D magneto-hydrodynamical simulations including for the first time an accurate non-equilibrium chemical network with 21 gas-phase species plus dust grains and 213 reactions. With this network we model the deuteration process in fully depleted prestellar cores in great detail and determine its response to variations in the initial conditions. We explore the dependence on the initial gas column density, the turbulent Mach number, the mass-to-magnetic flux ratio and the distribution of the magnetic field, as well as the initial ortho-to-para ratio of H2. We find excellent agreement with recent observations of deuterium fractions in quiescent sources. Our results show that deuteration is rather efficient, even when assuming a conservative ortho-to-para ratio of 3 and highly sub-virial initial conditions, leading to large deuterium fractions already within roughly a free-fall time. We discuss the implications of our results and give an outlook to relevant future investigations.

The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".

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 from their infrared and radio properties to be representative for the whole range of evolutionary stages. 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. We determine flux densities from mid-infrared to submm wavelength (8-870 micron) images and use these values to fit their spectral energy distributions (SEDs) 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. 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. We show that most of the sample has the ability to form massive stars (including the most massive O-type stars) and that the majority is gravitationally unstable and hence likely to be collapsing. The highest column density ATLASGAL sources presented cover the whole range of evolutionary stages from the youngest to the most evolved high-mass star forming clumps. Their study 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.

The mechanisms leading to the formation of disks around young stellar objects (YSOs) and to the launching of the associated jets are crucial to the understanding of the earliest stages of star and planet formation. HH 212 is a privileged laboratory to study a pristine jet-disk system. Therefore we investigate the innermost region ($<100$ AU) around the HH 212-MM1 protostar through ALMA band\,7 observations of methanol. The 8 GHz bandwidth spectrum towards the peak of the continuum emission of the HH 212 system reveals at least 19 transitions of methanol. Several of these lines (among which several vibrationally excited lines in the v$_{\rm t}=1,2$ states) have upper energies above 500 K. They originate from a compact ($<135$ AU in diameter), hot ($\sim 295$ K) region elongated along the direction of the SiO jet. We performed a fit in the $uv$ plane of various velocity channels of the strongest high-excitation lines. The blue- and red-shifted velocity centroids are shifted roughly symmetrically on either side of the jet axis, indicating that the line-of-sight velocity beyond 0.7 km s$^{-1}$ from systemic is dominated by rotational motions. The velocity increases moving away from the protostar further indicating that the emission of methanol is not associated with a Keplerian disk or rotating-infalling cavity, and it is more likely associated with outflowing gas. We speculate that CH$_3$OH traces a disk wind gas accelerated at the base. The launching region would be at a radius of a few astronomical units from the YSO.

Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and the way these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations show that magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ~0.25'' resolution of the thermal dust continuum emission at ~278 GHz towards a turbulent, dense, and massive clump, IRAS 16061-5048c1, in a very early evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure. We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominated by a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in the fragmentation process of massive turbulent clump.

Context: Infrared dark clouds are the coldest and densest portions of giant molecular clouds. The most massive ones represent some of the most likely birthplaces for the next generation of massive stars in the Milky Way. Because a strong mid-IR background is needed to make them appear in absorption, they are usually assumed to be nearby. Aims: We use THz absorption spectroscopy to solve the distance ambiguity associated with kinematic distances for the IR-dark clouds in the TOP100 ATLASGAL sample, a flux-limited selection of massive clumps in different evolutionary phases of star formation. Methods: The para-H2O ground state transition at 1113.343 GHz, observed with Herschel/HIFI, was used to investigate the occurrence of foreground absorption along the line of sight directly towards infrared-dark clouds. Additional consistency checks were performed using MALT90 and HiGAL archival data and targeted Mopra and APEX spectroscopic observations. Results: We report the first discovery of five IRDCs in the TOP100 lying conclusively at the far kinematic distance, showing that the mere presence of low-contrast mid-IR absorption is not sufficient to unequivocally resolve the near/far ambiguity in favour of the former. All IRDCs are massive and actively forming high-mass stars; four of them also show infall signatures. Conclusions: We give a first estimate of the fraction of dark sources at the far distance (~11% in the TOP100) and describe their appearance and properties. The assumption that all dark clouds lie at the near distance may lead, in some cases, to underestimating masses, sizes, and luminosities, possibly causing clouds to be missed that will form very massive stars and clusters.

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.

NGC6357 is an active star forming region with very young massive open clusters (OC). These clusters contain some of the most massive stars in the Galaxy and strongly interact with nearby giant molecular clouds (GMC). We study the young stellar populations of the region and of the OC Pismis24, focusing on their relationship with the nearby GMCs. We seek evidence of triggered star formation propagating from the clusters. We used new deep JHKs photometry, along with unpublished deep IRAC/Spitzer MIR photometry, complemented with optical HST/WFPC2 high spatial resolution photometry and X-ray Chandra observations, to constrain age, initial mass function, and star formation modes in progress. We carefully examine and discuss all sources of bias (saturation, confusion, different sensitivities, extinction). NGC6357 hosts three large young stellar clusters, of which Pismis24 is the most prominent. We found that Pismis24 is a very young (~1-3 Myr) OC with a Salpeter-like IMF and a few thousand members. A comparison between optical and IR photometry indicates that the fraction of members with a NIR excess (i. e., with a circumstellar disk) is in the range 0.3-0.6, consistent with its photometrically derived age. We also find that Pismis24 is likely subdivided into a few different sub-clusters, one of which contains almost all the massive members. There are indications of current star formation triggered by these massive stars, but clear age trends could not be derived (although the fraction of stars with a NIR excess does increase towards the HII region associated with the cluster). The gas out of which Pismis24 formed must have been distributed in dense clumps within a cloud of less dense gas ~1 pc in radius. Our findings provide some new insight into how young stellar populations and massive stars emerge, and evolve in the first few Myr after birth, from a giant molecular cloud complex.

In the low-mass regime, it is found that the gas-phase abundances of C-bearing molecules in cold starless cores rapidly decrease with increasing density, as the molecules form mantles on dust grains. We study CO depletion in 102 massive clumps selected from the ATLASGAL 870 micron survey, and investigate its correlation with evolutionary stage and with the physical parameters of the sources. Moreover, we study the gradients in [12C]/[13C] and [18O]/[17O] isotopic ratios across the inner Galaxy, and the virial stability of the clumps. We use low-J emission lines of CO isotopologues and the dust continuum emission to infer the depletion factor fD. RATRAN one-dimensional models were also used to determine fD and to investigate the presence of depletion above a density threshold. The isotopic ratios and optical depth were derived with a Bayesian approach. We find a significant number of clumps with a large fD, up to ~20. Larger values are found for colder clumps, thus for earlier evolutionary phases. For massive clumps in the earliest stages of evolution we estimate the radius of the region where CO depletion is important to be a few tenths of a pc. Clumps are found with total masses derived from dust continuum emission up to ~20 times higher than the virial mass, especially among the less evolved sources. These large values may in part be explained by the presence of depletion: if the CO emission comes mainly from the low-density outer layers, the molecules may be subthermally excited, leading to an overestimate of the dust masses. CO depletion in high-mass clumps seems to behave as in the low-mass regime, with less evolved clumps showing larger values for the depletion than their more evolved counterparts, and increasing for denser sources. The C and O isotopic ratios are consistent with previous determinations, and show a large intrinsic scatter.

(Abridged) Aims. To investigate the first stages of the process of high-mass star formation, we selected a sample of massive clumps previously observed with the SEST at 1.2 mm and with the ATNF ATCA at 1.3 cm. We want to characterize the physical conditions in such sources, and test whether their properties depend on the evolutionary stage of the clump. Methods. With ATCA we observed the selected sources in the NH3(1,1) and (2,2) transitions and in the 22 GHz H2O maser line. Ammonia lines are a good temperature probe that allow us to accurately determine the mass and the column-, volume-, and surface densities of the clumps. We also collected all data available to construct the spectral energy distribution of the individual clumps and to determine if star formation is already occurring, through observations of its most common signposts, thus putting constraints on the evolutionary stage of the source. We fitted the spectral energy distribution between 1.2 mm and 70 microns with a modified black body to derive the dust temperature and independently determine the mass. Results. The clumps are cold (T~10-30 K), massive (M~10^2-10^3 Mo), and dense (n(H2)>~10^5 cm^-3) and they have high column densities (N(H2)~10^23 cm^-2). All clumps appear to be potentially able to form high-mass stars. The most massive clumps appear to be gravitationally unstable, if the only sources of support against collapse are turbulence and thermal pressure, which possibly indicates that the magnetic field is important in stabilizing them. Conclusions. After investigating how the average properties depend on the evolutionary phase of the source, we find that the temperature and central density progressively increase with time. Sources likely hosting a ZAMS star show a steeper radial dependence of the volume density and tend to be more compact than starless clumps.

Infrared-dark high-mass clumps are among the most promising objects to study the initial conditions of the formation process of high-mass stars and rich stellar clusters. In this work, we have observed the (3-2) rotational transition of C18O with the APEX telescope, and the (1,1) and (2,2) inversion transitions of NH3 with the Australia Telescope Compact Array in 21 infrared-dark clouds already mapped in the 1.2 mm continuum, with the aim of measuring basic chemical and physical parameters such as the CO depletion factor (fD), the gas kinetic temperature and the gas mass. In particular, the C18O (3-2) line allows us to derive fD in gas at densities higher than that traced by the (1-0) and (2-1) lines, typically used in previous works. We have detected NH3 and C18O in all targets. The clumps possess mass, H2 column and surface densities consistent with being potentially the birthplace of high-mass stars. We have measured fD in between 5 and 78, with a mean value of 32 and a median of 29. These values are, to our knowledge, larger than the typical CO depletion factors measured towards infrared-dark clouds and high-mass dense cores, and are comparable to or larger than the values measured in low-mass pre-stellar cores close to the onset of the gravitational collapse. This result suggests that the earliest phases of the high-mass star and stellar cluster formation process are characterised by fD larger than in low-mass pre-stellar cores. Thirteen out of 21 clumps are undetected in the 24 {\mu}m Spitzer images, and have slightly lower kinetic temperatures, masses and H2 column densities with respect to the eight Spitzer-bright sources. This could indicate that the Spitzer-dark clumps are either less evolved or are going to form less massive objects.

The Galactic HII region G353.2+0.9 is excited by the massive open cluster Pismis-24. By analyzing (sub-)mm molecular-line and -continuum we study the detailed morphology of the gas and dust, as well as their physical parameters and their variation across the PDR. We observed various molecules and transitions to derive the physical properties of the molecular gas through line ratios, and both LTE and non-LTE analyses. The physical properties of the gas were derived with a Bayesian approach for the non-LTE analysis. Based on the continuum data at 870 micron, we derived the dust mass and the column density of H2, and thus the molecular abundances. The total mass of the gas in the region is ca. 2000 Mo, while that of the dust is ca. 21 Mo. A velocity gradient in the region suggests that the expansion of the ionized gas is pushing the molecular gas away from the observer. We unambiguously identify the ionization front, at the location of which we detect an increase in gas density and temperature. We find at least 14 clumps at different positions and LSR velocities. We derive kinetic temperatures in the ranges 11-45 K (CS) and 20-45 K (CN). The H2 number density is typically around 1e5 cm^-3 from CS and few 1e5 cm^-3 from CN, with maxima above 1e6 cm^-3. The abundances of the molecules observed vary across the region, and appear to be higher in regions further away from the ionization front.