Higher-order effects in the dynamics of hierarchical triple systems. III. Astrophysical implications of second-order and dotriacontapole terms

We study the long-term evolution of selected hierarchical triple systems in Newtonian gravity. We employ analytic equations derived in Paper II for the evolution of orbit-averaged orbital elements for both inner and outer orbits, which include two classes of contributions. One class consists of linear-order contributions, including quadrupole, octupole, hexadecapole and dotriacontapole orders, the latter scaling as ${\ensuremath{\epsilon}}^{6}$, where $\ensuremath{\epsilon}=a/A$, the ratio of the semimajor axes of the inner and outer orbits. The second class consists of contributions at second order in the fundamental perturbation parameter; they contribute at orders ${\ensuremath{\epsilon}}^{9/2}$, ${\ensuremath{\epsilon}}^{5}$, ${\ensuremath{\epsilon}}^{11/2}$ and ${\ensuremath{\epsilon}}^{6}$. For well studied triples such as star-planet systems perturbed by a low-mass third body (``hot Jupiters''), second-order and dotriacontapole (SOD) effects induce only small corrections. For stellar-mass binaries orbiting supermassive black holes, SOD corrections can suppress orbital flips that are generated by purely first-order effects. Planets orbiting binary star systems are susceptible to significant variations in the planetary semimajor axis, an effect that does not occur at first perturbative order. SOD effects in triple black hole systems can induce migrations of the eccentricity to significantly larger values than predicted by first-order perturbations, with implications for the gravitational-wave induced inspiral of the inner binary. We also show that in most cases, evolutions using our SOD equations are in better agreement with those from direct integration of the N-body equations of motion than those from first-order perturbations through hexadecapole order.