In 2017, an Atacama Large Millimeter/submillimeter Array (ALMA) high-frequency long baseline campaign was organized to test image capabilities with baselines up to 16 km at submillimeter (submm) wavelengths. We investigated image qualities using ALMA receiver Bands 7, 8, 9, and 10 (285-875 GHz) by adopting band-to-band (B2B) phase referencing in which a phase calibrator is tracked at a lower frequency. For B2B phase referencing, it is expected that a closer phase calibrator to a target can be used, comparing to standard in-band phase referencing. In the first step, it is ensured that an instrumental phase offset difference between low- and high-frequency Bands can be removed using a differential gain calibration in which a phase calibrator is certainly detected while frequency switching. In the next step, comparative experiments are arranged to investigate the image quality between B2B and in-band phase referencing with phase calibrators at various separation angles. In the final step, we conducted long baseline imaging tests for a quasar at 289 GHz in Band 7 and 405 GHz in Band 8 and complex structure sources of HL Tau and VY CMa at ~670 GHz in Band 9. The B2B phase referencing was successfully applied, allowing us to achieve an angular resolution of 14x11 and 10x8 mas for HL Tau and VY CMa, respectively. There is a high probability of finding a low-frequency calibrator within 5.4 deg in B2B phase referencing, bright enough to use an 8 s scan length combined with a 7.5 GHz bandwidth.

Water is a key volatile that provides insights into the initial stages of planet formation. The low water abundances inferred from water observations toward low-mass protostellar objects may point to a rapid locking of water as ice by large dust grains during star and planet formation. However, little is known about the water vapor abundance in newly formed planet-forming disks. We aim to determine the water abundance in embedded Keplerian disks through spatially-resolved observations of H$_2^{18}$O lines to understand the evolution of water during star and planet formation. We present H$_2^{18}$O line observations with ALMA and NOEMA millimeter interferometers toward five young stellar objects. NOEMA observed the 3$_{1,3}$ - $2_{2,0}$ line (E$_{\rm up}$ = 203.7 K) while ALMA targeted the $4_{1,4}$ - $3_{2,1}$ line (E$_{\rm up}$ = 322.0 K). Water column densities are derived considering optically thin and thermalized emission. Our observations are sensitive to the emission from the known Keplerian disks around three out of the five Class I objects in the sample. No H$_2^{18}$O emission is detected toward any of our five Class I disks. We report upper limits to the integrated line intensities. The inferred water column densities in Class I disks are N < 10$^{15}$ cm$^{-2}$ on 100 au scales which include both disk and envelope. The upper limits imply a disk-averaged water abundance of $\lesssim 10^{-6}$ with respect to H$_2$ for Class I objects. After taking into account the physical structure of the disk, the upper limit to the water abundance averaged over the inner warm disk with $T>$ 100 K is between 10$^{-7}$ up to 10$^{-5}$. Water vapor is not abundant in warm protostellar envelopes around Class I protostars. Upper limits to the water vapor column densities in Class I disks are at least two orders magnitude lower than values found in Class 0 disk-like structures.

We present high-resolution (30 mas or 130 au at 4.2 kpc) Atacama Large Millimeter/submillimeter Array observations at 1.2 mm of the disc around the forming O-type star AFGL 4176 mm1. The disc (AFGL 4176 mm1-main) has a radius of ~1000 au and contains significant structure, most notably a spiral arm on its redshifted side. We fitted the observed spiral with logarithmic and Archimedean spiral models. We find that both models can describe its structure, but the Archimedean spiral with a varying pitch angle fits its morphology marginally better. As well as signatures of rotation across the disc, we observe gas arcs in CH$_3$CN that connect to other millimetre continuum sources in the field, supporting the picture of interactions within a small cluster around AFGL 4176 mm1-main. Using local thermodynamic equilibrium modelling of the CH$_3$CN K-ladder, we determine the temperature and velocity field across the disc, and thus produce a map of the Toomre stability parameter. Our results indicate that the outer disc is gravitationally unstable and has already fragmented or is likely to fragment in the future, possibly producing further companions. These observations provide evidence that disc fragmentation is one possible pathway towards explaining the high fraction of multiple systems around high-mass stars.

We have conducted a survey of 328 protostars in the Orion molecular clouds with ALMA at 0.87 mm at a resolution of $\sim$0.1" (40 au), including observations with the VLA at 9 mm toward 148 protostars at a resolution of $\sim$0.08" (32 au). This is the largest multi-wavelength survey of protostars at this resolution by an order of magnitude. We use the dust continuum emission at 0.87 mm and 9 mm to measure the dust disk radii and masses toward the Class 0, Class I, and Flat Spectrum protostars, characterizing the evolution of these disk properties in the protostellar phase. The mean dust disk radii for the Class 0, Class I, and Flat Spectrum protostars are 44.9$^{+5.8}_{-3.4}$, 37.0$^{+4.9}_{-3.0}$, and 28.5$^{+3.7}_{-2.3}$ au, respectively, and the mean protostellar dust disk masses are 25.9$^{+7.7}_{-4.0}$, 14.9$^{+3.8}_{-2.2}$, 11.6$^{+3.5}_{-1.9}$ Earth masses, respectively. The decrease in dust disk masses is expected from disk evolution and accretion, but the decrease in disk radii may point to the initial conditions of star formation not leading to the systematic growth of disk radii or that radial drift is keeping the dust disk sizes small. At least 146 protostellar disks (35% out of 379 detected 0.87 mm continuum sources plus 42 non-detections) have disk radii greater than 50 au in our sample. These properties are not found to vary significantly between different regions within Orion. The protostellar dust disk mass distributions are systematically larger than that of Class II disks by a factor of $>$4, providing evidence that the cores of giant planets may need to at least begin their formation during the protostellar phase.

We present a detailed observational and modeling study of the hot core VLA 3 in the high-mass star-forming region AFGL 2591, which is a target region of the NOrthern Extended Millimeter Array (NOEMA) large program CORE. Using NOEMA observations at 1.37 mm with an angular resolution of ~0."42 (1 400 au at 3.33 kpc), we derived the physical and chemical structure of the source. We modeled the observed molecular abundances with the chemical evolution code MUSCLE (MUlti Stage ChemicaL codE). Results. With the kinetic temperature tracers CH3CN and H2CO we observe a temperature distribution with a power-law index of q = 0.41+-0.08. Using the visibilities of the continuum emission we derive a density structure with a power-law index of p = 1.7+-0.1. The hot core spectra reveal high molecular abundances and a rich diversity in complex molecules. The majority of the molecules have an asymmetric spatial distribution around the forming protostar(s), which indicates a complex physical structure on scales < 1 400 au. Using MUSCLE, we are able to explain the observed molecular abundance of 10 out of 14 modeled species at an estimated hot core chemical age of ~21 100 years. In contrast to the observational analysis, our chemical modeling predicts a lower density power-law index of p < 1.4. Reasons for this discrepancy are discussed. Conclusions. Combining high spatial resolution observations with detailed chemical modeling allows us to derive a concise picture of the physical and chemical structure of the famous AFGL 2591 hot core. The next steps are to conduct a similar analysis for the whole CORE sample, and then use this analysis to constrain the chemical diversity in high-mass star formation to a much greater depth.

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.

The large majority of protoplanetary disks have very compact ($\lesssim15\,$AU) continuum emission at mm wavelengths. However, high angular resolution observations that resolve these small disks are still lacking, due to their intrinsically fainter emission compared with large bright disks. In this letter, we present $1.3\,$mm ALMA data of the faint ($\sim10\,$mJy) disk orbiting the TTauri star CX Tau at a resolution of $\sim40\,$mas, $\sim5\,$AU in diameter. The mm-dust disk is compact, with a 68$\%$ enclosing flux radius of 14$\,$AU, and the intensity profile exhibits a sharp drop between 10-20$\,$AU, and a shallow tail between 20-40$\,$AU. No clear signatures of substructure in the dust continuum are observed, down to the same sensitivity level of the DSHARP large program. However, the angular resolution does not allow to detect substructures at a scale of the disk aspect ratio in the inner regions. The radial intensity profile resembles well the inner regions of more extended disks imaged at the same resolution in DSHARP, but with no rings present in the outer disk. No inner cavity is detected, even though the disk has been classified as a transition disk from the spectral energy distribution in the near infrared. The emission of $^{12}$CO is much more extended, with a 68$\%$ enclosing flux radius of 75$\,$AU. The large difference of the mm dust and gas extents ($>5$) strongly points to the occurrence of radial drift, and matches well the predictions of theoretical models.

The formation of planets occurs within protoplanetary disks surrounding young stars, resulting in perturbation of the gas and dust surface densities. Here, we report the first evidence of spatially resolved gas surface density ($\Sigma_{g}$) perturbation towards the AS~209 protoplanetary disk from the optically thin C$^{18}$O ($J=2-1$) emission. The observations were carried out at 1.3~mm with ALMA at a spatial resolution of about 0.3$\arcsec$ $\times$ 0.2$\arcsec$ (corresponding to $\sim$ 38 $\times$ 25 au). The C$^{18}$O emission shows a compact ($\le$60~au), centrally peaked emission and an outer ring peaking at 140~au, consistent with that observed in the continuum emission and, its azimuthally averaged radial intensity profile presents a deficit that is spatially coincident with the previously reported dust map. This deficit can only be reproduced with our physico-thermochemical disk model by lowering $\Sigma_{gas}$ by nearly an order of magnitude in the dust gaps. Another salient result is that contrary to C$^{18}$O, the DCO$^{+}$ ($J=3-2$) emission peaks between the two dust gaps. We infer that the best scenario to explain our observations (C$^{18}$O deficit and DCO$^{+}$ enhancement) is a gas perturbation due to forming-planet(s), that is commensurate with previous continuum observations of the source along with hydrodynamical simulations. Our findings confirm that the previously observed dust gaps are very likely due to perturbation of the gas surface density that is induced by a planet of at least 0.2~M$\rm_{Jupiter}$ in formation. Finally, our observations also show the potential of using CO isotopologues to probe the presence of saturn mass planet(s).

Methods: We observed the high-mass hot core region G351.77-0.54 with ALMA and more than 16km baselines. Results: At a spatial resolution of 18/40au (depending on the distance), we identify twelve sub-structures within the inner few thousand au of the region. The brightness temperatures are high, reaching values greater 1000K, signposting high optical depth toward the peak positions. Core separations vary between sub-100au to several 100 and 1000au. The core separations and approximate masses are largely consistent with thermal Jeans fragmentation of a dense gas core. Due to the high continuum optical depth, most spectral lines are seen in absorption. However, a few exceptional emission lines are found that most likely stem from transitions with excitation conditions above1000K. Toward the main continuum source, these emission lines exhibit a velocity gradient across scales of 100-200au aligned with the molecular outflow and perpendicular to the previously inferred disk orientation. While we cannot exclude that these observational features stem from an inner hot accretion disk, the alignment with the outflow rather suggests that it stems from the inner jet and outflow region. The highest-velocity features are found toward the peak position, and no Hubble-like velocity structure can be identified. Therefore, these data are consistent with steady-state turbulent entrainment of the hot molecular gas via Kelvin-Helmholtz instabilities at the interface between the jet and the outflow. Conclusions: Resolving this high-mass star-forming region at sub-50au scales indicates that the hierarchical fragmentation process in the framework of thermal Jeans fragmentation can continue down to the smallest accessible spatial scales. Velocity gradients on these small scales have to be treated cautiously and do not necessarily stem from disks, but may be better explained with outflow emission.

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.

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.

The growth of dust grains in protoplanetary disks is a necessary first step towards planet formation. This growth has been inferred via observations of thermal dust emission towards mature protoplanetary systems (age >2 million years) with masses that are, on average, similar to Neptune3. In contrast, the majority of confirmed exoplanets are heavier than Neptune. Given that young protoplanetary disks are more massive than their mature counterparts, this suggests that planet formation starts early, but evidence for grain growth that is spatially and temporally coincident with a massive reservoir in young disks remains scarce. Here, we report observations on a lack of emission of carbon monoxide isotopologues within the inner ~15 au of a very young (age ~100,000 years) disk around the Solar-type protostar TMC1A. By using the absence of spatially resolved molecular line emission to infer the gas and dust content of the disk, we conclude that shielding by millimeter-size grains is responsible for the lack of emission. This suggests that grain growth and millimeter-size dust grains can be spatially and temporally coincident with a mass reservoir sufficient for giant planet formation. Hence, planet formation starts during the earliest, embedded phases in the life of young stars.

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.

Recent measurements carried out at comet 67P/C-G with the ${\it Rosetta}$ probe revealed that molecular oxygen, O$_2$, is the fourth most abundant molecule in comets. Models show that O$_2$ is likely of primordial nature, coming from the interstellar cloud from which our Solar System was formed. However, gaseous O$_2$ is an elusive molecule in the interstellar medium with only one detection towards quiescent molecular clouds, in the $\rho$ Oph A core. We perform a deep search for molecular oxygen, through the $2_1 - 0_1$ rotational transition at 234 GHz of its $^{16}$O$^{18}$O isotopologue, towards the warm compact gas surrounding the nearby Class 0 protostar IRAS 16293--2422 B with the ALMA interferometer. The targeted $^{16}$O$^{18}$O transition is surrounded by two brighter transitions at $\pm 1$ km s$^{-1}$ relative to the expected $^{16}$O$^{18}$O transition frequency. After subtraction of these two transitions, residual emission at a 3$\sigma$ level remains, but with a velocity offset of $0.3 - 0.5$ km s$^{-1}$ relative to the source velocity, rendering the detection "tentative". We derive the O$_2$ column density for two excitation temperatures $T_{\rm ex}$ of 125 and 300 K, as indicated by other molecules, in order to compare the O$_2$ abundance between IRAS16293 and comet 67P/C-G. Assuming that $^{16}$O$^{18}$O is not detected and using methanol CH$_3$OH as a reference species, we obtain a [O$_2$]/[CH$_3$OH] abundance ratio lower than $2-5$, depending on the assumed $T_{\rm ex}$, a three to four times lower abundance than the [O$_2$]/[CH$_3$OH] ratio of $5-15$ found in comet 67P/C-G. Such a low O$_2$ abundance could be explained by the lower temperature of the dense cloud precursor of IRAS16293 with respect to the one at the origin of our Solar System that prevented an efficient formation of O$_2$ in interstellar ices.

Understanding the formation of massive stars is one of the unsolved problems in modern astronomy. The main difficulty is that the intense radiation from the high-luminosity stars and the thermal pressure from the resulting ionized gas (both insignificant for low-mass stars) may be able to reverse the accretion flow and prevent the star from accreting fresh material. Such feedback effects can naturally be mitigated if accretion proceeds through discs, which is the established mechanism to form sun-like stars. However, recent 3D MHD simulations have shown that accretion on 1000 au scales is through filaments rather than a large disc. This theoretical prediction has never been confirmed via observations owing to the poor linear resolution of previous studies (>1000 au). Here we present the first observational evidence that mass assembly in young high-mass stars forming in protoclusters is predominantly asymmetric and disordered. In particular, we observed the innermost regions around three deeply embedded high-mass protostars with very high spatial resolution (~100 au). We identified multiple massive (several solar masses), warm (50-150 Kelvin) filamentary streamers pointing onto the central sources, which we interpret as multi-directional accretion channels. These structures inhibit the formation of a large, steady disc. Nevertheless, the identification of fast collimated outflows in the three observed systems indicates that (non-steady) compact discs may be present (we measure upper limits on their radii of <80 for one object and <350 astronomical units for the remaining two objects). Our finding contrasts with the simplified classic paradigm of an ordered (and stable) disc/jet system and provides an experimental confirmation of a multi-directional and unsteady accretion model for massive star formation supported by recent 3D (magneto)hydrodynamic simulations.

A wide variety of ring-like dust structures has been detected in protoplanetary disks, but their origin and frequency are still unclear. We characterize the structure of an extended, multi-ringed disk discovered serendipitously in the ALMA Lupus disk survey and put it in the context of the Lupus disk population. ALMA observations in Band 6 at 234 GHz and Band 7 at 328 GHz at 0.3" resolution toward the K6 star V1094 Sco in Lupus III are presented, and its disk structure is analyzed. The spectral index $\alpha_{mm}$ is determined in the inner 150 AU of the disk. The ALMA continuum data show a very extended disk with two gap/ring pairs. The gaps are located at 100 AU and 170 AU, the bright rings at 130 AU and 220 AU. Continuum emission is detected out to a 300 AU distance, similar to IM Lup but a factor of 5 larger than typically found for Lupus disks at this sensitivity and resolution. The bright central region of the disk (within 35 AU) is possibly optically thick at 1 mm wavelengths, and has a brightness temperature of only 13 K. The spectral index increases between the inner disk and the first ring, at the location of the first gap. Due to the low temperature of the disk midplane, snow lines can be excluded as the drivers behind the ring and gap formation in this disk. Disks the size of V1094 Sco are rare, and only 2.1+-1.5% of disks in Lupus show continuum emission beyond 200 AU. Possible connections between the large primordial disk population, transition disks, and exoplanets are discussed.

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.

We present a composite model and radiative transfer simulations of the massive star forming core W33A MM1. The model was tailored to reproduce the complex features observed with ALMA at $\approx 0.2$ arcsec resolution in CH$_3$CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within $\sim 1000$ au. In our models, three of these compact sources are better represented as disc-envelope systems around a central (proto)star, two as envelopes with a central object, and one as a pure envelope. The model of the most prominent object (Main) contains the most massive (proto)star ($M_\star\approx7~M_\odot$) and disc+envelope ($M_\mathrm{gas}\approx0.4~M_\odot$), and is the most luminous ($L_\mathrm{Main} \sim 10^4~L_\odot$). The model discs are small (a few hundred au) for all sources. The composite model shows that the elongated spiral-like feature converging to the MM1 core can be convincingly interpreted as a filamentary accretion flow that feeds the rising stellar system. The kinematics of this filament is reproduced by a parabolic trajectory with focus at the center of mass of the region. Radial collapse and fragmentation within this filament, as well as smaller filamentary flows between pairs of sources are proposed to exist. Our modelling supports an interpretation where what was once considered as a single massive star with a $\sim 10^3$ au disc and envelope, is instead a forming stellar association which appears to be virialized and to form several low-mass stars per high-mass object.

We present the results of an outflow and infall survey towards a distance limited sample of 31 massive star forming regions drawn from the RMS survey. The presence of young, active outflows is identified from SiO (8-7) emission and the infall dynamics are explored using HCO$^+$/H$^{13}$CO$^+$ (4-3) emission. We investigate if the infall and outflow parameters vary with source properties, exploring whether regions hosting potentially young active outflows show similarities or differences with regions harbouring more evolved, possibly momentum driven, "fossil" outflows. SiO emission is detected towards approximately 46% of the sources. When considering sources with and without an SiO detection (i.e. potentially active and fossil outflows respectively), only the $^{12}$CO outflow velocity shows a significant difference between samples, indicating SiO is more prevalent towards sources with higher outflow velocities. Furthermore, we find the SiO luminosity increases as a function of the Herschel 70$\mu$m to WISE 22$\mu$m flux ratio, suggesting the production of SiO is prevalent in younger, more embedded regions. Similarly, we find tentative evidence that sources with an SiO detection have a smaller bolometric luminosity-to-mass ratio, indicating SiO (8-7) emission is associated with potentially younger regions. We do not find a prevalence towards sources displaying signatures of infall in our sample. However, the higher energy HCO$^+$ transitions may not be the best suited tracer of infall at this spatial resolution in these regions.

Sulphur-bearing volatiles are observed to be significantly depleted in interstellar and circumstellar regions. This missing sulphur is postulated to be mostly locked up in refractory form. With ALMA we have detected sulphur monoxide (SO), a known shock tracer, in the HD 100546 protoplanetary disk. Two rotational transitions: $J=7_{7}-6_{6}$ (301.286 GHz) and $J=7_8-6_7$ (304.078 GHz) are detected in their respective integrated intensity maps. The stacking of these transitions results in a clear 5$\sigma$ detection in the stacked line profile. The emission is compact but is spectrally resolved and the line profile has two components. One component peaks at the source velocity and the other is blue-shifted by 5 km s$^{-1}$. The kinematics and spatial distribution of the SO emission are not consistent with that expected from a purely Keplerian disk. We detect additional blue-shifted emission that we attribute to a disk wind. The disk component was simulated using LIME and a physical disk structure. The disk emission is asymmetric and best fit by a wedge of emission in the north east region of the disk coincident with a `hot-spot' observed in the CO $J=3-2$ line. The favoured hypothesis is that a possible inner disk warp (seen in CO emission) directly exposes the north-east side of the disk to heating by the central star, creating locally the conditions to launch a disk wind. Chemical models of a disk wind will help to elucidate why the wind is particularly highlighted in SO emission and whether a refractory source of sulphur is needed. An alternative explanation is that the SO is tracing an accretion shock from a circumplanetary disk associated with the proposed protoplanet embedded in the disk at 50 au. We also report a non-detection of SO in the protoplanetary disk around HD 97048.

The Atacama Large millimetre/submillimetre Array (ALMA) makes use of water vapour radiometers (WVR), which monitor the atmospheric water vapour line at 183 GHz along the line of sight above each antenna to correct for phase delays introduced by the wet component of the troposphere. The application of WVR derived phase corrections improve the image quality and facilitate successful observations in weather conditions that were classically marginal or poor. We present work to indicate that a scaling factor applied to the WVR solutions can act to further improve the phase stability and image quality of ALMA data. We find reduced phase noise statistics for 62 out of 75 datasets from the long-baseline science verification campaign after a WVR scaling factor is applied. The improvement of phase noise translates to an expected coherence improvement in 39 datasets. When imaging the bandpass source, we find 33 of the 39 datasets show an improvement in the signal-to-noise ratio (S/N) between a few to ~30 percent. There are 23 datasets where the S/N of the science image is improved: 6 by <1%, 11 between 1 and 5%, and 6 above 5%. The higher frequencies studied (band 6 and band 7) are those most improved, specifically datasets with low precipitable water vapour (PWV), <1mm, where the dominance of the wet component is reduced. Although these improvements are not profound, phase stability improvements via the WVR scaling factor come into play for the higher frequency (>450 GHz) and long-baseline (>5km) observations. These inherently have poorer phase stability and are taken in low PWV (<1mm) conditions for which we find the scaling to be most effective. A promising explanation for the scaling factor is the mixing of dry and wet air components, although other origins are discussed. We have produced a python code to allow ALMA users to undertake WVR scaling tests and make improvements to their data.

We targeted the massive star forming region W33A using the Atacama Large Sub/Millimeter Array (ALMA) in band 6 (230 GHz) and 7 (345 GHz) to search for a sub-1000au disc around the central O-type massive young stellar object (MYSO) W33A MM1-Main. Our data achieves a resolution of ~0.2" (~500au) and resolves the central core, MM1, into multiple components and reveals complex and filamentary structures. There is strong molecular line emission covering the entire MM1 region. The kinematic signatures are inconsistent with only Keplerian rotation although we propose that the shift in the emission line centroids within ~1000au of MM1-Main could hint at an underlying compact disc with Keplerian rotation. We cannot however rule out the possibility of an unresolved binary or multiple system. A putative smaller disc could be fed by the large scale spiral `feeding filament' we detect in both gas and dust emission. We also discuss the nature of the now-resolved continuum sources.

This paper presents the first detailed investigation of the characteristics of mm/submm phase fluctuation and phase correction methods obtained using ALMA with baseline lengths up to ~15 km. Most of the spatial structure functions (SSFs) show that the phase fluctuation increases as a function of baseline length, with a power-law slope of ~0.6. In many cases, we find that the slope becomes shallower (average of ~0.2-0.3) at baseline lengths longer than ~1 km, namely showing a turn-over in SSF. The phase correction method using water vapor radiometers (WVRs) works well, especially for the cases where PWV >1 mm, which reduces the degree of phase fluctuations by a factor of two in many cases. However, phase fluctuations still remain after the WVR phase correction, suggesting the existence of other turbulent constituent that cause the phase fluctuation. This is supported by occasional SSFs that do not exhibit any turn-over; these are only seen when the PWV is low or after WVR phase correction. This means that the phase fluctuation caused by this turbulent constituent is inherently smaller than that caused by water vapor. Since there is no turn-over in the SSF up to the maximum baseline length of ~15 km, this turbulent constituent must have scale height of 10 km or more, and thus cannot be water vapor, whose scale height is around 1 km. This large scale height turbulent constituent is likely to be water ice or a dry component. Excess path length fluctuation after the WVR phase correction at a baseline length of 10 km is large (>200 micron), which is significant for high frequency (>450 GHz or <700 micron) observations. These results suggest the need for an additional phase correction method, such as fast switching, in addition to the WVR phase correction. We simulated the fast switching, and the result suggests that it works well, with shorter cycle times linearly improving the coherence.

Transitional disks show a lack of excess emission at infrared wavelengths due to a large dust cavity, that is often corroborated by spatially resolved observations at ~ mm wavelengths. We present the first spatially resolved ~ mm-wavelength images of the disk around the Herbig Ae/Be star, HD 97048. Scattered light images show that the disk extends to ~640 au. The ALMA data reveal a circular-symmetric dusty disk extending to ~350 au, and a molecular disk traced in CO J=3-2 emission, extending to ~750 au. The CO emission arises from a flared layer with an opening angle ~ 30 deg - 40 deg. HD 97048 is another source for which the large (~ mm-sized) dust grains are more centrally concentrated than the small (~ {\mu}m-sized) grains and molecular gas, likely due to radial drift. The images and visibility data modelling suggests a decrement in continuum emission within ~50 au, consistent with the cavity size determined from mid-infrared imaging (34 +/- 4 au). The extracted continuum intensity profiles show ring-like structures with peaks at ~50, 150, and 300 au, with associated gaps at ~100 and 250 au. This structure should be confirmed in higher-resolution images (FWHM ~ 10 - 20 au). These data confirm the classification of HD 97048 as a transitional disk that also possesses multiple ring-like structures in the dust continuum emission. Additional data are required at multiple and well-separated frequencies to fully characterise the disk structure, and thereby constrain the mechanism(s) responsible for sculpting the HD 97048 disk.

Outflowing motions, whether a wind launched from the disk, a jet launched from the protostar, or the entrained molecular outflow, appear to be an ubiquitous feature of star formation. These outwards motions have a number of root causes, and how they manifest is intricately linked to their environment as well as the process of star formation itself. Using the ALMA Science Verification data of HL Tau, we investigate the high velocity molecular gas being removed from the system as a result of the star formation process. We aim to place these motions in context with the optically detected jet, and the disk. With these high resolution ($\sim 1"$) ALMA observations of CO (J=1-0), we quantify the outwards motions of the molecular gas. We find evidence for a bipolar outwards flow, with an opening angle, as measured in the red-shifted lobe, starting off at 90$^\circ$, and narrowing to 60$^\circ$ further from the disk, likely because of magnetic collimation. Its outwards velocity, corrected for inclination angle is of order 2.4 km s$^{-1}$.

We present 1.3mm Submillimeter Array (SMA) observations at $\sim$3$^{\prime\prime}$ resolution towards the brightest section of the intermediate/massive star forming cluster NGC 2264-C. The millimetre continuum emission reveals ten 1.3mm continuum peaks, of which four are new detections. The observed frequency range includes the known molecular jet/outflow tracer SiO (5-4), thus providing the first high resolution observations of SiO towards NGC 2264-C. We also detect molecular lines of twelve additional species towards this region, including CH$_3$CN, CH$_3$OH, SO, H$_2$CO, DCN, HC$_3$N, and $^{12}$CO. The SiO (5-4) emission reveals the presence of two collimated, high velocity (up to 30kms$^{-1}$ with respect to the systemic velocity) bi-polar outflows in NGC 2264-C. In addition, the outflows are traced by emission from $^{12}$CO, SO, H$_2$CO, and CH$_3$OH. We find an evolutionary spread between cores residing in the same parent cloud. The two unambiguous outflows are driven by the brightest mm continuum cores, which are IR-dark, molecular line weak, and likely the youngest cores in the region. Furthermore, towards the RMS source AFGL 989-IRS1, the IR-bright and most evolved source in NGC 2264-C, we observe no molecular outflow emission. A molecular line rich ridge feature, with no obvious directly associated continuum source, lies on the edge of a low density cavity and may be formed from a wind driven by AFGL 989-IRS1. In addition, 229GHz class I maser emission is detected towards this feature.

We have observed 99 mid-infrared-bright, massive young stellar objects and compact HII regions drawn from the Red MSX source (RMS) survey in the J=3$-$2 transition of $^{12}$CO and $^{13}$CO, using the James Clerk Maxwell Telescope. 89 targets are within 6 kpc of the Sun, covering a representative range of luminosities and core masses. These constitute a relatively unbiased sample of bipolar molecular outflows associated with massive star formation. Of these, 59, 17 and 13 sources (66, 19 and 15 percent) are found to have outflows, show some evidence of outflow, and have no evidence of outflow, respectively. The time-dependent parameters of the high-velocity molecular flows are calculated using a spatially variable dynamic timescale. The canonical correlations between the outflow parameters and source luminosity are recovered and shown to scale with those of low-mass sources. For coeval star formation we find the scaling is consistent with all the protostars in an embedded cluster providing the outflow force, with massive stars up to $\sim$30 M$_{\odot}$ generating outflows. Taken at face value, the results support the model of a scaled-up version of the accretion-related outflow-generation mechanism associated with discs and jets in low-mass objects with time-averaged accretion rates of $\sim$10$^{-3}$ M$_{\odot}$ yr$^{-1}$ onto the cores. However, we also suggest an alternative model, in which the molecular outflow dynamics are dominated by the entrained mass and are unrelated to the details of the acceleration mechanism. We find no evidence that outflows contribute significantly to the turbulent kinetic energy of the surrounding dense cores.

We analyse C$^{18}$O ($J=3-$2) data from a sample of 99 infrared-bright massive young stellar objects (MYSOs) and compact HII regions that were identified as potential molecular-outflow sources in the Red MSX source (RMS) survey. We extract a distance limited (D $<$ 6 kpc) sample shown to be representative of star formation covering the transition between the source types. At the spatial resolution probed, Larson-like relationships are found for these cores, though the alternative explanation, that Larson's relations arise where surface-density-limited samples are considered, is also consistent with our data. There are no significant differences found between source properties for the MYSOs and HII regions, suggesting that the core properties are established prior to the formation of massive stars, which subsequently have little impact at the later evolutionary stages investigated. There is a strong correlation between dust-continuum and C$^{18}$O-gas masses, supporting the interpretation that both trace the same material in these IR-bright sources. A clear linear relationship is seen between the independently established core masses and luminosities. The position of MYSOs and compact HII regions in the mass-luminosity plane is consistent with the luminosity expected a cluster of protostars when using a $\sim$40 percent star-formation efficiency and indicates that they are at a similar evolutionary stage, near the end of the accretion phase.

A major goal of the Atacama Large Millimeter/submillimeter Array (ALMA) is to make accurate images with resolutions of tens of milliarcseconds, which at submillimeter (submm) wavelengths requires baselines up to ~15 km. To develop and test this capability, a Long Baseline Campaign (LBC) was carried out from September to late November 2014, culminating in end-to-end observations, calibrations, and imaging of selected Science Verification (SV) targets. This paper presents an overview of the campaign and its main results, including an investigation of the short-term coherence properties and systematic phase errors over the long baselines at the ALMA site, a summary of the SV targets and observations, and recommendations for science observing strategies at long baselines. Deep ALMA images of the quasar 3C138 at 97 and 241 GHz are also compared to VLA 43 GHz results, demonstrating an agreement at a level of a few percent. As a result of the extensive program of LBC testing, the highly successful SV imaging at long baselines achieved angular resolutions as fine as 19 mas at ~350 GHz. Observing with ALMA on baselines of up to 15 km is now possible, and opens up new parameter space for submm astronomy.

We study the evolution of circumstellar disks in 22 young (1 to 100 Myr) nearby (within 500 pc) associations over the entire mass spectrum using photometry covering from the optical to the mid-infrared. We compiled a catalog of 2340 spectroscopically-confirmed members of these nearby associations. We analyzed their spectral energy distributions and searched for excess related to the presence of protoplanetary disks in a homogeneous way. Sensitivity limits and spatial completeness were also considered. We derive disk fractions as probed by mid-infrared excess in these regions. The unprecedented size of our sample allows us to confirm the timescale of disk decay reported in the literature and to find new trends. The fraction of excess sources increases systematically if measured at longer wavelengths. Disk percentages derived using different wavelength ranges should therefore be compared with caution. The dust probed at 22-24 um evolves slower than that probed at shorter wavelengths (3.4-12 um). Assuming an exponential decay, we derive a timescale tau=4.2-5.8 Myr at 22-24 um for primordial disks, compared to 2-3 Myr at shorter wavelength (3.4-12 um). Primordial disks disappear around 10 Myr, matching in time a brief increase of the number of 'evolved' disks. The increase in timescale of excess decay at longer wavelength is compatible with inside-out disk clearing scenarios. The increased timescale of decay and larger dispersion in the distribution of disk fractions at 22-24 um suggest that the inner and outer zones evolve differently, the latter potentially following a variety of evolutionary paths. The drop of primordial disks and the coincident rise of evolved disks at 10 Myr are compatible with planet formation theories suggesting that the disappearance of the gas is immediately followed by the dynamical stirring of the disk.

We investigate the circumstellar dust properties of the oxygen-rich bipolar proto-planetary nebula IRAS 18276-1431 by means of two-dimensional radiative transfer simulations of the circumstellar dust shell. The model geometry is assumed to have a torus and an envelope. The parameters of the dust and the dust shell are constrained by comparing the SED and NIR intensity and polarisation data with the models. The polarisation in the envelope reaches 50 -- 60 % and is nearly constant in the H and K_S bands in the observations. This weak wavelength dependence of the polarisation can be reproduced with a grain size distribution function for the torus: 0.05 micron <= a with n(a)=a^{-(p=5.5)}exp(-a/{a_c=0.3 micron}). The power index p is significantly steeper than that for interstellar dust. Similar results have also been found in some other PPNs and suggest that mechanisms that grind down large particles may also have acted when the dust particles formed. The spectral opacity index beta is found to be 0.6+/-0.5 from the millimeter fluxes. This low value indicates the presence of large dust grains in the torus. We discuss two possible dust models for the torus. One has a size distribution function of 1.0 micron <= a <= a_max=5,000.0 micron with n(a)=a^{-(p=2.5)} and the other is 1.0 micron <= a <= a_max=10,000.0 micron with n(a)=a^{-(p=3.5)}. The former has beta of 0.633, but we are not able to find reasonable geometry parameters to fit the SED in the infrared. The latter has beta of 1.12, but reproduces the SED better over a wide wavelength range. With this dust model, the geometric parameters are estimated as follows: the inner and outer radii are 30 AU and 1000 AU and the torus mass is 3.0 M_sun. Assuming an expansion velocity of 15 kms^{-1}, the torus formation time and mass-loss rate are found to be \sim300 yrs and \sim10^{-2}M_sun yr^{-1} respectively.

Massive stars play an important role in many areas of astrophysics, but numerous details regarding their formation remain unclear. In this paper we present and analyse high resolution (R ~ 30,000) near-infrared 2.3 micron spectra of 20 massive young stellar objects from the RMS database, in the largest such study of CO first overtone bandhead emission to date. We fit the emission under the assumption it originates from a circumstellar disc in Keplerian rotation. We explore three approaches to modelling the physical conditions within the disc - a disc heated mainly via irradiation from the central star, a disc heated mainly via viscosity, and a disc in which the temperature and density are described analytically. We find that the models described by heating mechanisms are inappropriate because they do not provide good fits to the CO emission spectra. We therefore restrict our analysis to the analytic model, and obtain good fits to all objects that possess sufficiently strong CO emission, suggesting circumstellar discs are the source of this emission. On average, the temperature and density structure of the discs correspond to geometrically thin discs, spread across a wide range of inclinations. Essentially all the discs are located within the dust sublimation radius, providing strong evidence that the CO emission originates close to the central protostar, on astronomical unit scales. In addition, we show that the objects in our sample appear no different to the general population of MYSOs in the RMS database, based on their near- and mid-infrared colours. The combination of observations of a large sample of MYSOs with CO bandhead emission and our detailed modelling provide compelling evidence of the presence of small scale gaseous discs around such objects, supporting the scenario in which massive stars form via disc accretion.