Planet-planet scattering in planetesimal disks

We study the final architecture of planetary systems that evolve under the combined effects of planet-planet and planetesimal scattering. Using N-body simulations we investigate the dynamics of marginally unstable systems of gas and ice giants both in isolation and when the planets form interior to a planetesimal belt. The unstable isolated systems evolve under planet-planet scattering to yield an eccentricity distribution that matches that observed for extrasolar planets. When planetesimals are included the outcome depends upon the total mass of the planets. For system masses exceeding about one Jupiter mass the final eccentricity distribution remains broad, whereas for lower mass planetary systems a combination of divergent orbital evolution and recircularization of scattered planets results in a preponderance of nearly circular final orbits. We also study the fate of marginally stable multiple planet systems in the presence of planetesimal disks, and find that for high planet masses the majority of such systems evolve into resonance. A significant fraction lead to resonant chains that are planetary analogs of Jupiter’s Galilean satellites. We predict that a transition from eccentric to near-circular orbits will be observed once extrasolar planet surveys detect sub-Jovian mass planets at orbital radii of 5-10 AU.

💡 Research Summary

The paper investigates how the combined effects of planet‑planet scattering and interactions with a surrounding planetesimal disk shape the final architecture of planetary systems containing gas and ice giants. Using a large suite of N‑body integrations, the authors model marginally unstable systems of three to five giant planets both in isolation and embedded interior to a planetesimal belt. In the “isolated” runs, where only mutual gravitational perturbations act, the systems quickly become dynamically unstable, leading to close encounters, ejections, and collisions. The resulting eccentricity distribution of the surviving planets is broad (e ≈ 0.2–0.8) and matches the observed eccentricities of known extrasolar giant planets, confirming earlier work that planet‑planet scattering alone can explain the high‑e tail of the exoplanet population.

When a planetesimal disk is added, the outcome depends critically on the total planetary mass. The disk is modeled with a mass equal to roughly ten percent of the total planetary mass, extending from just beyond the outermost planet to ~30 AU, and composed of low‑mass test particles on near‑circular orbits. For systems whose combined mass exceeds about one Jupiter mass, the gravitational stirring from the planets dominates over the damping provided by the disk. Consequently, the high eccentricities generated by scattering are only modestly reduced, and the final eccentricity distribution remains broad, similar to the isolated case. In addition, the massive planets can capture each other into mean‑motion resonances (MMRs) as the disk slowly migrates them, often forming resonant chains reminiscent of the Laplace resonance among the Galilean moons.

Conversely, for lower‑mass systems (total mass < 1 MJup), the planetesimal disk exerts a strong dynamical friction on the planets. After an initial scattering episode that may raise eccentricities, the disk efficiently damps these values, while also causing divergent migration that spreads the planets apart. The net effect is a population of planets on nearly circular orbits (e < 0.1) at modest separations. The divergent migration also reduces the likelihood of resonant capture, so low‑mass systems tend to remain non‑resonant.

The authors further explore “marginally stable” configurations—systems that would remain stable in the absence of a disk. When such systems are placed in a planetesimal belt, massive planets again tend to be driven into resonances, forming long‑lived resonant chains. Low‑mass planets, however, are more likely to be scattered outward or to have their orbits circularized without entering resonance.

A key observational prediction emerges from these results. Current radial‑velocity and transit surveys are most sensitive to Jupiter‑mass planets inside ~5 AU, where the eccentricity distribution is indeed broad. As future instruments (e.g., direct imaging, astrometry, and next‑generation radial‑velocity spectrographs) become capable of detecting sub‑Jovian planets at 5–10 AU, the authors predict a clear transition: the eccentricity distribution should shift toward low values, reflecting the dominant role of planetesimal‑driven damping in lower‑mass systems. Detecting such a transition would provide strong evidence that both scattering and planetesimal disk interactions are essential components of planetary system evolution.

In summary, the study demonstrates that (1) planet‑planet scattering alone reproduces the observed high‑eccentricity tail of giant exoplanets; (2) the presence of a planetesimal disk can either preserve or erase this signature depending on the total planetary mass; (3) massive systems are prone to resonant chain formation, while low‑mass systems evolve toward circular, widely spaced orbits; and (4) an observable shift from eccentric to near‑circular orbits is expected as surveys probe lower‑mass planets at larger orbital radii. This integrated framework advances our understanding of the dynamical histories that produce the diverse architectures seen in extrasolar planetary systems.