Computer-assisted proofs of existence of KAM tori in planetary dynamical models of $\upsilon$-And $\mathbf{b}$

Rita Mastroianni,Ugo Locatelli

2023-11-03

Abstract:We reconsider the problem of the orbital dynamics of the innermost exoplanet of the $\upsilon$-Andromedae system (i.e., $\upsilon$-And $\mathbf{b}$) into the framework of a Secular Quasi-Periodic Restricted Hamiltonian model. This means that we preassign the orbits of the planets that are expected to be the biggest ones in that extrasolar system (namely, $\upsilon$-And $\mathbf{c}$ and $\upsilon$-And $\mathbf{d}$). The Fourier decompositions of their secular motions are injected in the equations describing the orbital dynamics of $\upsilon$-And $\mathbf{b}$ under the gravitational effects exerted by those two exoplanets. By a computer-assisted procedure, we prove the existence of KAM tori corresponding to orbital motions that we consider to be very robust configurations, according to the analysis and the numerical explorations made in our previous article. The computer-assisted assisted proofs are successfully performed for two variants of the Secular Quasi-Periodic Restricted Hamiltonian model, which differs for what concerns the effects of the relativistic corrections on the orbital motion of $\upsilon$-And $\mathbf{b}$, depending on whether they are considered or not.

Mathematical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is related to the orbital stability of the innermost planet υ - And b in the υ - Andromedæ system. Specifically, the authors verified the existence of KAM tori through computer - aided proof methods, both with and without considering general - relativistic corrections. These tori correspond to very stable orbital configurations. The main objective of the paper is to prove the existence of these KAM tori in a simplified but realistic Hamiltonian model through rigorous mathematical methods, thereby providing theoretical support for understanding the long - term orbital stability of υ - And b.

### Background and Problem Description
The υ - Andromedæ system is a star system containing four known exoplanets. These planets are named υ - And b, υ - And c, υ - And d, and υ - And e in order of their distances from the host star. Due to the gravitational influence of the other larger planets on the innermost planet υ - And b, its orbital stability has become an important research problem. Numerical simulations show that unstable motions are common in this system, so more rigorous mathematical methods are required to analyze its orbital dynamics.

### Research Method
The authors adopted a simplified model, namely the **Secular Quasi - Periodic Restricted (SQPR) Hamiltonian model**. This model simplifies the orbital - dynamics problem of υ - And b into a Hamiltonian system with 2 + 3/2 degrees of freedom. In this model, the orbital motions of the larger planets υ - And c and υ - And d are predetermined and injected into the equations describing the orbital dynamics of υ - And b through Fourier decomposition.

### Computer - Aided Proof
To prove the existence of KAM tori, the authors used a computer - aided proof (CAP) technique based on the normal - form method. The specific steps are as follows:
1. **Selection of Initial Conditions**: First, a most stable orbital configuration is selected as the initial condition through numerical methods.
2. **Construction of Normal Form**: Apply the normal - form algorithm to gradually eliminate the perturbation terms in the Hamiltonian function and finally transform the Hamiltonian function into an approximate normal form.
3. **Application of KAM Theorem**: Use the KAM theorem to prove the existence of invariant tori under specific conditions. These tori correspond to stable orbital configurations.

### Results and Conclusions
The paper successfully proved the existence of KAM tori, both with and without considering general - relativistic corrections. These results not only provide a theoretical basis for understanding the orbital stability of υ - And b but also offer a new method for the dynamic analysis of similar systems.

### Formula Presentation
Some of the key formulas involved in the paper are as follows:

1. **Hamiltonian Function**:
\[
H=\sum_{j = 2}^{3}\left(\frac{\mathbf{p}_j\cdot\mathbf{p}_j}{2\beta_j}-\frac{Gm_0m_j}{r_j}\right)+\frac{\mathbf{p}_2\cdot\mathbf{p}_3}{m_0}-\frac{Gm_2m_3}{|\mathbf{r}_2-\mathbf{r}_3|}
\]
where \(m_0\) is the mass of the star, \(m_j\) is the mass of the planet, \(\mathbf{r}_j\) and \(\mathbf{p}_j\) are the position vector and conjugate momentum of the planet respectively, \(G\) is the gravitational constant, and \(\beta_j=\frac{m_0m_j}{m_0 + m_j}\) is the reduced mass.

2. **Quasi - Periodic Approximation**:
\[
\xi_j(t)+i\eta_j(t)\approx\sum_{s = 1}^{N_C}A_{j,s}e^{i(\mathbf{k}_{j,s}\cdot\boldsymbol{\theta}(t)+\vartheta_{j,s})}
\]
\[
P_j(t)+iQ_j(t)