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 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.

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.

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} $.

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.

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.

The IAU Commission 4 Working Group on Standardizing Access to Ephemerides recommends the use of the Spacecraft and Planet Kernel (SPK) format as a standard format for the position ephemerides of planets and other natural solar system bodies, and the use of the Planetary Constants Kernel (PCK) format for the orientation of these bodies. It further recommends that other supporting data be stored in a text PCK. These formats were developed for use by the SPICE Toolkit by the Navigation and Ancillary Information Facility of NASA's Jet Propulsion Laboratory (JPL). The CALCEPH library developed by the Institut de mecanique celeste de calcul des ephemerides (IMCCE) is also able to make use of these files. High accuracy ephemerides available in files conforming to the SPK and PCK formats include: the Development Ephemerides (DE) from JPL, Integrateur Numerique Planetaire de l'Observatoire de Paris (INPOP) from IMCCE, and the Ephemerides Planets and the Moon (EPM), developed by the Institute for Applied Astronomy (IAA). The bulk of this report is a description of the portion of PCK and SPK formats required for these ephemerides. New SPK and PCK data types, both called Type 20: Chebyshev (Velocity Only), have been added. Other changes to the specification are (i) a new object identification number for coordinate time ephemerides and (ii) a set of three new data types that use the TCB rather than the TDB time scale for the ephemerides, but are otherwise identical to their TDB versions.

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.

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.

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.

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

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.

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}.

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.