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.

Aims: The inner regions of high-mass protostars are often invisible in the near-infrared. We aim to investigate the inner gaseous disc of IRAS11101-5829 through scattered light from the outflow cavity walls. Methods: We observed the environment of the high-mass young stellar object IRAS11101-5829 and the closest knots of its jet, HH135-136, with the VLT/SINFONI. We also retrieved archival data from the high-resolution long-slit spectrograph VLT/X-shooter. Results: We detect the first three bandheads of the $\upsilon=2-0$ CO vibrational emission for the first time in this object. It is coincident with continuum and Br$\gamma$ emission and extends up to $\sim10000$ au to the north-east and $\sim10 000$ au to the south-west. The line profiles have been modelled as a Keplerian rotating disc assuming a single ring in LTE. The model output gives a temperature of $\sim3000$ K, a CO column density of $\sim1\times10^{22}\mathrm{ cm^{-2}}$, and a projected Keplerian velocity $v_\mathrm{K}\sin i_\mathrm{disc} \sim 25\mathrm{ km s^{-1}}$, which is consistent with previous modelling in other high-mass protostars. In particular, the low value of $v_\mathrm{K}\sin i_\mathrm{disc}$ suggests that the disc is observed almost face-on, whereas the well-constrained geometry of the jet imposes that the disc must be close to edge-on. This apparent discrepancy is interpreted as the CO seen reflected in the mirror of the outflow cavity wall. Conclusions: From both jet geometry and disc modelling, we conclude that all the CO emission is seen through reflection by the cavity walls and not directly. This result implies that in the case of highly embedded objects, as for many high-mass protostars, line profile modelling alone might be deceptive and the observed emission could affect the derived physical and geometrical properties; in particular the inclination of the system can be incorrectly interpreted.

Context. As a building block for amino acids, formamide (NH$_2$CHO) is an important molecule in astrobiology and astrochemistry, but its formation path in the interstellar medium is not understood well. Aims. We aim to find empirical evidence to support the chemical relationships of formamide to HNCO and H$_2$CO. Methods. We examine high angular resolution (~0.2") Atacama Large Millimeter/submillimeter Array (ALMA) maps of six sources in three high-mass star-forming regions and compare the spatial extent, integrated emission peak position, and velocity structure of HNCO and H$_2$CO line emission with that of NH$_2$CHO by using moment maps. Through spectral modeling, we compare the abundances of these three species. Results. In these sources, the emission peak separation and velocity dispersion of formamide emission is most often similar to HNCO emission, while the velocity structure is generally just as similar to H$_2$CO and HNCO (within errors). From the spectral modeling, we see that the abundances between all three of our focus species are correlated, and the relationship between NH$_2$CHO and HNCO reproduces the previously demonstrated abundance relationship. Conclusions. In this first interferometric study, which compares two potential parent species to NH$_2$CHO, we find that all moment maps for HNCO are more similar to NH$_2$CHO than H$_2$CO in one of our six sources (G24 A1). For the other five sources, the relationship between NH$_2$CHO, HNCO, and H$_2$CO is unclear as the different moment maps for each source are not consistently more similar to one species as opposed to the other.

To understand how Phosphorus-bearing molecules are formed in star-forming regions, we have analysed ALMA observations of PN and PO towards the massive star-forming region AFGL 5142, combined with a new analysis of the data of the comet 67P/Churyumov-Gerasimenko taken with the ROSINA instrument onboard Rosetta. The ALMA maps show that the emission of PN and PO arises from several spots associated with low-velocity gas with narrow linewidths in the cavity walls of a bipolar outflow. PO is more abundant than PN in most of the spots, with the PO/PN ratio increasing as a function of the distance to the protostar. Our data favor a formation scenario in which shocks sputter phosphorus from the surface of dust grains, and gas-phase photochemistry induced by UV photons from the protostar allows efficient formation of the two species in the cavity walls. Our analysis of the ROSINA data has revealed that PO is the main carrier of P in the comet, with PO/PN>10. Since comets may have delivered a significant amount of prebiotic material to the early Earth, this finding suggests that PO could contribute significantly to the phosphorus reservoir during the dawn of our planet. There is evidence that PO was already in the cometary ices prior to the birth of the Sun, so the chemical budget of the comet might be inherited from the natal environment of the Solar System, which is thought to be a stellar cluster including also massive stars.

We consider a dusty clump in the two cases of spherical and cylindrical symmetry to investigate the effect of temperature and density gradients on the observed flux density. Conversely, we evaluate how the presence of such gradients affects the calculation of the clump mass from the observed flux. We provide the reader with approximate expressions relating flux density and mass in the optically thick and thin limits, in the Rayleigh-Jeans regime, and discuss the reliability of these expressions by comparing them to the outcome of a numerical code. Finally, we present an application of our calculations to three examples taken from the literature, which shows how the correction introduced after taking into account temperature and density gradients may affect our conclusions on the stability of the clumps.

Context. Submillimeter Array (SMA) 870 micron polarization observations of the hot molecular core G31.41+0.31 revealed one of the clearest examples up to date of an hourglass-shaped magnetic field morphology in a high-mass star-forming region. Aims. To better establish the role that the magnetic field plays in the collapse of G31.41+0.31, we carried out Atacama Large Millimeter/submillimeter Array (ALMA) observations of the polarized dust continuum emission at 1.3 mm with an angular resolution four times higher than that of the previous (sub)millimeter observations to achieve an unprecedented image of the magnetic field morphology. Methods. We used ALMA to perform full polarization observations at 233 GHz (Band 6). The resulting synthesized beam is 0.28"x0"20 which, at the distance of the source, corresponds to a spatial resolution of ~875 au. Results. The observations resolve the structure of the magnetic field in G31.41+0.31 and allow us to study the field in detail. The polarized emission in the Main core of G31.41+0.41is successfully fit with a semi-analytical magnetostatic model of a toroid supported by magnetic fields. The best fit model suggests that the magnetic field is well represented by a poloidal field with a possible contribution of a toroidal component of ~10% of the poloidal component, oriented southeast to northwest at ~ -44 deg and with an inclination of ~-45 degr. The magnetic field is oriented perpendicular to the northeast to southwest velocity gradient detected in this core on scales from 1E3-1E4 au. This supports the hypothesis that the velocity gradient is due to rotation and suggests that such a rotation has little effect on the magnetic field. The strength of the magnetic field estimated in the central region of the core with the Davis-Chandrasekhar-Fermi method is ~8-13 mG and implies that the mass-to-flux ratio in this region is slightly supercritical ...

The formation process of high-mass stars (>8M$_\odot$) is poorly constrained, particularly, the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500M$_\odot$) star-forming region IRAS 23033+5951. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30-m single-dish telescope at 220GHz, we probe the gas and dust emission at an angular resolution of ~0.45arcsec, corresponding to 1900au. In the mm continuum emission, we identify a protostellar cluster with at least four mm-sources, where three of them show a significantly higher peak intensity well above a signal-to-noise ratio of 100. Hierarchical fragmentation from large to small spatial scales is discussed. Two fragments are embedded in rotating structures and drive molecular outflows, traced by $^{13}$CO (2-1) emission. The velocity profiles across two of the cores are similar to Keplerian but are missing the highest velocity components close to the center of rotation, which is a common phenomena from observations like these, and other rotation scenarios are not excluded entirely. Position-velocity diagrams suggest protostellar masses of ~6 and 19M$_\sun$. Rotational temperatures from fitting CH$_3$CN ($12_K-11_K$) spectra are used for estimating the gas temperature and by that the disk stability against gravitational fragmentation, utilizing Toomre's $Q$ parameter. [We] identify only one candidate disk to be unstable against gravitational instability caused by axisymmetric perturbations. The dominant sources cover different evolutionary stages within the same maternal gas clump. The appearance of rotation and outflows of the cores are similar to those found in low-mass star-forming regions.

We present the highest angular resolution (20x15mas - 44x33au) Atacama Large Millimeter/sub-millimeter Array (ALMA) observations currently possible of the proto-O-star G17.64+0.16 in Band 6. The Cycle 5 observations with baselines out to 16km probes scales <50au and reveal the rotating disc around G17.64+0.16, a massive forming O-type star. The disc has a ring-like enhancement in the dust emission, especially visible as arc structures to the north and south. The Keplerian kinematics are most prominently seen in the vibrationally excited water line, H2O (Eu=3461.9K). The mass of the central source found by modelling the Keplerian rotation is consistent with 45+/-10Mo. The H30alpha (231.9GHz) radio-recombination line and the SiO (5-4) molecular line were detected at up to the 10 sigma$ level. The estimated disc mass is 0.6-2.6Mo under the optically thin assumption. Analysis of the Toomre Q parameter, in the optically thin regime, indicates that the disc stability is highly dependent on temperature. The disc currently appears stable for temperatures >150K, this does not preclude that the substructures formed earlier through disc fragmentation.

This study is part of the project ``CORE'', an IRAM/NOEMA large program consisting of observations of the millimeter continuum and molecular line emission towards 20 selected high-mass star forming regions. We focus on IRAS23385+6053, which is believed to be the least evolved source of the CORE sample. The observations were performed at ~1.4 mm and employed three configurations of NOEMA and additional single-dish maps, merged with the interferometric data to recover the extended emission. Our correlator setup covered a number of lines from well-known hot core tracers and a few outflow tracers. The angular (~0.45"$-$0.9") and spectral (0.5 km/s) resolutions were sufficient to resolve the clump in IRAS23385+6053 and investigate the existence of large-scale motions due to rotation, infall, or expansion. We find that the clump splits into six distinct cores when observed at sub-arcsecond resolution. These are identified through their 1.4 mm continuum and molecular line emission. We produce maps of the velocity, line width, and rotational temperature from the methanol and methyl cyanide lines, which allow us to investigate the cores and reveal a velocity and temperature gradient in the most massive core. We also find evidence of a bipolar outflow, possibly powered by a low-mass star. We present the tentative detection of a circumstellar self-gravitating disk lying in the most massive core and powering a large-scale outflow previously known in the literature. In our scenario, the star powering the flow is responsible for most of the luminosity of IRAS23385+6053 (~$3000~L_\odot$). The other cores, albeit with masses below the corresponding virial masses, appear to be accreting material from their molecular surroundings and are possibly collapsing or on the verge of collapse. We conclude that we are observing a sample of star-forming cores that is bound to turn into a cluster of massive stars.

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.

It is well established that Solar-mass stars gain mass via disk accretion, until the mass reservoir of the disk is exhausted and dispersed, or condenses into planetesimals. Accretion disks are intimately coupled with mass ejection via polar cavities, in the form of jets and less collimated winds, which allow mass accretion through the disk by removing a substantial fraction of its angular momentum. Whether disk accretion is the mechanism leading to the formation of stars with much higher masses is still unclear. Here, we are able to build a comprehensive picture for the formation of an O-type star, by directly imaging a molecular disk which rotates and undergoes infall around the central star, and drives a molecular jet which arises from the inner disk regions. The accretion disk is truncated between 2000-3000au, it has a mass of about a tenth of the central star mass, and is infalling towards the central star at a high rate (6x10^-4 Msun/yr), as to build up a very massive object. These findings, obtained with the Atacama Large Millimeter/submillimeter Array at 700au resolution, provide observational proof that young massive stars can form via disk accretion much like Solar-mass stars.

We previously observed the star-forming region G16.59$-$0.05 through interferometric observations of both thermal and maser lines, and identified a high-mass young stellar object (YSO) which is surrounded by an accretion disk and drives a nonthermal radio jet. We performed high-angular-resolution (beam FWHM ~0.15") 1.2-mm continuum and line observations towards G16.59$-$0.05 with the Atacama Large Millimeter Array (ALMA). The main dust clump, with size ~10$^4$ au, is resolved into four relatively compact (diameter ~2000 au) millimeter (mm) sources. The source harboring the high-mass YSO is the most prominent in molecular emission. By fitting the emission profiles of several unblended and optically thin transitions of CH$_3$OCH$_3$ and CH$_3$OH, we derived gas temperatures inside the mm-sources in the range 42--131 K, and calculated masses of 1--5 $M_{\odot}$. A well-defined Local Standard of Rest velocity (Vlsr) gradient is detected in most of the high-density molecular tracers at the position of the high-mass YSO, pinpointed by compact 22-GHz free-free emission. This gradient is oriented along a direction forming a large (~70 degree) angle with the radio jet, traced by elongated 13-GHz continuum emission. The butterfly-like shapes of the P-V plots and the linear pattern of the emission peaks of the molecular lines at high velocity confirm that this Vlsr gradient is due to rotation of the gas in the disk surrounding the high-mass YSO. The disk radius is ~500 au, and the Vlsr distribution along the major axis of the disk is well reproduced by a Keplerian profile around a central mass of 10$\pm$2 $M_{\odot}$. The position of the YSO is offset by >~ 0.1" from the axis of the radio jet and the dust emission peak. To explain this displacement we argue that the high-mass YSO could have moved from the center of the parental mm source owing to dynamical interaction with one or more companions.

The formyl radical HCO has been proposed as the basic precursor of many complex organic molecules such as methanol (CH$_3$OH) or glycolaldehyde (CH$_2$OHCHO). Using ALMA, we have mapped, for the first time at high angular resolution ($\sim$1$^{\prime\prime}$, $\sim$140 au), HCO towards the Solar-type protostellar binary IRAS 16293$-$2422, where numerous complex organic molecules have been previously detected. We also detected several lines of the chemically related species H$_2$CO, CH$_3$OH and CH$_2$OHCHO. The observations revealed compact HCO emission arising from the two protostars. The line profiles also show redshifted absorption produced by foreground material of the circumbinary envelope that is infalling towards the protostars. Additionally, IRAM 30m single-dish data revealed a more extended HCO component arising from the common circumbinary envelope. The comparison between the observed molecular abundances and our chemical model suggests that whereas the extended HCO from the envelope can be formed via gas-phase reactions during the cold collapse of the natal core, the HCO in the hot corinos surrounding the protostars is predominantly formed by the hydrogenation of CO on the surface of dust grains and subsequent thermal desorption during the protostellar phase. The derived abundance of HCO in the dust grains is high enough to produce efficiently more complex species such as H$_2$CO, CH$_3$OH, and CH$_2$OHCHO by surface chemistry. We found that the main formation route of CH$_2$OHCHO is the reaction between HCO and CH$_2$OH.

We present high angular resolution 0.2 arcsec continuum and molecular emission line Atacama Large Millimeter/sub-millimeter Array (ALMA) observations of G17.64+0.16 in Band 6 (220GHz) taken as part of a campaign in search of circumstellar discs around (proto)-O-stars. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. At a resolution of 400au the main continuum core is essentially unresolved and isolated from other strong and compact emission peaks. We detect SiO (5-4) emission that is marginally resolved and elongated in a direction perpendicular to the large-scale outflow seen in the 13CO (2-1) line using the main ALMA array in conjunction with the Atacama Compact Array (ACA). Morphologically, the SiO appears to represent a disc-like structure. Using parametric models we show that the position-velocity profile of the SiO is consistent with the Keplerian rotation of a disc around an object between 10-30Mo in mass, only if there is also radial expansion from a separate structure. The radial motion component can be interpreted as a disc wind from the disc surface. Models with a central stellar object mass between 20 and 30Mo are the most consistent with the stellar luminosity (100000 Lo) and indicative of an O-type star. The H30a millimetre recombination line (231.9GHz) is also detected, but spatially unresolved, and is indicative of a very compact, hot, ionised region co-spatial with the dust continuum core. Accounting for all observables, we suggest that G17.64 is consistent with a O-type young stellar object in the final stages of protostellar assembly, driving a wind, but that has not yet developed into a compact HII region. The existance and detection of the disc in G17.64 is likely related to its isolated and possibly more evolved nature, traits which may underpin discs in similar sources.

Phosphorus is a crucial element in prebiotic chemistry, especially the P$-$O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO had only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detection of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches $\sim$35 K. The observed molecular abundances of 10$^{-10}$ are predicted by the model if a relatively high initial abundance of phosphorus, 5$\times$10$^{-9}$, is assumed.

The fragmentation mode of high-mass molecular clumps and the properties of the central rotating structures surrounding the most luminous objects have yet to be comprehensively characterised. Using the IRAM NOrthern Extended Millimeter Array (NOEMA) and the IRAM 30-m telescope, the CORE survey has obtained high-resolution observations of 20 well-known highly luminous star-forming regions in the 1.37 mm wavelength regime in both line and dust continuum emission. We present the spectral line setup of the CORE survey and a case study for W3(H2O). At ~0.35" (700 AU at 2 kpc) resolution, the W3(H2O) clump fragments into two cores (West and East), separated by ~2300 AU. Velocity shifts of a few km/s are observed in the dense-gas tracer, CH3CN, across both cores, consistent with rotation and perpendicular to the directions of two bipolar outflows, one emanating from each core. The kinematics of the rotating structure about W3(H2O) W shows signs of differential rotation of material, possibly in a disk-like object. The observed rotational signature around W3(H2O) E may be due to a disk-like object, an unresolved binary (or multiple) system, or a combination of both. We fit the emission of CH3CN (12-11) K = 4-6 and derive a gas temperature map with a median temperature of ~165 K across W3(H2O). We create a Toomre Q map to study the stability of the rotating structures against gravitational instability. The rotating structures appear to be Toomre unstable close to their outer boundaries, with a possibility of further fragmentation in the differentially-rotating core W3(H2O) W. Rapid cooling in the Toomre-unstable regions supports the fragmentation scenario. Combining millimeter dust continuum and spectral line data toward the famous high-mass star-forming region W3(H2O), we identify core fragmentation on large scales, and indications for possible disk fragmentation on smaller spatial scales.

As part of our effort to search for circumstellar disks around high-mass stellar objects, we observed the well-known core G31.41+0.31 with ALMA at 1.4 mm with an angular resolution of~0.22" (~1700 au). The dust continuum emission has been resolved into two cores namely Main and NE. The Main core, which has the stronger emission and is the more chemically rich, has a diameter of ~5300 au, and is associated with two free-free continuum sources. The Main core looks featureless and homogeneous in dust continuum emission and does not present any hint of fragmentation. Each transition of CH3CN and CH3OCHO, both ground and vibrationally excited, as well as those of CH3CN isotopologues, shows a clear velocity gradient along the NE-SW direction, with velocity linearly increasing with distance from the center, consistent with solid-body rotation. However, when comparing the velocity field of transitions with different upper level energies, the rotation velocity increases with increasing energy of the transition, which suggests that the rotation speeds up towards the center. Spectral lines towards the dust continuum peak show an inverse P-Cygni profile that supports the existence of infall in the core. The infall velocity increases with the energy of the transition suggesting that the infall is accelerating towards the center of the core, consistent with gravitational collapse. Despite the monolithic appearance of the Main core, the presence of red-shifted absorption, the existence of two embedded free-free sources at the center, and the rotational spin-up are consistent with an unstable core undergoing fragmentation with infall and differential rotation due to conservation of angular momentum. Therefore, the most likely explanation for the monolithic morphology is that the large opacity of the dust emission prevents the detection of any inhomogeneity in the core.

Aims: We aim to understand the fragmentation as well as the disk formation, outflow generation and chemical processes during high-mass star formation on spatial scales of individual cores. Methods: Using the IRAM Northern Extended Millimeter Array (NOEMA) in combination with the 30m telescope, we have observed in the IRAM large program CORE the 1.37mm continuum and spectral line emission at high angular resolution (~0.4'') for a sample of 20 well-known high-mass star-forming regions with distances below 5.5kpc and luminosities larger than 10^4Lsun. Results: We present the overall survey scope, the selected sample, the observational setup and the main goals of CORE. Scientifically, we concentrate on the mm continuum emission on scales on the order of 1000AU. We detect strong mm continuum emission from all regions, mostly due to the emission from cold dust. The fragmentation properties of the sample are diverse. We see extremes where some regions are dominated by a single high-mass core whereas others fragment into as many as 20 cores. A minimum-spanning-tree analysis finds fragmentation at scales on the order of the thermal Jeans length or smaller suggesting that turbulent fragmentation is less important than thermal gravitational fragmentation. The diversity of highly fragmented versus singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions: The smallest observed separations between cores are found around the angular resolution limit which indicates that further fragmentation likely takes place on even smaller spatial scales. The CORE project with its numerous spectral line detections will address a diverse set of important physical and chemical questions in the field of high-mass star formation.

Context. Recent observations of the massive young stellar object S255 NIRS 3 have revealed a large increase in both methanol maser flux density and IR emission, which have been interpreted as the result of an accretion outburst, possibly due to instabilities in a circumstellar disk. This indicates that this type of accretion event could be common in young/forming early-type stars and in their lower mass siblings, and supports the idea that accretion onto the star may occur in a non-continuous way. Aims. As accretion and ejection are believed to be tightly associated phenomena, we wanted to confirm the accretion interpretation of the outburst in S255 NIRS 3 by detecting the corresponding burst of the associated thermal jet. Methods. We monitored the radio continuum emission from S255 NIRS 3 at four bands using the Karl G. Jansky Very Large Array. The millimetre continuum emission was also observed with both the Northern Extended Millimeter Array of IRAM and the Atacama Large Millimeter/submillimeter Array. Results. We have detected an exponential increase in the radio flux density from 6 to 45 GHz starting right after July 10, 2016, namely about 13 months after the estimated onset of the IR outburst. This is the first ever detection of a radio burst associated with an IR accretion outburst from a young stellar object. The flux density at all observed centimetre bands can be reproduced with a simple expanding jet model. At millimetre wavelengths we infer a marginal flux increase with respect to the literature values and we show this is due to free-free emission from the radio jet. Abridged.

STATCONT is a python-based tool designed to determine the continuum emission level in spectral data, in particular for sources with a line-rich spectrum. The tool inspects the intensity distribution of a given spectrum and automatically determines the continuum level by using different statistical approaches. The different methods included in STATCONT are tested against synthetic data. We conclude that the sigma-clipping algorithm provides the most accurate continuum level determination, together with information on the uncertainty in its determination. This uncertainty can be used to correct the final continuum emission level, resulting in the here called `corrected sigma-clipping method' or c-SCM. The c-SCM has been tested against more than 750 different synthetic spectra reproducing typical conditions found towards astronomical sources. The continuum level is determined with a discrepancy of less than 1% in 50% of the cases, and less than 5% in 90% of the cases, provided at least 10% of the channels are line free. The main products of STATCONT are the continuum emission level, together with a conservative value of its uncertainty, and datacubes containing only spectral line emission, i.e., continuum-subtracted datacubes. STATCONT also includes the option to estimate the spectral index, when different files covering different frequency ranges are provided.

We present Submillimeter Array (SMA) observations toward the high-mass star-forming region IRAS 18566+0408. Observations at 1.3 mm continuum and in several molecular line transitions were performed in the compact (2."4 angular resolution) and very-extended (~0."4 angular resolution) configurations. The continuum emission from the compact configuration shows a dust core of 150 Msun, while the very-extended configuration reveals a dense (2.6 x 10^7 cm^-3) and compact (~4,000 AU) condensation of 8 Msun. We detect 31 molecular transitions from 14 species including CO isotopologues, SO, CH3OH, OCS, and CH3CN. Using the different k-ladders of the CH3CN line, we derive a rotational temperature at the location of the continuum peak of 240 K. The 12CO(2-1), 13CO(2-1), and SO(6_5-5_4) lines reveal a molecular outflow at PA ~135^o centered at the continuum peak. The extended 12CO(2-1) emission has been recovered with the IRAM 30 m telescope observations. Using the combined data set, we derive an outflow mass of 16.8 Msun. The chemically rich spectrum and the high rotational temperature confirm that IRAS 18566+0408 is harboring a hot molecular core. We find no clear velocity gradient that could suggest the presence of a rotational disk-like structure, even at the high resolution observations obtained with the very-extended configuration.

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.

In the study of high-mass star formation, hot cores are empirically defined stages where chemically rich emission is detected toward a massive YSO. It is unknown whether the physical origin of this emission is a disk, inner envelope, or outflow cavity wall and whether the hot core stage is common to all massive stars. We investigate the chemical make up of several hot molecular cores to determine physical and chemical structure. We use high spectral and spatial resolution Cycle 0 ALMA observations to determine how this stage fits into the formation sequence of a high mass star. We observed the G35.20-0.74N and G35.03+0.35 hot cores at 350 GHz. We analyzed spectra and maps from four continuum peaks (A, B1, B2 and B3) in G35.20, separated by 1000-2000 AU, and one continuum peak in G35.03. We made all possible line identifications across 8 GHz of spectral windows of molecular emission lines and determined column densities and temperatures for as many as 35 species assuming local thermodynamic equilibrium. In comparing the spectra of the four peaks, we find each has a distinct chemical composition expressed in over 400 different transitions. In G35.20, B1 and B2 contain oxygen- and sulfur-bearing organic and inorganic species but few nitrogen-bearing species whereas A and B3 are strong sources of O, S, and N-bearing species (especially those with the CN-bond). CH$_2$DCN is clearly detected in A and B3 with D/H ratios of 8 and 13$\%$, respectively, but is much weaker at B1 and undetected at B2. No deuterated species are detected in G35.03, but similar molecular abundances to G35.20 were found in other species. We also find co-spatial emission of HNCO and NH$_2$CHO in both sources indicating a strong chemical link between the two species. The chemical segregation between N-bearing organic species and others in G35.20 suggests the presence of multiple protostars, surrounded by a disk or torus.

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.

Solar-mass stars form via circumstellar disk accretion (disk-mediated accretion). Recent findings indicate that this process is likely episodic in the form of accretion bursts, possibly caused by disk fragmentation. Although it cannot be ruled out that high-mass young stellar objects (HMYSOs; $M>$8 M$_\odot$, $L_{bol}>$5$\times$10$^3$ L$_\odot$) arise from the coalescence of their low-mass brethren, latest results suggest that they more likely form via disks. Accordingly, disk-mediated accretion bursts should occur. Here we report on the discovery of the first disk-mediated accretion burst from a $\sim$20 M$_\odot$ HMYSO. Our near-infrared images show the brightening of the central source and its outflow cavities. Near-infrared spectroscopy reveals emission lines typical of accretion bursts in low-mass protostars, but orders of magnitude more luminous. Moreover, the energy released and the inferred mass-accretion rate are also orders of magnitude larger. Our results identify disk accretion as the common mechanism of star formation across the entire stellar mass spectrum.

The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission in line-rich sources, we use a new statistical method: STATCONT. We detect 27 continuum sources in SgrB2(M) and 20 in SgrB2(N). We study the continuum emission across the ALMA band 6, and compare it with previous SMA 345 GHz and VLA 40 GHz observations, to study the nature of the sources detected. The brightest sources are dominated by (partially optically thick) dust emission, while there is an important degree of contamination from ionized gas free-free emission in weaker sources. While the total mass in SgrB2(M) is distributed in many fragments, most of the mass in SgrB2(N) arises from a single object, with filamentary-like structures converging towards the center. There seems to be a lack of low-mass dense cores in both regions. We determine H2 volume densities for the cores of about 10^5-10^7 Msun pc^-3, one to two orders of magnitude higher than the stellar densities of super star clusters. In general, SgrB2(N) is chemically richer than SgrB2(M). There seems to be a correlation between the chemical richness and the mass of the fragments, with more massive clumps being more chemically rich. Both SgrB2(N) and SgrB2(M) harbour a cluster of hot molecular cores. We compare the continuum images with predictions from a detailed 3D radiative transfer model that reproduces the structure of SgrB2 from 45 pc down to 100 au. This dataset, together with ongoing projects in the range 5 to 200 GHz, better constrain the 3D structure of SgrB2, and allow us to understand its physical and chemical structure.

We report on subarcsecond observations of complex organic molecules (COMs) in the high-mass protostar IRAS20126+4104 with the Plateau de Bure Interferometer in its most extended configurations. In addition to the simple molecules SO, HNCO and H2-13CO, we detect emission from CH3CN, CH3OH, HCOOH, HCOOCH3, CH3OCH3, CH3CH2CN, CH3COCH3, NH2CN, and (CH2OH)2. SO and HNCO present a X-shaped morphology consistent with tracing the outflow cavity walls. Most of the COMs have their peak emission at the putative position of the protostar, but also show an extension towards the south(east), coinciding with an H2 knot from the jet at about 800-1000 au from the protostar. This is especially clear in the case of H2-13CO and CH3OCH3. We fitted the spectra at representative positions for the disc and the outflow, and found that the abundances of most COMs are comparable at both positions, suggesting that COMs are enhanced in shocks as a result of the passage of the outflow. By coupling a parametric shock model to a large gas-grain chemical network including COMs, we find that the observed COMs should survive in the gas phase for about 2000 yr, comparable to the shock lifetime estimated from the water masers at the outflow position. Overall, our data indicate that COMs in IRAS20126+4104 may arise not only from the disc, but also from dense and hot regions associated with the outflow.

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 detection of organic molecules with increasing complexity and potential biological relevance is opening the possibility to understand the formation of the building blocks of life in the interstellar medium. One of the families of molecules with astrobiological interest are the esters, whose simplest member, methyl formate, is rather abundant in star-forming regions. The next step in the chemical complexity of esters is ethyl formate, C$_2$H$_5$OCHO. Only two detections of this species have been reported so far, which strongly limits our understanding of how complex molecules are formed in the interstellar medium. We have searched for ethyl formate towards the W51 e2 hot molecular core, one of the most chemically rich sources in the Galaxy and one of the most promising regions to study prebiotic chemistry, especially after the recent discovery of the P$-$O bond, key in the formation of DNA. We have analyzed a spectral line survey towards the W51 e2 hot molecular core, which covers 44 GHz in the 1, 2 and 3 mm bands, carried out with the IRAM 30m telescope. We report the detection of the trans and gauche conformers of ethyl formate. A Local Thermodynamic Equilibrium analysis indicates that the excitation temperature is 78$\pm$10 K and that the two conformers have similar source-averaged column densities of (2.0$\pm$0.3)$\times$10$^{16}$ cm$^{-2}$ and an abundance of $\sim$10$^{-8}$. We compare the observed molecular abundances of ethyl formate with different competing chemical models based on grain surface and gas-phase chemistry. We propose that grain-surface chemistry may have a dominant role in the formation of ethyl formate (and other complex organic molecules) in hot molecular cores, rather than reactions in the gas phase.

Aims: We want to study the physical properties of the ionized jet emission in the vicinity of an O-type young stellar object (YSO), and estimate how efficient is the transfer of energy and momentum from small- to large-scale outflows. Methods: We conducted Karl G. Jansky Very Large Array (VLA) observations, at both 22 and 45 GHz, of the compact and faint radio continuum emission in the high-mass star-forming region G023.01-00.41, with an angular resolution between 0.3" and 0.1", and a thermal rms of the order of 10 uJy/beam. Results: We discovered a collimated thermal (bremsstrahlung) jet emission, with a radio luminosity (L_rad) of 24 mJy kpc^2 at 45 GHz, in the inner 1000 AU from an O-type YSO. The radio thermal jet has an opening angle of 44 degrees and brings a momentum rate of 8 10^-3 M_sun yr^-1 km/s. By combining the new data with previous observations of the molecular outflow and water maser shocks, we can trace the outflow emission from its driving source through the molecular clump, across more than two order of magnitude in length (500 AU-0.2 pc). We find that the momentum-transfer efficiency, between the inner jet emission and the extended outflow of entrained ambient gas, is near unity. This result suggests that the large-scale flow is swept-up by the mechanical force of the radio jet emission, which originates in the inner 1000 AU from the high-mass YSO.

We present a high sensitivity radio continuum survey at 6 and 1.3$\,$cm using the Karl G. Jansky Very Large Array towards a sample of 58 high-mass star forming regions. Our sample was chosen from dust clumps within infrared dark clouds with and without IR sources (CMC-IRs, CMCs, respectively), and hot molecular cores (HMCs), with no previous, or relatively weak radio continuum detection at the $1\,$mJy level. Due to the improvement in the continuum sensitivity of the VLA, this survey achieved map rms levels of $\sim$ 3-10 $\mu$Jy beam$^{-1}$ at sub-arcsecond angular resolution. We extracted 70 centimeter continuum sources associated with 1.2$\,$mm dust clumps. Most sources are weak, compact, and are prime candidates for high-mass protostars. Detection rates of radio sources associated with the mm dust clumps for CMCs, CMC-IRs and HMCs are 6$\%$, 53$\%$ and 100$\%$, respectively. This result is consistent with increasing high-mass star formation activity from CMCs to HMCs. The radio sources located within HMCs and CMC-IRs occur close to the dust clump centers with a median offset from it of 12,000$\,$AU and 4,000$\,$AU, respectively. We calculated 5 - 25$\,$GHz spectral indices using power law fits and obtain a median value of 0.5 (i.e., flux increasing with frequency), suggestive of thermal emission from ionized jets. In this paper we describe the sample, observations, and detections. The analysis and discussion will be presented in Paper II.

Massive stars, multiple stellar systems and clusters are born from the gravitational collapse of massive dense gaseous clumps, and the way these systems form strongly depends on how the parent clump fragments into cores during collapse. Numerical simulations show that magnetic fields may be the key ingredient in regulating fragmentation. Here we present ALMA observations at ~0.25'' resolution of the thermal dust continuum emission at ~278 GHz towards a turbulent, dense, and massive clump, IRAS 16061-5048c1, in a very early evolutionary stage. The ALMA image shows that the clump has fragmented into many cores along a filamentary structure. We find that the number, the total mass and the spatial distribution of the fragments are consistent with fragmentation dominated by a strong magnetic field. Our observations support the theoretical prediction that the magnetic field plays a dominant role in the fragmentation process of massive turbulent clump.

We study the molecular abundance and spatial distribution of the simplest sugar alcohol, ethylene glycol (EG), the simplest sugar glycoladehyde (GA), and other chemically related complex organic species towards the massive star-forming region G31.41+0.31. We have analyzed multiple single dish and interferometric data, and obtained excitation temperatures and column densities using an LTE analysis. We have reported for the first time the presence of EG towards G31.41+0.31, and we have also detected multiple transitions of other complex organic molecules such as GA, methyl formate (MF), dimethyl ether (DME) and ethanol (ET). The high angular resolution images show that the EG emission is very compact, peaking towards the maximum of the continuum. These observations suggest that low abundance complex organic molecules, like EG or GA, are good probes of the gas located closer to the forming stars. Our analysis confirms that EG is more abundant than GA in G31.41+0.31, as previously observed in other interstellar regions. Comparing different star-forming regions we find evidence of an increase of the EG/GA abundance ratio with the luminosity of the source. The DME/MF and EG/ET ratios are nearly constant with luminosity. We have also found that the abundance ratios of pairs of isomers GA/MF and ET/DME decrease with the luminosity of the sources. The most likely explanation for the behavior of the EG/GA ratio is that these molecules are formed by different chemical formation routes not directly linked; although warm-up timescales effects and different formation and destruction efficiencies in the gas phase cannot be ruled out. The most likely formation route of EG is by combination of two CH$_{2}$OH radicals on dust grains. We also favor that GA is formed via the solid-phase dimerization of the formyl radical HCO, and a chemical link between MF and DME.

Context. ALMA observations of the high-mass star-forming region G35.20-0.74N have revealed the presence of a Keplerian disk in core B rotating about a massive object of 18 Msun, as computed from the velocity field. The luminosity of such a massive star would be comparable to (or higher than) the luminosity of the whole star-forming region. To solve this problem it has been proposed that core B could harbor a binary system. This could also explain the possible precession of the radio jet associated with this core, which has been suggested by its S-shaped morphology. Aims. To establish the origin of the free-free emission from core B and investigate the existence of a binary system at the center of this massive core and the possible precession of the radio jet. Methods. We carried out VLA continuum observations of G35.20-0.74N at 2 cm in the B configuration and at 1.3 cm and 7 mm in the A and B configurations. The bandwidth at 7 mm covers the CH3OH maser line at 44.069 GHz. Continuum images at 6 and 3.6 cm in the A configuration were obtained from the VLA archive. We also carried out VERA observations of the H2O maser line at 22.235 GHz. Results. The observations have revealed the presence of a binary system of UC/HC HII regions at the geometrical center of the radio jet in G35.20-0.74N. This binary system, which is associated with a Keplerian rotating disk, consists of two B-type stars of 11 and 6 Msun. The S-shaped morphology of the radio jet has been successfully explained as due to precession produced by the binary system. The analysis of the precession of the radio jet has allowed us to better interpret the IR emission in the region, which would be not tracing a wide-angle cavity open by a single outflow with a position angle of ~55 deg but two different flows: a precessing one in the NE-SW direction associated with the radio jet, and a second one in an almost E--W direction ...

Phosphorus is a crucial element in biochemistry, especially the P-O bond, which is key for the formation of the backbone of the deoxyribonucleic acid. So far, PO has only been detected towards the envelope of evolved stars, and never towards star-forming regions. We report the first detections of PO towards two massive star-forming regions, W51 e1/e2 and W3(OH), using data from the IRAM 30m telescope. PN has also been detected towards the two regions. The abundance ratio PO/PN is 1.8 and 3 for W51 and W3(OH), respectively. Our chemical model indicates that the two molecules are chemically related and are formed via gas-phase ion-molecule and neutral-neutral reactions during the cold collapse. The molecules freeze out onto grains at the end of the collapse and desorb during the warm-up phase once the temperature reaches 35 K. Similar abundances of the two species are expected during a period of 5x10^{4} yr at the early stages of the warm-up phase, when the temperature is in the range 35-90 K. The observed molecular abundances of 10^{-10} are predicted by the model if a relatively high initial abundance of 5x10^{-9} of depleted phosphorus is assumed.

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.

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.

We want to study the velocity and magnetic field morphology in the vicinity (<1000 AU) of a massive young stellar object (YSO), at very high spatial resolution (10-100 AU). We performed milli-arcsecond polarimetric observations of the strong CH3OH maser emission observed in the vicinity of an O-type YSO, in G023.01-00.41. We have combined this information with the velocity field of the CH3OH masing gas previously measured at the same angular resolution. We analyse the velocity and magnetic fields in the reference system defined by the direction of the molecular outflow and the equatorial plane of the hot molecular core at its base, as recently observed on sub-arcsecond scales. We provide a first detailed picture of the gas dynamics and magnetic field configuration within a radius of 2000 AU from a massive YSO. We have been able to reproduce the magnetic field lines for the outer regions (>600 AU) of the molecular envelope, where the magnetic field orientation shows a smooth change with the maser cloudlets position (0.2 degree/AU). Overall, the velocity field vectors well accommodate with the local, magnetic field direction, but still show an average misalignment of 30 degrees. We interpret this finding as the contribution of a turbulent velocity field of about 3.5 km/s, responsible for braking up the alignment between the velocity and magnetic field vectors. We do resolve different gas flows which develop both along the outflow axis and across the disk plane, with an average speed of 7 km/s. In the direction of the outflow axis, we establish a collimation of the gas flow, at a distance of about 1000 AU from the disk plane. In the disk region, gas appears to stream outward along the disk plane for radii greater than 500-600 AU, and inward for shorter radii.

The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.

The high-mass star-forming site G009.62+00.20E hosts the 6.7 GHz methanol maser source with the greatest flux density in the Galaxy which has been flaring periodically over the last ten years. We performed high-resolution astrometric measurements of the CH3OH, H2O, and OH maser emission and 7 mm continuum in the region. The radio continuum emission was resolved in two sources separated by 1300 AU. The CH3OH maser cloudlets are distributed along two north-south ridges of emission to the east and west of the strongest radio continuum component. This component likely pinpoints a massive young stellar object which heats up its dusty envelope, providing a constant IR pumping for the Class II CH3OH maser transitions. We suggest that the periodic maser activity may be accounted for by an independent, pulsating, IR radiation field provided by a bloated protostar in the vicinity of the brightest masers. We also report about the discovery of an elliptical distribution of CH3OH maser emission in the region of periodic variability.

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 report 1.3 cm and 6 cm continuum observations toward the massive proto-stellar candidate G11.11-0.12P1 using the Karl G. Jansky Very Large Array (VLA). We detect a string of four unresolved radio continuum sources coincident with the mid-IR source in G11P1. The continuum sources have positive spectral indices consistent with a thermal (free-free) ionized jet. The most likely origin of the ionized gas are shocks due to the interaction of a stellar wind with the surrounding high-density material. We also present NIR United Kingdom Infrared Telescope (UKIRT) archival data which show an extended structure detected only at K-band (2.2 $\mu$m), which is oriented perpendicular to the jet, and that may be scattered light from a circumstellar disk around the massive protostar. Our observations plus the UKIRT archival data thus provide new evidence that a disk/jet system is present in the massive protostellar candidate located in the G11.11-0.12P1 core.

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.

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

We wish to investigate the structure, velocity field, and stellar content of the G16.59-0.05 high-mass star-forming region, where previous studies have established the presence of two almost perpendicular (NE-SW and SE-NW), massive outflows, and a rotating disk traced by methanol maser emission. We performed Very Large Array observations of the radio continuum and ammonia line emission, complemented by COMICS/Subaru and Hi-GAL/Herschel images in the mid- and far-infrared (IR). Our centimeter continuum maps reveal a collimated radio jet that is oriented E-W and centered on the methanol maser disk, placed at the SE border of a compact molecular core. The spectral index of the jet is negative, indicating non-thermal emission over most of the jet, except the peak close to the maser disk, where thermal free-free emission is observed. We find that the ammonia emission presents a bipolar structure consistent (on a smaller scale) in direction and velocity with that of the NE-SW bipolar outflow detected in previous CO observations. After analyzing our previous N2H+(1-0) observations again, we conclude that two scenarios are possible. In one case both the radio jet and the ammonia emission would trace the root of the large-scale CO bipolar outflow. The different orientation of the jet and the ammonia flow could be explained by precession and/or a non-isotropic density distribution around the star. In the other case, the N2H+(1-0) and ammonia bipolarity is interpreted as two overlapping clumps moving with different velocities along the line of sight. The ammonia gas also seems to undergo rotation consistent with the maser disk. Our IR images complemented by archival data allow us to derive a bolometric luminosity of about 10^4 L_sun and to conclude that most of the luminosity is due to the young stellar object associated with the maser disk.

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

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 report on ALMA observations of continuum and molecular line emission with 0.4" resolution towards the high-mass star forming region G35.20-0.74 N. Two dense cores are detected in typical hot-core tracers, such as CH3CN, which reveal velocity gradients. In one of these cores, the velocity field can be fitted with an almost edge-on Keplerian disk rotating about a central mass of 18 Msun. This finding is consistent with the results of a recent study of the CO first overtone bandhead emission at 2.3mum towards G35.20-0.74 N. The disk radius and mass are >2500 au and 3 Msun. To reconcile the observed bolometric luminosity (3x10^4 Lsun) with the estimated stellar mass of 18 Msun, we propose that the latter is the total mass of a binary system.

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 imaged the molecular outflows towards the cluster of high-mass young stellar objects G24.78+0.08 at high-angular resolution using SiO emission, which is considered the classical tracer of protostellar jets. We performed SiO observations with the VLA interferometer in the J = 1-0 v=0 transition and with the SMA array in the 5-4 transition. A complementary IRAM 30-m single-dish survey in the (2-1), (3-2), (5-4), and (6-5) SiO lines was also carried out. Two collimated SiO high-velocity outflows driven by the A2 and C millimeter continuum massive cores have been imaged. On the other hand, we detected no SiO outflow driven by the young stellar objects in more evolved evolutionary phases that are associated with ultracompact (B) or hypercompact (A1) HII regions. The LVG analysis reveals high-density gas (10^3-10^4 cm-3), with well constrained SiO column densities (0.5-1 10^15 cm-2). The driving source of the A2 outflow is associated with typical hot core tracers such as methyl formate, vinyl cyanide, cyanoacetilene, and acetone. The driving source of the main SiO outflow in G24 has an estimated luminosity of a few 10^4 Lsun (typical of a late O-type star) and is embedded in the 1.3 mm continuum core A2, which in turn is located at the centre of a hot core that rotates on a plane perpendicular to the outflow main axis. The present SiO images support a scenario similar to the low-mass case for massive star formation, where jets that are clearly traced by SiO emission, create outflows of swept-up ambient gas usually traced by CO.

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed with the PdBI in the most extended configuration the continuum at 1.3 mm and the CO(2-1) emission of four massive cores. We detect dust condensations down to ~0.3 Msun and separate millimeter sources down to 0.4" or ~1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO(2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning 3 orders of magnitude. Among the 18 regions, ~30% show no signs of fragmentation, while 50% split up into ~4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of magnetic field, radiative feedback and turbulence in the fragmentation process, we compared our observations to radiation magneto-hydrodynamic simulations, and found that the low-fragmented regions are well reproduced in the magnetized core case, while the highly-fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

This work presents a detailed study of the gas kinematics towards the "Hot Molecular Core" (HMC) G31.41+0.31 via multi-epoch VLBI observations of the H2O 22 GHz and CH3OH 6.7 GHz masers, and single-epoch VLBI of the OH 1.6 GHz masers. Water masers present a symmetric spatial distribution with respect to the HMC center, where two nearby (0.2" apart), compact, VLA sources (labeled "A" and "B") are previously detected. The spatial distribution of a first group of water masers, named "J1", is well fit with an elliptical profile, and the maser proper motions mainly diverge from the ellipse center, with average speed of 36 km s-1. These findings strongly suggest that the "J1" water maser group traces the heads of a young (dynamical time of 1.3 10^3 yr), powerful (momentum rate of ~0.2 M_sun yr-1 km s-1), collimated (semi-opening angle ~10 deg) jet emerging from a MYSO located close (within 0.15") to the VLA source "B". Most of the water features not belonging to "J1" present an elongated (about 2" in size), NE--SW oriented (PA = 70 deg), S-shape distribution, which we denote with the label "J2". The elongated distribution of the "J2" group and the direction of motion, approximately parallel to the direction of elongation, of most "J2" water masers suggests the presence of another collimated outflow, emitted from a MYSO near the VLA source "A". The orientation of the "J2" jet agrees well with that (PA = 68 deg) of the well-defined V_LSR gradient across the HMC revealed by previous interferometric, thermal line observations. Furthermore, the "J2" jet is powerful enough to sustain the large momentum rate, 0.3 M_sun yr-1 km s-1, estimated assuming that the V_LSR gradient represents a collimated outflow. These two facts lead us to favour the interpretation of the V_LSR gradient across the G31.41+0.31 HMC in terms of a compact and collimated outflow.

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

Infrared-dark high-mass clumps are among the most promising objects to study the initial conditions of the formation process of high-mass stars and rich stellar clusters. In this work, we have observed the (3-2) rotational transition of C18O with the APEX telescope, and the (1,1) and (2,2) inversion transitions of NH3 with the Australia Telescope Compact Array in 21 infrared-dark clouds already mapped in the 1.2 mm continuum, with the aim of measuring basic chemical and physical parameters such as the CO depletion factor (fD), the gas kinetic temperature and the gas mass. In particular, the C18O (3-2) line allows us to derive fD in gas at densities higher than that traced by the (1-0) and (2-1) lines, typically used in previous works. We have detected NH3 and C18O in all targets. The clumps possess mass, H2 column and surface densities consistent with being potentially the birthplace of high-mass stars. We have measured fD in between 5 and 78, with a mean value of 32 and a median of 29. These values are, to our knowledge, larger than the typical CO depletion factors measured towards infrared-dark clouds and high-mass dense cores, and are comparable to or larger than the values measured in low-mass pre-stellar cores close to the onset of the gravitational collapse. This result suggests that the earliest phases of the high-mass star and stellar cluster formation process are characterised by fD larger than in low-mass pre-stellar cores. Thirteen out of 21 clumps are undetected in the 24 {\mu}m Spitzer images, and have slightly lower kinetic temperatures, masses and H2 column densities with respect to the eight Spitzer-bright sources. This could indicate that the Spitzer-dark clumps are either less evolved or are going to form less massive objects.

We used the Very Long Baseline Array (VLBA) and the European VLBI Network (EVN) to perform phase-referenced VLBI observations of the three most powerful maser transitions associated with the high-mass star-forming region G28.87+0.07: the 22.2 GHz H$_{2}$O, 6.7 GHz CH$_{3}$OH, and 1.665 GHz OH lines. We also performed VLA observations of the radio continuum emission at 1.3 and 3.6 cm and Subaru observations of the continuum emission at 24.5 $\mu$m. Two centimeter continuum sources are detected and one of them (named "HMC") is compact and placed at the center of the observed distribution of H$_{2}$O, CH$_{3}$OH and OH masers. The bipolar distribution of line-of-sight (l.o.s) velocities and the pattern of the proper motions suggest that the water masers are driven by a (proto)stellar jet interacting with the dense circumstellar gas. The same jet could both excite the centimeter continuum source named "HMC" (interpreted as free-free emission from shocked gas) and power the molecular outflow observed at larger scales -- although one cannot exclude that the free-free continuum is rather originating from a hypercompact \ion{H}{2} region. At 24.5 $\mu$m, we identify two objects separated along the north-south direction, whose absolute positions agree with those of the two VLA continuum sources. We establish that $\sim$90% of the luminosity of the region ($\sim$\times10^{5} L_\sun$) is coming from the radio source "HMC", which confirms the existence of an embedded massive young stellar object (MYSO) exciting the masers and possibly still undergoing heavy accretion from the surrounding envelope.

We made a detailed observational analysis of a well known hot molecular core lying in the high-mass star-forming region G31.41+0.31. This core is believed to contain deeply embedded massive stars and presents a velocity gradient that has been interpreted either as rotation or as expansion, depending on the authors. Our aim was to shed light on this question and possibly prepare the ground for higher resolution ALMA observations which could directly detect circumstellar disks around the embedded massive stars. Observations at sub-arcsecond resolution were performed with the Submillimeter Array in methyl cyanide, a typical hot molecular core tracer, and 12CO and 13CO, well known outflow tracers. We also obtained sensitive continuum maps at 1.3 mm. Our findings confirm the existence of a sharp velocity gradient across the core, but cannot confirm the existence of a bipolar outflow perpendicular to it. The improved angular resolution and sampling of the uv plane allow us to attain higher quality channel maps of the CH3CN lines with respect to previous studies and thus significantly improve our knowledge of the structure and kinematics of the hot molecular core. While no conclusive argument can rule out any of the two interpretations (rotation or expansion) proposed to explain the velocity gradient observed in the core, in our opinion the observational evidence collected so far indicates the rotating toroid as the most likely scenario. The outflow hypothesis appears less plausible, because the dynamical time scale is too short compared to that needed to form species such as CH3CN, and the mass loss and momentum rates estimated from our measurements appear too high.

Context. This study is part of a large project to study the physics of accretion and molecular outflows towards a selected sample of high-mass star-forming regions that show evidence of infall and rotation from previous studies. Aims. We wish to make a thorough study at high-angular resolution of the structure and kinematics of the HMCs and corresponding molecular outflows in the high-mass star-forming region G24.78+0.08. Methods. We carried out SMA and IRAM PdBI observations at 1.3 and 1.4 mm, respectively, of dust and of typical high-density and molecular outflow tracers with resolutions of <1". Complementary IRAM 30-m 12CO and 13CO observations were carried out to recover the short spacing information of the molecular outflows. Results. The millimeter continuum emission towards cores G24 A1 and A2 has been resolved into 3 and 2 cores, respectively, and named A1, A1b, A1c, A2, and A2b. All these cores are aligned in a southeast-northwest direction coincident with that of the molecular outflows detected in the region, which suggests a preferential direction for star formation in this region. The masses of the cores range from 7 to 22 Msun, and the rotational temperatures from 128 to 180 K. The high-density tracers have revealed the existence of 2 velocity components towards A1, one of them peaks close to the position of the millimeter continuum peak and of the HC HII region, and is associated with the velocity gradient seen in CH3CN towards this core, while the other one peaks southwest of core A1 and is not associated with any millimeter continuum emission peak. The position-velocity plots along outflow A and the 13CO averaged blueshifted and redshifted emission indicate that this outflow is driven by core A2.

Direct observations of accretion disks around high-mass young stellar objects would help to discriminate between different models of formation of massive stars. However, given the complexity of massive star forming regions, such studies are still limited in number. Additionally, there is still no general consensus on the molecular tracers to be used for such investigations. Because of its close distance and high luminosity, IRAS 17233-3606 is a potential good laboratory to search for traces of rotation in the inner gas around the protostar(s). Therefore, we selected the source for a detailed analysis of its molecular emission at 230 GHz with the SMA. We systematically investigated the velocity fields of transitions in the SMA spectra which are not affected by overlap with other transitions, and searched for coherent velocity gradients to compare them to the distribution of outflows in the region. Beside CO emission we also used high-angular H2 images to trace the outflow motions driven by the IRAS 17233-3606 cluster. We find linear velocity gradients in many transitions of the same molecular species and in several molecules. We report the first detection of HNCO in molecular outflows from massive YSOs. We discuss the CH3CN velocity gradient taking into account various scenarios: rotation, presence of multiple unresolved sources with different velocities, and outflow(s). Although other interpretations cannot be ruled out, we propose that the CH3CN emission might be affected by the outflows of the region. Higher angular observations are needed to discriminate between the different scenarios. The present observations, with the possible association of CH3CN with outflows in a few thousands AU around the YSOs' cluster, (i) question the choice of the tracer to probe rotating structures, and (ii) show the importance of the use of H2 images for detailed studies of kinematics.

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.

IRAS 20126+4104 is one of the best candidates for a high-mass (proto)star surrounded by an accretion disk. Such a fact may be used to set constraints on theories of high-mass star formation, but requires confirmation that the mass and luminosity of IRAS 20126+4104 are indeed typical of a B0.5 star, which in turn requires an accurate estimate of the distance. We used the Very Long Baseline Array and the European VLBI Network to observe the 22.2 GHz water and 6.7 GHz methanol masers in IRAS 20126+4104 at a number of epochs suitably distributed in time. The absolute positions of the maser features were established with respect to reference quasars, which allowed us to derive absolute proper motions. From the parallax of the water masers we obtain a distance of 1.64 \pm 0.05 kpc, which is very similar to the value adopted so far in the literature (1.7 kpc) and confirms that IRAS 20126+4104 is a high-mass (proto)star. From the methanol masers we derive the component in the plane of the sky of the systemic velocity of the disk+star system (-16 km/s in right-ascension and +7.6 km/s in declination). Accurate knowledge of the distance and systemic velocity allows us to improve on the model fit to the water maser jet presented in a previous study. Finally, we identify two groups of methanol maser features, one undergoing rotation in the disk and possibly distributed along a narrow ring centered on the star, the other characterised by relative proper motions indicating that the features are moving away from the disk, perpendicular to it. We speculate that the latter group might be tracing the disk material marginally entrained by the jet.

Context. In recent years, we have detected clear evidence of rotation in more than 5 hot molecular cores (HMCs). Their identification is confirmed by the fact that the rotation axes are parallel to the axes of the associated bipolar outflows. We have now pursued our investigation by extending the sample to 3 known massive cores, G10.62-0.38, G19.61-0.23, and G29.96-0.02. Aims. We wish to make a thorough study of the structure and kinematics of HMCs and corresponding molecular outflows to reveal possible velocity gradients indicative of rotation of the cores. Methods. We carried out PdBI observations at 2.7 and 1.4~mm of gas and dust with angular resolutions of 2"-3", and 1"-2", respectively. To trace both rotation and expansion, we simultaneously observed CH3CN, a typical HMC tracer, and 13CO, a typical outflow tracer. Results. The CH3CN(12-11) observations have revealed the existence of clear velocity gradients in the three HMCs oriented perpendicular to the direction of the bipolar outflows. For G19 and G29 the molecular outflows have been mapped in 13CO. The gradients have been interpreted as rotating toroids. The rotation temperatures, used to derive the mass of the cores, have been obtained by means of the rotational diagram method, and lie in the range of 87-244 K. The diameters and masses of the toroids lie in the range of 4550-12600 AU, and 28-415 Msun, respectively. Given that the dynamical masses are 2 to 30 times smaller than the masses of the cores (if the inclination of the toroids with respect to the plane of the sky is not much smaller than 45 degrees), we suggest that the toroids could be accreting onto the embedded cluster. For G19 and G29, the collapse is also suggested by the redshifted absorption seen in the 13CO(2-1) line. We infer that infall onto the embedded (proto)stars must proceed with rates of 1E-2 Msun/yr, and on timescales of the order of 4E3-1E4yr...

Aims: The main goal of this study is to perform a sub-arcsecond resolution analysis of the high-mass star formation region G19.61-0.23, both in the continuum and molecular line emission. While the centimeter continuum images will be discussed in detail in a forthcoming paper, here we focus on the (sub)mm emission, devoting special attention to the hot molecular core. Results: Our observations resolve the HMC into three cores whose masses are on the order of 10^1-10^3 Msun. No submm core presents detectable free-free emission in the centimeter regime, but they appear to be associated with masers and thermal line emission from complex organic molecules. Towards the most massive core, SMA1, the CH3CN (18_K-17_K) lines reveal hints of rotation about the axis of a jet/outflow traced by H2O maser and H13CO+ (1--0) line emission. Inverse P-Cygni profiles of the 13CO (3--2) and C18O (3--2) lines seen towards SMA1 indicate that the central high-mass (proto)star(s) is (are) still gaining mass with an accretion rate $ge 3 ~10^{-3}$ Msun/yr. Due to the linear scales and the large values of the accretion rate, we hypothesize that we are observing an accretion flow towards a cluster in the making, rather than towards a single massive star.

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

The present paper focuses on the high-mass star-forming region G23.01-0.41. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G23.01-0.41: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the HMC at 1.3 and 3.6 cm. Results: We have detected H2O, CH3OH, and OH maser emission clustered within 2000 AU from the center of a flattened HMC, oriented SE-NW, from which emerges a massive 12CO outflow, elongated NE-SW, extended up to the pc-scale. Although the three maser species show a clearly different spatial and velocity distribution and sample distinct environments around the massive YSO, the spatial symmetry and velocity field of each maser specie can be explained in terms of expansion from a common center, which possibly denotes the position of the YSO driving the maser motion. Water masers trace both a fast shock (up to 50 km/s) closer to the YSO, powered by a wide-angle wind, and a slower (20 km/s) bipolar jet, at the base of the large-scale outflow. Since the compact free-free emission is found offset from the putative location of the YSO along a direction consistent with that of the maser jet axis, we interpret the radio continuum in terms of a thermal jet. The velocity field of methanol masers can be explained in terms of a composition of slow (4 km/s in amplitude) motions of radial expansion and rotation about an axis approximately parallel to the maser jet. Finally, the distribution of line of sight velocities of the hydroxyl masers suggests that they can trace gas less dense (n(H2) < 10^6 cm^-3) and more distant from the YSO than that traced by the water and methanol masers, which is expanding toward the observer. (Abridged)

Among the tracers of the earliest phases in the massive star formation process, methanol masers have gained increasing importance. The phenomenological distinction between Class I and II methanol masers is based on their spatial association with objects such as jets, cores, and ultracompact HII regions, but is also believed to correspond to different pumping mechanisms: radiation for Class II masers, collisions for Class I masers. In this work, we have surveyed a large sample of massive star forming regions - 296 objects divided into two groups named 'High' and 'Low' according to their [25-12] and [60-12] IRAS colours - in Class I and II methanol masers. Previous studies indicate that the High sources are likely more evolved. Therefore, the sample can be used to assess the existence of a sequence for the occurrence of Class I and II methanol masers during the evolution of a massive star forming region. We observed the 6 GHz (Class II) CH3OH maser with the Effelsberg 100-m telescope, and the 44 GHz and 95 GHz (Class I) CH3OH masers with the Nobeyama 45-m telescope. We have detected: 55 sources in the Class II line (12 new detections); 27 sources in the 44 GHz Class I line (17 new detections); 11 sources in the 95 GHz Class I line (all except one are new detections). Our statistical analysis shows that the ratio between the detection rates of Class II and Class I methanol masers is basically the same in High and Low sources. Therefore, both masers are equally associated with each evolutionary phase. In contrast, all maser species have about 3 times higher detection rates in High than in Low sources. This might indicate that the phenomena that originate all masers become progressively more active with time, during the earliest evolutionary phases of a high-mass star forming region.

The present paper focuses on the high-mass star-forming region G16.59-0.05. Methods: Using the VLBA and the EVN arrays, we conducted phase-referenced observations of the three most powerful maser species in G16.59-0.05: H2O at 22.2 GHz (4 epochs), CH3OH at 6.7 GHz (3 epochs), and OH at 1.665 GHz (1 epoch). In addition, we performed high-resolution (> 0".1), high-sensitivity (< 0.1 mJy) VLA observations of the radio continuum emission from the star-forming region at 1.3 and 3.6 cm. Results: This is the first work to report accurate measurements of the "relative" proper motions of the 6.7 GHz CH3OH masers. The different spatial and 3-D velocity distribution clearly indicate that the 22 GHz water and 6.7 GHz methanol masers are tracing different kinematic environments. The bipolar distribution of 6.7 GHz maser l.o.s. velocities and the regular pattern of observed proper motions suggest that these masers are tracing rotation around a central mass of about 35 solar masses. The flattened spatial distribution of the 6.7 GHz masers, oriented NW-SE, suggests that they can originate in a disk/toroid rotating around the massive YSO which drives the 12CO(2-1) outflow, oriented NE-SW, observed on arcsec scale. The extended, radio continuum source observed close to the 6.7 GHz masers could be excited by a wide-angle wind emitted from the YSO associated with the methanol masers, and such a wind is proven to be sufficiently energetic to drive the NE-SW 12CO(2-1) outflow. The H2O masers distribute across a region offset about 0".5 to the NW of the CH3OH masers, in the same area where emission of high-density molecular tracers, typical of HMCs, was detected. We postulate that a distinct YSO, possibly in an earlier evolutionary phase than that exciting the methanol masers, is responsible for the excitation of the water masers and the HMC molecular lines. (Abridged)

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

We have searched for massive molecular outflows in a sample of high-mass star forming regions, and we have characterised both the outflow properties and those of their associated molecular clumps. With a sample composed largely of more luminous objects than previous ones, this work complements analogous surveys performed by other authors by adding the missing highest luminosity sources. The sample under study has been selected so as to favour the earliest evolutionary phases of star formation, and is composed of very luminous objects (L_bol > 2x10^4 L_sun and up to ~10^6 L_sun), possibly containing O-type stars. Each source has been mapped in 13CO(2-1) and C18O(2-1) with the IRAM-30m telescope on Pico Veleta (Spain). The whole sample shows high-velocity wings in the 13CO(2-1) spectra, indicative of outflowing motions. In addition, we have obtained outflow maps in 9 of our 11 sources, which display well-defined blue and/or red lobes. For these sources, the outflow parameters have been derived from the line wing 13CO(2-1) emission. An estimate of the clump masses from the C18O(2-1) emission is also provided and found to be comparable to the virial masses. From a comparison between our results and those found by other authors at lower masses, it is clear that the outflow mechanical force increases with the bolometric luminosity of the clump and with the ionising photon rate of the associated HII regions, indicating that high-mass stars drive more powerful outflows. A tight correlation between outflow mass and clump mass is also found. Molecular outflows are found to be as common in massive star forming regions as in low-mass star forming regions. This, added to the detection of a few tentative large-scale rotating structures suggests that high-mass stars may generally form via accretion, as low-mass stars.

Glycolaldehyde is the simplest of the monosaccharide sugars and is directly linked to the origin of life. We report on the detection of glycolaldehyde (CH2OHCHO) towards the hot molecular core G31.41+0.31 through IRAM PdBI observations at 1.4, 2.1, and 2.9 mm. The CH2OHCHO emission comes from the hottest (> 300 K) and densest (>2E8 cm^-3) region closest (< 10^4 AU) to the (proto)stars. The comparison of data with gas-grain chemical models of hot cores suggests for G31.41+0.31 an age of a few 10^5 yr. We also show that only small amounts of CO need to be processed on grains in order for existing hot core gas-grain chemical models to reproduce the observed column densities of glycolaldehyde, making surface reactions the most feasible route to its formation.

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.

We used the CSO 10.4 meter telescope to image the 350 micron and 450 micron continuum and CO J=6-5 line emission of the IRAS 20126+4104 clump. The continuum and line observations show that the clump is isolated over a 4 pc region and has a radius of ~ 0.5 pc. Our analysis shows that the clump has a radial density profile propto r ^{-1.2} for r <~ 0.1 pc and has propto r^{-2.3} for r >~ 0.1 pc which suggests the inner region is infalling, while the infall wave has not yet reached the outer region. Assuming temperature gradient of r^{-0.35}, the power law indices become propto r ^{-0.9} for r < ~0.1 pc and propto r^{-2.0} for r >~ 0.1 pc. Based on a map of the flux ratio of 350micron/450micron, we identify three distinct regions: a bipolar feature that coincides with the large scale CO bipolar outflow; a cocoon-like region that encases the bipolar feature and has a warm surface; and a cold layer outside of the cocoon region. The complex patterns of the flux ratio map indicates that the clump is no longer uniform in terms of temperature as well as dust properties. The CO emission near the systemic velocity traces the dense clump and the outer layer of the clump shows narrow line widths (< ~3 km/s). The clump has a velocity gradient of ~ 2 km/s pc^{-1}, which we interpret as due to rotation of the clump, as the equilibrium mass (~ 200 Msun) is comparable to the LTE mass obtained from the CO line. Over a scale of ~ 1 pc, the clump rotates in the opposite sense with respect to the >~ 0.03 pc disk associated with the (proto)star. This is one of four objects in high-mass and low-mass star forming regions for which a discrepancy between the rotation sense of the envelope and the core has been found, suggesting that such a complex kinematics may not be unusual in star forming regions.

We have observed the deuterated gas in the high-mass star formation region IRAS 05345+3157 at high-angular resolution, in order to determine the morphology and the nature of such gas. We have mapped the N2H+ (1-0) line with the Plateau de Bure Interferometer, and the N2D+ (3-2) and N2H+ (3-2) lines with the Submillimeter Array. The N2D+ (3-2) integrated emission is concentrated in two condensations, with masses of 2-3 and 9 M_sun and diameters of 0.05 and 0.09 pc, respectively. The high deuterium fractionation (0.1) and the line parameters in the N2D+ condensations indicate that they are likely low- to intermediate-mass pre-stellar cores, even though other scenarios are possible.

Context. Water vapor emission at 22 GHz from masers associated with star-forming regions is highly variable. Aims. We present a database of up to 20 years of monitoring of a sample of 43 masers within star-forming regions. The sample covers a large range of luminosities of the associated IRAS source and is representative of the entire population of H2O masers of this type. The database forms a good starting point for any further study of H2O maser variability. Methods. The observations were obtained with the Medicina 32-m radiotelescope, at a rate of 4-5 observations per year. Results. To provide a database that can be easily accessed through the web, we give for each source: plots of the calibrated spectra, the velocity-time-flux density plot, the light curve of the integrated flux, the lower and upper envelopes of the maser emission, the mean spectrum, and the rate of the maser occurrence as a function of velocity. Figures for just one source are given in the text for representative purposes. Figures for all the sources are given in electronic form in the on-line appendix. A discussion of the main properties of the H2O variability in our sample will be presented in a forthcoming paper.

Using the OVRO, Nobeyama, and IRAM mm-arrays, we searched for ``disk''-outflow systems in three high-mass (proto)star forming regions: G16.59-0.05, G23.01-0.41, and G28.87+0.07. These were selected from a sample of NH3 cores associated with OH and H2O maser emission and with no or very faint continuum emission. Our imaging of molecular line (including rotational transitions of CH3CN and 3mm dust continuum emission revealed that these are compact, massive, and hot molecular cores (HMCs), that is likely sites of high-mass star formation prior to the appearance of UCHII regions. All three sources turn out to be associated with molecular outflows from CO and/or HCO+ J=1--0 line imaging. In addition, velocity gradients of 10 -- 100 km/s per pc in the innermost densest regions of the G23.01 and G28.87 HMCs are identified along directions roughly perpendicular to the axes of the corresponding outflows. All the results suggest that these cores might be rotating about the outflow axis, although the contribution of rotation to gravitational equilibrium of the HMCs appears to be negligible. Our analysis indicates that the 3 HMCs are close to virial equilibrium due to turbulent pressure support. Comparison with other similar objects where rotating toroids have been identified so far shows that in our case rotation appears to be much less prominent; this can be explained by the combined effect of unfavorable projection, large distance, and limited angular resolution with the current interferometers.

Context. G24.78+0.08 A1 is a 20 Msun star surrounded by a hypercompact (HC) HII region, driving a CO bipolar outflow, and located at the center of a massive rotating toroid undergoing infall towards the HC region. Recent water maser observations suggest that the HC region is expanding and accretion onto the star is halted. Aims. This study aims to confirm the expansion scenario proposed for the HC region on the basis of recent H2O maser observations. Methods. We carried out continuum VLA observations at 1.3cm and 7mm with the A array plus Pie Town configuration to map the HC region towards G24 A1. Results. The emission of the HC region has been resolved and shows a ring shape structure. The profiles of the emission obtained by taking slices at different angles passing through the barycenter of the HC region confirm the shell structure of the emission. The ratio between the inner and the outer radius of the shell, Ri/Ro, derived fitting the normalized brightness temperature profile passing through the peak of the 7mm emission, is 0.9, which indicates that the shell is thin. The deconvolved outer radius estimated from the fit is 590 AU. These results imply that the HC region in G24 A1 cannot be described in terms of a classical, homogeneous HII region but is instead an ionized shell. This gives support to the model of an expanding wind-driven, ionized shell suggested by the kinematics and distribution of the H2O masers associated with the HC region. According to this model, the HC region is expanding on very short times scales, 21-66 yr.

We have observed several emission lines of two Nitrogen-bearing (C2H5CN and C2H3CN) and two Oxygen-bearing (CH3OCH3 and HCOOCH3) molecules towards a sample of well-known hot molecular cores (HMCs) in order to check whether the chemical differentiation seen in the Orion-HMC and W3(H_2O) between O- and N-bearing molecules is a general property of HMCs. With the IRAM-30m telescope we have observed 12 HMCs in 21 bands, centered at frequencies from 86250 to 258280 MHz. The rotational temperatures obtained range from ~100 to ~150 K in these HMCs. Single Gaussian fits performed to unblended lines show a marginal difference in the line peak velocities of the C2H5CN and CH3OCH3 lines, indicating a possible spatial separation between the region traced by the two molecules. On the other hand, neither the linewidths nor the rotational temperatures and column densities confirm such a result. By comparing the abundance ratio of the pair C2H5CN/C2H3CN with the predictions of theoretical models, we derive that the age of our cores ranges between 3.7 and 5.9x10^{4} yrs. The abundances of C2H5CN and C2H3CN are strongly correlated, as expected from theory which predicts that C2H3CN is formed through gas phase reactions involving C2H5CN. A correlation is also found between the abundances of C2H3CN and CH3OCH3, and C2H5CN and CH3OCH3. In all tracers the fractional abundances increase with the H_2 column density while they are not correlated with the gas temperature.

We mapped the HH212 Class 0 outflow in SiO(2--1, 5--4) and continuum using the PdBI in its extended configurations. The unprecedented angular resolution (down to 0.34") allows accurate comparison with a new, deep H2 image obtained at the VLT. The SiO emission is confined to a highly-collimated bipolar jet (width 0.35") along the outflow axis. The jet can be traced down to within 500 AU of the protostar, in a region that is heavily obscured in H2 images. Where both species are detected, SiO shows the same overall kinematics and structure as H2, indicating that both molecules are tracing the same material. We find that the high-velocity SiO gas near the protostar is not tracing a wide-angle wind but is already confined to a flow inside a narrow cone of half-opening angle < 6 deg.

Most previous high-angular (<0.1 arcsec) resolution studies of molecular masers in high-mass star forming regions (SFRs) have concentrated mainly on either water or methanol masers. While high-angular resolution observations have clarified that water masers originate from shocks associated with protostellar jets, different environments have been proposed in several sources to explain the origin of methanol masers. Tha aim of the paper is to investigate the nature of the methanol maser birthplace in SFRs and the association between the water and methanol maser emission in the same young stellar object. We have conducted phase-reference Very Long Baseline Interferometry (VLBI) observations of water and methanol masers toward two high-mass SFRs, Sh 2-255 IR and AFGL 5142. In Sh 2-255 IR water masers are aligned along a direction close to the orientation of the molecular outflow observed on angular scales of 1-10 arcsec, tracing possibly the disk-wind emerging from the disk atmosphere. In AFGL 5142 water masers trace expansion at the base of a protostellar jet, whilst methanol masers are more probably tracing infalling than outflowing gas. The results for AFGL 5142 suggest that water and methanol masers trace different kinematic structures in the circumstellar gas.

Theory predicts and observations confirm that low-mass stars (like the Sun) in their early life grow by accreting gas from the surrounding material. But for stars ~ 10 times more massive than the Sun (~10 M_sun), the powerful stellar radiation is expected to inhibit accretion and thus limit the growth of their mass. Clearly, stars with masses >10 M_sun exist, so there must be a way for them to form. The problem may be solved by non-spherical accretion, which allows some of the stellar photons to escape along the symmetry axis where the density is lower. The recent detection of rotating disks and toroids around very young massive stars has lent support to the idea that high-mass (> 8 M_sun) stars could form in this way. Here we report observations of an ammonia line towards a high-mass star forming region. We conclude from the data that the gas is falling inwards towards a very young star of ~20 M_sun, in line with theoretical predictions of non-spherical accretion.

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.

We report on new aspects of the star-forming region S235AB revealed through high-resolution observations at radio and mid-infrared wavelengths. Using the Very Large Array, we carried out sensitive observations of S235AB in the cm continuum (6, 3.6, 1.3, and 0.7) and in the 22 GHz water maser line. These were complemented with Spitzer Space Telescope Infrared Array Camera archive data to clarify the correspondence between radio and IR sources. We made also use of newly presented data from the Medicina water maser patrol, started in 1987, to study the variability of the water masers found in the region. S235A is a classical HII region whose structure is now well resolved. To the south, no radio continuum emission is detected either from the compact molecular core or from the jet-like structure observed at 3.3 mm, suggesting emission from dust in both cases. We find two new compact radio continuum sources (VLA-1 and VLA-2) and three separate maser spots. VLA-1 coincides with one of the maser spots and with a previously identified IR source (M1). VLA-2 lies towards S235B and represents the first radio detection from this peculiar nebula that may represent an ionized wind from a more evolved star. The two other maser spots coincide with an elongated structure previously observed within the molecular core in the C34S line. This structure is perpendicular to a bipolar molecular outflow observed in HCO+(1-0) and may trace the associated equatorial disk. The Spitzer images reveal a red object towards the molecular core. This is the most viable candidate for the embedded source originating the outflow and maser phenomenology. The picture emerging from these and previous data shows the extreme complexity of a small (< 0.5 pc) star-forming region where widely different stages of stellar evolution are present.

Disks are a natural outcome of the star formation process in which they play a crucial role. Luminous, massive stars of spectral type earlier than B4 are likely to be those that benefit most from the existence of accretion disks, which may significantly reduce the effect of radiation pressure on the accreting material. The scope of the present contribution is to review the current knowledge about disks in young high-mass (proto)stars and discuss their implications. The issues of disk stability and lifetime are also discussed. We conclude that for protostars of less than ~20 solar masses, disks with mass comparable to that of the central star are common. Above this limit the situation is unclear and there are no good examples of proto O4-O8 stars surrounded by accretion disks: in these objects only huge, massive, toroidal, non-equilibrium rotating structures are seen. It is clear on the other hand that the observed disks in stars of 10-20 solar masses are likely to be unstable and with short lifetimes.

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.

We report the results of a multi-epoch survey of water maser observations at 22.2 GHz with the Medicina radiotelescope from 44 bright rimmed clouds (BRCs) of the northern hemisphere identified by Sugitani et al. (1989) as potential sites of star formation. The data span 16 years of observations and allow to draw conclusions about the maser detection rate in this class of objects. In spite of the relatively high far-infrared luminosities of the embedded sources ($L_{\rm FIR}\ga 10^2$ L$_\odot$), H$_2$O maser emission was detected towards three globules only. Since the occurrence of water masers is higher towards bright IRAS sources, the lack of frequent H$_2$O maser emission is somewhat surprising if the suggestion of induced intermediate- and high-mass star formation within these globules is correct. The maser properties of two BRCs are characteristic of exciting sources of low-mass, while the last one (BRC~38) is consistent with an intermediate-mass object. We argue that most BRCs host young stellar objects of low-luminosity, likely in an evolutionary phase later than the protostellar Class 0 sources, and that a significant contribution to the observd IRAS luminosity comes from warm dust heated by the radiation from the bright rim.

We report on the first multi-epoch, phase referenced VLBI observations of the water maser emission in a high-mass protostar associated with a disk-jet system. The source under study, IRAS 20126+4104, has been extensively investigated in a large variety of tracers, including water maser VLBA data acquired by us three years before the present observations. The new findings fully confirm the interpretation proposed in our previous study, namely that the maser spots are expanding from a common origin coincident with the protostar. We also demonstrate that the observed 3-D velocities of the maser spots can be fitted with a model assuming that the spots are moving along the surface of a conical jet, with speed increasing for increasing distance from the cone vertex. We also present the results of single-dish monitoring of the water maser spectra in IRAS 20126+4104. These reveal that the peak velocity of some maser lines decreases linearly with time. We speculate that such a deceleration could be due to braking of the shocks from which the maser emission originates, due to mass loading at the shock front or dissipation of the shock energy.

We present the results of high angular resolution millimeter observations of gas and dust toward G31.41+0.31 and G24.78+0.08, two high-mass star forming regions where four rotating massive toroids have been previously detected by Beltran et al. (2004). The CH3CN (12-11) emission of the toroids in G31.41+0.31 and core A1 in G24.78+0.08 has been modeled assuming that it arises from a disk-like structure seen edge-on, with a radial velocity field. For G31.41+0.31 the model properly fits the data for a velocity v_rot~1.7 km/s at the outer radius R_out~13400 AU and an inner radius R_inn~1340 AU, while for core A1 in G24.78+0.08 the best fit is obtained for v_rot~2.0 km/s at R_out~7700 AU and R_inn~2300 AU. Unlike the rotating disks detected around less luminous stars, these toroids are not undergoing Keplerian rotation. From the modeling itself, however, it is not possible to distinguish between constant rotation or constant angular velocity, since both velocity fields suitably fit the data. The best fit models have been computed adopting a temperature gradient of the type T proportional R^{-3/4}, with a temperature at the outer radius T_out~100 K for both cores. The M_dyn needed for equilibrium derived from the models is much smaller than the mass of the cores, suggesting that such toroids are unstable and undergoing gravitational collapse. The collapse is also supported by the CH3^{13}CN or CH3CN line width measured in the cores, which increases toward the center of the toroids. The estimates of v_inf and \dot M_acc are ~2 km/s and 3x10^{-2} M_sun/yr for G31.41+0.31, and ~1.2 km/s and ~9x10^{-3} M_sun/yr for G24.78+0.08 A1. Such large accretion rates could weaken the effect of stellar winds and radiation pressure and allow further accretion on the star.

Using the Owens Valley and Nobeyama Radio Observatory interferometers, we carried out an unbiased search for hot molecular cores and ultracompact UC HII regions toward the high-mass star forming region G19.61--0.23. In addition, we performed 1.2 mm imaging with SIMBA, and retrieved 3.5 and 2 cm images from the VLA archive data base. The newly obtained 3 mm image brings information on a cluster of high-mass (proto)stars located in the innermost and densest part of the parsec scale clump detected in the 1.2 mm continuum. We identify a total of 10 high-mass young stellar objects: one hot core (HC) and 9 UC HII regions, whose physical parameters are obtained from model fits to their continuum spectra. The ratio between the current and expected final radii of the UC \HII regions ranges from 0.3 to 0.9, which leaves the possibility that all O-B stars formed simultaneously. Under the opposite assumption -- namely that star formation occurred randomly -- we estimate that HC lifetime is less than $\sim$1/3 of that of UCHII regions on the basis of the source number ratio between them.

We present VLBI observations of the H$_2$O maser emission towards a selection of massive young stellar objects (YSOs). In one of these sources, IRAS20126+4104, the measured proper motions confirm that the water masers spots are tracing the very root (around 100 AU) of a conical bipolar outflow emerging near the position of the embedded YSO, and are an excellent tool to investigate the structure and kinematics of the outflow/jet. We also present preliminary results of a first epoch EVN observations of a selection of 5 high-mass YSO to assess the precise association of maser spots and molecular outflows in these sources.

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

An overview is given of the analysis of more than a decade of H2O maser data from our monitoring program. We find the maser emission to generally depend on the luminosity of the YSO as well as on the geometry of the SFR. There appears to be a threshold luminosity of a few times 10**4 Lsol above and below which we find different maser characteristics.

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

The S235A-B star forming region has been extensively observed in the past from the radio to the near-IR, but what was happening in the immediate surroundings of the water maser, placed in between the two nebulosities, was still unclear because of insufficient resolution especially in the spectral range from the Far-IR to the mm, even though there were sound indications that new young stellar objects (YSOs) are being formed there. We present here new high resolution maps at mm wavelengths in different molecules (HCO+, C34S, H2CS, SO2 and CH3CN), as well as in the 1.2 and 3.3 mm continuum obtained with the Plateau de Bure interferometer, and JCMT observations at 450 micron and 850 micron that unambiguously reveal the presence of new YSOs placed in between the two HII regions S235A and S235B and associated with the water maser. A molecular core and an unresolved source in the mm and in the sub-mm are centred on the maser, with indication of mass infall onto the core. Two molecular bipolar outflows and a jet originate from the same position. Weak evidence is found for a molecular rotating disk perpendicular to the direction of the main bipolar outflow. The derived parameters indicate that one of the YSOs is an intermediate luminosity object (L~1000 Lsun) in a very early evolutionary phase, embedded in a molecular core of ~100 Msun, with a temperature of 30 K. The main source of energy for the YSO could come from gravitational infall, thus making of this YSO a rare example of intermediate luminosity protostar representing a link between the earliest evolutionary phases of massive stars and low mass protostars of class 0-I.

We report on the detection of four rotating massive disks in two regions of high-mass star formation. The disks are perpendicular to known bipolar outflows and turn out to be unstable but long lived. We infer that accretion onto the embedded (proto)stars must proceed through the disks with rates of ~10E-2 Msun/yr.

We report on a water maser survey towards a sample of 27 planetary nebulae (PNe) using the Robledo de Chavela and Medicina single-dish antennas, as well as the Very Large Array (VLA). Two detections have been obtained: the already known water maser emission in K 3-35, and a new cluster of masers in IRAS 17347-3139. This low rate of detections is compatible with the short life-time of water molecules in PNe (~100 yr). The water maser cluster at IRAS 17347-3139 are distributed on a ellipse of size ~ 0.2" x 0.1", spatially associated with compact 1.3 cm continuum emission (simultaneously observed with the VLA). From archive VLA continuum data at 4.9, 8.4, and 14.9 GHz, a spectral index alpha = 0.76 +- 0.03 is derived for this radio source, which is consistent with either a partially optically thick ionized region or with an ionized wind. However, the latter scenario can be ruled out on mass-loss considerations, thus indicating that this source is probably a young PN. The spatial distribution and the radial velocities of the water masers are suggestive of a rotating and expanding maser ring, tracing the innermost regions of a torus formed at the end of the AGB phase. Given that the 1.3 cm continuum emission peak is located near one of the tips of the major axis of the ellipse of masers, we speculate on a possible binary nature of IRAS 17347-3139, where the radio continuum emission could belong to one of the components and the water masers would be associated with a companion.

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}$.

(Abridged) H2O masers in 14 SFRs have been monitored once every 2-3 months for up to 13 years. We investigate the dependence of the overall spectral morphology of the maser emission and its variability on the luminosity of the YSO, and look for bursts and gradients in individual components. We find that higher-luminosity sources tend to be associated with stronger and more stable masers.Higher-luminosity YSOs can excite more emission components over a larger range in velocity, yet the emission that dominates the spectra is at a velocity very near that of the molecular cloud in which the objects are embedded. For Lfir > 3E+04 Lo the maser emission becomes increasingly structured and more extended in velocity with increasing Lfir. Below this limit the maser emission shows the same variety of morphologies, but without a clear dependence on Lfir and with a smaller velocity extent. Also, for sources with Lfir above this limit, the water maser is always present above the 5sigma-level; below it, the typical 5sigma detection rate is 75-80%. We conclude that there must be a lower limit to Lfir, below which the associated maser is below the detection level most of the time. These results can be understood in terms of scaled versions of similar SFRs with different YSO luminosities, each with many potential sites of maser amplification, which can be excited provided there is sufficient energy to pump them, i.e. the basic pumping process is identical regardless of the YSO luminosity. In SFRs with lower input energies, the conditions of maser amplification are much closer to the threshold conditions, and consequently more unstable. We find indications that the properties of the maser emission may be determined also by the geometry of the SFR, specifically by the beaming and collimation properties of the outflow driven by the YSO.

We have used the IRAM 30-m and FCRAO 14-m telescopes to observe the molecular clumps associated with 12 ultracompact (UC) HII regions in the J=6-5, 8-7 and 13-12 rotational transitions of methyl-acetylene (CH3C2H). Under the assumption of LTE and optically thin emission, we have derived temperature estimates ranging from 30 to 56 K. We estimate that the clumps have diameters of 0.2-1.6 pc, H_2 densities of 10^5-10^6 {cm^{-3}}, and masses of 10^2-2 10^4 M_\odot. We compare these values with those obtained by other authors from different molecular tracers and find that the H_2 density and the temperature inside the clumps vary respectively like n_{H_2} ~ R^{-2.6} and T ~ R^{-0.5}, with R distance from the centre. We also find that the virial masses of the clumps are ~3 times less than those derived from the CH3C2H column densities: we show that a plausible explanation is that magnetic fields play an important role to stabilise the clumps, which are on the verge of gravitational collapse. Finally, we show that the CH3C2H line width increases for decreasing distance from the clump centre: this effect is consistent with infall in the inner regions of the clumps. We conclude that the clumps around UC HII regions are likely to be transient (~10^(5) yr) entities, remnants of isothermal spheres currently undergoing gravitational collapse: the high mass accretion rates (~10^{-2} M_\odot yr^{-1}) lead to massive star formation at the centre of such clumps.

We present the results of 1.3 mm observations of the Crab Nebula, performed with the MPIfR bolometer arrays at the IRAM 30-m telescope. The maps obtained, of unprecedented quality at these wavelengths, allow a direct comparison with high-resolution radio maps. Although the spatial structure of the Crab Nebula does not change much from radio to millimetre wavelengths, we have detected significant spatial variations of the spectral index between 20 cm and 1.3 mm. The main effect is a spectral flattening in the inner region, which can be hardly explained just in terms of the evolution of a single population of synchrotron emitting electrons. We propose instead that this is the result of the emergence of a second synchrotron component, that we have tried to extract from the data. Shape and size of this component resemble those of the Crab Nebula in X rays. However, while the more compact structure of the Crab Nebula in X rays is commonly regarded as an effect of synchrotron downgrading, it cannot be explained why a similar structure is present also at mm wavelengths, where the electron lifetimes far exceed the nebular age. Our data, combined with published upper limits on spatial variations of the radio spectral index, also imply a low-energy cutoff for the distribution of electrons responsible for this additional synchrotron component. Although no model has been developed so far to explain the details of this component, one may verify that the total number of the electrons responsible for it is in agreement with what predicted by the classical pulsar-wind models, which otherwise are known to fail in accounting for the number of radio emitting electrons. We have also detected a spectral steepening at mm wavelengths in some elongated regions, whose positions match those of radio synchrotron filaments.

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

We present high angular resolution interferometric observations of the 3 and 1.3mm continuum emission, and HCO+(1-0) and SiO(2-1)v=0 lines, obtained with the Owens Valley Radio Observatory millimeter-wave array, toward two hot cores (HCs) associated with two well known ultracompact (UC) HII regions: G29.96-0.02 and G31.41+0.31. These HCs are believed to host young forming massive stars which have been suggested to be surrounded by massive rotating accretion disks. The aim of these new observations is to study the structure and kinematics of the molecular clumps surrounding the HCs and nearby UCHII regions at moderately high angular resolution. Our observations reveal that the clumps within which the HCs and UCHII regions are embedded have a complex kinematical structure. The total mass of the clumps is estimated to be in the range 1000-3000 Msun, consistent with previous findings. Our observations also show compelling evidence that the clump in G29.96-0.02 is contracting onto the HC position, suggesting that the accretion process onto the massive young stellar object embedded in the HC is still ongoing. In these objects the kinematical structure that we observe is also compatible with the presence of a massive rotating disk within the HC, even though we cannot prove this suggestion with our data. The case of G31.41+0.31 is more complicated, and our data, although consistent with the presence of an inner disk and an infalling envelope around it, do not have the required spatial resolution to resolve the different structures.

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.

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.