Gravitational N-body problem on the accretion process of terrestrial planets

Numerical integration of the gravitational N-body problem has been carried out for a variety of photoplanetary clusters in the range N = 100 to 200. Particles are assumed to coagulate at collisions irrespective of relative velocity and mass ratio of the particles. It is shown graphically how the dispersed N-bodies accumulate to a single planet through mutual collisions. The velocity distribution and size distribution of bodies are also investigated as functions of time in the accretion process. The root mean square velocity of bodies in a cluster increases with time in an early stage of accretion but decreases with time in a late stage of accretion. Accretion rates of planets are found to be dependent strongly on the initial number density distribution, the initial size distribution, and the initial velocity distribution of bodies. Formation of satellites of about 10% in the planet mass is common to most cases in the present study. A substantial mass of bodies also escapes from the cluster. Many satellites and escapers formed during the accretion process of planets may be source materials of heavy bombardment in the early history of planets.