Abstract
The JUICE and Europa Clipper missions will soon be launched to investigate the presence of liquid subsurface oceans underneath the icy moons of Jupiter. Both missions rely on complex moon tour trajectories whose design can benefit from a preliminary understanding of the moons' relative geometry when candidate flyby sequences are researched. To aid with this objective, a new time-periodic system of equation of motion is introduced based on the 4:2:1 Laplace resonance of the Io-Europa-Ganymede system. The dynamical system is investigated with numerical continuation techniques that seek to grow equilibria and periodic orbits from the circular restricted three-body problem into periodic and quasi-periodic solutions that populate the phase space when fourth-and fifth-body perturbations are taken into account. The stability of these solutions is further analyzed in the new time-periodic model, searching for stable and unstable manifolds that can help identify fuel-efficient transfer opportunities between the Galilean moons of Jupiter.