Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userHigh-resolution simulations of planetesimal formation in turbulent protoplanetary discs
109
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present high-resolution computer simulations of dust dynamics and planetesimal formation in turbulence generated by the magnetorotational instability. We show that the turbulent viscosity associated with magnetorotational turbulence in a non-stratified shearing box increases when going from 256^3 to 512^3 grid points in the presence of a weak imposed magnetic field, yielding a turbulent viscosity of $alphaapprox0.003$ at high resolution. Particles representing approximately meter-sized boulders concentrate in large-scale high-pressure regions in the simulation box. The appearance of zonal flows and particle concentration in pressure bumps is relatively similar at moderate (256^3) and high (512^3) resolution. In the moderate-resolution simulation we activate particle self-gravity at a time when there is little particle concentration, in contrast with previous simulations where particle self-gravity was activated during a concentration event. We observe that bound clumps form over the next ten orbits, with initial birth masses of a few times the dwarf planet Ceres. At high resolution we activate self-gravity during a particle concentration event, leading to a burst of planetesimal formation, with clump masses ranging from a significant fraction of to several times the mass of Ceres. We present a new domain decomposition algorithm for particle-mesh schemes. Particles are spread evenly among the processors and the local gas velocity field and assigned drag forces are exchanged between a domain-decomposed mesh and discrete blocks of particles. We obtain good load balancing on up to 4096 cores even in simulations where particles sediment to the mid-plane and concentrate in pressure bumps.
rate research
Read More
The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies (Safronov 1969). How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem (Dominik et al. 2007): boulders stick together poorly (Benz 2000), and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas (Weidenschilling 1977). Gravitational collapse of the solid component has been suggested to overcome this barrier (Safronov 1969, Goldreich & Ward 1973, Youdin & Shu 2002). Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer (Weidenschilling & Cuzzi 1993, Dominik et al. 2007), but turbulence must be present to explain observed gas accretion in protostellar discs (Hartmann 1998). Here we report the discovery of efficient gravitational collapse of boulders in locally overdense regions in the midplane. The boulders concentrate initially in transient high pressures in the turbulent gas (Johansen, Klahr, & Henning 2006), and these concentrations are augmented a further order of magnitude by a streaming instability (Youdin & Goodman 2005, Johansen, Henning, & Klahr 2006, Johansen & Youdin 2007) driven by the relative flow of gas and solids. We find that gravitationally bound clusters form with masses comparable to dwarf planets and containing a distribution of boulder sizes. Gravitational collapse happens much faster than radial drift, offering a possible path to planetesimal formation in accreting circumstellar discs.
This document contains refereed supplementary information for the paper ``Rapid planetesimal formation in turbulent circumstellar discs. It contains 15 sections (S1.1 -- S1.15) that address a number of subjects related to the main paper. We describe in detail the Poisson solver used to find the self-potential of the solid particles, including a linear and a non-linear test problem (S1.3). Dissipative collisions remove energy from the motion of the particles by collisional cooling (S1.4), an effect that allows gravitational collapse to occur in somewhat less massive discs (S1.7). A resolution study of the gravitational collapse of the boulders is presented in S1.6. We find that gravitational collapse can occur in progressively less massive discs as the grid resolution is increased, likely due to the decreased smoothing of the particle-mesh self-gravity solver with increasing resolution. In S1.10 we show that it is in good agreement with the Goldreich & Ward (1973) stability analysis to form several-hundred-km-sized bodies, when the analysis is applied to 5 AU and to regions of increased boulder column density. S11 is devoted to the measurement of random speeds and collision speeds between boulders. We find good agreement between our measurements and analytical theory for the random speeds, but the measured collision speeds are 3 times lower than expected from analytical theory. Higher resolution studies, and an improved analytical theory of collision speeds that takes into account epicyclic motion, will be needed to determine whether collision speeds have converged. In S1.12 we present models with no magnetic fields. The boulder layer still exhibits strong clumping, due to the streaming instability, if the global solids-to-gas ratio is increased by a factor 3. Gravitational collapse occurs as readily as in magnetised discs.
Large-scale vertical magnetic fields are believed to play a key role in the evolution of protoplanetary discs. Associated with non-ideal effects, such as ambipolar diffusion, they are known to launch a wind that could drive accretion in the outer part of the disc ($R> 1$ AU). They also potentially lead to self-organisation of the disc into large-scale axisymmetric structures, similar to the rings recently imaged by sub-millimetre or near-infrared instruments (ALMA and SPHERE). The aim of this paper is to investigate the mechanism behind the formation of these gaseous rings, but also to understand the dust dynamics and its emission in discs threaded by a large-scale magnetic field. To this end, we performed global magneto-hydrodynamics (MHD) axisymmetric simulations with ambipolar diffusion using a modified version of the PLUTO code. We explored different magnetisations with the midplane $beta$ parameter ranging from $10^5$ to $10^3$ and included dust grains -- treated in the fluid approximation -- ranging from $100 mu$m to 1 cm in size. We first show that the gaseous rings (associated with zonal flows) are tightly linked to the existence of MHD winds. Secondly, we find that millimetre-size dust is highly sedimented, with a typical scale height of 1 AU at $R=100$ AU for $beta=10^4$, compatible with recent ALMA observations. We also show that these grains concentrate into pressure maxima associated with zonal flows, leading to the formation of dusty rings. Using the radiative transfer code MCFOST, we computed the dust emission and make predictions on the ring-gap contrast and the spectral index that one might observe with interferometers like ALMA.
The mechanism through which meter-sized boulders grow to km-sized planetesimals in protoplanetary discs is a subject of active research, since it is critical for planet formation. To avoid spiralling into the protostar due to aerodynamic drag, objects must rapidly grow from cm-sized pebbles, which are tightly coupled to the gas, to large boulders of 1-100m in diameter. It is already well known that over-densities in the gaseous component of the disc provide potential sites for the collection of solids, and that significant density structures in the gaseous component of the disc (e.g., spiral density waves) can trap solids efficiently enough for the solid component of the disc to undergo further gravitational collapse due to their own self-gravity. In this work, we employ the PENCIL CODE to conduct local shearing sheet simulations of massive self-gravitating protoplanetary discs, to study the effect of anticyclonic transient vortices, or eddies, on the evolution of solids in these discs. We find that these types of structures are extremely efficient at concentrating small and intermediate-sized dust particles with friction times comparable to, or less than, the local orbital period of the disc. This can lead to significant over-densities in the solid component of the disc, with density enhancements comparable to, and even higher, than those within spiral density waves; increasing the rate of gravitational collapse of solids into bound structures.
We investigated whether outward Planetesimal Driven Migration (PDM) takes place or not in simulations when the self gravity of planetesimals is included. We performed $N$-body simulations of planetesimal disks with large width (0.7 - 4AU) which ranges over the ice line. The simulations consisted of two stages. The first stage simulations were carried out to see the runaway growth phase using the planetesimals of initially the same mass. The runaway growth took place both at the inner edge of the disk and at the region just outside the ice line. This result was utilized for the initial setup of the second stage simulations in which the runaway bodies just outside the ice line were replaced by the protoplanets with about the isolation mass. In the second stage simulations, the outward migration of the protoplanet was followed by the stopping of the migration due to the increase of the random velocity of the planetesimals. Due to this increase of random velocities, one of the PDM criteria derived in Minton and Levison (2014) was broken. In the current simulations, the effect of the gas disk is not considered. It is likely that the gas disk plays an important role in planetesimal driven migration, and we plan to study its effect in future papers.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
