The distribution of hot Jupiters, for which star-planet interactions can be significant, questions the evolution of exosystems. We aim to follow the orbital evolution of a planet along the rotational and structural evolution of the host star by taking into account the coupled effects of tidal and magnetic torques from ab initio prescriptions. It allows us to better understand the evolution of star-planet systems and to explain some properties of the distribution of observed close-in planets. To this end we use a numerical model of a coplanar circular star-planet system taking into account stellar structural changes, wind braking and star-planet interactions, called ESPEM (Benbakoura et al. 2019). We find that depending on the initial configuration of the system, magnetic effects can dominate tidal effects during the various phases of the evolution, leading to an important migration of the planet and to significant changes on the rotational evolution of the star. Both kinds of interactions thus have to be taken into account to predict the evolution of compact star-planet systems.

We present the scientific motivation for future space tests of the equivalence principle, and in particular the universality of free fall, at the $10^{-17}$ level or better. Two possible mission scenarios, one based on quantum technologies, the other on electrostatic accelerometers, that could reach that goal are briefly discussed.

We investigate the performance of the upcoming ACES (Atomic Clock Ensemble in Space) space mission in terms of its primary scientific objective, the test of the gravitational redshift. Whilst the ultimate performance of that test is determined by the systematic uncertainty of the on-board clock at 2-3 ppm, we determine whether, and under which conditions, that limit can be reached in the presence of colored realistic noise, data gaps and orbit determination uncertainties. To do so we have developed several methods and software tools to simulate and analyse ACES data. Using those we find that the target uncertainty of 2-3 ppm can be reached after only a few measurement sessions of 10-20 days each, with a relatively modest requirement on orbit determination of around 300 m.

With the discovery over the last two decades of a large diversity of exoplanetary systems, it is now of prime importance to characterize star-planet interactions and how such systems evolve. We address this question by studying systems formed by a solar-like star and a close-in planet. We focus on the stellar wind spinning down the star along its main sequence phase and tidal interaction causing orbital evolution of the systems. Despite recent significant advances in these fields, all current models use parametric descriptions to study at least one of these effects. Our objective is to introduce simultaneously ab-initio prescriptions of the tidal and braking torques, so as to improve our understanding of the underlying physics. We develop a 1D numerical model of coplanar circular star-planet systems taking into account stellar structural changes, wind braking and tidal interaction and implement it in a code called ESPEM. We follow the secular evolution of the stellar rotation assuming a bi-layer internal structure, and of the semi-major axis of the orbit. After comparing our predictions to recent observations and models, we perform tests to emphasize the contribution of ab-initio prescriptions. Our secular model of stellar wind braking reproduces well the recent observations of stellar rotation in open clusters. Our results show that a planet can affect the rotation of its host star and that the resulting spin-up or spin-down depends on the orbital semi-major axis and on the joint influence of magnetic and tidal effects. The ab-initio prescription for tidal dissipation that we used predicts fast outward migration of massive planet orbiting fast-rotating young stars. Finally, we provide the reader with a criterion based on the system's characteristics that allows us to assess whether or not the planet will undergo orbital decay due to tidal interaction.

In this proceeding, we show how observations of Solar System Objects with Gaia can be used to test General Relativity and to constrain modified gravitational theories. The high number of Solar System objects observed and the variety of their orbital parameters associated with the impressive astrometric accuracy will allow us to perform local tests of General Relativity. In this communication, we present a preliminary sensitivity study of the Gaia observations on dynamical parameters such as the Sun quadrupolar moment and on various extensions to general relativity such as the parametrized post-Newtonian parameters, the fifth force formalism and a violation of Lorentz symmetry parametrized by the Standard-Model extension framework. We take into account the time sequences and the geometry of the observations that are particular to Gaia for its nominal mission (5 years) and for an extended mission (10 years).

Accurate tracking of probes is one of the key points of space exploration. Range and Doppler techniques are the most commonly used. In this paper we analyze the impact of the transponder delay, $i.e.$ the processing time between reception and re-emission of a two-way tracking link at the satellite, on tracking observables and on spacecraft orbits. We show that this term, only partially accounted for in the standard formulation of computed space observables, can actually be relevant for future missions with high nominal tracking accuracies or for the re-processing of old missions. We present several applications of our formulation to Earth flybys, the NASA GRAIL and the ESA BepiColombo missions.

A key objective of the ESA Gaia satellite is the realization of a quasi-inertial reference frame at visual wavelengths by means of global astrometric techniques. This requires an accurate mathematical and numerical modeling of relativistic light propagation, as well as double-blind-like procedures for the internal validation of the results, before they are released to the scientific community at large. Aim of this work is to specialize the Time Transfer Functions (TTF) formalism to the case of the Gaia observer and prove its applicability to the task of Global Sphere Reconstruction (GSR), in anticipation of its inclusion in the GSR system, already featuring the suite of RAMOD models, as an additional semi-external validation of the forthcoming Gaia baseline astrometric solutions. We extend the current GSR framework and software infrastructure (GSR2) to include TTF relativistic observation equations compatible with Gaia's operations. We use simulated data generated by the Gaia Data Reduction and Analysis Consortium (DPAC) to obtain different least-squares estimations of the full stellar spheres and gauge results. These are compared to analogous solutions obtained with the current RAMOD model in GSR2 and to the catalog generated with GREM, the model baselined for Gaia and used to generate the DPAC synthetic data. Linearized least-squares TTF solutions are based on spheres of about 132,000 primary stars uniformly distributed on the sky and simulated observations spanning the entire 5-yr range of Gaia's nominal operational lifetime. The statistical properties of the results compare well with those of GREM. Finally, comparisons to RAMOD@GSR2 solutions confirmed the known lower accuracy of that model and allowed us to establish firm limits on the quality of the linearization point outside of which an iteration for non-linearity is required for its proper convergence.

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.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Lorentz symmetry violations can be described by an effective field theory framework that contains both General Relativity and the Standard Model of particle physics called the Standard-Model extension (SME). Recently, post-fit analysis of Gravity Probe B and binary pulsars lead to an upper limit at the $10^{-4}$ level on the time-time coefficient $\bar s^{TT}$ of the pure-gravity sector of the minimal SME. In this work, we derive the observable of Very Long Baseline Interferometry (VLBI) in SME and then we implement it into a real data analysis code of geodetic VLBI observations. Analyzing all available observations recorded since 1979, we compare estimates of $\bar s^{TT}$ and errors obtained with various analysis schemes, including global estimations over several time spans and with various Sun elongation cut-off angles, and with analysis of radio source coordinate time series. We obtain a constraint on $\bar s^{TT}=(-5\pm 8)\times 10^{-5}$, directly fitted to the observations and improving by a factor 5 previous post-fit analysis estimates.

Lorentz symmetry violations can be parametrized by an effective field theory framework that contains both General Relativity and the Standard Model of particle physics, called the Standard-Model Extension or SME. We consider in this work only the pure gravitational sector of the minimal SME. We present new constraints on the SME coefficients obtained from lunar laser ranging, very long baseline interferometry, and planetary motions.

Turbulent friction in convective regions in stars and planets is one of the key physical mechanisms that drive the dissipation of the kinetic energy of tidal flows in their interiors and the evolution of their systems. This friction acts both on the equilibrium/non-wave like tide and on tidal inertial waves in these layers. It is thus necessary to obtain a robust prescription for this friction. In the current state-of-the-art, it is modeled by a turbulent eddy-viscosity coefficient, based on mixing-length theory, applied on velocities of tides. However, none of the current prescriptions take into account the action of rotation that can strongly affects turbulent convection. Therefore, we use theoretical scaling laws for convective velocities and characteristic lengthscales in rotating stars and planets that have been recently confirmed by 3-D high-resolution nonlinear Cartesian numerical simulations to derive a new prescription. A corresponding local model of tidal waves is used to understand the consequences for the linear tidal dissipation. Finally, new grids of rotating stellar models and published values of planetary convective Rossby numbers are used to discuss astrophysical consequences. The action of rotation on convection deeply modifies the turbulent friction applied on tides. In the regime of rapid rotation (with a convective Rossby number below 0.25), the eddy-viscosity may be decreased by several ordres of magnitude. It may lead to a loss of efficiency of the viscous dissipation of the equilibrium tide and to a more efficient complex and resonant dissipation of tidal inertial waves in the bulk of convective regions. Therefore, it is necessary to have a completely coupled treatment of the tidal/rotational evolution of star-planet systems and multiple stars with a coherent treatment of the variations of tidal flows and of their dissipation as a function of rotation.

Using astrometric observations spanning more than a century and including a large set of Cassini data, we determine Saturn's tidal parameters through their current effects on the orbits of the eight main and four coorbital moons. We have used the latter to make the first determination of Saturn's Love number, $k_2=0.390 \pm 0.024$, a value larger than the commonly used theoretical value of 0.341 (Gavrilov & Zharkov, 1977), but compatible with more recent models (Helled & Guillot, 2013) for which $k_2$ ranges from 0.355 to 0.382. Depending on the assumed spin for Saturn's interior, the new constraint can lead to a reduction of up to 80% in the number of potential models, offering great opportunities to probe the planet's interior. In addition, significant tidal dissipation within Saturn is confirmed (Lainey et al., 2012) corresponding to a high present-day tidal ratio $k_2/Q=(1.59 \pm 0.74) \times 10^{-4}$ and implying fast orbital expansions of the moons. This high dissipation, with no obvious variations for tidal frequencies corresponding to those of Enceladus and Dione, may be explained by viscous friction in a solid core, implying a core viscosity typically ranging between $10^{14}$ and $10^{16}$ Pa.s (Remus et al., 2012). However, a dissipation increase by one order of magnitude at Rhea's frequency could suggest the existence of an additional, frequency-dependent, dissipation process, possibly from turbulent friction acting on tidal waves in the fluid envelope of Saturn (Ogilvie & Li, 2004). Alternatively, a few of Saturn's moons might themselves experience large tidal dissipation.

In this communication, we show how asteroids observations from the Gaia mission can be used to perform local tests of General Relativity (GR). This ESA mission, launched in December 2013, will observe --in addition to the stars-- a large number of small Solar System Objects (SSOs) with unprecedented astrometric precision. Indeed, it is expected that about 360,000 asteroids will be observed with a nominal sub-mas precision. Here, we show how these observations can be used to constrain some extensions to General Relativity. We present results of SSOs simulations that take into account the time sequences over 5 years and geometry of the observations that are particular to Gaia. We present a sensitivity study on various GR extensions and dynamical parameters including: the Sun quadrupolar moment $J_2$, the parametrized post-Newtonian parameter $\beta$, the Nordtvedt parameter $\eta$, the fifth force formalism, the Lense-Thirring effect, a temporal variation of the gravitational parameter $GM_\textrm{sun}$ (a linear variation as well as a periodic variation), the Standard Model Extension formalism,... Some implications for planetary ephemerides analysis are also briefly discussed.

Planetary systems evolve over secular time scales. One of the key mechanisms that drive this evolution is tidal dissipation. Submitted to tides, stellar and planetary fluid layers do not behave like rocky ones. Indeed, they are the place of resonant gravito-inertial waves. Therefore, tidal dissipation in fluid bodies strongly depends on the excitation frequency while this dependence is smooth in solid ones. Thus, the impact of the internal structure of celestial bodies must be taken into account when studying tidal dynamics. The purpose of this work is to present a local model of tidal gravito-inertial waves allowing us to quantify analytically the internal dissipation due to viscous friction and thermal diffusion, and to study the properties of the resonant frequency spectrum of the dissipated energy. We derive from this model scaling laws characterizing tidal dissipation as a function of fluid parameters (rotation, stratification, diffusivities) and discuss them in the context of star-planet systems.

Planetary ephemerides are a very powerful tool to constrain deviations from the theory of General Relativity using orbital dynamics. The effective field theory framework called the Standard-Model Extension (SME) has been developed in order to systematically parametrize hypothetical violations of Lorentz symmetry (in the Standard Model and in the gravitational sector). In this communication, we use the latest determinations of the supplementary advances of the perihelia and of the nodes obtained by planetary ephemerides analysis to constrain SME coefficients from the pure gravity sector and also from gravity-matter couplings. Our results do not show any deviation from GR and they improve current constraints. Moreover, combinations with existing constraints from Lunar Laser Ranging and from atom interferometry gravimetry allow us to disentangle contributions from the pure gravity sector from the gravity-matter couplings.

Tidal dissipation in planets and stars is one of the key physical mechanisms driving the evolution of star-planet and planet-moon systems. Several signatures of its action are observed in planetary systems thanks to their orbital architecture and the rotational state of their components. Tidal dissipation inside the fluid layers of celestial bodies are intrinsically linked to the dynamics and the physical properties of the latter. This complex dependence must be characterized. We compute the tidal kinetic energy dissipated by viscous friction and thermal diffusion in a rotating local fluid Cartesian section of a star/planet/moon submitted to a periodic tidal forcing. The properties of tidal gravito-inertial waves excited by the perturbation are derived analytically as explicit functions of the tidal frequency and local fluid parameters (i.e. the rotation, the buoyancy frequency characterizing the entropy stratification, viscous and thermal diffusivities) for periodic normal modes. The sensitivity of the resulting possibly highly resonant dissipation frequency-spectra to a control parameter of the system is either important or negligible depending on the position in the regime diagram relevant for planetary and stellar interiors. For corresponding asymptotic behaviors of tidal gravito-inertial waves dissipated by viscous friction and thermal diffusion, scaling laws for the frequencies, number, width, height and contrast with the non-resonant background of resonances are derived to quantify these variations. We characterize the strong impact of the internal physics and dynamics of fluid planetary layers and stars on the dissipation of tidal kinetic energy in their bulk. We point out the key control parameters that really play a role and demonstrate how it is now necessary to develop ab-initio modeling for tidal dissipation in celestial bodies.

Determining range, Doppler and astrometric observables is of crucial interest for modelling and analyzing space observations. We recall how these observables can be computed when the travel time of a light ray is known as a function of the positions of the emitter and the receiver for a given instant of reception (or emission). For a long time, such a function--called a reception (or emission) time transfer function--has been almost exclusively calculated by integrating the null geodesic equations describing the light rays. However, other methods avoiding such an integration have been considerably developped in the last twelve years. We give a survey of the analytical results obtained with these new methods up to the third order in the gravitational constant $G$ for a mass monopole. We briefly discuss the case of quasi-conjunctions, where higher-order enhanced terms must be taken into account for correctly calculating the effects. We summarize the results obtained at the first order in $G$ when the multipole structure and the motion of an axisymmetric body is taken into account. We present some applications to on-going or future missions like Gaia and Juno. We give a short review of the recent works devoted to the numerical estimates of the time transfer functions and their derivatives.

Tidal dissipation is known as one of the main drivers of the secular evolution of planetary systems. It directly results from dissipative mechanisms that occur in planets and stars' interiors and strongly depends on the structure and dynamics of the bodies. This work focuses on the mechanism of viscous friction in stars and planetary layers. A local model is used to study tidal dissipation. It provides general scaling laws that give a qualitative overview of the different possible behaviors of fluid tidal waves. Furthermore, it highlights the sensitivity of dissipation to the tidal frequency and the roles played by the internal parameters of the fluid such as rotation, stratification, viscosity and thermal diffusivity that will impact the spins/orbital architecture in planetary systems.

Tidal dissipation in stars is one of the key physical mechanisms that drive the evolution of binary and multiple stars. As in the Earth oceans, it corresponds to the resonant excitation of their eigenmodes of oscillation and their damping. Therefore, it strongly depends on the internal structure, rotation, and dissipative mechanisms in each component. In this work, we present a local analytical modeling of tidal gravito-inertial waves excited in stellar convective and radiative regions respectively. This model allows us to understand in details the properties of the resonant tidal dissipation as a function of the excitation frequencies, the rotation, the stratification, and the viscous and thermal properties of the studied fluid regions. Then, the frequencies, height, width at half-height, and number of resonances as well as the non-resonant equilibrium tide are derived analytically in asymptotic regimes that are relevant in stellar interiors. Finally, we demonstrate how viscous dissipation of tidal waves leads to a strongly erratic orbital evolution in the case of a coplanar binary system. We characterize such a non-regular dynamics as a function of the height and width of resonances, which have been previously characterized thanks to our local fluid model.

Given the extreme accuracy of modern space science, a precise relativistic modeling of observations is required. We use the Time Transfer Functions formalism to study light propagation in the field of uniformly moving axisymmetric bodies, which extends the field of application of previous works. We first present a space-time metric adapted to describe the geometry of an ensemble of uniformly moving bodies. Then, we show that the expression of the Time Transfer Functions in the field of a uniformly moving body can be easily derived from its well-known expression in a stationary field by using a change of variables. We also give a general expression of the Time Transfer Function in the case of an ensemble of arbitrarily moving point masses. This result is given in the form of an integral easily computable numerically. We also provide the derivatives of the Time Transfer Function in this case, which are mandatory to compute Doppler and astrometric observables. We particularize our results in the case of moving axisymmetric bodies. Finally, we apply our results to study the different relativistic contributions to the range and Doppler tracking for the JUNO mission in the Jovian system.

During the last decade, the precision in the tracking of spacecraft has constantly improved. The discovery of few astrometric anomalies, such as the Pioneer and Earth flyby anomalies, stimulated further analysis of the operative modeling currently adopted in Deep Space Navigation (DSN). Our study shows that some traditional approximations lead to neglect tiny terms that could have consequences in the orbit determination of a probe in specific configurations such as during an Earth flyby. Therefore, we suggest here a way to improve the light time calculation used for probe tracking.

In this communication, the current tests of gravitation available at Solar System scales are recalled. These tests rely mainly on two frameworks: the PPN framework and the search for a fifth force. Some motivations are given to look for deviations from General Relativity in other frameworks than the two extensively considered. A recent analysis of Cassini data in a MOND framework is presented. Furthermore, possibilities to constrain Standard Model Extension parameters using Solar System data are developed.

Given the extreme accuracy of modern space science, a precise relativistic modeling of observations is required. In particular, it is important to describe properly light propagation through the Solar System. For two decades, several modeling efforts based on the solution of the null geodesic equations have been proposed but they are mainly valid only for the first order Post-Newtonian approximation. However, with the increasing precision of ongoing space missions as Gaia, GAME, BepiColombo, JUNO or JUICE, we know that some corrections up to the second order have to be taken into account for future experiments. We present a procedure to compute the relativistic coordinate time delay, Doppler and astrometric observables avoiding the integration of the null geodesic equation. This is possible using the Time Transfer Function formalism, a powerful tool providing key quantities such as the time of flight of a light signal between two point-events and the tangent vector to its null-geodesic. Indeed we show how to compute the Time Transfer Functions and their derivatives (and thus range, Doppler and astrometric observables) up to the second post-Minkowskian order. We express these quantities as quadratures of some functions that depend only on the metric and its derivatives evaluated along a Minkowskian straight line. This method is particularly well adapted for numerical estimations. As an illustration, we provide explicit expressions in static and spherically symmetric space-time up to second post-Minkowskian order. Then we give the order of magnitude of these corrections for the range/Doppler on the BepiColombo mission and for astrometry in a GAME-like observation.

Given the extreme accuracy of modern space astrometry, a precise relativistic modeling of observations is required. Moreover, the availability of several models formulated in different and independent ways is a security against the presence of systematic errors in the analysis of future experimental results, like in the case of the Gaia mission. In this work, we simulate a series of observations using the two models to be used for the data analysis of Gaia, the Gaia RElativistic Model (GREM) and the Relativistic Astrometric MODel (RAMOD), and we compare them with the results of our astrometric model based on the Time Transfer Functions.

Context. Tidal dissipation in planets and in stars is one of the key physical mechanisms that drive the evolution of planetary systems. Aims. Tidal dissipation properties are intrisically linked to the internal structure and the rheology of studied celestial bodies. The resulting dependence of the dissipation upon the tidal frequency is strongly different in the cases of solids and fluids. Methods. We compute the tidal evolution of a two-body coplanar system, using the tidal quality factor's frequency-dependencies appropriate to rocks and to convective fluids. Results. The ensuing orbital dynamics comes out smooth or strongly erratic, dependent on how the tidal dissipation depends upon frequency. Conclusions. We demonstrate the strong impact of the internal structure and of the rheology of the central body on the orbital evolution of the tidal perturber. A smooth frequency-dependence of the tidal dissipation renders a smooth orbital evolution while a peaked dissipation can furnish erratic orbital behaviour.

In this communication, we focus on possibilities to constrain SME coefficients using Cassini and Messenger data. We present simulations of radioscience observables within the framework of the SME, identify the linear combinations of SME coefficients the observations depend on and determine the sensitivity of these measurements to the SME coefficients. We show that these datasets are very powerful for constraining SME coefficients.

Given the extreme accuracy reached by future global space astrometry, one needs a global relativistic modeling of observations. A relativistic definition of astrometric observables is then essential to find uniquely coordinates, parallax and proper motion of a stellar object. The standard procedure is to explicitly solve the null geodesic equations to describe the trajectory of a photon emitted by a celestial object and received by a moving observing satellite. However, this task can be avoided if one builds an astrometric set up using the recent formalism of time transfer functions.

In this paper, we focus on the possibility to test General Relativity in the Solar System with radioscience measurements. To this aim, we present a new software that simulates Range and Doppler signals directly from the space-time metric. This flexible approach allows one to perform simulations in General Relativity and in alternative metric theories of gravity. In a second step, a least-squares fit of the different initial conditions involved in the situation is performed in order to compare anomalous signals produced by a given alternative theory with the ones obtained in General Relativity. This software provides orders of magnitude and signatures stemming from hypothetical alternative theories of gravity on radioscience signals. As an application, we present some simulations done for the Cassini mission in Post-Einsteinian Gravity and in the context of MOND External Field Effect. We deduce constraints on the Post-Einsteinian parameters but find that the considered arc of the Cassini mission is not useful to constrain the MOND External Field Effect.

A lot of fundamental tests of gravitational theories rely on highly precise measurements of the travel time and/or the frequency shift of electromagnetic signals propagating through the gravitational field of the Solar System. In practically all of the previous studies, the explicit expressions of such travel times and frequency shifts as predicted by various metric theories of gravity are derived from an integration of the null geodesic differential equations. However, the solution of the geodesic equations requires heavy calculations when one has to take into account the presence of mass multipoles in the gravitational field or the tidal effects due to the planetary motions, and the calculations become quite complicated in the post-post-Minkowskian approximation. This difficult task can be avoided using the time transfer function's formalism. We present here our last advances in the formulation of the one-way frequency shift using this formalism up to the post-post-Minkowskian approximation.

Tidal interactions between Saturn and its satellites play a crucial role in both the orbital migration of the satellites and the heating of their interiors. Therefore constraining the tidal dissipation of Saturn (here the ratio k2/Q) opens the door to the past evolution of the whole system. If Saturn's tidal ratio can be determined at different frequencies, it may also be possible to constrain the giant planet's interior structure, which is still uncertain. Here, we try to determine Saturn's tidal ratio through its current effect on the orbits of the main moons, using astrometric data spanning more than a century. We find an intense tidal dissipation (k2/Q= (2.3 \pm 0.7) \times 10-4), which is about ten times higher than the usual value estimated from theoretical arguments. As a consequence, eccentricity equilibrium for Enceladus can now account for the huge heat emitted from Enceladus' south pole. Moreover, the measured k2/Q is found to be poorly sensitive to the tidal frequency, on the short frequency interval considered. This suggests that Saturn's dissipation may not be controlled by turbulent friction in the fluid envelope as commonly believed. If correct, the large tidal expansion of the moon orbits due to this strong Saturnian dissipation would be inconsistent with the moon formations 4.5 Byr ago above the synchronous orbit in the Saturnian subnebulae. But it would be compatible with a new model of satellite formation in which the Saturnian satellites formed possibly over longer time scale at the outer edge of the main rings. In an attempt to take into account for possible significant torques exerted by the rings on Mimas, we fitted a constant rate da/dt on Mimas semi-major axis, also. We obtained an unexpected large acceleration related to a negative value of da/dt= -(15.7 \pm 4.4) \times 10-15 au/day.

The laws of gravitation have been tested for a long time with steadily improving precision, leading at some moment of time to paradigmatic evolutions. Pursuing this continual effort is of great importance for science. In this communication, we focus on Solar System tests of gravity and more precisely on possible tests that can be performed with radio science observations (Range and Doppler). After briefly reviewing the current tests of gravitation at Solar System scales, we give motivations to continue such experiments. In order to obtain signature and estimate the amplitude of anomalous signals that could show up in radio science observables because of modified gravitational laws, we developed a new software that simulates Range/Doppler signals. We present this new tool that simulates radio science observables directly from the space-time metric. We apply this tool to the Cassini mission during its cruise from Jupiter to Saturn and derive constraints on the parameters entering alternative theories of gravity beyond the standard Parametrized Post Newtonian theory.

The origin of Saturn's inner mid-sized moons (Mimas, Enceladus, Tethys, Dione and Rhea) and Saturn's rings is debated. Charnoz et al. (2010) introduced the idea that the smallest inner moons could form from the spreading of the rings' edge while Salmon et al. (2010) showed that the rings could have been initially massive, and so was the ring's progenitor itself. One may wonder if the mid-sized moons may have formed also from the debris of a massive ring progenitor, as also suggested in Canup (2010). However, the process driving mid-sized moons accretion from the icy debris disks has not been investigated in details. In particular, this process does not seem able to explain the varying silicate contents of the mid-sized moons (from 6% to 57% in mass). Here, we explore the formation of large objects from a massive ice-rich ring (a few times Rhea's mass) and describe the fundamental properties and implications of this new process. Using a hybrid computer model, we show that accretion within massive icy rings can form all mid-sized moons from Mimas to Rhea. However in order to explain their current locations, intense dissipation within Saturn (with Qp<2000) would be required. Our results are consistent with a satellite origin tied to the rings formation at least 2.5 Gy ago, both compatible with either a formation concurrent to Saturn or during the Late Heavy Bombardment. Tidal heating related to high-eccentricity post-accretional episodes may induce early geological activity. If some massive irregular chunks of silicates were initially present within the rings, they would be present today inside the satellites' cores which would have accreted icy shells while being tidally expelled from the rings (via a heterogeneous accretion process)while those still present in the rings are interpreted as today Saturn's rings' propellers and ring-moons (like Pan or Daphnis).

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.

With the discovery during the past decade of a large number of extrasolar planets orbiting their parent stars at a distance lower than 0.1 astronomical unit (and the launch and the preparation of dedicated space missions such as CoRoT and KEPLER), with the position of inner natural satellites around giant planets in our Solar System and with the existence of very closed but separated binary stars, tidal interaction has to be carefully studied. In particular, a question arises about the validity of usual approximations used in the modelling of this interaction. The purpose of this paper is to examine the step beyond the ponctual approximation for the tidal perturber. To achieve this aim, the gravitational interaction between two extended bodies and more precisely the interaction between mass multipole moments of their gravitational fields and the associated tidal phenomena are studied. Use of Cartesian Symmetric Trace Free (STF) tensors, of their relation with spherical harmonics and of the Kaula's transform enables to derive analytically the tidal and mutual interaction potentials as well as the associated disturbing functions in extended bodies systems. The tidal and mutual interaction potentials of two extended bodies are derived. Next, the external gravitational potential of such tidally disturbed extended body is obtained, using the Love's number theory, as well as the associated disturbing function. Finally, the dynamical evolution equations for such systems are given in their more general form without any linearization. The dynamical equations for the gravitational and tidal interactions between extended bodies and associated dynamics are derived in a form where they could be directly implemented to perform coherent numerical simulations of planetary systems or multiple stars tidal evolution.

With the discoveries of very close star-planet systems with planet orbiting sometimes at several star radius but also with well-known situations in our solar system where natural satellites are very close to their parent planet the validity of the ponctual mass approximation for the tidal perturber (respectively the parent star or planet when we study the close planet or natural satellite dynamics) has to be examined. In this short paper, we consider this problematic using results coming from a complete formalism that allows to treat the tidal interaction between extended bodies. We focus on its application to a simplified configuration.

Future space astrometry missions are planned to measure positions and/or parallaxes of celestial objects with an accuracy of the order of the microarcsecond. At such a level of accuracy, it will be indispensable to take into account the influence of the mass multipole structure of the giant planets on the bending of light rays. Within the parametrized post-Newtonian formalism, we present an algorithmic procedure enabling to determine explicitly this influence on a light ray connecting two points located at a finite distance. Then we specialize our formulae in the cases where 1) the light source is located at space infinity, 2) both the light source and the observer are located at space infinity. We examine in detail the cases where the unperturbed ray is in the equatorial plane or in a meridian plane.

The International Celestial Reference Frame (ICRF, Ma et al. 1998) is currently the best realization of a quasi-inertial reference system. It is based on more than 10 years of cumulated geodetic and astrometric VLBI observations of compact extragalactic objects at centimetric wavelengths. In the perspective of the realization of an accurate optical counterpart of the ICRF using future space astrometry missions like GAIA or SIM, this paper investigates the consistency of celestial reference frames realized through the same subset of compact extragalactic radio sources at optical wavelengths. Celestial reference frames realized in radio wavelengths with the VLBA Calibrator Survey (VCS) data and in optical wavelengths with the Sloan Digital Sky Survey (SDSS) data (DR3 quasar catalogue and DR5) are compared in terms of radio-optical distances between the common sources, global rotation of the axes and offset of the equator.

The new generation of atomic clocks will reach unprecedented uncertainties in frequency of $10^{-18}$. In order to prepare space missions such as ACES, we compute all relativistic frequency shifts detectable during this mission in the case of a clock aboard the International Space Station.

Almost all of the studies devoted to relativistic astrometry are based on the integration of the null geodesic differential equations. However, the gravitational deflection of light rays can be calculated by a different method, based on the determination of the bifunction giving half the squared geodesic distance between two arbitrary points-events, the so-called Synge's world function. We give a brief review of the main results obtained by this method.

In suitably chosen domains of space-time, the world function may be a powerful tool for modelling the deflection of light and the time/frequency transfer. In this paper we work out a recursive procedure for expanding the world function into a perturbative series of ascending powers of the Newtonian gravitational constant G. We show rigorously that each perturbation term is given by a line integral taken along the unperturbed geodesic between two points. Once the world function is known, it becomes possible to determine the time transfer functions giving the propagation time of a photon between its emission and its reception. We establish that the direction of a light ray as measured in the 3-space relative to a given observer can be derived from these time transfer functions, even if the metric is not stationary. We show how to deduce these functions up to any given order in G from the perturbative expansion of the world function. To illustrate the method, we carry out the calculation of the world function and of the time transfer function outside a static, spherically symmetric body up to the order G^2, the metric containing three arbitrary parameters.