We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

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.

CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines.

Massive stars and their associated ionized (HII) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H regions are still poorly known. To investigate the impact of HII regions on the formation of filaments, we imaged the Galactic HII region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied. We used the ArT\'eMiS camera on the APEX telescope and combined the ArT\'eMiS data at 350 and 450 microns with Herschel-SPIRE/HOBYS. We studied the dense gas distribution around RCW 120 with a resolution of 8 arcsec (0.05 pc at a distance of 1.34 kpc). Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the HII region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer like profile with a deconvolved Gaussian FWHM of 0.09 pc, as observed for filaments in low-mass star-forming regions. This study suggests that compression exerted by HII regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArT\'eMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure

The recent data collected by {\it Herschel} have confirmed that interstellar structures with filamentary shape are ubiquitously present in the Milky Way. Filaments are thought to be formed by several physical mechanisms acting from the large Galactic scales down to the sub-pc fractions of molecular clouds, and they might represent a possible link between star formation and the large-scale structure of the Galaxy. In order to study this potential link, a statistically significant sample of filaments spread throughout the Galaxy is required. In this work we present the first catalogue of $32,059$ candidate filaments automatically identified in the Hi-GAL survey of the entire Galactic Plane. For these objects we determined morphological (length, $l^{a}$, and geometrical shape) and physical (average column density, $N_{\rm H_{2}}$, and average temperature, $T$) properties. We identified filaments with a wide range of properties: 2$'$\,$\leq l^{a}\leq$\, 100$'$, $10^{20} \leq N_{\rm H_{2}} \leq 10^{23}$\,cm$^{-2}$ and $10 \leq T\leq$ 35\,K. We discuss their association with the Hi-GAL compact sources, finding that the most tenuous (and stable) structures do not host any major condensation and we also assign a distance to $\sim 18,400$ filaments for which we determine mass, physical size, stability conditions and Galactic distribution. When compared to the spiral arms structure, we find no significant difference between the physical properties of on-arm and inter-arm filaments. We compared our sample with previous studies, finding that our Hi-GAL filament catalogue represents a significant extension in terms of Galactic coverage and sensitivity. This catalogue represents an unique and important tool for future studies devoted to understanding the filament life-cycle.

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.

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.

Understanding where and when the mass of stars is determined is one of the fundamental, mostly unsolved, questions in astronomy. Here, we present the first results of GASTON, the Galactic Star Formation with NIKA2 large programme on the IRAM 30m telescope, that aims to identify new populations of low-brightness sources to tackle the question of stellar mass determination across all masses. In this paper, we focus on the high-mass star formation part of the project, for which we map a $\sim2$ deg$^2$ region of the Galactic plane around $l=24^\circ$ in both 1.2 mm and 2.0 mm continuum. Half-way through the project, we reach a sensitivity of 3.7 mJy/beam at 1.2mm. Even though larger than our target sensitivity of 2 mJy, the current sensitivity already allows the identification of a new population of cold, compact sources that remained undetected in any (sub-)mm Galactic plane survey so far. In fact, about 25% of the $\sim 1600$ compact sources identified in the 1.2 mm GASTON image are new detections. We present a preliminary analysis of the physical properties of the GASTON sources as a function of their evolutionary stage, arguing for a potential evolution of the mass distribution of these sources with time.

(Abridged) The initial physical conditions of high-mass stars and protoclusters remain poorly characterized. To this end we present the first targeted ALMA 1.3mm continuum and spectral line survey towards high-mass starless clump candidates, selecting a sample of 12 of the most massive candidates ($400-4000\, M_\odot$) within 5 kpc. The joint 12+7m array maps have a high spatial resolution of $\sim 3000\, \mathrm{au}$ ($\sim 0.8^{\prime\prime}$) and have point source mass-completeness down to $\sim 0.3\, M_\odot$ at $6\sigma$ (or $1\sigma$ column density sensitivity of $1.1\times10^{22}\, \mathrm{cm^{-2}}$). We discover previously undetected signposts of low-luminosity star formation from CO (2-1) and SiO (5-4) bipolar outflows and other signatures towards 11 out of 12 clumps, showing that current MIR/FIR Galactic Plane surveys are incomplete to low- and intermediate-mass protostars ($\lesssim 50\, L_\odot$). We compare a subset of the observed cores with a suite of radiative transfer models of starless cores. We find a high-mass starless core candidate with a model-derived mass consistent with $29^{52}_{15}\, M_\odot$ when integrated over size scales of $2\times10^4\, \mathrm{au}$. Unresolved cores are poorly fit by starless core models, supporting the interpretation that they are protostellar even without detection of outflows. Substantial fragmentation is observed towards 10 out of 12 clumps. We extract sources from the maps using a dendrogram to study the characteristic fragmentation length scale. Nearest neighbor separations when corrected for projection are consistent with being equal to the clump average thermal Jeans length. Our findings support a hierarchical fragmentation process, where the highest density regions are not strongly supported against thermal gravitational fragmentation by turbulence or magnetic fields.

Feedback from high-mass stars shapes the ISM of galaxies and thereby impacts gas that will form future generations of stars. However, due to our inability to track the time evolution of individual molecular clouds, quantifying the exact role of feedback on their star formation history is an observationally challenging task. In the present study, we take advantage of the unique properties of the G316.75-00.00 high-mass star-forming ridge to determine how feedback from O-stars impacts the dynamical stability of filaments. G316.75 is a 13.6pc long ridge containing 18,900Msun of H2 gas which is half IR dark, half IR bright. The IR bright half has already formed 4 O-type stars over the past 2Myr, while the IR dark half is still quiescent. Therefore, by assuming the gas properties of the dark part represents an earlier evolutionary stage of the bright part, we can quantify how feedback impacts these properties by contrasting them. We use archive Herschel and molecular line data, tracing both dense (NH3 and N2H+) and more diffuse (13CO) gas, to measure the ratio of kinetic to gravitational energy per-unit-length, virial-ratio-per-line, across the ridge for a range of gas volume densities. We show that despite the presence of 4 embedded O-stars, feedback cannot unbind the ridge's mass except for some small gas pockets near the O-stars. In fact, the virial-ratio-per-line is almost indistinguishable for both parts of the ridge. These results are at odds with most simulations where O-star-forming clouds are dispersed by feedback within a few cloud free-fall time. We conclude that the discrepancy between such simulations and our observations originates from different cloud morphologies and average densities when the first O-star forms. These results have important implications regarding, for instance, how feedback is implemented within cosmological and galactic-scale simulations.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

(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 the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

Filaments are a key step on the path that leads from molecular clouds to star formation. However, their characteristics are heavily debated, and the exact processes that lead to their formation and fragmentation into dense cores still remain to be fully understood. We aim at characterising the mass, kinematics, and stability against gravitational collapse of a statistically significant sample of filaments in the Orion B molecular cloud, which is renown for its very low star formation efficiency. We characterise the gas column densities and kinematics over a field of 1.9 deg$^2$, using C$^{18}$O(J=1-0) data from the IRAM-30m large programme ORION-B. Using two different Hessian-based filters, we extract and compare two filamentary networks, each containing over 100 filaments. Independent of the extraction method, the filaments have widths of 0.12$\pm$0.04 pc, and show a wide range of linear (1 - 100 $M_{\odot}$pc$^{-1}$) and volume densities (2.10$^3$ - 2.10$^5$ cm$^{-3}$). Compared to previous studies, the filament population is dominated by low-density, thermally sub-critical structures, suggesting that most of the identified filaments are not collapsing to form stars. In fact, only ~1% of the Orion B cloud mass covered by our observations can be found in super-critical, star-forming filaments, explaining the low star formation efficiency of the region. The velocity profiles observed across the filaments show quiescence in the centre, and coherency in the plane of the sky, despite being mostly supersonic. The filaments in Orion B apparently belong to a continuum which contains a few elements comparable to already studied star-forming filaments as well as many lower-density, gravitationally unbound structures. This comprehensive study of the Orion B filaments shows that the mass fraction in super-critical filaments is a key factor in determining star formation efficiency.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

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.

The thermal emission of dust grains is a powerful tool for probing cold, dense regions of molecular gas in the ISM, and so constraining dust properties is key to obtaining accurate measurements of dust mass and temperature. By placing constraints on the dust emissivity spectral index, beta, towards two star-forming infrared dark clouds, SDC18.888 and SDC24.489, we evaluate the role of mass concentration in the associated star-formation activity. We exploit the simultaneous 1.2mm and 2.0mm imaging capability of NIKA on the IRAM 30m telescope to construct maps of beta for both clouds, and by incorporating Herschel observations, we create H2 column density maps with 13" resolution. While we find no significant systematic radial variations around the most massive clumps in either cloud on >0.1 pc scales, their mean beta values are significantly different, with beta = 2.07 +/- 0.09 (rand) +/- 0.25 (syst) for SDC18.888 and beta = 1.71 +/- 0.09 (rand) +/- 0.25 (syst) for SDC24.489. These differences could be a consequence of the very different environments in which both clouds lie, and we suggest that the proximity of SDC18.888 to the W39 HII region may raise beta on scales of 1 pc. We also find that the mass in SDC24.489 is more centrally concentrated and circularly symmetric than in SDC18.888, and is consistent with a scenario in which spherical globally-collapsing clouds concentrate a higher fraction of their mass into a single core than elongated clouds that will more easily fragment, distributing their mass into many cores. We demonstrate that beta variations towards interstellar clouds can be robustly constrained with high-SNR NIKA observations, providing more accurate estimates of their masses. The methods presented here will be applied to the Galactic Star Formation with NIKA2 (GASTON) large programme, extending our analysis to a statistically significant sample of star-forming clouds.

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]

We present far-infrared observations of Monoceros R2 (a giant molecular cloud at approximately 830 pc distance, containing several sites of active star formation), as observed at 70 {\mu}m, 160 {\mu}m, 250 {\mu}m, 350 {\mu}m, and 500 {\mu}m by the Photodetector Array Camera and Spectrometer (PACS) and Spectral and Photometric Imaging Receiver (SPIRE) instruments on the Herschel Space Observatory as part of the Herschel imaging survey of OB young stellar objects (HOBYS) Key programme. The Herschel data are complemented by SCUBA-2 data in the submillimetre range, and WISE and Spitzer data in the mid-infrared. In addition, C18O data from the IRAM 30-m Telescope are presented, and used for kinematic information. Sources were extracted from the maps with getsources, and from the fluxes measured, spectral energy distributions were constructed, allowing measurements of source mass and dust temperature. Of 177 Herschel sources robustly detected in the region (a detection with high signal-to-noise and low axis ratio at multiple wavelengths), including protostars and starless cores, 29 are found in a filamentary hub at the centre of the region (a little over 1% of the observed area). These objects are on average smaller, more massive, and more luminous than those in the surrounding regions (which together suggest that they are at a later stage of evolution), a result that cannot be explained entirely by selection effects. These results suggest a picture in which the hub may have begun star formation at a point significantly earlier than the outer regions, possibly forming as a result of feedback from earlier star formation. Furthermore, the hub may be sustaining its star formation by accreting material from the surrounding filaments.

Previous attempts at segmenting molecular line maps of molecular clouds have focused on using position-position-velocity data cubes of a single line to separate the spatial components of the cloud. In contrast, wide field spectral imaging with large spectral bandwidth in the (sub)mm domain now allows to combine multiple molecular tracers to understand the different physical and chemical phases that constitute giant molecular clouds. We aim at using multiple tracers (sensitive to different physical processes) to segment a molecular cloud into physically/chemically similar regions (rather than spatially connected components). We use a machine learning clustering method (the Meanshift algorithm) to cluster pixels with similar molecular emission, ignoring spatial information. Simple radiative transfer models are used to interpret the astrophysical information uncovered by the clustering. A clustering analysis based only on the J=1-0 lines of 12CO, 13CO and C18O reveals distinct density/column density regimes (nH~100, 500, and >1000 cm-3), closely related to the usual definitions of diffuse, translucent and high-column-density regions. Adding two UV-sensitive tracers, the (1-0) lines of HCO+ and CN, allows us to distinguish two clearly distinct chemical regimes, characteristic of UV-illuminated and UV-shielded gas. The UV-illuminated regime shows overbright HCO+ and CN emission, which we relate to photochemical enrichment. We also find a tail of high CN/HCO+ intensity ratio in UV-illuminated regions. Finer distinctions in density classes (nH~7E3, and 4E4 cm-3) for the densest regions are also identified, likely related to the higher critical density of the CN and HCO+ (1-0) lines. The association of simultaneous multi-line, wide-field mapping and powerful machine learning methods such as the Meanshift algorithm reveals how to decode the complex information available in molecular tracers.

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.

The characterization of dust properties in the interstellar medium (ISM) is key for star formation. Mass estimates are crucial to determine gravitational collapse conditions for the birth of new stellar objects in molecular clouds. However, most of these estimates rely on dust models that need further observational constraints from clouds to prestellar and protostellar cores. We present results of a study of dust emissivity changes based on mm-continuum data obtained with the NIKA camera at the IRAM-30m telescope. Observing dust emission at 1.15 mm and 2 mm allows us to constrain the dust emissivity index ($\beta$) in the Rayleigh-Jeans tail of the dust spectral energy distribution (SED) far from its peak emission, where the contribution of other parameters (i.e. dust temperature) is important. Focusing on the Taurus molecular cloud, a low-mass star-forming regions in the Gould Belt, we analyze the emission properties of several distinct objects in the B213 filament: three prestellar cores, two Class-0/I protostellar cores and one Class-II object. By means of the ratio of the two NIKA channel-maps, we show that in the Rayleigh-Jeans approximation the dust emissivity index varies among the objects. For one prestellar and two protostellar cores, we produce a robust study using Herschel data to constrain the dust temperature of the sources. By using the Abel transform inversion technique we get accurate radial $\beta$ profiles. We find systematic spatial variations of $\beta$ in the protostellar cores that is not observed in the prestellar core. While in the former case $\beta$ decreases toward the center, in the latter it remains constant. Moreover, $\beta$ appears anticorrelated with the dust temperature. We discuss the implication of these results in terms of dust grain evolution between pre- and protostellar cores.

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.

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 first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We present the first data release of the James Clerk Maxwell Telescope (JCMT) Plane Survey (JPS), the JPS Public Release 1 (JPSPR1). JPS is an 850-um continuum survey of six fields in the northern inner Galactic Plane in a longitude range of l=7-63, made with the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2). This first data release consists of emission maps of the six JPS regions with an average pixel-to-pixel noise of 7.19 mJy beam^-1, when smoothed over the beam, and a compact-source catalogue containing 7,813 sources. The 95 per cent completeness limits of the catalogue are estimated at 0.04 Jy beam^-1 and 0.3 Jy for the peak and integrated flux densities, respectively. The emission contained in the compact-source catalogue is 42 +- 5 per cent of the total and, apart from the large-scale (greater than 8') emission, there is excellent correspondence with features in the 500-um Herschel maps. We find that, with two-dimensional matching, 98 +- 2 per cent of sources within the fields centred at l=20, 30, 40 and 50 are associated with molecular clouds, with 91 +- 3 per cent of the l=30 and 40 sources associated with dense molecular clumps. Matching the JPS catalogue to Herschel 70-um sources, we find that 38 +- 1 per cent of sources show evidence of ongoing star formation. The images and catalogue will be a valuable resource for studies of star formation in the Galaxy and the role of environment and spiral arms in the star formation process.

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.

Context. The combination of wideband receivers and spectrometers currently available in (sub-)millimeter observatories deliver wide- field hyperspectral imaging of the interstellar medium. Tens of spectral lines can be observed over degree wide fields in about fifty hours. This wealth of data calls for restating the physical questions about the interstellar medium in statistical terms. Aims. We aim at gaining information on the physical structure of the interstellar medium from a statistical analysis of many lines from different species over a large field of view, without requiring detailed radiative transfer or astrochemical modeling. Methods. We coupled a nonlinear rescaling of the data with one of the simplest multivariate analysis methods, namely the Principal Component Analysis, to decompose the observed signal into components that we interpret first qualitatively and then quantitatively based on our deep knowledge of the observed region and of the astrochemistry at play. Results. We identify 3 principal components, linear compositions of line brightness temperatures, that are correlated at various levels with the column density, the volume density and the UV radiation field. Conclusions. When sampling a sufficiently diverse mixture of physical parameters, it is possible to decompose the molecular emission in order to gain physical insight on the observed interstellar medium. This opens a new avenue for future studies of the interstellar medium.

The nature of turbulence in molecular clouds is one of the key parameters that control star formation efficiency: compressive motions, as opposed to solenoidal motions, can trigger the collapse of cores, or mark the expansion of Hii regions. We try to observationally derive the fractions of momentum density ($\rho v$) contained in the solenoidal and compressive modes of turbulence in the Orion B molecular cloud and relate these fractions to the star formation efficiency in the cloud. The implementation of a statistical method developed by Brunt & Federrath (2014), applied to a $^{13}$CO(J=1-0) datacube obtained with the IRAM-30m telescope, allows us to retrieve 3-dimensional quantities from the projected quantities provided by the observations, yielding an estimate of the compressive versus solenoidal ratio in various regions of the cloud. Despite the Orion B molecular cloud being highly supersonic (mean Mach number $\sim$ 6), the fractions of motion in each mode diverge significantly from equipartition. The cloud's motions are on average mostly solenoidal (excess > 8 % with respect to equipartition), which is consistent with its low star formation rate. On the other hand, the motions around the main star-forming regions (NGC 2023 and NGC 2024) prove to be strongly compressive. We have successfully applied to observational data a method that was so far only tested on simulations, and have shown that there can be a strong intra-cloud variability of the compressive and solenoidal fractions, these fractions being in turn related to the star formation efficiency. This opens a new possibility for star-formation diagnostics in galactic molecular clouds.

Investigations into the substructure of massive star forming regions are essential for understanding the observed relationships between core mass distributions and mass distributions in stellar clusters, differentiating between proposed mechanisms of massive star formation. We study the substructure in the two largest fragments (i.e. cores) MM1 and MM2, in the infrared dark cloud complex SDC13. As MM1 appears to be in a later stage of evolution than MM2, comparing their substructure provides an insight in to the early evolution of massive clumps. We report the results of high resolution SMA dust continuum observations towards MM1 and MM2. Combining these data with Herschel observations, we carry out RADMC-3D radiative transfer modelling to characterise the observed substructure. SMA continuum data indicates 4 sub-fragments in the SDC13 region. The nature of the second brightest sub-fragment (B) is uncertain as it does not appear as prominent at the lower MAMBO resolution or at radio wavelengths. Statistical analysis indicates that it is unlikely to be a background source, an AGB star, or the free-free emission of a HII region. It is plausible that B is a runaway object ejected from MM1. MM1, which is actively forming stars, consists of two sub-fragments A and C. This is confirmed by 70 micron Herschel data. While MM1 and MM2 appear quite similar in previous low resolution observations, at high resolution, the sub-fragment at the centre of MM2 (D) is much fainter than sub-fragment at the centre of MM1 (A). RADMC-3D models of MM1 and MM2 are able to reproduce these results, modelling MM2 with a steeper density profile and higher mass than is required for MM1. The relatively steep density profile of MM2 depends on a significant temperature decrease in its centre, justified by the lack of star formation in MM2. A final stellar population for MM1 was extrapolated, indicating a ...

Herschel observations of nearby molecular clouds suggest that interstellar filaments and prestellar cores represent two fundamental steps in the star formation process. The observations support a picture of low-mass star formation according to which ~ 0.1 pc-wide filaments form first in the cold interstellar medium, probably as a result of large-scale compression of interstellar matter by supersonic turbulent flows, and then prestellar cores arise from gravitational fragmentation of the densest filaments. Whether this scenario also applies to regions of high-mass star formation is an open question, in part because Herschel data cannot resolve the inner width of filaments in the nearest regions of massive star formation. We used the bolometer camera ArTeMiS on the APEX telescope to map the central part of the NGC6334 complex at a factor of > 3 higher resolution than Herschel at 350 microns. Combining ArTeMiS data with Herschel data allowed us to study the structure of the main filament of the complex with a resolution of 8" or < 0.07 pc at d ~ 1.7 kpc. Our study confirms that this filament is a very dense, massive linear structure with a line mass ranging from ~ 500 Msun/pc to ~ 2000 Msun/pc over nearly 10 pc. It also demonstrates that its inner width remains as narrow as W ~ 0.15 +- 0.05 pc all along the filament length, within a factor of < 2 of the characteristic 0.1 pc value found with Herschel for lower-mass filaments in the Gould Belt. While it is not completely clear whether the NGC 6334 filament will form massive stars or not in the future, it is two to three orders of magnitude denser than the majority of filaments observed in Gould Belt clouds, and yet has a very similar inner width. This points to a common physical mechanism for setting the filament width and suggests that some important structural properties of nearby clouds also hold in high-mass star forming regions.

(Abridged) We present the first public release of high-quality data products (DR1) from Hi-GAL, the {\em Herschel} infrared Galactic Plane Survey. Hi-GAL is the keystone of a suite of continuum Galactic Plane surveys from the near-IR to the radio, and covers five wavebands at 70, 160, 250, 350 and 500 micron, encompassing the peak of the spectral energy distribution of cold dust for 8 < T < 50K. This first Hi-GAL data release covers the inner Milky Way in the longitude range 68{\deg} > l > -70{\deg} in a |b|<1{\deg} latitude strip. Photometric maps have been produced with the ROMAGAL pipeline, that optimally capitalizes on the excellent sensitivity and stability of the bolometer arrays of the {\em Herschel} PACS and SPIRE photometric cameras, to deliver images of exquisite quality and dynamical range, absolutely calibrated with {\em Planck} and {\em IRAS}, and recovering extended emission at all wavelengths and all spatial scales. The compact source catalogues have been generated with the CuTEx algorithm, specifically developed to optimize source detection and extraction in the extreme conditions of intense and spatially varying background that are found in the Galactic Plane in the thermal infrared. Hi-GAL DR1 images will be accessible via a dedicated web-based image cutout service. The DR1 Compact Source Catalogues are delivered as single-band photometric lists containing, in addition to source position, peak and integrated flux and source sizes, a variety of parameters useful to assess the quality and reliability of the extracted sources, caveats and hints to help this assessment are provided. Flux completeness limits in all bands are determined from extensive synthetic source experiments and depend on the specific line of sight along the Galactic Plane. Hi-GAL DR1 catalogues contain 123210, 308509, 280685, 160972 and 85460 compact sources in the five bands, respectively.

We present a catalogue of dense cores in a $\sim 4^\circ\times2^\circ$ field of the Taurus star-forming region, inclusive of the L1495 cloud, derived from Herschel SPIRE and PACS observations in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, and 500 $\mu$m continuum bands. Estimates of mean dust temperature and total mass are derived using modified blackbody fits to the spectral energy distributions. We detect 525 starless cores of which $\sim10$-20% are gravitationally bound and therefore presumably prestellar. Our census of unbound objects is $\sim85$% complete for $M>0.015\,M_\odot$ in low density regions ($A_V\stackrel{<}{_\sim}5$ mag), while the bound (prestellar) subset is $\sim85$% complete for $M>0.1\,M_\odot$ overall. The prestellar core mass function (CMF) is consistent with lognormal form, resembling the stellar system initial mass function, as has been reported previously. All of the inferred prestellar cores lie on filamentary structures whose column densities exceed the expected threshold for filamentary collapse, in agreement with previous reports. Unlike the prestellar CMF, the unbound starless CMF is not lognormal, but instead is consistent with a power-law form below $0.3\,M_\odot$ and shows no evidence for a low-mass turnover. It resembles previously reported mass distributions for CO clumps at low masses ($M\stackrel{<}{_\sim}0.3\,M_\odot$). The volume density PDF, however, is accurately lognormal except at high densities. It is consistent with the effects of self-gravity on magnetized supersonic turbulence. The only significant deviation from lognormality is a high-density tail which can be attributed unambiguously to prestellar cores.

Galactic plane surveys of pristine molecular clouds are key for establishing a Galactic-scale view of the earliest stages of star formation. For this reason Peretto & Fuller (2009) built an unbiased sample of IRDCs in the 10 deg < |l| < 65 deg, |b|<1 deg region of the Galactic plane using Spitzer 8micron extinction. However, in absorption studies, intrinsic fluctuations in the mid-infrared background can be mis-interpreted as foreground clouds. The main goal of the study presented here is to disentangle real clouds in the Spitzer Dark Cloud (SDC) catalogue from artefacts due to fluctuations in the mid-infrared background. We constructed H_2 column density maps at ~18" resolution using the 160micron and 250micron data from the Herschel Galactic plane survey Hi-GAL. We also developed an automated detection scheme that confirms the existence of a SDC through its association with a peak on these Herschel column density maps. Detection simulations, along with visual inspection of a small sub-sample of SDCs, have been performed to get better insight into the limitations of our automated identification scheme. Our analysis shows that 76(+/-19)% of the catalogued SDCs are real. This fraction drops to 55(+/-12)% for clouds with angular diameters larger than ~1 arcminute. The contamination of the PF09 catalogue by large spurious sources reflect the large uncertainties associated to the construction of the 8micron background emission, a key stage towards the identification of SDCs. A comparison of the Herschel confirmed SDC sample with the BGPS and ATLASGAL samples shows that SDCs probe a unique range of cloud properties, reaching down to more compact and lower column density clouds than any of these two (sub-)millimetre Galactic plane surveys.

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.

Much of the dense gas in molecular clouds has a filamentary structure but the detailed structure and evolution of this gas is poorly known. We have observed 54 cores in infrared dark clouds (IRDCs) using N$_2$H$^+$ (1-0) and (3-2) to determine the kinematics of the densest material, where stars will form. We also observed N$_2$D$^+$ (3-2) towards 29 of the brightest peaks to analyse the level of deuteration which is an excellent probe of the quiescent of the early stages of star formation. There were 13 detections of N$_2$D$^+$ (3-2). This is one of the largest samples of IRDCs yet observed in these species. The deuteration ratio in these sources ranges between 0.003 and 0.14. For most of the sources the material traced by N$_2$D$^+$ and N$_2$H$^+$ (3-2) still has significant turbulent motions, however three objects show subthermal N$_2$D$^+$ velocity dispersion. Surprisingly the presence or absence of an embedded 70 $\mu$m source shows no correlation with the detection of N$_2$D$^+$ (3-2), nor does it correlate with any change in velocity dispersion or excitation temperature. Comparison with recent models of deuteration suggest evolutionary time-scales of these regions of several free-fall times or less.

The surroundings of HII regions can have a profound influence on their development, morphology, and evolution. This paper explores the effect of the environment on H II regions in the MonR2 molecular cloud. We aim to investigate the density structure of envelopes surrounding HII regions and to determine their collapse and ionisation expansion ages. The Mon R2 molecular cloud is an ideal target since it hosts an H II region association. Column density and temperature images derived from Herschel data were used together to model the structure of HII bubbles and their surrounding envelopes. The resulting observational constraints were used to follow the development of the Mon R2 ionised regions with analytical calculations and numerical simulations. The four hot bubbles associated with H II regions are surrounded by dense, cold, and neutral gas envelopes. The radial density profiles are reminiscent of those of low-mass protostellar envelopes. The inner parts of envelopes of all four HII regions could be free-falling because they display shallow density profiles. As for their outer parts, the two compact HII regions show a density profile, which is typical of the equilibrium structure of an isothermal sphere. In contrast, the central UCHii region shows a steeper outer profile, that could be interpreted as material being forced to collapse. The size of the heated bubbles, the spectral type of the irradiating stars, and the mean initial neutral gas density are used to estimate the ionisation expansion time, texp, 0.1Myr,for the dense UCHII and compact HII regions and 0.35 Myr for the extended HII region. The envelope transition radii between the shallow and steeper density profiles are used to estimate the time elapsed since the formation of the first proto stellar embryo, Tinf : 1Myr, for the ultra-compact, 1.5 / 3Myr for the compact, and greater than 6Myr for the extended HII regions.

Two major features of the prestellar CMF are: 1) a broad peak below 1 Msun, presumably corresponding to a mean gravitational fragmentation scale, and 2) a characteristic power-law slope, very similar to the Salpeter slope of the stellar initial mass function (IMF) at the high-mass end. While recent Herschel observations have shown that the peak of the prestellar CMF is close to the thermal Jeans mass in marginally supercritical filaments, the origin of the power-law tail of the CMF/IMF at the high-mass end is less clear. Inutsuka (2001) proposed a theoretical scenario in which the origin of the power-law tail can be understood as resulting from the growth of an initial spectrum of density perturbations seeded along the long axis of filaments by interstellar turbulence. Here, we report the statistical properties of the line-mass fluctuations of filaments in nearby molecular clouds observed with Herschel using a 1-D power spectrum analysis. The observed filament power spectra were fitted by a power-law function $(P_{true}(s) \propto s^{\alpha})$ after removing the effect of beam convolution at small scales. A Gaussian-like distribution of power-spectrum slopes was found centered at -1.6, close to that of the one-dimensional velocity power spectrum generated by subsonic Kolomogorov turbulence (-1.67). An empirical correlation, $P^{0.5}(s_0) \propto <N_{\rm H_2}>^{1.4 \pm 0.1} $, was also found between the amplitude of each filament power spectrum $P(s_0)$ and the mean column density along the filament $<N_{\rm H_2}>$. Finally, the dispersion of line-mass fluctuations along each filament $\sigma_{\rm M_{line}}$ was found to scale with the physical length $L$ of the filament, roughly as $\sigma_{M_{line}} \propto L^{0.7}$. Overall, our results are consistent with the suggestion that the bulk of the CMF/IMF results from the gravitational fragmentation of filaments.

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-{\mu}m continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-{\sigma} surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

We present and discuss the results of the Herschel Gould Belt survey observations in a ~11 deg^2 area of the Aquila molecular cloud complex at d~260 pc, imaged with the SPIRE/PACS cameras from 70 to 500 micron. We identify a complete sample of starless dense cores and embedded protostars in this region, and analyze their global properties and spatial distributions. We find a total of 651 starless cores, ~60% of which are gravitationally bound prestellar cores, and they will likely form stars in the future. We also detect 58 protostellar cores. The core mass function (CMF) derived for the prestellar cores is very similar in shape to the stellar initial mass function (IMF), supporting the earlier view that there is a close physical link between the IMF and the CMF. The global shift in mass scale observed between the CMF and the IMF is consistent with a typical star formation efficiency of ~40%. By comparing the numbers of starless cores to the number of young stellar objects, we estimate that the lifetime of prestellar cores is ~1 Myr. We find a strong correlation between the spatial distribution of prestellar cores and the densest filaments. About 90% of the Herschel-identified prestellar cores are located above a background column density corresponding to A_V~7, and ~75% of them lie within filamentary structures with supercritical masses per unit length >~16 M_sun/pc. These findings support a picture wherein the cores making up the peak of the CMF (and probably responsible for the base of the IMF) result primarily from the gravitational fragmentation of marginally supercritical filaments. Given that filaments appear to dominate the mass budget of dense gas at A_V>7, our findings also suggest that the physics of prestellar core formation within filaments is responsible for a characteristic "efficiency" SFR/M_dense ~5+-2 x 10^-8 yr^-1 for the star formation process in dense gas.

We present a catalogue of starless and protostellar clumps associated with infrared dark clouds (IRDCs) in a 40 degrees wide region of the inner Galactic Plane (b<1). We have extracted the far-infrared (FIR) counterparts of 3493 IRDCs with known distance in the Galactic longitude range 15<l<55 and searched for the young clumps using Hi-GAL, the survey of the Galactic Plane carried out with the Herschel satellite. Each clump is identified as a compact source detected at 160, 250 and 350 mum. The clumps have been classified as protostellar or starless, based on their emission (or lack of emission) at 70 mum. We identify 1723 clumps, 1056 (61%) of which are protostellar and 667 (39%) starless. These clumps are found within 764 different IRDCs, 375 (49%) of which are only associated with protostellar clumps, 178 (23%) only with starless clumps, and 211 (28%) with both categories of clumps. The clumps have a median mass of 250 M_sun and range up to >10^4$ M_sun in mass and up to 10^5 L_sun in luminosity. The mass-radius distribution shows that almost 30% of the starless clumps identified in this survey could form high-mass stars, however these massive clumps are confined in only ~4% of the IRDCs. Assuming a minimum mass surface density threshold for the formation of high-mass stars, the comparison of the numbers of massive starless clumps and those already containing embedded sources suggests an upper limit lifetime for the starless phase of 10^5 years for clumps with a mass M>500 M_sun.

We analyse column density and temperature maps derived from Herschel dust continuum observations of a sample of massive infrared dark clouds (G11.11-0.12, G18.82-0.28, G28.37+0.07, G28.53-0.25). We disentangle the velocity structure of the clouds using 13CO 1-0 and 12CO 3-2 data, showing that these IRDCs are the densest regions in massive giant molecular clouds and not isolated features. The probability distribution function (PDF) of column densities for all clouds have a power-law distribution over all (high) column densities, regardless of the evolutionary stage of the cloud: G11.11-0.12, G18.82-0.28, and G28.37+0.07 contain (proto)-stars, while G28.53-0.25 shows no signs of star formation. This is in contrast to the purely log-normal PDFs reported for near/mid-IR extinction maps. We only find a log-normal distribution for lower column densities, if we perform PDFs of the column density maps of the whole GMC in which the IRDCs are embedded. By comparing the PDF slope and the radial column density profile, we attribute the power law to the effect of large-scale gravitational collapse and to local free-fall collapse of pre- and protostellar cores. Independent from the PDF analysis, we find infall signatures in the spectral profiles of 12CO for G28.37+0.07 and G11.11-0.12, supporting the scenario of gravitational collapse. IRDCs are the densest regions within GMCs, which may be the progenitors of massive stars or clusters. At least some of the IRDCs are probably the same features as ridges (high column density regions with N>1e23 cm-2 over small areas), which were defined for nearby IR-bright GMCs. Because IRDCs are only confined to the densest (gravity dominated) cloud regions, the PDF constructed from this kind of a clipped image does not represent the (turbulence dominated) low column density regime of the cloud.

We present the first results from the MALT-45 (Millimetre Astronomer's Legacy Team - 45 GHz) Galactic Plane survey. We have observed 5 square-degrees ($l = 330 - 335$, $b = \pm0.5$) for spectral lines in the 7 mm band (42-44 and 48-49 GHz), including $\text{CS}$ $(1-0)$, class I $\text{CH}_3\text{OH}$ masers in the $7(0,7)-6(1,6)$ $\text{A}^{+}$ transition and $\text{SiO}$ $(1-0)$ $v=0,1,2,3$. MALT-45 is the first unbiased, large-scale, sensitive spectral line survey in this frequency range. In this paper, we present data from the survey as well as a few intriguing results; rigorous analyses of these science cases are reserved for future publications. Across the survey region, we detected 77 class I $\text{CH}_3\text{OH}$ masers, of which 58 are new detections, along with many sites of thermal and maser $\text{SiO}$ emission and thermal $\text{CS}$. We found that 35 class I $\text{CH}_3\text{OH}$ masers were associated with the published locations of class II $\text{CH}_3\text{OH}$, $\text{H}_2\text{O}$ and $\text{OH}$ masers but 42 have no known masers within 60 arcsec. We compared the MALT-45 $\text{CS}$ with $\text{NH}_3$ (1,1) to reveal regions of $\text{CS}$ depletion and high opacity, as well as evolved star-forming regions with a high ratio of $\text{CS}$ to $\text{NH}_3$. All $\text{SiO}$ masers are new detections, and appear to be associated with evolved stars from the $\it{Spitzer}$ Galactic Legacy Infrared Mid-Plane Survey Extraordinaire (GLIMPSE). Generally, within $\text{SiO}$ regions of multiple vibrational modes, the intensity decreases as $v=1,2,3$, but there are a few exceptions where $v=2$ is stronger than $v=1$.

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

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.

We use PACS and SPIRE continuum data at 160 um, 250 um, 350 um, and 500 um from the Herschel Gould Belt Survey to sample seven clumps in Perseus: B1, B1-E, B5, IC348, L1448, L1455, and NGC1333. Additionally, we identify and characterize the embedded Class 0 protostars using detections of compact Herschel sources at 70 um as well as archival Spitzer catalogues and SCUBA 850 um photometric data. We identify 28 candidate Class 0 protostars, four of which are newly discovered sources not identified with Spitzer. We find that the star formation efficiency of clumps, as traced by Class 0 protostars, correlates strongly with the flatness of their respective column density distributions at high values. This correlation suggests that the fraction of high column density material in a clump reflects only its youngest protostellar population rather than its entire source population. We propose that feedback from either the formation or evolution of protostars changes the local density structure of clumps.

The density and temperature structures of dense cores in the L1495 cloud of the Taurus star-forming region are investigated using Herschel SPIRE and PACS images in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m and 500 $\mu$m continuum bands. A sample consisting of 20 cores, selected using spectral and spatial criteria, is analysed using a new maximum likelihood technique, COREFIT, which takes full account of the instrumental point spread functions. We obtain central dust temperatures, $T_0$, in the range 6-12 K and find that, in the majority of cases, the radial density falloff at large radial distances is consistent with the $r^{-2}$ variation expected for Bonnor-Ebert spheres. Two of our cores exhibit a significantly steeper falloff, however, and since both appear to be gravitationally unstable, such behaviour may have implications for collapse models. We find a strong negative correlation between $T_0$ and peak column density, as expected if the dust is heated predominantly by the interstellar radiation field. At the temperatures we estimate for the core centres, carbon-bearing molecules freeze out as ice mantles on dust grains, and this behaviour is supported here by the lack of correspondence between our estimated core locations and the previously-published positions of H$^{13}$CO$^+$ peaks. On this basis, our observations suggest a sublimation-zone radius typically $\sim 10^4$ AU. Comparison with previously-published N$_2$H$^+$ data at 8400 AU resolution, however, shows no evidence for N$_2$H$^+$ depletion at that resolution.

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $\kappa_{\lambda} =0.1(\lambda/300 \mu m)^{-2}$, agrees to within 50% with the dust extinction constraints

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

Formation of stars is now believed to be tightly linked to the dynamical evolution of interstellar filaments in which they form. In this paper we analyze the density structure and kinematics of a small network of infrared dark filaments, SDC13, observed in both dust continuum and molecular line emission with the IRAM 30m telescope. These observations reveal the presence of 18 compact sources amongst which the two most massive, MM1 and MM2, are located at the intersection point of the parsec-long filaments. The dense gas velocity and velocity dispersion observed along these filaments show smooth, strongly correlated, gradients. We discuss the origin of the SDC13 velocity field in the context of filament longitudinal collapse. We show that the collapse timescale of the SDC13 filaments (from 1 Myr to 4 Myr depending on the model parameters) is consistent with the presence of Class I sources in them, and argue that, on top of bringing more material to the centre of the system, collapse could generate additional kinematic support against local fragmentation, helping the formation of starless super-Jeans cores.

We present the Herschel Gould Belt survey maps of the L1641 molecular clouds in Orion A. We extracted both the filaments and dense cores in the region. We identified which of dense sources are proto- or pre-stellar, and studied their association with the identified filaments. We find that although most (71%) of the pre-stellar sources are located on filaments there is still a significant fraction of sources not associated with such structures. We find that these two populations (on and off the identified filaments) have distinctly different mass distributions. The mass distribution of the sources on the filaments is found to peak at 4 Solar masses and drives the shape of the CMF at higher masses, which we fit with a power law of the form dN/dlogM \propto M^{-1.4+/-0.4}. The mass distribution of the sources off the filaments, on the other hand, peaks at 0.8 Solar masses and leads to a flattening of the CMF at masses lower than ~4 Solar masses. We postulate that this difference between the mass distributions is due to the higher proportion of gas that is available in the filaments, rather than in the diffuse cloud.

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.

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.

The whole of the Taurus region (a total area of 52 sq. deg.) has been observed by the Herschel SPIRE and PACS instruments at wavelengths of 70, 160, 250, 350 and 500 {\mu}m as part of the Herschel Gould Belt Survey. In this paper we present the first results from the part of the Taurus region that includes the Barnard 18 and L1536 clouds. A new source-finding routine, the Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based on CLUMPFIND, but that also generates a structure tree, or dendrogram, which can be used to interpret hierarchical clump structure in a complex region. Sources were extracted from the data using the hierarchical version of CSAR and plotted on a mass-size diagram. We found a hierarchy of objects with sizes in the range 0.024-2.7 pc. Previous studies showed that gravitationally bound prestellar cores and unbound starless clumps appeared in different places on the mass-size diagram. However, it was unclear whether this was due to a lack of instrumental dynamic range or whether they were actually two distinct populations. The excellent sensitivity of Herschel shows that our sources fill the gap in the mass-size plane between starless and pre-stellar cores, and gives the first clear supporting observational evidence for the theory that unbound clumps and (gravitationally bound) prestellar cores are all part of the same population, and hence presumably part of the same evolutionary sequence (c.f. Simpson et al. 2011).

A key parameter to the description of all star formation processes is the density structure of the gas. In this letter, we make use of probability distribution functions (PDFs) of Herschel column density maps of Orion B, Aquila, and Polaris, obtained with the Herschel Gould Belt survey (HGBS). We aim to understand which physical processes influence the PDF shape, and with which signatures. The PDFs of Orion B (Aquila) show a lognormal distribution for low column densities until Av 3 (6), and a power-law tail for high column densities, consistent with a rho r^-2 profile for the equivalent spherical density distribution. The PDF of Orion B is broadened by external compression due to the nearby OB stellar aggregates. The PDF of a quiescent subregion of the non-star-forming Polaris cloud is nearly lognormal, indicating that supersonic turbulence governs the density distribution. But we also observe a deviation from the lognormal shape at Av>1 for a subregion in Polaris that includes a prominent filament. We conclude that (i) the point where the PDF deviates from the lognormal form does not trace a universal Av-threshold for star formation, (ii) statistical density fluctuations, intermittency and magnetic fields can cause excess from the lognormal PDF at an early cloud formation stage, (iii) core formation and/or global collapse of filaments and a non-isothermal gas distribution lead to a power-law tail, and (iv) external compression broadens the column density PDF, consistent with numerical simulations.

We investigate the gas velocity dispersions of a sample of filaments recently detected as part of the Herschel Gould Belt Survey in the IC5146, Aquila, and Polaris interstellar clouds. To measure these velocity dispersions, we use 13CO, C18O, and N2H+ line observations obtained with the IRAM 30m telescope. Correlating our velocity dispersion measurements with the filament column densities derived from Herschel data, we show that interstellar filaments can be divided into two regimes: thermally subcritical filaments, which have transonic velocity dispersions (c_s ~< \sigma_tot < 2 c_s) independent of column density, and are gravitationally unbound; and thermally supercritical filaments, which have higher velocity dispersions scaling roughly as the square root of column density (\sigma_tot ~ \Sigma^0.5), and are self-gravitating. The higher velocity dispersions of supercritical filaments may not directly arise from supersonic interstellar turbulence but may be driven by gravitational contraction/accretion. Based on our observational results, we propose an evolutionary scenario whereby supercritical filaments undergo gravitational contraction and increase in mass per unit length through accretion of background material while remaining in rough virial balance. We further suggest that this accretion process allows supercritical filaments to keep their approximately constant inner widths (~ 0.1 pc) while contracting.

We present first results from the Herschel Gould Belt survey for the B211/L1495 region in the Taurus molecular cloud. Thanks to their high sensitivity and dynamic range, the Herschel images reveal the structure of the dense, star-forming filament B211 with unprecedented detail, along with the presence of striations perpendicular to the filament and generally oriented along the magnetic field direction as traced by optical polarization vectors. Based on the column density and dust temperature maps derived from the Herschel data, we find that the radial density profile of the B211 filament approaches a power-law behavior {\rho} {\propto} r^(-2.0{\pm}0.4) at large radii and that the temperature profile exhibits a marked drop at small radii. The observed density and temperature profiles of the B211 filament are in good agreement with a theoretical model of a cylindrical filament undergoing gravitational contraction with a polytropic equation of state: P {\propto} {\rho}^{\gamma} and T {\propto} {\rho}^({\gamma}-1), with {\gamma}=0.97{\pm}0.01<1 (i.e. not strictly isothermal). The morphology of the column density map, where some of the perpendicular striations are apparently connected to the B211 filament, further suggests that the material may be accreting along the striations onto the main filament. The typical velocities expected for the infalling material in this picture are ~0.5-1 km/s, which are consistent with the existing kinematical constraints from previous CO observations.

We present a study of the star formation histories of the Lupus I, III, and IV clouds using the Herschel 70-500 micron maps obtained by the Herschel Gould Belt Survey Key Project. By combining the new Herschel data with the existing Spitzer catalog we obtained an unprecedented census of prestellar sources and young stellar objects in the Lupus clouds, which allowed us to study the overall star formation rate (SFR) and efficiency (SFE). The high SFE of Lupus III, its decreasing SFR, and its large number of pre-main sequence stars with respect to proto- and prestellar sources, suggest that Lupus III is the most evolved cloud, and after having experienced a major star formation event in the past, is now approaching the end of its current star-forming cycle. Lupus I is currently undergoing a large star formation event, apparent by the increasing SFR, the large number of prestellar objects with respect to more evolved objects, and the high percentage of material at high extinction (e.g., above A_V=8 mag). Also Lupus IV has an increasing SFR; however, the relative number of prestellar sources is much lower, suggesting that its star formation has not yet reached its peak.

We present APEX/P-ArT\'eMiS 450\mu m continuum observations of RCW 36 and the adjacent ridge, a high-mass high-column density filamentary structure at the centre of the Vela C molecular cloud. These observations, at higher resolution than Herschel's SPIRE camera, reveal clear fragmentation of the central star-forming ridge. Combined with PACS far-infrared and SPIRE sub-millimetre observations from the Herschel HOBYS project we build a high resolution column density map of the region mapped with P-ArT\'eMiS. We extract the radial density profile of the Vela C ridge which with a ~ 0.1pc central width is consistent with that measured for low-mass star-forming filaments in the Herschel Gould Belt survey. Direct comparison with Serpens South, of the Gould Belt Aquila complex, reveals many similarities between the two regions. Despite likely different formation mechanisms and histories, the Vela C ridge and Serpens South filament share common characteristics, including their filament central widths.

We report far-IR Herschel observations obtained between 70 $\mu$m and 500 $\mu$m of two star-forming dusty condensations, B1-bS and B1-bN, in the B1 region of the Perseus star-forming cloud. In the Western part of the Perseus cloud, B1-bS is the only source detected in all of the 6 PACS and SPIRE photometric bands without being visible in the Spitzer map at 24 $\mu$m. B1-bN is clearly detected between 100 $\mu$m and 250 $\mu$m. We have fitted the spectral energy distributions of these sources to derive their physical properties, and find that a simple greybody model fails to reproduce the observed SEDs. At least a two-component model, consisting of a central source surrounded by a dusty envelope is required. The properties derived from the fit, however, suggest that the central source is not a Class 0 object. We then conclude that while B1-bS and B1-bN appear to be more evolved than a pre-stellar core, the best-fit models suggest that their central objects are younger than a Class 0 source. Hence, they may be good candidates to be examples of the first hydrostatic core phase. The projected distance between B1-bS and B1-bN is a few Jeans lengths. If their physical separation is close to this value, this pair would allow the mutual interactions between two forming stars at a very early stage of their evolution to be studied.

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.

In order to characterise the cloud structures responsible for the formation of high-mass stars, we present Herschel observations of the DR21 environment. Maps of the column density and dust temperature unveil the structure of the DR21 ridge and several connected filaments. The ridge has column densities larger than 1e23/cm^2 over a region of 2.3 pc^2. It shows substructured column density profiles and branching into two major filaments in the north. The masses in the studied filaments range between 130 and 1400 Msun whereas the mass in the ridge is 15000 Msun. The accretion of these filaments onto the DR21 ridge, suggested by a previous molecular line study, could provide a continuous mass inflow to the ridge. In contrast to the striations seen in e.g., the Taurus region, these filaments are gravitationally unstable and form cores and protostars. These cores formed in the filaments potentially fall into the ridge. Both inflow and collisions of cores could be important to drive the observed high-mass star formation. The evolutionary gradient of star formation running from DR21 in the south to the northern branching is traced by decreasing dust temperature. This evolution and the ridge structure can be explained by two main filamentary components of the ridge that merged first in the south.

In a previous paper we identified cores within infrared dark clouds (IRDCs). We regarded those without embedded sources as the least evolved, and labelled them starless. Here we identify the most isolated starless cores and model them using a three-dimensional, multi-wavelength, Monte Carlo, radiative transfer code. We derive the cores' physical parameters and discuss the relation between the mass, temperature, density, size and the surrounding interstellar radiation field (ISRF) for the cores. The masses of the cores were found not to correlate with their radial size or central density. The temperature at the surface of a core was seen to depend almost entirely on the level of the ISRF surrounding the core. No correlation was found between the temperature at the centre of a core and its local ISRF. This was seen to depend, instead, on the density and mass of the core.

We present Herschel images from the HOBYS key program of the Eagle Nebula (M16) in the far-infrared and sub-millimetre, using the PACS and SPIRE cameras at 70{\mu}m, 160{\mu}m, 250{\mu}m, 350{\mu}m, 500{\mu}m. M16, home to the Pillars of Creation, is largely under the influence of the nearby NGC6611 high-mass star cluster. The Herschel images reveal a clear dust temperature gradient running away from the centre of the cavity carved by the OB cluster. We investigate the heating effect of NGC6611 on the entire M16 star-forming complex seen by Herschel including the diffuse cloud environment and the dense filamentary structures identified in this region. In addition, we interpret the three-dimensional geometry of M16 with respect to the nebula, its surrounding environment, and the NGC6611 cavity. The dust temperature and column density maps reveal a prominent eastern filament running north-south and away from the high-mass star-forming central region and the NGC6611 cluster, as well as a northern filament which extends around and away from the cluster. The dust temperature in each of these filaments decreases with increasing distance from the NGC6611 cluster, indicating a heating penetration depth of \sim 10 pc in each direction in 3 - 6 \times 10^{22} cm-2 column density filaments. We show that in high-mass star-forming regions OB clusters impact the temperature of future star-forming sites, modifying the initial conditions for collapse and effecting the evolutionary criteria of protostars developed from spectral energy distributions. Possible scenarios for the origin of the morphology seen in this region are discussed, including a western equivalent to the eastern filament, which was destroyed by the creation of the OB cluster and its subsequent winds and radiation.

We present continuum observations of the Perseus B1-E region from the Herschel Gould Belt Survey. These Herschel data reveal a loose grouping of substructures at 160 - 500 micron not seen in previous submillimetre observations. We measure temperature and column density from these data and select the nine densest and coolest substructures for follow-up spectral line observations with the Green Bank Telescope. We find that the B1-E clump has a mass of ~ 100 solar masses and appears to be gravitationally bound. Furthermore, of the nine substructures examined here, one substructure (B1-E2) appears to be itself bound. The substructures are typically less than a Jeans length from their nearest neighbour and thus, may interact on a timescale of ~ 1 Myr. We propose that B1-E may be forming a first generation of dense cores, which could provide important constraints on the initial conditions of prestellar core formation. Our results suggest that B1-E may be influenced by a strong, localized magnetic field, but further observations are still required.

For many years feedback processes generated by OB-stars in molecular clouds, including expanding ionization fronts, stellar winds, or UV-radiation, have been proposed to trigger subsequent star formation. However, hydrodynamic models including radiation and gravity show that UV-illumination has little or no impact on the global dynamical evolution of the cloud. The Rosette molecular cloud, irradiated by the NGC2244 cluster, is a template region for triggered star-formation, and we investigated its spatial and density structure by applying a curvelet analysis, a filament-tracing algorithm (DisPerSE), and probability density functions (PDFs) on Herschel column density maps, obtained within the HOBYS key program. The analysis reveals not only the filamentary structure of the cloud but also that all known infrared clusters except one lie at junctions of filaments, as predicted by turbulence simulations. The PDFs of sub-regions in the cloud show systematic differences. The two UV-exposed regions have a double-peaked PDF we interprete as caused by shock compression. The deviations of the PDF from the log-normal shape typically associated with low- and high-mass star-forming regions at Av~3-4m and 8-10m, respectively, are found here within the very same cloud. This shows that there is no fundamental difference in the density structure of low- and high-mass star-forming regions. We conclude that star-formation in Rosette - and probably in high-mass star-forming clouds in general - is not globally triggered by the impact of UV-radiation. Moreover, star formation takes place in filaments that arose from the primordial turbulent structure built up during the formation of the cloud. Clusters form at filament mergers, but star formation can be locally induced in the direct interaction zone between an expanding HII--region and the molecular cloud.

A growing body of evidence indicates that the formation of filaments in interstellar clouds is a key component of the star formation process. In this paper, we present new Herschel PACS and SPIRE observations of the B59 and Stem regions in the Pipe Nebula complex, revealing a rich, organized network of filaments. The asymmetric column density profiles observed for several filaments, along with the bow-like edge of B59, indicates that the Pipe Nebula is being compressed from its western side, most likely by the winds from the nearby Sco OB2 association. We suggest that this compressive flow has contributed to the formation of some of the observed filamentary structures. In B59, the only region of the entire Pipe complex showing star formation activity, the same compressive flow has likely enhanced the initial column density of the clump, allowing it to become globally gravitationally unstable. Although more speculative, we propose that gravity has also been responsible for shaping the converging filamentary pattern observed in B59. While the question of the relative impact of large-scale compression and gravity remains open in B59, large-scale compression appears to be a plausible mechanism for the initial formation of filamentary structures in the rest of the complex

We have used data taken as part of the Herschel infrared Galactic Plane survey (Hi-GAL) to study 3171 infrared-dark cloud (IRDC) candidates that were identified in the mid-infrared (8 {\mu}m) by Spitzer (we refer to these as 'Spitzer-dark' regions). They all lie in the range l=300 - 330 \circ and |b| 6 1 \circ. Of these, only 1205 were seen in emission in the far-infrared (250-500 {\mu}m) by Herschel (we call these 'Herschel-bright' clouds). It is predicted that a dense cloud will not only be seen in absorption in the mid-infrared, but will also be seen in emission in the far-infrared at the longest Herschel wavebands (250-500 {\mu}m). If a region is dark at all wavelengths throughout the mid-infrared and far-infrared, then it is most likely to be simply a region of lower background infrared emission (a 'hole in the sky'). Hence, it appears that previous surveys, based on Spitzer and other mid-infrared data alone, may have over-estimated the total IRDC population by a factor of 2. This has implications for estimates of the star formation rate in IRDCs in the Galaxy.We studied the 1205 Herschel-bright IRDCs at 250 {\mu}m, and found that 972 of them had at least one clearly defined 250-{\mu}m peak, indicating that they contained one or more dense cores. Of these, 653 (67 per cent) contained an 8-{\mu}m point source somewhere within the cloud, 149 (15 per cent) contained a 24-{\mu}m point source but no 8-{\mu}m source, and 170 (18 per cent) contained no 24-{\mu}m or 8-{\mu}m point sources. We use these statistics to make inferences about the lifetimes of the various evolutionary stages of IRDCs.

The filament IRDC G035.39--00.33 in the W48 molecular complex is one of the darkest infrared clouds observed by \textit{Spitzer}. It has been observed by the PACS (70 and 160\,$\micron$) and SPIRE (250, 350, and 500\,$\micron$) cameras of the \textit{Herschel} Space Observatory as part of the W48 molecular cloud complex in the framework of the HOBYS key programme. The observations reveal a sample of 28 compact sources (deconvolved FWHM sizes $<$0.3 pc) complete down to $\sim$$5 \msun$ in G035.39--00.33 and its surroundings. Among them, 13 compact sources are massive dense cores with masses $>$$20 \msun$. The cloud characteristics we derive from the analysis of their spectral energy distributions are masses of $20-50 \msun$, sizes of 0.1--0.2 pc, and average densities of $2-20 \times 10^{5} \cmc$, which make these massive dense cores excellent candidates to form intermediate- to high-mass stars. Most of the massive dense cores are located inside the G035.39--00.33 ridge and host IR-quiet high-mass protostars. The large number of protostars found in this filament suggests that we are witnessing a mini-burst of star formation with an efficiency of $\sim$15% and a rate density of $\sim$$40 \msun\,$yr$^{-1}\,$kpc$^{-2}$ within $\sim$8 pc$^2$, a large area covering the full ridge. Part of the extended SiO emission observed towards G035.39--00.33 is not associated with obvious protostars and may originate from low-velocity shocks within converging flows, as advocated by previous studies.

We describe a pilot survey conducted with the Mopra 22-m radio telescope in preparation for the Millimeter Astronomy Legacy Team Survey at 90 GHz (MALT90). We identified 182 candidate dense molecular clumps using six different selection criteria and mapped each source simultaneously in 16 different lines near 90 GHz. We present a summary of the data and describe how the results of the pilot survey shaped the design of the larger MALT90 survey. We motivate our selection of target sources for the main survey based on the pilot detection rates and demonstrate the value of mapping in multiple lines simultaneously at high spectral resolution.

We present the first Herschel PACS and SPIRE results of the Vela C molecular complex in the far-infrared and submillimetre regimes at 70, 160, 250, 350, and 500 um, spanning the peak of emission of cold prestellar or protostellar cores. Column density and multi-resolution analysis (MRA) differentiates the Vela C complex into five distinct sub-regions. Each sub-region displays differences in their column density and temperature probability distribution functions (PDFs), in particular, the PDFs of the `Centre-Ridge' and `South-Nest' sub-regions appear in stark contrast to each other. The Centre-Ridge displays a bimodal temperature PDF representative of hot gas surrounding the HII region RCW 36 and the cold neighbouring filaments, whilst the South-Nest is dominated by cold filamentary structure. The column density PDF of the Centre-Ridge is flatter than the South-Nest, with a high column density tail, consistent with formation through large-scale flows, and regulation by self-gravity. At small to intermediate scales MRA indicates the Centre-Ridge to be twice as concentrated as the South-Nest, whilst on larger scales, a greater portion of the gas in the South-Nest is dominated by turbulence than in the Centre-Ridge. In Vela C, high-mass stars appear to be preferentially forming in ridges, i.e., dominant high column density filaments.

We provide a first look at the results of the Herschel Gould Belt survey toward the IC5146 molecular cloud and present a preliminary analysis of the filamentary structure in this region. The column density map, derived from our 70-500 micron Herschel data, reveals a complex network of filaments, and confirms that these filaments are the main birth sites of prestellar cores. We analyze the column density profiles of 27 filaments and show that the underlying radial density profiles fall off as r^{-1.5} to r^{-2.5} at large radii. Our main result is that the filaments seem to be characterized by a narrow distribution of widths having a median value of 0.10 +- 0.03 pc, which is in stark contrast to a much broader distribution of central Jeans lengths. This characteristic width of ~0.1 pc corresponds to within a factor of ~2 to the sonic scale below which interstellar turbulence becomes subsonic in diffuse gas, supporting the argument that the filaments may form as a result of the dissipation of large-scale turbulence.

Context: Despite their profound effect on the universe, the formation of massive stars and stellar clusters remains elusive. In the past decade, compelling evidence has emerged that suggests Infrared Dark Clouds (IRDCs) may be precursors to stellar clusters. However, the usual method for identifying IRDCs is biased by the requirement that they are seen in absorption against bright mid-IR emission, whereas dust continuum observations allow cold, dense pre-stellar-clusters to be identified anywhere. Aims: We aim to understand what physical properties characterize IRDCs, to explore the population of dust continuum sources that are not IRDCs, and to roughly characterize the star formation activity in dust continuum sources. Results: We present temperature and column density maps in the Hi-GAL l=30 and l=59 SDP fields, as well as a robust algorithm for cirrus subtraction and source identification using Hi-GAL. We report on the fraction of Hi-GAL sources which are mid-IR-dark, -neutral, or -bright in both fields. We find significant trends in column density and temperature between mid-IR-dark and -bright pixels; mid-IR-dark pixels are about 10 K colder and have a factor of 2 higher column density on average than mid-IR-bright pixels. We find that Hi-GAL dust continuum sources span a range of evolutionary states from pre- to star-forming, and that warmer sources are associated with more star formation tracers. There is a trend of increasing temperature with tracer type from mid-IR-dark at the coldest, to outflow/maser sources in the middle, and finally to 8 and 24 micron bright sources at the warmest. Finally, we identify five candidate IRDC-like sources on the far-side of the Galaxy. These are cold (~ 20 K), high column density (N(H2) > 10^22 cm^-2) clouds identified with Hi-GAL which, despite bright surrounding mid-IR emission, show little to no absorption at 8 micron. (abridged)

The densest infrared dark clouds (IRDCs) may represent the earliest observable stage of high-mass star formation. These clouds are very cold, hence they emit mainly at far-infrared and sub-mm wavelengths. For the first time, Herschel has provided multi-wavelength, spatially resolved observations of cores within IRDCs, which, when combined with radiative transfer modelling, can constrain their properties, such as mass, density profile and dust temperature. We use a 3D, multi-wavelength Monte Carlo radiative transfer code to model in detail the emission from six cores in three typical IRDCs seen in the Hi-GAL survey (G030.50+00.95, G031.03+00.26 and G031.03+00.76), and thereby to determine the properties of these cores and compare them with their low-mass equivalents. We found masses ranging from 90 to 290 solar masses with temperatures from 8 to 11K at the centre of each core and 18 to 28K at the surface. The maximum luminosity of an embedded star within each core was calculated, and we rule out the possibility of significant high mass star formation having yet occurred in three of our cores.

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.

How and when the mass distribution of stars in the Galaxy is set is one of the main issues of modern astronomy. Here we present a statistical study of mass and density distributions of infrared dark clouds (IRDCs) and fragments within them. These regions are pristine molecular gas structures and progenitors of stars and so provide insights into the initial conditions of star formation. This study makes use of a IRDC catalogue (Peretto & Fuller 2009), the largest sample of IRDC column density maps to date, containing a total of ~11,000 IRDCs with column densities exceeding N_{H2} = 1 X10^{22} cm^{-2} and over 50,000 single peaked IRDC fragments. The large number of objects constitutes an important strength of this study, allowing detailed analysis of the completeness of the sample and so statistically robust conclusions. Using a statistical approach to assigning distances to clouds, the mass and density distributions of the clouds and the fragments within them are constructed. The mass distributions show a steepening of the slope when switching from IRDCs to fragments, in agreement with previous results of similar structures. IRDCs and fragments are divided into unbound/bound objects by assuming Larson's relation and calculating their virial parameter. IRDCs are mostly gravitationally bound, while a significant fraction of the fragments are not. The density distribution of gravitationally unbound fragments shows a steep characteristic slope. (see paper for full Abstract).

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.

We summarize the first results from the Gould Belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demonstration phase' of Herschel. Our 70-500 micron images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class 0 protostars can be identified in the Aquila field, while ~ unbound starless cores and no protostars are observed in the Polaris field. The prestellar core mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability.

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

As part of the science demonstration phase of the Herschel mission of the Gould Belt Key Program, the Aquila Rift molecular complex has been observed. The complete ~ 3.3deg x 3.3deg imaging with SPIRE 250/350/500 micron and PACS 70/160 micron allows a deep investigation of embedded protostellar phases, probing of the dust emission from warm inner regions at 70 and 160 micron to the bulk of the cold envelopes between 250 and 500 micron. We used a systematic detection technique operating simultaneously on all Herschel bands to build a sample of protostars. Spectral energy distributions are derived to measure luminosities and envelope masses, and to place the protostars in an M_env - L_bol evolutionary diagram. The spatial distribution of protostars indicates three star-forming sites in Aquila, with W40/Sh2-64 HII region by far the richest. Most of the detected protostars are newly discovered. For a reduced area around the Serpens South cluster, we could compare the Herschel census of protostars with Spitzer results. The Herschel protostars are younger than in Spitzer with 7 Class 0 YSOs newly revealed by Herschel. For the entire Aquila field, we find a total of ~ 45-60 Class 0 YSOs discovered by Herschel. This confirms the global statistics of several hundred Class~0 YSOs that should be found in the whole Gould Belt survey.

Infrared dark clouds (IRDCs) are cold and dense reservoirs of gas potentially available to form stars. Many of these clouds are likely to be pristine structures representing the initial conditions for star formation. The study presented here aims to construct and analyze accurate column density and dust temperature maps of IRDCs by using the first Herschel data from the Hi-GAL galactic plane survey. These fundamental quantities, are essential for understanding processes such as fragmentation in the early stages of the formation of stars in molecular clouds. We have developed a simple pixel-by-pixel SED fitting method, which accounts for the background emission. By fitting a grey-body function at each position, we recover the spatial variations in both the dust column density and temperature within the IRDCs. This method is applied to a sample of 22 IRDCs exhibiting a range of angular sizes and peak column densities. Our analysis shows that the dust temperature decreases significantly within IRDCs, from background temperatures of 20-30 K to minimum temperatures of 8-15 K within the clouds, showing that dense molecular clouds are not isothermal. Temperature gradients have most likely an important impact on the fragmentation of IRDCs. Local temperature minima are strongly correlated with column density peaks, which in a few cases reach NH2 = 1 x 10^{23} cm^{-2}, identifying these clouds as candidate massive prestellar cores. Applying this technique to the full Hi-GAL data set will provide important constraints on the fragmentation and thermal properties of IRDCs, and help identify hundreds of massive prestellar core candidates.

The majority of stars form in clusters. Therefore a comprehensive view of star formation requires understanding the initial conditions for cluster formation. The goal of our study is to shed light on the physical properties of infrared dark clouds (IRDCs) and the role they play in the formation of stellar clusters. This article, the first of a series dedicated to the study of IRDCs, describes techniques developed to establish a complete catalogue of Spitzer IRDCs in the Galaxy. We have analysed Spitzer GLIMPSE and MIPSGAL data to identify a complete sample of IRDCs in the region of Galactic longitude and latitude 10deg < |l|<65deg and |b|<1deg. From the 8micron observations we have constructed opacity maps and used a newly developed extraction algorithm to identify structures above a column density of N_{H2} > 1x10^{22} cm^{-2}. The 24micron data are then used to characterize the star formation activity of each extracted cloud. A total of 11303 clouds have been extracted. A comparison with the existing MSX based catalogue of IRDCs shows that 80% of these Spitzer dark clouds were previously unknown. The algorithm also extracts ~ 20000 to 50000 fragments within these clouds, depending on detection threshold used. A first look at the MIPSGAL data indicates that between 20% and 68% of these IRDCs show 24micron point-like association. This new database provides an important resource for future studies aiming to understand the initial conditions of star formation in the Galaxy.

The earliest phases of clustered star formation and the origin of the stellar initial mass function (IMF) are currently much debated. In order to constrain the origin of the IMF, we investigated the internal and relative motions of starless condensations and protostars previously detected by us in the dust continuum at 1.2mm in the L1688 protocluster of the Ophiuchus molecular cloud complex. The starless condensations have a mass spectrum resembling the IMF and are therefore likely representative of the initial stages of star formation in the protocluster. We carried out detailed molecular line observations, including some N2H+(1-0) mapping, of the Ophiuchus protocluster condensations using the IRAM 30m telescope. We measured subsonic or at most transonic levels of internal turbulence within the condensations, implying virial masses which generally agree within a factor of ~ 2 with the masses derived from the 1.2mm dust continuum. This supports the notion that most of the L1688 starless condensations are gravitationally bound and prestellar in nature. We measured a global one-dimensional velocity dispersion of less than 0.4 km/s between condensations. This small relative velocity dispersion implies that, in general, the condensations do not have time to interact with one another before evolving into pre-main sequence objects. Our observations support the view that the IMF is partly determined by cloud fragmentation at the prestellar stage. Competitive accretion is unlikely to be the dominant mechanism at the protostellar stage in the Ophiuchus protocluster, but it may possibly govern the growth of starless, self-gravitating condensations initially produced by gravoturbulent fragmentation toward an IMF, Salpeter-like mass spectrum.

We measure the mass and kinematics of the massive molecular torus in the planetary nebula NGC 6302. The nebula is the proto-typical butterfly nebula. The origin of the wing-like morphology is disputed: determining the mass-loss history of the confining torus is an important step in understanding the formation of this structure. We performed submillimeter observations with JCMT and the SMA interferometer. The continuum emission as well as the J=2-1 and 3-2 transitions of 12CO and 13CO are analysed at arcsecond resolution. The CO emission indicates a mass of the torus of ~ 2Msun +/- 1Msun. The 12CO and 13CO emission matches the dark lane seen in absorption in the Halpha image of the object. The CO torus is expanding with a velocity of ~ 8 km/s, centred at Vlsr=-31.5 km/s. The size and expansion velocity of the torus indicates that the torus was ejected from ~ 7500 yr to 2900 yr ago, with a mass-loss rate of 5x10^{-4}Msun/yr. We also see a ballistic component in the CO images with a velocity gradient of 140 km/s/pc. The derived mass-loss history of the torus favours binary interaction as the cause of the ejection of the torus. We predict the existence of a companion with an orbital period P < 1 month.

The earliest phases of massive star formation in clusters are still poorly understood. Here, we test the hypothesis for high-mass star formation proposed in our earlier paper (Peretto et al. 2006). In order to confirm the physical validity of this hypothesis, we carried out IRAM Plateau de Bure interferometer observations of NGC 2264-C and performed SPH numerical simulations of the collapse of a Jeans-unstable, prolate dense clump. Our Plateau de Bure observations reveal the presence of a new compact source (C-MM13) located only \~ 10000 AU away, but separated by ~ 1.1 km/s in (projected) velocity, from the most massive Class 0 object (C-MM3) lying at the very center of NGC 2264-C. Detailed comparison with our numerical SPH simulations supports the view that NGC 2264-C is an elongated cluster-forming clump in the process of collapsing and fragmenting along its long axis, leading to a strong dynamical interaction and possible protostar merger in the central region of the clump. The present study also sets several quantitative constraints on the initial conditions of large-scale collapse in NGC 2264-C. Our hydrodynamic simulations indicate that the observed velocity pattern characterizes an early phase of protocluster collapse which survives for an only short period of time (i.e., < 10^5 yr). To provide a good match to the observations the simulations require an initial ratio of turbulent to gravitational energy of only ~ 5 %, which strongly suggests that the NGC 2264-C clump is structured primarily by gravity rather than turbulence. The required "cold'' initial conditions may result from rapid compression by an external trigger.

Classical T Tauri stars are encircled by accretion discs most of the time unresolved by conventional imaging observation. However, numerical simulations show that unresolved aperture linear polarimetry can be used to extract information about the geometry of the immediate circumstellar medium that scatter the starlight. Monin, Menard & Duchene (1998) previously suggested that polarimetry can be used to trace the relative orientation of discs in young binary systems in order to shed light on the stellar and planet formation process. In this paper, we report on new VLT/FORS1 optical linear polarisation measurements of 23 southern binaries spanning a range of separation from 0.8'' to 10''. In each field, the polarisation of the central binary is extracted, as well as the polarisation of nearby stars in order to estimate the local interstellar polarisation. We find that, in general, the linear polarisation vectors of individual components in binary systems tend to be parallel to each other. The amplitude of their polarisations are also correlated. These findings are in agreement with our previous work and extend the trend to smaller separations. They are also similar to other studies, e.g., Donar et al. 1999; Jensen et al. 2000, 2004; Wolf et al. 2001. However, we also find a few systems showing large differences in polarisation level, possibly indicating different inclinations to the line-of-sight for their discs.

We present the results of dust continuum and molecular line observations of two massive cluster-forming clumps, NGC 2264-C and NGC 2264-D, including extensive mapping performed with the MAMBO bolometer array and the HERA heterodyne array on the IRAM 30m telescope. Both NGC 2264 clumps are located in the Mon OB1 giant molecular cloud complex, adjacent to one another. Twelve and fifteen compact millimeter continuum sources (i.e. MMSs) are identified in clumps C and D, respectively. Evidence for widespread infall motions is found in, e.g., HCO+(3-2) or CS(3-2) in both NGC 2264-C and NGC 2264-D. A sharp velocity discontinuity ~ 2 km/s in amplitude is observed in N_2H+(1-0) and H^{13}CO+(1-0) in the central, innermost part of NGC 2264-C, which we interpret as the signature of a strong dynamical interaction between two MMSs and their possible merging with the central MMS C-MM3. Radiative transfer modelling supports the idea that NGC 2264-C is a highly unstable prolate clump in the process of collapsing along its long axis on a near free-fall dynamical timescale ~ 1.7x10^5 yr. Our model fit of this large-scale collapse suggests a maximum mass inflow rate ~ 3x10^{-3} Msun/yr toward the central protostellar object C-MM3. Such infall rates are sufficiently high to overcome radiation pressure and allow the formation of ~ 20 Msun stars by accretion in ~ 1.7x10^5 yr, i.e., a time similar to the global dynamical timescale of the central part of NGC 2264-C. We conclude that we are likely witnessing the formation of a high-mass (> 10 Msun) protostar in the central part of NGC 2264-C. Our results suggest a picture of massive star formation intermediate between the scenario of stellar mergers of Bonnell et al. (1998) and the massive turbulent core model of McKee & Tan (2003).