We present the discovery and characterization of WASP-148, a new extrasolar system that includes at least two giant planets. The host star is a slowly rotating inactive late-G dwarf with a V=12 magnitude. The planet WASP-148b is a hot Jupiter of 0.72 R_Jup and 0.29 M_Jup that transits its host with an orbital period of 8.80 days. We found the planetary candidate with the SuperWASP photometric survey, then characterized it with the SOPHIE spectrograph. Our radial velocity measurements subsequently revealed a second planet in the system, WASP-148c, with an orbital period of 34.5 days and a minimum mass of 0.40 M_Jup. No transits of this outer planet were detected. The orbits of both planets are eccentric and fall near the 4:1 mean-motion resonances. This configuration is stable on long timescales, but induces dynamical interactions so that the orbits differ slightly from purely Keplerian orbits. In particular, WASP-148b shows transit-timing variations of typically 15 minutes, making it the first interacting system with transit-timing variations that is detected on ground-based light curves. We establish that the mutual inclination of the orbital plane of the two planets cannot be higher than 35 degrees, and the true mass of WASP-148c is below 0.60 M_Jup. We present photometric and spectroscopic observations of this system that cover a time span of ten years. We also provide their Keplerian and Newtonian analyses; these analyses should be significantly improved through future TESS~observations.

Sky survey telescopes and powerful targeted telescopes play complementary roles in astronomy. In order to investigate the nature and characteristics of the motions of very faint objects, a flexibly-pointed instrument capable of high astrometric accuracy is an ideal complement to current astrometric surveys and a unique tool for precision astrophysics. Such a space-based mission will push the frontier of precision astrometry from evidence of earth-massed habitable worlds around the nearest starts, and also into distant Milky way objects up to the Local Group of galaxies. As we enter the era of the James Webb Space Telescope and the new ground-based, adaptive-optics-enabled giant telescopes, by obtaining these high precision measurements on key objects that Gaia could not reach, a mission that focuses on high precision astrometry science can consolidate our theoretical understanding of the local universe, enable extrapolation of physical processes to remote redshifts, and derive a much more consistent picture of cosmological evolution and the likely fate of our cosmos. Already several missions have been proposed to address the science case of faint objects in motion using high precision astrometry ESA missions: NEAT for M3, micro-NEAT for S1 mission, and Theia for M4 and M5. Additional new mission configurations adapted with technological innovations could be envisioned to pursue accurate measurements of these extremely small motions. The goal of this white paper is to address the fundamental science questions that are at stake when we focus on the motions of faint sky objects and to briefly review quickly instrumentation and mission profiles.

In some planetary systems, the orbital periods of two of its members present a commensurability, usually known by mean-motion resonance. These resonances greatly enhance the mutual gravitational influence of the planets. As a consequence, these systems present uncommon behaviors, and their motions need to be studied with specific methods. Some features are unique and allow us a better understanding and characterization of these systems. Moreover, mean-motion resonances are a result of an early migration of the orbits in an accretion disk, so it is possible to derive constraints on their formation. Here we review the dynamics of a pair of resonant planets and explain how their orbits evolve in time. We apply our results to the HD 45365 planetary system.

We use the planetary ephemeris INPOP17b to constrain the mass of the graviton in the Newtonian limit. We also give an interpretation of this result for a specific case of fifth force framework. We find that the residuals for the Cassini spacecraft significantly (90\% C.L.) degrade for Compton wavelengths of the graviton smaller than $1.83\times 10^{13}$ km, corresponding to a graviton mass bigger than $6.76\times 10^{-23} eV/c^2$. This limit is comparable in magnitude to the one obtained by the LIGO-Virgo collaboration in the radiative regime. We also use this specific example to illustrate why constraints on alternative theories of gravity obtained from postfit residuals are generically overestimated.

We present the obliquity measurement, that is, the angle between the normal angle of the orbital plane and the stellar spin axis, of the sub-Neptune planet HD3167 c, which transits a bright nearby K0 star. We study the orbital architecture of this multi-planet system to understand its dynamical history. We also place constraints on the obliquity of planet d based on the geometry of the planetary system and the dynamical study of the system. New observations obtained with HARPS-N at the Telescopio Nazionale Galileo (TNG) were employed for our analysis. The sky-projected obliquity was measured using three different methods: the Rossiter-McLaughlin anomaly, Doppler tomography, and reloaded Rossiter-McLaughlin techniques. We performed the stability analysis of the system and investigated the dynamical interactions between the planets and the star. HD3167 c is found to be nearly polar with sky-projected obliquity, $\lambda$ = -97$^{\circ} \pm$ 23$^{\circ}$. This misalignment of the orbit of planet c with the spin axis of the host star is detected with 97\% confidence. The analysis of the dynamics of this system yields coplanar orbits of planets c and d. It also shows that it is unlikely that the currently observed system can generate this high obliquity for planets c and d by itself. However, the polar orbits of planets c and d could be explained by the presence of an outer companion in the system. Follow-up observations of the system are required to confirm such a long-period companion.

Eccentricity is a parameter of particular interest as it is an informative indicator of the past of planetary systems. It is however not always clear whether the eccentricity fitted on radial velocity data is real or if it is an artefact of an inappropriate modelling. In this work, we address this question in two steps: we first assume that the model used for inference is correct and present interesting features of classical estimators. Secondly, we study whether the eccentricity estimates are to be trusted when the data contain incorrectly modelled signals, such as missed planetary companions, non Gaussian noises, correlated noises with unknown covariance, etc. Our main conclusion is that data analysis via posterior distributions, with a model including a free error term gives reliable results provided two conditions. First, convergence of the numerical methods needs to be ascertained. Secondly, the noise power spectrum should not have a particularly strong peak at the semi period of the planet of interest. As a consequence, it is difficult to determine if the signal of an apparently eccentric planet might be due to another inner companion in 2:1 mean motion resonance. We study the use of Bayes factors to disentangle these cases. Finally, we suggest methods to check if there are hints of an incorrect model in the residuals. We show on simulated data the performance of our methods and comment on the eccentricities of Proxima b and 55 Cnc f.

We present a new mixed variable symplectic (MVS) integrator for planetary systems, that fully resolve close encounters. The method is based on a time regularisation that allows keeping the stability properties of the symplectic integrators, while also reducing the effective step size whenever two planets encounter. We use a high order MVS scheme such that it is possible to integrate with large time steps far away from close encounters. We show that this algorithm is able to resolve almost exact collisions (i.e with a mutual separation of a fraction of the physical radius) while using the same time-step as in weakly perturbed problem such as the Solar System. We demonstrate the long term behaviour on systems of six super-Earths experiencing strong scattering for 50 kyr. We compare our algorithm to hybrid methods such as MERCURY and show that for an equivalent cost we obtain much better energy conservation.

Lunar laser ranging (LLR) data and Apollo seismic data analyses, revealed independent evidence for the presence of a fluid lunar core. However, the size of the lunar fluid core remained uncertain by $\pm55$ km (encompassing two contrasting 2011 Apollo seismic data analyses). Here we show that a new description of the lunar interior's dynamical model provides a determination of the radius and geometry of the lunar core-mantle boundary (CMB) from the LLR observations. We compare the present-day lunar core oblateness obtained from LLR analysis with the expected hydrostatic model values, over a range of previously expected CMB radii. The findings suggest a core oblateness ($f_c=(2.2\pm0.6)\times10^{-4}$) that satisfies the assumption of hydrostatic equilibrium over a tight range of lunar CMB radii ($\mathcal{R}_{CMB}=381\pm12$ km). Our estimates of a presently-relaxed lunar CMB translates to a core mass fraction in the range of $1.59-1.77\%$ with a present-day Free Core Nutation (FCN) within $(367\pm100)$ years.

We propose to use the properties of the Lie algebra of the angular momentum to build symplectic integrators dedicated to the Hamiltonian of the free rigid body. By introducing a dependence of the coefficients of integrators on the moments of inertia of the integrated body, we can construct symplectic dedicated integrators with fewer stages than in the general case for a splitting in three parts of the Hamiltonian. We perform numerical tests to compare the developed dedicated 4th-order integrators to the existing reference integrators for the water molecule. We also estimate analytically the accuracy of these new integrators for the set of the rigid bodies and conclude that they are more accurate than the existing ones only for very asymmetric bodies.

Context. The dwarf planet Ceres and the asteroid Vesta have been studied by the Dawn space mission. They are the two heaviest bodies of the main asteroid belt and have different characteristics. Notably, Vesta appears to be dry and inactive with two large basins at its south pole. Ceres is an ice-rich body with signs of cryovolcanic activity. Aims. The aim of this paper is to determine the obliquity variations of Ceres and Vesta and to study their rotational stability. Methods. The orbital and rotational motions have been integrated by symplectic integration. The rotational stability has been studied by integrating secular equations and by computing the diffusion of the precession frequency. Results. The obliquity variations of Ceres over $[-20:0]\,\mathrm{Myr}$ are between $2$ and $20{\deg}$ and the obliquity variations of Vesta are between $21$ and $45{\deg}$. The two giant impacts suffered by Vesta modified the precession constant and could have put Vesta closer to the resonance with the orbital frequency $2s_6-s_V$. Given the uncertainty on the polar moment of inertia, the present Vesta could be in this resonance where the obliquity variations can vary between $17$ and $48{\deg}$. Conclusions. Although Ceres and Vesta have precession frequencies close to the secular orbital frequencies of the inner planets, their long-term rotations are relatively stable. The perturbations of Jupiter and Saturn dominate the secular orbital dynamics of Ceres and Vesta and the perturbations of the inner planets are much weaker. The secular resonances with the inner planets also have smaller widths and do not overlap, contrary to the case of the inner planets.

Context: Seasonal variations and climate stability of a planet are very sensitive to the planet obliquity and its evolution. This is of particular interest for the emergence and sustainability of land-based life, but orbital and rotational parameters of exoplanets are still poorly constrained. Numerical explorations usually realised in this situation are thus in heavy contrast with the uncertain nature of the available data. Aims: We aim to provide an analytical formulation of the long-term spin-axis dynamics of exoplanets, linking it directly to physical and dynamical parameters, but still giving precise quantitative results if the parameters are well known. Together with bounds for the poorly constrained parameters of exoplanets, this analysis is designed to allow a quick and straightforward exploration of the spin-axis dynamics. Methods: The long-term orbital solution is decomposed in quasi-periodic series and the spin-axis Hamiltonian is expanded in powers of eccentricity and inclination. Chaotic zones are measured by the resonance overlap criterion. Bounds for the poorly known parameters of exoplanets are obtained from physical grounds (rotational breakup) and dynamical considerations (equipartition of AMD). Results: This method gives accurate results when the orbital evolution is well known. The chaotic zones for planets of the Solar System can be retrieved in details from simple analytical formulas. For less constrained planetary systems, the maximal extent of the chaotic regions can be computed, requiring only the mass, the semi-major axis and the eccentricity of the planets present in the system. Additionally, some estimated bounds of the precession constant allow to classify which observed exoplanets are necessarily out of major spin-orbit secular resonances (unless the precession rate is affected by the presence of massive satellites).

Short-period planets suffer from extreme tidal forces from their parent stars causing them to deform and attain non-spherical shapes. The planet shapes, modeled here as triaxial ellipsoids, can have an impact on the observed transit light curves and the planetary parameters derived. We investigate the detectability of tidal deformation in short-period planets from their transit light curves and the instrumental precision needed. We show how the detection of deformation from the light curve allows us to obtain an observational estimate of the second fluid Love number which gives valuable insight about the planet's internal structure. We adopted a model to calculate the shape of a planet due to the external potentials acting on it and used this model to modify the ellc transit tool. Our model is parameterized by the Love number, hence for a given light curve we can derive the value of the Love number that best matches the observations. We simulated the known cases of WASP-103b and WASP-121b expected to be highly deformed. Our analyses showed that instrumental precision $\leq$50ppm/min is needed to reliably estimate the Love number and detect tidal deformation. This precision can be achieved for WASP-103b in ~40 transits using HST and in ~300 transits using the forthcoming CHEOPS. However, fewer transits will be required for short-period planets that may be found around bright stars in the TESS and PLATO survey missions. The unprecedented precisions expected from PLATO and JWST can permit the detection of deformation with a single transit. However, the effects of instrumental and astrophysical noise must be well-considered as they can increase the number of transits required to reach the 50 ppm/min detection limit. We also show that improper modeling of limb darkening can act to bury signals related to the planet's shape thereby leading us to infer sphericity for a deformed planet.

Oceanic tides are a major source of tidal dissipation. They are a key actor for the orbital and rotational evolution of planetary systems, and contribute to the heating of icy satellites hosting a subsurface ocean. Oceanic tides are characterized by a highly frequency-resonant behavior, which is mainly due to the propagation of surface gravity waves in the case of thin oceans, and internal waves when they are deeper. In this work, we derive self-consistent ab initio expressions of the oceanic tidal torque as a function of the key physical parameters of the system (the ocean depth, the Brunt-V\"ais\"al\"a stratification frequency, the rotation rate, the tidal frequency, the Rayleigh friction). These solutions include the coupled mechanisms of internal and surface gravito-inertial waves, which allows us to study the case of planets hosting deep oceans and offer interesting prospects for the coupling between subsurface oceans and ice shells in the case of icy satellites.

In a two-planet system, due to Sundman (1912) inequality, a topological boundary can forbid close encounters between the two planets for infinite time. A system is said Hill stable if it verifies this topological condition. Hill stability is widely used in the study of extra solar planets dynamics. However people often use the coplanar and circular orbits approximation. In this paper, we explain how the Hill stability can be understood in the framework of Angular Momentum Deficit (AMD). In the secular approximation, the AMD allows to discriminate between a priori stable systems and systems for which a more in depth dynamical analysis is required. We show that the general Hill stability criterion can be expressed as a function of only the semi major axes, the masses and the total AMD of the system. The proposed criterion is only expanded in the planets-to-star mass ratio $\epsilon$ and not in the semi-major axis ratio, in eccentricities nor in the mutual inclination. Moreover the expansion in $\epsilon$ remains excellent up to values of about $10^{-3}$ even for two planets with very different mass values. We performed numerical simulations in order to highlight the sharp change of behaviour between Hill stable and Hill unstable systems. We show that Hill stable systems tend to be very regular whereas Hill unstable ones often lead to rapid planet collisions. We also remind that Hill stability does not protect from the ejection of the outer planet.

Cyclic sedimentation has varied at several timescales and this variability has been geologically well documented at Milankovitch timescales, controlled in part by climatically (insolation) driven sea-level changes. At the longer (tens of Myr) timescales connection between astronomical parameters and sedimentation via cyclic solar-system motions within the Milky Way has also been proposed, but this hypothesis remains controversial because of the lack of long geological records. The absence of a physical mechanism that could explain the connection between climate and astronomy at these longer timescales led to the explanation of plate motions as the main driver of climate on Earth. Here we statistically show a prominent and persistent ~36 Myr sedimentary cyclicity superimposed on two megacycles (~250 Myr) in a relatively well-constrained sea-level (SL) record of the past 542 Myr (Phanerozoic eon). Given the possible link between amplitudes of the ~36 and ~250 Myr cyclicities in SL record and the potential that these periodicities fall into the frequency band of solar system motions, we suggest an astronomical origin, and model these periodicities as originating from the path of the solar system in the Milky Way as vertical and radial periods that modulate the flux of cosmic rays on Earth. Our finding of the ~36 Myr SL cyclicity lends credibility to the existing hypothesis about the imprint of solar-system vertical period on the geological record. The ~250 Myr megacycles are tentatively attributed to a radial period. However, the tectonic drivers also remain potentially plausible. The potential existence of a correlation between the modeled astronomical signal and the geological record may offer an indirect proxy to understand the structure and history of the Milky Way by providing a 542 Myr long record of the path of the Sun in our Galaxy.

Oceanic tides are a major source of tidal dissipation. They drive the evolution of planetary systems and the rotational dynamics of planets. However, 2D models commonly used for the Earth cannot be applied to extrasolar telluric planets hosting potentially deep oceans because they ignore the three-dimensional effects related to the ocean vertical structure. Our goal is to investigate in a consistant way the importance of the contribution of internal gravity waves in the oceanic tidal response and to propose a modeling allowing to treat a wide range of cases from shallow to deep oceans. A 3D ab initio model is developed to study the dynamics of a global planetary ocean. This model takes into account compressibility, stratification and sphericity terms, which are usually ignored in 2D approaches. An analytic solution is computed and used to study the dependence of the tidal response on the tidal frequency and on the ocean depth and stratification. In the 2D asymptotic limit, we recover the frequency-resonant behaviour due to surface inertial-gravity waves identified by early studies. As the ocean depth and Brunt-V\"ais\"al\"a frequency increase, the contribution of internal gravity waves grows in importance and the tidal response become three-dimensional. In the case of deep oceans, the stable stratification induces resonances that can increase the tidal dissipation rate by several orders of magnitude. It is thus able to affect significantly the evolution time scale of the planetary rotation.

Semidiurnal atmospheric thermal tides are important for terrestrial exoplanets in the habitable zone of their host stars. With solid tides, they torque these planets, thus contributing to determine their rotation states as well as their climate. Given the complex dynamics of thermal tides, analytical models are essential to understand its dependence on the structure and rotation of planetary atmospheres and the tidal frequency. In this context, the state of the art model proposed in the 60's by Lindzen and Chapman explains well the properties of thermal tides in the asymptotic regime of Earth-like rapid rotators but predicts a non-physical diverging tidal torque in the vicinity of the spin-orbit synchronization. In this work, we present a new model that addresses this issue by taking into account dissipative processes through a Newtonian cooling. First, we recover the tidal torque recently obtained with numerical simulations using General Circulation Models (GCM). Second, we show that the tidal response is very sensitive to the atmospheric structure, particularly to the stability with respect to convection. A strong stable stratification is able to annihilate the atmospheric tidal torque, leading to synchronization, while a convective atmosphere will be submitted to a strong torque, leading to a non-synchronized rotation state.

(in english): Urbain-Jean-Joseph Le Verrier was born in Saint-L\^o on March 11, 1811. He entered the Ecole Polytechnique in 1831, from which he was to emerge 8th two years later. After first devoting himself to chemistry, in 1836 he obtained a position as an astronomy assistant at the Ecole Polytechnique. This choice will decide his future career which culminates with the discovery of Neptune in 1846. Le Verrier wrote more than 200 contributions in the Comptes Rendus de l'Acad\'emie des Sciences. These contributions are very varied: some original articles but also reports on publications published elsewhere, sometimes even simple notes of a single page. The whole set gives a very vivid vision of the development of the science of the XIXth century. At that time, the Comptes Rendus are really a reflection of the debates of the sessions of the Academy. They are published very quickly, and leave a large freedom of speech to the authors. They are therefore a snapshot of the sometimes lively polemics which animated the sessions of the Academy of Sciences. In this limited essay, we will look to the first years of the career of Le Verrier until the discovery of Neptune. (in french) Le Verrier a \'ecrit plus de 200 contributions dans les Comptes Rendus de l'Acad\'emie des Sciences . Ces contributions sont le reflet des d\'ebats des s\'eances de l'Acad\'emie. Ils sont publi\'es tr\`es rapidement, et laissent une tr\`es grande libert\'e de parole aux auteurs. On y retrouve de mani\`ere instantan\'ee les pol\'emiques parfois vives qui ont anim\'e les s\'eances de l'Acad\'emie des sciences. Dans cet essai nous survolerons les premi\`eres ann\'ees de la carri\`ere de Le Verrier jusqu'\`a la d\'ecouverte de Neptune.

Thermal atmospheric tides can torque telluric planets away from spin-orbit synchronous rotation, as observed in the case of Venus. They thus participate to determine the possible climates and general circulations of the atmospheres of these planets. In this work, we write the equations governing the dynamics of thermal tides in a local vertically-stratified section of a rotating planetary atmosphere by taking into account the effects of the complete Coriolis acceleration on tidal waves. This allows us to derive analytically the tidal torque and the tidally dissipated energy, which we use to discuss the possible regimes of tidal dissipation and examine the key role played by stratification. In agreement with early studies, we find that the frequency dependence of the thermal atmospheric tidal torque in the vicinity of synchronization can be approximated by a Maxwell model. This behaviour corresponds to weakly stably stratified or convective fluid layers, as observed in ADLM2016a. A strong stable stratification allows gravity waves to propagate, which makes the tidal torque become negligible. The transition is continuous between these two regimes. The traditional approximation appears to be valid in thin atmospheres and in regimes where the rotation frequency is dominated by the forcing or the buoyancy frequencies. Depending on the stability of their atmospheres with respect to convection, observed exoplanets can be tidally driven toward synchronous or asynchronous final rotation rates. The domain of applicability of the traditional approximation is rigorously constrained by calculations.

Atmospheric tides can have a strong impact on the rotational dynamics of planets. They are of most importance for terrestrial planets located in the habitable zone of their host star, where their competition with solid tides is likely to drive the body towards non-synchronized rotation states of equilibrium, as observed in the case of Venus. Contrary to other planetary layers, the atmosphere is sensitive to both gravitational and thermal forcings, through a complex dynamical coupling involving the effects of Coriolis acceleration and characteristics of the atmospheric structure. These key physics are usually not taken into account in modelings used to compute the evolution of planetary systems, where tides are described with parametrised prescriptions. In this work, we present a new ab initio modeling of atmospheric tides adapting the theory of the Earth's atmospheric tides (Chapman & Lindzen 1970) to other terrestrial planets. We derive analytic expressions of the tidal torque, as a function of the tidal frequency and parameters characterizing the internal structure (e.g. the Brunt-V\"ais\"al\"a frequency, the radiative frequency, the pressure heigh scale). We show that stratification plays a key role, the tidal torque being strong in the case of convective atmospheres (i.e. with a neutral stratification) and weak in case of atmosphere convectively stable. In a second step, the model is used to determine the non-synchronized rotation states of equilibrium of Venus-like planets as functions of the physical parameters of the system. These results are detailed in Auclair-Desrotour et al. (2017a) and Auclair-Desrotour et al. (2017b).

The AMD-stability criterion allows to discriminate between a-priori stable planetary systems and systems for which the stability is not granted and needs further investigations. AMD-stability is based on the conservation of the Angular Momentum Deficit (AMD) in the averaged system at all orders of averaging. While the AMD criterion is rigorous, the conservation of the AMD is only granted in absence of mean-motion resonances (MMR). Here we extend the AMD-stability criterion to take into account mean-motion resonances, and more specifically the overlap of first order MMR. If the MMR islands overlap, the system will experience generalized chaos leading to instability. The Hamiltonian of two massive planets on coplanar quasi-circular orbits can be reduced to an integrable one degree of freedom problem for period ratios close to a first order MMR. We use the reduced Hamiltonian to derive a new overlap criterion for first order MMR. This stability criterion unifies the previous criteria proposed in the literature and admits the criteria obtained for initially circular and eccentric orbits as limit cases. We then improve the definition of AMD-stability to take into account the short term chaos generated by MMR overlap. We analyze the outcome of this improved definition of AMD-stability on selected multi-planet systems from the Extrasolar Planets Encyclopeadia.

In the context of the ESA M5 (medium mission) call we proposed a new satellite mission, Theia, based on relative astrometry and extreme precision to study the motion of very faint objects in the Universe. Theia is primarily designed to study the local dark matter properties, the existence of Earth-like exoplanets in our nearest star systems and the physics of compact objects. Furthermore, about 15 $\%$ of the mission time was dedicated to an open observatory for the wider community to propose complementary science cases. With its unique metrology system and "point and stare" strategy, Theia's precision would have reached the sub micro-arcsecond level. This is about 1000 times better than ESA/Gaia's accuracy for the brightest objects and represents a factor 10-30 improvement for the faintest stars (depending on the exact observational program). In the version submitted to ESA, we proposed an optical (350-1000nm) on-axis TMA telescope. Due to ESA Technology readiness level, the camera's focal plane would have been made of CCD detectors but we anticipated an upgrade with CMOS detectors. Photometric measurements would have been performed during slew time and stabilisation phases needed for reaching the required astrometric precision.

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.

Over short time intervals planetary ephemerides have been traditionally represented in analytical form as finite sums of periodic terms or sums of Poisson terms that are periodic terms with polynomial amplitudes. Nevertheless, this representation is not well adapted for the evolution of the planetary orbits in the solar system over million of years as they present drifts in their main frequencies, due to the chaotic nature of their dynamics. The aim of the present paper is to develop a numerical algorithm for slowly diffusing solutions of a perturbed integrable Hamiltonian system that will apply to the representation of the chaotic planetary motions with varying frequencies. By simple analytical considerations, we first argue that it is possible to recover exactly a single varying frequency. Then, a function basis involving time-dependent fundamental frequencies is formulated in a semi-analytical way. Finally, starting from a numerical solution, a recursive algorithm is used to numerically decompose the solution on the significant elements of the function basis. Simple examples show that this algorithm can be used to give compact representations of different types of slowly diffusing solutions. As a test example, we show how this algorithm can be successfully applied to obtain a very compact approximation of the La2004 solution of the orbital motion of the Earth over 40 Myr ([-35Myr,5Myr]). This example has been chosen as this solution is widely used for the reconstruction of the climates of the past.

In this paper, we present a formalism designed to model tidal interaction with a viscoelastic body made of Maxwell material. Our approach remains regular for any spin rate and orientation, and for any orbital configuration including high eccentricities and close encounters. The method is to integrate simultaneously the rotation and the position of the planet as well as its deformation. We provide the equations of motion both in the body frame and in the inertial frame. With this study, we generalize preexisting models to the spatial case and to arbitrary multipole orders using a formalism taken from quantum theory. We also provide the vectorial expression of the secular tidal torque expanded in Fourier series. Applying this model to close-in exoplanets, we observe that if the relaxation time is longer than the revolution period, the phase space of the system is characterized by the presence of several spin-orbit resonances, even in the circular case. As the system evolves, the planet spin can visit different spin-orbit configurations. The obliquity is decreasing along most of these resonances, but we observe a case where the planet tilt is instead growing. These conclusions derived from the secular torque are successfully tested with numerical integrations of the instantaneous equations of motion on HD 80606b. Our formalism is also well adapted to close-in super-Earths in multiplanet systems which are known to have non-zero mutual inclinations.

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.

The competition between the torques induced by solid and thermal tides drives the rotational dynamics of Venus-like planets and super-Earths orbiting in the habitable zone of low-mass stars. The tidal responses of the atmosphere and telluric core are related to their respective physical properties and strongly depend on the tidal frequency. The resulting torque determines the possible equilibrium states of the planet's spin. We compute here an analytic expression for the total tidal torque exerted on a Venus-like planet. This expression is used to characterize the equilibrium rotation of the body. Close to the star, the solid tide dominates. Far from it, the thermal tide drives the rotational dynamics of the planet. The transition regime corresponds to the habitable zone, where prograde and retrograde equilibrium states appear. We demonstrate the strong impact of the atmospheric properties and of the rheology of the solid part on the rotational dynamics of Venus-like planets, highlighting the key role played by dissipative mechanisms in the stability of equilibrium configurations.

Atmospheric tides can strongly affect the rotational dynamics of planets. In the family of Earth-like planets, such as Venus, this physical mechanism coupled with solid tides makes the angular velocity evolve over long timescales and determines the equilibrium configurations of their spin. Contrary to the solid core, the atmosphere is submitted to both tidal gravitational potential and insolation flux coming from the star. The complex response of the gas is intrinsically linked to its physical properties. This dependence has to be characterized and quantified to study the large variety of extrasolar planetary systems. We develop a theoretical global model where radiative losses, which are predominant in slowly rotating atmospheres, are taken into account. We analytically compute the tidal perturbation of pressure, density, temperature and velocity field from which we deduce the expressions of atmospheric Love numbers and tidal torque exerted by the star. The dynamics of atmospheric tides depends on the frequency regime of the tidal perturbation: the thermal regime near synchronization and the dynamical regime characterizing fast-rotating planets. The dependence of the torque on the tidal frequency is quantified for Earth-like and Venus-like exoplanets and is in good agreement with the results given by Global Climate Models (GCM) simulations. Introducing dissipative processes such as radiation regularizes the tidal response of the atmosphere, otherwise it is singular at synchronization. We demonstrate the important role played by the physical and dynamical properties of a super-Earth atmosphere (e.g. Coriolis, stratification, background pressure, density, temperature, radiative emission) and point out the key parameters defining tidal regimes (e.g. inertia, Brunt-V\"ais\"al\"a, radiative frequencies, tidal frequency) and characterize the behaviour of the fluid shell with dissipation.

We present a novel approach for analysing radial velocity data that combines two features: all the planets are searched at once and the algorithm is fast. This is achieved by utilizing compressed sensing techniques, which are modified to be compatible with the Gaussian processes framework. The resulting tool can be used like a Lomb-Scargle periodogram and has the same aspect but with much fewer peaks due to aliasing. The method is applied to five systems with published radial velocity data sets: HD 69830, HD 10180, 55 Cnc, GJ 876 and a simulated very active star. The results are fully compatible with previous analysis, though obtained more straightforwardly. We further show that 55 Cnc e and f could have been respectively detected and suspected in early measurements from the Lick observatory and Hobby-Eberly Telescope available in 2004, and that frequencies due to dynamical interactions in GJ 876 can be seen.

We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin-orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star-planet-satellite systems and triple stellar systems.

In this brief communication we provide the rationale for, and the outcome of the International Astronomical Union (IAU) resolution vote at the XXIX-th General Assembly in Honolulu, Hawaii, in 2015, on recommended nominal conversion constants for selected solar and planetary properties. The problem addressed by the resolution is a lack of established conversion constants between solar and planetary values and SI units: a missing standard has caused a proliferation of solar values (e.g., solar radius, solar irradiance, solar luminosity, solar effective temperature and solar mass parameter) in the literature, with cited solar values typically based on best estimates at the time of paper writing. As precision of observations increases, a set of consistent values becomes increasingly important. To address this, an IAU Working Group on Nominal Units for Stellar and Planetary Astronomy formed in 2011, uniting experts from the solar, stellar, planetary, exoplanetary and fundamental astronomy, as well as from general standards fields to converge on optimal values for nominal conversion constants. The effort resulted in the IAU 2015 Resolution B3, passed at the IAU General Assembly by a large majority. The resolution recommends the use of nominal solar and planetary values, which are by definition exact and are expressed in SI units. These nominal values should be understood as conversion factors only, not as the true solar/planetary properties or current best estimates. Authors and journal editors are urged to join in using the standard values set forth by this resolution in future work and publications to help minimize further confusion.

With about 2000 extrasolar planets confirmed, the results show that planetary systems have a whole range of unexpected properties. We present a full investigation of the HD219828 system, a bright metal-rich star for which a hot neptune has previously been detected. We used a set of HARPS, SOPHIE, and ELODIE radial velocities to search for the existence of orbiting companions to HD219828. A dynamical analysis is also performed to study the stability of the system and to constrain the orbital parameters and planet masses. We announce the discovery of a long period (P=13.1years) massive (msini=15.1MJup) companion (HD219828c) in a very eccentric orbit (e=0.81). The same data confirms the existence of a hot-neptune, HD219828b, with a minimum mass of 21 MEarth and a period of 3.83days. The dynamical analysis shows that the system is stable. The HD219828 system is extreme and unique in several aspects. First, among all known exoplanet systems it presents an unusually high mass ratio. We also show that systems like HD219828, with a hot neptune and a long-period massive companion are more frequent than similar systems with a hot jupiter instead. This suggests that the formation of hot neptunes follows a different path than the formation of their hot jovian counterparts. The high mass, long period, and eccentricity of HD219828c also make it a good target for Gaia astrometry as well as a potential target for atmospheric characterisation, using direct imaging or high-resolution spectroscopy. Astrometric observations will allow us to derive its real mass and orbital configuration. If a transit of HD219828b is detected, we will be able to fully characterise the system, including the relative orbital inclinations. With a clearly known mass, HD219828c may become a benchmark object for the range in between giant planets and brown dwarfs.

To explain the unusual distribution of Kuiper Belt objects, several authors have advocated the existence of a super-Earth planet in the outer solar system. It has recently been proposed that a 10 M$_{\oplus}$ object with an orbit of 700 AU semi major axis and 0.6 eccentricity can explain the observed distribution of Kuiper Belt objects around Sedna. Here we use the INPOP planetary ephemerides model as a sensor for testing for an additional body in the solar system. We test the possibility of adding the proposed planet without increasing the residuals of the planetary ephemerides, fitted over the whole INPOP planetary data sample. We demonstrate that the presence of such an object is not compatible with the most sensitive data set, the Cassini radio ranging data, if its true anomaly is in the intervals $[-130^\circ:-100^\circ]$ or $[-65^\circ : 85^\circ]$. Moreover, we find that the addition of this object can reduce the Cassini residuals, with a most probable position given by a true anomaly $v = {117.8^\circ}^{ + 11^\circ}_{ - 10^\circ} $.

We present new radial velocity measurements of eight stars secured with the spectrograph SOPHIE at the 193-cm telescope of the Haute-Provence Observatory allowing the detection and characterization of new giant extrasolar planets. The host stars are dwarfs of spectral types between F5 and K0 and magnitudes between 6.7 and 9.6; the planets have minimum masses M_p sin i between 0.4 to 3.8 M_Jup and orbital periods of several days to several months. The data allow only single planets to be discovered around the first six stars (HD143105, HIP109600, HD35759, HIP109384, HD220842, and HD12484), but one of them shows the signature of an additional substellar companion in the system. The seventh star, HIP65407, allows the discovery of two giant planets, just outside the 12:5 resonance in weak mutual interaction. The last star, HD141399, was already known to host a four-planetary system; our additional data and analyses allow new constraints to be put on it. We present Keplerian orbits of all systems, together with dynamical analyses of the two multi-planetary systems. HD143105 is one of the brightest stars known to host a hot Jupiter, which could allow numerous follow-up studies to be conducted despite this is not a transiting system. The giant planets HIP109600b, HIP109384b, and HD141399c are located in the habitable zone of their host star.

Planetary ephemerides are a good tool for studying general relativity at the scale of our solar system. We present here new evaluations of advances of perihelia for Mercury and Saturn.

Astronomers commonly quote the properties of celestial objects in units of parameters for the Sun, Jupiter, or the Earth. The resolution presented here was proposed by the IAU Inter-Division Working Group on Nominal Units for Stellar and Planetary Astronomy and passed by the XXIXth IAU General Assembly in Honolulu. IAU 2015 Resolution B3 adopts a set of nominal solar, terrestrial, and jovian conversion constants for stellar and (exo)planetary astronomy which are defined to be exact SI values. While the nominal constants are based on current best estimates (CBEs; which have uncertainties, are not secularly constant, and are updated regularly using new observations), they should be interpreted as standard values and not as CBEs. IAU 2015 Resolution B3 adopts five solar conversion constants (nominal solar radius, nominal total solar irradiance, nominal solar luminosity, nominal solar effective temperature, and nominal solar mass parameter) and six planetary conversion constants (nominal terrestrial equatorial radius, nominal terrestrial polar radius, nominal jovian equatorial radius, nominal jovian polar radius, nominal terrestrial mass parameter, and nominal jovian mass parameter).

The XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on recommended zero points for the absolute and apparent bolometric magnitude scales. The resolution was proposed by the IAU Inter-Division A-G Working Group on Nominal Units for Stellar and Planetary Astronomy after consulting with a broad spectrum of researchers from the astronomical community. Resolution B2 resolves the long-standing absence of an internationally-adopted zero point for the absolute and apparent bolometric magnitude scales. Resolution B2 defines the zero point of the absolute bolometric magnitude scale such that a radiation source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0) corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU 2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the value most commonly adopted in recent literature. The nominal total solar irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832. Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based on the IAU 2012 Resolution B2 definition of the astronomical unit.

Based on the use of MESSENGER radiotracking data in the construction of new Mercury ephemerides (Verma et al. 2014) a new planetary ephemerides INPOP13b was built including Mercury improvements but also improvements on the Mars orbit and on the tie of INPOP planetary ephemerides to ICRF in general.

We study the stability of mean-motion resonances (MMR) between two planets during their migration in a protoplanetary disk. We use an analytical model of resonances, and describe the effect of the disk by a migration timescale (T_{m,i}) and an eccentricity damping timescale (T_{e,i}) for each planet (i=1,2 respectively for the inner and outer planet). We show that the resonant configuration is stable if T_{e,1}/T_{e,2} > (e_1/e_2)^2. This general result can be used to put constraints on specific models of disk-planet interactions. For instance, using classical prescriptions for type I migration, we show that when the angular momentum deficit (AMD) of the inner orbit is larger than the outer's orbit AMD, resonant systems must have a locally inverted disk density profile to stay locked in resonance during the migration. This inversion is very untypical of type I migration and our criterion can thus provide an evidence against classical type I migration. That is indeed the case for the Jupiter-mass resonant systems HD 60532b, c (3:1 MMR), GJ 876b, c (2:1 MMR), and HD 45364b, c (3:2 MMR). This result may be an evidence for type II migration (gap opening planets), which is compatible with the large masses of these planets.

We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We derive in a systematic way an independent set of necessary and sufficient conditions to be satisfied by the coefficients of splitting methods to achieve a prescribed order of accuracy. Splitting methods satisfying such (generalized) order conditions are appropriate in particular for the numerical simulation of the Solar System described in Jacobi coordinates. We show that, when using Poincar\'e Heliocentric coordinates, the same order of accuracy may be obtained by imposing an additional polynomial equation on the coefficients of the splitting method. We construct several splitting methods appropriate for each of the two sets of coordinates by solving the corresponding systems of polynomial equations and finding the optimal solutions. The experiments reported here indicate that the efficiency of our new schemes is clearly superior to previous integrators when high accuracy is required.

Paper in honour of Michel Henon, based on a talk presented at Institut Henri Poincare, Paris, 5 Decembre 2013.

In this paper we present a new approach to tidal theory. Assuming a Maxwell viscoelastic rheology, we compute the instantaneous deformation of celestial bodies using a differential equation for the gravity field coefficients. This method allows large eccentricities and it is not limited to quasi-periodic perturbations. It can take into account an extended class of perturbations, including chaotic motions and transient events. We apply our model to some already detected eccentric hot Jupiters and super-Earths in planar configurations. We show that when the relaxation time of the deformation is larger than the orbital period, spin-orbit equilibria arise naturally at half-integers of the mean motion, even for gaseous planets. In the case of super-Earths, these equilibria can be maintained for very low values of eccentricity. Our method can also be used to study planets with complex internal structures and other rheologies.

Based on the new developped planetary ephemerides INPOP13c, determinations of acceptable intervals of General Relativity violation in considering simultaneously the PPN parameters $\beta$, PPN $\gamma$, the flattening of the sun $J_{2}^\odot$ and time variation of the gravitational mass of the sun $\mu$ are obtained in using Monte Carlo simulation coupled with basic genetic algorithm. Possible time variations of the gravitational constant G are also deduced. Discussions are lead about the better choice of indicators for the goodness-of-fit for each run and limits consistent with general relativity are obtained simultaneously.

Planets with masses between 0.1 - 10 M_earth are believed to host dense atmospheres. These atmospheres can play an important role on the planet's spin evolution, since thermal atmospheric tides, driven by the host star, may counterbalance gravitational tides. In this work we study the long-term spin evolution of Earth-sized exoplanets. We generalize previous works by including the effect of eccentric orbits and obliquity. We show that under the effect of tides and core-mantle friction, the obliquity of the planets evolve either to 0 or 180 degrees. The rotation of these planets is also expected to evolve into a very restricted number of equilibrium configurations. In general, none of this equilibria is synchronous with the orbital mean motion. The role of thermal atmospheric tides becomes more important for Earth-sized planets in the habitable zones of their systems, so they cannot be neglected when we search for their potential habitability.

Multi-planetary systems detected by the Kepler mission present an excess of planets close to first-order mean-motion resonances (2:1 and 3:2) but with a period ratio slightly higher than the resonant value. Several mechanisms have been proposed to explain this observation. Here we provide some clues that these near-resonant systems were initially in resonance and reached their current configuration through tidal dissipation. The argument that has been opposed to this scenario is that it only applies to the close-in systems and not to the farthest ones for which the tidal effect is too weak. Using the catalog of KOI of the Kepler mission, we show that the distributions of period ratio among the most close-in planetary systems and the farthest ones differ significantly. This distance dependent repartition is a strong argument in favor of the tidal dissipation scenario.

We analyze the impact of the installation of a lunar laser ranging device on the NTT 3.6m ESO telescope. With such an instrument, the scientific communities of fundamental physics and solar system formation will highly benefit of the only LLR station in the Southern Hemisphere. The quality of the NTT 3.6 meter telescope will also greatly complement the LLR 3.5 meter Apache Point telescope (3.5 m) instrument in the Northern Hemisphere (USA) which is the best instrument for tracking the Moon since 2006. Finally, we also consider the technical characteristics of such installation including the observational constraints.

We study the evolution of two planets around a star, in mean-motion resonance and undergoing tidal effect. We derive an integrable analytical model of mean-motion resonances of any order which reproduce the main features of the resonant dynamics. Using this simplified model, we obtain a criterion showing that depending on the balance of the tidal dissipation in both planets, their final period ratio may stay at the resonant value, increase above, or decrease below the resonant value. Applying this criterion to the two inner planets orbiting GJ163, we deduce that the current period ratio (2.97) could be the outcome of dissipation in the 3:1 MMR provided that the innermost planet is gaseous (slow dissipation) while the second one is rocky (faster dissipation). We perform N-body simulations with tidal dissipation to confirm the results of our analytical model. We also apply our criterion on GJ581b, c (5:2 MMR) and reproduce the current period ratio (2.4) if the inner planet is gaseous and the outer is rocky (as for GJ163). Finally, we apply our model to the Kepler mission's statistics. We show that the excess of planets pairs close to first order MMR but in external circulation, i.e., with period ratios P_out/P_in > (p+1)/p for the resonance (p+1):p, can be reproduced by tidal dissipation in the inner planet. There is no need for any other dissipative mechanism, provided that these systems left the resonance with non-negligible eccentricities.

The current knowledge of Mercury orbit has mainly been gained by direct radar ranging obtained from the 60s to 1998 and by five Mercury flybys made by Mariner 10 in the 70s, and MESSENGER made in 2008 and 2009. On March 18, 2011, MESSENGER became the first spacecraft to orbit Mercury. The radioscience observations acquired during the orbital phase of MESSENGER drastically improved our knowledge of the orbit of Mercury. An accurate MESSENGER orbit is obtained by fitting one-and-half years of tracking data using GINS orbit determination software. The systematic error in the Earth-Mercury geometric positions, also called range bias, obtained from GINS are then used to fit the INPOP dynamical modeling of the planet motions. An improved ephemeris of the planets is then obtained, INPOP13a, and used to perform general relativity tests of PPN-formalism. Our estimations of PPN parameters (beta and gamma?) are more stringent than previous results.

The meter-per-second precision achieved by today velocimeters enables the search for 1-10 M_Earth planets in the habitable zone of cool stars. This paper reports on the detection of 3 planets orbiting GJ163 (HIP19394), a M3 dwarf monitored by our ESO/HARPS search for planets. We made use of the HARPS spectrograph to collect 150 radial velocities of GJ163 over a period of 8 years. We searched the RV time series for coherent signals and found 5 distinct periodic variabilities. We investigated the stellar activity and casted doubts on the planetary interpretation for 2 signals. Before more data can be acquired we concluded that at least 3 planets are orbiting GJ163. They have orbital periods of P_b=8.632+-0.002, P_c=25.63+-0.03 and P_d=604+-8 days and minimum masses msini = 10.6+-0.6, 6.8+-0.9, and 29+-3 M_Earth, respectively. We hold our interpretations for the 2 additional signals with periods P_(e)=19.4 and P_(f)=108 days. The inner pair presents an orbital period ratio of 2.97, but a dynamical analysis of the system shows that it lays outside the 3:1 mean motion resonance. GJ163c, in particular, is a super-Earth with an equilibrium temperature of T_eq = (302+-10) (1-A)^(1/4) K and may lie in the so called habitable zone for albedo values (A=0.34-0.89) moderately higher than that of Earth (A_Earth=0.2-0.3).

Hierarchical two-planet systems, in which the inner body's semi-major axis is between 0.1 and 0.5 AU, usually present high eccentricity values, at least for one of the orbits. As a result of the formation process, one may expect that planetary systems with high eccentricities also have high mutual inclinations. However, here we show that tidal effects combined with gravitational interactions damp the initial mutual inclination to modest values in timescales that are shorter than the age of the system. This effect is not a direct consequence of tides on the orbits, but it results from a secular forcing of the inner planet's flattening. We then conclude that these hierarchical planetary systems are unlikely to present very high mutual inclinations, at least as long as the orbits remain outside the Lidov-Kozai libration areas. The present study can also be extended to systems of binary stars and to planet-satellite systems.

We present a highly scalable algorithm for multiplying sparse multivariate polynomials represented in a distributed format. This algo- rithm targets not only the shared memory multicore computers, but also computers clusters or specialized hardware attached to a host computer, such as graphics processing units or many-core coprocessors. The scal- ability on the large number of cores is ensured by the lacks of synchro- nizations, locks and false-sharing during the main parallel step.

The INPOP ephemerides have known several improvements and evolutions since the first INPOP06 release (Fienga et al. 2008) in 2008. In 2010, anticipating the IAU 2012 resolutions, adjustement of the gravitational solar mass with a fixed astronomical unit (AU) has been for the first time implemented in INPOP10a (Fienga et al. 2011) together with improvements in the asteroid mass determinations. With the latest INPOP10e version, such advancements have been enhanced and studies about solar corona have also been investigated (Verma et al. 2012). The use of planetary ephemerides for several physical applications are presented here from electronic densities of solar slow and fast winds to asteroid mass determinations and tests of general relativity operated with INPOP10a. Perspectives will also be drawn especially related to the analysis of the Messenger spacecraft data for the planetary orbits and future computation of the time variations of the gravitational mass of the sun.

Mercury's eccentricity is chaotic and can increase so much that collisions with Venus or the Sun become possible (Laskar, 1989, 1990, 1994, 2008, Batygin & Laughlin, 2008, Laskar & Gastineau, 2009). This chaotic behavior results from an intricate network of secular resonances, but in this paper, we show that a simple integrable model with only one degree of freedom is actually able to reproduce the large variations in Mercury's eccentricity, with the correct amplitude and timescale. We show that this behavior occurs in the vicinity of the separatrices of the resonance g1-g5 between the precession frequencies of Mercury and Jupiter. However, the main contribution does not come from the direct interaction between these two planets. It is due to the excitation of Venus' orbit at Jupiter's precession frequency g5. We use a multipolar model that is not expanded with respect to Mercury's eccentricity, but because of the proximity of Mercury and Venus, the Hamiltonian is expanded up to order 20 and more in the ratio of semimajor axis. When the effects of Venus' inclination are added, the system becomes nonintegrable and a chaotic zone appears in the vicinity of the separatrices. In that case, Mercury's eccentricity can chaotically switch between two regimes characterized by either low-amplitude circulations or high-amplitude librations.

Since the formulation of the problem by Newton, and during three centuries, astronomers and mathematicians have sought to demonstrate the stability of the Solar System. Thanks to the numerical experiments of the last two decades, we know now that the motion of the planets in the Solar System is chaotic, which prohibits any accurate prediction of their trajectories beyond a few tens of millions of years. The recent simulations even show that planetary collisions or ejections are possible on a period of less than 5 billion years, before the end of the life of the Sun.

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.

Using a Newtonian model of the Solar System with all 8 planets, we perform extensive tests on various symplectic integrators of high orders, searching for the best splitting scheme for long term studies in the Solar System. These comparisons are made in Jacobi and Heliocentric coordinates and the implementation of the algorithms is fully detailed for practical use. We conclude that high order integrators should be privileged, with a preference for the new $(10,6,4)$ method of (Blanes et al., 2012)

Impact basins identified by Mariner 10 and Messenger flyby images provide us a fossilized record of the impactor flux of asteroids on Mercury during the last stages of the early Solar System. The distribution of these basins is not uniform across the surface, and is consistent with a primordial synchronous rotation (Wieczorek et al. 2012). By analyzing the size of the impacts, we show that the distribution for asteroid diameters D < 110 km is compatible with an index power law of 1.2, a value that matches the predicted primordial distribution of the main-belt. We then derive a simple collisional model coherent with the observations, and when combining it with the secular evolution of the spin of Mercury, we are able to reproduce the present 3/2 spin-orbit resonance (about 50% of chances), as well as a primordial synchronous rotation. This result is robust with respect to variations in the dissipation and collisional models, or in the initial spin state of the planet.

The planet Mercury is locked in a spin-orbit resonance where it rotates three times about its spin axis for every two orbits about the Sun. The current explanation for this unique state assumes that the initial rotation of this planet was prograde and rapid, and that tidal torques decelerated the planetary spin to this resonance. When core-mantle boundary friction is accounted for, capture into the 3/2 resonance occurs with a 26% probability, but the most probable outcome is capture into one of the higher-order resonances. Here we show that if the initial rotation of Mercury were retrograde, this planet would be captured into synchronous rotation with a 68% probability. Strong spatial variations of the impact cratering rate would have existed at this time, and these are shown to be consistent with the distribution of pre-Calorian impact basins observed by Mariner 10 and MESSENGER. Escape from this highly stable resonance is made possible by the momentum imparted by large basin-forming impact events, and capture into the 3/2 resonance occurs subsequently under favourable conditions.

Planets close to their host stars are believed to undergo significant tidal interactions, leading to a progressive damping of the orbital eccentricity. Here we show that, when the orbit of the planet is excited by an outer companion, tidal effects combined with gravitational interactions may give rise to a secular increasing drift on the eccentricity. As long as this secular drift counterbalances the damping effect, the eccentricity can increase to high values. This mechanism may explain why some of the moderate close-in exoplanets are observed with substantial eccentricity values.

We numerically estimate the semi-major axis chaotic diffusion of the Vesta family asteroids induced by close encounters with 11 massive main belt asteroids : (1) Ceres, (2) Pallas, (3) Juno, (4) Vesta, (7) Iris, (10) Hygiea, (15) Eunomia, (19) Fortuna, (324) Bamberga, (532) Herculina, (704) Interamnia. We find that most of the diffusion is due to Ceres and Vesta. By extrapolating our results, we are able to constrain the global effect of close encounters with all the main belt asteroids. A comparison of this drift estimate with the one expected for the Yarkovsky effect shows that for asteroids whose diameter is larger than about 40 km, close encounters dominate the Yarkovsky effect. Overall, our findings confirm the standard scenario for the history of the Vesta family.

We present here in full details the linear secular theory with tidal damping that was used to constraint the fit of the HD10180 planetary system in (Lovis et al. 2011). The theory is very general and can provide some intuitive understanding of the final state of a planetary system when one or more planets are close to their central star. We globally recover the results of (Mardling 2007), but we show that in the HD209458 planetary system, the consideration of the tides raised by the central star on the planet lead to believe that the eccentricity of HD209458b is most probably much smaller than 0.01.

Compared to the previous INPOP versions, the INPOP10a planetary and lunar ephemeris has several improvements. For the planets of our solar system, no big change was brought in the dynamics but improvements were implemented in the fitting process, the data sets used in the fit and in the selection of fitted parameters. We report here the main characteristics of the planetary part of INPOP10a like the fit of the product of the Solar mass with the gravitational constant (GM$_{\odot}$) instead of the astronomical unit. Determinations of PPN parameters as well as adjustments of the Sun J2 and of asteroid masses are also presented. New advances of nodes and perihelia of planets were also estimated and are given here. As for INPOP08, INPOP10a provides to the user, positions and velocities of the planets, the moon, the rotation angles of the Earth and the Moon as well as TT-TDB chebychev polynomials at http://www.imcce.fr/inpop

(abridged) A complete census of planetary systems around a volume-limited sample of solar-type stars (FGK dwarfs) in the Solar neighborhood with uniform sensitivity down to Earth-mass planets within their Habitable Zones out to several AUs would be a major milestone in extrasolar planets astrophysics. This fundamental goal can be achieved with a mission concept such as NEAT - the Nearby Earth Astrometric Telescope. NEAT is designed to carry out space-borne extremely-high-precision astrometric measurements sufficient to detect dynamical effects due to orbiting planets of mass even lower than Earth's around the nearest stars. Such a survey mission would provide the actual planetary masses and the full orbital geometry for all the components of the detected planetary systems down to the Earth-mass limit. The NEAT performance limits can be achieved by carrying out differential astrometry between the targets and a set of suitable reference stars in the field. The NEAT instrument design consists of an off-axis parabola single-mirror telescope, a detector with a large field of view made of small movable CCDs located around a fixed central CCD, and an interferometric calibration system originating from metrology fibers located at the primary mirror. The proposed mission architecture relies on the use of two satellites operating at L2 for 5 years, flying in formation and offering a capability of more than 20,000 reconfigurations (alternative option uses deployable boom). The NEAT primary science program will encompass an astrometric survey of our 200 closest F-, G- and K-type stellar neighbors, with an average of 50 visits. The remaining time might be allocated to improve the characterization of the architecture of selected planetary systems around nearby targets of specific interest (low-mass stars, young stars, etc.) discovered by Gaia, ground-based high-precision radial-velocity surveys.

We investigate the dynamical evolution of hierarchical three-body systems under the effect of tides, when the ratio of the orbital semi-major axes is small and the mutual inclination is relatively large (greater than 20 degrees). Using the quadrupolar non-restricted approximation for the gravitational interactions and the viscous linear model for tides, we derive the averaged equations of motion in a vectorial formalism which is suitable to model the long-term evolution of a large variety of exoplanetary systems in very eccentric and inclined orbits. In particular, it can be used to derive constraints for stellar spin-orbit misalignment, capture in Cassini states, tidal-Kozai migration, or damping of the mutual inclination. Because our model is valid for the non-restricted problem, it can be used to study systems of identical mass or for the outer restricted problem, such as the evolution of a planet around a binary of stars. Here, we apply our model to three distinct situations: 1) the HD80606 planetary system, for which we obtain the probability density function distribution for the misalignment angle, with two pronounced peaks of higher probability around 53 and 109 degrees; 2) the HD98800 binary system, for which we show that initial prograde orbits inside the observed disc may become retrograde and vice-versa, only because of tidal migration within the binary stars; 3) the HD11964 planetary system, for which we show that tidal dissipation combined with gravitational perturbations may lead to a decrease in the mutual inclination, and a fast circularization of the inner orbit.

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.

Context. Low-mass extrasolar planets are presently being discovered at an increased pace by radial velocity and transit surveys, opening a new window on planetary systems. Aims. We are conducting a high-precision radial velocity survey with the HARPS spectrograph which aims at characterizing the population of ice giants and super-Earths around nearby solar-type stars. This will lead to a better understanding of their formation and evolution, and yield a global picture of planetary systems from gas giants down to telluric planets. Methods. Progress has been possible in this field thanks in particular to the sub-m/s radial velocity precision achieved by HARPS. We present here new high-quality measurements from this instrument. Results. We report the discovery of a planetary system comprising at least five Neptune-like planets with minimum masses ranging from 12 to 25 M_Earth, orbiting the solar-type star HD 10180 at separations between 0.06 and 1.4 AU. A sixth radial velocity signal is present at a longer period, probably due to a 65-M_Earth object. Moreover, another body with a minimum mass as low as 1.4 M_Earth may be present at 0.02 AU from the star. This is the most populated exoplanetary system known to date. The planets are in a dense but still well-separated configuration, with significant secular interactions. Some of the orbital period ratios are fairly close to integer or half-integer values, but the system does not exhibit any mean-motion resonances. General relativity effects and tidal dissipation play an important role to stabilize the innermost planet and the system as a whole. Numerical integrations show long-term dynamical stability provided true masses are within a factor ~3 from minimum masses. We further note that several low-mass planetary systems exhibit a rather "packed" orbital architecture with little or no space left for additional planets. (Abridged)

The Planetary and Lunar ephemerides INPOP10a version has several improvements in the fitting process, the data sets used in the fit and in the general features of the solution. No big change was brought in the dynamics. As a consequence of these changes, the extrapolation capabilities of INPOP10a are improved compared to INPOP08, especially for the Earth, Mars, Mercury and Saturn orbits. As for INPOP08, INPOP10a provides to the user, positions, velocities of the planets and the moon and TT-TDB chebychev polynomials at \url{http://www.imcce.fr/inpop}.

Tidal effects arise from differential and inelastic deformation of a planet by a perturbing body. The continuous action of tides modify the rotation of the planet together with its orbit until an equilibrium situation is reached. It is often believed that synchronous motion is the most probable outcome of the tidal evolution process, since synchronous rotation is observed for the majority of the satellites in the Solar System. However, in the 19th century, Schiaparelli also assumed synchronous motion for the rotations of Mercury and Venus, and was later shown to be wrong. Rather, for planets in eccentric orbits synchronous rotation is very unlikely. The rotation period and axial tilt of exoplanets is still unknown, but a large number of planets have been detected close to the parent star and should have evolved to a final equilibrium situation. Therefore, based on the Solar System well studied cases, we can make some predictions for exoplanets. Here we describe in detail the main tidal effects that modify the secular evolution of the spin and the orbit of a planet. We then apply our knowledge acquired from Solar System situations to exoplanet cases. In particular, we will focus on two classes of planets, "Hot-Jupiters" (fluid) and "Super-Earths" (rocky with atmosphere).

Since the original work of Hansen and Tisserand in the XIXth century, there have been many variations in the analytical expansion of the three-body disturbing function in series of the semi-major axis ratio. With the increasing number of planetary systems of large eccentricity, these expansions are even more interesting as they allow us to obtain for the secular systems finite expressions that are valid for all eccentricities and inclinations. We revisited the derivation of the disturbing function in Legendre polynomial, with a special focus on the secular system. We provide here expressions of the disturbing function for the planar and spatial case at any order with respect to the ratio of the semi-major axes. Moreover, for orders in the ratio of semi-major axis up to ten in the planar case and five in the spatial case, we provide explicit expansions of the secular system, and simple algorithms with minimal computation to extend this to higher order, as well as the algorithms for the computation of non secular terms.

We assess the ability of a solid ring to model a global perturbation induced by several thousands of main-belt asteroids. The ring is first studied in an analytical framework that provides an estimate of all the ring's parameters excepting mass. In the second part, numerically estimated perturbations on the Earth-Mars, Earth-Venus, and Earth-Mercury distances induced by various subsets of the main-belt population are compared with perturbations induced by a ring. To account for large uncertainties in the asteroid masses, we obtain results from Monte Carlo experiments based on asteroid masses randomly generated according to available data and the statistical asteroid model. The radius of the ring is analytically estimated at 2.8 AU. A systematic comparison of the ring with subsets of the main belt shows that, after removing the 300 most perturbing asteroids, the total main-belt perturbation of the Earth-Mars distance reaches on average 246 m on the 1969-2010 time interval. A ring with appropriate mass is able to reduce this effect to 38 m. We show that, by removing from the main belt ~240 asteroids that are not necessarily the most perturbing ones, the corresponding total perturbation reaches on average 472 m, but the ring is able to reduce it down to a few meters, thus accounting for more than 99% of the total effect.

The origin of the high inclination of Uranus' spin-axis (Uranus' obliquity) is one of the great unanswered questions about the Solar system. Giant planets are believed to form with nearly zero obliquity, and it has been shown that the present behaviour of Uranus' spin is essentially stable. Several attempts were made in order to solve this problem. Here we report numerical simulations showing that Uranus' axis can be tilted during the planetary migration, without the need of a giant impact, provided that the planet had an additional satellite and a temporary large inclination. This might have happened during the giant planet instability phase described in the Nice model. In our scenario, the satellite is ejected after the tilt by a close encounter at the end of the migration. This model can both explain Uranus' large obliquity and bring new constraints on the planet orbital evolution.

Precise radial-velocity measurements for data acquired with the HARPS spectrograph infer that three planets orbit the M4 dwarf star GJ876. In particular, we confirm the existence of planet "d", which orbits every 1.93785 days. We find that its orbit may have significant eccentricity (e=0.14), and deduce a more accurate estimate of its minimum mass of 6.3 Earth masses. Dynamical modeling of the HARPS measurements combined with literature velocities from the Keck Observatory strongly constrain the orbital inclinations of the "b" and "c" planets. We find that i_b = 48.9 degrees and i_c = 48.1 degrees, which infers the true planet masses of M_b = 2.64 Jupiter masses and M_c = 0.83 Jupiter masses, respectively. Radial velocities alone, in this favorable case, can therefore fully determine the orbital architecture of a multi-planet system, without the input from astrometry or transits. The orbits of the two giant planets are nearly coplanar, and their 2:1 mean motion resonance ensures stability over at least 5 Gyr. The libration amplitude is smaller than 2 degrees, suggesting that it was damped by some dissipative process during planet formation. The system has space for a stable fourth planet in a 4:1 mean motion resonance with planet "b", with a period around 15 days. The radial velocity measurements constrain the mass of this possible additional planet to be at most that of the Earth.

Long-term precise Doppler measurements with the CORALIE spectrograph revealed the presence of two massive companions to the solar-type star HD202206. Although the three-body fit of the system is unstable, it was shown that a 5:1 mean motion resonance exists close to the best fit, where the system is stable. We present here an extensive dynamical study of the HD202206 system aiming at constraining the inclinations of the two known companions, from which we derive possible ranges of value for the companion masses. We study the long term stability of the system in a small neighborhood of the best fit using Laskar's frequency map analysis. We also introduce a numerical method based on frequency analysis to determine the center of libration mode inside a mean motion resonance. We find that acceptable coplanar configurations are limited to inclinations to the line of sight between 30 and 90 degrees. This limits the masses of both companions to roughly twice the minimum. Non coplanar configurations are possible for a wide range of mutual inclinations from 0 to 90 degrees, although $\Delta\Omega = 0 [\pi]$ configurations seem to be favored. We also confirm the 5:1 mean motion resonance to be most likely. In the coplanar edge-on case, we provide a very good stable solution in the resonance, whose $\chi^2$ does not differ significantly from the best fit. Using our method to determine the center of libration, we further refine this solution to obtain an orbit with a very low amplitude of libration, as we expect dissipative effects to have dampened the libration.

The Lidov-Kozai mechanism allows a body to periodically exchange its eccentricity with inclination. It was first discussed in the framework of the quadrupolar secular restricted three-body problem, where the massless particle is the inner body, and later extended to the quadrupolar secular nonrestricted three body problem. In this paper, we propose a different point of view on the problem by looking first at the restricted problem where the massless particle is the outer body. In this situation, equilibria at high mutual inclination appear, which correspond to the population of stable particles that Verrier & Evans (2008,2009) find in stable, high inclination circumbinary orbits around one of the components of the quadruple star HD 98800. We provide a simple analytical framework using a vectorial formalism for these situations. We also look at the evolution of these high inclination equilibria in the non restricted case.

In this Letter, we give new constraints on planet migration. They were obtained under the assumption that Saturn's current obliquity is due to a capture in resonance with Neptune's ascending node. If planet migration is too fast, then Saturn crosses the resonance without being captured and it keeps a small obliquity. This scenario thus gives a lower limit on the migration time scale tau. We found that this boundary depends strongly on Neptune's initial inclination. For two different migration types, we found that tau should be at least greater than 7 Myr. This limit increases rapidly as Neptune's initial inclination decreases from 10 to 1 degree. We also give an algorithm to know if Saturn can be tilted for any migration law.

The present obliquity of Mercury is very low (less than 0.1 degree), which led previous studies to always adopt a nearly zero obliquity during the planet's past evolution. However, the initial orientation of Mercury's rotation axis is unknown and probably much different than today. As a consequence, we believe that the obliquity could have been significant when the rotation rate of the planet first encountered spin-orbit resonances. In order to compute the capture probabilities in resonance for any evolutionary scenario, we present in full detail the dynamical equations governing the long term evolution of the spin, including the obliquity contribution. The secular spin evolution of Mercury results from tidal interactions with the Sun, but also from viscous friction at the core-mantle boundary. Here, this effect is also regarded with particular attention. Previous studies show that a liquid core enhances drastically the chances of capture in spin-orbit resonances. We confirm these results for null obliquity, but we find that the capture probability generally decreases as the obliquity increases. We finally show that, when core-mantle friction is combined with obliquity evolution, the spin can evolve into some unexpected configurations as the synchronous or the 1/2 spin-orbit resonance.

The latest version of the planetary ephemerides developed at the Paris Observatory and at the Besancon Observatory is presented here. INPOP08 is a 4-dimension ephemeris since it provides to users positions and velocities of planets and the relation between TT and TDB. Investigations leading to improve the modeling of asteroids are described as well as the new sets of observations used for the fit of INPOP08. New observations provided by the European Space Agency (ESA) deduced from the tracking of the Mars Express (MEX) and Venus Express (VEX) missions are presented as well as the normal point deduced from the Cassini mission. We show the huge impact brought by these observations in the fit of INPOP08, especially in terms of Venus, Saturn and Earth-Moon barycenter orbits.

We consider the solid-solid interactions in the two body problem. The relative equilibria have been previously studied analytically and general motions were numerically analyzed using some expansion of the gravitational potential up to the second order, but only when there are no direct interactions between the orientation of the bodies. Here we expand the potential up to the fourth order and we show that the secular problem obtained after averaging over fast angles, as for the precession model of Boue and Laskar [Boue, G., Laskar, J., 2006. Icarus 185, 312-330], is integrable, but not trivially. We describe the general features of the motions and we provide explicit analytical approximations for the solutions. We demonstrate that the general solution of the secular system can be decomposed as a uniform precession around the total angular momentum and a periodic symmetric orbit in the precessing frame. More generally, we show that for a general n-body system of rigid bodies in gravitational interaction, the regular quasiperiodic solutions can be decomposed into a uniform precession around the total angular momentum, and a quasiperiodic motion with one frequency less in the precessing frame.

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.

Precise radial-velocity measurements with the HARPS spectrograph reveal the presence of two planets orbiting the solar-type star HD45364. The companion masses are 0.187 Mjup and 0.658 Mjup, with semi-major axes of 0.681 AU and 0.897 AU, and eccentricities of 0.168 and 0.097, respectively. A dynamical analysis of the system further shows a 3:2 mean motion resonance between the two planets, which prevents close encounters and ensures the stability of the system over 5 Gyr. This is the first time that such a resonant configuration has been observed for extra-solar planets, although there is an analogue in our Solar System formed by Neptune and Pluto. This singular planetary system may provide important constraints on planetary formation and migration scenarios.

The rotation of Mercury is presently captured in a 3/2 spin-orbit resonance with the orbital mean motion. The capture mechanism is well understood as the result of tidal interactions with the Sun combined with planetary perturbations. However, it is now almost certain that Mercury has a liquid core, which should induce a contribution of viscous friction at the core-mantle boundary to the spin evolution. This last effect greatly increases the chances of capture in all spin-orbit resonances, being 100% for the 2/1 resonance, and thus preventing the planet from evolving to the presently observed configuration. Here we show that for a given resonance, as the chaotic evolution of Mercury's orbit can drive its eccentricity to very low values during the planet's history, any previous capture can be destabilized whenever the eccentricity becomes lower than a critical value. In our numerical integrations of 1000 orbits of Mercury over 4 Gyr, the spin ends 99.8% of the time captured in a spin-orbit resonance, in particular in one of the following three configurations: 5/2 (22%), 2/1 (32%) and 3/2 (26%). Although the present 3/2 spin-orbit resonance is not the most probable outcome, we also show that the capture probability in this resonance can be increased up to 55% or 73%, if the eccentricity of Mercury in the past has descended below the critical values 0.025 or 0.005, respectively.

We summarise the scientific and technological aspects of the SAGAS (Search for Anomalous Gravitation using Atomic Sensors) project, submitted to ESA in June 2007 in response to the Cosmic Vision 2015-2025 call for proposals. The proposed mission aims at flying highly sensitive atomic sensors (optical clock, cold atom accelerometer, optical link) on a Solar System escape trajectory in the 2020 to 2030 time-frame. SAGAS has numerous science objectives in fundamental physics and Solar System science, for example numerous tests of general relativity and the exploration of the Kuiper belt. The combination of highly sensitive atomic sensors and of the laser link well adapted for large distances will allow measurements with unprecedented accuracy and on scales never reached before. We present the proposed mission in some detail, with particular emphasis on the science goals and associated measurements.

The equilibrium rotation of tidally evolved "Earth-like" extra-solar planets is often assumed to be synchronous with their orbital mean motion. The same assumption persisted for Mercury and Venus until radar observations revealed their true spin rates. As many of these planets follow eccentric orbits and are believed to host dense atmospheres, we expect the equilibrium rotation to differ from the synchronous motion. Here we provide a general description of the allowed final equilibrium rotation states of these planets, and apply this to already discovered cases in which the mass is lower than twelve Earth-masses. At low obliquity and moderate eccentricity, it is shown that there are at most four distinct equilibrium possibilities, one of which can be retrograde. Because most presently known "Earth-like" planets present eccentric orbits, their equilibrium rotation is unlikely to be synchronous.

A statistical analysis is performed over more than 1001 different integrations of the secular equations of the Solar system over 5 Gyr. With this secular system, the probability of the eccentricity of Mercury to reach 0.6 in 5 Gyr is about 1 to 2 %. In order to compare with (Ito and Tanikawa, 2002), we have performed the same analysis without general relativity, and obtained even more orbits of large eccentricity for Mercury. We have performed as well a direct integration of the planetary orbits, without averaging, for a dynamical model that do not include the Moon or general relativity with 10 very close initial conditions over 3 Gyr. The statistics obtained with this reduced set are comparable to the statistics of the secular equations, and in particular we obtain two trajectories for which the eccentricity of Mercury increases beyond 0.8 in less than 1.3 Gyr and 2.8 Gyr respectively. These strong instabilities in the orbital motion of Mecury results from secular resonance beween the perihelion of Jupiter and Mercury that are facilitated by the absence of general relativity. The statistical analysis of the 1001 orbits of the secular equations also provides probability density functions (PDF) for the eccentricity and inclination of the terrestrial planets.

Over the past two years, the search for low-mass extrasolar planets has led to the detection of seven so-called 'hot Neptunes' or 'super-Earths' around Sun-like stars. These planets have masses 5-20 times larger than the Earth and are mainly found on close-in orbits with periods of 2-15 days. Here we report a system of three Neptune-mass planets with periods of 8.67, 31.6 and 197 days, orbiting the nearby star HD 69830. This star was already known to show an infrared excess possibly caused by an asteroid belt within 1 AU (the Sun-Earth distance). Simulations show that the system is in a dynamically stable configuration. Theoretical calculations favour a mainly rocky composition for both inner planets, while the outer planet probably has a significant gaseous envelope surrounding its rocky/icy core; the outer planet orbits within the habitable zone of this star.

Eccentricity or obliquity tides have been proposed as the missing energy source that may explain the anomalously large radius of some transiting ``hot Jupiters''. To maintain a non-zero and large obliquity, it was argued that the planets can be locked in a Cassini state, i.e. a resonance between spin and orbital precessions. We compute the tidal heating within ``inflated'' close-in giant planets with a non-zero eccentricity or obliquity. We further inspect whether the spin of a ``hot Jupiter'' could have been trapped and maintained in a Cassini state during its early despinning and migration. We estimate the capture probability in a spin-orbit resonance between $\sim$ 0.5 AU (a distance where tidal effects become significant) and 0.05 AU for a wide range of secular orbital frequencies and amplitudes of gravitational perturbations. Numerical simulations of the spin evolution are performed to explore the influence of tidal despinning and migration processes on the resonance stability. We find that tidal heating within a non-synchronous giant planet is about twice larger than previous estimates based on the hypothesis of synchronization. Chances of capture in a spin-orbit resonance are very good around 0.5 AU but they decrease dramatically with the semi-major axis. Furthermore, even if captured, both tidal despinning and migration processes cause the tidal torque to become large enough that the obliquity ultimately leaves the resonance and switches to near $0^{\circ}$. Locking a ``hot Jupiter'' in an isolated spin-orbit resonance is unlikely at 0.05 AU but could be possible at larger distances. Another mechanism is then required to maintain a large obliquity and create internal heating through obliquity tides

The mu Ara planetary system is rather complex: It contains two already known planets, mu Ara b with P=640 days and mu Ara c with P=9.64 days, and a third companion on a wide but still poorly defined orbit. Even with three planets in the system, the data points keep anomalously high dispersion around the fitted solution. The high residuals are only partially due to the strong p-mode oscillations of the host star. We have therefore studied in this paper the possible presence of a fourth planet in the system. During the past years we have carried out additional and extremely precise radial-velocity measurements with the HARPS spectrograph. We provide in this paper a full orbital solution of the planetary system around mu Ara. It turns out to be the second system known to harbor 4 planetary companions. Thanks to the new data points acquired with HARPS we can confirm the presence of mu Ara c at P=9.64 days, which produces a coherent RV signal over more than two years. The new orbital fit sets the mass of mu Ara c to 10.5 M_Earth. Furthermore, we present the discovery of mu Ara d, a new planet on an almost circular 310 days-period and with a mass of 0.52 M_Jup. Finally, we give completely new orbital parameters for the longest-period planet, mu Ara e. It is the first time that this companion is constrained by radial-velocity data into a dynamical stable orbit, which leaves no doubt about its planetary nature. (Abridged).

Long-term precise Doppler measurements with the CORALIE spectrograph reveal the presence of a second planet orbiting the solar-type star HD202206. The radial-velocity combined fit yields companion masses of m_2\sini = 17.4 M_Jup and 2.44 M_Jup, semi-major axes of a = 0.83 AU and 2.55 AU, and eccentricities of e = 0.43 and 0.27, respectively. A dynamical analysis of the system further shows a 5/1 mean motion resonance between the two planets. This system is of particular interest since the inner planet is within the brown-dwarf limits while the outer one is much less massive. Therefore, either the inner planet formed simultaneously in the protoplanetary disk as a superplanet, or the outer Jupiter-like planet formed in a circumbinary disk. We believe this singular planetary system will provide important constraints on planetary formation and migration scenarios.

Recently, Breiter et al (2004) reported the computation of Hansen coefficients $X_k^{\gamma,m}$ for non integer values of $\gamma$. In fact, the Hansen coefficients are closely related to the Laplace $b_{s}^{(m)}$, and generalized Laplace coefficients $b_{s,r}^{(m)}$ (Laskar and Robutel, 1995) that do not require $s,r$ to be integers. In particular, the coefficients $X_0^{\g,m}$ have very simple expressions in terms of the usual Laplace coefficients $b_{\g+2}^{(m)}$, and all their properties derive easily from the known properties of the Laplace coefficients.

In a recent paper, Edvardsson etal (2002) propose a new solution for the spin evolution of the Earth and Mars. Their results differ significantly with respect to previous studies, as they found a large contribution on the precession of the planet axis from the tidal effects of Phobos and Deimos. In fact, this probably results from the omission by the authors of the torques exerted on the satellites orbits by the planet's equatorial bulge, as otherwise the average torque exerted by the satellites on the planet is null.

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.