Computer simulations of planetary accretion dynamics - Sensitivity to initial conditions

A computer simulation program which generates mature planetary systems using only Newtonian physics and accretion with unit sticking efficiency is examined. The dependence of the simulation on a variety of radial and vertical density distribution laws, on the ratio of gas to dust in the solar nebula, on the total nebular mass, and on the orbital eccentricity of the accreting grains is explored. To approach physically realistic models, these parameters are altered according to various empirically- and theoretically-suggested criteria. The program is found to generate many recognizable planetary systems, including multiple-star systems with accompanying planets; single stars surrounded only by asteroids; and systems with large Jovian planets, planets similar to Pluto, or objects of asteroidal mass, as well as the usual terrestrial planets. Planetary spacing in the simulation obeys a relation of the Titius Bode variety; in general, the number of planets per system is inversely proportional to the orbital eccentricity. Some reservations about the simulation program are expressed; however, it does support the idea of abundant and diverse planetary systems throughout the galaxy.