Speaker
Description
In this talk we revisit the problem of the stability of circumbinary planetary orbits and how to identify stable and unstable motion in such cases. In the first part of the talk, we discuss some past results and how the problem has been dealt with so far. We present some stability criteria along with their advantages and disadvantages. In the second part of the talk, we present the latest developments in the problem of circumbinary stability. In that context, we carry out more than 3x10^8 numerical simulations of planets between the size of Mercury and the lower fusion boundary (13 Jupiter masses) which revolve around the center of mass of a stellar binary over long timescales. For the first time, three dimensional and eccentric planetary orbits are considered. The results of the numerical integrations provide us with two critical borders: an outer border beyond which all planetary orbits are stable and an inner border closer to the binary below which all planetary orbits are unstable. In between the two borders, a mixture of stable and unstable planetary orbits is observed. We provide empirical expressions in the form of multidimensional, parameterized fits for the two borders that separate the three dynamical domains. Moreover, we train a machine learning model on our data set in order to have an additional tool for predicting stable and unstable motion. Both the empirical fits and the machine learning model are tested against randomly generated circumbinary systems. The empirical formulae are also applied to the Kepler and TESS circumbinary systems, confirming the stability of the planets in these systems. Finally, the empirical fits are compared against previously derived stability criteria.