Existence proof of librational invariant tori in an averaged model of HD60532 planetary system

We investigate the long-term dynamics of HD60532, an extrasolar system hosting two giant planets orbiting in a 3:1 mean motion resonance. We consider an average approximation at order one in the masses which results (after the reduction in the constants of motion) in a resonant Hamiltonian with two libration angles. In this framework, the usual algorithms constructing the Kolmogorov normal form approach do not easily apply and we need to perform some untrivial preliminary operations, in order to adapt the method to this kind of problems. First, we perform an average over the fast angle of libration which provides an integrable approximation of the Hamiltonian. Then, we introduce action-angle variables that are adapted to such an integrable approximation. This sequence of preliminary operations brings the Hamiltonian in a suitable form to successfully start the Kolmogorov normalization scheme. The convergence of the KAM algorithm is proved by applying a technique based on a computer-assisted proof. This allows us to reconstruct the quasi-periodic motion of the system, with initial conditions that are compatible with the observations.