HR 8799 is a four planet system that also hosts a debris disk. By numerically integrating both planets and a planetesimal disk, we find interactions between an exterior planetesimal disk and the planets can influence the lifetime of the system. We first consider resonant planetary configurations that remained stable for at least 7 Myrs sans debris disk. An exterior debris disk with only ∼ 1% the mass of the outermost planet (approximately a Neptune mass) was sufficiently large enough to pull the system out of resonance after 2 to 6 Myrs. Secondly, we consider configurations which are unstable in less than a few hundred thousand years. We find that these can be stabilized by a debris disk with a mass of more than ∼ 10% that of the outermost planet. Our two sets of simulations suggest that estimates of the long term stability of a planetary system should take into account the role of the debris disk.

We describe a parallel hybrid symplectic integrator for planetary system integration that runs on a graphics processing unit (GPU). The integrator identifies close approaches between particles and switches from symplectic to Hermite algorithms for particles that require higher resolution integrations. The integrator is approximately as accurate as other hybrid symplectic integrators but is GPU accelerated.

We explore scenarios for the origin of two different density planets in the Kepler 36 system in adjacent orbits near the 7:6 mean motion resonance. We find that fine tuning is required in the stochastic forcing amplitude, the migration rate and planet eccentricities to allow two convergently migrating planets to bypass mean motion resonances such as the 4:3, 5:4 and 6:5, and yet allow capture into the 7:6 resonance. Stochastic forcing can eject the system from resonance causing a collision between the planets, unless the disk causing migration and the stochastic forcing is depleted soon after resonance capture.We explore a scenario with approximately Mars mass embryos originating exterior to the two planets and migrating inwards toward two planets. We find that gravitational interactions with embryos can nudge the system out of resonances. Numerical integrations with about a half dozen embryos can leave the two planets in the 7:6 resonance. Collisions between planets and embryos have a wide distribution of impact angles and velocities ranging from accretionary to disruptive. We find that impacts can occur at sufficiently high impact angle and velocity that the envelope of a planet could have been stripped, leaving behind a dense core. Some of our integrations show the two planets exchanging locations, allowing the outer planet that had experienced multiple collisions with embryos to become the innermost planet. A scenario involving gravitational interactions and collisions with embryos may account for both the proximity of the Kepler 36 planets and their large density contrast. 1 The mutual Hill radius r mH ≡ m b +mc 3M * a b +ac 2 .c 0000 RAS

A fraction of multiple planet candidate systems discovered from transits by the Kepler mission contain pairs of planet candidates that are in orbital resonance or are spaced slightly too far apart to be in resonance. We focus here on the four planet system, KOI 730, that has planet periods satisfying the ratios 8:6:4:3. By numerically integrating four planets initially in this resonant configuration in proximity to an initially exterior cold planetesimal disk, we find that of the order of a Mars mass of planet-orbit-crossing planetesimals is sufficient to pull this system out of resonance. Approximately one Earth mass of planet-orbit-crossing planetesimals increases the interplanetary spacings sufficiently to resemble the multiple planet candidate Kepler systems that lie just outside of resonance. This suggests that the closely spaced multiple planet Kepler systems, host only low mass debris disks or their debris disks have been extremely stable. We find that the planetary inclinations increase as a function of the mass in planetesimals that have crossed the orbits of the planets. If systems are left at zero inclination and in resonant chains after depletion of the gas disk then we would expect a correlation between distance to resonance and mutual planetary inclinations. This may make it possible to differentiate between dynamical mechanisms that account for the fraction of multiple planet systems just outside of resonance. c 0000 RAS