Disentangling the parameter space: The role of planet multiplicity in triggering dynamical instabilities on planetary systems around white dwarfs

TL;DR

This study uses N-body simulations to explore how the number of planets influences the stability of planetary systems after their host star becomes a white dwarf, revealing increased instability with more planets and varying outcomes based on planet mass.

Contribution

It provides a systematic analysis of the impact of planet multiplicity on post-main sequence system stability through controlled N-body simulations.

Findings

Higher planet multiplicity leads to increased instability during the white dwarf phase.

Systems with low-mass planets often undergo orbit crossing without losing planets.

High-mass planet systems tend to lose multiple planets quickly after stellar evolution.

Abstract

Planets orbiting intermediate and low-mass stars are in jeopardy as their stellar hosts evolve to white dwarfs (WDs) because the dynamics of the planetary system changes due to the increase of the planet:star mass ratio after stellar mass-loss. In order to understand how the planet multiplicity affects the dynamical stability of post-main sequence (MS) systems, we perform thousands of N-body simulations involving planetary multiplicity as the variable and with a controlled physical and orbital parameter space: equal-mass planets; the same orbital spacing between adjacent planet's pairs; and orbits with small eccentricities and inclinations. We evolve the host star from the MS to the WD phase following the system dynamics for 10 Gyr. We find that the fraction of dynamically active simulations on the WD phase for two-planet systems is $10.2^{+1.2}_{-1.0}$-$25.2^{+2.5}_{-2.2}$ $\%$ and increases to $33.6^{+2.3}_{-2.2}$-$74.1^{+3.7}_{-4.6}$ $\%$ for the six-planet systems, where the ranges cover different ranges of initial orbital separations. Our simulations show that the more planets the system has, the more systems become unstable when the star becomes a WD, regardless of the planet masses and range of separations. Additional results evince that simulations with low-mass planets (1, 10 $\mathrm{M_\oplus}$) lose at most two planets, have a large fraction of systems undergoing orbit crossing without planet losses, and are dynamically active for Gyr time-scales on the WD's cooling track. On the other hand, systems with high-mass planets (100, 1000 $\mathrm{M_\oplus}$) lose up to five planets, preferably by ejections, and become unstable in the first few hundred Myr after the formation of the WD.