The Stability of Protoplanet Systems

Abstract We investigated the stability of 10 protoplanet systems using three-dimensional N-body simulations. We found that the time scale of instability T depends strongly on the initial random velocities v (eccentricities e and inclinations i) and orbital separations Δa. For zero initial random velocities, we confirmed the result of Chambers et al. (1996, Icarus119, 261–268) that T is proportional to exp(Δa). For finite random velocities, we found that T depends strongly on the initial random velocities. The relation between T and Δa is still expressed as logT=b+cΔa. However, both b and c depend on initial random velocities and the slope, b, becomes smaller for larger v. Even for relatively small initial eccentricities such as e∼2rH/a, where rH is the Hill radius, the time scale can be reduced by a factor of 10 compared with the case of the zero random velocity. Therefore, the time scale of the formation of inner planets might be much shorter than what implied by Chambers et al.