The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

We have used the ratio of column densities derived independently from the 850-um continuum JCMT Plane Survey (JPS) and the $^{13}$CO/C$^{18}$O J=3-2 Heterodyne Inner Milky Way Plane Survey (CHIMPS) to produce maps of the dense-gas mass fraction (DGMF) in two slices of the Galactic Plane centred at l=30 and l=40. The observed DGMF is a metric for the instantaneous clump-formation efficiency (CFE) in the molecular gas and a two-dimensional power-spectrum analysis of the DGMF reveals a break in slope at the approximate size scale of molecular clouds. We interpret this as the characteristic scale of the amplitude of variations in the CFE and a constraint on the dominant mechanism regulating the CFE and, hence, the star-formation efficiency in CO-traced clouds. By splitting the two fields into velocity components corresponding to the spiral arms that cross them, we find that individual molecular clouds have mean CFE values of 40, 41, and 46 per cent for the Scutum--Centaurus, Sagittarius, and Perseus spiral arms, respectively

Recent surveys of the Galactic plane in the dust continuum and CO emission lines reveal that large ($\gtrsim 50$~pc) and massive ($\gtrsim 10^5$~$M_\odot$) filaments, know as giant molecular filaments (GMFs), may be linked to galactic dynamics and trace the mid-plane of the gravitational potential in the Milky Way. We have imaged one entire GMF located at $l\sim$52--54$^\circ$ longitude, GMF54 ($\sim$68~pc long), in the empirical dense gas tracers using the HCN(1--0), HNC(1--0), HCO$^+$(1--0) lines, and their $^{13}$C isotopologue transitions, as well as the N$_2$H$^+$(1--0) line. We study the dense gas distribution, the column density probability density functions (N-PDFs) and the line ratios within the GMF. The dense gas molecular transitions follow the extended structure of the filament with area filling factors between 0.06 and 0.28 with respect to $^{13}$CO(1--0). We constructed the N-PDFs of H$_2$ for each of the dense gas tracers based on their column densities and assumed uniform abundance. The N-PDFs of the dense gas tracers appear curved in log-log representation, and the HCO$^+$ N-PDF has the largest log-normal width and flattest power-law slope index. Studying the N-PDFs for sub-regions of GMF54, we found an evolutionary trend in the N-PDFs that high-mass star forming and Photon-Dominate Regions (PDRs) have flatter power-law indices. The integrated intensity ratios of the molecular lines in GMF54 are comparable to those in nearby galaxies. In particular, the N$_2$H$^+$/$^{13}$CO ratio, which traces the dense gas fraction, has similar values in GMF54 and all nearby galaxies except ULIRGs.

Context: Atomic and molecular cloud formation is a dynamical process. However, kinematic signatures of these processes are still observationally poorly constrained. Methods: Targeting the cloud-scale environment of the prototypical infrared dark cloud G28.3, we employ spectral line imaging observations of the two atomic lines HI and [CI] as well as molecular lines observations in 13CO in the 1--0 and 3--2 transitions. The analysis comprises investigations of the kinematic properties of the different tracers, estimates of the mass flow rates, velocity structure functions, a Histogram of Oriented Gradients (HOG) study as well as comparisons to simulations. Results: The central IRDC is embedded in a more diffuse envelope of cold neutral medium (CNM) traced by HI self-absorption (HISA) and molecular gas. The spectral line data as well as the HOG and structure function analysis indicate a possible kinematic decoupling of the HI from the other gas compounds. Spectral analysis and position-velocity diagrams reveal two velocity components that converge at the position of the IRDC. Estimated mass flow rates appear rather constant from the cloud edge toward the center. The velocity structure function analysis is consistent with gas flows being dominated by the formation of hierarchical structures. Conclusions: The observations and analysis are consistent with a picture where the IRDC G28 is formed at the center of two converging gas flows. While the approximately constant mass flow rates are consistent with a self-similar, gravitationally driven collapse of the cloud, external compression by, e.g., spiral arm shocks or supernovae explosions cannot be excluded yet. Future investigations should aim at differentiating the origin of such converging gas flows.

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

We present a CO(3-2) survey of selected regions in the M31 disc as part of the JCMT large programme, HARP and SCUBA-2 High-Resolution Terahertz Andromeda Galaxy Survey (HASHTAG). The 12 CO(3-2) fields in this survey cover a total area of 60 square arcminutes, spanning a deprojected radial range of 2 - 14 kpc across the M31 disc. Combining these observations with existing IRAM 30m CO(1-0) observations and JCMT CO(3-2) maps of the nuclear region of M31, as well as dust temperature and star formation rate surface density maps, we are able to explore the radial distribution of the CO(3-2)/CO(1-0) integrated intensity ratio (R31) and its relationship with dust temperature and star formation. We find that the value of R31 between 2 - 9 kpc galactocentric radius is 0.14, significantly lower than what is seen in the nuclear ring at ~1 kpc (R31 ~ 0.8), only to rise again to 0.27 for the fields centred on the 10 kpc star forming ring. We also found that R31 is positively correlated with dust temperature, with Spearman's rank correlation coefficient $\rho$ = 0.55. The correlation between star formation rate surface density and CO(3--2) intensity is much stronger than with CO(1-0), with $\rho$ = 0.54 compared to -0.05, suggesting that the CO(3-2) line traces warmer and denser star forming gas better. We also find that R31 correlates well with star formation rate surface density, with $\rho$ = 0.69.

The latest generation of high-angular-resolution unbiased Galactic plane surveys in molecular-gas tracers are enabling the interiors of molecular clouds to be studied across a range of environments. The CHIMPS survey simultaneously mapped a sector of the inner Galactic plane, within 27.8 < l < 46.2 deg and |b| < 0.5 deg, in 13CO and C18O (3-2) at 15 arcsec resolution. The combination of CHIMPS data with 12CO (3-2) data from the COHRS survey has enabled us to perform a voxel-by-voxel local-thermodynamic-equilibrium analysis, determining the excitation temperature, optical depth, and column density of 13CO at each l,b,v position. Distances to discrete sources identified by FellWalker in the 13CO (3-2) emission maps were determined, allowing the calculation of numerous physical properties of the sources, and we present the first source catalogues in this paper. We find that, in terms of size and density, the CHIMPS sources represent an intermediate population between large-scale molecular clouds identified by CO and dense clumps seen in dust emission, and therefore represent the bulk transition from the diffuse to the dense phase of molecular gas. We do not find any significant systematic variations in the masses, column densities, virial parameters, excitation temperature, or the turbulent pressure over the range of Galactocentric distance probed, but we do find a shallow increase in the mean volume density with increasing Galactocentric distance. We find that inter-arm clumps have significantly narrower linewidths, and lower virial parameters and excitation temperatures than clumps located in spiral arms. When considering the most reliable distance-limited subsamples, the largest variations occur on the clump-to-clump scale, echoing similar recent studies that suggest that the star-forming process is largely insensitive to the Galactic-scale environment, at least within the inner disc.

Identifying the structure of our Galaxy has always been fraught with difficulties, and while modern surveys continue to make progress building a map of the Milky Way, there is still much to understand. The arm and bar features are important drivers in shaping the interstellar medium, but their exact nature and influence still require attention. We present results of smoothed particle hydrodynamic simulations of gas in the Milky Way including star formation, stellar feedback, and ISM cooling, when exposed to different arm and bar features, with the aim of better understanding how well newly formed stars trace out the underlying structure of the Galaxy. The bar is given a faster pattern speed than the arms, resulting in a complex, time-dependent morphology and star formation. Inter-arm branches and spurs are easily influenced by the bar, especially in the two-armed spiral models where there is a wide region of resonance overlap in the disc. As the bar over-takes the spiral arms it induces small boosts in star formation and enhances spiral features, which occur at regularly spaced beat-like intervals. The locations of star formation events are similar to those seen in observational data, and do not show a perfect 1:1 correspondence with the underlying spiral potential, though arm tangencies are generally well traced by young stars. Stellar velocity fields from the newly formed stars are compared to data from Gaia DR2, showing that the spiral and bar features can reproduce many of the non-axisymmetric features seen in the data. A simple analytical model is used to show many of these feature are a natural response of gas to rigidly rotating spiral and bar potentials.

We present images of radio recombination lines (RRLs) at wavelengths around 17 cm from the star-forming region W49A to determine the kinematics of ionized gas in the THOR survey (The HI/OH/Recombination line survey of the inner Milky Way) at an angular resolution of 16.8"x13.8". The distribution of ionized gas appears to be affected by feedback processes from the star clusters in W49A. The velocity structure of the RRLs shows a complex behavior with respect to the molecular gas. We find a shell-like distribution of ionized gas as traced by RRL emission surrounding the central cluster of OB stars in W49A. We describe the evolution of the shell with the recent feedback model code WARPFIELD that includes the important physical processes and has previously been applied to the 30 Doradus region in the Large Magellanic Cloud. The cloud structure and dynamics of W49A are in agreement with a feedback-driven shell that is re-collapsing. The shell may have triggered star formation in other parts of W49A. We suggest that W49A is a potential candidate for star formation regulated by feedback-driven and re-collapsing shells.

Our understanding of how molecular clouds form in the interstellar medium (ISM) would be greatly helped if we had a reliable observational tracer of the gas flows responsible for forming the clouds. Fine structure emission from singly ionised and neutral carbon ([CII], [CI]) and rotational line emission from CO are all observed to be associated with molecular clouds. However, it remains unclear whether any of these tracers can be used to study the inflow of gas onto an assembling cloud, or whether they primarily trace the cloud once it has already assembled. In this paper, we address this issue with the help of high resolution simulations of molecular cloud formation that include a sophisticated treatment of the chemistry and thermal physics of the ISM. Our simulations demonstrate that both [CI] and CO emission trace gas that is predominantly molecular, with a density n ~ 500-1000 cm$^{-3}$, much larger than the density of the inflowing gas. [CII] traces lower density material (n ~ 100 cm$^{-3}$) that is mainly atomic at early times. A large fraction of the [CII] emission traces the same structures as the [CI] or CO emission, but some arises in the inflowing gas. Unfortunately, this emission is very faint and will be difficult to detect with current observational facilities, even for clouds situated in regions with an elevated interstellar radiation field.

In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO $J=4(0,4)-3(0,3)$ data at 3 mm, using two new pieces of software which we make available to the community. First, scousepy, a Python implementation of the spectral line fitting algorithm scouse. Secondly, acorns (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line of sight velocity dispersion in this cloud, $\sigma_{v_{los}, {\rm 1D}}=4.4\pm2.1$ kms$^{-1}$, which is somewhat larger than predicted by velocity dispersion-size relations for the Central Molecular Zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding $\sigma_{v_{los}, {\rm 1D}}/\sigma_{v_{pos}, {\rm 1D}}\sim1.2\pm0.3$. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling we conclude that G0.253+0.016 is not a single, coherent, and centrally-condensed molecular cloud; 'the Brick' is not a \emph{brick}. Instead, G0.253+0.016 is a dynamically complex and hierarchically-structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ.

We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall.

The hydroxyl radical (OH) is present in the diffuse molecular and partially atomic phases of the interstellar medium (ISM), but its abundance relative to hydrogen is not clear. We aim to evaluate the abundance of OH with respect to molecular hydrogen using OH absorption against cm-continuum sources over the first Galactic quadrant. This OH study is part of the HI/OH/Recombination line survey (THOR). THOR is a Karl G. Jansky Very Large Array large program of atomic, molecular and ionized gas in the range 15{\deg}$\leq$l$\leq$67{\deg} and |b|$\leq$1{\deg}. It is the highest-resolution unbiased OH absorption survey to date towards this region. We combine the derived optical depths with literature 13CO(1-0) and HI observations to determine the OH abundance. We detect absorption in the 1665 and 1667 MHz transitions for continuum sources stronger than $F_{\rm cont}\geq$0.1 Jy/beam. OH absorption is found against $\sim$15% of these continuum sources with increasing fractions for stronger sources. Most of the absorption is associated with Galactic HII regions. We find OH and 13CO gas to have similar kinematic properties. The OH abundance decreases with increasing hydrogen column density. The OH abundance with respect to the total hydrogen nuclei column density (atomic and molecular phase) is in agreement with a constant abundance for $A_V$ < 10-20. Towards the lowest column densities, we find sources that exhibit OH absorption but no 13CO emission, indicating that OH is a well suited tracer of the low column density molecular gas. We present spatially resolved OH absorption towards W43. The unbiased nature of the THOR survey opens a new window onto the gas properties of the ISM. The characterization of the OH abundance over a large range of hydrogen gas column densities contributes to the understanding of OH as a molecular gas tracer and provides a starting point for future investigations.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

What role does Galactic structure play in star formation? We have used the Herschel Hi-GAL compact-clump catalogue to examine trends in evolutionary stage over large spatial scales in the inner Galaxy. We examine the relationship between the fraction of clumps with embedded star formation (the star-forming fraction, or SFF) and other measures of star-formation activity. Based on a positive correlation between SFF and evolutionary indicators such as the luminosity-to-mass ratio, we assert that the SFF principally traces the average evolutionary state of a sample and must depend on the local fraction of rapidly-evolving, high-mass young stellar objects. The spiral-arm tangent point longitudes show small excesses in the SFF, though these can be accounted for by a small number of the most massive clusters, just 7.6% of the total number of clumps in the catalogue. This suggests that while the arms tend to be home to the Galaxy's massive clusters, the remaining 92.4% of Hi-GAL clumps in our catalogue do not show an enhancement of star formation within arms. Globally, the SFF is highest at the Galactic midplane and inner longitudes. We find no significant trend in evolutionary stage as a function of position across spiral arms at the tangent-point longitudes. This indicates that the angular offset observed between gas and stars, if coordinated by a density wave, is not evident at the clump phase; alternatively, the onset of star formation is not triggered by the spiral density wave.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

Infrared dark clouds (IRDCs) are thought to be potential hosts of the elusive early phases of high-mass star formation. Here we conduct an in-depth kinematic analysis of one such IRDC, G034.43+00.24 (Cloud F), using high sensitivity and high spectral resolution IRAM-30m N$_2$H$^+$ ($1-0$) and C$^{18}$O ($1-0$) observations. To disentangle the complex velocity structure within this cloud we use Gaussian decomposition and hierarchical clustering algorithms. We find that four distinct coherent velocity components are present within Cloud F. The properties of these components are compared to those found in a similar IRDC, G035.39-00.33 (Cloud H). We find that the components in both clouds have: high densities (inferred by their identification in N$_2$H$^+$), trans-to-supersonic non-thermal velocity dispersions with Mach numbers of $\sim$ $1.5-4$, a separation in velocity of $\sim$3 km s$^{-1}$, and a mean red-shift of $\sim$ 0.3 km s$^{-1}$ between the N$_2$H$^+$ (dense gas) and C$^{18}$O emission (envelope gas). The latter of these could suggest that these clouds share a common formation scenario. We investigate the kinematics of the larger-scale Cloud F structures, using lower-density-tracing $^{13}$CO ($1-0$) observations. A good correspondence is found between the components identified in the IRAM-30m observations and the most prominent component in the$^{13}$CO data. We find that the IRDC Cloud F is only a small part of a much larger structure, which appears to be an inter-arm filament of the Milky Way.

We describe new Hi-GAL based maps of the entire Galactic Plane, obtained using continuum data in the wavelength range 70-500 $\mu$m. These maps are derived with the PPMAP procedure, and therefore represent a significant improvement over those obtained with standard analysis techniques. Specifically they have greatly improved resolution (12 arcsec) and, in addition to more accurate integrated column densities and mean dust temperatures, they give temperature-differential column densities, i.e., separate column density maps in twelve distinct dust temperature intervals, along with the corresponding uncertainty maps. The complete set of maps is available on-line. We briefly describe PPMAP and present some illustrative examples of the results. These include (a) multi-temperature maps of the Galactic HII region W5-E, (b) the temperature decomposition of molecular cloud column-density probability distribution functions, and (c) the global variation of mean dust temperature as a function of Galactocentric distance. Amongst our findings are: (i) a strong localised temperature gradient in W5-E in a direction orthogonal to that towards the ionising star, suggesting an alternative heating source and providing possible guidance for models of the formation of the bubble complex, and (ii) the overall radial profile of dust temperature in the Galaxy shows a monotonic decrease, broadly consistent both with models of the interstellar radiation field and with previous estimates at lower resolution. However, we also find a central temperature plateau within ~ 6 kpc of the Galactic centre, outside of which is a pronounced steepening of the radial profile. This behaviour may reflect the greater proportion of molecular (as opposed to atomic) gas in the central region of the Galaxy.

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions.

The dense, cold regions where high-mass stars form are poorly characterised, yet they represent an ideal opportunity to learn more about the initial conditions of high-mass star formation (HMSF), since high-mass starless cores (HMSCs) lack the violent feedback seen at later evolutionary stages. We present continuum maps obtained from the Submillimeter Array (SMA) interferometry at 1.1 mm for four infrared dark clouds (IRDCs, G28.34S, IRDC 18530, IRDC 18306, and IRDC 18308). We also present 1 mm/3 mm line surveys using IRAM 30 m single-dish observations. Our results are: (1) At a spatial resolution of 10^4 AU, the 1.1 mm SMA observations resolve each source into several fragments. The mass of each fragment is on average >10 Msun, which exceeds the predicted thermal Jeans mass of the whole clump by a factor of up to 30, indicating that thermal pressure does not dominate the fragmentation process. Our measured velocity dispersions in the 30 m lines imply that non-thermal motions provides the extra support against gravity in the fragments. (2) Both non-detection of high-J transitions and the hyperfine multiplet fit of N2H+(1-0), C2H(1-0), HCN(1-0), and H13CN(1-0) indicate that our sources are cold and young. However, obvious detection of SiO and the asymmetric line profile of HCO+(1-0) in G28.34S indicate a potential protostellar object and probable infall motion. (3) With a large number of N-bearing species, the existence of carbon rings and molecular ions, and the anti-correlated spatial distributions between N2H+/NH2D and CO, our large-scale high-mass clumps exhibit similar chemical features as small-scale low-mass prestellar objects. This study of a small sample of IRDCs illustrates that thermal Jeans instability alone cannot explain the fragmentation of the clump into cold (~15 K), dense (>10^5 cm-3) cores and that these IRDCs are not completely quiescent.

We present large-scale trends in the distribution of star-forming objects revealed by the Hi-GAL survey. As a simple metric probing the prevalence of star formation in Hi-GAL sources, we define the fraction of the total number of Hi-GAL sources with a 70-micron counterpart as the "star-forming fraction" or SFF. The mean SFF in the inner galactic disc (3.1 kpc < R_GC < 8.6 kpc) is 25%. Despite an apparent pile-up of source numbers at radii associated with spiral arms, the SFF shows no significant deviations at these radii, indicating that the arms do not affect the star-forming productivity of dense clumps either via physical triggering processes or through the statistical effects of larger source samples associated with the arms. Within this range of Galactocentric radii, we find that the SFF declines with R_GC at a rate of -0.026 +/- 0.002 per kiloparsec, despite the dense gas mass fraction having been observed to be constant in the inner Galaxy. This suggests that the SFF may be weakly dependent on one or more large-scale physical properties of the Galaxy, such as metallicity, radiation field, pressure or shear, such that the dense sub-structures of molecular clouds acquire some internal properties inherited from their environment.

Constraining the temperature and density structure of dense molecular cloud cores is fundamental for understanding the initial conditions of star formation. We use Herschel observations of the thermal FIR dust emission from nearby isolated molecular cloud cores and combine them with ground-based submillimeter continuum data to derive observational constraints on their temperature and density structure. The aim of this study is to verify the validity of a ray-tracing inversion technique developed to derive the dust temperature and density structure of isolated starless cores directly from the dust emission maps and to test if the resulting temperature and density profiles are consistent with physical models. Using this ray-tracing inversion technique, we derive the dust temperature and density structure of six isolated starless cloud cores. We employ self-consistent radiative transfer modeling to the derived density profiles, treating the ISRF as the only heating source. The best-fit values of local strength of the ISRF and the extinction by the outer envelope are derived by comparing the self-consistently calculated temperature profiles with those derived by the ray-tracing method. We find that all starless cores are significantly colder inside than outside, with the core temperatures showing a strong negative correlation with peak column density. This suggests that their thermal structure is dominated by external heating from the ISRF and shielding by dusty envelopes. The temperature profiles derived with the ray-tracing inversion method can be well-reproduced with self-consistent radiative transfer models.

Image cubes of differential column density as a function of dust temperature are constructed for Galactic Centre molecular cloud G0.253+0.016 ("The Brick") using the recently described PPMAP procedure. The input data consist of continuum images from the Herschel Space Telescope in the wavelength range 70-500 $\mu$m, supplemented by previously published interferometric data at 1.3 mm wavelength. While the bulk of the dust in the molecular cloud is consistent with being heated externally by the local interstellar radiation field, our image cube shows the presence, near one edge of the cloud, of a filamentary structure whose temperature profile suggests internal heating. The structure appears as a cool ($\sim 14$ K) tadpole-like feature, $\sim 6$ pc in length, in which is embedded a thin spine of much hotter ($\sim$ 40-50 K) material. We interpret these findings in terms of a cool filament whose hot central region is undergoing gravitational collapse and fragmentation to form a line of protostars. If confirmed, this would represent the first evidence of widespread star formation having started within this cloud.

Filamentary structures are common morphological features of the cold, molecular interstellar medium (ISM). Recent studies have discovered massive, hundred-parsec-scale filaments that may be connected to the large-scale, Galactic spiral arm structure. Addressing the nature of these Giant Molecular Filaments (GMFs) requires a census of their occurrence and properties. We perform a systematic search of GMFs in the fourth Galactic quadrant and determine their basic physical properties. We identify GMFs based on their dust extinction signatures in near- and mid-infrared and velocity structure probed by ^{13}CO line emission. We use the ^{13}CO line emission and ATLASGAL dust emission data to estimate the total and dense gas masses of the GMFs. We combine our sample with an earlier sample from literature and study the Galactic environment of the GMFs. We identify nine GMFs in the fourth Galactic quadrant; six are located in the Centaurus spiral arm and three in inter-arm regions. Combining this sample with an earlier study using the same identification criteria in the first Galactic quadrant results in 16 GMFs, nine of which are located within spiral arms. The GMFs have sizes of 80-160 pc and ^{13}CO-derived masses between 5-90 x 10^{4} Msun. Their dense gas mass fractions are between 1.5-37%, being higher in the GMFs connected to spiral arms. We also compare the different GMF-identification methods and find that emission and extinction based techniques overlap only partially, highlighting the need to use both to achieve a complete census.

We run numerical simulations of molecular clouds (MCs), adopting properties similar to those found in the Central Molecular Zone (CMZ) of the Milky Way. For this, we employ the moving mesh code Arepo and perform simulations which account for a simplified treatment of time-dependent chemistry and the non-isothermal nature of gas and dust. We perform simulations using an initial density of n_0 = 10^3 cm^{-3} and a mass of 1.3x10^5 M_sun. Furthermore, we vary the virial parameter, defined as the ratio of kinetic and potential energy, alpha = E_{kin} / |E_{pot}|, by adjusting the velocity dispersion. We set it to alpha = 0.5, 2.0 and 8.0, in order to analyze the impact of the kinetic energy on our results. We account for the extreme conditions in the CMZ and increase both the interstellar radiation field (ISRF) and the cosmic-ray flux (CRF) by a factor of 1000 compared to the values found in the solar neighbourhood. We use the radiative transfer code RADMC-3D to compute synthetic images in various diagnostic lines. These are [CII] at 158 micron, [OI] (145 micron), [OI] (63 micron), 12CO (J = 1 -> 0) and 13CO (J = 1 -> 0) at 2600 micron and 2720 micron, respectively. When alpha is large, the turbulence disperses much of the gas in the cloud, reducing its mean density and allowing the ISRF to penetrate more deeply into the cloud's interior. This significantly alters the chemical composition of the cloud, leading to the dissociation of a significant amount of the molecular gas. On the other hand, when alpha is small, the cloud remains compact, allowing more of the molecular gas to survive. We show that in each case the atomic tracers accurately reflect most of the physical properties of both the H2 and the total gas of the cloud and that they provide a useful alternative to molecular lines when studying the ISM in the CMZ.

We report on observations of the hydroxyl radical (OH) within The H{\sc I}, OH Recombination line survey (THOR) pilot region. The region is bounded approximately between Galactic coordinates l=29.2 to 31.5$^\circ$ and b=-1.0 to +1.0$^\circ$ and includes the high-mass star forming region W43. We identify 103 maser sites, including 72 with 1612\,MHz masers, 42 showing masers in either of the main line transitions at 1665 and 1667\,MHz and four showing 1720\,MHz masers. Most maser sites with either main-line or 1720\,MHz emission are associated with star formation, whereas most of the 1612\,MHz masers are associated with evolved stars. We find that nearly all of the main-line maser sites are co-spatial with an infrared source, detected by GLIMPSE. We also find diffuse OH emission, as well as OH in absorption towards selected unresolved or partially resolved sites. Extended OH absorption is found towards the well known star forming complex W43 Main.

Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU). Results: The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(+-0.18)pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, gravitationally bound gas cylinder, a high mass-to-length ratio of M/l~1000M_sun/pc requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1-0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent sub-filaments. Furthermore, this velocity gradient perpendicular to the filament resembles recent results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers ([CI] and C18O) reveal a similar red/blueshifted velocity structure on scales around 60'' east and west of the IRDC18223 filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller-scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may either be due to no significant gas flows along the filamentary axis, but it may partly also be caused by a low inclination angle of the filament with respect to the plane of the sky that could minimize such signature.

Aims: Understanding the fragmentation and collapse properties of the dense gas during the onset of high-mass star formation. Methods: We observed the massive (~800M_sun) starless gas clump IRDC18310-4 with the Plateau de Bure Interferometer (PdBI) at sub-arcsecond resolution in the 1.07mm continuum andN2H+(3-2) line emission. Results: Zooming from a single-dish low-resolution map to previous 3mm PdBI data, and now the new 1.07mm continuum observations, the sub-structures hierarchically fragment on the increasingly smaller spatial scales. While the fragment separations may still be roughly consistent with pure thermal Jeans fragmentation, the derived core masses are almost two orders of magnitude larger than the typical Jeans mass at the given densities and temperatures. However, the data can be reconciled with models using non-homogeneous initial density structures, turbulence and/or magnetic fields. While most sub-cores remain (far-)infrared dark even at 70mum, we identify weak 70mum emission toward one core with a comparably low luminosity of ~16L_sun, re-enforcing the general youth of the region. The spectral line data always exhibit multiple spectral components toward each core with comparably small line widths for the individual components (in the 0.3 to 1.0km/s regime). Based on single-dish C18O(2-1) data we estimate a low virial-to-gas-mass ratio <=0.25. We discuss that the likely origin of these spectral properties may be the global collapse of the original gas clump that results in multiple spectral components along each line of sight. Even within this dynamic picture the individual collapsing gas cores appear to have very low levels of internal turbulence.

To study the atomic, molecular and ionized emission of Giant Molecular Clouds (GMCs), we have initiated a Large Program with the VLA: 'THOR - The HI, OH, Recombination Line survey of the Milky Way'. We map the 21cm HI line, 4 OH lines, 19 H_alpha recombination lines and the continuum from 1 to 2 GHz of a significant fraction of the Milky Way (l=15-67deg, |b|<1deg) at ~20" resolution. In this paper, we focus on the HI emission from the W43 star-formation complex. Classically, the HI 21cm line is treated as optically thin with column densities calculated under this assumption. This might give reasonable results for regions of low-mass star-formation, however, it is not sufficient to describe GMCs. We analyzed strong continuum sources to measure the optical depth, and thus correct the HI 21cm emission for optical depth effects and weak diffuse continuum emission. Hence, we are able to measure the HI mass of W43 more accurately and our analysis reveals a lower limit of M=6.6x10^6 M_sun, which is a factor of 2.4 larger than the mass estimated with the assumption of optically thin emission. The HI column densities are as high as N(HI)~150 M_sun/pc^2 ~ 1.9x10^22 cm^-2, which is an order of magnitude higher than for low mass star formation regions. This result challenges theoretical models that predict a threshold for the HI column density of ~10 M_sun/pc^2, at which the formation of molecular hydrogen should set in. By assuming an elliptical layered structure for W43, we estimate the particle density profiles. While at the cloud edge atomic and molecular hydrogen are well mixed, the center of the cloud is strongly dominated by molecular hydrogen. We do not identify a sharp transition between hydrogen in atomic and molecular form. Our results are an important characterization of the atomic to molecular hydrogen transition in an extreme environment and challenges current theoretical models.

We want to understand the kinematic and thermal properties of young massive gas clumps prior to and at the earliest evolutionary stages of high-mass star formation. Do we find signatures of gravitational collapse? Do we find temperature gradients in the vicinity or absence of infrared emission sources? Do we find coherent velocity structures toward the center of the dense and cold gas clumps? To determine kinematics and gas temperatures, we used ammonia, because it is known to be a good tracer and thermometer of dense gas. We observed the NH$_3$(1,1) and (2,2) lines within seven very young high-mass star-forming regions with the VLA and the Effelsberg 100m telescope. This allows us to study velocity structures, linewidths, and gas temperatures at high spatial resolution of 3-5$"$, corresponding to $\sim$0.05 pc. We find on average cold gas clumps with temperatures in the range between 10 K and 30 K. The observations do not reveal a clear correlation between infrared emission peaks and ammonia temperature peaks. We report an upper limit for the linewidth of $\sim$1.3 km s$^{-1}$, at the spectral resolution limit of our VLA observation. This indicates a relatively low level of turbulence on the scale of the observations. Velocity gradients are present in almost all regions with typical velocity differences of 1 to 2 km s$^{-1}$ and gradients of 5 to 10 km s$^{-1}$ pc$^{-1}$. These velocity gradients are smooth in most cases, but there is one exceptional source (ISOSS23053), for which we find several velocity components with a steep velocity gradient toward the clump centers that is larger than 30 km s$^{-1}$ pc$^{-1}$. This steep velocity gradient is consistent with recent models of cloud collapse. Furthermore, we report a spatial correlation of ammonia and cold dust, but we also find decreasing ammonia emission close to infrared emission sources.

The massive infrared dark cloud G0.253+0.016 projected 45pc from the Galactic centre contains ~10^5Msun of dense gas whilst being mostly devoid of observed star-formation tracers. To scrutinise the physical properties, dynamics and structure of this cloud with reference to its star-forming potential, we have carried out a concerted SMA and IRAM 30m study of this cloud in dust continuum, CO isotopologues, shock tracing molecules, as well as H$_2$CO to trace the gas temperature. We detect and characterise the dust cores within G0.253+0.016 at ~1.3 mm and find that the kinetic temperature of the gas is >320K on size-scales of ~0.15 pc. Analysis of the position-velocity diagrams of our observed lines show broad linewidths and strong shock emission in the south of the cloud, indicating that G0.253+0.016 is colliding with another cloud at v(LSR)~70 km/s. We confirm via an analysis of the observed dynamics in the CMZ that it is an elongated structure, orientated with Sgr B2 closer to the Sun than Sgr A*, however our results suggest that the actual geometry may be more complex than an elliptical ring. We find that the column density PDF of G0.253+0.016 is log-normal with no discernible power-law tail, consistent with little star formation, and that its width can be explained in the framework of theory predicting the density structure of clouds created by supersonic, magnetised turbulence. We also present the delta-variance spectrum of this region, and show it is consistent with that expected for clouds with no star formation. Using G0.253+0.016 as a test-bed of the conditions required for star formation in a different physical environment to that of nearby clouds, we also conclude that there is not one column density threshold for star formation, but instead this value is dependant on the local physical conditions. [Abbrv.]

We present new Plateau de Bure Interferometer observations of a region in the filamentary infrared-dark cloud (IRDC) G011.11-0.12 containing young, star-forming cores. In addition to the 3.2mm continuum emission from cold dust, we map this region in the N$_2$H$^+$(1-0) line to trace the core kinematics with an angular resolution of 2" and velocity resolution of 0.2km s$^{-1}$. These data are presented in concert with recent {\em Herschel} results, single-dish N$_2$H$^+$(1-0) data, SABOCA 350$\mu$m continuum data, and maps of the C$^{18}$O (2-1) transition obtained with the IRAM 30m telescope. We recover the star-forming cores at 3.2mm continuum, while in N$_2$H$^+$ they appear at the peaks of extended structures. The mean projected spacing between N$_2$H$^+$ emission peaks is 0.18pc, consistent with simple isothermal Jeans fragmentation. The 0.1pc-sized cores have low virial parameters on the criticality borderline, while on the scale of the whole region, we infer that it is undergoing large-scale collapse. The N$_2$H$^+$ linewidth increases with evolutionary stage, while CO isotopologues show no linewidth variation with core evolution. Centroid velocities of all tracers are in excellent agreement, except in the starless region where two N$_2$H$^+$ velocity components are detected, one of which has no counterpart in C$^{18}$O. We suggest that gas along this line of sight may be falling into the quiescent core, giving rise to the second velocity component, possibly connected to the global collapse of the region.

We present CARMA 2.9 mm dust continuum emission observations of a sample of 14 Herschel-detected Class 0 protostars in the Orion A and B molecular clouds, drawn from the PACS Bright Red Sources (PBRS) sample (Stutz et al.). These objects are characterized by very red 24 \micron\ to 70 \micron\ colors and prominent submillimeter emission, suggesting that they are very young Class 0 protostars embedded in dense envelopes. We detect all of the PBRS in 2.9 mm continuum emission and emission from 4 protostars and 1 starless core in the fields toward the PBRS; we also report 1 new PBRS source. The ratio of 2.9 mm luminosity to bolometric luminosity is higher by a factor of $\sim$5 on average, compared to other well-studied protostars in the Perseus and Ophiuchus clouds. The 2.9 mm visibility amplitudes for 6 of the 14 PBRS are very flat as a function of uv-distance, with more than 50\% of the source emission arising from radii $<$ 1500 AU. These flat visibility amplitudes are most consistent with spherically symmetric envelope density profiles with $\rho$~$\propto$~R$^{-2.5}$. Alternatively, there could be a massive unresolved structure like a disk or a high-density inner envelope departing from a smooth power-law. The large amount of mass on scales $<$ 1500 AU (implying high average central densities) leads us to suggest that that the PBRS with flat visibility amplitude profiles are the youngest PBRS and may be undergoing a brief phase of high mass infall/accretion and are possibly among the youngest Class 0 protostars. The PBRS with more rapidly declining visibility amplitudes still have large envelope masses, but could be slightly more evolved.

Context: How do molecular clouds form out of the atomic phase? And what are the relative fractions of carbon in the ionized, atomic and molecular phase? These are questions at the heart of cloud and star formation. Methods: Using multiple observatories from Herschel and SOFIA to APEX and the IRAM 30m telescope, we mapped the ionized, atomic and molecular carbon ([CII]@1900GHz, [CI]@492GHz and C18O(2-1)@220GHz) at high spatial resolution (12"-25") in four young massive infrared dark clouds (IRDCs). Results: The three carbon phases were successfully mapped in all four regions, only in one source the [CII] line remained a non-detection. Both the molecular and atomic phases trace the dense structures well, with [CI] also tracing material at lower column densities. [CII] exhibits diverse morphologies in our sample, from compact to diffuse structures probing the cloud environment. In at least two out of the four regions, we find kinematic signatures strongly indicating that the dense gas filaments have formed out of a dynamically active and turbulent atomic/molecular cloud, potentially from converging gas flows. The atomic-to-molecular carbon gas mass ratios are low between 7% and 12% with the lowest values found toward the most quiescent region. In the three regions where [CII] is detected, its mass is always higher by a factor of a few than that of the atomic carbon. The ionized carbon emission depends as well on the radiation field, however, we also find strong [CII] emission in a region without significant external sources, indicating that other processes, e.g., energetic gas flows can contribute to the [CII] excitation as well.

Throughout the Milky Way, molecular clouds typically appear filamentary, and mounting evidence indicates that this morphology plays an important role in star formation. What is not known is to what extent the dense filaments most closely associated with star formation are connected to the surrounding diffuse clouds up to arbitrarily large scales. How are these cradles of star formation linked to the Milky Way's spiral structure? Using archival Galactic plane survey data, we have used multiple datasets in search of large-scale, velocity-coherent filaments in the Galactic plane. In this paper, we present our methods employed to identify coherent filamentary structures first in extinction and confirmed using Galactic Ring Survey data. We present a sample of seven Giant Molecular Filaments (GMFs) that have lengths of order $\sim$100 pc, total masses of 10$^4$ - 10$^5$ M$_{\odot}$, and exhibit velocity coherence over their full length. The GMFs we study appear to be inter-arm clouds and may be the Milky Way analogues to spurs observed in nearby spiral galaxies. We find that between 2 and 12% of the total mass (above $\sim$10$^{20}$ cm$^{-2}$) is "dense" (above 10$^{22}$ cm$^{-2}$), where filaments near spiral arms in the Galactic midplane tend to have higher dense gas mass fractions than those further from the arms.

We present Submillimeter Array (SMA) $\lambda =$ 0.88 and 1.3 mm broad band observations, and the Jansky Very Large Array (VLA) observations in $\rm{NH_3}$ $(J,K) = (1,1)$ up to $(5,5)$, $\rm{H_2O}$ and $\rm{CH_3OH}$ maser lines toward the two most massive molecular clumps in infrared dark cloud (IRDC) G11.11-0.12. Sensitive high-resolution images reveal hierarchical fragmentation in dense molecular gas from the $\sim 1$ pc clump scale down to $\sim 0.01$ pc condensation scale. At each scale, the mass of the fragments is orders of magnitude larger than the Jeans mass. This is common to all four IRDC clumps we studied, suggesting that turbulence plays an important role in the early stages of clustered star formation. Masers, shock heated $\rm{NH_3}$ gas, and outflows indicate intense ongoing star formation in some cores while no such signatures are found in others. Furthermore, chemical differentiation may reflect the difference in evolutionary stages among these star formation seeds. We find $\rm{NH_3}$ ortho/para ratios of $1.1\pm0.4$, $2.0\pm0.4$, and $3.0\pm0.7$ associated with three outflows, and the ratio tends to increase along the outflows downstream. Our combined SMA and VLA observations of several IRDC clumps present the most in depth view so far of the early stages prior to the hot core phase, revealing snapshots of physical and chemical properties at various stages along an apparent evolutionary sequence.

The mid- and far-infrared view on high-mass star formation, in particular with the results from the Herschel space observatory, has shed light on many aspects of massive star formation. However, these continuum studies lack kinematic information. We study the kinematics of the molecular gas in high-mass star-forming regions. We complemented the PACS and SPIRE far-infrared data of 16 high-mass star-forming regions from the Herschel key project EPoS with N2H+ molecular line data from the MOPRA and Nobeyama 45m telescope. Using the full N2H+ hyperfine structure, we produced column density, velocity, and linewidth maps. These were correlated with PACS 70micron images and PACS point sources. In addition, we searched for velocity gradients. For several regions, the data suggest that the linewidth on the scale of clumps is dominated by outflows or unresolved velocity gradients. IRDC18454 and G11.11 show two velocity components along several lines of sight. We find that all regions with a diameter larger than 1pc show either velocity gradients or fragment into independent structures with distinct velocities. The velocity profiles of three regions with a smooth gradient are consistent with gas flows along the filament, suggesting accretion flows onto the densest regions. We show that the kinematics of several regions have a significant and complex velocity structure. For three filaments, we suggest that gas flows toward the more massive clumps are present.

Infrared-dark clouds (IRDCs) harbor the early phases of cluster and high-mass star formation and are comprised of cold (~20 K), dense (n > 10$^4$ cm$^{-3}$) gas. The spectral energy distribution (SED) of IRDCs is dominated by the far-infrared and millimeter wavelength regime, and our initial Herschel study examined IRDCs at the peak of the SED with high angular resolution. Here we present a follow-up study using the SABOCA instrument on APEX which delivers 7.8" angular resolution at 350 micron, matching the resolution we achieved with Herschel/PACS, and allowing us to characterize substructure on ~0.1pc scales. Our sample of 11 nearby IRDCs are a mix of filamentary and clumpy morphologies, and the filamentary clouds show significant hierarchical structure, while the clumpy IRDCs exhibit little hierarchical structure. All IRDCs, regardless of morphology, have about 14% of their total mass in small scale core-like structures which roughly follow a trend of constant volume density over all size scales. Out of the 89 protostellar cores we identified in this sample with Herschel, we recover 40 of the brightest and re-fit their SEDs and find their properties agree fairly well with our previous estimates (<T> ~ 19K). We detect a new population of "cold cores" which have no 70 micron counterpart, but are 100 and 160 micron-bright, with colder temperatures (<T> ~ 16K). This latter population, along with SABOCA-only detections, are predominantly low-mass objects, but their evolutionary diagnostics are consistent with the earliest starless or prestellar phase of cores in IRDCs.

Star formation in molecular clouds is intimately linked to their internal mass distribution. We present an unprecedentedly detailed analysis of the column density structure of a high-mass, filamentary molecular cloud, namely IRDC G11.11-0.12 (G11). We use two novel column density mapping techniques: high-resolution (FWHM=2", or ~0.035 pc) dust extinction mapping in near- and mid-infrared, and dust emission mapping with the Herschel satellite. These two completely independent techniques yield a strikingly good agreement, highlighting their complementarity and robustness. We first analyze the dense gas mass fraction and linear mass density of G11. We show that G11 has a top heavy mass distribution and has a linear mass density (M_l ~ 600 Msun pc^{-1}) that greatly exceeds the critical value of a self-gravitating, non-turbulent cylinder. These properties make G11 analogous to the Orion A cloud, despite its low star-forming activity. This suggests that the amount of dense gas in molecular clouds is more closely connected to environmental parameters or global processes than to the star-forming efficiency of the cloud. We then examine hierarchical fragmentation in G11 over a wide range of size-scales and densities. We show that at scales 0.5 pc > l > 8 pc, the fragmentation of G11 is in agreement with that of a self-gravitating cylinder. At scales smaller than l < 0.5 pc, the results agree better with spherical Jeans' fragmentation. One possible explanation for the change in fragmentation characteristics is the size-scale-dependent collapse time-scale that results from the finite size of real molecular clouds: at scales l < 0.5 pc, fragmentation becomes sufficiently rapid to be unaffected by global instabilities.

Aims: We study the fragmentation and dynamical properties of a massive starless gas clump at the onset of high-mass star formation. Methods: Based on Herschel continuum data we identify a massive gas clump that remains far-infrared dark up to 100mum wavelengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 45m single-dish spectral line and continuum observations. Results: The massive gas reservoir fragments at spatial scales of ~18000AU in four cores. Comparing the spatial extent of this high-mass region with intermediate- to low-mass starless cores from the literature, we find that linear sizes do not vary significantly over the whole mass regime. However, the high-mass regions squeeze much more gas into these similar volumes and hence have orders of magnitude larger densities. The fragmentation properties of the presented low-to high-mass regions are consistent with gravitational instable Jeans fragmentation. Furthermore, we find multiple velocity components associated with the resolved cores. Recent radiative transfer hydrodynamic simulations of the dynamic collapse of massive gas clumps also result in multiple velocity components along the line of sight because of the clumpy structure of the regions. This result is supported by a ratio between viral and total gas mass for the whole region <1. Conclusions: This apparently still starless high-mass gas clump exhibits clear signatures of early fragmentation and dynamic collapse prior to the formation of an embedded heating source. A comparison with regions of lower mass reveals that the linear size of star-forming regions does not necessarily have to vary much for different masses, however, the mass reservoirs and gas densities are orders of magnitude enhanced for high-mass regions compared to their lower-mass siblings.

We present a series of smoothed particle hydrodynamical models of G0.253+0.016 (also known as 'The Brick'), a very dense molecular cloud that lies close to the Galactic Centre. We explore how its gas and dust temperatures react as we vary the strength of both the interstellar radiation field (ISRF) and the cosmic ray ionisation rate (CRIR). As the physical extent of G0.253+0.016 along our line-of-sight is unknown, we consider two possibilities: one in which the longest axis is that measured in the plane of the sky (9.4 pc in length), and one in which it is along the line of sight, in which case we take it to be 17 pc. To recover the observed gas and dust temperatures, we find find that the ISRF must be around 1000 times the solar neighbourhood value, and the CRIR must be roughly 1E-14 /s, regardless of the geometries studied. For such high values of the CRIR, we find that cooling in the cloud's interior is dominated by neutral oxygen, in contrast to standard molecular clouds, which at the same densities are mainly cooled via CO. Our results suggest that the conditions near G0.253+0.016 are more extreme than those generally accepted for the inner 500 pc of the galaxy.

The temperature and density structure of molecular cloud cores are the most important physical quantities that determine the course of the protostellar collapse and the properties of the stars they form. Nevertheless, density profiles often rely either on the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, we initiated the Guaranteed Time Key Project (GTKP) "The Earliest Phases of Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, including all observations, the analysis method, and the initial results of the survey. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm continuum maps at up to eight wavelengths between 100 micron and 1.2 mm. Our sample contains both globules with starless cores and embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which are then fit independently with modified blackbody curves to obtain line-of-sight-averaged dust temperature and column density maps. We find that the thermal structure of all globules is dominated by external heating from the interstellar radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense interiors (8-12 K). The protostars embedded in some of the globules raise the local temperature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the globules.

We present continuum observations of the infrared dark cloud (IRDC) G48.66-0.22 (G48) obtained with Herschel, Spitzer, and APEX, in addition to several molecular line observations. The Herschel maps are used to derive temperature and column density maps of G48 using a model based on a modified blackbody. We find that G48 has a relatively simple structure and is relatively isolated; thus this IRDC provides an excellent target to study the collapse and fragmentation of a filamentary structure in the absence of complicating factors such as strong external feedback. The derived temperature structure of G48 is clearly non-isothermal from cloud to core scale. The column density peaks are spatially coincident with the lowest temperatures (~ 17.5 K) in G48. A total cloud mass of ~390Msun is derived from the column density maps. By comparing the luminosity-to-mass ratio of 13 point sources detected in the Herschel/PACS bands to evolutionary models, we find that two cores are likely to evolve into high-mass stars (M>8 Msun). The derived mean projected separation of point sources is smaller than in other IRDCs but in good agreement with theoretical predications for cylindrical collapse. We detect several molecular species such as CO, HCO+, HCN, HNC and N2H+. CO is depleted by a factor of ~3.5 compared to the expected interstellar abundance, from which we conclude that CO freezes out in the central region. Furthermore, the molecular clumps, associated with the sub-millimeter peaks in G48, appear to be gravitationally unbound or just pressure confined. The analysis of critical line masses in G48 show that the entire filament is collapsing, overcoming any internal support.

(Abriged) In the framework of the Herschel GTKP "The earliest phases of star formation", we have imaged B68 between 100 and 500 um. Ancillary (sub)mm data, spectral line maps of the 12/13CO(2-1) transitions as well as a NIR extinction map were added to the analysis. We employed a ray-tracing algorithm to derive the 2D mid-plane dust temperature and volume density distribution without suffering from LoS averaging effects of simple SED fitting procedures. Additional 3D radiative transfer calculations were employed to investigate the connection between the external irradiation and the peculiar crescent shaped morphology found in the FIR maps. For the first time, we spatially resolve the dust temperature and density distribution of B68. We find T_dust dropping from 16.7 K at the edge to 8.2 K in the centre, which is about 4 K lower than the result of the simple SED fitting approach. N_H peaks at 4.3x10^22 cm^-2 and n_H at 3.4x10^5 cm^-3 in the centre. B68 has a mass of 3.1 M_sun of material with A_K > 0.2 mag for an assumed distance of 150 pc. We detect a compact source in the southeastern trunk, which is also seen in extinction and CO. We find the radial density distribution from the edge of the inner plateau outward to be n_H ~ r^-3.5. Such a steep profile can arise from either or both of the following: external irradiation with a significant UV contribution or the fragmentation of filamentary structures. Our 3D radiative transfer model of an externally irradiated core by an anisotropic ISRF reproduces the crescent morphology. Our CO observations show that B68 is part of a chain of globules in both space and velocity, which may indicate that it was once part of a filament which dispersed. We also resolve a new compact source in the SE trunk and find that it is slightly shifted in centroid velocity from B68, lending qualitative support to core collision scenarios.

(Abridged) We present an overview of the sample of high-mass star and cluster forming regions observed as part of the Earliest Phases of Star Formation (EPoS) Herschel Guaranteed Time Key Program. A sample of 45 infrared-dark clouds (IRDCs) were mapped at PACS 70, 100, and 160 micron and SPIRE 250, 350, and 500 micron. In this paper, we characterize a population of cores which appear in the PACS bands and place them into context with their host cloud and investigate their evolutionary stage. We construct spectral energy distributions (SEDs) of 496 cores which appear in all PACS bands, 34% of which lack counterparts at 24 micron. From single-temperature modified blackbody fits of the SEDs, we derive the temperature, luminosity, and mass of each core. These properties predominantly reflect the conditions in the cold, outer regions. Taking into account optical depth effects and performing simple radiative transfer models, we explore the origin of emission at PACS wavelengths. The core population has a median temperature of 20K and has masses and luminosities that span four to five orders of magnitude. Cores with a counterpart at 24 micron are warmer and bluer on average than cores without a 24 micron counterpart. We conclude that cores bright at 24 micron are on average more advanced in their evolution, where a central protostar(s) have heated the outer bulk of the core, than 24 micron-dark cores. The 24 micron emission itself can arise in instances where our line of sight aligns with an exposed part of the warm inner core. About 10% of the total cloud mass is found in a given cloud's core population. We uncover over 300 further candidate cores which are dark until 100 micron. These are candidate starless objects, and further observations will help us determine the nature of these very cold cores.

Infrared-dark clouds (IRDCs) are believed to be the birthplaces of rich clusters and thus contain the earliest phases of high-mass star formation. We use the Green Bank Telescope (GBT) and Very Large Array (VLA) maps of ammonia (NH3) in six IRDCs to measure their column density and temperature structure (Paper 1), and here, we investigate the kinematic structure and energy content. We find that IRDCs overall display organized velocity fields, with only localized disruptions due to embedded star formation. The local effects seen in NH3 emission are not high velocity outflows but rather moderate (few km/s) increases in the line width that exhibit maxima near or coincident with the mid-infrared emission tracing protostars. These line width enhancements could be the result of infall or (hidden in NH3 emission) outflow. Not only is the kinetic energy content insufficient to support the IRDCs against collapse, but also the spatial energy distribution is inconsistent with a scenario of turbulent cloud support. We conclude that the velocity signatures of the IRDCs in our sample are due to active collapse and fragmentation, in some cases augmented by local feedback from stars.

Context: The earliest stages of high-mass star formation are still poorly characterized. It is also unknown whether the initial conditions vary with environment. Aims: We want to investigate the youngest massive gas clumps in the environment of extremely active star formation. Methods: We selected the IRDC18454 complex associated with the W43 Galactic mini-starburst, and observed it in the continuum emission between 70mu and 1.2mm, and in spectral line emission of N2H+ and 13CO with single-dish instruments and an interferometer. Results: The multi-wavelength continuum study allows us to identify clumps that are infrared dark even at 70mu and hence the best candidates to be genuine high-mass starless gas clumps. The spectral energy distributions reveal elevated temperatures and luminosities compared to more quiescent environments. Furthermore, we identify a temperature gradient from the W43 mini-starburst toward the starless clumps. We discuss whether the radiation impact of the nearby mini-starburst changes the fragmentation properties of the gas clumps and by that maybe favors more high-mass star formation in such an environment. The spectral line data reveal two different velocity components of the gas at 100 and 50km/s. While chance projection is a possibility to explain these components, the projected associations of the emission sources as well as the prominent location at the Galactic bar -- spiral arm interface also allow the possibility that these two components may be spatially associated and even interacting. Conclusions: High-mass starless gas clumps can exist in the close environment of very active star formation without being destroyed. The impact of the active star formation sites may even allow for more high-mass stars to form in these 2nd generation gas clumps.

We present Very Large Array observations of NH3 (1,1) and (2,2) in a sample of six infrared-dark clouds (IRDCs) with distances from 2 to 5 kpc. We find that ammonia serves as an excellent tracer of dense gas in IRDCs, showing no evidence of depletion, and the average abundance in these clouds is 8.1 x 10^-7 relative to H2. Our sample consists of four IRDCs with 24 micron embedded protostars and two that appear starless. We calculate the kinetic temperature of the gas in IRDCs and find no significant difference between starless and star-forming IRDCs. We find that the bulk of the gas is between 8 and 13 K, indicating that any embedded or nearby stars or clusters do not affect the gas temperature dramatically. Though IRDCs have temperatures and volume densities on par with local star formation regions of lower mass, they consist of much more mass which induces very high internal pressures. In order for IRDCs to survive as coherent structures, the internal pressure must be balanced by a confining pressure provided by the high concentration of molecular clouds in the spiral arm in which they reside. The high molecular concentration and pressure is roughly consistent with gas dynamics of a bar galaxy.

Infrared-dark clouds (IRDCs) are the precursors to massive stars and stellar clusters. G011.11-0.12 is a well-studied filamentary IRDC, though, to date, the absence of far-infrared data with sufficient spatial resolution has limited the understanding of the structure and star-formation activity. We use Herschel to study the embedded population of young pre- and protostellar cores in this IRDC. We examine the cloud structure, which appears in absorption at short wavelength and in emission at longer wavelength. We derive the properties of the massive cores from the spectral energy distributions of bright far-infrared point sources detected with the PACS instrument aboard Herschel. We report on the detection and characterization of pre- and protostellar cores in a massive filamentary infrared-dark cloud G011.11-0.12 using PACS. We characterize 18 cores directly associated with the filament, two of which have masses over 50 Msun, making them the best candidates to become massive stars in G011.11-0.12. These cores are likely at various stages of protostar formation, showing elevated temperature (<T> ~ 22 K) with respect to the ambient gas reservoir. The core masses (<M> ~ 24 Msun) are small compared to that in the cold filament. The mean core separation is 0.9 pc, well in excess of the Jeans length in the filament. We confirm that star formation in IRDCs is underway and diverse, and IRDCs have the capability of forming massive stars and clusters.

Using a source selection biased towards high mass star forming regions, we used a Large Velocity Gradient (LVG) code to calculate the H2 densities and CS column densities for a sample of Midcourse Space Experiment (MSX) 8 micron infrared dark cores. Our average H2 density and CS column density were 1.14 x 10e6 cm-3 and 1.21 x 10e13 cm-2 respectively. In addition, we have calculated the Jeans mass and Virial mass for each core to get a better understanding of their gravitational stability. We found that core masses calculated from observations of N2H+ J = 1-0 and C18O J = 1-0 by Ragan et al. 2006 (Paper 1) were sufficient for collapse, though most regions are likely to form protoclusters. We have explored the star-forming properties of the molecular gas within our sample and find some diversity which extends the range of infrared dark clouds from very the massive clouds that will create large clusters, to clouds that are similar to some of our local counterparts (e.g. Serpens, Ophiuchus).

We have conducted a survey of a sample of infrared-dark clouds (IRDCs) with the Spitzer Space Telescope in order to explore their mass distribution. We present a method for tracing mass using dust absorption against the bright Galactic background at 8 microns. The IRDCs in this sample are comprised of tens of clumps, ranging in sizes from 0.02 to 0.3 pc in diameter and masses from 0.5 to a few 10 Msun, the broadest dynamic range in any clump mass spectrum study to date. Structure with this range in scales confirms that IRDCs are the the precursors to stellar clusters in an early phase of fragmentation. Young stars are distributed in the vicinity of the IRDCs, but the clumps are typically not associated with stars and appear pre-stellar in nature. We find an IRDC clump mass spectrum with a slope of 1.76 +/- 0.05 for masses from 30 to 3000 Msun. This slope is consistent with numerous studies, culled from a variety of observational techniques, of massive star formation regions and is close to the mass function of Galactic stellar clusters and star clusters in other galaxies. We assert that the shape of the mass function is an intrinsic and universal feature of massive star formation regions, that are the birth sites of stellar clusters. As these clouds evolve and their constituent clumps fragment, the mass spectrum will steepen and eventually assume the form of the core mass function that is observed locally.

We have identified 41 infrared dark clouds from the 8 micron maps of the Midcourse Space Experiment (MSX), selected to be found within one square degree areas centered on known ultracompact HII regions. We have mapped these infrared dark clouds in N2H+(1-0), CS(2-1) and C18O(1-0) emission using the Five College Radio Astronomy Observatory. The maps of the different species often show striking differences in morphologies, indicating differences in evolutionary state and/or the presence of undetected, deeply embedded protostars. We derive an average mass for these clouds using N2H+ column densities of ~2500 solar masses, a value comparable to that found in previous studies of high mass star forming cores using other mass tracers. The linewidths of these clouds are typically ~2.0 - 2.9 km/s. Based on the fact that they are dark at 8 micron, compact, massive, and have large velocity dispersions, we suggest that these clouds may be the precursor sites of intermediate and high mass star formation.