Dynamical Stability of Earth-Like Planetary Orbits in Binary Systems

This paper explores the stability of an Earth-like planet orbiting a solar mass star in the presence of an outer-lying intermediate mass companion. The overall goal is to estimate the fraction of binary systems that allow Earth-like planets to remain stable over long time scales. We numerically determine the planet's ejection time $\tauej$ over a range of companion masses ($M_C$ = 0.001 -- 0.5 $M_\odot$), orbital eccentricities $\epsilon$, and semi-major axes $a$. This suite of $\sim40,000$ numerical experiments suggests that the most important variables are the companion's mass $M_C$ and periastron distance $\rmin$ = $a(1-\epsilon)$ to the primary star. At fixed $M_C$, the ejection time is a steeply increasing function of $\rmin$ over the range of parameter space considered here (although the ejection time has a distribution of values for a given $\rmin$). Most of the integration times are limited to 10 Myr, but a small set of integrations extend to 500 Myr. For each companion mass, we find fitting formulae that approximate the mean ejection time as a function of $\rmin$. These functions can then be extrapolated to longer time scales. By combining the numerically determined ejection times with the observed distributions of orbital parameters for binary systems, we estimate that (at least) 50 percent of binaries allow an Earth-like planet to remain stable over the 4.6 Gyr age of our solar system.