Utilizing multiwavelength dust emission maps acquired with Herschel, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using an 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 were previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in the 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 ± 0.5 K and 9.8 ± 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 and 11 M , respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, κ λ = 0.1 × agrees to within 50% with the dust extinction constraints.

Context. For the project Galactic cold cores, Herschel photometric observations were carried out as a follow-up of cold regions of interstellar clouds previously identified with the Planck satellite. The aim of the project is to derive the physical properties of the population of cold sources and to study its connection to ongoing and future star formation. Aims. We build a catalogue of cold sources within the clouds in 116 fields observed with the Herschel PACS and SPIRE instruments. We wish to determine the general physical characteristics of the cold sources and to examine the correlations with their host cloud properties. Methods. From Herschel data, we computed colour temperature and column density maps of the fields. We estimated the distance to the target clouds and provide both uncertainties and reliability flags for the distances. The getsources multiwavelength source extraction algorithm was employed to build a catalogue of several thousand cold sources. Mid-infrared data were used, along with colour and position criteria, to separate starless and protostellar sources. We also propose another classification method based on submillimetre temperature profiles. We analysed the statistical distributions of the physical properties of the source samples. Results. We provide a catalogue of ∼4000 cold sources within or near star forming clouds, most of which are located either in nearby molecular complexes ( 1 kpc) or in star forming regions of the nearby galactic arms (∼2 kpc). About 70% of the sources have a size compatible with an individual core, and 35% of those sources are likely to be gravitationally bound. Significant statistical differences in physical properties are found between starless and protostellar sources, in column density versus dust temperature, mass versus size, and mass versus dust temperature diagrams. The core mass functions are very similar to those previously reported for other regions. On statistical grounds we find that gravitationally bound sources have higher background column densities (median N bg (H 2 ) ∼ 5 × 10 21 cm −2 ) than unbound sources (median N bg (H 2 ) ∼ 3 × 10 21 cm −2 ). These values of N bg (H 2 ) are higher for higher dust temperatures of the external layers of the parent cloud. However, only in a few cases do we find clear N bg (H 2 ) thresholds for the presence of cores. The dust temperatures of cloud external layers show clear variations with galactic location, as may the source temperatures. Conclusions. Our data support a more complex view of star formation than in the simple idea of a column density threshold. They show a clear influence of the surrounding UV-visible radiation on how cores distribute in their host clouds with possible variations on the Galactic scale.

Context. The Galactic Cold Cores project has carried out Herschel photometric observations of 116 fields where the Planck survey has found signs of cold dust emission. The fields contain sources in different environments and different phases of star formation. Previous studies have revealed variations in their dust submillimetre opacity. Aims. The aim is to measure the value of dust opacity spectral index and to understand its variations spatially and with respect to other parameters, such as temperature, column density, and Galactic location. Methods. The dust opacity spectral index β and the dust colour temperature T are derived using Herschel and Planck data. The relation between β and T is examined for the whole sample and inside individual fields. Results. Based on IRAS and Planck data, the fields are characterised by a median colour temperature of 16.1 K and a median opacity spectral index of β = 1.84. The values are not correlated with Galactic longitude. We observe a clear T -β anti-correlation. In Herschel observations, constrained at lower resolution by Planck data, the variations follow the column density structure and β FIR can rise to ∼2.2 in individual clumps. The highest values are found in starless clumps. The Planck 217 GHz band shows a systematic excess that is not restricted to cold clumps and is thus consistent with a general flattening of the dust emission spectrum at millimetre wavelengths. When fitted separately below and above 700 µm, the median spectral index values are β FIR ∼ 1.91 and β(mm) ∼ 1.66. Conclusions. The spectral index changes as a function of column density and wavelength. The comparison of different data sets and the examination of possible error sources show that our results are robust. However, β variations are partly masked by temperature gradients and the changes in the intrinsic grain properties may be even greater.

In this work we have carried out an in-depth analysis of the young stellar content in the W3 GMC. The YSO population was identified and classified in the IRAC/MIPS color-magnitude space according to the 'Class' scheme and compared to other classifications based on intrinsic properties. Class 0/I and II candidates were also compared to low/intermediate-mass pre-main-sequence stars selected through their colors and magnitudes in 2MASS. We find that a reliable color/magnitude selection of low-mass PMS stars in the infrared requires prior knowledge of the protostar population, while intermediate mass objects can be more reliably identified. By means of the MST algorithm and our YSO spatial distribution and age maps we investigated the YSO groups and the star formation history in W3. We find signatures of clustered and distributed star formation in both triggered and quiescent environments. The central/western parts of the GMC are dominated by large scale turbulence likely powered by isolated bursts of star formation that triggered secondary star formation events. Star formation in the eastern high density layer also shows signs of quiescent and triggered stellar activity, as well as extended periods of star formation. While our findings support triggering as a key factor for inducing and enhancing some of the major star forming activity in the HDL (e.g., W3 Main/W3(OH)), we argue that some degree of quiescent or spontaneous star formation is required to explain the observed YSO population. Our results also support previous studies claiming an spontaneous origin for the isolated massive star(s) powering KR 140.

We report the novel detection of complex high-column density tails in the probability distribution functions (PDFs) for three high-mass star-forming regions (CepOB3, MonR2, NGC6334), obtained from dust emission observed with Herschel. The low column density range can be fit with a lognormal distribution. A first power-law tail starts above an extinction (A V ) of ∼6-14. It has a slope of α=1.3-2 for the ρ ∝ r −α profile for an equivalent density distribution (spherical or cylindrical geometry), and is thus consistent with free-fall gravitational collapse. Above A V ∼40, 60, and 140, we detect an excess that can be fitted by a flatter power law tail with α >2. It correlates with the central regions of the cloud (ridges/hubs) of size ∼1 pc and densities above 10 4 cm −3 . This excess may be caused by physical processes that slow down collapse and reduce the flow of mass towards higher densities. Possible are: 1. rotation, which introduces an angular momentum barrier, 2. increasing optical depth and weaker cooling, 3. magnetic fields, 4. geometrical effects, and 5. protostellar feedback. The excess/second power-law tail is closely linked to high-mass star-formation though it does not imply a universal column density threshold for the formation of (high-mass) stars.

Aims. Ionization feedback should impact the probability distribution function (PDF) of the column density of cold dust around the ionized gas. We aim to quantify this effect and discuss its potential link to the core and initial mass function (CMF/IMF). Methods. We used Herschel column density maps of several regions observed within the HOBYS key program in a systematic way: M 16, the Rosette and Vela C molecular clouds, and the RCW 120 H ii region. We computed the PDFs in concentric disks around the main ionizing sources, determined their properties, and discuss the effect of ionization pressure on the distribution of the column density. Results. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a "double-peak" or an 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. Such a double peak is not visible for all clouds associated with ionization fronts, but it depends on the relative importance of ionization pressure and turbulent ram pressure. A power-law tail is present for higher column densities, which are generally ascribed to the effect of gravity. 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 that is able to disentangle triggered star formation from pre-existing star formation. Conclusions. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double-peaked or enlarged shape of the PDF may affect the formation of objects at both the low-mass and the high-mass ends of the CMF/IMF. In particular, a broader PDF is required by the gravo-turbulent scenario to fit the IMF properly 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 said to broaden the PDF, the relative importance of the different effects remains an open question.

We have studied the filaments extracted from the column density maps of the nearby Lupus 1, 3, and 4 molecular clouds, derived from photometric maps observed with the Herschel satellite. Filaments in the Lupus clouds have quite low column densities, with a median value of ∼1.5×10 21 cm −2 and most have masses per unit length lower than the maximum critical value for radial gravitational collapse. Indeed, no evidence of filament contraction has been seen in the gas kinematics. We find that some filaments, that on average are thermally subcritical, contain dense cores that may eventually form stars. This is an indication that in the low column density regime, the critical condition for the formation of stars may be reached only locally and this condition is not a global property of the filament. Finally, in Lupus we find multiple observational evidences of the key role that the magnetic field plays in forming filaments, and determining their confinement and dynamical evolution.

Aims. To constrain models of high-mass star formation, the Herschel-HOBYS key program aims at discovering massive dense cores (MDCs) able to host the high-mass analogs of low-mass prestellar cores, which have been searched for over the past decade. We here focus on NGC 6334, one of the best-studied HOBYS molecular cloud complexes. Methods. We used Herschel/PACS and SPIRE 70−500 µm images of the NGC 6334 complex complemented with (sub)millimeter and mid-infrared data. We built a complete procedure to extract ∼0.1 pc dense cores with the getsources software, which simultaneously measures their far-infrared to millimeter fluxes. We carefully estimated the temperatures and masses of these dense cores from their spectral energy distributions (SEDs). We also identified the densest pc-scale cloud structures of NGC 6334, one 2 pc × 1 pc ridge and two 0.8 pc × 0.8 pc hubs, with volume-averaged densities of ∼10 5 cm −3 . Results. A cross-correlation with high-mass star formation signposts suggests a mass threshold of 75 M for MDCs in NGC 6334. MDCs have temperatures of 9.5−40 K, masses of 75−1000 M , and densities of 1 × 10 5 −7 × 10 7 cm −3 . Their mid-infrared emission is used to separate 6 IR-bright and 10 IR-quiet protostellar MDCs while their 70 µm emission strength, with respect to fitted SEDs, helps identify 16 starless MDC candidates. The ability of the latter to host high-mass prestellar cores is investigated here and remains questionable. An increase in mass and density from the starless to the IR-quiet and IR-bright phases suggests that the protostars and MDCs simultaneously grow in mass. The statistical lifetimes of the high-mass prestellar and protostellar core phases, estimated to be 1−7 × 10 4 yr and at most 3 × 10 5 yr respectively, suggest a dynamical scenario of high-mass star formation. Conclusions. The present study provides good mass estimates for a statistically significant sample, covering the earliest phases of high-mass star formation. High-mass prestellar cores may not exist in NGC 6334, favoring a scenario presented here, which simultaneously forms clouds, ridges, MDCs, and high-mass protostars.