A distance-limited sample of massive molecular outflows

We have observed 99 mid-infrared-bright, massive young stellar objects and compact H ii regions drawn from the Red MSX source survey in the J = 3−2 transition of 12CO and 13CO, 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 per cent) 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 time-scale. 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 ∼30 M⊙ 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 ∼10−3 M⊙ yr−1 on to 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

Sulphur monoxide exposes a potential molecular disk wind from the planet-hosting disk around HD100546

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₇-6₆ (301.286 GHz) and J=7₈-6₇ (304.078 GHz) are detected in their respective integrated intensity maps. The stacking of these transitions results in a clear 5σ 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‾¹. 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

Spiral arms and instability within the AFGL 4176 mm1 disc

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

Radiative transfer modelling of W33A MM1: 3D structureand dynamics of a complex massive star-forming region

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 Atacama Large Millimeter/submillimeter Array at ≈0.2 arcsec resolution in CH₃CN and dust emission. The MM1 core is fragmented into six compact sources coexisting within ∼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⋆ ≈ 7 M⊙) and disc + envelope (Mgas ≈ 0.4 M⊙), and is the most luminous (LMain ∼ 10⁴ L⊙). 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 centre 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 ∼10³ 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

High-mass star formation at sub-50AU scales

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

ALMA Long Baseline Campaigns: Phase Characteristics of Atmosphere at Long Baselines in the Millimeter and Submillimeter Wavelengths

We present millimeter- and submillimeter-wave phase characteristics measured between 2012 and 2014 of Atacama Large Millimeter/submillimeter Array long baseline campaigns. This paper presents the first detailed investigation of the characteristics of phase fluctuation and phase correction methods obtained with baseline lengths up to ~15 km. The basic phase fluctuation characteristics can be expressed with the spatial structure function (SSF). Most of the 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. These power law slopes do not change with the amount of precipitable water vapor (PWV), but the fitted constants have a weak correlation with PWV, so that the phase fluctuation at a baseline length of 10 km also increases as a function of PWV. The phase correction method using water vapor radiometers (WVRs) works well, especially for the cases where PWV > 1mm, 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 (i.e., when the WVR phase correction works less effectively) 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 in these rare cases 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. Based on the characteristics, 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 ($\gtrsim$200 μm), which is significant for high frequency (>450 GHz or <700 μm) observations. These results suggest the need for an additional phase correction method to reduce the degree of phase fluctuation, such as fast switching, in addition to the WVR phase correction. We simulated the fast switching phase correction method using observations of single quasars, and the result suggests that it works well, with shorter cycle times linearly improving the coherence

From clump to disc scales in W3 IRS4 A case study of the IRAM NOEMA large programme CORE

Context. High-mass star formation typically takes place in a crowded environment, with a higher likelihood of young forming stars affecting and being affected by their surroundings and neighbours, as well as links between different physical scales affecting the outcome. However, observational studies are often focused on either clump or disc scales exclusively. Aims. We explore the physical and chemical links between clump and disc scales in the high-mass star formation region W3 IRS4, a region that contains a number of different evolutionary phases in the high-mass star formation process, as a case-study for what can be achieved as part of the IRAM NOrthern Extended Millimeter Array (NOEMA) large programme named CORE: “Fragmentation and disc formation in high-mass star formation”. Methods. We present 1.4 mm continuum and molecular line observations with the IRAM NOEMA interferometer and 30 m telescope, which together probe spatial scales from ~0.3−20′′ (600−40 000 AU or 0.003−0.2 pc at 2 kpc, the distance to W3). As part of our analysis, we used XCLASS to constrain the temperature, column density, velocity, and line-width of the molecular emission lines. Results. The W3 IRS4 region includes a cold filament and cold cores, a massive young stellar object (MYSO) embedded in a hot core, and a more evolved ultra-compact (UC)H II region, with some degree of interaction between all components of the region that affects their evolution. A large velocity gradient is seen in the filament, suggesting infall of material towards the hot core at a rate of 10−3−10−4 M⊙ yr−1, while the swept up gas ring in the photodissociation region around the UCH II region may be squeezing the hot core from the other side. There are no clear indications of a disc around the MYSO down to the resolution of the observations (600 AU). A total of 21 molecules are detected, with the abundances and abundance ratios indicating that many molecules were formed in the ice mantles of dust grains at cooler temperatures, below the freeze-out temperature of CO (≲35 K). This contrasts with the current bulk temperature of ~50 K, which was obtained from H2CO. Conclusions. CORE observations allow us to comprehensively link the different structures in the W3 IRS4 region for the first time. Our results argue that the dynamics and environment around the MYSO W3 IRS4 have a significant impact on its evolution. This context would be missing if only high resolution or continuum observations were available

Core fragmentation and Toomre stability analysis of W3(H2O): A case study of the IRAM NOEMA large program CORE

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

Fragmentation and disk formation during high-mass star formation: The IRAM NOEMA (Northern Extended Millimeter Array) large program 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

The feedback of an HC HII region on its parental molecular core: The case of core A1 in the star-forming region G24.78+0.08

Context: G24.78+0.08 is a well known high-mass star-forming region, where several molecular cores harboring OB young stellar objects are found inside a clump of size ≈1 pc. This article focuses on the most prominent of these cores, A1, where an intense hypercompact (HC) HII region has been discovered by previous observations. Aims: Our aim is to determine the physical conditions and the kinematics of core A1, and study the interaction of the HII region with the parental molecular core. Methods: We combine ALMA 1.4 mm high-angular resolution (≈0.′′2) observations of continuum and line emission with multi-epoch Very Long Baseline Interferometry data of water 22 GHz and methanol 6.7 GHz masers. These observations allow us to study the gas kinematics on linear scales from 10 to 104 au, and to accurately map the physical conditions of the gas over core A1. Results: The 1.4 mm continuum is dominated by free-free emission from the intense HC HII region (size ≈1000 au) observed to the North of core A1 (region A1N). Analyzing the H30α line, we reveal a fast bipolar flow in the ionized gas, covering a range of LSR velocities (VLSR) of ≈60 km s−1. The amplitude of the VLSR gradient, 22 km s−1 mpc−1, is one of the highest so far observed towards HC HII regions. Water and methanol masers are distributed around the HC HII region in A1N, and the maser three-dimensional (3D) velocities clearly indicate that the ionized gas is expanding at high speed (≥200 km s−1) into the surrounding molecular gas. The temperature distribution (in the range 100–400 K) over core A1, traced with molecular (CH3OCHO, 13CH3CN, 13CH3OH, and CH3CH2CN) transitions with level energy in the range 30 K ≤ Eu/k ≤ 300 K, reflects the distribution of shocks produced by the fast-expansion of the ionized gas of the HII region. The high-energy (550 K ≤ Eu/k ≤ 800 K) transitions of vibrationally excited CH3CN are likely radiatively pumped, and their rotational temperature can significantly differ from the kinetic temperature of the gas. Over core A1, the VLSR maps from both the 1.4 mm molecular lines and the 6.7 GHz methanol masers consistently show a VLSR gradient (amplitude ≈0.3 km s−1 mpc−1) directed approximately S–N. Rather than gravitationally supported rotation of a massive toroid, we interpret this velocity gradient as a relatively slow expansion of core A1