The far-infrared (FIR) regime is one of the few wavelength ranges where no astronomical data with sub-arcsecond spatial resolution exist. Neither of the medium-term satellite projects like SPICA, Millimetron nor O.S.T. will resolve this malady. For many research areas, however, information at high spatial and spectral resolution in the FIR, taken from atomic fine-structure lines, from highly excited carbon monoxide (CO), light hydrids, and especially from water lines would open the door for transformative science. A main theme will be to trace the role of water in proto-planetary disks, to observationally advance our understanding of the planet formation process and, intimately related to that, the pathways to habitable planets and the emergence of life. Furthermore, key observations will zoom into the physics and chemistry of the star-formation process in our own Galaxy, as well as in external galaxies. The FIR provides unique tools to investigate in particular the energetics of heating, cooling and shocks. The velocity-resolved data in these tracers will reveal the detailed dynamics engrained in these processes in a spatially resolved fashion, and will deliver the perfect synergy with ground-based molecular line data for the colder dense gas.

The Methanol MultiBeam survey (MMB) provides the most complete sample of Galactic massive young stellar objects (MYSOs) hosting 6.7GHz class II methanol masers. We characterise the properties of these maser sources using dust emission detected by the Herschel Infrared Galactic Plane Survey (Hi-GAL) to assess their evolutionary state. Associating 731 (73%) of MMB sources with compact emission at four Hi-GAL wavelengths, we derive clump properties and define the requirements of a MYSO to host a 6.7GHz maser. The median far-infrared (FIR) mass and luminosity are 630M$_{\odot}$ and 2500L$_{\odot}$ for sources on the near side of Galactic centre and 3200M$_{\odot}$ and 10000L$_{\odot}$ for more distant sources. The median luminosity-to-mass ratio is similar for both at $\sim$4.2L$_{\odot}/$M$_{\odot}$. We identify an apparent minimum 70$\mu$m luminosity required to sustain a methanol maser of a given luminosity (with $L_{70} \propto L_{6.7}^{0.6}$). The maser host clumps have higher mass and higher FIR luminosities than the general Galactic population of protostellar MYSOs. Using principal component analysis, we find 896 protostellar clumps satisfy the requirements to host a methanol maser but lack a detection in the MMB. Finding a 70$\mu$m flux density deficiency in these objects, we favour the scenario in which these objects are evolved beyond the age where a luminous 6.7GHz maser can be sustained. Finally, segregation by association with secondary maser species identifies evolutionary differences within the population of 6.7GHz sources.

We present the Forgotten Quadrant Survey (FQS), an ESO large project that used the 12m antenna of the Arizona Radio Observatory to map the Galactic Plane in the range 220\deg$<l<$240\deg and -2.5\deg$<b<$0\deg, both in $^{12}$CO(1-0) and $^{13}$CO(1-0), at a spectral resolution of 0.65 km s$^{-1}$ and 0.26 km s$^{-1}$. Our dataset allows us to easily identify how the molecular dense gas is organised at different spatial scales: from the giant clouds with their denser filamentary networks, down to the clumps and cores that host the newborn stars and to obtain reliable estimates of their key physical parameters. We present the first release of the FQS data and discuss their quality. Spectra with 0.65 km s$^{-1}$ velocity channels have a noise ranging from 0.8 K to 1.3 K for $^{12}$CO(1-0) and from 0.3 K to 0.6 K for $^{13}$CO(1-0). In this paper, we used the $^{12}$CO(1-0) spectral cubes to produce a catalogue of 263 molecular clouds. This is the first selfconsistent, statistical catalogue of molecular clouds of the outer Galaxy, obtained with a subarcminute spatial resolution and therefore able to detect not only the classical giant molecular clouds, but also the small clouds and to resolve the cloud structure at the subparsec scale up to a distance of a few kpc. We found two classes of objects: structures with size above a few parsecs that are typical molecular clouds and may be self-gravitating, and subparsec structures that cannot be in gravitational equilibrium and are likely transient or confined by external pressure. We used the ratio between the Herschel H$_2$ column density and the integrated intensity of the CO lines to calculate the CO conversion factor and we found mean values of (3.3$\pm$1.4)$\times 10^{20}$ cm$^{-2}$(K km s$^{-1})^{-1}$ and (1.2$\pm$0.4)$\times 10^{21}$ cm$^{-2}$(K km s$^{-1})^{-1}$, for $^{12}$CO(1-0) and $^{13}$CO(1-0), respectively.

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.

In this study we analyze 70 radio continuum sources associated with dust clumps and considered to be candidates for the earliest stages of high-mass star formation. The detection of these sources was reported by Rosero et al. (2016), who found most of them to show weak (${\scriptstyle <}$1 mJy) and compact (${\scriptstyle <}\,$0.6$^{\prime \prime}$) radio emission. Herein, we used the observed parameters of these sources to investigate the origin of the radio continuum emission. We found that at least $\sim 30\%$ of these radio detections are most likely ionized jets associated with high-mass protostars, but for the most compact sources we cannot discard the scenario that they represent pressure-confined HII regions. This result is highly relevant for recent theoretical models based on core accretion that predict the first stages of ionization from high-mass stars to be in the form of jets. Additionally, we found that properties such as the radio luminosity as a function of the bolometric luminosity of ionized jets from low and high-mass stars are extremely well-correlated. Our data improve upon previous studies by providing further evidence of a common origin for jets independently of luminosity.

A grid of 20 millions 3-1100$\mu$m SED models is presented for synthetic young clusters embedded in dense clumps. The models depend on four primary parameters: the clump mass M$_{clump}$ and dust temperature T$_{dust}$, the fraction of mass f$_{core}$ locked in dense cores, and the age of the clump t$_{SF}$. We populate the YSO clusters using the IMF from Kroupa(2001) and the YSOs SED models grid of Robitaille et al. (2006). We conduct extensive testing of SED fitting using a simulated dataset and we find that M$_{clump}$ essentially depends on the submillimeter portion of the SED, while T$_{dust}$ is mostly determined from the shape of the SED in the 70-350$\mu$m range. Thanks to the large number of models computed we verify that the combined analysis of L/M, [8-24] and [24-70] colours removes much of the SEDs f$_{core}$-t$_{SF}$ degeneracy. The L/M values are particularly useful to diagnose f$_{core}$. L/M$\leq$1 identifies protoclusters with f$_{core}\leq$0.1 and t$_{SF} \leq 10^5$ years, while L/M$\geq$10 excludes f$_{core}\leq$0.1. We characterize lower limits of L/M where ZAMS stars are not found in models, and we also find models with L/M $\geq$10 and no ZAMS stars, in which [8-24]$\geq0.8\pm 0.1$ independently from M$_{clump}$, temperature and luminosity. This is the first set of synthesis SED models suited to model for embedded and unresolved clusters of YSOs. A set of new evolutionary tracks in the L/M diagram is also presented.

We present a morphological and physical analysis of a Giant Molecular Cloud (GMC) using the carbon monoxide isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $^{3}P_{2}\rightarrow$ $^{3}P_{1}$) survey of the Galactic Plane (Mopra CO Southern Galactic Plane Survey), supplemented with neutral carbon maps from the HEAT telescope in Antarctica. The giant molecular cloud structure (hereinafter the ring) covers the sky region $332^\circ$ < $\ell$ < $333^\circ$ and $\mathit{b}$ = $\pm 0.5^\circ$ (hereinafter the G332 region). The mass of the ring and its distance are determined to be respectively ~2$\times10^{5}\mathrm{M_{\odot}}$ and ~3.7 kpc from Sun. The dark molecular gas fraction, estimated from the $^{13}$CO and [CI] lines, is $\sim17\%$ for a CO T$_{\mathrm{ex}}$ between [10,20 K]. Comparing the [CI] integrated intensity and N(H$_{2}$) traced by $^{13}$CO and $^{12}$CO, we define an X$\mathrm{_{CI}^{809}}$ factor, analogous to the usual X$_{\mathrm{co}}$, through the [CI] line. X$\mathrm{_{CI}^{809}}$ ranges between [1.8,2.0]$\times10^{21}\mathrm{cm}^{-2}\mathrm{K}^{-1}\mathrm{km}^{-1}\mathrm{s}$. We examined local variation in X$_{\mathrm{co}}$ and T$_{\mathrm{ex}}$ across the cloud, and find in regions where the star formation activity is not in an advanced state, an increase in the mean and dispersion of the X$_{\mathrm{co}}$ factor as the excitation temperature decreases. We present a catalogue of C$^{18}$O clumps within the cloud. The star formation (SF) activity ongoing in the cloud shows a correlation with T$_{\mathrm{ex}}$, [CI] and CO emissions, and anti-correlation with X$_{\mathrm{co}}$, suggesting a North-South spatial gradient in the SF activity. We propose a method to disentangle dust emission across the Galaxy, using HI and $^{13}$CO data. We describe Virtual Reality (VR) and Augmented Reality (AR) data visualisation techniques for the analysis of radio astronomy data.

Context. In recent years, an in-depth gamma-ray analysis of the Orion region has been carried out by the AGILE and Fermi-LAT (Large Area Telescope) teams with the aim of estimating the H2-CO conversion factor, XCO. The comparison of the data from both satellites with models of diffuse gamma-ray Galactic emission unveiled an excess at (l,b)=[213.9, -19.5], in a region at a short angular distance from the OB star k-Ori. Possible explanations of this excess are scattering of the so-called "dark gas", non-linearity in the H2-CO relation, or Cosmic-Ray (CR) energization at the k-Ori wind shock. Aims. Concerning this last hypothesis, we want to verify whether cosmic-ray acceleration or re-acceleration could be triggered at the k-Ori forward shock, which we suppose to be interacting with a star-forming shell detected in several wavebands and probably triggered by high energy particles. Methods. Starting from the AGILE spectrum of the detected gamma-ray excess, showed here for the first time, we developed a valid physical model for cosmic-ray energization, taking into account re-acceleration, acceleration, energy losses, and secondary electron contribution. Results. Despite the characteristic low velocity of an OB star forward shock during its "snowplow" expansion phase, we find that the Orion gamma-ray excess could be explained by re-acceleration of pre-existing cosmic rays in the interaction between the forward shock of k-Ori and the CO-detected, star-forming shell swept-up by the star expansion. According to our calculations, a possible contribution from freshly accelerated particles is sub-dominant with respect the re-acceleration contribution. However, a simple adiabatic compression of the shell could also explain the detected gamma-ray emission. Futher GeV and TeV observations of this region are highly recommended in order to correctly identify the real physical scenario.

We present a comparative study of physical properties derived from gas and dust emission in a sample of 1068 dense Galactic clumps. The sources are selected from the crossmatch of the Hi-GAL survey with 16 catalogues of NH$_3$ line emission in its lowest inversion (1,1) and (2,2) transitions. The sample covers a large range in masses and bolometric luminosities, with surface densities above $\Sigma=0.1$ g cm$^{-2}$ and with low virial parameters $\alpha<1$. The comparison between dust and gas properties shows an overall agreement between $T_{\textit{kin}}$ and $T_{\textit{dust}}$ at volumetric densities $n\gtrsim1.2\times10^{4}$ cm$^{-3}$, and a median fractional abundance $\chi$(NH$_3$)$=1.46\times10^{-8}$. While the protostellar clumps in the sample have small differences between $T_{\textit{kin}}$ and $T_{\textit{dust}}$, prestellar clumps have a median ratio $T_{\textit{kin}}/T_{\textit{dust}}=1.24$, suggesting that these sources are thermally decoupled. A correlation is found between the evolutionary tracer $L/M$ and the parameters $T_{\textit{kin}}/T_{\textit{dust}}$ and $\chi$(NH$_3$) in prestellar sources and protostellar clumps with $L/M<1$ L$_\odot$ M$_\odot^{-1}$. In addition, a weak correlation is found between non-thermal velocity dispersion and the $L/M$ parameter, possibly indicating an increase of turbulence with protostellar evolution in the interior of clumps. Finally, different processes are discussed to explain the differences between gas and dust temperatures in prestellar candidates, and the origin of non-thermal motions observed in the clumps.

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

The dynamics of massive clumps, the environment where massive stars originate, is still unclear. Many theories predict that these regions are in a state of near-virial equilibrium, or near energy equi-partition, while others predict that clumps are in a sub-virial state. Observationally, the majority of the massive clumps are in a sub-virial state with a clear anti-correlation between the virial parameter $\alpha_{vir}$ and the mass of the clumps $M_{c}$, which suggests that the more massive objects are also the more gravitationally bound. Although this trend is observed at all scales, from massive clouds down to star-forming cores, theories do not predict it. In this work we show how, starting from virialized clumps, an observational bias is introduced in the specific case where the kinetic and the gravitational energies are estimated in different volumes within clumps and how it can contribute to the spurious $\alpha_{vir}-M_{c}$ anti-correlation in these data. As a result, the observed effective virial parameter $\tilde{\alpha}_{eff}<\alpha_{vir}$, and in some circumstances it might not be representative of the virial state of the observed clumps.

The multifractal geometry remains an under-exploited approach to describe and quantify the large-scale structure of interstellar clouds. In this paper, the typical tools of multifractal analysis are applied to Herschel far-infrared (70-500 $\mu$m) dust continuum maps, which represent an ideal case of study. This dust component is a relatively optically thin tracer at these wavelengths and the size in pixel of the maps is well suitable for this statistical analysis. We investigate the so-called multifractal spectrum and generalised fractal dimension of six Hi-GAL regions in the third Galactic quadrant. We measure systematic variations of the spectrum at increasing wavelength, which generally correspond to an increasing complexity of the image, and we observe peculiar behaviours of the investigated fields, strictly related to the presence of high-emission regions, which in turn are connected to star formation activity. The same analysis is applied to synthetic column density maps, generated from numerical turbulent molecular cloud models and from fractal Brownian motion (fBm), allowing for the confrontation of the observations with models with well controlled physical parameters. This comparison shows that, in terms of multifractal descriptors of the structure, fBm images exhibit a very different, quite unrealistic behaviour when compared with Hi-GAL observations, whereas the numerical simulations appear in some cases (depending on the specific model) more similar to the observations. Finally, the link between mono-fractal parameters (commonly found in the literature) and multifractal indicators is investigated: the former appear to be only partially connected with the latter, with the multifractal analysis offering a more extended and complete characterization of the cloud structure.

We present a visual analytics tool, based on the VisIVO suite, to exploit a combination of all new-generation surveys of the Galactic Plane to study the star formation process of the Milky Way. The tool has been developed within the VIALACTEA project, founded by the 7th Framework Programme of the European Union, that creates a common forum for the major new-generation surveys of the Milky Way Galactic Plane from the near infrared to the radio, both in thermal continuum and molecular lines. Massive volumes of data are produced by space missions and ground-based facilities and the ability to collect and store them is increasing at a higher pace than the ability to analyze them. This gap leads to new challenges in the analysis pipeline to discover information contained in the data. Visual analytics focuses on handling these massive, heterogeneous, and dynamic volumes of information accessing the data previously processed by data mining algorithms and advanced analysis techniques with highly interactive visual interfaces offering scientists the opportunity for in-depth understanding of massive, noisy, and high-dimensional data.

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

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.

Diffuse galactic gamma-ray emission is produced by the interaction of cosmic rays (CRs) with the interstellar environment. The study of gamma-ray emission is therefore a powerful tool to investigate the origin of CRs and the processes through which they are accelerated. We aim to gain deeper insights of the nature of gamma-ray emission in the region of Orion, which is one of the best studied sites of on-going star formation, by analysing data from the AGILE satellite. The diffuse gamma-ray emission expected from the Orion region is relatively high. Its separation from the galactic plane also ensures a very low contribution from foreground or background emission, which makes it an ideal site for studying the processes of particle acceleration in star forming environments. The AGILE data are modelled through a template that quantifies the gamma-ray diffuse emission expected from atomic and molecular hydrogen. Other sources of emission are modelled as an isotropic contribution. Gamma-ray emission exceeding the amount expected by the diffuse emission model is detected with high level of significance. The main excess is in the high-longitude part of Orion A. A thorough analysis of this feature suggests a connection between the observed gamma-ray emission and the B0.5 Ia star k Orionis. The location of the gamma-ray excess is compatible with the site where stellar wind collides with the ISM. Both scattering on dark gas and cosmic-ray acceleration at the shock between the two environments are discussed as possible explanations, with the latter hypothesis being supported by the hardness of the energy spectrum of the emission. If confirmed, this would be the first direct detection of gamma-ray emission from the interaction between ISM and a single star's stellar wind.

We present a catalogue of prestellar and starless cores within the Corona Australis molecular cloud using photometric data from the Herschel Space Observatory. At a distance of d~130 pc, Corona Australis is one of the closest star-forming regions. Herschel has taken multi-wavelength data of Corona Australis with both the SPIRE and PACS photometric cameras in a parallel mode with wavelengths in the range 70 {\mu}m to 500 {\mu}m. A complete sample of starless and prestellar cores and embedded protostars is identified. Other results from the Herschel Gould Belt Survey have shown spatial correlation between the distribution of dense cores and the filamentary structure within the molecular clouds. We go further and show correlations between the properties of these cores and their spatial distribution within the clouds, with a particular focus on the mass distribution of the dense cores with respect to their filamentary proximity. We find that only lower-mass starless cores form away from filaments, while all of the higher-mass prestellar cores form in close proximity to, or directly on the filamentary structure. This result supports the paradigm that prestellar cores mostly form on filaments. We analyse the mass distribution across the molecular cloud, finding evidence that the region around the Coronet appears to be at a more dynamically advanced evolutionary stage to the rest of the clumps within the cloud.

In this paper, we present the first extended catalogue of far-infrared fluxes of Galactic bubbles. Fluxes were estimated for 1814 bubbles, defined here as the `golden sample', and were selected from the Milky Way Project First Data Release (Simpson et al.) The golden sample was comprised of bubbles identified within the Wide-field Infrared Survey Explorer (WISE) dataset (using 12- and 22-$\mu$m images) and Herschel data (using 70-, 160-, 250-, 350- and 500-$\mu$m wavelength images). Flux estimation was achieved initially via classical aperture photometry and then by an alternative image analysis algorithm that used active contours. The accuracy of the two methods was tested by comparing the estimated fluxes for a sample of bubbles, made up of 126 H II regions and 43 planetary nebulae, which were identified by Anderson et al. The results of this paper demonstrate that a good agreement between the two was found. This is by far the largest and most homogeneous catalogue of infrared fluxes measured for Galactic bubbles and it is a step towards the fully automated analysis of astronomical datasets.

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.

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

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

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 a comprehensive statistical analysis of star-forming objects located in the vicinities of 1 360 bubble structures throughout the Galactic Plane and their local environments. The compilation of ~70 000 star-forming sources, found in the proximity of the ionized (Hii) regions and detected in both Hi-GAL and GLIMPSE surveys, provided a broad overview of the different evolutionary stages of star-formation in bubbles, from prestellar objects to more evolved young stellar objects (YSOs). Surface density maps of star-forming objects clearly reveal an evolutionary trend where more evolved star-forming objects are found spatially located near the center, while younger star-forming objects are found at the edge of the bubbles. We derived dynamic ages for a subsample of 182 Hii regions for which kinematic distances and radio continuum flux measurements were available. We detect ~80% more star-forming sources per unit area in the direction of bubbles than in the surrounding fields. We estimate ~10% clump formation efficiency (CFE) of Hi-GAL clumps in bubbles, twice the CFE in fields not affected by feedback. We find higher CFE of protostellar clumps in younger bubbles, whose density of the bubble shells is higher. We argue that the formation rate from prestellar to protostellar phase is probably higher during the early stages of the bubble expansion. Evaluation of the fragmentation time inside the shell of bubbles advocates the preexistence of clumps in the medium before the bubble, as supported by numerical simulations. Approximately 23% of the Hi-GAL clumps are found located in the direction of a bubble, with 15% for prestellar clumps and 41% for protostellar clumps. We argue that the high fraction of protostellar clumps may be due to the acceleration of the star-formation process cause by the feedback of the (Hii) bubbles.

We report 3 mm continuum, CH3CN(5-4) and 13CS(2-1) line observations with CARMA, in conjunction with 6 and 1.3 cm continuum VLA data, and 12 and 25 micron broadband data from the Subaru Telescope toward the massive proto-star IRAS18566+0408. The VLA data resolve the ionized jet into 4 components aligned in the E-W direction. Radio components A, C, and D have flat cm SEDs indicative of optically thin emission from ionized gas, and component B has a spectral index alpha = 1.0, and a decreasing size with frequency proportional to frequency to the -0.5 power. Emission from the CARMA 3 mm continuum, and from the 13CS(2-1), and CH3CN(5-4) spectral lines is compact (i.e. < 6700 AU), and peaks near the position of VLA cm source, component B. Analysis of these lines indicates hot, and dense molecular gas, typical for HMCs. Our Subaru telescope observations detect a single compact source, coincident with radio component B, demonstrating that most of the energy in IRAS18566+0408 originates from a region of size < 2400 AU. We also present UKIRT near-infrared archival data for IRAS18566+0408 which show extended K-band emission along the jet direction. We detect an E-W velocity shift of about 10 km/sec over the HMC in the CH3CN lines possibly tracing the interface of the ionized jet with the surrounding core gas. Our data demonstrate the presence of an ionized jet at the base of the molecular outflow, and support the hypothesis that massive protostars with O-type luminosity form with a mechanism similar to lower mass stars.

To constrain models of high-mass star formation, the Herschel/HOBYS KP 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 NGC6334, one of the best-studied HOBYS molecular cloud complexes. We used Herschel PACS and SPIRE 70-500mu images of the NGC6334 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 SEDs. A cross-correlation with high-mass star formation signposts suggests a mass threshold of 75Msun for MDCs in NGC6334. MDCs have temperatures of 9.5-40K, masses of 75-1000Msun, and densities of 10^5-10^8cm-3. Their mid-IR emission is used to separate 6 IR-bright and 10 IR-quiet protostellar MDCs while their 70mu 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-7x10^4yr and at most 3x10^5yr respectively, suggest a dynamical scenario of high-mass star formation. 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 NGC6334, favoring a scenario presented here, which simultaneously forms clouds and high-mass protostars.

In the current data-driven science era, it is needed that data analysis techniques has to quickly evolve to face with data whose dimensions has increased up to the Petabyte scale. In particular, being modern astrophysics based on multi-wavelength data organized into large catalogues, it is crucial that the astronomical catalog cross-matching methods, strongly dependant from the catalogues size, must ensure efficiency, reliability and scalability. Furthermore, multi-band data are archived and reduced in different ways, so that the resulting catalogues may differ each other in formats, resolution, data structure, etc, thus requiring the highest generality of cross-matching features. We present $C^{3}$ (Command-line Catalogue Cross-match), a multi-platform application designed to efficiently cross-match massive catalogues from modern surveys. Conceived as a stand-alone command-line process or a module within generic data reduction/analysis pipeline, it provides the maximum flexibility, in terms of portability, configuration, coordinates and cross-matching types, ensuring high performance capabilities by using a multi-core parallel processing paradigm and a sky partitioning algorithm.

The degradation of spatial resolution in star-forming regions observed at large distances ($d\gtrsim1$ kpc) with Herschel,can lead to estimates of the physical parameters of the detected compact sources (clumps) which do not necessarily mirror the properties of the original population of cores. This paper aims at quantifying the bias introduced in the estimation of these parameters by the distance effect. To do so, we consider Herschel maps of nearby star-forming regions taken from the Herschel-Gould-Belt survey, and simulate the effect of increased distance to understand what amount of information is lost when a distant star-forming region is observed with Herschel resolution. In the maps displaced to different distances we extract compact sources, and we derive their physical parameters as if they were original Hi-GAL maps of the extracted source samples. In this way, we are able to discuss how the main physical properties change with distance. In particular, we discuss the ability of clumps to form massive stars: we estimate the fraction of distant sources that are classified as high-mass stars-forming objects due to their position in the mass vs radius diagram, that are only "false positives". We give also a threshold for high-mass star-formation $M>1282 \ \left(\frac{r}{[\mathrm{pc}]}\right)^{1.42} M_{\odot}$. In conclusion, this paper provides the astronomer dealing with Herschel maps of distant star-forming regions with a set of prescriptions to partially recover the character of the core population in unresolved clumps.

We present the physical and evolutionary properties of prestellar and protostellar clumps in the Herschel Infrared GALactic plane survey (Hi-GAL) in two large areas centered in the Galactic plane and covering the tips of the long Galactic bar at the intersection with the spiral arms. The areas fall in the longitude ranges 19 < l < 33 and 340 < l < 350, while latitude is -1 < b < 1. Newly formed high mass stars and prestellar objects are identified and their properties derived and compared. A study is also presented on five giant molecular complexes at the further edge of the bar. The star-formation rate was estimated from the quantity of proto-stars expected to form during the collapse of massive turbulent clumps into star clusters. This new method was developed by applying a Monte Carlo procedure to an evolutionary model of turbulent cores and takes into account the wide multiplicity of sources produced during the collapse. The star-formation rate density values at the tips are 1.2 +- 0.3 10-3 Msol/yr/kpc2 and 1.5+-0.3 10-3 Msol/yr/kpc2 in the first and fourth quadrant, respectively. The same values estimated on the entire field of view, that is including the tips of the bar and background and foreground regions, are 0.9+-0.2 10-3 Msol/yr/kpc2 and 0.8+-0.2 10-3 Msol/yr/kpc2. The conversion efficiency is approximately 0.8% in the first quadrant and 0.5% in the fourth quadrant, and does not show a significant difference in proximity of the bar. The star forming regions identified through CO contours at the further edge of the bar show star-formation rate densities larger than the surrounding regions but their conversion efficiencies are comparable. Our results suggest that the star-formation activity at the bar is due to a large amount of dust and molecular material rather than being due to a triggering process.

We present a multi-wavelength study of the Galactic Luminous Blue Variable HR Carinae, based on new high resolution mid-infrared (IR) and radio images obtained with the Very Large Telescope (VLT) and the Australia Telescope Compact Array (ATCA), which have been complemented by far-infrared Herschel-PACS observations and ATCA archive data. The Herschel images reveal the large-scale distribution of the dusty emitting nebula, which extends mainly to the North-East direction, up to 70 arcsec from the central star, and is oriented along the direction of the space motion of the star. In the mid-infrared images, the brightness distribution is characterized by two arc-shaped structures, tracing an inner envelope surrounding the central star more closely. At radio wavelengths, the ionized gas emission lies on the opposite side of the cold dust with respect to the position of the star, as if the ionized front was confined by the surrounding medium in the North-South direction. The comparison with previous data indicates significant changes in the radio nebula morphology and in the mass-loss rate from the central star, which has increased from 6.1$\times10^{-6}$ M$_{\odot}$yr$^{-1}$ in 1994-1995 to $1.17\times10^{-5}$ M$_{\odot}$yr$^{-1}$ in 2014. We investigate possible scenarios that could have generated the complex circumstellar environment revealed by our multi-wavelength data.

Modern Astrophysics is based on multi-wavelength data organized into large and heterogeneous catalogues. Hence, the need for efficient, reliable and scalable catalogue cross-matching methods plays a crucial role in the era of the petabyte scale. Furthermore, multi-band data have often very different angular resolution, requiring the highest generality of cross-matching features, mainly in terms of region shape and resolution. In this work we present $C^{3}$ (Command-line Catalogue Cross-match), a multi-platform application designed to efficiently cross-match massive catalogues. It is based on a multi-core parallel processing paradigm and conceived to be executed as a stand-alone command-line process or integrated within any generic data reduction/analysis pipeline, providing the maximum flexibility to the end-user, in terms of portability, parameter configuration, catalogue formats, angular resolution, region shapes, coordinate units and cross-matching types. Using real data, extracted from public surveys, we discuss the cross-matching capabilities and computing time efficiency also through a direct comparison with some publicly available tools, chosen among the most used within the community, and representative of different interface paradigms. We verified that the $C^{3}$ tool has excellent capabilities to perform an efficient and reliable cross-matching between large datasets. Although the elliptical cross-match and the parametric handling of angular orientation and offset are known concepts in the astrophysical context, their availability in the presented command-line tool makes $C^{3}$ competitive in the context of public astronomical tools.

The emerging need for efficient, reliable and scalable astronomical catalog cross-matching is becoming more pressing in the current data-driven science era, where the size of data has rapidly increased up to the Petabyte scale. C3 (Command-line Catalogue Cross-matching) is a multi-platform tool designed to efficiently cross-match massive catalogues from modern astronomical surveys, ensuring high-performance capabilities through the use of a multi-core parallel processing paradigm. The tool has been conceived to be executed as a stand-alone command-line process or integrated within any generic data reduction/analysis pipeline, providing the maximum flexibility to the end user, in terms of parameter configuration, coordinates and cross-matching types. In this work we present the architecture and the features of the tool. Moreover, since the modular design of the tool enables an easy customization to specific use cases and requirements, we present also an example of a customized C3 version designed and used in the FP7 project ViaLactea, dedicated to cross-correlate Hi-GAL clumps with multi-band compact sources.

IRAS 18511+0146 is a young embedded (proto)cluster located at 3.5 kpc surrounding what appears to be an intermediate mass protostar. In this paper, we investigate the nature of cluster members (two of which are believed to be the most massive and luminous) using imaging and spectroscopy in the near and mid-infrared. The brightest point-like object associated with IRAS 18511+0146 is referred to as S7 in the present work (designated UGPS J185337.88+015030.5 in the UKIRT Galactic Plane survey). Seven of the nine objects show rising spectral energy distributions (SED) in the near-infrared, with four objects showing Br-gamma emission. Three members: S7, S10 (also UGPS J185338.37+015015.3) and S11 (also UGPS J185338.72+015013.5) are bright in mid-infrared with diffuse emission being detected in the vicinity of S11 in PAH bands. Silicate absorption is detected towards these three objects, with an absorption maximum between 9.6 and 9.7 um, large optical depths (1.8-3.2), and profile widths of 1.6-2.1 um. The silicate profiles of S7 and S10 are similar, in contrast to S11 (which has the largest width and optical depth). The cold dust emission investigated using Herschel HiGal peaks at S7, with temperature at 26 K and column density N(H2) ~ 7 x 10^(22) cm^(-2). The bolometric luminosity of IRAS 18511 region is L ~ 1.8 x 10^4 L_sun. S7 is the main contributor to the bolometric luminosity, with L (S7) > 10^4 L_sun. S7 is a high mass protostellar object with ionised stellar winds, evident from the correlation between radio and bolometric luminosity as well as the asymmetric Br-gamma profile. The differences in silicate profiles of S7 and S11 could be due to different radiation environment as we believe the former to be more massive and in an earlier phase than the latter.

The aims are to investigate the clustering of the far-infrared sources from the Herschel infrared Galactic Plane Survey (Hi-GAL) in the Galactic longitude range of -71 to 67 deg. These clumps, and their spatial distribution, are an imprint of the original conditions within a molecular cloud. This will produce a catalogue of over-densities. The minimum spanning tree (MST) method was used to identify the over-densities in two dimensions. The catalogue was further refined by folding in heliocentric distances, resulting in more reliable over-densities, which are cluster candidates. We found 1,633 over-densities with more than ten members. Of these, 496 are defined as cluster candidates because of the reliability of the distances, with a further 1,137 potential cluster candidates. The spatial distributions of the cluster candidates are different in the first and fourth quadrants, with all clusters following the spiral structure of the Milky Way. The cluster candidates are fractal. The clump mass functions of the clustered and isolated are statistically indistinguishable from each other and are consistent with Kroupa's initial mass function.

We present an innovative method called FilExSeC (Filaments Extraction, Selection and Classification), a data mining tool developed to investigate the possibility to refine and optimize the shape reconstruction of filamentary structures detected with a consolidated method based on the flux derivative analysis, through the column-density maps computed from Herschel infrared Galactic Plane Survey (Hi-GAL) observations of the Galactic plane. The present methodology is based on a feature extraction module followed by a machine learning model (Random Forest) dedicated to select features and to classify the pixels of the input images. From tests on both simulations and real observations the method appears reliable and robust with respect to the variability of shape and distribution of filaments. In the cases of highly defined filament structures, the presented method is able to bridge the gaps among the detected fragments, thus improving their shape reconstruction. From a preliminary "a posteriori" analysis of derived filament physical parameters, the method appears potentially able to add a sufficient contribution to complete and refine the filament reconstruction.

The VIALACTEA project has a work package dedicated to Tools and Infrastructure and, inside it, a task for the Database and Virtual Observatory Infrastructure. This task aims at providing an infrastructure to store all the resources needed by the, more purposely, scientific work packages of the project itself. This infrastructure includes a combination of: storage facilities, relational databases and web services on top of them, and has taken, as a whole, the name of VIALACTEA Knowledge Base (VLKB). This contribution illustrates the current status of this VLKB. It details the set of data resources put together; describes the database that allows data discovery through VO inspired metadata maintenance; illustrates the discovery, cutout and access services built on top of the former two for the users to exploit the data content.

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

The evolutionary classification of massive clumps that are candidate progenitors of high-mass young stars and clusters relies on a variety of independent diagnostics based on observables from the near-infrared to the radio. A promising evolutionary indicator for massive and dense cluster-progenitor clumps is the L/M ratio between the bolometric luminosity and the mass of the clumps. With the aim of providing a quantitative calibration for this indicator we used SEPIA/APEX to obtain CH3C2H(12-11) observations, that is an excellent thermometer molecule probing densities > 10^5 cm^-3 , toward 51 dense clumps with M>1000 solar masses, and uniformly spanning -2 < Log(L/M) < 2.3. We identify three distinct ranges of L/M that can be associated to three distinct phases of star formation in massive clumps. For L/M <1 no clump is detected in CH3C2H , suggesting an inner envelope temperature below 30K. For 1< L/M < 10 we detect 58% of the clumps, with a temperature between 30 and 35 K independently from the exact value of L/M; such clumps are building up luminosity due to the formation of stars, but no star is yet able to significantly heat the inner clump regions. For L/M> 10 we detect all the clumps, with a gas temperature rising with Log(L/M), marking the appearance of a qualitatively different heating source within the clumps; such values are found towards clumps with UCHII counterparts, suggesting that the quantitative difference in T - L/M behaviour above L/M >10 is due to the first appearance of ZAMS stars in the clumps.

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

Previous studies of the initial conditions of massive star formation have mainly targeted Infrared-Dark Clouds (IRDCs) toward the inner Galaxy. This is due to the fact that IRDCs were first detected in absorption against the bright mid-IR background, requiring a favourable location to be observed. By selection, IRDCs represent only a fraction of the Galactic clouds capable of forming massive stars and star clusters. Due to their low dust temperatures, IRDCs are bright in the far-IR and millimeter and thus, observations at these wavelengths have the potential to provide a complete sample of star-forming massive clouds across the Galaxy. Our aim is to identify the clouds at the initial conditions of massive star formation across the Galaxy and compare their physical properties as a function of their Galactic location. We have examined the physical properties of a homogeneous galactic cold core sample obtained with the Planck satellite across the Galactic Plane. With the use of Herschel Hi-GAL observations, we have characterized the internal structure of them. By using background-subtracted Herschel images, we have derived the H2 column density and dust temperature maps for 48 Planck clumps. Their basic physical parameters have been calculated and analyzed as a function of location within the Galaxy. These properties have also been compared with the empirical relation for massive star formation derived by Kauffmann & Pillai (2010). Most of the Planck clumps contain signs of star formation. About 25% of them are massive enough to form high mass stars. Planck clumps toward the Galactic center region show higher peak column densities and higher average dust temperatures than those of the clumps in the outer Galaxy. Although we only have seven clumps without associated YSOs, the Hi-GAL data show no apparent differences in the properties of Planck cold clumps with and without star formation.

Ionized regions around early-type stars are believed to be well-known objects, but until recently, our knowledge of the relation between the free-free radio emission and the IR emission has been observationally hindered by the limited angular resolution in the far-IR. The advent of Herschel has now made it possible to obtain a more precise comparison between the two regimes, and it has been found that about a third of the young HII regions emit more Lyman continuum photons than expected, thus presenting a Lyman excess. With the present study we wish to distinguish between two scenarios that have been proposed to explain the existence of the Lyman excess: (i) underestimation of the bolometric luminosity, or (ii) additional emission of Lyman-continuum photons from an accretion shock. We observed an outflow (SiO) and an infall (HCO+) tracer toward a complete sample of 200 HII regions, 67 of which present the Lyman excess. Our goal was to search for any systematic difference between sources with Lyman excess and those without. While the outflow tracer does not reveal any significant difference between the two subsamples of HII regions, the infall tracer indicates that the Lyman-excess sources are more associated with infall signposts than the other objects. Our findings indicate that the most plausible explanation for the Lyman excess is that in addition to the Lyman continuum emission from the early-type star, UV photons are emitted from accretion shocks in the stellar neighbourhood. This result suggests that high-mass stars and/or stellar clusters containing young massive stars may continue to accrete for a long time, even after the development of a compact HII region.

Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, SCOUSE accurately describes the {l, b, v_LSR} distribution of CMZ gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km/s<\sigma<53.1 km/s. A median dispersion of 9.8 km/s, translates to a Mach number, M_3D>28. The gas is distributed throughout several "streams", with projected lengths ~100-250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km/s clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l ,b, v_LSR} space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km/s and 50 km/s clouds and Sgr C.

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.

We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consistently, provided a given initial density distribution. This density distribution of the entire cloud complex is then recursively reconstructed based on available continuum maps, including both single-dish and high-resolution interferometric maps covering a wide frequency range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100 au, i.e. a spatial dynamic range of 10^5. We find that the density distribution of Sagittarius B2 can be reasonably well fitted by applying a superposition of spherical cores with Plummer-like density profiles. In order to reproduce the spectral energy distribution, we position Sgr B2(N) along the line of sight behind the plane containing Sgr B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0 x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700 Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2 for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S), we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that most of the gas content in Sgr B2(M) has already been converted to stars or dispersed.

The physical mechanisms that induce the transformation of a certain mass of gas in new stars are far from being well understood. Infrared bubbles associated with HII regions have been considered to be good samples of investigating triggered star formation. In this paper we report on the investigation of the dust properties of the infrared bubble N4 around the HII region G11.898+0.747, analyzing its interaction with its surroundings and star formation histories therein, with the aim of determining the possibility of star formation triggered by the expansion of the bubble. Using Herschel PACS and SPIRE images with a wide wavelength coverage, we reveal the dust properties over the entire bubble. Meanwhile, we are able to identify six dust clumps surrounding the bubble, with a mean size of 0.50 pc, temperature of about 22 K, mean column density of 1.7 $\times10^{22}$ cm$^{-2}$, mean volume density of about 4.4 $\times10^{4}$ cm$^{-3}$, and a mean mass of 320 $M_{\odot}$. In addition, from PAH emission seen at 8 $\mu$m, free-free emission detected at 20 cm and a probability density function in special regions, we could identify clear signatures of the influence of the HII region on the surroundings. There are hints of star formation, though further investigation is required to demonstrate that N4 is the triggering source.

The VIALACTEA project aims at building a predictive model of star formation in our galaxy. We present the innovative integrated framework and the main technologies and methodologies to reach this ambitious goal.

We report high resolution observations of the $^{12}$CO$(1\rightarrow0)$ and $^{13}$CO$(1\rightarrow0)$ molecular lines in the Carina Nebula and the Gum 31 region obtained with the 22-m Mopra telescope as part of the The Mopra Southern Galactic Plane CO Survey. We cover 8 deg$^2$ from $l = 285^{\circ}$ to 290$^{\circ}$, and from $b = -1.5^{\circ}$ to +0.5$^{\circ}$. The molecular gas column density distributions from both tracers have a similar range of values. By fitting a grey-body function to the observed infrared spectral energy distribution from Herschel maps, we derive gas column densities and dust temperatures. The gas column density has values in the range from $6.3\times\ 10^{20}$ to $1.4\times 10^{23}$ cm$^{-2}$, while the dust temperature has values in the range from 17 to 43 K. The gas column density derived from the dust emission is approximately described by a log-normal function for a limited range of column densities. A high-column density tail is clearly evident for the gas column density distribution, which appears to be a common feature in regions with active star formation. There are regional variations in the fraction of the mass recovered by the CO emission lines with respect to the total mass traced by the dust emission. These variations may be related to changes in the radiation field strength, variation of the atomic to molecular gas fraction across the observed region, differences in the CO molecule abundance with respect to H$_{2}$, and evolutionary stage differences of the molecular clouds that compose the Carina Nebula-Gum 31 complex.

We use the Herschel Hi-GAL survey data to study the spatial distribution in Galactic longitude and latitude of the interstellar medium and of dense, star-forming clumps in the inner Galaxy. The peak position and width of the latitude distribution of the dust column density as well as of number density of compact sources from the band-merged Hi-GAL photometric catalogues are analysed as a function of longitude. The width of the diffuse dust column density traced by the Hi-GAL 500 micron emission varies across the inner Galaxy, with a mean value of 1{\deg}.2-1{\deg}.3, similar to that of the 250um Hi-GAL sources. 70um Hi-GAL sources define a much thinner disk, with a mean FWHM of 0{\deg}.75, and an average latitude of b=0{\deg}.06, coincident with the results from ATLASGAL. The GLAT distribution as a function of GLON shows modulations, both for the diffuse emission and for the compact sources, with ~0{\deg}.2 displacements mostly toward negative latitudes at l~ +40{\deg}, +12{\deg}, -25{\deg} and -40{\deg}. No such modulations can be found in the MIPSGAL 24 or WISE 22 um data when the entire source samples are considered. The distortions revealed by Herschel are interpreted as large-scale bending modes of the Plane. The lack of similar distortions in tracers of more evolved YSOs or stars rules out gravitational instabilities or satellite-induced perturbations, as they should act on both the diffuse and stellar disk components. We propose that the observed bends are caused by incoming flows of extra-planar gas interacting with the gaseous disk. Stars decouple from the gaseous ISM and relax into the stellar disk potential. The time required for the disappearance of the distortions from the diffuse ISM to the relatively evolved YSO stages are compatible with star-formation timescales.

The Evolutionary Map of the Universe (EMU) is a proposed radio continuum survey of the Southern Hemisphere up to declination +30 deg., with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU will use an automated source identification and measurement approach that is demonstrably optimal, to maximise the reliability, utility and robustness of the resulting radio source catalogues. As part of the process of achieving this aim, a "Data Challenge" has been conducted, providing international teams the opportunity to test a variety of source finders on a set of simulated images. The aim is to quantify the accuracy of existing automated source finding and measurement approaches, and to identify potential limitations. The Challenge attracted nine independent teams, who tested eleven different source finding tools. In addition, the Challenge initiators also tested the current ASKAPsoft source-finding tool to establish how it could benefit from incorporating successful features of the other tools. Here we present the results of the Data Challenge, identifying the successes and limitations for this broad variety of the current generation of radio source finding tools. As expected, most finders demonstrate completeness levels close to 100% at 10sigma dropping to levels around 10% by 5sigma. The reliability is typically close to 100% at 10sigma, with performance to lower sensitivities varying greatly between finders. All finders demonstrate the usual trade-off between completeness and reliability, whereby maintaining a high completeness at low signal-to-noise comes at the expense of reduced reliability, and vice-versa. We conclude with a series of recommendations for improving the performance of the ASKAPsoft source-finding tool.

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 $\sim$1.5$\times$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.

The ubiquity of filamentary structure at various scales through out the Galaxy has triggered a renewed interest in their formation, evolution, and role in star formation. The largest filaments can reach up to Galactic scale as part of the spiral arm structure. However, such large scale filaments are hard to identify systematically due to limitations in identifying methodology (i.e., as extinction features). We present a new approach to directly search for the largest, coldest, and densest filaments in the Galaxy, making use of sensitive Herschel Hi-GAL data complemented by spectral line cubes. We present a sample of the 9 most prominent Herschel filaments, including 6 identified from a pilot search field plus 3 from outside the field. These filaments measure 37-99 pc long and 0.6-3.0 pc wide with masses (0.5-8.3)$\times10^4 \, M_\odot$, and beam-averaged ($28"$, or 0.4-0.7 pc) peak H$_2$ column densities of (1.7-9.3)$\times 10^{22} \, \rm{cm^{-2}}$. The bulk of the filaments are relatively cold (17-21 K), while some local clumps have a dust temperature up to 25-47 K. All the filaments are located within <~60 pc from the Galactic mid-plane. Comparing the filaments to a recent spiral arm model incorporating the latest parallax measurements, we find that 7/9 of them reside within arms, but most are close to arm edges. These filaments are comparable in length to the Galactic scale height and therefore are not simply part of a grander turbulent cascade.

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

We report the discovery of a new young stellar cluster in the outer Galaxy located at the position of an IRAS PSC source that has been previously mis-identified as an external galaxy. The cluster is seen in our near-infrared imaging towards IRAS 04186+5143 and in archive Spitzer images confirming the young stellar nature of the sources detected. There is also evidence of sub-clustering seen in the spatial distributions of young stars and of gas and dust. Near- and mid-infrared photometry indicates that the stars exhibit colours compatible with reddening by interstellar and circumstellar dust and are likely to be low- and intermediate-mass YSOs with a large proportion of Class I YSOs. Ammonia and CO lines were detected, with the CO emission well centred near the position of the richest part of the cluster. The velocity of the CO and NH$_3$ lines indicates that the gas is Galactic and located at a distance of about 5.5 kpc, in the outer Galaxy. Herschel data of this region characterise the dust environment of this molecular cloud core where the young cluster is embedded. We derive masses, luminosities and temperatures of the molecular clumps where the young stars reside and discuss their evolutionary stages.

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.

Over the next decade, observations conducted with ALMA and the SKA will reveal the process of mass assembly and accretion onto young stars and will be revolutionary for studies of star formation. Here we summarise the capabilities of ALMA and discuss recent results from its early science observations. We then review infrared and radio variability observations of both young low-mass and high-mass stars. A time domain SKA radio continuum survey of star forming regions is then outlined. This survey will produce radio light-curves for hundreds of young sources, providing for the first time a systematic survey of radio variability across the full range of stellar masses. These light-curves will probe the magnetospheric interactions of young binary systems, the origins of outflows, trace episodic accretion on the central sources and potentially constrain the rotation rates of embedded sources.

We used Herschel Hi-GAL survey data to determine whether massive young stellar objects (MYSOs) are resolved at 70$\mu$m and to study their envelope density distribution. Our analysis of three relatively isolated sources in the l=30{\deg} and l=59{\deg} Galactic fields show that the objects are partially resolved at 70$\mu$m. The Herschel Hi-GAL survey data have a high scan velocity which makes unresolved and partially resolved sources appear elongated in the 70$\mu$m images. We analysed the two scan directions separately and examine the intensity profile perpendicular to the scan direction. Spherically symmetric radiative transfer models with a power law density distribution were used to study the circumstellar matter distribution. Single dish sub-mm data were also included to study how different spatial information affects the fitted density distribution. The density distribution which best fits both the 70$\mu$m intensity profile and SED has an average index of ~0.5. This index is shallower than expected and is probably due to the dust emission from bipolar outflow cavity walls not accounted for in the spherical models. We conclude that 2D axisymmetric models and Herschel images at low scan speeds are needed to better constrain the matter distribution around MYSOs.

The SKA will be a state of the art radiotelescope optimized for both large area surveys as well as for deep pointed observations. In this paper we analyze the impact that the SKA will have on Galactic studies, starting from the immense legacy value of the all-sky survey proposed by the continuum SWG but also presenting some areas of Galactic Science that particularly benefit from SKA observations both surveys and pointed. The planned all-sky survey will be characterized by unique spatial resolution, sensitivity and survey speed, providing us with a wide-field atlas of the Galactic continuum emission. Synergies with existing, current and planned radio Galactic Plane surveys will be discussed. SKA will give the opportunity to create a sensitive catalog of discrete Galactic radio sources, most of them representing the interaction of stars at various stages of their evolution with the environment: complete census of all stage of HII regions evolution; complete census of late stages of stellar evolution such as PNe and SNRs; detection of stellar winds, thermal jets, Symbiotic systems, Chemically Peculiar and dMe stars, active binary systems in both flaring and quiescent states. Coherent emission events like Cyclotron Maser in the magnetospheres of different classes of stars can be detected. Pointed, deep observations will allow new insights into the physics of the coronae and plasma processes in active stellar systems and single stars, enabling the detection of flaring activity in larger stellar population for a better comprehension of the mechanism of energy release in the atmospheres of stars with different masses and age.

We investigate the young stellar population in the Vela Molecular Ridge, Cloud-D (VMR-D), a star forming (SF) region observed by both Spitzer/NASA and Herschel/ESA space telescope. The point source, band-merged, Spitzer-IRAC catalog complemented with MIPS photometry previously obtained is used to search for candidate young stellar objects (YSO), also including sources detected in less than four IRAC bands. Bona fide YSO are selected by using appropriate color-color and color-magnitude criteria aimed to exclude both Galatic and extragalactic contaminants. The derived star formation rate and efficiency are compared with the same quantities characterizing other SF clouds. Additional photometric data, spanning from the near-IR to the submillimeter, are used to evaluate both bolometric luminosity and temperature for 33 YSOs located in a region of the cloud observed by both Spitzer and Herschel. The luminosity-temperature diagram suggests that some of these sources are representative of Class 0 objects with bolometric temperatures below 70 K and luminosities of the order of the solar luminosity. Far IR observations from the Herschel/Hi-GAL key project for a survey of the Galactic plane are also used to obtain a band-merged photometric catalog of Herschel sources aimed to independently search for protostars. We find 122 Herschel cores located on the molecular cloud, 30 of which are protostellar and 92 starless. The global protostellar luminosity function is obtained by merging the Spitzer and Herschel protostars. Considering that 10 protostars are found in both Spitzer and Herschel list it follows that in the investigated region we find 53 protostars and that the Spitzer selected protostars account for approximately two-thirds of the total.

Observations of molecular clouds reveal a complex structure, with gas and dust often arranged in filamentary rather than spherical geometries. The associations of pre- and proto- stellar cores with the filaments suggest a direct link with the process of star formation. Any study of the properties of such filaments requires a representative samples from different enviroments and so an unbiased detection method. We developed such an approach using the Hessian matrix of a surface-brightness distribution to identify filaments and determine their physical and morphological properties. After testing the method on simulated, but realistic filaments, we apply the algorithms to column-density maps computed from Herschel observations of the Galactic Plane obtained by the Hi-GAL project. We identified ~500 filaments, in the longitude range of l=216.5 to l=225.5, with lengths from ~1 pc up to ~30 pc and widths between 0.1 pc and 2.5 pc. Average column densities are between 10^20 cm^-2 and 10^22 cm^-2. Filaments include the majority of dense material with N_H_2 > 6x10^21cm^-2. We find that the pre- and proto-stellar compact sources already identified in the same region are mostly associated with filaments. However, surface densities in excess of the expected critical values for high-mass star formation are only found on the filaments, indicating that these structures are necessary to channel material into the clumps. Furthermore, we analyze the gravitational stability of filaments and discuss their relationship with star formation.

The present study aims at characterizing the massive star forming region G35.20N, which is found associated with at least one massive outflow and contains multiple dense cores, one of them recently found associated with a Keplerian rotating disk. We used ALMA to observe the G35.20N region in the continuum and line emission at 350 GHz. The observed frequency range covers tracers of dense gas (e.g. H13CO+, C17O), molecular outflows (e.g. SiO), and hot cores (e.g. CH3CN, CH3OH). The ALMA 870 um continuum emission map reveals an elongated dust structure (0.15 pc long and 0.013 pc wide) perpendicular to the large-scale molecular outflow detected in the region, and fragmented into a number of cores with masses 1-10 Msun and sizes 1600 AU. The cores appear regularly spaced with a separation of 0.023 pc. The emission of dense gas tracers such as H13CO+ or C17O is extended and coincident with the dust elongated structure. The three strongest dust cores show emission of complex organic molecules characteristic of hot cores, with temperatures around 200 K, and relative abundances 0.2-2x10^(-8) for CH3CN and 0.6-5x10^(-6) for CH3OH. The two cores with highest mass (cores A and B) show coherent velocity fields, with gradients almost aligned with the dust elongated structure. Those velocity gradients are consistent with Keplerian disks rotating about central masses of 4-18 Msun. Perpendicular to the velocity gradients we have identified a large-scale precessing jet/outflow associated with core B, and hints of an east-west jet/outflow associated with core A. The elongated dust structure in G35.20N is fragmented into a number of dense cores that may form massive stars. Based on the velocity field of the dense gas, the orientation of the magnetic field, and the regularly spaced fragmentation, we interpret this elongated structure as the densest part of a 1D filament fragmenting and forming massive stars.

High-energy phenomena in the cosmos, and in particular processes leading to the emission of gamma- rays in the energy range 10 MeV - 100 GeV, play a very special role in the understanding of our Universe. This energy range is indeed associated with non-thermal phenomena and challenging particle acceleration processes. The technology involved in detecting gamma-rays is challenging and drives our ability to develop improved instruments for a large variety of applications. GAMMA-LIGHT is a Small Mission which aims at an unprecedented advance of our knowledge in many sectors of astrophysical and Earth studies research. The Mission will open a new observational window in the low-energy gamma-ray range 10-50 MeV, and is configured to make substantial advances compared with the previous and current gamma-ray experiments (AGILE and Fermi). The improvement is based on an exquisite angular resolution achieved by GAMMA-LIGHT using state-of-the-art Silicon technology with innovative data acquisition. GAMMA-LIGHT will address all astrophysics issues left open by the current generation of instruments. In particular, the breakthrough angular resolution in the energy range 100 MeV - 1 GeV is crucial to resolve patchy and complex features of diffuse sources in the Galaxy as well as increasing the point source sensitivity. This proposal addresses scientific topics of great interest to the community, with particular emphasis on multifrequency correlation studies involving radio, optical, IR, X-ray, soft gamma-ray and TeV emission. At the end of this decade several new observatories will be operational including LOFAR, SKA, ALMA, HAWK, CTA. GAMMA-LIGHT will "fill the vacuum" in the 10 MeV-10 GeV band, and will provide invaluable data for the understanding of cosmic and terrestrial high-energy sources.

The cycling of material from the interstellar medium (ISM) into stars and the return of stellar ejecta into the ISM is the engine that drives the "galactic ecology" in normal spirals, a cornerstone in the formation and evolution of galaxies through cosmic time. Major observational and theoretical challenges need to be addressed in determining the processes responsible for converting the low-density ISM into dense molecular clouds, forming dense filaments and clumps, fragmenting them into stars, OB associations and bound clusters, and characterizing the feedback that limits the rate and efficiency of star formation. This formidable task can be now effectively attacked thanks to the combination of new global-scale surveys of the Milky Way Galactic Plane from infrared to radio wavelengths, offering the possibility of bridging the gap between local and extragalactic star formation studies. The Herschel, Spitzer and WISE mid to far infrared continuum surveys, complemented by analogue surveys from ground-based facilities in the millimetre and radio wavelengths, enables us to measure the Galactic distribution and physical properties of dust on all scales and in all components of the ISM from diffuse clouds to filamentary complexes and tens of thousands of dense clumps. A complementary suite of spectroscopic surveys in various atomic and molecular tracers is providing the chemical fingerprinting of dense clumps and filaments, as well as essential kinematic information to derive distances and thus transform panoramic data into a 3D representation. The latest results emerging from these Galaxy-scale surveys are reviewed. New insights into cloud formation and evolution, filaments and their relationship to channeling gas onto gravitationally-bound clumps, the properties of these clumps, density thresholds for gravitational collapse, and star and cluster formation rates are discussed.

We present a study of the structure of the Galactic interstellar medium through the Delta-variance technique, related to the power spectrum and the fractal properties of infrared/sub-mm maps. Through this method, it is possible to provide quantitative parameters which are useful to characterize different morphological and physical conditions, and to better constrain the theoretical models. In this respect, the Herschel Infrared Galactic Plane Survey carried out at five photometric bands from 70 to 500 \mu m constitutes an unique database for applying statistical tools to a variety of regions across the Milky Way. In this paper, we derive a robust estimate of the power-law portion of the power spectrum of four contiguous 2{\deg}x2{\deg} Hi-GAL tiles located in the third Galactic quadrant (217{\deg} < l < 225{\deg}, -2{\deg} < b < 0{\deg}). The low level of confusion along the line of sight testified by CO observations makes this region an ideal case. We find very different values of the power spectrum slope from tile to tile but also from wavelength to wavelength (2 < \beta < 3), with similarities between fields attributable to components located at the same distance. Thanks to the comparison with models of turbulence, an explanation of the determined slopes in terms of the fractal geometry is also provided, and possible relations with the underlying physics are investigated. In particular, an anti-correlation between ISM fractal dimension and star formation efficiency is found for the two main distance components observed in these fields. A possible link between the fractal properties of the diffuse emission and the resulting clump mass function is discussed.

We investigate the star formation activity in the molecular complex associated with the Galactic HII region Sh2-90, using radio-continuum maps obtained at 1280 MHz and 610 MHz, Herschel Hi-GAL observations at 70 -- 500 microns, and deep near-infrared observation at JHK bands, along with Spitzer observations. Sh2-90 presents a bubble morphology in the mid-IR (size ~ 0.9 pc x 1.6 pc). Radio observations suggest it is an evolved HII region with an electron density ~ 144 cm^-3, emission measure ~ 6.7 x 10^4 cm^-6 pc and a ionized mass ~ 55 Msun. From Hi-GAL observations it is found that the HII region is part of an elongated extended molecular cloud (size ~ 5.6 pc x 9.7 pc, H_2 column density >= 3 x 10^21 cm^-2 and dust temperature 18 -- 27 K) of total mass >= 1 x 10^4 Msun. We identify the ionizing cluster of Sh2-90, the main exciting star being an O8--O9 V star. Five cold dust clumps (mass ~ 8 -- 95 Msun), four mid-IR blobs around B stars, and a compact HII region are found at the edge of the bubble.The velocity information derived from CO (J=3-2) data cubes suggests that most of them are associated with the Sh2-90 region. 129 YSOs are identified (Class I, Class II, and near-IR excess sources). The majority of the YSOs are low mass (<= 3 Msun) sources and they are distributed mostly in the regions of high column density. Four candidate Class 0/I MYSOs have been found; they will possibly evolve to stars of mass >= 15 Msun. We suggest multi-generation star formation is present in the complex. From the evidences of interaction, the time scales involved and the evolutionary status of stellar/protostellar sources, we argue that the star formation at the immediate border/edges of Sh2-90 might have been triggered by the expanding HII region. However, several young sources in this complex are probably formed by some other processes.

Herschel FIR continuum data obtained as part of the Hi-GAL survey have been used, together with the GLIMPSE 8 $\mu$m and MIPSGAL 24 $\mu$m data, to attempt the first 3D-decomposition of dust emission associated with atomic, molecular and ionized gas at 15 arcmin angular resolution. Our initial test case is a 2$\times$2 square degrees region centred on (l,b)=(30$^{\circ}$,0$^{\circ}$), a direction that encompasses the origin point of the Scutum-Crux Arm at the tip of the Galactic Bar. Coupling the IR maps with velocity maps specific for different gas phases (HI 21cm, $^{12}$CO and $^{13}$CO, and RRLs), we estimate the properties of dust blended with each of the gas components and at different Galactocentric distances along the LOS. A statistical Pearson's coefficients analysis is used to study the correlation between the column densities and the intensity of the IR emission. This analysis provides evidence that the 2$\times$2 square degree field under consideration is characterized by the presence of a gas component not accounted for by the standard tracers, possibly associated with warm H$_{2}$ and cold HI. We demonstrate that the IR radiation in the range 8 $\mu$m $<$ $\lambda$ $<$ 500 $\mu$m is systematically dominated by emission originating within the Scutum-Crux Arm. By applying an inversion method, we recover the dust emissivities associated with atomic, molecular and ionized gas. Using the DustEM model we obtain an indication for PAHs depletion in the diffuse ionized gas. However, the main goal of this work is to discuss the impact of the missing column density associated with the dark gas component on the accurate evaluation of the dust properties, and to shed light on the limitations of the inversion method approach when this is applied to a small section of the Galactic Plane and when the working resolution allows sufficient de-blending of the gas components along the LOS.

We present the results of Herschel HOBYS photometric mapping combined with BIMA observations and additional archival data, and perform an in-depth study of the evolutionary phases of the star-forming clumps in W 48A and their surroundings. Age estimates for the compact sources were derived from bolometric luminosities and envelope masses, which were obtained from the dust continuum emission, and agree within an order of magnitude with age estimates from molecular line and radio data. The clumps in W 48A are linearly aligned by age (east-old to west-young): we find a ultra compact (UC) HII region, a young stellar object (YSO) with class II methanol maser emission, a YSO with a massive outflow, and finally the NH_2D prestellar cores from Pillai et al. This remarkable positioning reflects the (star) formation history of the region. We find that it is unlikely that the star formation in the W 48A molecular cloud was triggered by the UCHII region and discuss the Aquila supershell expansion as a mayor influence on the evolution of W 48A. We conclude that the combination of Herschel continuum data with interferometric molecular line and radio continuum data is important to derive trustworthy age estimates and interpret the origin of large scale structures through kinematic information.

(Abridged) METHODS: We performed SMA observations at 1.3 mm with both the most extended and compact array configurations, providing sub-arcsecond and high sensitivity maps of various molecular lines, including both hot-core and outflow tracers. We also reconstruct the spectral energy distribution of the region from millimeter to near infrared wavelengths, using the Herschel/Hi-GAL maps, as well as archival data. RESULTS: From the spectral energy distribution, we derive a bolometric luminosity of about 4x10^4 Lsun. Our interferometric observations reveal that the distribution of dense gas and dust in the HMC is significantly flattened and extends up to a radius of 8000 AU from the center of radio continuum and maser emission in the region. The equatorial plane of this HMC is strictly perpendicular to the elongation of the collimated bipolar outflow, as imaged on scales of about 0.1-0.5 pc in the main CO isotopomers as well as in the SiO(5-4) line. In the innermost HMC regions (ca. 1000 AU), the velocity field traced by the CH3CN(12_K-11_K) line emission shows that molecular gas is both expanding along the outflow direction following a Hubble-law, and rotating about the outflow axis, in agreement with the (3-D) velocity field traced by methanol masers. The velocity field associated with rotation indicates a dynamical mass of 19 Msun at the center of the core. The latter is likely to be concentrated in a single O9.5 ZAMS star, consistent with the estimated bolometric luminosity of G023.01-00.41. The physical properties of the CO(2-1) outflow emission, such as its momentum rate 6x10^-3 Msun km/s /yr and its outflow rate 2x10^-4 Msun/yr, support our estimates of the luminosity (and mass) of the embedded young stellar object.

Context. Stars form in dense, dusty clumps of molecular clouds, but little is known about their origin, their evolution and their detailed physical properties. In particular, the relationship between the mass distribution of these clumps (also known as the "clump mass function", or CMF) and the stellar initial mass function (IMF), is still poorly understood. Aims. In order to better understand how the CMF evolve toward the IMF, and to discern the "true" shape of the CMF, large samples of bona-fide pre- and proto-stellar clumps are required. Two such datasets obtained from the Herschel infrared GALactic Plane Survey (Hi-GAL) have been described in paper I. Robust statistical methods are needed in order to infer the parameters describing the models used to fit the CMF, and to compare the competing models themselves. Methods. In this paper we apply Bayesian inference to the analysis of the CMF of the two regions discussed in Paper I. First, we determine the Bayesian posterior probability distribution for each of the fitted parameters. Then, we carry out a quantitative comparison of the models used to fit the CMF. Results. We have compared the results from several methods implementing Bayesian inference, and we have also analyzed the impact of the choice of priors and the influence of various constraints on the statistical conclusions for the preferred values of the parameters. We find that both parameter estimation and model comparison depend on the choice of parameter priors. Conclusions. Our results confirm our earlier conclusion that the CMFs of the two Hi-GAL regions studied here have very similar shapes but different mass scales. Furthermore, the lognormal model appears to better describe the CMF measured in the two Hi-GAL regions studied here. However, this preliminary conclusion is dependent on the choice of parameters priors.

Maser lines from different molecular species, including water, hydroxyl, and methanol, are common observational phenomena associated with massive star forming regions. In particular, the methanol maser appears as an ideal tool to study the early phases of massive star formation. However, it is difficult to establish the exact start of the methanol maser phase, and it would then be interesting to detect and study low-flux density methanol masers (i.e., < 0.1 Jy or even << 0.1 Jy), in order to determine if they can effectively be used to mark a specific evolutionary phase in high-mass star formation. Past surveys have been unable to systematically detect many low-flux density methanol masers, and thus we do not yet know how many such masers exist in the Galaxy and what is their physical nature. Out of a sample of 107 observed Hi-GAL sources we detected a total of 32 methanol masers, with 22 sources being new and weak (median peak flux density 0.07 Jy) detections, in the Galactic longitude range [32.0, 59.8]deg. We also detected 12 6.035-GHz OH maser, with 9 objects being new detections. Our survey covers a similar range of source distances as the "Arecibo Methanol Maser Galactic Plane Survey", but the methanol masers detected by us are clearly shifted towards lower integrated flux densities. The newly detected methanol masers are mostly of low-luminosity and, except for some sources, their weakness is not due to distance effects or positional offsets. No specific correlation is found with the physical parameters of the Hi-GAL clumps, except for sources with both CH3OH and OH masers which tend to have higher mass and luminosity. The intensity of the methanol masers correlates well with the velocity range of the maser emission, which suggests that the low brightness of these masers is related to the number of maser spots in the emitting region and their evolution with time.

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.

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

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.

Sagittarius B2 (Sgr B2) is one of the most massive and luminous star-forming regions in the Galaxy and shows chemical and physical conditions similar to those in distant extragalactic starbursts. We present large-scale far-IR/submm photometric images and spectroscopic maps taken with the PACS and SPIRE instruments onboard Herschel. The spectra towards the Sgr B2 star-forming cores, B2(M) and B2(N), are characterized by strong CO line emission, emission lines from high-density tracers (HCN, HCO+, and H2S), [N II] 205 um emission from ionized gas, and absorption lines from hydride molecules (OH+, H2O+, H2O, CH+, CH, SH+, HF, NH, NH2, and NH3). The rotational population diagrams of CO suggest the presence of two gas temperature components: an extended warm component, which is associated with the extended envelope, and a hotter component, which is seen towards the B2(M) and B2(N) cores. As observed in other Galactic Center clouds, the gas temperatures are significantly higher than the dust temperatures inferred from photometric images. We determined far-IR and total dust masses in the cores. Non-local thermodynamic equilibrium models of the CO excitation were used to constrain the averaged gas density in the cores. A uniform luminosity ratio is measured along the extended envelope, suggesting that the same mechanism dominates the heating of the molecular gas at large scales. The detection of high-density molecular tracers and of strong [N II] 205 um line emission towards the cores suggests that their morphology must be clumpy to allow UV radiation to escape from the inner HII regions. Together with shocks, the strong UV radiation field is likely responsible for the heating of the hot CO component. At larger scales, photodissociation regions models can explain both the observed CO line ratios and the uniform L(CO)/LFIR luminosity ratios.

We present the first overview of the Herschel observations of the nearby high-mass star-forming region NGC 7538, taken as part of the Herschel imaging study of OB Young Stellar objects (HOBYS) Key Programme. These PACS and SPIRE maps cover an approximate area of one square degree at five submillimeter and far-infrared wavebands. We have identified 780 dense sources and classified 224 of those. With the intention of investigating the existence of cold massive starless or class 0-like clumps that would have the potential to form intermediate- to high-mass stars, we further isolate 13 clumps as the most likely candidates for followup studies. These 13 clumps have masses in excess of 40 M_sun and temperatures below 15 K. They range in size from 0.4 pc to 2.5 pc and have densities between 3x10^3 cm^-3 to 4x10^4 cm^-3. Spectral energy distributions are then used to characterize their energetics and evolutionary state through a luminosity-mass diagram. NGC 7538 has a highly filamentary structure, previously unseen in the dust continuum of existing submillimeter surveys. We report the most complete imaging to date of a large, evacuated ring of material in NGC 7538 which is bordered by many cool sources.

Theoretical models suggest that massive stars form via disk-mediated accretion, with bipolar outflows playing a fundamental role. A recent study toward massive molecular outflows has revealed a decrease of the SiO line intensity as the object evolves. The present study aims at characterizing the variation of the molecular outflow properties with time, and at studying the SiO excitation conditions in outflows associated with massive YSOs. We used the IRAM30m telescope to map 14 massive star-forming regions in the SiO(2-1), SiO(5-4) and HCO+(1-0) outflow lines, and in several dense gas and hot core tracers. Hi-GAL data was used to improve the spectral energy distributions and the L/M ratio, which is believed to be a good indicator of the evolutionary stage of the YSO. We detect SiO and HCO+ outflow emission in all the sources, and bipolar structures in six of them. The outflow parameters are similar to those found toward other massive YSOs. We find an increase of the HCO+ outflow energetics as the object evolve, and a decrease of the SiO abundance with time, from 10^(-8) to 10^(-9). The SiO(5-4) to (2-1) line ratio is found to be low at the ambient gas velocity, and increases as we move to high velocities, indicating that the excitation conditions of the SiO change with the velocity of the gas (with larger densities and/or temperatures for the high-velocity gas component). The properties of the SiO and HCO+ outflow emission suggest a scenario in which SiO is largely enhanced in the first evolutionary stages, probably due to strong shocks produced by the protostellar jet. As the object evolves, the power of the jet would decrease and so does the SiO abundance. During this process, however, the material surrounding the protostar would have been been swept up by the jet, and the outflow activity, traced by entrained molecular material (HCO+), would increase with time.

We present the first Herschel PACS and SPIRE photometric observations in a portion of the outer Galaxy ($216.5^{\circ} \lesssim \ell \lesssim 225.5^{\circ}$ and $-2^{\circ} \lesssim b \lesssim 0^{\circ}$) as a part of the Hi-GAL survey. The maps between 70 and 500 $\mu$m, the derived column density and temperature maps, and the compact source catalog are presented. NANTEN CO(1-0) line observations are used to derive cloud kinematics and distances, so that we can estimate distance-dependent physical parameters of the compact sources (cores and clumps) having a reliable spectral energy distribution, that we separate in 255 proto-stellar and 688 starless. Both typologies are found in association with all the distance components observed in the field, up to $\sim 5.8$ kpc, testifying the presence of star formation beyond the Perseus arm at these longitudes. Selecting the starless gravitationally bound sources we identify 590 pre-stellar candidates. Several sources of both proto- and pre-stellar nature are found to exceed the minimum requirement for being compatible with massive star formation, based on the mass-radius relation. For the pre-stellar sources belonging to the Local arm ($d\lesssim1.5$ kpc) we study the mass function, whose high-mass end shows a power-law $N(\log M) \propto M^{-1.0 \pm 0.2}$. Finally, we use a luminosity vs mass diagram to infer the evolutionary status of the sources, finding that most of the proto-stellar are in the early accretion phase (with some cases compatible with a Class I stage), while for pre-stellar sources, in general, accretion has not started yet.

Super star clusters are the end product of star formation under the most extreme conditions. As such, studying how their final stellar populations are assembled from their natal progenitor gas clouds can provide strong constraints on star formation theories. An obvious place to look for the initial conditions of such extreme stellar clusters are gas clouds of comparable mass and density, with no star formation activity. We present a method to identify such progenitor gas clouds and demonstrate the technique for the gas in the inner few hundred pc of our Galaxy. The method highlights three clouds in the region with similar global physical properties to the previously identified extreme cloud, G0.253+0.016, as potential young massive cluster (YMC) precursors. The fact that four potential YMC progenitor clouds have been identified in the inner 100 pc of the Galaxy, but no clouds with similar properties have been found in the whole first quadrant despite extensive observational efforts, has implications for cluster formation/destruction rates across the Galaxy. We put forward a scenario to explain how such dense gas clouds can arise in the Galactic centre environment, in which YMC formation is triggered by gas streams passing close to the minimum of the global Galactic gravitational potential at the location of the central supermassive black hole, Sgr A*. If this triggering mechanism can be verified, we can use the known time interval since closest approach to Sgr A* to study the physics of stellar mass assembly in an extreme environment as a function of absolute time.

Context. G29.96-0.02 is a high-mass star-forming cloud observed at 70, 160, 250, 350, and 500 microns as part of the Herschel survey of the Galactic Plane during the Science Demonstration Phase. Aims. We wish to conduct a far-infrared study of the sources associated with this star-forming region by estimating their physical properties and evolutionary stage, and investigating the clump mass function, the star formation efficiency and rate in the cloud. Methods. We have identified the Hi-GAL sources associated with the cloud, searched for possible counterparts at centimeter and infrared wavelengths, fitted their spectral energy distribution and estimated their physical parameters. Results. A total of 198 sources have been detected in all 5 Hi-GAL bands, 117 of which are associated with 24 microns emission and 87 of which are not associated with 24 microns emission. We called the former sources 24 microns-bright and the latter ones 24 microns-dark. The [70-160] color of the 24 microns-dark sources is smaller than that of the 24 microns-bright ones. The 24 microns-dark sources have lower L_bol and L_bol/M_env than the 24 microns-bright ones for similar M_env, which suggests that they are in an earlier evolutionary phase. The G29-SFR cloud is associated with 10 NVSS sources and with extended centimeter continuum emission well correlated with the 70 microns emission. Most of the NVSS sources appear to be early B or late O-type stars. The most massive and luminous Hi-GAL sources in the cloud are located close to the G29-UC region, which suggests that there is a privileged area for massive star formation towards the center of the G29-SFR cloud. Almost all the Hi-GAL sources have masses well above the Jeans mass but only 5% have masses above the virial mass, which indicates that most of the sources are stable against gravitational collapse. The sources with M_env > M_virial and that ...

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.

The Herschel survey of the Galactic Plane (Hi-GAL) provides a unique opportunity to study star formation over large areas of the sky and different environments in the Milky Way. We use the best studied Hi-GAL fields to date, two 2x2 tiles centered on (l, b) = (30, 0) deg and (l, b) = (59, 0) deg, to study the star formation activity using a large sample of well selected young stellar objects (YSOs). We estimate the star formation rate (SFR) for these fields using the number of candidate YSOs and their average time scale to reach the Zero Age Main Sequence, and compare it with the rate estimated using their integrated luminosity at 70 micron combined with an extragalactic star formation indicator. We measure a SFR of (9.5 +- 4.3)*10^{-4} Msol/yr and (1.6 +- 0.7)*10^{-4} Msol/yr with the source counting method, in l=30 deg and l=59 deg, respectively. Results with the 70 micron estimator are (2.4 +- 0.4)*10^{-4} Msol/yr and (2.6 +- 1.1)*10^{-6} Msol/yr. Since the 70 micron indicator is derived from averaging extragalactic star forming complexes, we perform an extrapolation of these values to the whole Milky Way and obtain SFR_{MW} = (0.71 +- 0.13) Msol/yr from l = 30 deg and SFR_{MW} = (0.10 +- 0.04) Msol/yr from l=59 deg. The estimates in l=30 deg are in agreement with the most recent results on the Galactic star formation activity, indicating that the characteristics of this field are likely close to those of the star-formation dominated galaxies used for its derivation. Since the sky coverage is limited, this analysis will improve when the full Hi-GAL survey will be available.

We present ATCA observations of the H2O maser line and radio continuum at 18.0GHz and 22.8GHz, toward a sample of 192 massive star forming regions containing several clumps already imaged at 1.2mm. The main aim of this study is to investigate the water maser and centimeter continuum emission (likely tracing thermal free-free emission) in sources at different evolutionary stages, using the evolutionary classifications proposed by Palla et al (1991) and Molinari et al (2008). We used the recently comissioned CABB backend at ATCA obtaining images with 20arcsec resolution in the 1.3cm continuum and H2O maser emission, in all targets. For the evolutionary analysis of the sources we used the millimeter continuum emission from Beltran et al (2006) and the infrared emission from the MSX Point Source Catalogue. We detect centimeter continuum emission in 88% of the observed fields with a typical rms noise level of 0.45mJy/beam. Most of the fields show a single radio continuum source, while in 20% of them we identify multiple components. A total of 214 centimeter continuum sources have been identified, likely tracing optically thin HII regions, with physical parameters typical of both extended and compact HII regions. Water maser emission was detected in 41% of the regions, resulting in a total of 85 distinct components. The low angular (20arcsec) and spectral (14km/s) resolutions do not allow a proper analysis of the water maser emission, but suffice to investigate its association with the continuum sources. We have also studied the detection rate of HII regions in the two types of IRAS sources defined by Palla et (1991) on the basis of the IRAS colours: High and Low. No significant differences are found, with large detection rates (>90%) for both High and Low sources. We classify the millimeter and infrared sources in our fields in three evolutionary stages following the scheme presented by ...

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.

Stars form in dense, dusty clumps of molecular clouds, but little is known about their origin and evolution. In particular, the relationship between the mass distribution of these clumps (also known as the "clump mass function", or CMF) and the stellar initial mass function (IMF), is still poorly understood. In order to discern the "true" shape of the CMF and to better understand how the CMF may evolve toward the IMF, large samples of bona-fide pre- and proto-stellar clumps are required. The sensitive observations of the Herschel Space Observatory (HSO) are now allowing us to look at large clump populations in various clouds with different physical conditions. We analyse two fields in the Galactic plane mapped by HSO during its science demonstration phase, as part of the more complete and unbiased Herschel infrared GALactic Plane Survey (Hi-GAL). These fields undergo a source-extraction and flux-estimation pipeline, which allows us to obtain a sample with thousands of clumps. Starless and proto-stellar clumps are separated using both color and positional criteria to find those coincident with MIPS 24 micron sources. We describe the probability density functions of the power-law and lognormal models that are used to fit the CMFs, and we then find their best-fit parameters. For the lognormal model we apply several statistical techniques to the data and compare their results. The CMFs of the two SDP fields show very similar shapes, but very different mass scales. This similarity is confirmed by the values of the best-fit parameters of either the power-law or lognormal model. The power-law model leads to almost identical CMF slopes, whereas the lognormal model shows that the CMFs have similar widths. The similar CMF shape but different mass scale represents an evidence that the overall process of star formation in the two regions is very different.

The conversion of gas into stars is a fundamental process in astrophysics and cosmology. Stars are known to form from the gravitational collapse of dense clumps in interstellar molecular clouds, and it has been proposed that the resulting star formation rate is proportional to either the amount of mass above a threshold gas surface density, or the gas volume density. These star-formation prescriptions appear to hold in nearby molecular clouds in our Milky Way Galaxy's disk as well as in distant galaxies where the star formation rates are often much larger. The inner 500 pc of our Galaxy, the Central Molecular Zone (CMZ), contains the largest concentration of dense, high-surface density molecular gas in the Milky Way, providing an environment where the validity of star-formation prescriptions can be tested. Here we show that by several measures, the current star formation rate in the CMZ is an order-of-magnitude lower than the rates predicted by the currently accepted prescriptions. In particular, the region 1 deg < l < 3.5 deg, |b| < 0.5 deg contains ~10^7 Msun of dense molecular gas -- enough to form 1000 Orion-like clusters -- but the present-day star formation rate within this gas is only equivalent to that in Orion. In addition to density, another property of molecular clouds, such as the amplitude of turbulent motions, must be included in the star-formation prescription to predict the star formation rate in a given mass of molecular gas.

We present a Herschel far-infrared study towards the rich massive star- forming complex G305, utilising PACS 70, 160 {\mu}m and SPIRE 250, 350, and 500 {\mu}m observations from the Hi-GAL survey of the Galactic plane. The focus of this study is to identify the embedded massive star-forming population within G305, by combining far-infrared data with radio continuum, H2O maser, methanol maser, MIPS, and Red MSX Source survey data available from previous studies. By applying a frequentist technique we are able to identify a sample of the most likely associations within our multi-wavelength dataset, that can then be identified from the derived properties obtained from fitted spectral energy distributions (SEDs). By SED modelling using both a simple modified blackbody and fitting to a comprehensive grid of model SEDs, some 16 candidate associations are identified as embedded massive star-forming regions. We derive a two-selection colour criterion from this sample of log(F70/F500)\geq 1 and log(F160/F350)\geq 1.6 to identify an additional 31 embedded massive star candidates with no associated star-formation tracers. Using this result we can build a picture of the present day star-formation of the complex, and by extrapolating an initial mass function, suggest a current population of \approx 2 \times 10^4 young stellar objects (YSOs) present, corresponding to a star formation rate (SFR) of 0.01-0.02 M\odot yr^-1. Comparing this resolved star formation rate, to extragalactic star formation rate tracers (based on the Kennicutt-Schmidt relation), we find the star formation activity is underestimated by a factor of \geq 2 in comparison to the SFR derived from the YSO population.

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.

As part of the Herschel guaranteed time key program 'HOBYS', we present the photometric survey of the star forming region Vela-C, one of the nearest sites of low-to-high-mass star formation in the Galactic plane. Vela-C has been observed with PACS and SPIRE in parallel mode between 70 um and 500 um over an area of about 3 square degrees. A photometric catalogue has been extracted from the detections in each band, using a threshold of 5 sigma over the local background. Out of this catalogue we have selected a robust sub-sample of 268 sources, of which 75% are cloud clumps and 25% are cores. Their Spectral Energy Distributions (SEDs) have been fitted with a modified black body function. We classify 48 sources as protostellar and 218 as starless. For two further sources, we do not provide a secure classification, but suggest they are Class 0 protostars. From SED fitting we have derived key physical parameters. Protostellar sources are in general warmer and more compact than starless sources. Both these evidences can be ascribed to the presence of an internal source(s) of moderate heating, which also causes a temperature gradient and hence a more peaked intensity distribution. Moreover, the reduced dimensions of protostellar sources may indicate that they will not fragment further. A virial analysis of the starless sources gives an upper limit of 90% for the sources gravitationally bound and therefore prestellar. We fit a power law N(logM) prop M^-1.1 to the linear portion of the mass distribution of prestellar sources. This is in between that typical of CO clumps and those of cores in nearby star-forming regions. We interpret this as a result of the inhomogeneity of our sample, which is composed of comparable fractions of clumps and cores.

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.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

Past and recent observations have revealed unexpected variations in the FIR-mm dust emissivity. In the Herschel spectral range, those are often referred to as a 500{\mu}m emission excess. Several dust emission models have been developed to interpret astrophysical data in the FIR-mm domain. However, these are commonly unable to fully reconcile theoretical predictions with observations. In contrast, the recently revised two level system (TLS) model seems to provide a promising way of interpreting the existing data. The newly available Herschel Hi-GAL data which covers most of the inner Milky-Way offers a unique opportunity to investigate possible variations in the dust emission properties both with wavelength and environment. By combining the IRIS 100 {\mu}m with the Hi-GAL 160, 250, 350 and 500 {\mu}m data, we model the dust emission spectra in each pixel of the Hi-GAL maps, using both the TLS model and, for comparison, a single modified black-body fit. The effect of temperature mixing along the line of sight is investigated. We find a slight decrease in the dust temperature with distance from the Galactic center. We also report the detection of a significant 500 {\mu}m emissivity excess in the peripheral regions of the plane (35\circ<|l|<70\circ) of about 13-15% of the emissivity, that can reach up to 20% in some HII regions. We present the spatial distribution of the best-fit values for the two main parameters of the TLS model, i.e. the charge correlation length, lc, used to characterize the disordered charge distribution (DCD) part of the model, and the amplitude A of the TLS processes, with respect to the DCD effect. They highlight the plausible existence of an overall gradient with distance to the Galactic center. A comparison with previous findings in the solar neighborhood shows that the local value of the excess is less than expected from the Galactic gradient observed here.

Young massive clusters (YMCs) with stellar masses of 10^4 - 10^5 Msun and core stellar densities of 10^4 - 10^5 stars per cubic pc are thought to be the `missing link' between open clusters and extreme extragalactic super star clusters and globular clusters. As such, studying the initial conditions of YMCs offers an opportunity to test cluster formation models across the full cluster mass range. G0.253+0.016 is an excellent candidate YMC progenitor. We make use of existing multi-wavelength data including recently available far-IR continuum (Herschel/Hi-GAL) and mm spectral line (HOPS and MALT90) data and present new, deep, multiple-filter, near-IR (VLT/NACO) observations to study G0.253+0.016. These data show G0.253+0.016 is a high mass (1.3x10^5 Msun), low temperature (T_dust~20K), high volume and column density (n ~ 8x10^4 cm^-3; N_{H_2} ~ 4x10^23 cm^-2) molecular clump which is close to virial equilibrium (M_dust ~ M_virial) so is likely to be gravitationally-bound. It is almost devoid of star formation and, thus, has exactly the properties expected for the initial conditions of a clump that may form an Arches-like massive cluster. We compare the properties of G0.253+0.016 to typical Galactic cluster-forming molecular clumps and find it is extreme, and possibly unique in the Galaxy. This uniqueness makes detailed studies of G0.253+0.016 extremely important for testing massive cluster formation models.

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.

The molecular clouds Lupus 1, 3 and 4 were mapped with the Mopra telescope at 3 and 12 mm. Emission lines from high density molecular tracers were detected, i.e. NH$_3$ (1,1), NH$_3$ (2,2), N$_2$H$^+$ (1-0), HC$_3$N (3-2), HC$_3$N (10-9), CS (2-1), CH$_3$OH (2$_0-1_0$)A$^+$ and CH$_3$OH (2$_{-1}-1_{-1}$)E. Velocity gradients of more than 1 km s$^{-1}$ are present in Lupus 1 and 3 and multiple gas components are present in these clouds along some lines of sight. Lupus 1 is the cloud richest in high density cores, 8 cores were detected in it, 5 cores were detected in Lupus 3 and only 2 in Lupus 4. The intensity of the three species HC$_3$N, NH$_3$ and N$_2$H$^+$ changes significantly in the various cores: cores that are brighter in HC$_3$N are fainter or undetected in NH$_3$ and N$_2$H$^+$ and vice versa. We found that the column density ratios HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ change by one order of magnitude between the cores, indicating that also the chemical abundance of these species is different. The time dependent chemical code that we used to model our cores shows that the HC$_3$N/N$_2$H$^+$ and HC$_3$N/NH$_3$ ratios decrease with time therefore the observed column density of these species can be used as an indicator of the chemical evolution of dense cores. On this base we classified 5 out of 8 cores in Lupus 1 and 1 out of 5 cores in Lupus 3 as very young protostars or prestellar cores. Comparing the millimetre cores population with the population of the more evolved young stellar objects identified in the Spitzer surveys, we conclude that in Lupus 3 the bulk of the star formation activity has already passed and only a moderate number of stars are still forming. On the contrary, in Lupus 1 star formation is on-going and several dense cores are still in the pre--/proto--stellar phase. Lupus 4 is at an intermediate stage, with a smaller number of individual objects.

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 present the data reduction pipeline for the Hi-GAL survey. Hi-GAL is a key project of the Herschel satellite which is mapping the inner part of the Galactic plane (|l| <= 70\cdot and |b| <= 1\cdot), using 2 PACS and 3 SPIRE frequency bands, from 70{\mu}m to 500{\mu}m. Our pipeline relies only partially on the Herschel Interactive Standard Environment (HIPE) and features several newly developed routines to perform data reduction, including accurate data culling, noise estimation and minimum variance map-making, the latter performed with the ROMAGAL algorithm, a deep modification of the ROMA code already tested on cosmological surveys. We discuss in depth the properties of the Hi-GAL Science Demonstration Phase (SDP) data.

Thermal images of cold dust in the Central Molecular Zone of the Milky Way, obtained with the far-infrared cameras on-board the Herschel satellite, reveal a 3x10^7 solar masses ring of dense and cold clouds orbiting the Galactic Center. Using a simple toy-model, an elliptical shape having semi-major axes of 100 and 60 parsecs is deduced. The major axis of this 100-pc ring is inclined by about 40 degrees with respect to the plane-of-the-sky and is oriented perpendicular to the major axes of the Galactic Bar. The 100-pc ring appears to trace the system of stable x_2 orbits predicted for the barred Galactic potential. Sgr A* is displaced with respect to the geometrical center of symmetry of the ring. The ring is twisted and its morphology suggests a flattening-ratio of 2 for the Galactic potential, which is in good agreement with the bulge flattening ratio derived from the 2MASS data.

We use a Minimum Spanning Tree algorithm to characterize the spatial distribution of Galactic Far-IR sources and derive their clustering properties. We aim to reveal the spatial imprint of different types of star forming processes, e.g. isolated spontaneous fragmentation of dense molecular clouds, or events of triggered star formation around HII regions, and highlight global properties of star formation in the Galaxy. We plan to exploit the entire Hi-GAL survey of the inner Galactic plane to gather significant statistics on the clustering properties of star forming regions, and to look for possible correlations with source properties such as mass, temperature or evolutionary stage. In this paper we present a pilot study based on the two 2x2 square degree fields centered at longitudes l=30 and l=59 obtained during the Science Demonstration Phase (SDP) of the Herschel mission. We find that over half of the clustered sources are associated with HII regions and infrared dark clouds. Our analysis also reveals a smooth chromatic evolution of the spatial distribution where sources detected at short-wavelengths, likely proto-stars surrounded by warm circumstellar material emitting in the far-infrared, tend to be clustered in dense and compact groups around HII regions while sources detected at long-wavelengths, presumably cold and dusty density enhancements of the ISM emitting in the sub-millimeter, are distributed in larger and looser groups.

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.

Some theoretical models propose that O-B stars form via accretion, in a similar fashion to low-mass stars. Jet-driven molecular outflows play an important role in this scenario, and their study can help to understand the process of high-mass star formation and the different evolutionary phases involved. Observations towards low-mass protostars so far favour an evolutionary picture in which jets are always associated with Class 0 objects while more evolved Class I/II objects show less evidence of powerful jets. The present study aims at checking whether an analogous picture can be found in the high-mass case. The IRAM 30-m telescope (Spain) has been used to perform single-pointing SiO(2-1) and (3-2) observations towards a sample of 57 high-mass molecular clumps in different evolutionary stages. Continuum data at different wavelengths, from mid-IR to 1.2 mm, have been gathered to build the spectral energy distributions of all the clumps and estimate their bolometric luminosities. SiO emission at high velocities, characteristic of molecular jets, is detected in 88% of our sources, a very high detection rate indicating that there is ongoing star formation activity in most of the sources of our sample. The SiO(2-1) luminosity drops with L/M, which suggests that jet activity declines as time evolves. This represents the first clear evidence of a decrease of SiO outflow luminosity with time in a homogeneous sample of high-mass molecular clumps in different evolutionary stages. The SiO(3-2) to SiO(2-1) integrated intensity ratio shows only minor changes with evolutionary state.

(Abridged) We present a new method for detecting and measuring compact sources in conditions of intense, and highly variable, fore/background. While all most commonly used packages carry out the source detection over the signal image, our proposed method builds from the measured image a "curvature" image by double-differentiation in four different directions. In this way point-like as well as resolved, yet relatively compact, objects are easily revealed while the slower varying fore/background is greatly diminished. Candidate sources are then identified by looking for pixels where the curvature exceeds, in absolute terms, a given threshold; the methodology easily allows us to pinpoint breakpoints in the source brightness profile and then derive reliable guesses for the sources extent. Identified peaks are fit with 2D elliptical Gaussians plus an underlying planar inclined plateau, with mild constraints on size and orientation. Mutually contaminating sources are fit with multiple Gaussians simultaneously using flexible constraints. We ran our method on simulated large-scale fields with 1000 sources of different peak flux overlaid on a realistic realization of diffuse background. We find detection rates in excess of 90% for sources with peak fluxes above the 3-sigma signal noise limit; for about 80% of the sources the recovered peak fluxes are within 30% of their input values.

The first Herschel Hi-Gal images of the galactic plane unveil the far-infrared diffuse emission of the interstellar medium with an unprecedented angular resolution and sensitivity. In this paper, we present the first analysis of these data in combination with that of Spitzer Glimpse & Mipsgal. We selected a relatively diffuse and low excitation region of the l~59\,^{\circ} Hi-Gal Science Demonstration Phase field to perform a pixel by pixel fitting of the 8 to 500 microns SED using the DustEM dust emission model. We derived maps of the Very Small Grains (VSG) and PAH abundances from the model. Our analysis allows us to illustrate that the Aromatic Infrared Bands (AIB) intensity does not trace necessarily the PAH abundance but rather the product of "abundance x column density x intensity of the exciting radiation field". We show that the spatial structure of PACS70microns map resembles the shorter wavelengths (e.g. IRAC8microns) maps, because they trace both the intensity of exciting radiation field and column density. We also show that the modeled VSG contribution to PACS70microns (PACS160microns) band intensity can be up to 50% (7%). The interpretation of diffuse emission spectra at these wavelengths must take stochastically heated particles into account. Finally, this preliminary study emphasizes the potential of analyzing the full dust SED sampled by Herschel and Spitzer data, with a physical dust model (DustEM) to reach the properties of the dust at simultaneously large and small scales.

Variations in the dust emissivity are critical for gas mass determinations derived from far-infrared observations, but also for separating dust foreground emission from the Cosmic Microwave Background (CMB). Hi-GAL observations allow us for the first time to study the dust emissivity variations in the inner regions of the Galactic plane at resolution below 1 degree. We present maps of the emissivity spectral index derived from the combined Herschel PACS 160 \mu m, SPIRE 250 \mu m, 350 \mu m, and 500 \mu m data, and the IRIS 100 \mu m data, and we analyze the spatial variations of the spectral index as a function of dust temperature and wavelength in the two Science Demonstration Phase Hi-GAL fields, centered at l=30{\deg} and l=59{\deg}. Applying two different methods, we determine both dust temperature and emissivity spectral index between 100 and 500 \mu m, at an angular resolution of 4'. Combining both fields, the results show variations of the emissivity spectral index in the range 1.8-2.6 for temperatures between 14 and 23 K. The median values of the spectral index are similar in both fields, i.e. 2.3 in the range 100-500 \mu m, while the median dust temperatures are equal to 19.1 K and 16.0 K in the l=30{\deg} and l=59{\deg} field, respectively. Statistically, we do not see any significant deviations in the spectra from a power law emissivity between 100 and 500 \mu m. We confirm the existence of an inverse correlation between the emissivity spectral index and dust temperature, found in previous analyses.

We present single-dish (sub)millimeter observations of gas and dust in the Galactic high-mass star-forming region G19.61-0.23, with the aim of studying the large-scale properties and physical conditions of the molecular gas across the region. The final aim is to compare the large-scale (about 100 pc) properties with the small-scale (about 3 pc) properties and to consider possible implications for extragalactic studies. We have mapped CO isotopologues in the J=1-0 transition using the FCRAO-14m telescope and the J=2-1 transition using the IRAM-30m telescope. We have also used data from the ATLASGAL survey and from the BU-FCRAO Galactic Ring Survey, as well as the Spitzer infrared Galactic plane surveys GLIMPSE and MIPSGAL to characterize the star-formation activity within the molecular clouds. We reveal a population of molecular clumps in the 13CO(1-0) emission. Our analysis of the 13CO suggests that the virial parameter (ratio of kinetic to gravitational energy) varies over an order of magnitude between clumps that are unbound and some that are apparently "unstable". This conclusion is independent of whether they show evidence of ongoing star formation. We find that the majority of ATLASGAL sources have MIPSGAL counterparts with luminosities in the range 10^4 - 5 10^4 Lsun and are presumably forming relatively massive stars. We compare our results with previous extragalactic studies of the nearby spiral galaxies M31 and M33; and the blue compact dwarf galaxy Henize2-10. We find that the main giant molecular cloud surrounding G19.61-0.23 has physical properties typical for Galactic GMCs and which are comparable to the GMCs in M31 and M33. However, the GMC studied here shows smaller surface densities and masses than the clouds identified in Henize2-10 and associated with super star cluster formation.

We demonstrate the use of the 3D Monte Carlo radiative transfer code PHAETHON to model infrared-dark clouds (IRDCs) that are externally illuminated by the interstellar radiation field (ISRF). These clouds are believed to be the earliest observed phase of high-mass star formation, and may be the high-mass equivalent of lower-mass prestellar cores. We model three different cases as examples of the use of the code, in which we vary the mass, density, radius, morphology and internal velocity field of the IRDC. We show the predicted output of the models at different wavelengths chosen to match the observing wavebands of Herschel and Spitzer. For the wavebands of the long- wavelength SPIRE photometer on Herschel, we also pass the model output through the SPIRE simulator to generate output images that are as close as possible to the ones that would be seen using SPIRE. We then analyse the images as if they were real observations, and compare the results of this analysis with the results of the radiative transfer models. We find that detailed radiative transfer modelling is necessary to accurately determine the physical parameters of IRDCs (e.g. dust temperature, density profile). This method is applied to study G29.55+00.18, an IRDC observed by the Herschel Infrared Galactic Plane survey (Hi-GAL), and in the future it will be used to model a larger sample of IRDCs from the same survey.

The Spectral and Photometric Imaging Receiver (SPIRE), is the Herschel Space Observatory`s submillimetre camera and spectrometer. It contains a three-band imaging photometer operating at 250, 350 and 500 microns, and an imaging Fourier Transform Spectrometer (FTS) which covers simultaneously its whole operating range of 194-671 microns (447-1550 GHz). The SPIRE detectors are arrays of feedhorn-coupled bolometers cooled to 0.3 K. The photometer has a field of view of 4' x 8', observed simultaneously in the three spectral bands. Its main operating mode is scan-mapping, whereby the field of view is scanned across the sky to achieve full spatial sampling and to cover large areas if desired. The spectrometer has an approximately circular field of view with a diameter of 2.6'. The spectral resolution can be adjusted between 1.2 and 25 GHz by changing the stroke length of the FTS scan mirror. Its main operating mode involves a fixed telescope pointing with multiple scans of the FTS mirror to acquire spectral data. For extended source measurements, multiple position offsets are implemented by means of an internal beam steering mirror to achieve the desired spatial sampling and by rastering of the telescope pointing to map areas larger than the field of view. The SPIRE instrument consists of a cold focal plane unit located inside the Herschel cryostat and warm electronics units, located on the spacecraft Service Module, for instrument control and data handling. Science data are transmitted to Earth with no on-board data compression, and processed by automatic pipelines to produce calibrated science products. The in-flight performance of the instrument matches or exceeds predictions based on pre-launch testing and modelling: the photometer sensitivity is comparable to or slightly better than estimated pre-launch, and the spectrometer sensitivity is also better by a factor of 1.5-2.

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.

The Rosette molecular cloud is promoted as the archetype of a triggered star-formation site. This is mainly due to its morphology, because the central OB cluster NGC 2244 has blown a circular-shaped cavity into the cloud and the expanding HII-region now interacts with the cloud. Studying the spatial distribution of the different evolutionary states of all star-forming sites in Rosette and investigating possible gradients of the dust temperature will help to test the 'triggered star-formation' scenario in Rosette. We use continuum data obtained with the PACS (70 and 160 micron) and SPIRE instruments (250, 350, 500 micron) of the Herschel telescope during the Science Demonstration Phase of HOBYS. Three-color images of Rosette impressively show how the molecular gas is heated by the radiative impact of the NGC 2244 cluster. A clear negative temperature gradient and a positive density gradient (running from the HII-region/molecular cloud interface into the cloud) are detected. Studying the spatial distribution of the most massive dense cores (size scale 0.05 to 0.3 pc), we find an age-sequence (from more evolved to younger) with increasing distance to the cluster NGC 2244. No clear gradient is found for the clump (size-scale up to 1 pc) distribution. The existence of temperature and density gradients and the observed age-sequence imply that star formation in Rosette may indeed be influenced by the radiative impact of the central NGC 2244 cluster. A more complete overview of the prestellar and protostellar population in Rosette is required to obtain a firmer result.

Aims: To explore the infrared and radio properties of one of the closest Galactic starburst regions. Methods: Images obtained with the Herschel Space Observatory at wavelengths of 70, 160, 250, 350, and 500 microns using the PACS and SPIRE arrays are analyzed and compared with radio continuum VLA data and 8 micron images from the Spitzer Space Telescope. The morphology of the far-infrared emission is combined with radial velocity measurements of millimeter and centimeter wavelength transitions to identify features likely to be associated with the W43 complex. Results: The W43 star-forming complex is resolved into a dense cluster of protostars, infrared dark clouds, and ridges of warm dust heated by massive stars. The 4 brightest compact sources with L > 1.5 x 10^4 Lsun embedded within the Z-shaped ridge of bright dust emission in W43 remain single at 4" (0.1 pc) resolution. These objects, likely to be massive protostars or compact clusters in early stages of evolution are embedded in clumps with masses of 10^3 to 10^4 Msun, but contribute only 2% to the 3.6 x 10^6 Lsun far-IR luminosity of W43 measured in a 16 by 16 pc box. The total mass of gas derived from the far-IR dust emission inside this region is ~10^6 Msun. Cometary dust clouds, compact 6 cm radio sources, and warm dust mark the locations of older populations of massive stars. Energy release has created a cavity blowing-out below the Galactic plane. Compression of molecular gas in the plane by the older HII region near G30.684-0.260 and the bipolar structure of the resulting younger W43 HII region may have triggered the current mini-star burst.

We present a preliminary analysis of the small-scale structure found in new 70-520 micron continuum maps of the Rosette molecular cloud (RMC), obtained with the SPIRE and PACS instruments of the Herschel Space Observatory. We find 473 clumps within the RMC using a new structure identification algorithm, with sizes up to ~1.0 pc in diameter. A comparison with recent Spitzer maps reveals that 371 clumps are "starless" (without an associated young stellar object), while 102 are "protostellar." Using the respective values of dust temperature, we determine the clumps have masses (M_C) over the range -0.75 <= log (M_C/M_sun) <= 2.50. Linear fits to the high-mass tails of the resulting clump mass spectra (CMS) have slopes that are consistent with those found for high-mass clumps identified in CO emission by other groups.

We present a power spectrum analysis of the Herschel-SPIRE observations of the Polaris flare, a high Galactic latitude cirrus cloud midway between the diffuse and molecular phases. The SPIRE images of the Polaris flare reveal for the first time the structure of the diffuse interstellar medium down to 0.01 parsec over a 10 square degrees region. These exceptional observations highlight the highly filamentary and clumpy structure of the interstellar medium even in diffuse regions of the map. The power spectrum analysis shows that the structure of the interstellar medium is well described by a single power law with an exponent of -2.7 +- 0.1 at all scales from 30" to 8 degrees. That the power spectrum slope of the dust emission is constant down to the SPIRE angular resolution is an indication that the inertial range of turbulence extends down to the 0.01 pc scale. The power spectrum analysis also allows the identification of a Poissonian component at sub-arcminute scales in agreement with predictions of the cosmic infrared background level at SPIRE wavelengths. Finally, the comparison of the SPIRE and IRAS 100 micron data of the Polaris flare clearly assesses the capability of SPIRE in maping diffuse emission over large areas.

The Herschel OB young stellar objects survey (HOBYS) has observed the Rosette molecular cloud, providing an unprecedented view of its star formation activity. These new far-infrared data reveal a population of compact young stellar objects whose physical properties we aim to characterise. We compiled a sample of protostars and their spectral energy distributions that covers the near-infrared to submillimetre wavelength range. These were used to constrain key properties in the protostellar evolution, bolometric luminosity, and envelope mass and to build an evolutionary diagram. Several clusters are distinguished including the cloud centre, the embedded clusters in the vicinity of luminous infrared sources, and the interaction region. The analysed protostellar population in Rosette ranges from 0.1 to about 15 Msun with luminosities between 1 and 150 Lsun, which extends the evolutionary diagram from low-mass protostars into the high-mass regime. Some sources lack counterparts at near- to mid-infrared wavelengths, indicating extreme youth. The central cluster and the Phelps & Lada 7 cluster appear less evolved than the remainder of the analysed protostellar population. For the central cluster, we find indications that about 25% of the protostars classified as Class I from near- to mid-infrared data are actually candidate Class 0 objects. As a showcase for protostellar evolution, we analysed four protostars of low- to intermediate-mass in a single dense core, and they represent different evolutionary stages from Class 0 to Class I. Their mid- to far-infrared spectral slopes flatten towards the Class I stage, and the 160 to 70um flux ratio is greatest for the presumed Class 0 source. This shows that the Herschel observations characterise the earliest stages of protostellar evolution in detail.

Our PACS and SPIRE images of the Aquila Rift and part of the Polaris Flare regions, taken during the science demonstration phase of Herschel discovered fascinating, omnipresent filamentary structures that appear to be physically related to compact cores. We briefly describe a new multi-scale, multi-wavelength source extraction method used to detect objects and measure their parameters in our Herschel images. All of the extracted starless cores (541 in Aquila and 302 in Polaris) appear to form in the long and very narrow filaments. With its combination of the far-IR resolution and sensitivity, Herschel directly reveals the filaments in which the dense cores are embedded; the filaments are resolved and have deconvolved widths of 35 arcsec in Aquila and 59 arcsec in Polaris (9000 AU in both regions). Our first results of observations with Herschel enable us to suggest that in general dense cores may originate in a process of fragmentation of complex networks of long, thin filaments, likely formed as a result of an interplay between gravity, interstellar turbulence, and magnetic fields. To unravel the roles of the processes, one has to obtain additional kinematic and polarization information; these follow-up observations are planned.

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.

In Herschel images of the Galactic plane and many star forming regions, a major factor limiting our ability to extract faint compact sources is cirrus confusion noise, operationally defined as the "statistical error to be expected in photometric measurements due to confusion in a background of fluctuating surface brightness". The histogram of the flux densities of extracted sources shows a distinctive faint-end cutoff below which the catalog suffers from incompleteness and the flux densities become unreliable. This empirical cutoff should be closely related to the estimated cirrus noise and we show that this is the case. We compute the cirrus noise directly, both on Herschel images from which the bright sources have been removed and on simulated images of cirrus with statistically similar fluctuations. We connect these direct estimates with those from power spectrum analysis, which has been used extensively to predict the cirrus noise and provides insight into how it depends on various statistical properties and photometric operational parameters. We report multi-wavelength power spectra of diffuse Galactic dust emission from Hi-GAL observations at 70 to 500 microns within Galactic plane fields at l= 30 degrees and l= 59 degrees. We find that the exponent of the power spectrum is about -3. At 250 microns, the amplitude of the power spectrum increases roughly as the square of the median brightness of the map and so the expected cirrus noise scales linearly with the median brightness. Generally, the confusion noise will be a worse problem at longer wavelengths, because of the combination of lower angular resolution and the rising power spectrum of cirrus toward lower spatial frequencies, but the photometric signal to noise will also depend on the relative spectral energy distribution of the source compared to the cirrus.

Context: Sh2-104 is a Galactic H ii region with a bubble morphology, detected at optical and radio wavelengths. It is considered the first observational confirmation of the collect-and-collapse model of triggered star-formation. Aims: We aim to analyze the dust and gas properties of the Sh2-104 region to better constrain its effect on local future generations of stars. In addition, we investigate the relationship between the dust emissivity index {\beta} and the dust temperature, T_dust. Methods: Using Herschel PACS and SPIRE images at 100, 160, 250, 350 and 500 {\mu}m we determine T_dust and {\beta} throughout Sh2-104, fitting the spectral energy distributions (SEDs) obtained from aperture photometry. With the SPIRE Fourier transform spectrometer (FTS) we obtained spectra at different positions in the Sh2-104 region. We detect J-ladders of CO and 13CO, with which we derive the gas temperature and column density. We also detect proxies of ionizing flux as the [NII] 3P1-3P0 and [CI] 3P2-3P1 transitions. Results: We find an average value of {\beta} ~ 1.5 throughout Sh2-104, as well as a T dust difference between the photodissociation region (PDR, ~ 25 K) and the interior (~ 40 K) of the bubble. We recover the anti-correlation between {\beta} and dust temperature reported numerous times in the literature. The relative isotopologue abundances of CO appear to be enhanced above the standard ISM values, but the obtained value is very preliminary and is still affected by large uncertainties.

We present far-infrared spectra and maps of the DR21 molecular cloud core between 196 and 671 microns, using the Herschel-SPIRE spectrometer. Nineteen molecular lines originating from CO, 13CO, HCO+ and H2O, plus lines of [N II] and [CI] were recorded, including several transitions not previously detected. The CO lines are excited in warm gas with Tkin ~ 125 K and nH2 ~ 7 x 10^4 cm-3, CO column density N(CO) ~ 3.5 x 10^18 cm^-2 and a filling factor of ~ 12%, and appear to trace gas associated with an outflow. The rotational temperature analysis incorporating observations from ground-based telescopes reveals an additional lower excitation CO compoment which has a temperature ~ 78 K and N(CO) ~ 4.5 x 10^21 cm^-2. Astronomy & Astrophysics HERSCHEL special Issue, in press.

The origin and possible universality of the stellar initial mass function (IMF) is a major issue in astrophysics. One of the main objectives of the Herschel Gould Belt Survey is to clarify the link between the prestellar core mass function (CMF) and the IMF. We present and discuss the core mass function derived from Herschel data for the large population of prestellar cores discovered with SPIRE and PACS in the Aquila Rift cloud complex at d ~ 260 pc. We detect a total of 541 starless cores in the entire ~11 deg^2 area of the field imaged at 70-500 micron with SPIRE/PACS. Most of these cores appear to be gravitationally bound, and thus prestellar in nature. Our Herschel results confirm that the shape of the prestellar CMF resembles the stellar IMF, with much higher quality statistics than earlier submillimeter continuum ground-based surveys.

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.

The Polaris Flare cloud region contains a great deal of extended emission. It is at high declination and high Galactic latitude. It was previously seen strongly in IRAS Cirrus emission at 100 microns. We have detected it with both PACS and SPIRE on Herschel. We see filamentary and low-level structure. We identify the five densest cores within this structure. We present the results of a temperature, mass and density analysis of these cores. We compare their observed masses to their virial masses, and see that in all cases the observed masses lie close to the lower end of the range of estimated virial masses. Therefore, we cannot say whether they are gravitationally bound prestellar cores. Nevertheless, these are the best candidates to be potentialprestellar cores in the Polaris cloud region.

We use the SPIRE Fourier-Transform Spectrometer (FTS) on-board the ESA Herschel Space Telescope to analyse the submillimetre spectrum of the Ultra-compact HII region G29.96-0.02. Spectral lines from species including 13CO, CO, [CI], and [NII] are detected. A sparse map of the [NII] emission shows at least one other HII region neighbouring the clump containing the UCHII. The FTS spectra are combined with ISO SWS and LWS spectra and fluxes from the literature to present a detailed spectrum of the source spanning three orders of magnitude in wavelength. The quality of the spectrum longwards of 100 {\mu}m allows us to fit a single temperature greybody with temperature 80.3\pm0.6K and dust emissivity index 1.73\pm0.02, an accuracy rarely obtained with previous instruments. We estimate a mass of 1500 Msol for the clump containing the HII region. The clump's bolometeric luminosity of 4 x 10^6 Lsol is comparable to, or slightly greater than, the known O-star powering the UCHII region.

Aims. To follow the species chemistry arising in diverse sources of the Galaxy with Herschel. Methods. SPIRE FTS sparse sampled maps of the Orion bar & compact HII regions G29.96-0.02 and G32.80+0.19 have been analyzed. Results. Beyond the wealth of atomic and molecular lines detected in the high-resolution spectra obtained with the FTS of SPIRE in the Orion Bar, one emission line is found to lie at the position of the fundamental rotational transition of CH+ as measured precisely in the laboratory (Pearson & Drouion 2006). This coincidence suggests that it is the first detection of the fundamental rotational transition of CH+. This claim is strengthened by the observation of the lambda doublet transitions arising from its relative, CH, which are also observed in the same spectrum. The broad spectral coverage of the SPIRE FTS allows for the simultaneous measurement of these closely related chemically species, under the same observing conditions. The importance of these lines are discussed and a comparison with results obtained from models of the Photon Dominated Region (PDR) of Orion are presented. The CH+ line also appears in absorption in the spectra of the two galactic compact HII regions G29.96-0.02 and G32.80+0.19, which is likely due to the presence of CH+ in the the Cold Neutral Medium of the galactic plane. These detections will shed light on the formation processes and on the existence of CH+, which are still outstanding questions in astrophysics.

We present a first study of the star-forming compact dust condensations revealed by Herschel in the two 2 \times 2 \degr Galactic Plane fields centered at [l;b] = [30\degr; 0 \degr] and [l;b] = [59\degr; 0 \degr], respectively, and observed during the Science Demonstration Phase for the Herschel infrared Galactic Plane survey (Hi-GAL) Key-Project. Compact source catalogs extracted for the two fields in the five Hi-GAL bands (70, 160, 250, 350 and 500 $\mu$m) were merged based on simple criteria of positional association and spectral energy distribution (SED) consistency into a final catalog which contains only coherent SEDs with counterparts in at least three adjacent Herschel bands. These final source lists contain 528 entries for the l = 30\degr field, and 444 entries for the l = 59\degr field. The SED coverage has been augmented with ancillary data at 24 $\mu$m and 1.1 mm. SED modeling for the subset of 318 and 101 sources (in the two fields, respectively) for which the distance is known was carried out using both a structured star/disk/envelope radiative transfer model and a simple isothermal grey-body. Global parameters like mass, luminosity, temperature and dust properties have been estimated. The Lbol/Menv ratio spans four orders of magnitudes from values compatible with the pre-protostellar phase to embedded massive zero-age main sequence stars. Sources in the l = 59\degr field have on average lower L/M, possibly outlining an overall earlier evolutionary stage with respect to the sources in the l = 30\degr field. Many of these cores are actively forming high-mass stars, although the estimated core surface densities appear to be an order of magnitude below the 1 g cm$^{-2}$ critical threshold for high-mass star formation.

Context. RCW 120 is a well-studied, nearby Galactic HII region with ongoing star formation in its surroundings. Previous work has shown that it displays a bubble morphology at mid-infrared wavelengths and has a massive layer of collected neutral material seen at sub-mm wavelengths. Given the well-defined photo-dissociation region (PDR) boundary and collected layer, it is an excellent laboratory to study the "collect and collapse" process of triggered star formation. Using Herschel Space Observatory data at 100, 160, 250, 350, and 500 micron, in combination with Spitzer and APEX-LABOCA data, we can for the first time map the entire spectral energy distribution of an HII region at high angular resolution. Aims. We seek a better understanding of RCW120 and its local environment by analysing its dust temperature distribution. Additionally, we wish to understand how the dust emissivity index, beta, is related to the dust temperature. Methods. We determine dust temperatures in selected regions of the RCW 120 field by fitting their spectral energy distribution (SED), derived using aperture photometry. Additionally, we fit the SED extracted from a grid of positions to create a temperature map. Results. We find a gradient in dust temperature, ranging from >30 K in the interior of RCW 120, to ~20K for the material collected in the PDR, to ~10K toward local infrared dark clouds and cold filaments. Our results suggest that RCW 120 is in the process of destroying the PDR delineating its bubble morphology. The leaked radiation from its interior may influence the creation of the next generation of stars. We find support for an anti-correlation between the fitted temperature and beta, in rough agreement with what has been found previously. The extended wavelength coverage of the Herschel data greatly increases the reliability of this result.

It has been shown that by means of different physical mechanisms the expansion of HII regions can trigger the formation of new stars of all masses. This process may be important to the formation of massive stars but has never been quantified in the Galaxy. We use Herschel-PACS and -SPIRE images from the Herschel Infrared survey of the Galactic plane, Hi-GAL, to perform this study. We combine the Spitzer-GLIMPSE and -MIPSGAL, radio-continuum and sub-millimeter surveys such as ATLASGAL with Hi-GAL to study Young Stellar Objects (YSOs) observed towards Galactic HII regions. We select a representative HII region, N49, located in the field centered on l=30 degr observed as part of the Hi-GAL Science Demonstration Phase, to demonstrate the importance Hi-GAL will have to this field of research. Hi-GAL PACS and SPIRE images reveal a new population of embedded young stars, coincident with bright ATLASGAL condensations. The Hi-GAL images also allow us, for the first time, to constrain the physical properties of the newly formed stars by means of fits to their spectral energy distribution. Massive young stellar objects are observed at the borders of the N49 region and represent second generation massive stars whose formation has been triggered by the expansion of the ionized region. Hi-GAL enables us to detect a population of young stars at different evolutionary stages, cold condensations only being detected in the SPIRE wavelength range. The far IR coverage of Hi-GAL strongly constrains the physical properties of the YSOs. The large and unbiased spatial coverage of this survey offers us a unique opportunity to lead, for the first time, a global study of star formation triggered by HII regions in our Galaxy.

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.

Hi-GAL, the Herschel infrared Galactic Plane Survey, is an Open Time Key Project of the Herschel Space Observatory. It will make an unbiased photometric survey of the inner Galactic Plane by mapping a two-degree wide strip in the longitude range |l|<60 degrees in five wavebands between 70um and 500um. The aim of Hi-GAL is to detect the earliest phases of the formation of molecular clouds and high-mass stars and to use the optimum combination of Herschel wavelength coverage, sensitivity, mapping strategy and speed to deliver a homogeneous census of star-forming regions and cold structures in the interstellar medium. The resulting representative samples will yield the variation of source temperature, luminosity, mass and age in a wide range of Galactic environments at all scales from massive YSOs in protoclusters to entire spiral arms, providing an evolutionary sequence for the formation of intermediate and high-mass stars. This information is essential to the formulation of a predictive global model of the role of environment and feedback in regulating the star-formation process. Such a model is vital to understanding star formation on galactic scales and in the early Universe. Hi-GAL will also provide a science legacy for decades to come with incalculable potential for systematic and serendipitous science in a wide range of astronomical fields, enabling the optimum use of future major facilities such as JWST and ALMA.

(Abridged) Studying continuum emission from interstellar dust is essential to locating and characterizing the highest density regions in the interstellar medium. In particular, the early stages of massive star formation remain poorly understood. Our goal is to produce a large-scale, systematic database of massive pre- and proto-stellar clumps in the Galaxy, to understand how and under what conditions star formation takes place. A well characterized sample of star-forming sites will deliver an evolutionary sequence and a mass function of high-mass, star-forming clumps. This systematic survey at submm wavelengths also represents a preparatory work for Herschel and ALMA. The APEX telescope is ideally located to observe the inner Milky Way. The Large APEX Bolometer Camera (LABOCA) is a 295-element bolometer array observing at 870 microns, with a beam size of 19". Taking advantage of its large field of view (11.4') and excellent sensitivity, we started an unbiased survey of the Galactic Plane, with a noise level of 50-70 mJy/beam: the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL). As a first step, we covered 95 sq. deg. These data reveal 6000 compact sources brighter than 0.25 Jy, as well as extended structures, many of them filamentary. About two thirds of the compact sources have no bright infrared counterpart, and some of them are likely to correspond to the precursors of (high-mass) proto-stars or proto-clusters. Other compact sources harbor hot cores, compact HII regions or young embedded clusters. Assuming a typical distance of 5 kpc, most sources are clumps smaller than 1 pc with masses from a few 10 to a few 100 M_sun. In this introductory paper, we show preliminary results from these ongoing observations, and discuss the perspectives of the survey.

[Abridged] Twenty-six high-luminosity IRAS sources believed to be collection of stars in the early phases of high-mass star formation have been observed in the NIR (J, H, K) to characterize the clustering properties of their young stellar population and gain insight into the initial conditions of star formation in these clusters (Initial Mass Function [IMF], Star Formation History [SFH]), and to deduce mean values for cluster ages. K luminosity functions (KLFs) are compared with simulated ones from a model that generates populations of synthetic clusters starting from assumptions on the IMF, the SFH, and the Pre-MS evolution, and using the average properties of the observed clusters as boundary conditions Twenty-two sources show evidence of clustering from a few up to several tens of objects, and a median cluster radius of 0.7 pc. A considerable number of cluster members present an infrared excess characteristic of young Pre-Main-Sequence objects. We find that the median stellar age ranges between 2.5 10^5 and 5 10^6 years, with evidence of an age spread of the same entity within each cluster. We also find evidence that older clusters tend to be smaller in size, in line with the fact that our clusters are on average larger than those around relatively older Herbig Ae/Be stars. The relationship of the mass of the most massive star in the cluster with both the clusters richness and their total stellar mass suggest that our modeled clusters may not be consistent with them resulting from random sampling of the IMF. Our results are consistent with a star formation which takes place continuously over a period of time which is longer than a typical crossing time.

Determining the timeline for the formation of massive YSOs requires the identification and characterisation of all the phases that a massive forming YSO undergoes. It is of particular interest to verify the observability of the phase in which the object is rapidly accreting while not yet igniting the fusion of hydrogen that marks the arrival on the ZAMS. One of the candidate prototypical objects for this phase is Mol160/IRAS23385+6053, which previous studies suggest it could be in a pre-Hot Core stage. We further investigate this issue by means of Spitzer imaging and spectroscopy in the 5-70 micron range. The dense core of Mol160/IRAS23385+6053, which up to now had only been detected at submm and mm wavelenghts has been revealed for the first time at 24 and 70 micron by Spitzer. The complete 24 micron -3.4 mm continuum cannot be fitted with a standard model of a Zero-Age Main-Sequence (ZAMS) star embedded in an envelope. A simple greybody fit yields a mass of 220 solar masses. The luminosity is slightly in excess of 3000 solar luminosities, which is a factor of 5 less than previous estimates when only IRAS fluxes were available between 20 and 100 micron. The source is under-luminous by the same factor with respect to UCHII regions or Hot-Cores of similar circumstellar mass, and simple models show that this is compatible with an earlier evolutionary stage. Spectroscopy between 5 and 40 microns revelas typical PDR/PIR conditions, where the required UV illumination may be provided by other sources revealed at 24 microns in the same region, and which can be plausibly modeled as moderately embedded intermediate-mass ZAMS stars. Our results strengthen the suggestion that the central core in Mol160/IRAS23385+6053 is a massive YSO actively accreting from its circumstellar envelope and which did not yet begin hydrogen fusion.

Context: The evolution of a young protocluster depends on the relative spatial distributions and dynamics of both stars and gas. Aims: We study the distribution and properties of the gas and stars surrounding the luminous (10^4 L_sun) protocluster IRAS 18511+0146. Methods: IRAS 18511+0146 and the cluster associated with it has been investigated using the sub-millimetre (JCMT-SCUBA), infrared (Spitzer-MIPSGAL, Spitzer-GLIMPSE, Palomar) and radio (VLA) continuum data. Cluster simulations have been carried out in order to understand the properties of clusters as well as to compare with the observations. Results: The central most obscured part of the protocluster coincident with the compact sub-millimetre source found with SCUBA is responsible for at least 2/3 of the total luminosity. A number of cluster members have been identified which are bright in mid infrared and show rising (near to mid infrared) spectral energy distributions suggesting that these are very young stellar sources. In the mid infrared 8.0 micron image, a number of filamentary structures and clumps are detected in the vicinity of IRAS 18511+0146. Conclusions: Based on the luminosity and cluster size as well as on the evolutionary stages of the cluster members, IRAS 18511+0146 is likely to be protocluster with the most massive object being a precursor to a Herbig type star.

We have measured the deuterium fractionation (through the column density ratio N(N2D+)/N(N2H+)) and the CO depletion factor (ratio between expected and observed CO abundance) in a sample of 10 high-mass protostellar candidates, in order to understand whether the earliest evolutionary stages of high-mass stars have chemical characteristics similar to those of low-mass ones. The observations were carried out with the IRAM-30m telescope and the JCMT. We have detected N2D+ emission in 7 of the 10 sources of our sample, and found an average value N(N2D+)/N(N2H+)~0.015. This value is 3 orders of magnitude larger than the interstellar D/H ratio, indicating the presence of cold and dense gas, in which the physical-chemical conditions are similar to those observed in low-mass pre-stellar cores. Also, the integrated CO depletion factors show that in the majority of the sources the expected CO abundances are larger than the observed values, with a median ratio of 3.2. In principle, the cold gas that gives origin to the N2D+ emission can be the remnant of the massive molecular core in which the high-mass (proto-)star was born, not yet heated up by the central object. If so, our results indicate that the chemical properties of the clouds in which high-mass stars are born are similar to their low-mass counterparts. Alternatively, this cold gas can be located into one (or more) starless core (cores) near the protostellar object. Due to the poor angular resolution of our data, we cannot decide which is the correct scenario.

In an ongoing effort to identify and study high-mass protostellar candidates we have observed in various tracers a sample of 235 sources selected from the IRAS Point Source Catalog, mostly with dec < -30 deg, with the SEST antenna at millimeter wavelengths. The sample contains 142 Low sources and 93 High, which are believed to be in different evolutionary stages. Both sub-samples have been studied in detail by comparing their physical properties and morphologies. Massive dust clumps have been detected in all but 8 regions, with usually more than one clump per region. The dust emission shows a variety of complex morphologies, sometimes with multiple clumps forming filaments or clusters. The mean clump has a linear size of ~0.5 pc, a mass of ~320 Msolar for a dust temperature Td=30 K, an H_2 density of 9.5E5 cm-3, and a surface density of 0.4 g cm-2. The median values are 0.4 pc, 102 Msolar, 4E4 cm-3, and 0.14 g cm-2, respectively. The mean value of the luminosity-to-mass ratio, L/M ~99 Lsolar/Msolar, suggests that the sources are in a young, pre-ultracompact HII phase. We have compared the millimeter continuum maps with images of the mid-IR MSX emission, and have discovered 95 massive millimeter clumps non-MSX emitters, either diffuse or point-like, that are potential prestellar or precluster cores. The physical properties of these clumps are similar to those of the others, apart from the mass that is ~3 times lower than for clumps with MSX counterpart. Such a difference could be due to the potential prestellar clumps having a lower dust temperature. The mass spectrum of the clumps with masses above M ~100 Msolar is best fitted with a power-law dN/dM proportional to M-alpha with alpha=2.1, consistent with the Salpeter (1955) stellar IMF, with alpha=2.35.

(Abridged) We have observed CS and C17O lines, and 1.2 mm cont. emission towards a sample of 130 high-mass protostellar candidates with DEC<-30 deg. This is the first step of the southern extension of a project started more than a decade ago aimed at the identification of massive protostellar candidates. We selected from the IRAS PSC 429 sources which potentially are compact molecular clouds. The sample is divided into two groups: the 298 sources with [25-12]>0.57 and [60-12]>1.30 we call 'High' sources, the remaining 131 we call 'Low' sources. In this paper, we check the association with dense gas and dust in 130 'Low' sources. We find a detection rate of ca. 85% in CS, demonstrating a tight association with dense molecular clumps. Among the sources detected in CS, ca. 76% have also been detected in C17O and ca. 93% in the 1.2 mm cont. Mm-cont. maps show the presence of clumps with diameters 0.2-2 pc and masses from a few Msun to 10^5 Msun; H2 volume densities lie between ca. 10^{4.5} and 10^{5.5} cm^{-3}. The L(bol) are 10^3-10^6 Lsun, consistent with embedded high-mass objects. Based on our results and those found in the literature for other samples, we conclude that our sources are massive objects probably in a stage prior to the formation of an HII region. We propose a scenario in which 'High' and 'Low' sources are both made of a massive clump hosting a high-mass protostellar candidate and a nearby stellar cluster. The difference might be due to the fact that the IRAS 12mu flux, the best discriminant between the two groups, is dominated by the emission from the cluster in 'Lows' and from the massive protostellar object in 'Highs'.

We present observations of continuum and molecular lines towards the protostar candidates IRAS 21307+5049 and IRAS 22172+5549. Single-dish maps in the 12CO(2--1), C18O(2--1), HCO+(1--0) lines and sub-mm continuum are compared with interferometric maps in the 12CO(1--0) line and 3 mm continuum, and with mid- and near-infrared images. A plausible interpretation of our data based on the continuum maps and spectral energy distributions is that two components are present: a compact molecular core, responsible for the continuum emission at wavelengths longer than ~25 \mum, and a cluster of stars located close to the center of the core, but not spatially coincident with it, responsible for the emission at shorter wavelengths. The core is approximately located at the center of the associated molecular outflow, detected for both sources in the 12CO(1--0) and (2--1) lines. The cores have masses of ~50 M_{\odot}, and luminosities of ~10^{3} L_{\odot}. The outflows parameters are consistent with those typically found in high-mass young stellar objects. Our results support the hypothesis that in these sources the luminosity is dominated by accretion rather than by nuclear burning. We conclude that the sources embedded inside the cores are likely protostars with mass ~5-8 M_{\odot}.

We present the results of a multi-line and continuum study towards the source IRAS 23385+6053,performed with the IRAM-30m telescope, the Plateau de Bure Interferometer, the Very Large Array Interferometer and the James Clerk Maxwell Telescope. The new results confirm our earlier findings, namely that IRAS 23385+6053 is a good candidate high-mass protostellar object, precursor of an ultracompact H$_{II}$ region. The source is roughly composed of two regions: a molecular core $\sim0.03\div0.04$ pc in size, with a temperature of $\sim40$ K and an H$_{2}$ volume density of the order of 10$^{7}$ cm$^{-3}$, and an extended halo of diameter $\leq$0.4 pc, with an average kinetic temperature of $\sim 15$ K and H$_{2}$ volume density of the order of 10$^{5}$ cm$^{-3}$. The core temperature is much smaller than what is typically found in molecular cores of the same diameter surrounding massive ZAMS stars. We deduce that the core luminosity is between 150 and $1.6\times10^{4}L_{\odot}$, and we believe that the upper limit is near the ``true'' source luminosity. Moreover, by comparing the H$_{2}$ volume density obtained at different radii from the IRAS source, we find that the halo has a density profile of the type $n_{\rm H_{2}}\propto r^{-2.3}$. This suggests that the source is gravitationally unstable. Finally, we demonstrate that the temperature at the core surface is consistent with a core luminosity of $10^3 L_{\odot}$ and conclude that we might be observing a protostar still accreting material from its parental cloud, whose mass at present is $\sim 6 M_{\odot}$.

The HH 1/2 system has been observed with the two spectrometers on board the Infrared Space Observatory. Diffuse [CII] 158um emission indicate the presence of a Photo-Dissociation Region (PDR) which is only in part due to the external Far-UV irradiation from the nearby Orion Nebula. Additional irradiation must be originated locally and we show that the FUV field produced in the recombination regions behind the shocks traced by the HH objects is sufficient to induce a PDR at the flow cavity walls. The analysis of [OI], [SiII]34.8um and [NeII] 12.8um lines suggests shock velocities v=100-140 km/s with pre-shock densities between 100 and 1000 per cc. H2 pure rotational lines trace a T=600 K gas which is likely to be warmed up in slow, planar, molecular shocks. The coexistence in the areas subtended by the instrumental beamsizes of shocks at different velocities is consistent with the bow-shock morphology of the HH 1 and, to a lesser extent, HH 2 as traced by optical images.

We present molecular line observations made with the IRAM 30-m telescope of the immediate surroundings of a sample of 11 candidate high-mass protostars. These observations are part of an effort to clarify the evolutionary status of a set of objects which we consider to be precursors of UC HII regions. We find that the present data are consistent with such an interpretation. The sources have FIR luminosities indicating the presence of B2.5 to O8.5V0 stars, are associated with dense gas and dust, have (sub-)mm continuum spectra indicating temperatures of ca. 30K, and have no detectable radio continuum emission. All observed sources are associated with well-defined molecular clumps. Masses, sizes, and other parameters depend on the tracer used, but typically the cores have average diameters of ca. 0.5-1 pc, and masses of a few tens to a few thousand solar masses. Compared to a similar analysis of sources which are mostly identifyable with UC HII regions, the present sample has molecular clumps that are more massive, larger, cooler, and less turbulent. They also tend to have a smaller ratio of virial-to-luminous mass, indicating they are less dynamically stable than their counterparts embedded in the UC HII-like sources. The large sizes suggest these clumps should still undergo substantial contraction. The lower temperatures and small linewidths are also expected in objects in an earlier evolutionary state. We find that the masses of the molecular clumps associated with our objects increase with Lfir [M(clump) propto Lfir^(1.17)], and that there is a (weak) relation between the clump mass and the mass of the embedded protostellar object M(proto) propto M(clump)^(0.30).

The Cepheus E outflow has been studied in the mid and far infrared using the ISO CAM and LWS instruments, and at millimetric wavelengths using OVRO. In the near and mid-IR, its morphology is similar to that expected for a jet driven outflow, where the leading bow shocks entrain and accelerate the surrounding molecular gas. As expected, fine structure atomic/ionic emission lines arise from the bow shocks, at both the Mach Disk and the stagnation tip, where J-shocks are dominant. The H2, H2O and CO molecular emission could arise further `downstream' at the bow shock wings where the shocks (v = 8-35 km/s) are oblique and more likely to be C-type. The 13CO emission arises from entrained molecular gas and a compact high velocity emission is observed, together with an extended low velocity component that almost coincides spatially with the H2 near-IR emission. The millimetric continuum emission shows two sources. We identify one of them with IRAS 23011+6126, postulating is the driver of the Cepheus E outflow; the other, also an embedded source, is likely to be driving one of other outflows observed in the region.

The two spectrometers on board the Infrared Space Observatory were used to observe the Herbig-Haro objects HH 80, 81 and 80N, as well as their candidate exciting source IRAS 18162-2048. The fine structure lines of [OI]63um, [OI]145um and [CII]158um are detected everywhere, while [NII]122um and [OIII]88um are only detected toward the HH objects; line ratios confirm for the first time the collisionally excited HH nature of HH 80N. No molecular line is detected in any of the observed positions. We use a full shock code to diagnose shock velocities vs~100 km/s toward the HH objects, as expected from the optical spectroscopy. Since proper motions suggest velocities in excess of 600 km/s, the HH objects probably represent the interface between two flow components with velocity differing by ~vs. Aside from the flow exciting source, the [CII]158um line is everywhere brighter than the [OI]63um line, indicating the presence of a Photo-Dissociation Region (PDR) all along the flow. Continuum emission from the HH objects and from other positions along the flow is only detected longword of ~50 micron, and its proportionality to the [CII]158um line flux suggests it is PDR in origin. We propose that the FUV continuum irradiated by the HH objects and the jet is responsible for the generation of a PDR at the walls of the flow cavity. We develop a very simple model which strengthens the plausibility of this hypothesis.

ISO-LWS mapping observations of the Serpens molecular cloud core are presented. The spectral range is 50 - 200 micron and the map size is 8' X 8'. These observations suffer from severe source confusion at FIR wavelengths and we employ a Maximum Likelihood Method for the spectro-spatial deconvolution. The strong and fairly isolated source SMM1 FIRS1 presented a test case, whose modelled spectral energy distribution (SED), within observational errors, is identical to the observed one. The model results for the other infrared and submillimetre sources are therefore likely to represent their correct SEDs. Simulations demonstrating the reliability and potential of the developed method support this view. For the majority of sources the peak of the SEDs is found within the spectral range of the LWS and derived temperatures are generally higher (>30K) than have been found by earlier deconvolution attempts using IRAS data. SMM sizes are found to be only a few arcsec in diameter. In addition, the SMMs are generally optically thick even at LWS wavelengths. Self-consistent radiative transfer calculations indicates that the SMMs are optically thick out to even longer wavelengths. Models were calculated for five sources, for which sufficient data were available, viz. SMM 1, 2, 3, 4 and 9. These models are optically thick out to millimetre wavelengths (wavelength of unit optical depth 900 to 1400 micron). Envelope masses for these SMMs are in the range 2-6 Msun, which is of course considerably more massive than estimates based on the optically thin assumption. The luminosities are in the range 10-70 Lsun, suggesting the formation of low-mass to intermediate mass stars, so that the existence of such massive envelopes argues for extreme youth of the SMMs in the Serpens cloud core.

We have used the two spectrometers on the Infrared Space Observatory to observe the HH 7-11 flow, its red-shifted counterpart, and the candidate exciting source SVS 13, in the star formation region NGC 1333. We detect atomic ([OI]63um and 145um, [SiII]34.8um, [CII]158um) and molecular (H_2, CO, H_2O) lines at various positions along the bipolar flow. Most of the observed lines can be explained in terms of shock-excited emission. In particular, our analysis shows that dissociative (J-type) and non-dissociative (C-type) shocks are simultaneously present everywhere along both lobes of the flow. We confirm the low-excitation nature of the Herbig-Haro nebulosities, with shock velocities <40-50 km/s. Toward both lobes of the outflow we find pre-shock densities of n_0~10^4 cm-3 for both the J and C components, implying B_0~100uG. In the central region of the flow, close to the exciting source, the pre-shock density deduced for the C-shock component is n_0~10^5 cm-3, suggesting a magnetic field ~3 times stronger. We propose that the deficiency of gas-phase water in the post C-shock regions is due to freezing onto warm grains processed through the J-shock front and traveling along the magnetic field lines. The total observed cooling from the dissociative shock components is consistent with the power lost by a slow molecular outflow accelerated by a fast neutral HI wind. Finally, the skin of the cloud seen in projection toward the flow appears to be weakly photo-ionised by BD +30 549, the dominant illuminating source of the NGC 1333 reflection nebula.

The James Clerk Maxwell Telescope has been used to obtain submillimeter and millimeter continuum photometry of a sample of 30 IRAS sources previously studied in molecular lines and centimeter radio continuum. All the sources have IRAS colours typical of very young stellar objects (YSOs) and are associated with dense gas. In spite of their high luminosities (L>10000 solar units), only ten of these sources are also associated with a radio counterpart. In 17 cases we could identify a clear peak of millimeter emission associated with the IRAS source, while in 9 sources the millimeter emission was extended or faint and a clear peak could not be identified; upper limits were found in 4 cases only. Using simple greybody fitting model to the observed SED, we derive global properties of the circumstellar dust. The dust temperature varies from 24 K to 45 K, while the exponent of the dust emissivity vs frequency power-law spans a range 1.56<beta<2.38, characteristic of silicate dust; total circumstellar masses range up to more than 500 solar masses. We find that for sources with comparable luminosities, the total column densities derived from the dust masses do not distinguish between sources with and without radio counterpart. We interpret this result as an indication that dust does not play a dominant role in inhibiting the formation of the HII region. We examine several scenarios for their origin in terms of newborn ZAMS stars and although most of these fail to explain the observations, we cannot exclude that these sources are young stars already on the ZAMS with modest residual accretion that quenches the expansion of the HII region. Finally, we consider the possibility that the IRAS sources are high-mass pre-ZAMS (or pre-H-burning) objects deriving most of the emitted luminosity from accretion.

We report the detection of the 62 micron feature of crystalline water ice in emission towards the bow-shaped Herbig-Haro object HH 7. Significant amounts of far infrared continuum emission are also detected between 10 and 200 micron, so that Herbig-Haro objects cease to be pure emission-line objects at FIR wavelengths. The formation of crystalline water ice mantles requires grain temperatures T > 100 K at the time of mantle formation, suggesting that we are seeing material processed by the HH 7 shock front. The deduced ice mass is \~2e-5 solar masses corresponding to a water column density N~ 1e18 cm-3; an estimate of the [H2O]/[H] abundance yields values close to the interstellar gas-phase oxygen abundance. The relatively high dust temperature and the copious amounts of gas-phase water needed to produce the observed quantity of crystalline water ice, suggest a scenario where both dissociative and non-dissociative shocks co-exist. The timescale for ice mantle formation is of the order of ~400 years, so that the importance of gas-phase water cooling as a shock diagnostic may be greatly diminished.

IRAS 23385+6053 is a Young Stellar Object with luminosity 16,000 solar luminosities at a kinematic distance of 4.9 kpc. This candidate precursor of an ultracompact HII region is associated with a millimeter source detected at the JCMT but is undetected at centimeter wavelengths with the VLA. We observed this source with the OVRO millimeter array at 3.4 mm in the continuum, HCO+(1-0), H13CO+(1-0) and SiO(v=0, 2-1) line emission, and with CAM aboard ISO at 6.75 and 15 um. The IRAS source is coincident with a 3.4 mm compact (r~0.048 pc) and massive (M~370 solar masses) core, which is undetected at 15 um to a 3 sigma level of 6 mJy; this is compatible with the derived H2 column density of 2x10(24) and the estimated visual extinction of about 2000 mag. We find L(submm)/L(bol)>10(-3) and M(env)/M(star)>>1, typical of Class 0 objects. The source is also associated with a compact outflow characterized by a size similar to the core radius, a dynamical timescale of about 7,000 years, and a mass loss rate of about 10(-3) solar masses/year. The axis of the outflow is oriented nearly perpendicular to the plane of the sky, ruling out the possibility that the non-detection at 15 um is the result of a geometric effect. All these properties suggest that IRAS 23385+6053 is the first example of a bona fide massive Class 0 object.

The physical characteristics of Cepheus E (Cep E) `embedded' outflow are analyzed using ISOCAM images in the v=0-0 S(5) 6.91 um and S(3) 9.66 um molecular hydrogen lines. We find that the morphology of the Cep E outflow in the ground vibrational H2 lines is similar to that of the near infrared v=1-0 2.12 um line. At these mid-IR wavelengths, we do not detect the second H2 outflow which is almost perpendicular to Cep E 2.121 um flow or traces of H2 emission along the second 12CO J = 2-1 outflow at 52 degrees angle, down to a surface brightness of 12 - 46 uJy/arcsec square. We do detect at 6.91 um the likely source of the main H2 and CO outflows, IRAS 23011+6126, and show that the source is easily seen in all IRAS bands using HiRes images. The source is not detected at 9.66 um, but we think this agrees with the interstellar extinction curve which has a minimum at 7 um, but rises a 9.7 um due to the strong absorption silicate feature, enhanced in this case by a cocoon surrounding the Class 0 object. This idea is supported by our models of the spectral energy distribution (SED) of the central object. The models assume that the main source of opacity is due to bare silicates and our best fit for the SED yields a total mass of envelope of 17 solar masses and a dust temperature of 18 K.