The hierarchical three-body problem at $\mathcal{O} (G^2)$

Employing techniques from scattering amplitudes and effective field theory, we model the dynamics of hierarchical triples, which are three-body systems composed of two bodies separated by a distance $r$ and a third body a distance $\rho$ away, with $r \ll \rho$. We apply the method of regions to systematically expand in the small ratio $r/\rho$ and illustrate this approach for evaluating Fourier transform integrals, which have been the bottleneck for deriving complete results in position space. In the limit where the distant third body is much heavier than the other two, we derive new analytic results in position space for the three-body conservative potential at $\mathcal{O}(G^2)$ and at leading and next-to-leading order in $r/\rho$. We also derive new results for arbitrary masses in the rest frame of the distant particle. Our results are exact in velocity, and can be used in analyses involving both bound and unbound hierarchical triples in astrophysical systems.