Numerical approach to second-order canonical perturbation theory in the planetary 3-body problem: Application to exoplanets

Abstract

Extrasolar planetary systems commonly exhibit planets on eccentric orbits, with many systems located near or within mean-motion resonances, showcasing a wide diversity of orbital architectures. Such complex systems challenge traditional secular theories, which are limited to first-order approximations in planetary masses or rely on expansions in orbital elements--eccentricities, inclinations, and semi-major axis ratios--that are subject to convergence issues, especially in highly eccentric, inclined, or tightly-packed systems. To overcome these limitations, we develop a numerical approach to second-order perturbation theory based on the Lie transform formalism. Our method avoids the need for expansions in orbital elements, ensuring broader applicability and more robust convergence. We first outline the Hamiltonian framework for the 3-body planetary problem, and apply a canonical transformation to eliminate fast angle dependencies, deriving the secular Hamiltonian up to second order in the mass ratio. We then use the fast Fourier transform algorithm to numerically simulate, in an accurate way, the long-term evolution of planetary systems near or away from mean-motion resonances. Finally, we validate our methods against well-known planetary configurations, such as the Sun-Jupiter-Saturn system, as well as to exoplanetary systems like WASP-148, TIC 279401253 and GJ 876, demonstrating the applicability of our models across a wide range of planetary configurations.

arXiv comment17 pages, 9 figures. Published in Physical Review D

References in corpus (17)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• The Occurrence and Architecture of Exoplanetary Systems
• Architecture and Dynamics of Kepler's Candidate Multiple Transiting Planet Systems
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• Secular Evolution of Hierarchical Planetary Systems
• Analytical Treatment of Planetary Resonances
• Explicit expansion of the three-body disturbing function for arbitrary eccentricities and inclinations
• High precision Symplectic Integrators for the Solar System
• Dynamical stability analysis of the HD202206 system and constraints to the planetary orbits
• Dissipation in planar resonant planetary systems
• On the extension of the Laplace-Lagrange secular theory to order two in the masses for extrasolar systems
• Discovery and characterization of the exoplanets WASP-148b and c. A transiting system with two interacting giant planets
• Long-term dynamics of the solar system inner planets
• The origin of chaos in the Solar System through computer algebra
• Resonant Laplace-Lagrange theory for extrasolar systems in mean-motion resonance
• A 2:1 Mean-Motion Resonance Super-Jovian pair revealed by TESS, FEROS, and HARPS
• Photodynamical analysis of the nearly resonant planetary system WASP-148: Accurate transit-timing variations and mutual orbital inclination