N-Body Simulation of Planetesimal Formation through Gravitational Instability of a Dust Layer in Laminar Gas Disk (original) (raw)
Related papers
N ‐Body Simulation of Planetesimal Formation through Gravitational Instability of a Dust Layer
The Astrophysical Journal, 2007
We performed N-body simulations of a dust layer without a gas component and examined the formation process of planetesimals. We found that the formation process of planetesimals can be divided into three stages: the formation of non-axisymmetric wake-like structures, the creation of aggregates, and the collisional growth of the aggregates. Finally, a few large aggregates and many small aggregates are formed. The mass of the largest aggregate is larger than the mass predicted by the linear perturbation theory. We examined the dependence of system parameters on the planetesimal formation. We found that the mass of the largest aggregates increase as the size of the computational domain increases. However the ratio of the aggregate mass to the total mass M aggr /M total is almost constant 0.8 − 0.9. The mass of the largest aggregate increases with the optical depth and the Hill radius of particles.
2016
We performed N-body simulations of a dust layer without a gas component and examined the formation process of planetesimals. We found that the formation process of planetesimals can be divided into three stages: the formation of non-axisymmetric wake-like structures, the creation of aggregates, and the collisional growth of the aggregates. Finally, a few large aggregates and many small aggregates are formed. The mass of the largest aggregate is larger than the mass predicted by the linear perturbation theory. We examined the dependence of system parameters on the planetesimal formation. We found that the mass of the largest aggregates increase as the size of the computational domain increases. However the ratio of the aggregate mass to the total mass M aggr /M total is almost constant 0.8 − 0.9. The mass of the largest aggregate increases with the optical depth and the Hill radius of particles.
2009
The gravitational instability of a dust layer is one of the scenarios for planetesimal formation. If the density of a dust layer becomes sufficiently high as a result of the sedimentation of dust grains toward the midplane of a protoplanetary disk, the layer becomes gravitationally unstable and spontaneously fragments into planetesimals. Using a shearing box method, we performed local N-body simulations of gravitational instability of a dust layer and subsequent coagulation without gas and investigated the basic formation process of planetesimals. In this paper, we adopted the accretion model as a collision model. A gravitationally bound pair of particles is replaced by a single particle with the total mass of the pair. This accretion model enables us to perform long-term and large-scale calculations. We confirmed that the formation process of planetesimals is the same as that in the previous paper with the rubble pile models. The formation process is divided into three stages: the ...
Turbulence effects in planetesimal formation
1998
The formation of planetesimals is investigated by studying the transport of dust particles in a local three- dimensionalsimulationofaccretiondiscturbulence.Heavypar- ticlesfallrapidlytowardsthemidplane,whereaslighterparticles are strongly advected by the flow. For light particles the turbu- lence leads to a rapid redistribution of particles such that their density per unit mass is approximately constant with height. There is no pronounced concentration of particles in vortices or anticyclones, as was suggested previously. This is partly be- cause of the adverse effect of keplerian shear and also because in our simulation vortices are only short lived. However, if we assume the gas velocity to be frozen in time, there is a concen- tration of dust in ring-like structures after a few orbits. This is caused mainly by a convergence of the gas flow in those loca- tions, rather than the presence of vortices or anticyclones.
Planetesimal Formation by Gravitational Instability
The Astrophysical Journal, 2002
We investigate the formation of planetesimals via the gravitational instability of solids that have settled to the midplane of a circumstellar disk. Vertical shear between the gas and a subdisk of solids induces turbulent mixing that inhibits gravitational instability. Working in the limit of small, well-coupled particles, we find that the mixing becomes ineffective when the surface density ratio of solids to gas exceeds a critical value. Solids in excess of this precipitation limit can undergo midplane gravitational instability and form planetesimals. However, this saturation effect typically requires increasing the local ratio of solid to gaseous surface density by factors of 2-10 times cosmic abundances, depending on the exact properties of the gas disk. We discuss existing astrophysical mechanisms for augmenting the ratio of solids to gas in protoplanetary disks by such factors and investigate a particular process that depends on the radial variations of orbital drift speeds induced by gas drag. This mechanism can concentrate millimeter-sized chondrules to the supercritical surface density in few  10 6 yr, a suggestive timescale for the disappearance of dusty disks around T Tauri stars. We discuss the relevance of our results to some outstanding puzzles in planet formation theory-the size of the observed solar system and the rapid type I migration of Earth-mass bodies.
Turbulent Clustering of Protoplanetary Dust and Planetesimal Formation
The Astrophysical Journal, 2011
We study the clustering of inertial particles in turbulent flows and discuss its applications to dust particles in protoplanetary disks. Using numerical simulations, we compute the radial distribution function (RDF), which measures the probability of finding particle pairs at given distances, and the probability density function of the particle concentration. The clustering statistics depend on the Stokes number, St, defined as the ratio of the particle friction timescale, τ p , to the Kolmogorov timescale in the flow. In agreement with previous studies, we find that, in the dissipation range, the clustering intensity strongly peaks at St 1, and the RDF for St ∼ 1 shows a fast power-law increase toward small scales, suggesting that turbulent clustering may considerably enhance the particle collision rate. Clustering at inertial-range scales is of particular interest to the problem of planetesimal formation. At these large scales, the strongest clustering is from particles with τ p in the inertial range. Clustering of these particles occurs primarily around a scale where the eddy turnover time is ∼τ p. We find that particles of different sizes tend to cluster at different locations, leading to flat RDFs between different particles at small scales. In the presence of multiple particle sizes, the overall clustering strength decreases as the particle size distribution broadens. We discuss particle clustering in two recent models for planetesimal formation. We argue that, in the model based on turbulent clustering of chondrule-size particles, the probability of finding strong clusters that can seed planetesimals may have been significantly overestimated. We discuss various clustering mechanisms in simulations of planetesimal formation by gravitational collapse of dense clumps of meter-size particles, in particular the contribution from turbulent clustering due to the limited numerical resolution.
Gravoturbulent Formation of Planetesimals
The Astrophysical Journal, 2006
We explore the effect of magnetorotational turbulence on the dynamics and concentrations of boulders in local box simulations of a sub-Keplerian protoplanetary disc. The solids are treated as particles each with an independent space coordinate and velocity. We find that the turbulence has two effects on the solids. 1) Meter and decameter bodies are strongly concentrated, locally up to a factor 100 times the average dust density, whereas decimeter bodies only experience a moderate density increase. The concentrations are located in large scale radial gas density enhancements that arise from a combination of turbulence and shear. 2) For meter-sized boulders, the concentrations cause the average radial drift speed to be reduced by 40%. We find that the densest clumps of solids are gravitationally unstable under physically reasonable values for the gas column density and for the dust-to-gas ratio due to sedimentation. We speculate that planetesimals can form in a dust layer that is not in itself dense enough to undergo gravitational fragmentation, and that fragmentation happens in turbulent density fluctuations in this sublayer. Subject headings: instabilities -MHD -planetary systems: formation -planetary systems: protoplanetary disks -turbulence 1 The code is available at
Planetesimal formation by the gravitational instability of dust ring structures
arXiv (Cornell University), 2022
We investigate the gravitational instability (GI) of dust-ring structures and the formation of planetesimals by their gravitational collapse. The normalized dispersion relation of a self-gravitating ring structure includes two parameters that are related to its width and line mass (the mass per unit length). We survey these parameters and calculate the growth rate and wavenumber. Additionally, we investigate the planetesimal formation by growth of the GI of the ring that is formed by the growth of the secular GI of the protoplanetary disk. We adopt a massive, dust rich disk as a disk model. We find the range of radii for the fragmentation by the ring GI as a function of the width of the ring. The innermost radius for the ring GI is smaller for the smaller ring width. We also determine the range of the initial planetesimal mass resulting from the fragmentation of the ring GI. Our results indicate that the planetesimal mass can be as large as 10 28 g at its birth after the fragmentation. It can be as low as about 10 25 g if the ring width is 0.1% of the ring radius and the lower limit increases with the ring width. Furthermore, we obtain approximate formulas for the upper and lower limits of the planetesimal mass. We predict that the planetesimals formed by the ring GI have prograde rotations because of the Coriolis force acting on the contracting dust. This is consistent with the fact that many trans-Neptunian binaries exhibit prograde rotation.
Rapid planetesimal formation in turbulent circumstellar disks
Nature, 2007
The initial stages of planet formation in circumstellar gas discs proceed via dust grains that collide and build up larger and larger bodies 1 . How this process continues from metre-sized boulders to kilometre-scale planetesimals is a major unsolved problem 2 : boulders stick together poorly 3 , and spiral into the protostar in a few hundred orbits due to a head wind from the slower rotating gas 4 . Gravitational collapse of the solid component has been suggested to overcome this barrier 1, 5, 6 . Even low levels of turbulence, however, inhibit sedimentation of solids to a sufficiently dense midplane layer 2, 7 , but turbulence must be present to explain observed gas accretion in protostellar discs 8 . 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 9 , and these concentra-1 arXiv:0708.3890v1 [astro-ph]