Parking planets in circumbinary discs

¹ Institut für Astronomie und Astrophysik, Universität Tübingen, Auf der Morgenstelle 10 72076, Germany
e-mail: anna.penzlin@uni-tuebingen.de
² Astronomy Unit, School of Physics and Astronomy, Queen Mary University of London, London E1 4NS, UK

Received: 1 September 2020
Accepted: 5 December 2020

Abstract

The Kepler space mission has discovered about a dozen planets orbiting around binary stars systems. Most of these circumbinary planets lie near their instability boundaries, at about three to five binary separations. Past attempts to match these final locations through an inward migration process were only shown to be successful for the Kepler-16 system. Here, we study ten circumbinary systems and attempt to match the final parking locations and orbital parameters of the planets with a disc-driven migration scenario. We performed 2D locally isothermal hydrodynamical simulations of circumbinary discs with embedded planets and followed their migration evolution using different values for the disc viscosity and aspect ratio. We found that for the six systems with intermediate binary eccentricities (0.1 ≤ e_bin ≤ 0.21), the final planetary orbits matched the observations closely for a single set of disc parameters, specifically, a disc viscosity of α = 10⁻⁴ and an aspect ratio of H∕r ~ 0.04. For these systems the planet masses are large enough to open at least a partial gap in their discs as they approach the binary, forcing the discs to become circularised and allowing for further migration towards the binary – ultimately leading to a good agreement with the observed planetary orbital parameters. For systems with very small or large binary eccentricities, the match was not as good as the very eccentric discs and the large inner cavities in these cases prevented close-in planet migration. In test simulations with higher than observed planet masses, a better agreement was found for those systems. The good agreement for six out of the ten modelled systems, where the relative difference between observed and simulated final planet orbit is ≤10% strongly supports the idea that planet migration in the disc brought the planets to their present locations.