The mass growth of protostars is a central element to the determination of fundamental stellar population properties such as the initial mass function. Constraining the accretion history of individual protostars is therefore an important aspect of star formation research. The goal of the study presented here is to determine whether high-mass (proto)stars gain their mass from a compact (<0.1pc) fixed-mass reservoir of gas, often referred to as dense cores, in which they are embedded, or whether the mass growth of high-mass stars is governed by the dynamical evolution of the parsec-scale clump that typically surrounds them. To achieve this goal, we performed a 350micron continuum mapping of 11 infrared dark clouds, along side some of their neighbouring clumps, with the ArT\'eMiS camera on APEX. By identifying about 200 compact ArT\'eMiS sources, and matching them with Herschel Hi-GAL 70micron sources, we have been able to produce mass vs. temperature diagrams. We compare the nature (i.e. starless or protostellar) and location of the ArT\'eMiS sources in these diagrams with modelled evolutionary tracks of both core-fed and clump-fed accretion scenarios. We argue that the latter provide a better agreement with the observed distribution of high-mass star-forming cores. However, a robust and definitive conclusion on the question of the accretion history of high-mass stars requires larger number statistics.

In this work we investigate the interplay between gravity and turbulence at different spatial scales and in different density regimes. We analyze a sample of 70 $\mu$m quiet clumps that are divided into three surface density bins and we compare the dynamics of each group with the dynamics of their respective filaments. The densest clumps form within the densest filaments on average, and they have the highest value of the velocity dispersion. The kinetic energy is transferred from the filaments down to the clumps most likely through a turbulent cascade, but we identify a critical value of the surface density, $\Sigma\simeq0.1$ g cm$^{2}$, above which the dynamics changes from being mostly turbulent-driven to mostly gravity-driven. The scenario we obtain from our data is a continuous interplay between turbulence and gravity, where the former creates structures at all scales and the latter takes the lead when the critical surface density threshold is reached. In the densest filaments this transition can occur at the parsec, or even larger scales, leading to a global collapse of the whole region and most likely to the formation of the massive objects.

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.

(Abbreviated) We investigate the interplay between magnetic (B) field, gravity, and turbulence in the fragmentation process of cores within the filamentary infrared dark cloud G34.43+00.24. We observe the magnetic field (B) morphology across G34.43 and compare with the kinematics obtained from N2H+ across the filament. We derive local velocity gradients from N2H+, tracing motion in the plane of sky, and compare with the observed local B field orientations in the plane of sky. Besides a large-scale east-west velocity gradient, we find a close alignment between local B field orientations and local velocity gradients toward the MM1/MM2 ridge. This local correlation in alignment suggests that gas motions are influenced by the magnetic field morphology or vice versa. Additionally, this alignment seems to be getting even closer with increasing integrated emission in N2H+, possibly indicating that a growing gravitational pull is more and more aligning B field and gas motion. We analyze and quantify B field, gravity, turbulence, and their relative importance toward the MM1, MM2 and MM3 regions with various techniques over two scales, a larger clump area at 2 pc scale and the smaller core area at 0.6 pc scale. While gravitational energy, B field, and turbulent pressure all grow systematically from large to small scale, the ratios among the three constituents develop clearly differently over scale. We propose that this varying relative importance between B field, gravity, and turbulence over scale drives and explains the different fragmentation types seen at sub-pc scale (no fragmentation in MM1; aligned fragmentation in MM2; clustered fragmentation in MM3). We discuss uncertainties, subtleties, and the robustness of our conclusion, and we stress the need of a multi-scale joint analysis to understand the dynamics in these systems.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

We define the molecular cloud properties of the Milky Way first quadrant using data from the JCMT CO(3-2) High Resolution Survey. We apply the Spectral Clustering for Interstellar Molecular Emission Segmentation (SCIMES) algorithm to extract objects from the full-resolution dataset, creating the first catalog of molecular clouds with a large dynamic range in spatial scale. We identify $>85\,000$ clouds with two clear sub-samples: $\sim35\,500$ well-resolved objects and $\sim540$ clouds with well-defined distance estimations. Only 35% of the cataloged clouds (as well as the total flux encompassed by them) appear enclosed within the Milky Way spiral arms. The scaling relationships between clouds with known distances are comparable to the characteristics of the clouds identified in previous surveys. However, these relations between integrated properties, especially from the full catalog, show a large intrinsic scatter ($\sim0.5$ dex), comparable to other cloud catalogs of the Milky Way and nearby galaxies. The mass distribution of molecular clouds follows a truncated-power law relationship over three orders of magnitude in mass with a form $dN/dM \propto M^{-1.7}$ with a clearly defined truncation at an upper mass of $M_0 \sim 3 \times 10^6~M_\odot$, consistent with theoretical models of cloud formation controlled by stellar feedback and shear. Similarly, the cloud population shows a power-law distribution of size with $dN/dR \propto R^{-2.8}$ with a truncation at $R_0 = 70$ pc.

Aims. To constrain the physical processes that lead to the birth of high-mass stars it is mandatory to study the very first stages of their formation. We search for high-mass analogs of low-mass prestellar cores in W43-MM1. Methods. We conducted a 1.3 mm ALMA mosaic of the complete W43-MM1 cloud, which has revealed numerous cores with ~ 2000 au FWHM sizes. We investigated the nature of cores located at the tip of the main filament, where the clustering is minimum. We used the continuum emission to measure the core masses and the $^{13}$CS(5-4) line emission to estimate their turbulence level. We also investigated the prestellar or protostellar nature of these cores by searching for outflow signatures traced by CO(2-1) and SiO(5-4) line emission, and for molecular complexity typical of embedded hot cores. Results. Two high-mass cores of ~ 1300 au diameter and ~ $60~M_\odot$ mass are observed to be turbulent but gravitationally bound. One drives outflows and is associated with a hot core. The other core, W43-MM1\#6, does not yet reveal any star formation activity and thus is an excellent high-mass prestellar core candidate.

The JCMT Gould Belt Survey was one of the first Legacy Surveys with the James Clerk Maxwell Telescope in Hawaii, mapping 47 square degrees of nearby (< 500 pc) molecular clouds in both dust continuum emission at 850 $\mu$m and 450 $\mu$m, as well as a more-limited area in lines of various CO isotopologues. While molecular clouds and the material that forms stars have structures on many size scales, their larger-scale structures are difficult to observe reliably in the submillimetre regime using ground-based facilities. In this paper, we quantify the extent to which three subsequent data-reduction methods employed by the JCMT GBS accurately recover emission structures of various size scales, in particular, dense cores which are the focus of many GBS science goals. With our current best data-reduction procedure, we expect to recover 100% of structures with Gaussian sigma sizes of $\le$30" and intensity peaks of at least five times the local noise for isolated peaks of emission. The measured sizes and peak fluxes of these compact structures are reliable (within 15% of the input values), but source recovery and reliability both decrease significantly for larger emission structures and for fainter peaks. Additional factors such as source crowding have not been tested in our analysis. The most recent JCMT GBS data release includes pointing corrections, and we demonstrate that these tend to decrease the sizes and increase the peak intensities of compact sources in our dataset, mostly at a low level (several percent), but occasionally with notable improvement.

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.

Molecular-line observations of filaments in star-forming regions have revealed the existence of elongated coherent features within the filaments; these features are termed fibres. Here we caution that, since fibres are traced in PPV space, there is no guarantee that they represent coherent features in PPP space. We illustrate this contention using simulations of the growth of a filament from a turbulent medium. Synthetic C$^{18}$O observations of the simulated filaments reveal the existence of fibres very similar to the observed ones, i.e. elongated coherent features in the resulting PPV data-cubes. Analysis of the PPP data-cubes (i.e. 3D density fields) also reveals elongated coherent features, which we term sub-filaments. Unfortunately there is very poor correspondence between the fibres and the sub-filaments in the simulations. Both fibres and sub-filaments derive from inhomogeneities in the turbulent accretion flow onto the main filament. As a consequence, fibres are often affected by line-of-sight confusion. Similarly, sub-filaments are often affected by large velocity gradients, and even velocity discontinuities. These results suggest that extreme care should be taken when using velocity coherent features to constrain the underlying substructure within a filament.

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

We tested the validity of the three Larson relations in a sample of 213 massive clumps selected from the Herschel Hi-GAL survey and combined with data from the MALT90 survey of 3mm emission lines. The clumps have been divided in 5 evolutionary stages to discuss the Larson relations also as function of evolution. We show that this ensemble does not follow the three Larson relations, regardless of clump evolutionary phase. A consequence of this breakdown is that the virial parameter $\alpha_{vir}$ dependence with mass (and radius) is only a function of the gravitational energy, independent of the kinetic energy of the system, and $\alpha_{vir}$ is not a good descriptor of clump dynamics. Our results suggest that clumps with clear signatures of infall motions are statistically indistinguishable from clumps with no such signatures. The observed non-thermal motions are not necessarily ascribed to turbulence acting to sustain the gravity, but they may be due to the gravitational collapse at the clump scales. This seems particularly true for the most massive (M$\geq$1000 M$_{\odot}$) clumps in the sample, where also exceptionally high magnetic fields may not be enough to stabilize the collapse.

We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.

How important is the magnetic (B-) field when compared to gravity and turbulence in the star-formation process? Does its importance depend on scale and location? We summarize submm dust polarization observations towards the large filamentary infrared dark cloud G34 and towards a dense core in the high-mass star-forming region W51. We detect B-field orientations that are either perpendicular or parallel to the G34 filament axis. These B-field orientations further correlate with local velocity gradients. Towards three cores in G34 we find a varying importance between B-field, gravity, and turbulence that seems to dictate varying types of fragmentation. At highest resolution towards the gravity-dominated collapsing core W51 e2 we resolve new B-field features, such as converging B-field lines and possibly magnetic channels.

Converging networks of interstellar filaments i.e. hubs, have been recently linked to the formation of stellar clusters and massive stars. The goal is to study the kinematic and density structure of the SDC13 hub at high angular resolution to determine what drives its evolution and fragmentation. We have mapped SDC13, a 1000Msun infrared dark hub, in NH3(1,1) and NH3(2,2) emission lines, with both the JVLA and GBT down to 0.07pc. The mass-per-unit-lengths of all four hub filaments are thermally super-critical, consistent with the presence of tens of gravitationally bound cores along them. These cores exhibit regular separation of 0.37 +/- 0.16 pc suggesting gravitational instabilities running along these super-critical filaments are responsible for their fragmentation. The observed local increase of the dense gas velocity dispersion towards starless cores is believed to be a consequence of such fragmentation process. We see velocity gradient peaks towards 63% of the cores as expected during the early stages of filament fragmentation. The most massive cores are located at the filament junctions, where the velocity dispersion is the largest. We interpret this as the result of the hub morphology generating the largest acceleration gradients near the hub centre. We propose a scenario for the evolution of the SDC13 hub in which filaments first form as post-shock structures in a supersonic turbulent flow. Then gravity takes over and starts shaping the evolution of the hub, both fragmenting filaments and pulling the gas towards the centre of the gravitational well. By doing so, gravitational energy is converted into kinetic energy in both local (cores) and global (hub centre) potential well minima. Furthermore, the generation of larger gravitational acceleration gradients at the filament junctions promotes the formation of more massive cores. [abridged]

The origin of non-thermal motions in massive star forming regions can be ascribed to turbulence acting against the gravitational collapse, or to the self-gravity itself driving the rapid global collapse. The dependence between velocity dispersion, radius and clouds surface density found by Heyer et al. (2009), $\sigma/R^{1/2}\propto \Sigma^{1/2}$, has been interpreted in terms of global collapse of clouds. In this work we demonstrate that this relation is an expression of a more general relation between accelerations. We introduce the gravo-turbulent acceleration, a$_k$, which describe the non-thermal motions in each region, and the acceleration generated by the gravitational field a$_G$, which is proportional to $\Sigma$. We also introduce a new coefficient, the force partition coefficient $\alpha_{for}$ which is equivalent to the virial parameter but does not distinguish between collapsing and non-collapsing regions. In this work we use the a$_k$ - a$_G$ formalism in the analysis of a new sample of 16 massive starless clumps (MSCls) combined with data from the literature. We show that a$_k$ and a$_G$ are not independent. The non-thermal motions in each region can originate from both local turbulence and self-gravity but overall the data in the a$_k$ vs. a$_G$ diagram demonstrate that the majority of the non-thermal motions originate from self-gravity. We further show that all the MSCls with $\Sigma\geq 0.1$ g cm$^{-2}$ show signs of infall motions, a strong indication that the denser regions are the first to collapse. Finally, we include in the formalism the contribution of an external pressure and the magnetic fields.

The dynamic activity in massive star forming regions prior to the formation of bright protostars is still not fully investigated. In this work we present observations of HCO+ J=1-0 and N2H+ J=1-0 made with the IRAM 30m telescope towards a sample of 16 Herschel-identified massive 70 micron quiet clumps associated with infrared dark clouds. The clumps span a mass range from 300 M_sun to 2000 M_sun. The N2H+ data show that the regions have significant non-thermal motions with velocity dispersion between 0.28 km s^-1 and 1.5 km s^-1, corresponding to Mach numbers between 2.6 and 11.5. The majority of the 70 micron quiet clumps have asymmetric HCO+ line profiles, indicative of significant dynamical activity. We show that there is a correlation between the degree of line asymmetry and the surface density Sigma of the clumps, with clumps of Sigma>=0.1 g cm^-2 having more asymmetric line profiles, and so are more dynamically active, than clumps with lower Sigma. We explore the relationship between velocity dispersion, radius and Sigma and show how it can be interpreted as a relationship between an acceleration generated by the gravitational field a_G, and the measured kinetic acceleration, a_k, consistent with the majority of the non-thermal motions originating from self-gravity. Finally, we consider the role of external pressure and magnetic fields in the interplay of forces.

In recent years there has been a growing interest in studying giant molecular filaments (GMFs), which are extremely elongated (> 100pc in length) giant molecular clouds (GMCs). They are often seen as inter-arm features in external spiral galaxies, but have been tentatively associated with spiral arms when viewed in the Milky Way. In this paper, we study the time evolution of GMFs in a high-resolution section of a spiral galaxy simulation, and their link with spiral arm GMCs and star formation, over a period of 11Myrs. The GMFs generally survive the inter-arm passage, although they are subject to a number of processes (e.g. star formation, stellar feedback and differential rotation) which can break the giant filamentary structure into smaller sections. The GMFs are not gravitationally bound clouds as a whole, but are, to some extent, confined by external pressure. Once they reach the spiral arms, the GMFs tend to evolve into more substructured spiral arm GMCs, suggesting that GMFs may be precursors to arm GMCs. Here, they become incorporated into the more complex and almost continuum molecular medium that makes up the gaseous spiral arm. Instead of retaining a clear filamentary shape, their shapes are distorted both by their climb up the spiral potential and their interaction with the gas within the spiral arm. The GMFs do tend to become aligned with the spiral arms just before they enter them (when they reach the minimum of the spiral potential), which could account for the observations of GMFs in the Milky Way.

Massive clumps, prior to the formation of any visible protostars, are the best candidates to search for the elusive massive starless cores. In this work we investigate the dust and gas properties of massive clumps selected to be 70 micron quiet, therefore good starless candidates. Our sample of 18 clumps has masses 300 < M < 3000 M_sun, radius 0.54 < R < 1.00 pc, surface densities Sigma > 0.05 g cm^-2 and luminosity/mass ratio L/M < 0.3. We show that half of these 70 micron quiet clumps embed faint 24 micron sources. Comparison with GLIMPSE counterparts shows that 5 clumps embed young stars of intermediate stellar mass up to ~5.5 M_sun. We study the clump dynamics with observations of N2H+ (1-0), HNC (1-0) and HCO+ (1-0) made with the IRAM 30m telescope. Seven clumps have blue-shifted spectra compatible with infall signatures, for which we estimate a mass accretion rate 0.04 < M_dot < 2.0 x 10^-3 M_sun yr^-1, comparable with values found in high-mass protostellar regions, and free-fall time of the order of t_ff = 3 x 10^5 yr. The only appreciable difference we find between objects with and without embedded 24 micron sources is that the infall rate appears to increase from 24 micron dark to 24 micron bright objects. We conclude that all 70 micron quiet objects have similar properties on clump scales, independently of the presence of an embedded protostar. Based on our data we speculate that the majority, if not all of these clumps may already embed faint, low-mass protostellar cores. If these clumps are to form massive stars, this must occur after the formation of these lower mass stars.

We present the results of smoothed particle hydrodynamic simulations investigating the evolution and fragmentation of filaments that are accreting from a turbulent medium. We show that the presence of turbulence, and the resulting inhomogeneities in the accretion flow, play a significant role in the fragmentation process. Filaments which experience a weakly turbulent accretion flow fragment in a two-tier hierarchical fashion, similar to the fragmentation pattern seen in the Orion Integral Shaped Filament. Increasing the energy in the turbulent velocity field results in more sub-structure within the filaments, and one sees a shift from gravity-dominated fragmentation to turbulence-dominated fragmentation. The sub-structure formed in the filaments is elongated and roughly parallel to the longitudinal axis of the filament, similar to the fibres seen in observations of Taurus, and suggests that the fray and fragment scenario is a possible mechanism for the production of fibres. We show that the formation of these fibre-like structures is linked to the vorticity of the velocity field inside the filament and the filament's accretion from an inhomogeneous medium. Moreover, we find that accretion is able to drive and sustain roughly sonic levels of turbulence inside the filaments, but is not able to prevent radial collapse once the filaments become supercritical. However, the supercritical filaments which contain fibre-like structures do not collapse radially, suggesting that fibrous filaments may not necessarily become radially unstable once they reach the critical line-density.

The early evolution of massive cluster progenitors is poorly understood. We investigate the fragmentation properties from 0.3 pc to 0.06 pc scales of a homogenous sample of infrared-quiet massive clumps within 4.5 kpc selected from the ATLASGAL survey. Using the ALMA 7m array we detect compact dust continuum emission towards all targets, and find that fragmentation, at these scales, is limited. The mass distribution of the fragments uncovers a large fraction of cores above 40 $M_\odot$, corresponding to massive dense cores (MDCs) with masses up to ~400 $M_\odot$. 77 % of the clumps contain at most 3 MDCs per clump, and we also reveal single clumps/MDCs. The most massive cores are formed within the more massive clumps, and a high concentration of mass on small scales reveals a high core formation efficiency. The mass of MDCs highly exceeds the local thermal Jeans-mass, and observational evidence is lacking for a sufficiently high level of turbulence or strong enough magnetic fields to keep the most massive MDCs in equilibrium. If already collapsing, the observed fragmentation properties with a high core formation efficiency are consistent with the collapse setting in at parsec scales.

We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.

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".

We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.

Previous literature suggests that the densest structures in the interstellar medium form through colliding flows, but patent evidence of this process is still missing. Recent literature proposes using SiO line emission to trace low-velocity shocks associated with cloud formation through collision. In this paper we investigate the bright and extended SiO(2-1) emission observed along the ~5 pc-long W43-MM1 ridge to determine its origin. We used high angular resolution images of the SiO(2-1) and HCN(1-0) emission lines obtained with the IRAM plateau de Bure (PdBI) interferometer and combined with data from the IRAM 30 m radiotelescope. These data were complemented by a Herschel column density map of the region. We performed spectral analysis of SiO and HCN emission line profiles to identify protostellar outflows and spatially disentangle two velocity components associated with low- and high-velocity shocks. Then, we compared the low-velocity shock component to a dedicated grid of one-dimensional (1D) radiative shock models. We find that the SiO emission originates from a mixture of high-velocity shocks caused by bipolar outflows and low-velocity shocks. Using SiO and HCN emission lines, we extract seven bipolar outflows associated with massive dense cores previously identified within the W43-MM1 mini-starburst cluster. Comparing observations with dedicated Paris-Durham shock models constrains the velocity of the low-velocity shock component from 7 to 12km/s. The SiO arising from low-velocity shocks spreads along the complete length of the ridge. Its contribution represents at least 45% and up to 100% of the total SiO emission depending on the area considered. The low-velocity component of SiO is most likely associated with the ridge formation through colliding flows or cloud-cloud collision.

We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.

The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.

We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.

The radiative feedback of massive stars on molecular clouds creates pillars, globules and other features at the interface between the HII region and molecular cloud. We present here Herschel observations between 70 and 500 micron of the immediate environment of the Cygnus OB2 association, performed within the HOBYS program. All structures were detected based on their appearance at 70 micron, and have been classified as pillars, globules, evaporating gasous globules (EGGs), proplyd-like objects, and condensations. From the 70 and 160 micron flux maps, we derive the local FUV field on the PDR surfaces. In parallel, we use a census of the O-stars to estimate the overall FUV-field, that is 10^3-10^4 G_0 close to the central OB cluster (within 10 pc) and decreases down to a few tens G_0, in a distance of 50 pc. From a SED fit to the four longest Herschel wavelengths, we determine column density and temperature maps and derive masses, volume densities and surface densities for these structures. We find that the morphological classification corresponds to distinct physical properties. Pillars and globules have the longest estimated photoevaporation lifetimes, a few 10^6 yr, while all other features should survive less than that. These lifetimes are consistent with that found in simulations of turbulent, UV-illuminated clouds. We propose a tentative evolutionary scheme in which pillars can evolve into globules, which in turn then evolve into EGGs, condensations and proplyd-like objects.

We study the properties of giant molecular clouds (GMCs) from a smoothed particle hydrodynamics simulation of a portion of a spiral galaxy, modelled at high resolution, with robust representations of the physics of the interstellar medium. We examine the global molecular gas content of clouds, and investigate the effect of using CO or H2 densities to define the GMCs. We find that CO can reliably trace the high-density H2 gas, but misses less dense H2 clouds. We also investigate the effect of using 3D CO densities versus CO emission with an observer's perspective, and find that CO-emission clouds trace well the peaks of the actual GMCs in 3D, but can miss the lower density molecular gas between density peaks which is often CO-dark. Thus, the CO emission typically traces smaller clouds within larger GMC complexes. We also investigate the effect of the galactic environment (in particular the presence of spiral arms), on the distribution of GMC properties, and we find that the mean properties are similar between arm and inter-arm clouds, but the tails of some distributions are indicative of intrinsic differences in the environment. We find highly filamentary clouds (similar to the giant molecular filaments of our Galaxy) exclusively in the inter-arm region, formed by galactic shear. We also find that the most massive GMC complexes are located in the arm, and that as a consequence of more frequent cloud interactions/mergers in the arm, arm clouds are more sub-structured and have higher velocity dispersions than inter-arm clouds.

The role played by gravity in the transfer of interstellar matter from molecular cloud scales to protostellar scales is still highly debated. Only detailed studies on the kinematics of large samples of star-forming clouds will settle the issue. We present new IRAM 30m observations of a sample of 27 infrared dark clouds covering a large range of sizes, masses, and aspect ratios. Preliminary results suggest that gravity is regulating the dynamical evolution of these clouds on a couple of parsec scales.

The origin of massive stars is a fundamental open issue in modern astrophysics. Pre-ALMA interferometric studies reveal precursors to early B to late O type stars with collapsing envelopes of 15-20 M$_\odot$ on 1000-3000 AU size-scales. To search for more massive envelopes we selected the most massive nearby young clumps from the ATLASGAL survey to study their protostellar content with ALMA. Our first results using the intermediate scales revealed by the ALMA ACA array providing 3-5" angular resolution, corresponding to $\sim$0.05-0.1 pc size-scales, reveals a sample of compact objects. These massive dense cores are on average two-times more massive than previous studies of similar types of objects. We expect that once the full survey is completed, it will provide a comprehensive view on the origin of the most massive stars.

We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.

The statistical description of Giant Molecular Cloud (GMC) properties relies heavily on the performance of automatic identification algorithms, which are often seriously affected by the survey design. The algorithm we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation), is able to overcome some of these limitations by considering the cloud segmentation problem in the broad framework of the graph theory. The application of the code on the CO(3-2) High Resolution Survey (COHRS) data allowed for a robust decomposition of more than 12,000 objects in the Galactic Plane. Together with the wealth of Galactic Plane surveys of the recent years, this approach will help to open the door to a future, systematic cataloging of all discrete molecular features of our own Galaxy.

We present a generalization of the Giant Molecular Cloud (GMC) identification problem based on cluster analysis. The method we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation) considers the dendrogram of emission in the broader framework of graph theory and utilizes spectral clustering to find discrete regions with similar emission properties. For Galactic molecular cloud structures, we show that the characteristic volume and/or integrated CO luminosity are useful criteria to define the clustering, yielding emission structures that closely reproduce "by-eye" identification results. SCIMES performs best on well-resolved, high-resolution data, making it complementary to other available algorithms. Using 12CO(1-0) data for the Orion-Monoceros complex, we demonstrate that SCIMES provides robust results against changes of the dendrogram-construction parameters, noise realizations and degraded resolution. By comparing SCIMES with other cloud decomposition approaches, we show that our method is able to identify all canonical clouds of the Orion-Monoceros region, avoiding the over-division within high resolution survey data that represents a common limitation of several decomposition algorithms. The Orion-Monoceros objects exhibit hierarchies and size-line width relationships typical to the turbulent gas in molecular clouds, although "the Scissors" region deviates from this common description. SCIMES represents a significant step forward in moving away from pixel-based cloud segmentation towards a more physical-oriented approach, where virtually all properties of the ISM can be used for the segmentation of discrete objects.

Our goal is to identify bipolar HII regions and to understand their morphology, their evolution, and the role they play in the formation of new generations of stars. We use the Spitzer and Herschel Hi-GAL surveys to identify bipolar HII regions. We search for their exciting star(s) and estimate their distances using near-IR data. Dense clumps are detected using Herschel-SPIRE data. MALT90 observations allow us to ascertain their association with the central HII region. We identify Class 0/I YSOs using their Spitzer and Herschel-PACS emissions. These methods will be applied to the entire sample of candidate bipolar HII regions. This paper focuses on two bipolar HII regions, one interesting in terms of its morphology, G319.88$+$00.79, and one in terms of its star formation, G010.32$-$00.15. Their exciting clusters are identified and their photometric distances estimated to be 2.6 kpc and 1.75 kpc, respectively. We suggest that these regions formed in dense and flat structures that contain filaments. They have a central ionized region and ionized lobes perpendicular to the parental cloud. The remains of the parental cloud appear as dense (more than 10^4 per cm^3) and cold (14-17 K) condensations. The dust in the PDR is warm (19-25 K). Dense massive clumps are present around the central ionized region. G010.32-00.14 is especially remarkable because five clumps of several hundred solar masses surround the central HII region; their peak column density is a few 10^23 per cm^2, and the mean density in their central regions reaches several 10^5 per cm^3. Four of them contain at least one massive YSO; these clumps also contain extended green objects and Class II methanol masers. This morphology suggests that the formation of a second generation of massive stars has been triggered by the central bipolar HII region. It occurs in the compressed material of the parental cloud.

We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.

In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.

We report the first detection of the J = 1 - 0 (102.6 GHz) rotational lines of CF+ (fluoromethylidynium ion) towards CygX-N63, a young and massive protostar of the Cygnus X region. This detection occurred as part of an unbiased spectral survey of this object in the 0.8-3 mm range, performed with the IRAM 30m telescope. The data were analyzed using a local thermodynamical equilibrium model (LTE model) and a population diagram in order to derive the column density. The line velocity (-4 km s-1) and line width (1.6 km s-1) indicate an origin from the collapsing envelope of the protostar. We obtain a CF+ column density of 4.10e11 cm-2. The CF+ ion is thought to be a good tracer for C+ and assuming a ratio of 10e-6 for CF+/C+, we derive a total number of C+ of 1.2x10e53 within the beam. There is no evidence of carbon ionization caused by an exterior source of UV photons suggesting that the protostar itself is the source of ionization. Ionization from the protostellar photosphere is not efficient enough. In contrast, X-ray ionization from the accretion shock(s) and UV ionization from outflow shocks could provide a large enough ionizing power to explain our CF+ detection. Surprisingly, CF+ has been detected towards a cold, massive protostar with no sign of an external photon dissociation region (PDR), which means that the only possibility is the existence of a significant inner source of C+. This is an important result that opens interesting perspectives to study the early development of ionized regions and to approach the issue of the evolution of the inner regions of collapsing envelopes of massive protostars. The existence of high energy radiations early in the evolution of massive protostars also has important implications for chemical evolution of dense collapsing gas and could trigger peculiar chemistry and early formation of a hot core.

Recent ALMA observations identified one of the most massive star-forming cores yet observed in the Milky Way; SDC335-MM1, within the infrared dark cloud SDC335.579-0.292. Along with an accompanying core MM2, SDC335 appears to be in the early stages of its star formation process. In this paper we aim to constrain the properties of the stars forming within these two massive millimetre sources. Observations of SDC335 at 6, 8, 23 and 25GHz were made with the ATCA. We report the results of these continuum measurements, which combined with archival data, allow us to build and analyse the spectral energy distributions (SEDs) of the compact sources in SDC335. Three HCHII regions within SDC335 are identified, two within the MM1 core. For each HCHII region, a free-free emission curve is fit to the data allowing the derivation of the sources' emission measure, ionising photon flux and electron density. Using these physical properties we assign each HCHII region a ZAMS spectral type, finding two protostars with characteristics of spectral type B1.5 and one with a lower limit of B1-B1.5. Ancillary data from infrared to mm wavelength are used to construct free-free component subtracted SEDs for the mm-cores, allowing calculation of the bolometric luminosities and revision of the previous gas mass estimates. The measured luminosities for the two mm-cores are lower than expected from accreting sources displaying characteristics of the ZAMS spectral type assigned to them. The protostars are still actively accreting, suggesting that a mechanism is limiting the accretion luminosity, we present the case for two different mechanisms capable of causing this. Finally, using the ZAMS mass values as lower limit constraints, a final stellar population for SDC335 was synthesised finding SDC335 is likely to be in the process of forming a stellar cluster comparable to the Trapezium Cluster and NGC6334 I(N).

We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.

We present CO, H2, HI and HISA distributions from a set of simulations of grand design spirals including stellar feedback, self-gravity, heating and cooling. We replicate the emission of the 2nd Galactic Quadrant by placing the observer inside the modelled galaxies and post process the simulations using a radiative transfer code, so as to create synthetic observations. We compare the synthetic datacubes to observations of the 2nd Quadrant of the Milky Way to test the ability of the current models to reproduce the basic chemistry of the Galactic ISM, as well as to test how sensitive such galaxy models are to different recipes of chemistry and/or feedback. We find that models which include feedback and self-gravity can reproduce the production of CO with respect to H2 as observed in our Galaxy, as well as the distribution of the material perpendicular to the Galactic plane. While changes in the chemistry/feedback recipes do not have a huge impact on the statistical properties of the chemistry in the simulated galaxies, we find that the inclusion of both feedback and self-gravity are crucial ingredients, as our test without feedback failed to reproduce all of the observables. Finally, even though the transition from H2 to CO seems to be robust, we find that all models seem to underproduce molecular gas, and have a lower molecular to atomic gas fraction than is observed. Nevertheless, our fiducial model with feedback and self-gravity has shown to be robust in reproducing the statistical properties of the basic molecular gas components of the ISM in our Galaxy.

We investigate cloud-cloud collisions, and GMC evolution, in hydrodynamic simulations of isolated galaxies. The simulations include heating and cooling of the ISM, self--gravity and stellar feedback. Over timescales $<5$ Myr most clouds undergo no change, and mergers and splits are found to be typically two body processes, but evolution over longer timescales is more complex and involves a greater fraction of intercloud material. We find that mergers, or collisions, occur every 8-10 Myr (1/15th of an orbit) in a simulation with spiral arms, and once every 28 Myr (1/5th of an orbit) with no imposed spiral arms. Both figures are higher than expected from analytic estimates, as clouds are not uniformly distributed in the galaxy. Thus clouds can be expected to undergo between zero and a few collisions over their lifetime. We present specific examples of cloud--cloud interactions in our results, including synthetic CO maps. We would expect cloud--cloud interactions to be observable, but find they appear to have little or no impact on the ISM. Due to a combination of the clouds' typical geometries, and moderate velocity dispersions, cloud--cloud interactions often better resemble a smaller cloud nudging a larger cloud. Our findings are consistent with the view that spiral arms make little difference to overall star formation rates in galaxies, and we see no evidence that collisions likely produce massive clusters. However, to confirm the outcome of such massive cloud collisions we ideally need higher resolution simulations.

To date, most numerical simulations of molecular clouds, and star formation within them, assume a uniform density sphere or box with an imposed turbulent velocity field. In this work, we select molecular clouds from galactic scale simulations as initial conditions, increase their resolution, and re-simulate them using the SPH code Gadget2. Our approach provides clouds with morphologies, internal structures, and kinematics that constitute more consistent and realistic initial conditions for simulations of star formation. We perform comparisons between molecular clouds derived from a galactic simulation, and spheres of turbulent gas of similar dimensions, mass and velocity dispersion. We focus on properties of the clouds such as their density, velocity structure and star formation rate. We find that the inherited velocity structure of the galactic clouds has a significant impact on the star formation rate and evolution of the cloud. Our results indicate that, although we can follow the time evolution of star formation in any simulated cloud, capturing the entire history is difficult as we ignore any star formation that might have occurred before initialisation. Overall, the turbulent spheres do not match the complexity of the galactic clouds.

We used PdBI observations of SiO (2-1) to investigate the morphology and profile of the SiO emission within several massive dense clumps (MDCs) in Cygnus-X. We find that most molecular outflows are detected in both SiO and CO, although there are some cases of CO outflows with no SiO counterpart. We find a significant amount of narrow line SiO emission that appears to be unrelated to outflows. The fraction of the total SiO luminosity that is not associated with outflows is highly variable in the different MDCs (from 10% to 90%); this might be a problem when extrapolating outflow properties from SiO luminosities without resolving individual outflows. The extent of the narrow SiO emission varies from rather compact (~ 0.03 pc) to widespread (~0.2 pc), and its kinematics often differs from those found by other high-density tracers such as H13CO+. We find that the least centrally concentrated clumps with the least massive protostellar cores have the most widespread narrow SiO emission. In line with previous evidence of SiO emission associated with low-velocity shocks, we propose an evolutionary picture to explain the existence and distribution of narrow SiO line profiles. In this scenario, the least centrally condensed MDCs are at an early stage where the SiO emission traces shocks from the large-scale collapse of material onto the MDC (e.g. CygX-N40). As the MDC collapses, the SiO emission becomes more confined to the close surroundings of cores, tracing the post-shock material from the infalling MDC against the dense cores (e.g. CygX-N3, N12, and N48). At later stages, when single massive protostars are formed, the SiO luminosity is largely dominated by powerful outflows, and the weaker narrow component shows perhaps the last remnants of the initial collapse (e.g. CygX-N53 and N63).

Context: Star formation efficiency (SFE) theories are currently based on statistical distributions of turbulent cloud structures and a simple model of star formation from cores. They remain poorly tested, especially at the highest densities. Aims: We investigate the effects of gas density on the SFE through measurements of the core formation efficiency (CFE). With a total mass of $\sim2\times10^4$ M$_\odot$, the W43-MM1 ridge is one of the most convincing candidate precursor of starburst clusters and thus one of the best place to investigate star formation. Methods: We used high-angular resolution maps obtained at 3 mm and 1 mm within W43-MM1 with the IRAM Plateau de Bure Interferometer to reveal a cluster of 11 massive dense cores (MDCs), and, one of the most massive protostellar cores known. An Herschel column density image provided the mass distribution of the cloud gas. We then measured the 'instantaneous' CFE and estimated the SFE and the star formation rate (SFR) within subregions of the W43-MM1 ridge. Results: The high SFE found in the ridge ($\sim$6% enclosed in $\sim$8 pc$^3$) confirms its ability to form a starburst cluster. There is however a clear lack of dense cores in the northern part of the ridge, which may be currently assembling. The CFE and the SFE are observed to increase with volume gas density while the SFR steeply decreases with the virial parameter, $\alpha_{vir}$. Statistical models of the SFR may well describe the outskirts of the W43-MM1 ridge but struggle to reproduce its inner part, which corresponds to measurements at low $\alpha_{vir}$. It may be that ridges do not follow the log-normal density distribution, Larson relations, and stationary conditions forced in the statistical SFR models.

The Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey aims to characterise the physical and chemical evolution of high-mass star-forming clumps. Exploiting the unique broad frequency range and on-the-fly mapping capabilities of the Australia Telescope National Facility Mopra 22 m single-dish telescope, MALT90 has obtained 3' x 3' maps toward ~2000 dense molecular clumps identified in the ATLASGAL 870 um Galactic plane survey. The clumps were selected to host the early stages of high-mass star formation and to span the complete range in their evolutionary states (from prestellar, to protostellar, and on to HII regions and photodissociation regions). Because MALT90 mapped 16 lines simultaneously with excellent spatial (38") and spectral (0.11 km/s) resolution, the data reveal a wealth of information about the clump's morphologies, chemistry, and kinematics. In this paper we outline the survey strategy, observing mode, data reduction procedure, and highlight some early science results. All MALT90 raw and processed data products are available to the community. With its unprecedented large sample of clumps, MALT90 is the largest survey of its type ever conducted and an excellent resource for identifying interesting candidates for high resolution studies with ALMA.

As natural consequences of the accretion process, outflows are one of the few (indirect) tracers of accretion. We used CO(2-1) PdBI observations towards 6 MDCs in Cygnus-X, containing 9 high-mass cores, to investigate what the accretion process and origin of the material feeding the precursors of high-mass stars are. We compared our sample to low-mass objects from the literature and developed simple evolutionary models to reproduce the observables. We find that 8/9 high-mass cores drive clear individual outflows as true equivalents of Class 0 protostars in the high-mass regime. The remaining core has only a tentative outflow detection. It could be amongst the first examples of a true individual high-mass prestellar core. We find that the momentum flux of high-mass objects has a linear relation to the envelope mass, as a scale-up of the relations found for low-mass protostars. This suggests a fundamental proportionality between accretion rates and mass reservoir, suggesting identical collapse timescales for all masses. We conclude that if the pre-collapse evolution is quasi-static, the fragmentation scale (which is similar for all masses) would limit the size of the initial mass reservoirs, leading to shorter free-fall times for higher mass stars. However, as we find identical collapse timescales for all masses, a significant turbulent/magnetic support is needed to slow down the collapse of the more massive envelopes in a quasi-static view. With this support still to be discovered, and with indications of large dynamics in pre-collapse gas for high-mass star formation, we propose that such an identical collapse timescale implies that the initial densities, which should set the duration of the collapse, should be similar for all masses. This suggests that the mass that incorporates massive stars has to have been accreted in a dynamical way from larger scales than those of low-mass stars.

The formation of high-mass stars is tightly linked to that of their parental clouds. We here focus on the high-density parts of W43, a molecular cloud undergoing an efficient event of formation. The cloud structure is studied with a column density image derived from Herschel continuum maps obtained at 70, 160, 250, 350, and 500 micron. We identify two high-column density filamentary clouds, quoted as the W43-MM1 and W43-MM2 ridges, which both account for 1.5x10^4 Msun gas mass above 10^23 cm-2 and within areas of 5 and 14pc^2, respectively. We used the N_2H^+ 1--0 line to confirm that the W43-MM1 and W43-MM2 ridges are structures coherent in velocity and gravitationally bound, despite their large velocity dispersion and ~5 kms line widths. The most intriguing result of the W43 large program is the bright wide-spread SiO 2--1 emission: 1--11 K kms$ stretching an area of ~28 pc^2. Concentrated toward the W43-MM1 and W43-MM2 ridges and their immediate surroundings, it leads to a total luminosity of L_SiO 2-1 ~4 10^4 K kms kpc^2pc^2. We measured a steep relation between the luminosity and velocity extent of the SiO~2--1 lines and propose to use it to distinguish the low-velocity shocks observed here from the more classical high-velocity ones associated with outflows of high-mass young stellar objects. We used state-of-the-art shock models to demonstrate that low-velocity (<10 kms^-1) shocks with a small amount (10%) of Si atoms initially in gas phase or in grain mantles can explain the observed SiO column density in W43. The spatial and velocity overlaps between the ridges high-density gas (n_H2>10^4-10^5 cm^-3) and the shocked SiO gas suggests that ridges could be forming via colliding flows driven by gravity and accompanied by low-velocity shocks. This mechanism may be the initial conditions for the formation of young massive clusters in these ridges.

We characterize the Millimeter Astronomy Legacy Team 90 GHz (MALT90) Survey and the Mopra telescope at 90 GHz. We combine repeated position-switched observations of the source G300.968+01.145 with a map of the same source in order to estimate the pointing reliability of the position-switched observations and, by extension, the MALT90 survey; we estimate our pointing uncertainty to be 8 arcseconds. We model the two strongest sources of systematic gain variability as functions of elevation and time-of-day and quantify the remaining absolute flux uncertainty. Corrections based on these two variables reduce the scatter in repeated observations from 12-25% down to 10-17%. We find no evidence for intrinsic source variability in G300.968+01.145. For certain applications, the corrections described herein will be integral for improving the absolute flux calibration of MALT90 maps and other observations using the Mopra telescope at 90 GHz.

The relative importance of primordial molecular cloud fragmentation versus large-scale accretion still remains to be assessed in the context of massive core/star formation. Studying the kinematics of the dense gas surrounding massive-star progenitors can tell us the extent to which large-scale flow of material impacts the growth in mass of star-forming cores. Here we present a comprehensive dataset of the 5500(+/-800) Msun infrared dark cloud SDC335.579-0.272 (hereafter SDC335) which exhibits a network of cold, dense, parsec-long filaments. Atacama Large Millimeter Array (ALMA) Cycle 0 observations reveal two massive star-forming cores, MM1 and MM2, sitting at the centre of SDC335 where the filaments intersect. With a gas mass of 545(+770,-385) Msun contained within a source diameter of 0.05pc, MM1 is one of the most massive, compact protostellar cores ever observed in the Galaxy. As a whole, SDC335 could potentially form an OB cluster similar to the Trapezium cluster in Orion. ALMA and Mopra single-dish observations of the SDC335 dense gas furthermore reveal that the kinematics of this hub-filament system are consistent with a global collapse of the cloud. These molecular-line data point towards an infall velocity V_{inf} =0.7(+/-0.2) km/s, and a total mass infall rate \dot{M}_{inf} = 2.5(+/-1.0) x 10^{-3} Msun/yr towards the central pc-size region of SDC335. This infall rate brings 750(+/-300) Msun of gas to the centre of the cloud per free-fall time (t_{ff}=3x10^5 yr). This is enough to double the mass already present in the central pc-size region in 3.5(+2.2,-1.0) x t_{ff}. These values suggest that the global collapse of SDC335 over the past million year resulted in the formation of an early O-type star progenitor at the centre of the cloud's gravitational potential well.

We present Herschel observations from the Herschel Gould Belt Survey and SCUBA-2 science verification observations from the JCMT Gould Belt Survey of the B1 clump in the Perseus molecular cloud. We determined the dust emissivity index using four different techniques to combine the Herschel PACS+SPIRE data at 160 - 500 microns with the SCUBA-2 data at 450 microns and 850 microns. Of our four techniques, we found the most robust method was to filter-out the large-scale emission in the Herschel bands to match the spatial scales recovered by the SCUBA-2 reduction pipeline. Using this method, we find beta ~ 2 towards the filament region and moderately dense material and lower beta values (beta > 1.6) towards the dense protostellar cores, possibly due to dust grain growth. We find that beta and temperature are more robust with the inclusion of the SCUBA-2 data, improving estimates from Herschel data alone by factors of ~ 2 for beta and by ~ 40% for temperature. Furthermore, we find core mass differences of < 30% compared to Herschel-only estimates with an adopted beta = 2, highlighting the necessity of long wavelength submillimeter data for deriving accurate masses of prestellar and protostellar cores.

The Pipe Nebula is a molecular cloud hosting the B59 region as its only active star-forming clump. While the particular importance of outflows in active star forming regions is subject of debate, the quiet nature of the gas in B59 makes it a good site to directly see the impact of protostellar feedback on the quiescent dense gas. Using HARP at the JCMT, we mapped the B59 region with the J=3-2 transition of 12CO to study the kinematics and energetics of the outflows, and 13CO and C18O to study the overall dynamics of the ambient cloud, the physical properties of the gas, and the hierarchical structure of the region. The B59 region has a total of 30Msun of cold and quiescent material, mostly gravitationally bound, with narrow line widths throughout. Such low levels of turbulence in non-star-forming sites of B59 are indicative of the intrinsic initial conditions of the cloud. On the other hand, close to the forming protostars the impact of the outflows is observed as a localised increase of both line widths from 0.3 to 1 km/s, and 13CO excitation temperatures by 2-3K. The impact of the outflows is also evident in the low column density material which shows signs of being pushed, shaped and carved by the outflow bow shocks as they pierce their way out of the cloud. Much of this structure is readily apparent in a dendrogram analysis of the cloud. The low mass of B59 together with its intrinsically quiescent gas and small number of protostars, allows the identification of specific regions where the outflows from the embedded sources interact the dense gas. Our study suggests that outflows are an important mechanism for injecting and sustaining supersonic turbulence at sub-parsec size scales. We find that less than half of the outflow energy is deposited as turbulent energy of the gas, however this turbulent energy is sufficient to slow down the collapse of the region.

The Gould Belt Legacy Survey will map star-forming regions within 500 pc, using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common-User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes HARP observations of the J = 3-2 transitions of 13CO and C18O towards Orion A. The 1500-resolution observations cover 5 pc of the Orion filament, including OMC1 (inc. BN-KL and Orion Bar), OMC 2/3 and OMC 4, and allow a comparative study of the molecular gas properties throughout the star-forming cloud. The filament shows a velocity gradient of ~1 km/s /pc between OMC 1, 2 and 3, and high velocity emission is detected in both isotopologues. The Orion Nebula and Bar have the largest masses and line widths, and dominate the mass and energetics of the high velocity material. Compact, spatially resolved emission from CH3CN, 13CH3OH, SO, HCOOCH3, C2H5OH, CH3CHO and CH3OCHO is detected towards the Orion Hot Core. The cloud is warm, with a median excitation temperature of ~24 K; the Orion Bar has the highest excitation temperature gas, at >80 K. The C18O excitation temperature correlates well with the dust temperature (to within 40%). The C18O emission is optically thin, and the 13CO emission is marginally optically thick; despite its high mass, OMC 1 shows the lowest opacities. A virial analysis indicates that Orion A is too massive for thermal or turbulent support, but is consistent with a model of a filamentary cloud that is threaded by helical magnetic fields. The variation of physical conditions across the cloud is reflected in the physical characteristics of the dust cores....continued

The complexity of the ISM is such that it is unlikely that star formation is initiated in the same way in all molecular clouds. While some clouds seem to collapse on their own, others may be triggered by an external event such as a cloud/flow collision forming a gravitationally unstable enhanced density layer. This work tests cloud-cloud collisions as the triggering mechanism for star formation in the Serpens Main Cluster as has been suggested by previous work. A set of smoothed particle hydrodynamics (SPH) simulations of the collision between two cylindrical clouds are performed and compared to (sub)millimetre observations of the Serpens Main Cluster. A configuration has been found which reproduces many of the observed characteristics of Serpens, including some of the main features of the peculiar velocity field. The evolution of the velocity with position throughout the model is similar to that observed and the column density and masses within the modeled cloud agree with those measured for the SE sub-cluster. Furthermore, our results also show that an asymmetric collision provides the ingredients to reproduce lower density filaments perpendicular to the main structure, similar to those observed. In this scenario, the formation of the NW sub-cluster of Serpens can be reproduced only if there is a pre-existing marginally gravitationally unstable region at the time the collision occurs. This work supports the interpretation that a collision between two clouds may have been the trigger of the most recent burst of star formation in Serpens. It not only explains the complicated velocity structure seen in the region, but also the temperature differences between the north (in "isolated" collapse) and the south (resulting from the shock between the clouds). In addition it provides an explanation for the sources in the south having a larger spread in age than those in the north.

The Gould Belt Legacy Survey (GBS) on the JCMT has observed a region of 260 square arcminutes in 12CO J=3--2 emission, and a 190 square arcminute subset of this in 13CO and C18O towards the Serpens molecular cloud. We examine the global velocity structure of the non-outflowing gas, and calculate excitation temperatures and opacities. The large scale mass and energetics of the region are evaluated, with special consideration for high velocity gas. We find the cloud to have a mass of 203 solar masses, and to be gravitationally bound, and that the kinetic energy of the outflowing gas is approximately seventy percent of the turbulent kinetic energy of the cloud. We identify compact outflows towards some of the submillimetre Class 0/I sources in the region

The Serpens North Cluster is a nearby low mass star forming region which is part of the Gould Belt. It contains a range of young stars thought to correspond to two different bursts of star formation and provides the opportunity to study different stages of cluster formation. This work aims to study the molecular gas in the Serpens North Cluster to probe the origin of the most recent burst of star formation in Serpens. Transitions of the C17O and C18O observed with the IRAM 30m telescope and JCMT are used to study the mass and velocity structure of the region while the physical properties of the gas are derived using LTE and non-LTE analyses of the three lowest transitions of C18O. The molecular emission traces the two centres of star formation which are seen in submillimetre dust continuum emission. In the ~40M_sun NW sub-cluster the gas and dust emission trace the same structures although there is evidence of some depletion of the gas phase C18O. The gas has a very uniform temperature (~10K) and velocity (~8.5km/s) throughout the region. This is in marked contrast to the SE sub-cluster. In this region the dust and the gas trace different features, with the temperature peaking between the submillimetre continuum sources, reaching up to ~14K. The gas in this region has double peaked line profiles which reveal the presence of a second cloud in the line of sight. The submillimetre dust continuum sources predominantly appear located in the interface region between the two clouds. Even though they are at a similar stage of evolution, the two Serpens sub-clusters have very different characteristics. We propose that these differences are linked to the initial trigger of the collapse in the regions and suggest that a cloud-cloud collision could explain the observed properties.

As part of a JCMT Legacy Survey of star formation in the Gould Belt, we present early science results for Taurus. CO J=3-2 maps have been secured along the north-west ridge and bowl, collectively known as L 1495, along with deep 13CO and C18O J=3-2 maps in two sub-regions. With these data we search for molecular outflows, and use the distribution of flows, HH objects and shocked H2 line emission features, together with the population of young stars, protostellar cores and starless condensations to map star formation across this extensive region. In total 21 outflows are identified. It is clear that the bowl is more evolved than the ridge, harbouring a greater population of T Tauri stars and a more diffuse, more turbulent ambient medium. By comparison, the ridge contains a much younger, less widely distributed population of protostars which, in turn, is associated with a greater number of molecular outflows. We estimate the ratio of the numbers of prestellar to protostellar cores in L 1495 to be ~ 1.3-2.3, and of gravitationally unbound starless cores to (gravitationally bound) prestellar cores to be ~ 1. If we take previous estimates of the protostellar lifetime of ~ 5 x 10^5 yrs, this indicates a prestellar lifetime of 9(+/-3) x 10^5 yrs. From the number of outflows we also crudely estimate the star formation efficiency in L 1495, finding it to be compatible with a canonical value of 10-15 %. We note that molecular outflow-driving sources have redder near-IR colours than their HH jet-driving counterparts. We also find that the smaller, denser cores are associated with the more massive outflows, as one might expect if mass build-up in the flow increases with the collapse and contraction of the protostellar envelope.

The Gould Belt Legacy Survey will survey nearby star-forming regions (within 500 pc), using HARP (Heterodyne Array Receiver Programme), SCUBA-2 (Submillimetre Common- User Bolometer Array 2) and POL-2 (Polarimeter 2) on the James Clerk Maxwell Telescope (JCMT). This paper describes the initial data obtained using HARP to observe 12CO, 13CO and C18O J = 3 - 2 towards two regions in Orion B, NGC 2024 and NGC 2071. We describe the physical characteristics of the two clouds, calculating temperatures and opacities utilizing all three isotopologues. We find good agreement between temperatures calculated from CO and from dust emission in the dense, energetic regions. We determine the mass and energetics of the clouds, and of the high-velocity material seen in 12CO emission, and compare the relative energetics of the high- and low-velocity material in the two clouds. We present a CLUMPFIND analysis of the 13CO condensations. The slope of the condensation mass functions, at the high-mass ends, is similar to the slope of the initial mass function.

The first images of 6.7-GHz methanol masers in the massive star-forming regions DR21(OH) and DR21(OH)N are presented. By measuring the shapes, radial velocities and polarization properties of these masers it is possible to map out the structure, kinematics and magnetic fields in the molecular gas that surrounds newly-formed massive stars. The intrinsic angular resolution of the observations was 43 mas (~100 AU at the distance of DR21), but structures far smaller than this were revealed by employing a non-standard mapping technique. This technique was used in an attempt to identify the physical structure (e.g. disc, outflow, shock) associated with the methanol masers. Two distinct star-forming centres were identified. In DR21(OH) the masers had a linear morphology, and the individual maser spots each displayed an internal velocity gradient in the same direction as the large-scale structure. They were detected at the same position as the OH 1.7-GHz ground-state masers, close to the centre of an outflow traced by CO and class I methanol masers. The shape and velocity gradients of the masers suggests that they probably delineate a shock. In DR21(OH)N the methanol masers trace an arc with a double-peaked profile and a complex velocity gradient. This velocity gradient closely resembles that of a Keplerian disc. The masers in the arc are 4.5% linearly polarized, with a polarization angle that indicates that the magnetic field direction is roughly perpendicular to the large-scale magnetic field in the region (indicated by lower angular resolution measurements of the CO and dust polarization). The suitability of channel-by-channel centroid mapping is discussed as an improved and viable means to maximise the information gained from the data.