Numerical experiment applicable to the latest stage of planet growth

Abstract

A three-dimensional numerical model was developed with the goal of studying limited dynamical problems relevant to the latest stage of planet growth in the accretion theory. A small number of large protoplanets (∼ Moon size) of different masses, moving around the Sun, are considered. The dynamical evolution and growth of the population is studied under mutual gravitational perturbations, accretion, and collisional fragmentation processes. Gravitational encounters are treated exactly by numerical integration of the N-body problem. Outcomes of collisional fragmentation are modeled according to the results of R. Greenberg et al. (1978, Icarus, 35, 1–26). In the present work, we consider 25 protoplanets with uniform mass distribution in the range 2 × 10 25−4 × 10 26 g on heliocentric orbits in the Earth zone. These bodies are initially confined to a small volume of space to permit gravitational perturbations by close approaches and collisions within a finite length of integration time. The dynamical evolution of the swarm is followed for four different sets of initial ranges in semimajor axis, eccentricity, and inclination: Δa=0.01, 0.02, 0.04, 0.08 AU; Δe= 0.005, 0.01, 0.02, 0.04; Δi=0°3, 0°6, 1°2, 2°4. Among other results, it is found that average eccentricities and inclinations evolve toward a steady state such that i ⋍ 1 2 , e ; it is also found that, whatever the initial conditions, the population evolves toward a quasi-equilibrium relative velocity distribution corresponding to a Safronov parameter value θ⋍10. Moreover, the growth process of the growing planet presents very similar behavior in the four cases considered, except for the time scale of evolution, which increases with the initial range of orbital elements. Earlier works of this kind have been presented by L.P. Cox and J.S. Lewis (1980, Icarus, 44, 706–721) and by G.N. Wetherill (1980b, In Geol. Soc. Canad. Spec. Publ., p. 20), although a number of differences exist between the three approaches.