Orbital evolution of potentially habitable planets of tidally interacting binary stars

¹ Institut für Theoretische Astrophysik (ITA), Zentrum für Astronomie (ZAH), Universität Heidelberg, Albert-Ueberle-Str. 2, 69120 Heidelberg, Germany
e-mail: dglezg7@hotmail.com
² Astronomy Department, University of Washington, Box 951580, Seattle, WA 98195, USA
³ NASA Virtual Planetary Laboratory, Seattle, WA 98195, USA

Received: 15 July 2020
Accepted: 2 November 2020

Abstract

We simulate the coupled stellar and tidal evolution of short-period binary stars (orbital period P_orb ≲ 8 days) to investigate the orbital oscillations, instellation cycles, and orbital stability of circumbinary planets (CBPs). We consider two tidal models and show that both predict an outward-then-inward evolution of the binary’s semi-major axis a_bin and eccentricity e_bin. This orbital evolution drives a similar evolution of the minimum CBP semi-major axis for orbital stability. By expanding on previous models to include the evolution of the mass concentration, we show that the maximum in the CBP orbital stability limit tends to occur 100 Myr after the planets form, a factor of 100 longer than previous investigations. This result provides further support for the hypothesis that the early stellar-tidal evolution of binary stars has removed CBPs from short-period binaries. We then apply the models to Kepler-47 b, a CBP orbiting close to its host stars’ stability limit, to show that if the binary’s initial e_bin ≳ 0.24, the planet would have been orbiting within the instability zone in the past and probably wouldn’t have survived. For stable, hypothetical cases in which the stability limit does not reach a planet’s orbit, we find that the amplitudes of a_bin and e_bin oscillations can damp by up to 10% and 50%, respectively. Finally, we consider equal-mass stars with P_orb = 7.5 days and compare the HZ to the stability limit. We find that for stellar masses ≲0.12 M_⊙, the HZ is completely unstable, even if the binary orbit is circular. For e_bin ≲ 0.5, that limit increases to 0.17 M_⊙, and the HZ is partially destabilized for stellar masses up to 0.45 M_⊙. These results may help guide searches for potentially habitable CBPs, as well as characterize their evolution and likelihood to support life after they are found.