\section{Introduction}

When the geometry of constant curvature spaces was brought to attention in the 19th century, mechanical problems in such spaces presented an opportunity for an interesting variation on familiar themes. In particular, studying the motion of point masses subject to extensions of Newton's usual (flat) inverse square force law to analogous `curved inverse square laws'  was soon undertaken, see e.g. the historical reviews in \cite{BMK2, DCS}.  Consequently, we speak of \textit{the} curved Kepler problem or \textit{the} curved $n$-body problem to refer to the mechanics of point masses in a space of constant curvature under such  nowadays customary extensions of Newton's inverse square force.

While the curved Kepler problem is quite similar to the flat Kepler problem --both being characterized by their conformance to analogues of Kepler's laws, or super-integrability-- there are striking differences between the curved 2-body problem and the usual 2-body problem. Namely, it has long been noticed that the absence of Galilean boosts prohibits straightforward reduction of a curved 2-body problem to a curved Kepler problem, and more recently it has been shown --as an application of Morales-Ramis theory in \cite{Shch}, and by numerically exhibiting complicated dynamics in \cite{BMK}-- that the curved 2-body problem, in contrast to the usual 2-body problem, is non-integrable. In this article we will apply classical averaging techniques to compare some features of certain orbits in these curved and flat 2-body dynamics.

When the two bodies are close, their motion is a small perturbation of the flat 2-body problem. Indeed, after reduction, we find their relative motion governed by a Kepler problem and additional perturbing terms which vanish as the distance between the bodies goes to zero. Consequently their motion at each instant is approximately along an \textit{osculating conic}: the conic section the two bodies would follow if these perturbing terms were ignored. The averaged, or \textit{secular}, dynamics take into account how the perturbing terms slowly deform these osculating conics.

Our main results, theorems \ref{thm:qper} and \ref{thm:per} below, describe quasi-periodic and periodic motions in the reduced curved 2-body problems (see fig.~\ref{fig:orbs}).

\begin{figure}[!htb]
\centering
\includegraphics[scale = .28]{Orbs1.png}
\caption{\footnotesize{(a) For two sufficiently close bodies on a sphere, there are motions (blue) along which the two bodies revolve nearly along conic sections which precess in the direction \textit{opposite} to the bodies (fast) motion around the conic. The angular momentum here is vertical, and the whole system rotates about the angular momentum axis. As the two bodies get closer to the equator, their rate of precession decreases and the eccentricity of their osculating conic increases until at the equator two types of (regularized) periodic collision orbits survive: one for which the two bodies bounce perpendicularly to the equator (red) and another for which they bounce along the equator. Similar results hold for two bodies in a negatively curved space, but now the precession is in the \textit{same} direction as the bodies motion around the conic (b) with vertical angular momentum and (c) for angular momentum 'to the right'. We expect some chaotic dynamics to occur near the collision orbits (remark \ref{rmrk:split}).}}\label{fig:orbs}
\end{figure}

We will begin in section \ref{sec:Kep} by recalling properties of the curved Kepler problem, with details in appendix \ref{app:DelPoi}.  As an aside, in seeking parametrizations of the curved Kepler orbits useful for averaging, we derive a 'curved Kepler equation' eq. (\ref{eq:KepEq}) that, as far as we can tell, is new.

The curved 2-body problem admits symmetries by the isometry group of the space of constant curvature. In section \ref{sec:2bdy}, we carry out a Meyer-Marsden-Weinstein (symplectic) reduction to formulate the equations of motion for the reduced curved 2-body problem (prop.~\ref{prop:red}), and define a scaling of the coordinates (prop.~\ref{prop:scal})  allowing us to focus our attention on the dynamics when the two bodies are close relative to the curvature of the space.

In section \ref{sec:SecDyn} we prove our two main results, establishing that the behaviour depicted in fig.~\ref{fig:orbs} holds for 'most' orbits when the two bodies in a given curved space are sufficiently close.

\begin{rmrk}
Our results are related to a remark (in section 4.1 of \cite{BMK}), in that determining the secular dynamics is a preliminary step towards a proof of non-integrability by 'Poincar\'e's methods'.   Upon defining curved analogues of the Delaunay and Poincar\'e coordinates, our methods are quite similar to those applied to the usual 3-body problem in \cite{FejQP}, although the setting here is simpler in the sense that we only need to average over one fast angle.
\end{rmrk}

\section{Curved Kepler problems}
\label{sec:Kep}
Newton's '$1/r$' potential may be characterized by either of the following properties:
\begin{itemize}
\item it is a (constant multiple) of the fundamental solution of  the 3-dimensional laplacian, $\Delta_{\mathbb{R}^3}$

\item the corresponding central force problem has orbits following Kepler's laws.
\end{itemize}
Generalizing these properties to spaces of constant positive (resp.~negative) curvature yields '$\cot\varphi$' (resp.~'$\coth\varphi$') potentials. We refer to \cite{Al} for a discussion of these long established potentials and their properties using central projection (fig. \ref{fig:CentProj}), presenting some details here in appendix \ref{app:DelPoi}. In particular, the curved Kepler problems are super integrable and on the open subset of initial conditions leading to bounded non-circular motions we have action-angle coordinates, $(L, \ell, G, g)$ where the energy depends only on $L$ (see prop. \ref{prop:DP}).

\begin{figure}[!htb]
\centering
\includegraphics[scale = .5]{gnomic.png}
\caption{\footnotesize{The curved Kepler dynamics --with \textit{fixed} center of attraction or 'sun' at $q_s$-- may be characterized by requiring that its \textit{unparametrized} trajectories centrally project to those of the flat Kepler problem in the tangent plane at $q_s$. Setting $\rho = |cq_s|$, the sphere has curvature $\kappa = 1/\rho^2$ and the hyperboloid (with the restricted Minkowski metric) has curvature $\kappa = - 1/\rho^2$.  We call $\varphi := \angle (q_scq)$ (resp.~$\tilde\varphi:= \angle(q_sc\tilde q)$ measured with the Minkowski inner product) the \textit{angular distance} from $q$ to $q_s$ (resp.~$\tilde q$ to $q_s$). The distance from $q$ on the sphere to $q_s$ is given by $\rho\varphi$, while the distance from $\tilde q$ to $q_s$ on the hyperboloid is given by $\rho \tilde\varphi$. We set $r:=|q_oq_s| = \rho \tan\varphi = \rho\tanh\tilde\varphi$.}}\label{fig:CentProj}
\end{figure}

In what follows, we will consider motions satisfying $\varphi = O({\mathbf e})$ for all time, where ${\mathbf e}$ is a small parameter. Evidently, for these Kepler problems, such orbits are central projections of elliptic orbits in the plane. Also, to avoid repetitive arguments and computations, we will give details for the positive curvature case and remark on the analogous results for negative curvature.

\section{Curved two body problems}\label{sec:2bdy}

Let $S^2$ be the sphere of curvature $\kappa$, with  $\|\cdot\|_\kappa$ the norm given by its constant curvature metric. Consider two point masses, $$(q_1,q_2) =:q\in (S^2\times S^2)\backslash \{|\cot \varphi| = \infty\} =: Q,$$ on this sphere of masses $m_1, m_2>0$ with $\varphi$ the angular distance between $q_1$ and $q_2$. Normalize the masses so that $m_1 + m_2 = 1$. We consider the Hamiltonian flow on $T^*Q$ of: $$F := \frac{\|p_1\|_\kappa^2}{2m_1} + \frac{\|p_2\|_\kappa^2}{2m_2} - \frac{m_1m_2}{\rho}\cot\varphi.$$

The tangential components of the forces are of equal magnitude and directed towards eachother along the arc $\stackrel{\frown}{q_1q_2}$. By letting $\vec q_j$ be the position vector of $q_j$ from the center of the sphere, the components of the \textit{angular momentum vector}: $$\vec C := m_1\vec q_1\times \dot{\vec q}_1 + m_2 \vec q_1\times \dot{\vec q}_2$$ are first integrals. They correspond to the symmetry by the diagonal action, $(q_1,q_2)\mapsto (gq_1, gq_2)$, of $\mathrm{SO}_3$ on $Q$. Since $\varphi\ne 0, \pi$ in $Q$, this $\mathrm{SO}_3$ action on each level set $\{\varphi = cst.\}\subset Q$ is free and transitive. By choosing a representative configuration (see fig.~\ref{fig:com}) in each level set of $\varphi$ we obtain a slice of this group action and have: \[
    Q\cong I\times \mathrm{SO}_3,~~I = (0,\pi)\ni \varphi,\]  such that the $\mathrm{SO}_3$ action is by left multiplication on the second factor.

\begin{figure}[!htb]
\centering
\includegraphics[scale = .5]{com.png}
\caption{\footnotesize{Our choice of representative configuration. Embed the sphere of curvature $\kappa$ in $\mathbb{R}^3$, centered at the origin and take an orthonormal basis $\hat i, \hat j, \hat k$ with $\vec C = C\hat k$. We take $q^o(\varphi):=( -\rho\sin m_2\varphi \hat j + \rho\cos m_2\varphi \hat k, \rho\sin m_1\varphi \hat j + \rho\cos m_1\varphi \hat k)$. Given $q\in Q$ with angular distance $\varphi$ between the two bodies, there is a unique element $g\in \mathrm{SO}_3$ with $q = gq^o(\varphi)$.}}\label{fig:com}
\end{figure}

\subsection{Reduced equations of motion} An application of symplectic reduction (see e.g. \cite{ACM} App. 5) gives:

\begin{prop}\label{prop:red} Fix the angular momentum $\vec C\ne 0$. The reduced 2-body dynamics on a sphere of curvature $\kappa = \frac{1}{\rho^2}$ takes place in $T^*I\times\mathcal{O}_\mu$, where $\mathcal{O}_\mu$ is a co-adjoint orbit in $\mathfrak{so}_3^*$. It is the Hamiltonian flow of:
$$F_{red} = Kep_\kappa + \kappa(\frac{\|\vec C\|^2}{2} - p_\theta^2) + O(\varphi),~~~dp_\varphi\wedge d\varphi + dp_\theta\wedge d\theta,$$ for $|p_\theta|<\|\vec C\|$, and with $Kep_\kappa = \frac{p_\varphi^2 + p_\theta^2/\sin^2\varphi}{2\rho^2m_1m_2} - \frac{m_1m_2}{\rho}\cot\varphi$, the Hamiltonian for a curved Kepler problem.
\end{prop}

\begin{proof} We first find the mass weighted metric, $\langle (\vec v_1, \vec v_2), (\vec u_1, \vec u_2)\rangle = m_1 \vec v_1\cdot \vec u_1 + m_2 \vec v_2\cdot \vec u_2$, in terms of our identification $Q\cong I\times \mathrm{SO}_3$ (fig.~\ref{fig:com}). Since the metric is invariant under left translations, it suffices to determine it at the identity (our representative configurations, $q^o$, in fig.~\ref{fig:com}).

Let $X_1, X_2, X_3$ be the infinitesimal symmetry vector fields generated by rotations about the $\hat i, \hat j, \hat k$ axes. Together with $\partial_\varphi$, they frame $Q$. We compute:
\begin{align*}
\langle \partial_\varphi, \partial_\varphi\rangle_{q^o} &= \rho^2m_1m_2 & \langle X_1, X_1\rangle_{q^o} &= m_1\|\hat i \times \vec q_1^o\|^2 + m_2\|\hat i \times \vec q_2^o\|^2 = \rho^2\\
\langle X_2, X_2\rangle_{q^o} &= \rho^2(m_1\cos^2m_2\varphi + m_2\cos^2m_2\varphi) & \langle X_3, X_3\rangle_{q^o} &= \rho^2(m_1\sin^2m_2\varphi + m_2\sin^2m_2\varphi)
\end{align*}
 \[   \langle X_2, X_3\rangle_{q^o} = \rho^2(m_1\cos m_2\varphi\sin m_2\varphi - m_2\cos m_1\varphi\sin m_1\varphi), \]
 all other inner products being zero. In terms of $Q$'s co-frame, $\{d\varphi, X^j\}$, dual to the above frame, the metric is given by inverting its $\{\partial_\varphi, X_j\}$ matrix representation:
 \begin{align*}
\langle d\varphi, d\varphi\rangle_{q^o} &= \frac{1}{\rho^2m_1m_2} & \langle X^1, X^1\rangle_{q^o} &= \frac{1}{\rho^2}\\
\langle X^2, X^2\rangle_{q^o} &= \frac{m_1\sin^2m_2\varphi + m_2\sin^2m_2\varphi}{m_1m_2\rho^2\sin^2\varphi} & \langle X^3, X^3\rangle_{q^o} &= \frac{m_1\cos^2m_2\varphi + m_2\cos^2m_2\varphi}{m_1m_2\rho^2\sin^2\varphi}
\end{align*}
 \[   \langle X^2, X^3\rangle_{q^o} = \frac{m_2\cos m_1\varphi\sin m_1\varphi - m_1\cos m_2\varphi\sin m_2\varphi }{m_1m_2\rho^2\sin^2\varphi}. \]

We identify $T^*Q\cong T^*I\times (\mathrm{SO}_3\times \mathfrak{so}_3^*)$ by \textit{right} translation: $\alpha_g\in T_g^*SO_3\mapsto (g, R_g^*\alpha_g)$. In these coordinates, the symplectic lift of $\mathrm{SO}_3$'s left action on $Q$ is given by: $g\cdot (\varphi, p_\varphi, h, \mu) = (\varphi, p_\varphi, gh, Ad_{g^{-1}}^*\mu)$ with associated moment map: $J(\varphi, p_\varphi, g, \mu) =  \mu$.  The reduced dynamics takes place on the quotient $P_\mu := J^{-1}(\mu)/G_\mu\cong T^*I\times\mathcal{O}_\mu$, realized by:  $(\varphi, p_\varphi, g, \mu)\mapsto (\varphi, p_\varphi, Ad_g^*\mu)$. The reduced symplectic form on $P_\mu$ is $dp_\varphi\wedge d\varphi + \Omega$, where $\Omega_\nu(ad_\xi^*\nu, ad_\eta^*\nu) = \nu([\xi, \eta])$ is the Kirillov-Kostant form on $\mathcal{O}_\mu$.

By left invariance of the two body Hamiltonian, $F$, we have:  $F(\varphi, p_\varphi, g, \mu) = F(\varphi, p_\varphi, e, Ad_g^*\mu)$, or letting $\nu = Ad_g^*\mu = \nu_j X^j(q^o)$, and $m = m_1m_2$, we have:  \begin{equation}\label{eq:Fred} F_{red} = \frac{p_\varphi^2}{2m\rho^2} - \frac{m}{\rho}\cot\varphi + \frac{\nu_1^2}{2\rho^2} + \nu_2^2\frac{m_1\sin^2m_2\varphi + m_2\sin^2m_1\varphi}{2m\rho^2\sin^2\varphi}\end{equation} \[ + \nu_3^2\frac{m_1\cos^2m_2\varphi + m_2\cos^2m_1\varphi}{2m\rho^2\sin^2\varphi} + \nu_2\nu_3\frac{m_2\cos m_1\varphi\sin m_1\varphi - m_1\cos m_2\varphi\sin m_2\varphi}{m\rho^2\sin^2\varphi}.\]

Now, using that $\frac{\sin^2ax}{\sin^2x}$ is an analytic function around $x=0$, with expansion $a^2 + O(x^2)$, and that the $\nu_2\nu_3$ coefficient is analytic around $\varphi = 0$ of $O(\varphi)$, and $\nu_1^2 + \nu_2^2 + \nu_3^2 = \|\vec C\|^2$, we have:

\[ F_{red} = \frac{p_\varphi^2}{2m\rho^2} + \frac{\nu_3^2}{2m\rho^2\sin^2\varphi} - \frac{m}{\rho}\cot\varphi + \kappa(\frac{\|\vec C\|^2}{2} - \nu_3^2) + O(\varphi).\]

Finally, recall that the Kirillov-Kostant form is given by $\Omega = dp_\theta \wedge d\theta$ when we parametrize $\mathcal{O}_\mu$ -- away from the poles $|p_\theta| = \|\vec C\|$ -- as $(\sqrt{\|\vec C\|^2 - p_\theta^2}\sin\theta, \sqrt{\|\vec C\|^2 - p_\theta^2}\cos\theta, p_\theta) = (\nu_1, \nu_2, \nu_3)$.
\end{proof}

\begin{rmrk}\label{rmrk:eqMass} For equal masses, the expression for $F_{red}$ may be simplified by taking $\psi = \varphi/2, p_\psi = 2p_\varphi$. One has: $F_{red} = Kep_\kappa + \frac{\tan\psi}{8\rho} + \kappa\frac{\|\vec C\|^2 - p_\theta^2}{2}(1 + (\tan\psi\cos\theta)^2)$, where $Kep_\kappa = \frac{1}{2\rho^2}(p_\psi^2 + \frac{p_\theta^2}{\sin^2\psi}) - \frac{\cot\psi}{8\rho}$.
\end{rmrk}

\begin{rmrk} In \cite{Shch}, symplectic reduction was also applied to reduce the curved 2-body problems. However, for purposes of comparison, it was not clear to us what representative configuration or slice of the group action was used, which led us to present a reduction here 'from scratch'.  Our choice of representative (fig.~\ref{fig:com}) was motivated by the mass metric being more diagonal, namely requiring: $\langle \partial_\varphi, X_j\rangle = 0$.
\end{rmrk}

Next, we introduce a small parameter, ${\mathbf e}$, to study the dynamics when $\varphi = O({\mathbf e})$:

\begin{prop}\label{prop:scal}
The dynamics on the open set $\varphi = O({\mathbf e})$ may be reparametrized as that of:
$$\hat F_{red} = Kep_{{\mathbf e}^2} +  Per,~~\hat\omega = d\hat L\wedge d\hat\ell + d\hat G\wedge d\hat g$$ where $Kep_{{\mathbf e}^2} = - \frac{m^3}{2\hat L^2} + \frac{{\mathbf e}^2 \hat L^2}{2m}$, $m = m_1m_2$, and $Per = {\mathbf e}^2\left(\frac{\hat C^2}{2} - \hat G^2 + (m_2-m_1)O({\mathbf e}) + O({\mathbf e}^2)\right)$.
\end{prop}

\begin{proof} Take Delaunay coordinates (app. \ref{prop:DP}), $(L,G,\ell, g)$,  for the $Kep_\kappa$ term of $F_{red}$.  Let $L^2 =: \rho{\mathbf e}\hat L^2, G^2 =: \rho{\mathbf e}\hat G^2, \|\vec C\|^2 =: \rho{\mathbf e} \hat C^2$, so that $\varphi = O({\mathbf e})$. The dynamics of $\hat F_{red} := \rho{\mathbf e} F_{red}$ with symplectic form $\hat\omega := d\hat L\wedge d\ell + d\hat G\wedge dg = \omega/\sqrt{\rho{\mathbf e}}$, corresponds to the time reparametrization: $(\rho{\mathbf e})^{3/2}\hat t = t$, where $\hat t$ is the new time. Note that the mean anomaly, $\hat\ell$, for $Kep_{{\mathbf e}^2}$ is the same as the mean anomaly, $\ell$, for the $Kep_{\kappa}$ term, as can be easily checked by comparing $d\ell = \partial_LKep_\kappa~dt$ and $d\hat\ell = \partial_{\hat L} Kep_{{\mathbf e}^2}~d\hat t$.
\end{proof}

\begin{rmrk} The scaling above includes the possibility of taking ${\mathbf e} = 1/\rho$, and in place of imagining the bodies as close, view their distance remaining bounded as the curvature of the space goes to zero.
\end{rmrk}

\begin{rmrk}
When the curvature is negative, eq. (\ref{eq:Fred}) has the trigonometric functions replaced by their hyperbolic counterparts, and with   $\mathfrak{so}_{2,1}^*$'s coadjoint orbits:  $\nu_1^2 + \nu_2^2 - \nu_3^2 = \|\vec C\|_{2,1}^2$ replacing the coadjoint orbits of $\mathfrak{so}_3^*$.  One arrives at $\hat F_{red} = Kep_{-{\mathbf e}^2} + {\mathbf e}^2(\frac{\|\hat C\|_{2,1}^2}{2} + \hat G^2 + (m_2-m_1)O({\mathbf e}) + O({\mathbf e}^2))$.
\end{rmrk}

\section{Secular dynamics}
\label{sec:SecDyn}

The \textit{secular} or averaged dynamics refers to the dynamics of the Hamiltonian: \[\langle F_{red}\rangle := \frac{1}{2\pi}\int_0^{2\pi} F_{red}(L,\ell, G, g)~d\ell. \] The secular dynamics is integrable, $L$ being a first integral. We will use the brackets to denote a functions averaged over $\ell$, namely $\langle f\rangle = \frac{1}{2\pi}\int_0^{2\pi} f~d\ell$. The dynamical relevance of the secular Hamiltonian is due to its appearance in a series of 'integrable normal forms' for the reduced curved 2-body dynamics:

\begin{prop} Consider the scaled Hamiltonian, $(\hat F_{red}, \hat\omega)$, of proposition \ref{prop:scal}. For each $k\in\mathbb N$, there is a symplectic change of coordinates, $(\psi^k)^*\hat\omega = \hat\omega$, with:
\[\hat F_{red}\circ\psi^k = Kep_{{\mathbf e}^2} +  \langle Per\rangle + \langle F^1\rangle + ... + \langle F^{k-1}\rangle + F^{k}\]
and each $F^j = O({\mathbf e}^{2j+3})$. When the masses are equal, $F^j = O({\mathbf e}^{2j + 4})$.
\end{prop}
\begin{proof} The iterative process to determine the sequence of $\psi^k$ and $F^k$'s is not new (see \cite{FejQP}). For the first step, let $\psi^1$ be the time 1-flow of a 'to be determined' Hamiltonian $\chi$. By Taylor's theorem:
\[ \hat F_{red}\circ \psi^1 = Kep_{{\mathbf e}^2} + Per + \{Kep_{{\mathbf e}^2}, \chi\} + \{Per, \chi\} + \int_0^1 (1-t)\frac{d^2}{dt^2}(\hat F_{red}\circ\psi^t)~dt.\]
Let $\tilde Per := Per - \langle Per\rangle$. Note that, while $Per = O({\mathbf e}^2)$, we have $\tilde Per = O({\mathbf e}^3)$. Taking $\hat n \chi := \int_0^\ell \tilde Per(\hat L, l, \hat G, g)~dl,$ where $\hat n = \partial_{\hat L}Kep_{{\mathbf e}^2}$, we have: \[ Per + \{Kep_{{\mathbf e}^2}, \chi\} = \langle Per\rangle,~~\chi = O({\mathbf e}^3),\] so that for ${\mathbf e}$ sufficiently small $\psi^1$ is defined. Finally, integrating by parts, one has \[ F^1 := \{Per, \chi\} + \int_0^1 (1-t)\frac{d^2}{dt^2}(\hat F_{red}\circ\psi^t)~dt = \int_0^1 \{\langle Per\rangle + t\tilde Per, \chi\}\circ \psi^t~dt = O({\mathbf e}^5).\]
\end{proof}

\begin{rmrk}\label{rmrk:kSec} The truncated, or \textit{$k$'th order secular system}: $Sec^k:=Kep_{{\mathbf e}^2} +  \langle Per\rangle + \langle F^1\rangle + ... + \langle F^{k-1}\rangle$ is integrable, since it does not depend on $\ell$. If the sequence $\psi^k$ converged, the curved 2-body problem would be integrable, although by \cite{Shch} we know there is no such convergence. For each fixed $k$ however, the $k$'th order secular system approximates the true motion over long time scales in a region of configurations with the 2-bodies sufficiently close.
\end{rmrk}

In appendix \ref{app:Per}, we compute an expansion of $\langle Per\rangle$ in powers of ${\mathbf e}$, eq. (\ref{eq:PerExp}). Via this expansion, one finds two unstable periodic orbits (red in fig.~\ref{fig:orbs}), whose stable and unstable manifolds coincide for the integrable secular dynamics. The complicated dynamics which would result from the splitting of these manifolds is a likely mechanism responsible for the non-integrability of the curved 2-body problem.

By ignoring the $O(\varphi)$ terms in $F_{red}$, the Keplerian orbits experience precession (see fig.~\ref{fig:orbs}): \[ \dot G = 0,~~\dot g = - 2\kappa G.\] When the curvature is positive, the conic precesses in a direction opposite to the particles (fast) motion around the conic, while when the curvature is negative the precession occurs in the same direction as the particles motion. We will check that for $\varphi$ sufficiently small, many such orbits survive the perturbation of including the $O(\varphi)$ terms: they may be continued to orbits of the reduced curved 2-body problem and then lifted to orbits of the curved 2-body problem.

\subsection{Continued orbits}

Applying a KAM theorem from \cite{Fej}, establishes certain quasi-periodic motions of the reduced problem. Introducing some notation, for $\gamma>0, \tau>1$, let \[ D_{\gamma,\tau}:= \{ v\in\mathbb{R}^2~:~|\mathbf{k}\cdot v|\ge \frac{\gamma}{(|k_1|+|k_2|)^\tau}, \forall \mathbf{k} = (k_1,k_2)\in\mathbb{Z}^2\backslash \mathbf{0}\} \] be the $(\gamma,\tau)$-\textit{Diophantine frequencies}. Take \[\mathbf{\hat B}:= \{ (\hat L, \hat G) ~:~ 0 < \hat G \le \hat L < 1\} \] and let \[(I,\theta) = (I_1,I_2,\theta_1, \theta_2) = (\hat L, \hat G, \ell, g) + O({\mathbf e}^3)\] be action-angle variables for some $k$th order secular system, with $\mathbf{B} = \{ (I_1, I_2) :  (\hat L, \hat G)\in \mathbf{\hat B}\}\subset \mathbb{R}^2$ equipped with the Lebesgue measure, $Leb$, restricted to $\mathbf{B}$.

The reduced curved 2-body problem has, for $\varphi$ sufficiently small,  a positive measure of invariant tori along each of which the motions are quasi-periodic. More precisely:

\begin{theorem}\label{thm:qper}
Fix $\tau >1, \hat\gamma >0, m>2$ and set $\gamma := {\mathbf e}^m\hat\gamma$. There exists ${\mathbf e}_0>0$ such that for each $0<{\mathbf e}<{\mathbf e}_0$, we have:

(i) A positive measure set, $d_{\gamma, \tau}\subset \mathbf{B}$,

(ii) for each $I^o\in d_{\gamma,\tau}$, there is a local symplectic change of coordinates, $\psi^o$, for which $\psi^o(I^o,  \theta)$ is an invariant torus of $\hat F_{red}$ with frequency in $D_{\gamma,\tau}$.

(iii) As ${\mathbf e}\to 0$, $Leb(d_{\gamma,\tau})\to Leb(\mathbf{B})$ and $\psi^o\to id$ (in the Whitney $C^\infty$ topology).
\end{theorem}
\begin{proof} We verify the hypotheses of \cite{Fej}.  Let $(I,\theta)$ be action-angle variables for a $k$'th order secular system (remark \ref{rmrk:kSec}) with $2k>2m-3$. Setting $s := I^o\in \mathbf{B}$, and $r := I - I^o$, we have by Taylor expansion:
\[ F_s^o := Sec^k = c_s^o + \alpha_s^o\cdot r + O(r^2) \] where $c_s^o := Sec_k(I^o), \alpha_s^o := dSec^k(I^o) = O(1, {\mathbf e}^2)$. Likewise, we take $F_s:=\hat F_{red}(r,\theta)$ around $s = I^o$. Since we have chosen $2k + 1>2m$, we have \[|F_s - F_s^o| = o(\gamma^2)\] for ${\mathbf e}$ sufficiently small. So  (\cite{Fej} Thm.~15 and remark 21) there exists a map $\mathbf{B}\ni s\mapsto \alpha_s\in \mathbb{R}^2$, with $|\alpha_s - \alpha_s^o|<<1$ in the $C^\infty$-Whitney topology such that \textit{provided} $\alpha_s\in D_{\gamma,\tau}$, we have $\hat F_{red}\circ\psi^o = c_s + \alpha_s\cdot r + O(r^2;\theta)$ for some local symplectic map $\psi^o$. It remains to check that $d_{\gamma,\tau}:=\{ s~:~\alpha_s\in D_{\gamma,\tau}\}$ is non-empty. Since $m>2$ and $\alpha_s = O(1, {\mathbf e}^2)$ we have (\cite{Fej} Cor.~29) \[ Leb(\mathbf{B}\backslash d_{\gamma,\tau}) = O({\mathbf e}^{\frac{m-2}{\mu}})\] for some constant $\mu>0$. So indeed, for ${\mathbf e}$ sufficiently small, $d_{\gamma,\tau}$ has positive measure and is non-empty.
\end{proof}

\begin{rmrk} As the reduced dynamics has two degrees of freedom, there is KAM-stability in the $\varphi<<1$ regime: the invariant tori form barriers in energy level sets.
\end{rmrk}

One may also apply the implicit function theorem to establish periodic orbits of long periods for the reduced dynamics:

\begin{theorem}\label{thm:per}
Let $m,n\in\mathbb N$. Then for $m$ sufficiently large, there exist periodic orbits of the reduced curved 2-body problem for which the bodies make $m$ revolutions about their osculating conics, while the osculating conic precesses around the origin $n$ times.
\end{theorem}
\begin{proof}
We will consider non-equal masses, a similar argument applies for equal masses. After iso-energetic reduction, the equations of motion for $\hat G, g$ are eq. (\ref{eq:PerExp}):
\[ \frac{d\hat G}{d\ell} = {\mathbf e}^3\mathfrak{m} \sin g + O({\mathbf e}^4),~~~\frac{dg}{d\ell} = -2{\mathbf e}^2\hat G + O({\mathbf e}^3),\] where $\mathfrak{m}(\hat L, \hat G)\ne 0$ for non-circular motions. Consider the 'long time return map':
\[ \overline P_{\mathbf e} (G^o, g^o) :=  \int_0^{2\pi/{\mathbf e}^2} (\frac{1}{{\mathbf e}}\frac{d\hat G}{d\ell}, ~\frac{dg}{d\ell})~d\ell,\] the integral taken over an orbit with initial condition $G^o, g^o$.  When $\overline P_{\mathbf e}(G^o, g^o) = (0, -2\pi n)$ and $1/{\mathbf e}^2 = m\in\mathbb N$, we have a periodic orbit as described in the theorem. The map $\overline P_{\mathbf e}$ is analytic in ${\mathbf e}$, with
 \[\overline P_0(G^o, g^o) = (\mathfrak{m}\sin g^o, ~ - 2G^o).\]
By the implicit function theorem, $g^o\equiv 0~\mod \pi, G^o = \pi n$, continues to solutions of $\overline P_{\mathbf e}(G^{\mathbf e}, g^{\mathbf e}) = (0, -2\pi n)$ for $0<{\mathbf e}<{\mathbf e}_o$, yielding periodic orbits of long periods for those ${\mathbf e}$ with  $ 1/{\mathbf e}^2 = m\in\mathbb N$.
\end{proof}

\begin{rmrk}
When the curvature is positive the precession is, as with the quasi-periodic motions, in the opposite direction to the particles motion around the conic, while for negative curvature in the same direction.
\end{rmrk}

\subsection{Lifted orbits}

The continued orbits are the image under a near identity symplectic transformation of certain orbits of $F_{red}$ with the $O(\varphi)$ terms neglected: \[ H_{red} :=  Kep_\kappa + \kappa(\frac{\|\vec C\|^2}{2} - G^2). \] We describe how the precessing Keplerian orbits of $H_{red}$ lift to $T^*Q$.

In the reduced space, $T^*I\times\mathcal{O}_\mu$, the orbits of $H_{red}$ are given by:
\[ (p_\varphi(t), \varphi(t), Ad_{g(t)}^*\mu)\]
where $g(t) = R_{\hat i}(\lambda)R_{\hat k}(\theta(t))$ is a rotation by $\theta(t)$ about the $\hat k$-axis followed by a rotation by $\lambda$ about the $\hat i$-axis. Moreover, $(\varphi(t), \theta(t))$ describe a precessing orbit of the curved Kepler problem having $G = \|\vec C\|\cos \lambda$ for its angular momentum.

We take $\vec C = C\hat k$ (fig.~\ref{fig:com}), so that the isotropy is by rotations about the $\hat k$-axis. The ambiguity in the reduced orbit is $Ad_{g(t)}^*Ad_{h(t)}^*\mu$ where $h(t) = R_{\hat k}(\omega(t))$. Hence the lifted orbit, on $T^*Q$, is of the form:
\[(*)~~~ (\varphi(t), p_\varphi(t),  h(t)g(t), \mu).\]

The Hamiltonian $H_{red}$ is the symplectic reduction of $H$ on $T^*Q$ given by $H(\varphi, p_\varphi, g, \mu) = H_{red}(\varphi, p_\varphi, Ad_g^*\mu)$. Requiring the lifted curve $(*) $ to satisfy the equations of motion of $H$, imposes the condition: \[\dot\omega = \kappa \|\vec C\|.\]

In summary, a precessing Keplerian orbit with angular momentum $G$ lifts to an orbit with angular momentum $\vec C = C\hat k$ on $Q$ where the two bodies follow precessing Keplerian orbits centered (we use the notion of 'center' from fig.~\ref{fig:com} here) on the colatitude, $\lambda$, satisfying $C\cos\lambda = G$. The whole system is rotated with angular speed $\kappa C$ about the $\hat k$ axis.

\begin{rmrk}
As $\lambda\to \pi/2$, the center of the system approaches the equator and the bodies move along more eccentric conic sections until at the equator we have collision orbits. The rate of precession decreases as one moves closer to the equator, where the collision orbits cease to precess (see remark \ref{rmrk:split}, for a finer description of the equatorial behaviour). When $C > L$, there is a maximal value of $\lambda$, and such orbits are constrained to bands around the equator.
\end{rmrk}

The true motions of the curved 2-body problem established in the previous section are qualitatively the same as such lifted orbits of $H_{red}$, performing small oscillations around them.