2022ApJ...939...81H

We examine the stability of hierarchical triple systems using direct N-body simulations without adopting a secular perturbation assumption. We estimate their disruption timescales in addition to the mere stable/unstable criterion, with particular attention to the mutual inclination between the inner and outer orbits. First, we improve the fit to the dynamical stability criterion by Mardling & Aarseth widely adopted in the previous literature. Especially, we find that that the stability boundary is very sensitive to the mutual inclination; coplanar retrograde triples and orthogonal triples are much more stable and unstable, respectively, than coplanar prograde triples. Next, we estimate the disruption timescales of triples satisfying the stability condition up to 10⁹ times the inner orbital period. The timescales follow the scaling predicted by Mushkin & Katz, especially at high e_out where their random walk model is most valid. We obtain an improved empirical fit to the disruption timescales, which indicates that the coplanar retrograde triples are significantly more stable than the previous prediction. We furthermore find that the dependence on the mutual inclination can be explained by the energy transfer model based on a parabolic encounter approximation. We also show that the disruption timescales of triples are highly sensitive to tiny changes of the initial parameters, reflecting the genuine chaotic nature of the dynamics of those systems.