Search

Context. The fragmentation of high-mass star-forming regions depends on a variety of physical parameters, including density, the magnetic field, and turbulent gas properties. Aims. We evaluate the importance of the density and magnetic field structures in relation to the fragmentation properties during high-mass star formation. Methods. Observing the large parsec-scale Stokes I millimeter dust continuum emission with the IRAM 30 m telescope and the intermediate-scale (<0.1 pc) polarized submillimeter dust emission with the Submillimeter Array toward a sample of 20 high-mass star-forming regions allows us to quantify the dependence of the fragmentation behavior of these regions on the density and magnetic field structures. Results. Based on the IRAM 30 m data, we infer density distributions n ∝ r−p of the regions with typical power-law slopes p around ~1.5. There is no obvious correlation between the power-law slopes of the density structures on larger clump scales (~1 pc) and the number of fragments on smaller core scales (<0.1 pc). Comparing the large-scale single-dish density profiles to those derived earlier from interferometric observations at smaller spatial scales, we find that the smaller-scale power-law slopes are steeper, typically around ~2.0. The flattening toward larger scales is consistent with the star-forming regions being embedded in larger cloud structures that do not decrease in density away from a particular core. The magnetic fields of several regions appear to be aligned with filamentary structures that lead toward the densest central cores. Furthermore, we find different polarization structures; some regions exhibit central polarization holes, whereas other regions show polarized emission also toward the central peak positions. Nevertheless, the polarized intensities are inversely related to the Stokes I intensities, following roughly a power-law slope of ∝ SI−0.62. We estimate magnetic field strengths between ~0.2 and ~4.5 mG, and we find no clear correlation between magnetic field strength and the fragmentation level of the regions. A comparison of the turbulent to magnetic energies shows that they are of roughly equal importance in this sample. The mass-to-flux ratios range between ~2 and ~7, consistent with collapsing star-forming regions. Conclusions. Finding no clear correlations between the present-day large-scale density structure, the magnetic field strength, and the smaller-scale fragmentation properties of the regions, indicates that the fragmentation of high-mass star-forming regions may not be affected strongly by the initial density profiles and magnetic field properties. However, considering the limited evolutionary range and spatial scales of the presented CORE analysis, future research directions should include density structure analysis of younger regions that better resemble the initial conditions, as well as connecting the observed intermediate-scale magnetic field structure with the larger-scale magnetic fields of the parental molecular clouds. [ABSTRACT FROM AUTHOR]

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. Most massive stars are located in multiple stellar systems. The modeling of disk fragmentation, a mechanism that may plausibly lead to stellar multiplicity, relies on parallel 3D simulation codes whose agreement remains to be evaluated. Aims. Cartesian adaptive-mesh refinement (AMR) and spherical codes have frequently been used in the past decade to study massive star formation. We aim to study how the details of collapse and disk fragmentation depend on these codes. Methods. Using the Cartesian AMR code RAMSES within its self-gravity radiation-hydrodynamical framework, we compared disk fragmentation in a centrally condensed protostellar system to the findings of earlier studies performed on a grid in spherical coordinates using PLUTO. Results. To perform the code comparison, two RAMSES runs were considered, effectively giving qualitatively distinct pictures. On the one hand, when allowing for unlimited sink particle creation with no initial sink, Toomre instability and subsequent gas fragmentation leads to a multiple stellar system whose multiplicity is affected by the grid when triggering fragmentation and via numerically assisted mergers. On the other hand, using a unique, central, fixed-sink particle, a centrally-condensed system forms that is similar to that reported by PLUTO. Hence, the RAMSES-PLUTO comparison was performed with the latter and an agreement between the two codes is found as to the first rotationally supported disk formation, the presence of an accretion shock onto it, and the first fragmentation phase. Gaseous fragments form. The properties of the fragments (i.e., number, mass, and temperature) are dictated by local thermodynamics and are in agreement between the two codes given that the system has entered a highly nonlinear phase. Over the simulations, the stellar accretion rate is made of accretion bursts and continuous accretion on the same order of magnitude. As a minor difference between both codes, the dynamics of the fragments causes the disk structure to be sub-Keplerian in RAMSES, whereas it is found to be Keplerian, thus reaching quiescence, in PLUTO. We attribute this discrepancy to the central star being twice less massive in RAMSES because of the different stellar accretion subgrid models in use - rather than the potential grid effects. Conclusions. In a centrally condensed system, the agreement between RAMSES and PLUTO regarding many of the collapse properties and fragmentation process is good. In contrast, fragmentation occurring in the innermost region and given specific numerical choices (use of sink particles, grid, etc.) have a crucial impact when similar but smooth initial conditions are employed. These aspects prove more crucial than the choice of code, with regard to the system being multiple or centrally condensed. [ABSTRACT FROM AUTHOR]