Analytical Models for Secular Descents in Hierarchical Triple Systems

Abstract

Three-body systems are prevalent in nature, from planetary to stellar to supermassive black hole scales. In a hierarchical triple system, oscillations of the inner orbit’s eccentricity and inclination can be induced on secular timescales. Over many cycles, the octupole-level terms in the secular equations of motion can drive the system to extremely high eccentricities via the eccentric Kozai–Lidov (EKL) mechanism. The overall decrease in the inner orbit’s pericenter distance has potentially dramatic effects for realistic systems, such as tidal disruption events. We present an analytical approximation in the test-particle limit to describe individual stepwise increases in eccentricity of the inner orbit. A second approximation, also in the test-particle limit, is obtained by integrating the equations of motion and calibrating to numerical simulations to estimate the overall octupole-level time evolution of the eccentricity. The latter approach is then extended beyond the test particle to the general case. The three novel analytical approximations are compared to numerical solutions to show that the models accurately describe the form and timescale of the secular descent from large distances to a close-encounter distance (e.g., the Roche limit). By circumventing the need for numerical simulations to obtain the long-term behavior, these approximations can be used to readily estimate properties of close encounters and descent timescales for populations of systems. We demonstrate this by calculating rates of EKL-driven migration for Hot Jupiters in stellar binaries.