Computation of periodic solutions of conservative systems with application to the 3-body problem

We show how to compute families of periodic solutions of conservative systems with two-point boundary value problem continuation software. The computations include detection of bifurcations and corresponding branch switching. A simple example is used to illustrate the main idea. Thereafter we compute families of periodic solutions of the circular restricted 3-body problem. We also continue the figure-8 orbit recently discovered by Chenciner and Montgomery, and numerically computed by Simo, as the mass of one of the bodies is allowed to vary. In particular, we show how the invariances (phase-shift, scaling law, and x, y, z translations and rotations) can be dealt with. Our numerical results show, among other things, that there exists a continuous path of periodic solutions from the figure-8 orbit to a periodic solution of the restricted 3-body problem.