STABILITY OF THE OUTER PLANETS IN MULTIRESONANT CONFIGURATIONS WITH A SELF-GRAVITATING PLANETESIMAL DISK

We study the effect of a massive planetesimal disk on the dynamical stability of the outer planets in a system representing the early solar system assuming, as has been suggested recently, that these planets were initially locked in a compact and multiresonant configuration as a result of gas-driven migration in a protoplanetary disk. The planetesimal disk is represented by an ensemble of 2000 lunar mass bodies for which the gravitational interaction is calculated self-consistently using the Mercury6.5 code. Several initial multiresonant configurations and planetesimal disk models are considered. Under such conditions a strong dynamical instability, manifested as a rapid giant planet migration and planetesimal disk dispersal, develops on a timescale of less than 40 Myr in most cases. Dynamical disk heating due to the gravitational interactions among planetesimals leads to more frequent interactions between the planetesimals and the ice giants, in comparison to models in which planetesimal–planetesimal interactions are neglected. The number of particles used to represent the planetesimal disk has implications for our results, and although our studies represent the first self-consistent calculations of unstable planetesimal-driven migration, our results point toward the need for using more realistic treatments of the planetesimal disk. Finally, in the framework of our model, we discuss the possible implications of our results on the early evolution of the solar system.