Context. The mechanism of formation of massive stars is still a matter of debate. It is not yet clear if it can be considered to be a scaled-up analogue of the low-mass star regime, or if there are additional agents like merging of lower-mass forming objects or accretion from initially unbound material. Most of the uncertainties come from the lack of diagnostic tools to evolutionarily classify large samples of candidate massive protostellar objects that can then be studied in more detail. Aims. We want to verify whether diagnostic tools like the SED shape and the relationship between envelope mass and bolometric luminosity can be extended to the study of high-mass star formation. Methods. The 8−1200 µm SED of YSOs in 42 regions of massive star formation has been reconstructed using MSX, IRAS, and submm data partly available from previous works. Apart from IRAS catalogue fluxes, the fluxes in the Mid-IR and sub-mm/mm were derived directly from the images. The SEDs were fitted to an extensive grid of envelope models with embedded ZAMS stars, available from the literature. Sources that could not be fitted with a single model were then fitted with a two-component model composed of an embedded ZAMS for the mid-IR part and a single-temperature optically thin greybody for the longer wavelength emitting component. Sources were classified as "IR" if they were fitted with an embedded ZAMS envelope, and "MM" if they could only be fitted with a greybody with a peak at high λ; further subclassification was based on being the most massive object in the field ("P", for primary) or not ("S", for secondary). Results. The different classes of sources identified in our analysis have very different SEDs and occupy distinct areas in the L bol −M env diagram; by analogy with the low-mass regime, we see that MM-P, IR-P and IR-S objects could be interpreted as the high-mass analogue of Class 0-I-II. Evolutionary tracks obtained from a simple model based on the turbulent core prescriptions show that the three classes of sources possibly mark different periods in the formation of a massive YSO. The IR-P objects are consistent with being at the end of the main accretion phase, while MM-P sources are probably in an earlier evolutionary stage. The timescales for the formation decrease from ∼4 × 10 5 to ∼1 × 10 5 years with stellar mass increasing from ∼6 to ∼40 M ; these timescales, and the association with young clusters with median stellar age of a few 10 6 years suggest that the most massive objects are among the last ones to form. Conclusions. Our results are consistent with the high-mass star formation being a scaled-up analogue of the traditional accretiondominated paradigm valid for the low-mass regime.

IRAS 23385ϩ6053 is a young stellar object with a luminosity of ∼ L at a kinematic distance of 4

We present a study of molecular outflows toward a sample of 69 luminous IRAS point sources. The sample is associated with dense molecular gas and has far-infrared luminosities ranging from 10 2 to 10 5 L , indicating these objects as regions likely forming high-mass stars. Mapping in the CO J ¼ 2 1 line shows that molecular outflows are ubiquitous in these regions. Most of the outflows have masses of tens of M . The typical dynamical timescale of the flow, without correcting for inclination of the flow axis, is a few times 10 4 yr. The typical energy in the outflows is 10 46 ergs, comparable to the turbulent energy in the core. Nearly half of the outflows show spatially resolved bipolar lobes. This indicates that low-mass young stars that coexist in the region are not responsible for the bipolar outflows observed. It is the more massive stars that drive the outflow. The large detection rate of outflows in the region favors an accretion process in the formation of massive stars. The maximum mass-loss rate in the wind is about 10 À4 M yr À1 . If these outflows are driven via accretion, the accretion rate should be as high as a few times 10 À4 M yr À1

We present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key program 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 • × 2 • tiles approximately centered at l = 30 • and l = 59 • . 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 ∼ 1 is exceeded for the regions in the l = 59 • field; a A V value between 5 and 10 is found for the l = 30 • field, likely due to the relatively higher 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 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 δ < −30• , 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 M for a dust temperature T d = 30 K, an H 2 density of 9.5 × 10 5 cm −3 , and a surface density of 0.4 g cm −2 . The median values are 0.4 pc, 102 M , 4 × 10 4 cm −3 , and 0.14 g cm −2 , respectively. The mean value of the luminosity-to-mass ratio, L/M 99 L /M , 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 pointlike, 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 M is best fitted with a power-law dN/dM ∝ M −α with α = 2.1, consistent with the Salpeter (1955) stellar IMF, with α = 2.35. On the other hand, the mass function of clumps with masses 10 M < ∼ M < ∼ 120 M is better fitted with a power law of slope α = 1.5, more consistent with the mass function of molecular clouds derived from gas observations.