On the nature of the motion of a test particle in the pseudo-Newtonian Hill system

The scope of this work is to perform a numerical investigation of the orbital dynamics for a test particle in the pseudo-Newtonian Hill problem. Large two-dimensional sets of initial conditions of prograde and retrograde orbits are investigated. The orbits are classified as bounded (chaotic, sticky or regular), escaping and collision orbits. The smaller alignment index (SALI) method is used to identify chaotic orbits. Additionally, the influence of the energy (or equivalently the value of the Jacobi constant) and of the Schwarzschild radius on the orbital structure of the system are determined. Our numerical results are compared with related previous ones, corresponding to the classical version of the Hill problem.