Long-term dynamics of the solar system inner planets

Although the discovery of the chaotic motion of the inner planets in the solar system dates back to more than thirty years ago, the secular chaos of their orbits still dares more analytical analyses. Apart from the high-dimensional structure of the motion, this is probably related to the lack of an adequately simple dynamical model. Here, we consider a new secular dynamics for the inner planets, with the aim of retaining a fundamental set of interactions responsible for their chaotic behaviour, while being consistent with the predictions of the most precise orbital solutions currently available. We exploit the regularity in the secular motion of the outer planets, to predetermine a quasi-periodic solution for their orbits. This reduces the secular phase space to the degrees of freedom dominated by the inner planets. On top of that, the smallness of the inner planet masses and the absence of strong mean-motion resonances permits to restrict ourselves to first-order secular averaging. The resulting dynamics can be integrated numerically in a very efficient way through Gauss's method, while computer algebra allows for analytical inspection of planet interactions, once the Hamiltonian is truncated at a given total degree in eccentricities and inclinations. The new model matches very satisfactorily reference orbital solutions of the solar system over timescales shorter than or comparable to the Lyapunov time. It correctly reproduces the maximum Lyapunov exponent of the inner system and the statistics of the high eccentricities of Mercury over the next five billion years. The destabilizing role of the $g_1-g_5$ secular resonance also arises. A numerical experiment, consisting of a thousand orbital solutions over one hundred billion years, reveals the essential properties of the stochastic process driving the destabilization of the inner solar system and clarifies its current metastable state.