Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userRadiation-Hydrodynamic Simulations of Collapse and Fragmentation in Massive Protostellar Cores
82
0
0.0
(
0
)
Authors
Mark R. Krumholz
Ask ChatGPT about the research
No Arabic abstract
We simulate the early stages of the evolution of turbulent, virialized, high-mass protostellar cores, with primary attention to how cores fragment, and whether they form a small or large number of protostars. Our simulations use the Orion adaptive mesh refinement code to follow the collapse from ~0.1 pc scales to ~10 AU scales, for durations that cover the main fragmentation phase, using three-dimensional gravito-radiation hydrodynamics. We find that for a wide range of initial conditions radiation feedback from accreting protostars inhibits the formation of fragments, so that the vast majority of the collapsed mass accretes onto one or a few objects. Most of the fragmentation that does occur takes place in massive, self-shielding disks. These are driven to gravitational instability by rapid accretion, producing rapid mass and angular momentum transport that allows most of the gas to accrete onto the central star rather than forming fragments. In contrast, a control run using the same initial conditions but an isothermal equation of state produces much more fragmentation, both in and out of the disk. We conclude that massive cores with observed properties are not likely to fragment into many stars, so that, at least at high masses, the core mass function probably determines the stellar initial mass function. Our results also demonstrate that simulations of massive star forming regions that do not include radiative transfer, and instead rely on a barotropic equation of state or optically thin heating and cooling curves, are likely to produce misleading results.
rate research
Read More
We present results from three-dimensional, self-gravitating, radiation-hydrodynamic simulations of low-mass protostellar outflows. We construct synthetic observations in 12CO in order to compare with observed outflows and evaluate the effects of beam resolution and outflow orientation on inferred outflow properties. To facilitate the comparison, we develop a quantitative prescription for measuring outflow opening angles. Using this prescription, we demonstrate that, in both simulations and synthetic observations, outflow opening angles broaden with time similarly to observed outflows. However, the interaction between the outflowing gas and the turbulent core envelope produces significant asymmetry between the red and blue shifted outflow lobes. We find that applying a velocity cutoff may result in outflow masses that are underestimated by a factor 5 or more, and masses derived from optically thick CO emission further underpredict the mass of the high-velocity gas by a factor of 5-10. Derived excitation temperatures indicate that outflowing gas is hotter than the ambient gas with temperature rising over time, which is in agreement with the simulation gas temperatures. However, excitation temperatures are otherwise not well correlated with the actual gas temperature.
Radiative transfer plays a major role in the process of star formation. Many simulations of gravitational collapse of a cold gas cloud followed by the formation of a protostellar core use a grey treatment of radiative transfer coupled to the hydrodynamics. However, dust opacities which dominate extinction show large variations as a function of frequency. In this paper, we used frequency-dependent radiative transfer to investigate the influence of the opacity variations on the properties of Larsons first core. We used a multigroup M1 moment model in a 1D radiation hydrodynamics code to simulate the spherically symmetric collapse of a 1 solar mass cloud core. Monochromatic dust opacities for five different temperature ranges were used to compute Planck and Rosseland means inside each frequency group. The results are very consistent with previous studies and only small differences were observed between the grey and multigroup simulations. For a same central density, the multigroup simulations tend to produce first cores with a slightly higher radius and central temperature. We also performed simulations of the collapse of a 10 and 0.1 solar mass cloud, which showed the properties of the first core to be independent of the initial cloud mass, with again no major differences between grey and multigroup models. For Larsons first collapse, where temperatures remain below 2000 K, the vast majority of the radiation energy lies in the IR regime and the system is optically thick. In this regime, the grey approximation does a good job reproducing the correct opacities, as long as there are no large opacity variations on scales much smaller than the width of the Planck function. The multigroup method is however expected to yield more important differences in the later stages of the collapse when high energy (UV and X-ray) radiation is present and matter and radiation are strongly decoupled.
(Abridged) Stars more massive than $20-30M_{odot}$ are so luminous that the radiation force on the cooler, more opaque outer layers can balance or exceed the force of gravity. These near or super-Eddington outer envelopes represent a long standing challenge for calculating the evolution of massive stars in one dimension, a situation that limits our understanding of the stellar progenitors of some of the most exciting and energetic explosions in the universe. In particular, the proximity to the Eddington limit has been the suspected cause for the variability, large mass loss rate and giant eruptions of an enigmatic class of massive stars: the luminous blue variables (LBVs). When in quiescence, LBVs are usually found on the hot ($T_{eff} approx 2 - 4 times 10^4$ K) S Dor instability strip. While in outburst, most LBVs stay on the cold S Dor instability strip with a $T_{eff} approx 9000$ K. Here we show that physically realistic three dimensional global radiation hydrodynamic simulations of radiation dominated massive stars with the largest supercomputers in the world naturally reproduce many observed properties of LBVs, specifically their location in the Hertzsprung-Russell (HR) diagram and their episodic mass loss with rates of $10^{-7}-10^{-5} M_{odot}/yr$. The helium opacity peak is found to play an important role to determine these properties, which is not realized in the traditional one dimensional models of massive stars. The simulations also predict that convection causes irregular envelope oscillations and 10-30% brightness variations on a typical timescale of a few days. The variability is more prominent in our models that are on the cool part of the S Dor instability. These calculations pave the way to a quantitative understanding of the structure, stability and the dominant mode of mass loss of massive stars.
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.
Subsequent to Paper I, the evolution and fragmentation of a rotating magnetized cloud are studied with use of three-dimensional MHD nested-grid simulations. After the isothermal runaway collapse, an adiabatic gas forms a protostellar first core at the center of the cloud. When the isothermal gas is stable for fragmentation in a contracting disk, the adiabatic core often breaks into several fragments. Conditions for fragmentation and binary formation are studied. All the cores which show fragmentations are geometrically thin, as the diameter-to-thickness ratio is larger than 3. Two patterns of fragmentation are found. (1) When a thin disk is supported by centrifugal force, the disk fragments through a ring configuration (ring fragmentation). This is realized in a fast rotating adiabatic core as Omega >0.2 tau_ff^-1, where Omega and tau_ff represent the angular rotation speed and the free-fall time of the core, respectively. (2) On the other hand, the disk is deformed to an elongated bar in the isothermal stage for a strongly magnetized or rapidly rotating cloud. The bar breaks into 2 - 4 fragments (bar fragmentation). Even if a disk is thin, the disk dominated by the magnetic force or thermal pressure is stable and forms a single compact body. In either ring or bar fragmentation mode, the fragments contract and a pair of outflows are ejected from the vicinities of the compact cores. The orbital angular momentum is larger than the spin angular momentum in the ring fragmentation. On the other hand, fragments often quickly merge in the bar fragmentation, since the orbital angular momentum is smaller than the spin angular momentum in this case. Comparison with observations is also shown.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
