Dynamical evolution and stability maps of the Proxima Centauri system

TL;DR

This study investigates the dynamical stability of the Proxima Centauri planetary system using numerical simulations and semi-analytical models, revealing stable orbital configurations and constraints on planetary eccentricities and inclinations.

Contribution

It provides the first comprehensive dynamical analysis of Proxima Centauri's planets, incorporating relativistic and tidal effects, and establishes stability ranges for orbital parameters.

Findings

Proxima b's orbit is stable for semi-major axes 0.02-0.1 au with eccentricity < 0.4.

Relativistic effects influence eccentricity evolution, while tidal effects affect long-term eccentricity.

Mutual inclinations below 50° are necessary for system stability.

Abstract

Proxima Centauri was recently discovered to host an Earth-mass planet of Proxima b, and a 215-day signal which is probably a potential planet c. In this work, we investigate the dynamical evolution of the Proxima Centauri system with the full equations of motion and semi-analytical models including relativistic and tidal effects. We adopt the modified Lagrange-Laplace secular equations to study the evolution of eccentricity of Proxima b, and find that the outcomes are consistent with those from the numerical simulations. The simulations show that relativistic effects have an influence on the evolution of eccentricities of planetary orbits, whereas tidal effects primarily affects the eccentricity of Proxima b over long timescale. Moreover, using the MEGNO (the Mean Exponential Growth factor of Nearby Orbits) technique, we place dynamical constraints on orbital parameters that result in stable or quasi-periodic motions for coplanar and non-coplanar configurations. In the coplanar case, we find that the orbit of Proxima b is stable for the semi-major axis ranging from 0.02 au to 0.1 au and the eccentricity being less than 0.4. This is where the best-fitting parameters for Proxima b exactly fall. Additional simulations show that the robust stability of this system would favor an eccentricity of Proxima b less than 0.45 and that of Proxima c below 0.65. In the non-coplanar case, we find that mutual inclinations of two planets must be lower than $50^{\circ}$ in order to provide stability. Finally, we estimate the mass of Proxima c to be $3.13~{M_ \oplus } \le {m_c} \le 70.7~{M_ \oplus }$ when $1.27~{M_ \oplus } \le {m_b} \le 1.6~{M_ \oplus }$, if ${i_{mutual}} \le {50^{\circ}}$ and $\Delta \Omega = 0^{\circ}$.