Search

Context. The fragmentation mode of high-mass molecular clumps and the accretion processes that form the most massive stars (M ≳ 8 M⊙) are still not well understood. A growing number of case studies have found massive young stellar objects (MYSOs) to harbour disk-like structures, painting a picture that the formation of high-mass stars may proceed through disk accretion, similar to that of lower-mass stars. However, the properties of such structures have yet to be uniformly and systematically characterised. Aims. The aim of this work is to uniformly study the kinematic properties of a large sample of MYSOs and characterise the stability of possible circumstellar disks against gravitational fragmentation. Methods. We have undertaken a large observational programme (CORE) making use of interferometric observations from the Northern Extended Millimetre Array (NOEMA) for a sample of 20 luminous (L > 104L⊙) protostellar objects in the 1.37 mm wavelength regime in both continuum and spectral line emission, reaching 0.4″ resolution (800 au at 2 kpc). Results. We present the gas kinematics of the full sample and detect dense gas emission surrounding 15 regions within the CORE sample. Using the dense gas tracer CH3CN, we find velocity gradients across 13 cores perpendicular to the directions of bipolar molecular outflows, making them excellent disk candidates. The extent of the CH3CN emission tracing the disk candidates varies from 1800 to 8500 au. Analysing the free-fall to rotational timescales, we find that the sources are rotationally supported. The rotation profiles of some disk candidates are well described by differential rotation while for others the profiles are poorly resolved. Fitting the velocity profiles with a Keplerian model, we find protostellar masses in the range of ~ 10–25 M⊙. Modelling the level population of CH3CN (12K–11K) K = 0–6 lines, we present temperature maps and find median temperature in the range 70–210 K with a diversity in distributions. Radial profiles of the specific angular momentum (j) for the best disk candidates span a range of 1–2 orders of magnitude, on average ~10−3 km s−1 pc, and they follow j ∝ r1.7, which is consistent with a poorly resolved rotating and infalling envelope-disk model. Studying the Toomre stability of the disk candidates, we find almost all (11 out of 13) disk candidates to be prone to fragmentation due to gravitational instabilities at the scales probed by our observations, as a result of their high disk to stellar mass ratio. In particular, disks with masses greater than ~ 10–20% of the mass of their host (proto)stars are Toomre unstable, and more luminous YSOs tend to have disks that are more massive compared to their host star and hence more prone to fragmentation. Conclusions. In this work, we show that most disk structures around high-mass YSOs are prone to disk fragmentation early in their formation due to their high disk to stellar mass ratio. This impacts the accretion evolution of high-mass protostars which will have significant implications for the formation of the most massive stars. [ABSTRACT FROM AUTHOR]

Context. There is growing evidence that massive stars grow by disk accretion in a similar way to their low-mass counterparts. Early in evolution, these disks can achieve masses that are comparable to the current stellar mass, and therefore the forming disks are highly susceptible to gravitational fragmentation. Aims. We investigate the formation and early evolution of an accretion disk around a forming massive protostar, focussing on its fragmentation physics. To this end, we follow the collapse of a molecular cloud of gas and dust, the formation of a massive protostar, the formation of its circumstellar disk, and the formation and evolution of the disk fragments. Methods. We used a grid-based, self-gravity radiation hydrodynamics code including a sub-grid module for stellar evolution and dust evolution. We purposely do not use a sub-grid module for fragmentation such as sink particles to allow for all paths of fragment formation and destruction, but instead we keep the spatial grid resolution high enough to properly resolve the physical length scales of the problem, namely the pressure scale height and Jeans length of the disk. Simulations are performed on a grid in spherical coordinates with a logarithmic spacing of the grid cells in the radial direction and a cosine distribution of the grid cells in the polar direction, focusing the spatial resolution on the disk midplane. As a consequence, roughly 25% of the total number of grid cells, corresponding to ~26 million grid cells, are used to model the disk physics. These constitute the highest resolution simulations performed up to now on disk fragmentation around a forming massive star with the physics considered here. For a better understanding of the effects of spatial resolution and to compare our high-resolution results with previous lower resolution studies in the literature, we perform the same simulation at five different resolutions, each run differing in resolution from its predecessor by a factor of two. Results. The cloud collapses and a massive (proto)star is formed in its center surrounded by a fragmenting Keplerian-like accretion disk with spiral arms. The fragments have masses of ~1 M⊙, and their continuous interactions with the disk, spiral arms, and other fragments result in eccentric orbits. Fragments form hydrostatic cores surrounded by secondary disks with spiral arms that also produce new fragments. We identified several mechanisms of fragment formation, interaction, and destruction. Central temperatures of the fragments can reach the hydrogen dissociation limit, form second Larson cores, and evolve into companion stars. Based on this, we study the multiplicity predicted by the simulations and find approximately six companions at different distances from the primary: from possible spectroscopic multiples, to companions at distances between 1000 and 2000 au. [ABSTRACT FROM AUTHOR]

We have studied the dust density, temperature, and velocity distributions of the archetypal massive young stellar object (MYSO) AFGL 2591. Given its high luminosity (⁠|$L=2\times 10^5\, \mbox{L$_{\odot}$}$|⁠) and distance (d  = 3.3 kpc), AFGL 2591 has one of the highest |$\sqrt{L}/d$| ratio, giving better resolved dust emission than any other MYSO. As such, this paper provides a template on how to use resolved multiwavelength data and radiative transfer to obtain a well-constrained 2D axisymmetric analytic rotating infall model. We show for the first time that the resolved dust continuum emission from Herschel 70-  μ m observations is extended along the outflow direction, whose origin is explained in part from warm dust in the outflow cavity walls. However, the model can only explain the kinematic features from CH3CN observations with unrealistically low stellar masses (<15 M⊙), indicating that additional physical processes may be playing a role in slowing down the envelope rotation. As part of our three-step continuum and line fitting, we have identified model parameters that can be further constrained by specific observations. High-resolution mm visibilities were fitted to obtain the disc mass (6 M⊙) and radius (2200 au). A combination of SED and near-infrared observations were used to estimate the luminosity and envelope mass together with the outflow cavity inclination and opening angles. [ABSTRACT FROM AUTHOR]

Context. The formation process of high-mass stars (>8 M⊙) is poorly constrained, particularly the effects of clump fragmentation creating multiple systems and the mechanism of mass accretion onto the cores. Aims. We study the fragmentation of dense gas clumps, and trace the circumstellar rotation and outflows by analyzing observations of the high-mass (~500 M⊙) star-forming region IRAS 23033+5951. Methods. Using the Northern Extended Millimeter Array (NOEMA) in three configurations and the IRAM 30 m single-dish telescope at 220 GHz, we probe the gas and dust emission at an angular resolution of ~0.45′′, corresponding to 1900 au. Results. In the millimeter (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 13CO (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 19 M⊙. Rotational temperatures from fitting CH3CN (12K− 11K) spectra are used for estimating the gas temperature and thereby also the disk stability against gravitational fragmentation, utilizing Toomre's Q parameter. Assuming that the candidate disk is in Keplerian rotation about the central stellar object and considering different disk inclination angles, we identify only one candidate disk as being unstable against gravitational instability caused by axisymmetric perturbations. Conclusions. 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. [ABSTRACT FROM AUTHOR]

Studies of the evolution of massive protostars and the evolution of their host molecular cloud cores are commonly treated as separate problems. However, interdependencies between the two can be significant. Here, we study the simultaneous evolution of massive protostars and their host molecular cores using a multi-dimensional radiation hydrodynamics code that incorporates the effects of the thermal pressure and radiative acceleration feedback of the centrally forming protostar. The evolution of the massive protostar is computed simultaneously using the stellar evolution code STELLAR, modified to include the effects of variable accretion. The interdependencies are studied in three different collapse scenarios. For comparison, stellar evolutionary tracks at constant accretion rates and the evolution of the host cores using pre-computed stellar evolutionary tracks are computed. The resulting interdependencies of the protostellar evolution and the evolution of the environment are extremely diverse and depend on the order of events, in particular the time of circumstellar accretion disk formation with respect to the onset of the bloating phase of the star. Feedback mechanisms affect the instantaneous accretion rate and the protostar's radius, temperature, and luminosity on timescales t ⩽ 5 kyr, corresponding to the accretion timescale and Kelvin-Helmholtz contraction timescale, respectively. Nevertheless, it is possible to approximate the overall protostellar evolution in many cases by pre-computed stellar evolutionary tracks assuming appropriate constant average accretion rates. [ABSTRACT FROM AUTHOR]