Radial drift of particles in the solar nebula: implications for planetesimal formation

For standard cosmic abundances of heavy elements, a layer of small particles in the central plane of the solar nebula cannot attain the critical density for gravitational instability. Youdin and Shu (2002, Astrophys. J. 580, 494-505) suggest that the local surface density of solids can be enhanced by radial migration of particles due to gas drag. However, they consider only motions of individual particles. Collective motion due to turbulent stress on the particle layer acts to inhibit such enhancement and may prevent gravitational instability.     

Planetesimal formation by turbulent concentration

The formation of 1-1000 km diameter planetesimals from dust grains in a protoplanetary disk is a key step in planet formation. Conventional models for planetesimal formation involve pairwise sticking of dust grains, or the sedimentation of dust grains to a thin layer at the disk midplane followed by gravitational instability. Each of these mechanisms is likely to be frustrated if the disk is turbulent. Particles with stopping times comparable to the turnover time of the smallest eddies in a turbulent disk can become concentrated into dense clumps that may be the precursors of planetesimals. Such particles are roughly millimeter-sized for a typical protoplanetary disk. To survive to become planetesimals, clumps need to form in regions of low vorticity to avoid rotational breakup. In addition, clumps must have sufficient self gravity to avoid break up due to the ram pressure of the surrounding gas. Given these constraints, the rate of planetesimal formation can be estimated using a cascade model for the distribution of particle concentration and vorticity within eddies of various sizes in a turbulent disk. We estimate planetesimal formation rates and planetesimal diameters as a function of distance from a star for a range of protoplanetary disk parameters. For material with a solar composition, the dust-to-gas ratio is too low to allow efficient planetesimal formation, and most solid material will remain in small particles. Enhancement of the dust-to-gas ratio by 1-2 orders of magnitude, either vertically or radially, allows most solid material to be converted into planetesimals within the typical lifetime of a disk. Such dust-to-gas ratios may occur near the disk midplane as a result of vertical settling of short-lived clumps prior to clump breakup. Planetesimal formation rates are sensitive to the assumed size and rotational speed of the largest eddies in the disk, and formation rates increase substantially if the largest eddies rotate more slowly than the disk itself. Planetesimal formation becomes more efficient with increasing distance from the star unless the disk surface density profile has a slope of −1.5 or steeper as a function of distance. Planetesimal formation rates typically increase by an order-of-magnitude or more moving outward across the snow line for a solid surface density increase of a factor of 2. In all cases considered, the modal planetesimal size increases with roughly the square root of distance from the star. Typical modal diameters are 100 km and 400 km in the regions corresponding to the asteroid belt and Kuiper belt in the Solar System, respectively.     

Metallicity of the massive protoplanets around HR 8799 If formed by gravitational instability

The final composition of giant planets formed as a result of gravitational instability in the disk gas depends on their ability to capture solid material (planetesimals) during their ‘pre-collapse’ stage, when they are extended and cold, and contracting quasi-statically. The duration of the pre-collapse stage is inversely proportional roughly to the square of the planetary mass, so massive protoplanets have shorter pre-collapse timescales and therefore limited opportunity for planetesimal capture. The available accretion time for protoplanets with masses of 3, 5, 7, and 10 Jupiter masses is found to be and 5.67×10³ years, respectively. The total mass that can be captured by the protoplanets depends on the planetary mass, planetesimal size, the radial distance of the protoplanet from the parent star, and the local solid surface density. We consider three radial distances, 24, 38, and 68 AU, similar to the radial distances of the planets in the system HR 8799, and estimate the mass of heavy elements that can be accreted. We find that for the planetary masses usually adopted for the HR 8799 system, the amount of heavy elements accreted by the planets is small, leaving them with nearly stellar compositions.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: I. Effects of interaction with planetesimals and other protoplanets

We have performed N-body simulations on the stage of protoplanet formation from planetesimals, taking into account so-called “type-I migration,” and damping of orbital eccentricities and inclinations, as a result of tidal interaction with a gas disk without gap formation. One of the most serious problems in formation of terrestrial planets and jovian planet cores is that the migration time scale predicted by the linear theory is shorter than the disk lifetime (10⁶-10⁷ years). In this paper, we investigate retardation of type-I migration of a protoplanet due to a torque from a planetesimal disk in which a gap is opened up by the protoplanet, and torques from other protoplanets which are formed in inner and outer regions. In the first series of runs, we carried out N-body simulations of the planetesimal disk, which ranges from 0.9 to 1.1 AU, with a protoplanet seed in order to clarify how much retardation can be induced by the planetesimal disk and how long such retardation can last. We simulated six cases with different migration speeds. We found that in all of our simulations, a clear gap is not maintained for more than 10⁵ years in the planetesimal disk. For very fast migration, a gap cannot be created in the planetesimal disk. For migration slower than some critical speed, a gap does form. However, because of the growth of the surrounding planetesimals, gravitational perturbation of the planetesimals eventually becomes so strong that the planetesimals diffuse into the vicinity of the protoplanets, resulting in destruction of the gap. After the gap is destroyed, close encounters with the planetesimals rather accelerate the protoplanet migration. In this way, the migration cannot be retarded by the torque from the planetesimal disk, regardless of the migration speed. In the second series of runs, we simulated accretion of planetesimals in wide range of semimajor axis, 0.5 to 2-5 AU, starting with equal mass planetesimals without a protoplanet seed. Since formation of comparable-mass multiple protoplanets (“oligarchic growth”) is expected, the interactions with other protoplanets have a potential to alter the migration speed. However, inner protoplanets migrate before outer ones are formed, so that the migration and the accretion process of a runaway protoplanet are not affected by the other protoplanets placed inner and outer regions of its orbit. From the results of these two series of simulations, we conclude that the existence of planetesimals and multiple protoplanets do not affect type-I migration and therefore the migration shall proceed as the linear theory has suggested.     

The gravitational instability in the dust layer of a protoplanetary disk:: axisymmetric solutions for nonuniform dust density distributions in the direction vertical to the midplane

The gravitational instability in the dust layer of a protoplanetary disk with nonuniform dust density distributions in the direction vertical to the midplane is investigated. The linear analysis of the gravitational instability is performed. The following assumptions are used: (1) One fluid model is adopted, that is, difference of velocities between dust and gas are neglected. (2) The gas is incompressible. (3) Models are axisymmetric with respect to the rotation axis of the disk. Numerical results show that the critical density at the midplane is higher than the one for the uniform dust density distribution by Sekiya (1983, Prog. Theor. Phys. 69, 1116-1130). For the Gaussian dust density distribution, the critical density is 1.3 times higher, although we do not consider this dust density distribution to be realistic because of the shear instability in the dust layer. For the dust density distribution with a constant Richardson number, which is considered to be realized due to the shear instability, the critical density is 2.85 times higher and is independent of the value of the Richardson number. Further, if a constant Richardson number could decrease to the order of 0.001, the gravitational instability would be realized even for the dust to gas surface density ratio with the solar abundance. Our results give a new restriction on planetesimal formation by the gravitational instability.     

Low-speed impacts between rubble piles modeled as collections of polyhedra

We present results of modeling rubble piles as collections of polyhedra. The use of polyhedra allows more realistic (irregular) shapes and interactions (e.g. collisions), particularly for objects of different sizes. Rotational degrees of freedom are included in the modeling, which may be important components of the motion. We solved the equations of rigid-body dynamics, including frictional/inelastic collisions, for collections of up to several hundred elements. As a demonstration of the methods and to compare with previous work by other researchers, we simulated low-speed collisions between km-scale bodies with the same general parameters as those simulated by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151]. High-speed collisions appropriate to present-day asteroid encounters require additional treatment of shock effects and fragmentation and are the subject of future work; here we study regimes appropriate to planetesimal accretion and re-accretion in the aftermath of catastrophic events. Collisions between equal-mass objects at low speeds () were simulated for both head-on and off-center collisions between rubble piles made of a power-law mass spectrum of sub-elements. Very low-speed head-on collisions produce single objects from the coalescence of the impactors. For slightly higher speeds, extensive disruption occurs, but re-accretion produces a single object with most of the total mass. For increasingly higher speeds, the re-accreted object has smaller mass, finally resulting in complete catastrophic disruption with all sub-elements on escape trajectories and only small amounts of mass in re-accreted bodies. Off-center collisions at moderately low speeds produce two re-accreted objects of approximately equal mass, separating at greater than escape speed. At high speed, complete disruption occurs as with the high-speed head-on collisions. Head-on collisions at low to moderate speeds result in objects of mostly oblate shape, while higher speed collisions produce mostly prolate objects, as do off-center collisions at moderate and high speeds. Collisions carried out with the same dissipative coefficients (coefficient of restitution ?n=0.8, zero friction) as used by Leinhardt et al. [Leinhardt, Z.M., Richardson, D.C., Quinn, T., 2000. Icarus 146, 133-151] result in a value for specific energy for disruption , somewhat lower than the value of 2 J/kg found by them, while collisions with a lower coefficient of restitution and friction [?n=0.5, ?t=0, μ=0.5, similar to those used by Michel, et al. [Michel, P., Benz, W., Richardson, D.C., 2004. Planet. Space Sci. 52, 1109-1117] for SPH + N-body calculations] yield .     

A gas-poor planetesimal capture model for the formation of giant planet satellite systems

Assuming that an unknown mechanism (e.g., gas turbulence) removes most of the subnebula gas disk in a timescale shorter than that for satellite formation, we develop a model for the formation of regular (and possibly at least some of the irregular) satellites around giant planets in a gas-poor environment. In this model, which follows along the lines of the work of Safronov et al. [1986. Satellites. Univ. of Arizona Press, Tucson, pp. 89-116], heliocentric planetesimals collide within the planet's Hill sphere and generate a circumplanetary disk of prograde and retrograde satellitesimals extending as far out as ∼RH/2. At first, the net angular momentum of this proto-satellite swarm is small, and collisions among satellitesimals leads to loss of mass from the outer disk, and delivers mass to the inner disk (where regular satellites form) in a timescale ?¹⁰⁵ years. This mass loss may be offset by continued collisional capture of sufficiently small <1 km interlopers resulting from the disruption of planetesimals in the feeding zone of the giant planet. As the planet's feeding zone is cleared in a timescale ?¹⁰⁵ years, enough angular momentum may be delivered to the proto-satellite swarm to account for the angular momentum of the regular satellites of Jupiter and Saturn. This feeding timescale is also roughly consistent with the independent constraint that the Galilean satellites formed in a timescale of ¹⁰⁵-¹⁰⁶ years, which may be long enough to accommodate Callisto's partially differentiated state [Anderson et al., 1998. Science 280, 1573; Anderson et al., 2001. Icarus 153, 157-161]. In turn, this formation timescale can be used to provide plausible constraints on the surface density of solids in the satellitesimal disk (excluding satellite embryos for satellitesimals of size ∼1 km), which yields a total disk mass smaller than the mass of the regular satellites, and means that the satellites must form in several ∼10 collisional cycles. However, much more work will need to be conducted concerning the collisional evolution both of the circumplanetary satellitesimals and of the heliocentric planetesimals following giant planet formation before one can assess the significance of this agreement. Furthermore, for enough mass to be delivered to form the regular satellites in the required timescale one may need to rely on (unproven) mechanisms to replenish the feeding zone of the giant planet. We compare this model to the solids-enhanced minimum mass (SEMM) model of Mosqueira and Estrada [2003a. Icarus 163, 198-231; 2003b. Icarus 163, 232-255], and discuss its main consequences for Cassini observations of the saturnian satellite system.     

Considerations on the magnitude distributions of the Kuiper belt and of the Jupiter Trojans

By examining the absolute magnitude (H) distributions (hereafter HD) of the cold and hot populations in the Kuiper belt and of the Trojans of Jupiter, we find evidence that the Trojans have been captured from the outer part of the primordial trans-neptunian planetesimal disk. We develop a sketch model of the HDs in the inner and outer parts of the disk that is consistent with the observed distributions and with the dynamical evolution scenario known as the ‘Nice model’. This leads us to predict that the HD of the hot population should have the same slope of the HD of the cold population for 6.5<H<9, both as steep as the slope of the Trojans' HD. Current data partially support this prediction, but future observations are needed to clarify this issue. Because the HD of the Trojans rolls over at H∼9 to a collisional equilibrium slope that should have been acquired when the Trojans were still embedded in the primordial trans-neptunian disk, our model implies that the same roll-over should characterize the HDs of the Kuiper belt populations, in agreement with the results of Bernstein et al. [Bernstein, G.M., and 5 colleagues, 2004. Astron. J. 128, 1364-1390] and Fuentes and Holman [Fuentes, C.I., Holman, M.J., 2008. Astron. J. 136, 83-97]. Finally, we show that the constraint on the total mass of the primordial trans-neptunian disk imposed by the Nice model implies that it is unlikely that the cold population formed beyond 35 AU.     

Asteroids were born big

How big were the first planetesimals? We attempt to answer this question by conducting coagulation simulations in which the planetesimals grow by mutual collisions and form larger bodies and planetary embryos. The size frequency distribution (SFD) of the initial planetesimals is considered a free parameter in these simulations, and we search for the one that produces at the end objects with a SFD that is consistent with Asteroid belt constraints. We find that, if the initial planetesimals were small (e.g. km-sized), the final SFD fails to fulfill these constraints. In particular, reproducing the bump observed at diameter in the current SFD of the asteroids requires that the minimal size of the initial planetesimals was also ∼100 km. This supports the idea that planetesimals formed big, namely that the size of solids in the proto-planetary disk “jumped” from sub-meter scale to multi-kilometer scale, without passing through intermediate values. Moreover, we find evidence that the initial planetesimals had to have sizes ranging from 100 to several 100 km, probably even 1000 km, and that their SFD had to have a slope over this interval that was similar to the one characterizing the current asteroids in the same size range. This result sets a new constraint on planetesimal formation models and opens new perspectives for the investigation of the collisional evolution in the Asteroid and Kuiper belts as well as of the accretion of the cores of the giant planets.     

Ejecta from impacts at 0.2-2.3 m/s in low gravity

Collisions between planetary ring particles and in some protoplanetary disk environments occur at speeds below 10 m/s. The particles involved in these low-velocity collisions have negligible gravity and may be made of or coated with smaller dust grains and aggregates. We undertook microgravity impact experiments to better understand the dissipation of energy and production of ejecta in these collisions. Here we report the results of impact experiments of solid projectiles into beds of granular material at impact velocities from 0.2 to 2.3 m/s performed under near-weightless conditions on the NASA KC-135 Weightless Wonder V. Impactors of various densities and radii of 1 and 2 cm were launched into targets of quartz sand, JSC-1 lunar regolith simulant, and JSC-Mars-1 martian regolith simulant. Most impacts were at normal or near-normal incidence angles, though some impacts were at oblique angles. Oblique impacts led to much higher ejection velocities and ejecta masses than normal impacts. For normal incidence impacts, characteristic ejecta velocities increase with impactor kinetic energy, KE, as approximately KE^0.5. Ejecta masses could not be measured accurately due to the nature of the experiment, but qualitatively also increased with impactor kinetic energy. Some experiments were near the threshold velocity of 0.2 m/s identified in previous microgravity impact experiments as the minimum velocity needed to produce ejecta [Colwell, J.E., 2003. Icarus 164, 188-196], and the experimental scatter is large at these low speeds in the airplane experiment. A more precise exploration of the transition from low-ejecta-mass impacts to high-ejecta-mass impacts requires a longer and smoother period of reduced gravity. Coefficient of restitution measurements are not possible due to the varying acceleration of the airplane throughout the experiment.