scholarly journals Gravitoviscous protoplanetary disks with a dust component

2019 ◽  
Vol 631 ◽  
pp. A1 ◽  
Author(s):  
Eduard I. Vorobyov ◽  
Vardan G. Elbakyan

Aims. Spatial distribution and growth of dust in a clumpy protoplanetary disk subject to vigorous gravitational instability and fragmentation is studied numerically with sub-au resolution using the FEOSAD code. Methods. Hydrodynamics equations describing the evolution of self-gravitating and viscous protoplanetary disks in the thin-disk limit were modified to include a dust component consisting of two parts: sub-micron-sized dust and grown dust with a variable maximum radius. The conversion of small to grown dust, dust growth, friction of dust with gas, and dust self-gravity were also considered. Results. We found that the disk appearance is notably time-variable with spiral arms, dusty rings, and clumps, constantly forming, evolving, and decaying. As a consequence, the total dust-to-gas mass ratio is highly non-homogeneous throughout the disk extent, showing order-of-magnitude local deviations from the canonical 1:100 value. Gravitationally bound clumps formed through gravitational fragmentation have a velocity pattern that deviates notably from the Keplerian rotation. Small dust is efficiently converted into grown dust in the clump interiors, reaching a maximum radius of several decimeters. Concurrently, grown dust drifts towards the clump center forming a massive compact central condensation (70–100 M⊕). We argue that protoplanets may form in the interiors of inward-migrating clumps before they disperse through the action of tidal torques. We foresee the formation of protoplanets at orbital distances of several tens of au with initial masses of gas and dust in the protoplanetary seed in the (0.25–1.6) MJup and (1.0–5.5) M⊕ limits, respectively. The final masses of gas and dust in the protoplanets may however be much higher due to accretion from surrounding massive metal-rich disks/envelopes. Conclusions. Dusty rings formed through tidal dispersal of inward-migrating clumps may have a connection to ring-like structures found in youngest and massive protoplanetary disks. Numerical disk models with a dust component that can follow the evolution of gravitationally bound clumps through their collapse phase to the formation of protoplanets are needed to make firm conclusions on the characteristics of planets forming through gravitational fragmentation.

2020 ◽  
Vol 637 ◽  
pp. A5 ◽  
Author(s):  
Vardan G. Elbakyan ◽  
Anders Johansen ◽  
Michiel Lambrechts ◽  
Vitaly Akimkin ◽  
Eduard I. Vorobyov

Aims. We study the dynamics and growth of dust particles in circumstellar disks of different masses that are prone to gravitational instability during the critical first Myr of their evolution. Methods. We solved the hydrodynamics equations for a self-gravitating and viscous circumstellar disk in a thin-disk limit using the FEOSAD numerical hydrodynamics code. The dust component is made up of two different components: micron-sized dust and grown dust of evolving size. For the dust component, we considered the dust coagulation, fragmentation, momentum exchange with the gas, and dust self-gravity. Results. We found that the micron-sized dust particles grow rapidly in the circumstellar disk, reaching a few cm in size in the inner 100 au of the disk during less than 100 kyr after the disk formation, provided that fragmentation velocity is 30 ms−1. Due to the accretion of micron dust particles from the surrounding envelope, which serves as a micron dust reservoir, the approximately cm-sized dust particles continue to be present in the disk for more than 900 kyr after the disk formation and maintain a dust-to-gas ratio close to 0.01. We show that a strong correlation exists between the gas and pebble fluxes in the disk. We find that radial surface density distribution of pebbles in the disk shows power-law distribution with an index similar to that of the Minimum-mass solar nebula regardless the disk mass. We also show that the gas surface density in our models agrees well with measurements of dust in protoplanetary disks of AS 209, HD 163296, and DoAr 25 systems. Conclusions. Pebbles are formed during the very early stages of protoplanetary disk evolution. They play a crucial role in the planet formation process. Our disc simulations reveal the early onset (<105 yr) of an inwards-drifting flux of pebble-sized particles that makes up approximately between one hundredth and one tenth of the gas mass flux, which appears consistent with mm-observations of discs. Such a pebble flux would allow for the formation of planetesimals by streaming instability and the early growth of embryos by pebble accretion. We conclude that unlike the more common studies of isolated steady-state protoplanetary disks, more sophisticated global numerical simulations of circumstellar disk formation and evolution, including the pebble formation from the micron dust particles, are needed for performing realistic planet formation studies.


2019 ◽  
Vol 627 ◽  
pp. A154 ◽  
Author(s):  
Eduard I. Vorobyov ◽  
Aleksandr M. Skliarevskii ◽  
Vardan G. Elbakyan ◽  
Yaroslav Pavlyuchenkov ◽  
Vitaly Akimkin ◽  
...  

Aims. The central region of a circumstellar disk is difficult to resolve in global numerical simulations of collapsing cloud cores, but its effect on the evolution of the entire disk can be significant. Methods. We used numerical hydrodynamics simulations to model the long-term evolution of self-gravitating and viscous circumstellar disks in the thin-disk limit. Simulations start from the gravitational collapse of pre-stellar cores of 0.5–1.0 M⊙ and both gaseous and dusty subsystems were considered, including a model for dust growth. The inner unresolved 1.0 au of the disk is replaced with a central smart cell (CSC), a simplified model that simulates physical processes that may occur in this region. Results. We found that the mass transport rate through the CSC has an appreciable effect on the evolution of the entire disk. Models with slow mass transport form more massive and warmer disks, and are more susceptible to gravitational instability and fragmentation, including a newly identified episodic mode of disk fragmentation in the T Tauri phase of disk evolution. Models with slow mass transport through the CSC feature episodic accretion and luminosity bursts in the early evolution, while models with fast transport are characterized by a steadily declining accretion rate with low-amplitude flickering. Dust grows to a larger, decimeter size in the slow transport models and efficiently drifts in the CSC, where it accumulates and reaches the limit where a streaming instability becomes operational. We argue that gravitational instability, together with a streaming instability likely operating in the inner disk regions, constitute two concurrent planet-forming mechanisms, which may explain the observed diversity of exoplanetary orbits. Conclusions. We conclude that sophisticated models of the inner unresolved disk regions should be used when modeling the formation and evolution of gaseous and dusty protoplanetary disks.


2019 ◽  
Vol 491 (1) ◽  
pp. L24-L28 ◽  
Author(s):  
Shigeki Inoue ◽  
Naoki Yoshida

ABSTRACT We study the gravitational instability and fragmentation of primordial protostellar discs by using high-resolution cosmological hydrodynamics simulations. We follow the formation and evolution of spiral arms in protostellar discs, examine the dynamical stability, and identify a physical mechanism of secondary protostar formation. We use linear perturbation theory based on the spiral-arm instability (SAI) analysis in our previous studies. We improve the analysis by incorporating the effects of finite thickness and shearing motion of arms, and derive the physical conditions for SAI in protostellar discs. Our analysis predicts accurately the stability and the onset of arm fragmentation that is determined by the balance between self-gravity and gas pressure plus the Coriolis force. Formation of secondary and multiple protostars in the discs is explained by the SAI, which is driven by self-gravity and thus can operate without rapid gas cooling. We can also predict the typical mass of the fragments, which is found to be in good agreement with the actual masses of secondary protostars formed in the simulation.


2019 ◽  
Vol 485 (1) ◽  
pp. L146-L150 ◽  
Author(s):  
Elad Steinberg ◽  
Eric R Coughlin ◽  
Nicholas C Stone ◽  
Brian D Metzger

ABSTRACT The tidal destruction of a star by a massive black hole, known as a tidal disruption event (TDE), is commonly modelled using the ‘frozen-in’ approximation. Under this approximation, the star maintains exact hydrostatic balance prior to entering the tidal sphere (radius rt), after which point its internal pressure and self-gravity become instantaneously negligible and the debris undergoes ballistic free fall. We present a suite of hydrodynamical simulations of TDEs with high penetration factors β ≡ rt/rp = 5−7, where rp is the pericentre of the stellar centre of mass, calculated using a Voronoi-based moving-mesh technique. We show that basic assumptions of the frozen-in model, such as the neglect of self-gravity inside rt, are violated. Indeed, roughly equal fractions of the final energy spread accumulate exiting and entering the tidal sphere, though the frozen-in prediction is correct at the order-of-magnitude level. We also show that an $\mathcal {O}(1)$ fraction of the debris mass remains transversely confined by self-gravity even for large β which has implications for the radio emission from the unbound debris and, potentially, for the circularization efficiency of the bound streams.


2019 ◽  
Vol 492 (1) ◽  
pp. 79-95 ◽  
Author(s):  
Davide Martizzi

ABSTRACT A suite of idealized, global, gravitationally unstable, star-forming galactic disc simulations with 2 pc spatial resolution, performed with the adaptive mesh refinement code ramses, is used in this paper to predict the emergent effects of supernova feedback. The simulations include a simplified prescription for the formation of single stellar populations of mass $\sim 100 \, {\rm M}_{\odot }$, radiative cooling, photoelectric heating, an external gravitational potential for a dark matter halo and an old stellar disc, self-gravity, and a novel implementation of supernova feedback. The results of these simulations show that gravitationally unstable discs can generate violent supersonic winds with mass-loading factors η ≳ 10, followed by a galactic fountain phase. These violent winds are generated by highly clustered supernovae exploding in dense environments created by gravitational instability, and they are not produced in simulation without self-gravity. The violent winds significantly perturb the vertical structure of the disc, which is later re-established during the galactic fountain phase. Gas resettles into a quasi-steady, highly turbulent disc with volume-weighted velocity dispersion $\sigma \gt 50 \, {\rm km\, s}^{-1}$. The new configuration drives weaker galactic winds with a mass-loading factor η ≤ 0.1. The whole cycle takes place in ≤10 dynamical times. Such high time variability needs to be taken into account when interpreting observations of galactic winds from starburst and post-starburst galaxies.


2018 ◽  
Vol 614 ◽  
pp. A98 ◽  
Author(s):  
Eduard I. Vorobyov ◽  
Vitaly Akimkin ◽  
Olga Stoyanovskaya ◽  
Yaroslav Pavlyuchenkov ◽  
Hauyu Baobab Liu

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.


2018 ◽  
Vol 609 ◽  
pp. A50 ◽  
Author(s):  
C. P. Dullemond ◽  
A. B. T. Penzlin

Protoplanetary disks often appear as multiple concentric rings in dust continuum emission maps and scattered light images. These features are often associated with possible young planets in these disks. Many non-planetary explanations have also been suggested, including snow lines, dead zones and secular gravitational instabilities in the dust. In this paper we suggest another potential origin. The presence of copious amounts of dust tends to strongly reduce the conductivity of the gas, thereby inhibiting the magneto-rotational instability, and thus reducing the turbulence in the disk. From viscous disk theory it is known that a disk tends to increase its surface density in regions where the viscosity (i.e. turbulence) is low. Local maxima in the gas pressure tend to attract dust through radial drift, increasing the dust content even more. We have investigated mathematically if this could potentially lead to a feedback loop in which a perturbation in the dust surface density could perturb the gas surface density, leading to increased dust drift and thus amplification of the dust perturbation and, as a consequence, the gas perturbation. We find that this is indeed possible, even for moderately small dust grain sizes, which drift less efficiently, but which are more likely to affect the gas ionization degree. We speculate that this instability could be triggered by the small dust population initially, and when the local pressure maxima are strong enough, the larger dust grains get trapped and lead to the familiar ring-like shapes. We also discuss the many uncertainties and limitations of this model.


Sign in / Sign up

Export Citation Format

Share Document