No Arabic abstract
Similar to their low-mass counterparts, massive stars likely form via the collapse of pre-stellar molecular cores. Recent observations suggest that most massive cores are subvirial (i.e., not supported by turbulence) and therefore are likely unstable to gravitational collapse. Here we perform radiation hydrodynamic simulations to follow the collapse of turbulent massive pre-stellar cores with subvirial and virialized initial conditions to explore how their dynamic state affects the formation of massive stars and core fragmentation into companion stars. We find that subvirial cores undergo rapid monolithic collapse resulting in higher accretion rates at early times as compared to the collapse of virialized cores that have the same physical properties. In contrast, we find that virialized cores undergo a slower, gradual collapse and significant turbulent fragmentation at early times resulting in numerous companion stars. In the absence of strong magnetic fields and protostellar outflows we find that the faster growth rate of massive stars that are born out of subvirial cores leads to an increase in the radiative heating of the core thereby further suppressing fragmentation at early times when turbulent fragmentation occurs for virialized cores. Regardless of initial condition, we find that the massive accretion disks that form around massive stars dominant the accretion flow onto the star at late times and eventually become gravitationally unstable and fragment to form companion stars at late times.
(Abridged) Context. Massive stars form in magnetized and turbulent environments, and are often located in stellar clusters. Their accretion mechanism, as well as the origin of their systems stellar multiplicity are poorly understood. Aims. We study the influence of both magnetic fields and turbulence on the accretion mechanism of massive protostars and their multiplicity. Methods. We present a series of four Radiation-MHD simulations of the collapse of a massive magnetized, turbulent core of 100 $M_odot$ with the AMR code Ramses, including a hybrid radiative transfer method for stellar irradiation and ambipolar diffusion. We vary the Mach and Alfvenic Mach numbers to probe sub- and superalfvenic turbulence as well as sub- and supersonic turbulence regimes. Results. Subalfvenic turbulence leads to single stellar systems while superalfvenic turbulence leads to binary formation from disk fragmentation following spiral arm collision, with mass ratios of 1.1-1.6. In those runs, infalling gas reaches the individual disks via a transient circumbinary structure. Magnetically-regulated, thermally-dominated (plasma beta $beta>1$), Keplerian disks form in all runs, with sizes 100-200 AU and masses 1-8 $M_odot$. The disks around primary and secondary sink particles share similar properties. We observe higher accretion rates onto the secondary stars than onto their primary star companion. The primary disk orientation is found to be set by the initial angular momentum carried by turbulence. Conclusions. Small (300 AU) massive protostellar disks as those frequently observed nowadays can only be reproduced so far in the presence of (moderate) magnetic fields with ambipolar diffusion, even in a turbulent medium. The interplay between magnetic fields and turbulence sets the multiplicity of stellar clusters. A plasma beta $beta>1$ is a good indicator of streamers and disks.
The formation of massive stars is a long standing problem. Although a number of theories of massive star formation exist, ideas appear to converge to a disk-mediated accretion scenario. Here we present radiative hydrodynamic simulations of a star accreting mass via a disk embedded in a torus. We use a Monte Carlo based radiation hydrodynamics code to investigate the impact that ionizing radiation has on the torus. Ionized regions in the torus midplane are found to be either gravitationally trapped or in pressure driven expansion depending on whether or not the size of the ionized region exceeds a critical radius. Trapped Hii regions in the torus plane allow accretion to progress, while expanding Hii regions disrupt the accretion torus preventing the central star from aggregating more mass, thereby setting the stars final mass. We obtain constraints for the luminosities and torus densities that lead to both scenarios.
Using the HPC ressources of the state of Baden-Wurttemberg, we modelled for the first time the luminous burst from a young massive star by accretion of material from its close environment. We found that the surroundings of young massive stars are shaped as a clumpy disk whose fragments provoke outbursts once they fall onto the protostar and concluded that similar strong luminous events observed in high-mass star forming regions may be a signature of the presence of such disks.
We present high angular resolution observations (0.5x0.3) carried out with the Submillimeter Array (SMA) toward the AFGL2591 high-mass star forming region. Our SMA images reveal a clear chemical segregation within the AFGL2591 VLA 3 hot core, where different molecular species (Type I, II and III) appear distributed in three concentric shells. This is the first time that such a chemical segregation is ever reported at linear scales <3000 AU within a hot core. While Type I species (H2S and 13CS) peak at the AFGL2591 VLA 3 protostar, Type II molecules (HC3N, OCS, SO and SO2) show a double-peaked structure circumventing the continuum peak. Type III species, represented by CH3OH, form a ring-like structure surrounding the continuum emission. The excitation temperatures of SO2, HC3N and CH3OH (185+-11 K, 150+-20 K and 124+-12 K, respectively) show a temperature gradient within the AFGL2591 VLA 3 envelope, consistent with previous observations and modeling of the source. By combining the H2S, SO2 and CH3OH images, representative of the three concentric shells, we find that the global kinematics of the molecular gas follow Keplerian-like rotation around a 40 Mo-star. The chemical segregation observed toward AFGL2591 VLA 3 is explained by the combination of molecular UV photo-dissociation and a high-temperature (~1000 K) gas-phase chemistry within the low extinction innermost region in the AFGL2591 VLA 3 hot core.
We investigate protostellar collapse of molecular cloud cores by numerical simulations, taking into account turbulence and magnetic fields. By using the adaptive mesh refinement technique, the collapse is followed over a wide dynamic range from the scale of a turbulent cloud core to that of the first core. The cloud core is lumpy in the low density region owing to the turbulence, while it has a smooth density distribution in the dense region produced by the collapse. The shape of the dense region depends mainly on the mass of the cloud core; a massive cloud core tends to be prolate while a less massive cloud core tends to be oblate. In both cases, anisotropy of the dense region increases during the isothermal collapse. The minor axis of the dense region is always oriented parallel to the local magnetic field. All the models eventually yield spherical first cores supported mainly by the thermal pressure. Most of turbulent cloud cores exhibit protostellar outflows around the first cores. These outflows are classified into two types, bipolar and spiral flows, according to the morphology of the associated magnetic field. Bipolar flow often appears in the less massive cloud core. The rotation axis of the first core is oriented parallel to the local magnetic field for bipolar flow, while the orientation of the rotation axis from the global magnetic field depends on the magnetic field strength. In spiral flow, the rotation axis is not aligned with the local magnetic field.