Search for periodic solutions with special symmetry in the Hill problem

Hill problem is a well-known celestial mechanics problem, which gives the description of dynamics of a satellite near the minor of two active masses in the uniformly rotating frame (sinodical system of coordinates). This problem is a limiting case of the well-known restricted three-body problem (RTBP), and its periodical solutions can be continued into the corresponding solutions of the RTBP.
We propose to embed the Hill problem in so-called generalized Hill problem, which Hamiltonian can be written in the form $H=H_0+\varepsilon R$, where
\begin{equation*} \begin{aligned} H_0&=\frac12\left(y_1^2+y_2^2\right)+x_2y_1-x_1y_2+\frac{\sigma}{r},\\ R&=-x_1^2+\frac12x_2^2,\quad r=\sqrt{x_1^2+x_2^2}. \end{aligned} \end{equation*}
Here $ x$ and $ y$ are vectors of canonically conjugate coordinates and momenta correspondingly, $\varepsilon\in[0;1]$ and $\sigma\in\{-1,0\}$ are some parameters.
Canonical equations of motions of the problem are invariant under the discrete group with two generators of linear transformations of the extended phase space. Therefore all periodic solutions drop into the four group of asymmetric, singly symmetric, doubly symmetric and centrally symmetric solutions.
Singly and doubly symmetric periodic solutions were intensively studied with the help of singular perturbation theory earlier. But some heuristics make it possible to state that centrally symmetric orbits exist. In this paper we compute generating solutions with central symmetry of the Hill problem.
The method of computing centrally symmetric generating solutions is based on iterative calculation of a normal form of the Hill problem Hamiltonian in the vicinity of periodic solutions of the Kepler problem in the uniformly rotating frame. The method of invariant normalization (symmetrization) proposed by academician V. F. Zhuravlev was applied. The method is based on integration of so called homological equation for obtaining both Lie generator $\mathcal G$ of normal transformation and normal form $\mathcal H$.
The normal form $\mathcal H_2$ of the second order of initial Hamiltonian $H$ was computed and condition on existence of centrally symmetric generating solutions was formulated. Their initial conditions and periods were computed as well. These generating solutions can be used for their numerical continuation along the parameter $\varepsilon$ up to periodic solutions with central symmetry of the Hill problem.