ABSTRACT We present the latest development of the disk gravitational instability and fragmentation model, originally introduced by us to explain episodic accretion bursts in the early stages of star ...formation. Using our numerical hydrodynamics model with improved disk thermal balance and star-disk interaction, we computed the evolution of protostellar disks formed from the gravitational collapse of prestellar cores. In agreement with our previous studies, we find that cores of higher initial mass and angular momentum produce disks that are more favorable to gravitational instability and fragmentation, while a higher background irradiation and magnetic fields moderate the disk tendency to fragment. The protostellar accretion in our models is time-variable, thanks to the nonlinear interaction between different spiral modes in the gravitationally unstable disk, and can undergo episodic bursts when fragments migrate onto the star owing to the gravitational interaction with other fragments or spiral arms. Most bursts occur in the partly embedded Class I phase, with a smaller fraction taking place in the deeply embedded Class 0 phase and a few possible bursts in the optically visible Class II phase. The average burst duration and mean luminosity are found to be in good agreement with those inferred from observations of FUors. The model predicts the existence of two types of bursts: the isolated ones, showing well-defined luminosity peaks separated with prolonged periods ( yr) of quiescent accretion, and clustered ones, demonstrating several bursts occurring one after another during just a few hundred years. Finally, we estimate that 40%-70% of the star-forming cores can display bursts after forming a star-disk system.
We revisit our original papers on the burst mode of accretion by incorporating a detailed energy balance equation into a thin-disk model for the formation and evolution of circumstellar disks around ...low-mass protostars. Our model includes the effect of radiative cooling, viscous and shock heating, and heating due to stellar and background irradiation. Following the collapse from the prestellar phase allows us to model the early embedded phase of disk formation and evolution. During this time, the disk is susceptible to fragmentation, depending upon the properties of the initial prestellar core. Globally, we find that higher initial core angular momentum and mass content favors more fragmentation, but higher levels of background radiation can moderate the tendency to fragment. A higher rate of mass infall onto the disk than that onto the star is a necessary but not a sufficient condition for disk fragmentation. More locally, both the Toomre Q-parameter needs to be below a critical value and the local cooling time needs to be shorter than a few times the local dynamical time. Fragments that form during the early embedded phase tend to be driven into the inner disk regions and likely trigger mass accretion and luminosity bursts that are similar in magnitude to FU-Orionis-type or EX-Lupi-like events. Disk accretion is shown to be an intrinsically variable process, thanks to disk fragmentation, nonaxisymmetric structure, and the effect of gravitational torques. The additional effect of a generic {alpha}-type viscosity acts to reduce burst frequency and accretion variability, and is likely to not be viable for values of {alpha} significantly greater than 0.01.
Aims.
Accretion and luminosity bursts can be triggered by three distinct mechanisms: the magnetorotational instability (MRI) in the inner disk regions, clump infall in gravitationally fragmented ...disks, and close encounters with an intruder star. We study all three of these burst mechanisms to determine the disk kinematic characteristics that can help to distinguish between them.
Methods.
Numerical hydrodynamics simulations in the thin-disk limit were employed to model the bursts in disk environments that are expected for each burst mechanism.
Results.
We found that the circumstellar disks featuring accretion bursts can bear kinematic features that are distinct for different burst mechanisms, which can be useful when identifying the origin of a particular burst. The disks in the stellar encounter and clump-infall models are characterized by deviations from the Keplerian rotation of tens of per cent, while the disks in the MRI models are characterized by deviations of only a few per cent, which is mostly caused by the gravitational instability that fuels the MRI bursts. Velocity channel maps also show distinct kinks and wiggles, which are caused by gas disk flows that are particular to each considered burst mechanism. The deviations of velocity channels in the burst-hosting disks from a symmetric pattern typical of Keplerian disks are strongest for the clump-infall and collision models, and carry individual features that may be useful for the identification of the corresponding burst mechanism. The considered burst mechanisms produce a variety of light curves with the burst amplitudes varying in the Δ
m
= 2.5−3.7 limits, except for the clump-infall model where Δ
m
can reach 5.4, although the derived numbers may be affected by a small sample and boundary conditions.
Conclusions.
Burst-triggering mechanisms are associated with distinct kinematic features in the burst-hosting disks that may be used for their identification. Further studies including a wider model parameter space and the construction of synthetic disk images in thermal dust and molecular line emission are needed to constrain the mechanisms that lead to FU Orionis bursts.
Context. Aims. The long-term evolution of a circumstellar disk starting from its formation and ending in the T Tauri phase was simulated numerically with the purpose of studying the evolution of dust ...in the disk with distinct values of the viscous α-parameter and dust fragmentation velocity vfrag. Methods. We solved numerical hydrodynamics equations in the thin-disk limit, which were modified to include a dust component consisting of two parts: sub-micron-sized dust, and grown dust with a maximum radius ar. The former is strictly coupled to the gas, while the latter interacts with the gas through friction. Dust growth, dust self-gravity, and the conversion of small to grown dust were also considered. Results. We found that the process of dust growth that is known for the older protoplanetary phase also holds for the embedded phase of the disk evolution. The dust growth efficiency depends on the radial distance from the star – ar is largest in the inner disk and gradually declines with radial distance. In the inner disk, ar is limited by the dust fragmentation barrier. The process of small-to-grown dust conversion is very fast once the disk is formed. The total mass of the grown dust in the disk (beyond 1 AU) reaches tens or even hundreds of Earth masses as soon as in the embedded phase of star formation, and an even greater amount of grown dust drifts in the inner, unresolved 1 AU of the disk. Dust does not usually grow to radii greater than a few cm. A notable exception are models with α ≤ 10−3, in which case a zone with reduced mass transport develops in the inner disk and dust can grow to meter-sized boulders in the inner 10 AU. Grown dust drifts inward and accumulates in the inner disk regions. This effect is most pronounced in the α ≤ 10−3 models, where several hundreds of Earth masses can be accumulated in a narrow region of several AU from the star by the end of embedded phase. The efficiency of grown dust accumulation in spiral arms is stronger near corotation where the azimuthal velocity of dust grains is closest to the local velocity of the spiral pattern. In the framework of the adopted dust growth model, the efficiency of small-to-grown dust conversion was found to increase for lower values of α and vfrag.
Aims. Migration of dense gaseous clumps that form in young protostellar disks via gravitational fragmentation is investigated to determine the likelihood of giant planet formation. Methods. ...High-resolution numerical hydrodynamics simulations in the thin-disk limit are employed to compute the formation and long-term evolution of a gravitationally unstable protostellar disk around a solar-mass star. Results. We show that gaseous clumps that form in the outer regions of the disk (>100 au) through disk fragmentation are often perturbed by other clumps or disk structures, such as spiral arms, and migrate toward the central star on timescales from a few thousand to few tens of thousands of years. The migration timescale is slowest when stellar motion in response to the disk gravity is considered. When approaching the star, the clumps first gain mass (up to several tens of MJup), but then quickly lose most of their diffuse envelopes through tidal torques. Part of the clump envelope can be accreted onto the central star causing an FU-Orionis-type accretion and luminosity outburst. The tidal mass loss helps the clumps to significantly slow down or even halt their inward migration at a distance of a few tens of au from the protostar. The resulting clumps are heavily truncated both in mass and size compared to their wider orbit counterparts, keeping only a dense and hot nucleus. During the inward migration, the temperature in the clump interiors may exceed the molecular hydrogen dissociation limit (2000 K) and the central region of the clump can collapse into a gas giant protoplanet. Moreover, migrating clumps may experience close encounters with other clumps, resulting in the ejection of the least massive (planetary-mass) clumps from the disk. We argue that FU-Orionis-type luminosity outbursts may be the end product of disk fragmentation and clump inward migration, preceding the formation of giant protoplanets on tens of au orbits in systems such as HR 8799.
We present a calculation of protostellar disk formation and evolution in which gaseous clumps (essentially, the first Larson cores formed via disk fragmentation) are ejected from the disk during the ...early stage of evolution. This is a universal process related to the phenomenon of ejection in multiple systems of point masses. However, it occurs in our model entirely due to the interaction of compact, gravitationally bound gaseous clumps and is free from the smoothing-length uncertainty that is characteristic of models using sink particles. Clumps that survive ejection span a mass range of 0.08-0.35 M sub(middot in circle), and have ejection velocities 0.8 + or - 0.35 km s super(-1), which are several times greater than the escape speed. We suggest that, upon contraction, these clumps can form substellar or low-mass stellar objects with notable disks, or even close-separation very low mass binaries. In this hybrid scenario, allowing for ejection of clumps rather than finished protostars/proto-brown-dwarfs, disk formation and the low velocity dispersion of low-mass objects are naturally explained, while it is also consistent with the observation of isolated low-mass clumps that are ejection products. We conclude that clump ejection and the formation of isolated low-mass stellar and substellar objects is a common occurrence, with important implications for understanding the initial mass function, the brown dwarf desert, and the formation of stars in all environments and epochs.
Abstract
Previous observations have shown that the ≲10 au, ≳400 K hot inner disk of the archetypal accretion outburst young stellar object, FU Ori, is dominated by viscous heating. To constrain dust ...properties in this region, we have performed radio observations toward this disk using the Karl G. Jansky Very Large Array in 2020 June–July, September, and November. We also performed complementary optical photometric monitoring observations. We found that the dust thermal emission from the hot inner disk mid-plane of FU Ori has been approximately stationary and the maximum dust grain size is ≳1.6 mm in this region. If the hot inner disk of FU Ori, which is inward of the 150–170 K water snowline, is turbulent (e.g., corresponding to a Sunyaev & Shakura viscous
α
t
≳ 0.1), or if the actual maximum grain size is still larger than the lower limit we presently constrain, then as suggested by the recent analytical calculations and the laboratory measurements, water-ice-free dust grains may be stickier than water-ice-coated dust grains in protoplanetary disks. Additionally, we find that the free–free emission and the Johnson
B-
and
V
-band magnitudes of these binary stars were brightening in 2016–2020. The optical and radio variability might be related to the dynamically evolving protostellar- or disk-accretion activities. Our results highlight that the hot inner disks of outbursting objects are important laboratories for testing models of dust grain growth. Given the active nature of such systems, to robustly diagnose the maximum dust grain sizes, it is important to carry out coordinated multiwavelength radio observations.
Aims.
Accretion bursts triggered by the magnetorotational instability (MRI) in the innermost disk regions were studied for protoplanetary gas-dust disks that formed from prestellar cores of a various ...mass
M
core
and mass-to-magnetic flux ratio
λ
.
Methods.
Numerical magnetohydrodynamics simulations in the thin-disk limit were employed to study the long-term (~1.0 Myr) evolution of protoplanetary disks with an adaptive turbulent
α
-parameter, which explicitly depends on the strength of the magnetic field and ionization fraction in the disk. The numerical models also feature the co-evolution of gas and dust, including the back-reaction of dust on gas and dust growth.
Results.
A dead zone with a low ionization fraction of
x
≲10
−13
and temperature on the order of several hundred Kelvin forms in the inner disk soon after its formation, extending from several to several tens of astronomical units depending on the model. The dead zone features pronounced dust rings that are formed due to the concentration of grown dust particles in the local pressure maxima. Thermal ionization of alkaline metals in the dead zone trigger the MRI and associated accretion burst, which is characterized by a sharp rise, small-scale variability in the active phase, and fast decline once the inner MRI-active region is depleted of matter. The burst occurrence frequency is highest in the initial stages of disk formation and is driven by gravitational instability (GI), but it declines with diminishing disk mass-loading from the infalling envelope. There is a causal link between the initial burst activity and the strength of GI in the disk fueled by mass infall from the envelope. We find that the MRI-driven burst phenomenon occurs for
λ
= 2–10, but diminishes in models with
M
core
≲
M
⊙
, suggesting a lower limit on the stellar mass for which the MRI-triggered burst can occur.
Conclusions.
The MRI-triggered bursts occur for a wide range of mass-to-magnetic flux ratios and initial cloud core masses. The burst occurrence frequency is highest in the initial disk formation stage and reduces as the disk evolves from a gravitationally unstable to a viscous-dominated state. The MRI-triggered bursts are intrinsically connected with the dust rings in the inner disk regions, and both can be a manifestation of the same phenomenon, that is to say the formation of a dead zone.
ABSTRACT
3D simulations of high mass young stellar object (HMYSO) growth show that their circumstellar discs fragment on to multiple self-gravitating objects. Accretion of these by HMYSO may explain ...episodic accretion bursts discovered recently. We post-process results of a previous 3D simulation of a HMYSO disc with a 1D code that resolves the disc and object dynamics down to the stellar surface. We find that burst-like deposition of material into the inner disc seen in 3D simulations by itself does not always signify powerful accretion bursts. Only high density post-collapse clumps crossing the inner computational boundary may result in observable bursts. The rich physics of the inner disc has a significant impact on the expected accretion bursts: (1) in the standard turbulent viscosity discs, migrating objects can stall at a migration trap at the distance of a few au from the star. However, in discs powered by magnetized winds, the objects are able to cross the trap and produce bursts akin to those observed so far. (2) Migrating objects may interact with and modify the thermal (hydrogen ionization) instability of the inner disc, which can be responsible for longer duration and lower luminosity bursts in HMYSOs. (3) If the central star is bloated to a fraction of an au by a previous episode of high accretion rate, or if the migrating object is particularly dense, a merger rather than a disc-mediated accretion burst results; (4) Object disruption bursts may be super-Eddington, leading to episodic feedback on HMYSO surroundings via powerful outflows.
ABSTRACT
Supermassive stars (SMSs) with ∼104–105 M⊙ are candidate objects for the origin of supermassive black holes observed at redshift z > 6. They are supposed to form in primordial-gas clouds ...that provide the central stars with gas at a high accretion rate, but their growth may be terminated in the middle due to the stellar ionizing radiation if the accretion is intermittent and its quiescent periods are longer than the Kelvin–Helmholtz (KH) time-scales at the stellar surfaces. In this paper, we examine the role of the ionizing radiation feedback based on the accretion history in two possible SMS-forming clouds extracted from cosmological simulations, following their evolution with vertically integrated two-dimensional hydrodynamic simulations with detailed thermal and chemical models. The consistent treatment of the gas thermal evolution is crucial for obtaining the realistic accretion history, as we demonstrate by performing an additional run with a barotropic equation of state, in which the fluctuation of the accretion rate is artificially suppressed. We find that although the accretion becomes intermittent due to the formation of spiral arms and clumps in gravitationally unstable discs, the quiescent periods are always shorter than the KH time-scales, implying that SMSs can form without affected by the ionizing radiation.