By matching infrared-selected, massive young stellar objects (MYSOs) and compact HII regions in the RMS survey to massive clumps found in the submillimetre ATLASGAL survey, we have identified ∼1000 embedded young massive stars between 280• < ℓ < 350• and 10Combined with an existing sample of radio-selected methanol masers and compact HII regions, the result is a catalogue of ∼1700 massive stars embedded within ∼1300 clumps located across the inner Galaxy, containing three observationally distinct subsamples, methanol-maser, MYSO and HII-region associations, covering the most important tracers of massive star formation, thought to represent key stages of evolution. We find that massive star formation is strongly correlated with the regions of highest column density in spherical, centrally condensed clumps. We find no significant differences between the three samples in clump structure or the relative location of the embedded stars, which suggests that the structure of a clump is set before the onset of star formation, and changes little as the embedded object evolves towards the main sequence. There is a strong linear correlation between clump mass and bolometric luminosity, with the most massive stars forming in the most massive clumps. We find that the MYSO and HII-region subsamples are likely to cover a similar range of evolutionary stages and that the majority are near the end of their main accretion phase. We find few infrared-bright MYSOs associated with the most massive clumps, probably due to very short pre-main sequence lifetimes in the most luminous sources.

Context. High-mass stars form in clusters, but neither the early fragmentation processes nor the detailed physical processes leading to the most massive stars are well understood. 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 30 m telescope, we have observed in the IRAM large program CORE the 1.37 mm 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.5 kpc and luminosities larger than 104 L⊙. Results. We present the overall survey scope, the selected sample, the observational setup, and the main goals of CORE. Scientifically, we concentrated on the mm continuum emission on scales on the order of 1000 AU. 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 vs. singular regions can be explained by varying initial density structures and/or different initial magnetic field strengths. Conclusions. A large sample of high-mass star-forming regions at high spatial resolution allows us to study the fragmentation properties of young cluster-forming regions. 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. Circumstellar discs around massive stars could mediate the accretion onto the star from the infalling envelope, and could minimize the effects of radiation pressure. Despite such a crucial role, only a few convincing candidates have been provided for discs around deeply embedded O-type (proto)stars. Aims. In order to establish whether disc-mediated accretion is the formation mechanism for the most massive stars, we have searched for circumstellar, rotating discs around a limited sample of six luminous (>10 5 L ) young stellar objects. These objects were selected on the basis of their IR and radio properties in order to maximize the likelihood of association with disc+jet systems. Methods. We used ALMA with ∼0 . 2 resolution to observe a large number of molecular lines typical of hot molecular cores. In this paper we limit our analysis to two disc tracers (methyl cyanide, CH 3 CN, and its isotopologue, 13 CH 3 CN), and an outflow tracer (silicon monoxide, SiO). Results. We reveal many cores, although their number depends dramatically on the target. We focus on the cores that present prominent molecular line emission. In six of these a velocity gradient is seen across the core, three of which show evidence of Keplerian-like rotation. The SiO data reveal clear but poorly collimated bipolar outflow signatures towards two objects only. This can be explained if real jets are rare (perhaps short-lived) in very massive objects and/or if stellar multiplicity significantly affects the outflow structure. For all cores with velocity gradients, the velocity field is analysed through position-velocity plots to establish whether the gas is undergoing rotation with rot ∝ R −α , as expected for Keplerian-like discs. Conclusions. Our results suggest that in three objects we are observing rotation in circumstellar discs, with three more tentative cases, and one core where no evidence for rotation is found. In all cases but one, we find that the gas mass is less than the mass of any embedded O-type star, consistent with the (putative) discs undergoing Keplerian-like rotation. With the caveat of low number statistics, we conclude that the disc detection rate could be sensitive to the evolutionary stage of the young stellar object. In young, deeply embedded sources, the evidence for discs could be weak because of confusion with the surrounding envelope, while in the most evolved sources the molecular component of the disc could have already been dispersed. Only in those objects that are at an intermediate stage of the evolution would the molecular disc be sufficiently prominent and relatively less embedded to be detectable by mm/submm observations.

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. Our aim is to determine the critical parameters in water chemistry and the contribution of water to the oxygen budget by observing and modelling water gas and ice for a sample of eleven low-mass protostars, for which both forms of water have been observed. Methods. A simplified chemistry network, which is benchmarked against more sophisticated chemical networks, is developed that includes the necessary ingredients to determine the water vapour and ice abundance profiles in the cold, outer envelope in which the temperature increases towards the protostar. Comparing the results from this chemical network to observations of water emission lines and previously published water ice column densities, allows us to probe the influence of various agents (e.g., far-ultraviolet (FUV) field, initial abundances, timescales, and kinematics). Results. The observed water ice abundances with respect to hydrogen nuclei in our sample are 30−80 ppm, and therefore contain only 10-30% of the volatile oxygen budget of 320 ppm. The keys to reproduce this result are a low initial water ice abundance after the pre-collapse phase together with the fact that atomic oxygen cannot freeze-out and form water ice in regions with T dust > ∼ 15 K. This requires short prestellar core lifetimes < ∼ 0.1 Myr. The water vapour profile is shaped through the interplay of FUV photodesorption, photodissociation, and freeze-out. The water vapour line profiles are an invaluable tracer for the FUV photon flux and envelope kinematics. Conclusions. The finding that only a fraction of the oxygen budget is locked in water ice can be explained either by a short precollapse time of < ∼ 0.1 Myr at densities of n H ∼ 10 4 cm −3 , or by some other process that resets the initial water ice abundance for the post-collapse phase. A key for the understanding of the water ice abundance is the binding energy of atomic oxygen on ice.