La2010: A new orbital solution for the long term motion of the Earth

We present here a new solution for the astronomical computation of the orbital motion of the Earth spanning from 0 to -250 Myr. The main improvement with respect to the previous numerical solution La2004 (Laskar et al. 2004) is an improved adjustment of the parameters and initial conditions through a fit over 1 Myr to a special version of the high accurate numerical ephemeris INPOP08 (Fienga et al. 2009). The precession equations have also been entirely revised and are no longer averaged over the orbital motion of the Earth and Moon. This new orbital solution is now valid over more than 50 Myr in the past or in the future with proper phases of the eccentricity variations. Due to chaotic behavior, the precision of the solution decreases rapidly beyond this time span, and we discuss the behavior of various solutions beyond 50 Myr. For paleoclimate calibrations, we provide several different solutions that are all compatible with the most precise planetary ephemeris. We have thus reached the time where geological data are now required to discriminate among planetary orbital solutions beyond 50 Myr.Comment: 17 pages, 14 figure

Statistics and Universality in Simplified Models of Planetary Formation

In this paper, we modify Laskar's simplified model of planetary evolution and accretion [J. Laskar, Phys. Rev. Lett, vol 84, p 3240 (2000)] to account for the full conservation of the total angular momentum of the system, and extend it to incorporate an accretion probability that depends on the mass and relative velocity of the colliding particles. We present statistical results for the mass and eccentricity of the planets formed, in terms of their semi-major axes, for a large number of realisations of different versions of the model. In particular, we find that by combining the mass-dependent accretion probability and the velocity-selection mechanism, the planets formed display a systematic occurrence at specific locations. By introducing properly scaled variables, our results are universal with respect to the total angular momentum of the system, the mass of the planetesimal disc, and the mass of the central star.Comment: 13 pages, 21 figures, some in colour. Accepted in MNRA

Andoyer construction for Hill and Delaunay variables

Andoyer variables are well known for the study of rotational dynamics. These variables were derived by Andoyer through a procedure that can be also used to obtain the Hill variables of the Kepler problem. Andoyer construction can also forecast the Delaunay variables which canonicity is then obtained without the use of a generating function.Comment: 8 pages, 2 figures, revised versio

HD60532, a planetary system in a 3:1 mean motion resonance

In a recent paper it was reported a planetary system around the star HD60532, composed by two giant planets in a possible 3:1 mean motion resonance, that should be confirmed within the next decade. Here we show that the analysis of the global dynamics of the system allows to confirm this resonance. The present best fit to data already corresponds to this resonant configuration and the system is stable for at least 5Gry. The 3:1 resonance is so robust that stability is still possible for a wide variety of orbital parameters around the best fit solution and also if the inclination of the system orbital plane with respect to the plane of the sky is as small as 15 deg. Moreover, if the inclination is taken as a free parameter in the adjustment to the observations, we find an inclination ~ 20 deg, which corresponds to M_b =3.1 M_Jup and M_c = 7.4 M_Jup for the planetary companions.Comment: 4 Pages, 4 Figures, accepted by A&

Stability analysis of the Martian obliquity during the Noachian era

We performed numerical simulations of the obliquity evolution of Mars during the Noachian era, at which time the giant planets were on drastically different orbits than today. For the preferred primordial configuration of the planets we find that there are two large zones where the Martian obliquity is stable and oscillates with an amplitude lower than 20$^\circ$. These zones occur at obliquities below 30$^\circ$ and above 60$^\circ$; intermediate values show either resonant or chaotic behaviour depending on the primordial orbits of the terrestrial planets

High order symplectic integrators for perturbed Hamiltonian systems

We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+\epsilon B$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps with a remainder of order $O(\tau^p\epsilon +\tau^2\epsilon^2)$, where $\tau$ is the stepsize of the integrator. The analytical expressions of the leading terms of the remainders are given at all orders. In many cases, a corrector step can be performed such that the remainder becomes $O(\tau^p\epsilon +\tau^4\epsilon^2)$. The performances of these integrators are compared for the simple pendulum and the planetary 3-Body problem of Sun-Jupiter-Saturn.Comment: 24 pages, 6 figurre

AMD-stability and the classification of planetary systems

We present here in full detail the evolution of the angular momentum deficit (AMD) during collisions as it was described in (Laskar, PRL,2000). Since then, the AMD has been revealed to be a key parameter for the understanding of the outcome of planetary formation models. We define here the AMD-stability criterion that can be easily verified on a newly discovered planetary system. We show how AMD-stability can be used to establish a classification of the multiplanet systems in order to exhibit the planetary systems that are long-term stable because they are AMD-stable, and those that are AMD-unstable which then require some additional dynamical studies to conclude on their stability. The AMD-stability classification is applied to the 131 multiplanet systems from The Extrasolar Planet Encyclopaedia database (exoplanet.eu) for which the orbital elements are sufficiently well known.Comment: 18 pages, 13 figures, A&A in pres

Dissipation in planar resonant planetary systems

Close-in planetary systems detected by the Kepler mission present an excess of periods ratio that are just slightly larger than some low order resonant values. This feature occurs naturally when resonant couples undergo dissipation that damps the eccentricities. However, the resonant angles appear to librate at the end of the migration process, which is often believed to be an evidence that the systems remain in resonance. Here we provide an analytical model for the dissipation in resonant planetary systems valid for low eccentricities. We confirm that dissipation accounts for an excess of pairs that lie just aside from the nominal periods ratios, as observed by the Kepler mission. In addition, by a global analysis of the phase space of the problem, we demonstrate that these final pairs are non-resonant. Indeed, the separatrices that exist in the resonant systems disappear with the dissipation, and remains only a circulation of the orbits around a single elliptical fixed point. Furthermore, the apparent libration of the resonant angles can be explained using the classical secular averaging method. We show that this artifact is only due to the severe damping of the amplitudes of the eigenmodes in the secular motion.Comment: 18 pages, 20 figures, accepted to A&