With an estimated diameter in the 320 to 350 km range, (704) Interamnia is the fifth largest main belt asteroid and one of the few bodies that fills the gap in size between the four largest bodies with $D$ > 400 km (Ceres, Vesta, Pallas and Hygiea) and the numerous smaller bodies with $D$ $\lesssim$ 200 km. However, despite its large size, little is known about the shape and spin state of Interamnia and, therefore, about its bulk composition and past collisional evolution. We aimed to test at what size and mass the shape of a small body departs from a nearly ellipsoidal equilibrium shape (as observed in the case of the four largest asteroids) to an irregular shape as routinely observed in the case of smaller ($D$ $\lesssim$ 200 km) bodies. We observed Interamnia as part of our ESO VLT/SPHERE large program (ID: 199.C-0074) at thirteen different epochs. In addition, several new optical lightcurves were recorded. These data, along with stellar occultation data from the literature, were fed to the All-Data Asteroid Modeling (ADAM) algorithm to reconstruct the 3D-shape model of Interamnia and to determine its spin state. Interamnia's volume-equivalent diameter of 332 $\pm$ 6 km implies a bulk density of $\rho$=1.98 $\pm$ 0.68 gcm$^{-3}$ , which suggests that Interamnia - like Ceres and Hygiea - contains a high fraction of water ice, consistent with the paucity of apparent craters. Our observations reveal a shape that can be well approximated by an ellipsoid, and that is compatible with a fluid hydrostatic equilibrium at the 2 $\sigma$ level. The rather regular shape of Interamnia implies that the size and mass limit, under which the shapes of minor bodies with a high amount of water ice in the subsurface become irregular, has to be searched among smaller ($D$ $\lesssim$ 300km) less massive ($m$ $\lesssim$ 3x10$^{19}$ kg) bodies.

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.

During the thirteen years in orbit around Saturn before its final plunge, the Cassini spacecraft provided more than ten thousand astrometric measurements. Such large amounts of accurate data enable the search for extremely faint signals in the orbital motion of the saturnian moons. Among these, the detection of the dynamical feedback of the rotation of the inner moons of Saturn on their respective orbits becomes possible. Using all the currently available astrometric data associated with Atlas, Prometheus, Pandora, Janus and Epimetheus, we first provide a detailed analysis of the Cassini Imaging Science Subsystem (ISS) data, with special emphasis on their statistical behavior and sources of bias. Then, we give updated estimates of the moons' averaged densities and try to infer more details about their interior properties by estimating the physical librations for Prometheus, Pandora, Epimetheus and Janus from anomalies in their apsidal precession. Our results are compatible with a homogeneous interior for Janus and Epimetheus, within the uncertainty of the measurements. On the other hand, we found some inconsistency for Pandora and Prometheus, which might result from a dynamical mismodeling of Saturn's gravity field. Last, we show how the synergistic introduction of libration measurements directly derived from imaging should allow the moons' moments of inertia to be better constrained.

Natural satellite systems present a large variety of orbital configurations in the solar system. While some are clearly the result of known processes, others still have largely unexplained eccentricity and inclination values. Iapetus has a still unexplained 3% orbital eccentricity and its orbital plane is tilted with respect to its local Laplace plane. On the other hand, astrometric measurements of saturnian moons have revealed high tidal migration rates, corresponding to a quality factor Q of Saturn of around 1600 for the mid-sized icy moons. We show how a past crossing of the 5:1 mean motion resonance between Titan and Iapetus may be a plausible scenario to explain Iapetus' orbit. We have carried out numerical simulations of the resonance crossing using an N-Body code as well as using averaged equations of motion. A large span of migration rates were explored for Titan and Iapetus was started on its local Laplace plane with a circular orbit. The resonance crossing can trigger a chaotic evolution of the eccentricity and the inclination of Iapetus. The outcome of the resonance is highly dependent on the migration rate. For a quality factor Q of over around 2000, the chaotic evolution of Iapetus in the resonance leads in most cases to its ejection, while simulations with a quality factor between 100 and 2000 show a departure from the resonance with post-resonant eccentricities spanning from 0 up to 15%, and free inclinations capable of reaching 11 degrees. Usually high inclinations come with high eccentricities but some simulations show elements compatible with Iapetus' current orbit. A quality factor between 100 and 2000 at the frequency of Titan would bring Titan and Iapetus into a 5:1 resonance, which would perturb Iapetus' eccentricity and inclination to values observed today. Such rapid tidal migration would have avoided Iapetus' ejection around 40 to 800 million years ago.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

We revisit the rotation dynamics of a rigid satellite with either a liquid core or a global sub-surface ocean. In both problems, the flow of the fluid component is assumed inviscid. The study of a hollow satellite with a liquid core is based on the Poincare-Hough model which provides exact equations of motion. We introduce an approximation when the ellipticity of the cavity is low. This simplification allows to model both types of satellite in the same manner. The analysis of their rotation is done in a non-canonical Hamiltonian formalism closely related to Poincare's "forme nouvelle des equations de la mecanique". In the case of a satellite with a global ocean, we obtain a seven-degree of freedom system. Six of them account for the motion of the two rigid components, and the last one is associated with the fluid layer. We apply our model to Titan for which the origin of the obliquity is still a debated question. We show that the observed value is compatible with Titan slightly departing from the hydrostatic equilibrium and being in a Cassini equilibrium state.

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.

We investigate the resonant rotation of circumbinary bodies in planar quasi-circular orbits. Denoting $n_b$ and $n$ the orbital mean motion of the inner binary and of the circumbinary body, respectively, we show that spin-orbit resonances exist at the frequencies $n\pm k\nu/2$, where $\nu = n_b - n$, and $k$ is an integer. Moreover, when the libration at natural frequency has the same magnitude as $\nu$, the resonances overlap and the rotation becomes chaotic. We apply these results to the small satellites in the Pluto-Charon system, and conclude that their rotations are likely chaotic. However, the rotation can also be stable and not synchronous for small axial asymmetries.

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.

The analyses of Titan's gravity field obtained by Cassini space mission suggest the presence of an internal ocean beneath its icy surface. The characterization of the geophysical parameters of the icy shell and the ocean is important to constrain the evolution models of Titan. The knowledge of the librations, that are periodic oscillations around a uniform rotational motion, can bring piece of information on the interior parameters. The objective of this paper is to study the librational response in longitude from an analytical approach for Titan composed of a deep atmosphere, an elastic icy shell, an internal ocean, and an elastic rocky core perturbed by the gravitational interactions with Saturn. We start from the librational equations developed for a rigid satellite in synchronous spin-orbit resonance. We introduce explicitly the atmospheric torque acting on the surface computed from the Titan IPSL GCM (Institut Pierre Simon Laplace General Circulation Model) and the periodic deformations of elastic solid layers due to the tides. We investigate the librational response for various interior models in order to compare and to identify the influence of the geophysical parameters and the impact of the elasticity. The main librations arise at two well-separated forcing frequency ranges: low forcing frequencies dominated by the Saturnian annual and semi-annual frequencies, and a high forcing frequency regime dominated by Titan's orbital frequency around Saturn. We find that internal structure models including an internal ocean with elastic solid layers lead to the same order of libration amplitude than the oceanless models, which makes more challenging to differentiate them by the interpretation of librational motion.

The Cassini spacecraft collects high resolution images of the saturnian satellites and reveals the surface of these new worlds. The shape and rotation of the satellites can be determined from the Cassini Imaging Science Subsystem data, employing limb coordinates and stereogrammetric control points. This is the case for Epimetheus (Tiscareno et al. 2009) that opens elaboration of new rotational models (Tiscareno et al. 2009; Noyelles 2010; Robutel et al. 2011). Especially, Epimetheus is characterized by its horseshoe shape orbit and the presence of the swap is essential to introduce explicitly into rotational models. During its journey in the saturnian system, Cassini spacecraft accumulates the observational data of the other satellites and it will be possible to determine the rotational parameters of several of them. To prepare these future observations, we built rotational models of the coorbital (also called Trojan) satellites Telesto, Calypso, Helene, and Polydeuces, in addition to Janus and Epimetheus. Indeed, Telesto and Calypso orbit around the L_4 and L_5 Lagrange points of Saturn-Tethys while Helene and Polydeuces are coorbital of Dione. The goal of this study is to understand how the departure from the Keplerian motion induced by the perturbations of the coorbital body, influences the rotation of these satellites. To this aim, we introduce explicitly the perturbation in the rotational equations by using the formalism developed by Erdi (1977) to represent the coorbital motions, and so we describe the rotational motion of the coorbitals, Janus and Epimetheus included, in compact form.

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.

The Cassini mission in the Saturnian system is an outstanding opportunity to improve our knowledge of the satellites of Saturn. The data obtained thanks to this mission must be confronted to theoretical models. This paper aims at modeling the rotation of Mimas, with respect to its possible internal structure. For that, we first build different interior models, in considering Mimas as composed of 2 rigid layers with different porosity. Then we simulate the rotational behavior of these models in a 3-degree of freedom numerical code, in considering complete ephemerides of a Mimas whose rotation is disturbed by Saturn. We also estimate the deviation of its longitudinal orientation due to tides. We expect a signature of the internal structure up to 0.53{\deg} in the longitudinal librations and an obliquity between 2 and 3 arcmin, the exact values depending on the interior. The longitudinal librations should be detectable, inverting them to get clues on the internal structure of Mimas is challenging.

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

Context. Ceres is the most massive body of the asteroid belt and contains about 25 wt.% (weight percent) of water. Understanding its thermal evolution and assessing its current state are major goals of the Dawn Mission. Constraints on internal structure can be inferred from various observations. Especially, detailed knowledge of the rotational motion can help constrain the mass distribution inside the body, which in turn can lead to information on its geophysical history. Aims. We investigate the signature of the interior on the rotational motion of Ceres and discuss possible future measurements performed by the spacecraft Dawn that will help to constrain Ceres' internal structure. Methods. We compute the polar motion, precession-nutation, and length-of-day variations. We estimate the amplitudes of the rigid and non-rigid response for these various motions for models of Ceres interior constrained by recent shape data and surface properties. Results. As a general result, the amplitudes of oscillations in the rotation appear to be small, and their determination from spaceborne techniques will be challenging. For example, the amplitudes of the semi-annual and annual nutations are around ~364 and ~140 milli-arcseconds, and they show little variation within the parametric space of interior models envisioned for Ceres. This, combined with the very long-period of the precession motion, requires very precise measurements. We also estimate the timescale for Ceres' orientation to relax to a generalized Cassini State, and we find that the tidal dissipation within that object was probably too small to drive any significant damping of its obliquity since formation. However, combining the shape and gravity observations by Dawn offers the prospect to identify departures of non-hydrostaticity at the global and regional scale, which will be instrumental in constraining Ceres' past and current thermal state. We also discuss the existence of a possible Chandler mode in the rotational motion of Ceres, whose potential excitation by endogenic and/or exogenic processes may help detect the presence of liquid reservoirs within the asteroid.

% context heading (optional) {Thanks to the Venus Express Mission, new data on the properties of Venus could be obtained in particular concerning its rotation.} % aims heading (mandatory) {In view of these upcoming results, the purpose of this paper is to determine and compare the major physical processes influencing the rotation of Venus, and more particularly the angular rotation rate.} % methods heading (mandatory) {Applying models already used for the Earth, the effect of the triaxiality of a rigid Venus on its period of rotation are computed. Then the variations of Venus rotation caused by the elasticity, the atmosphere and the core of the planet are evaluated.} % results heading (mandatory) {Although the largest irregularities of the rotation rate of the Earth at short time scales are caused by its atmosphere and elastic deformations, we show that the Venus ones are dominated by the tidal torque exerted by the Sun on its solid body. Indeed, as Venus has a slow rotation, these effects have a large amplitude of 2 minutes of time (mn). These variations of the rotation rate are larger than the one induced by atmospheric wind variations that can reach 25-50 seconds of time (s), depending on the simulation used. The variations due to the core effects which vary with its size between 3 and 20s are smaller. Compared to these effects, the influence of the elastic deformation cause by the zonal tidal potential is negligible.} % conclusions heading (optional), leave it empty if necessary {As the variations of the rotation of Venus reported here are of the order 3mn peak to peak, they should influence past, present and future observations providing further constraints on the planet internal structure and atmosphere.}

Janus and Epimetheus are famously known for their distinctive horseshoe-shaped orbits resulting from a 1:1 orbital resonance. Every four years these two satellites swap their orbits by a few tens of kilometers as a result of their close encounter. Recently Tiscareno et al. (2009) have proposed a model of rotation based on images from the Cassini orbiter. These authors inferred the amplitude of rotational librational motion in longitude at the orbital period by fitting a shape model to the recent Cassini ISS images. By a quasiperiodic approximation of the orbital motion, we describe how the orbital swap impacts the rotation of the satellites. To that purpose, we have developed a formalism based on quasi-periodic series with long and short-period librations. In this framework, the amplitude of the libration at the orbital period is found proportional to a term accounting for the orbital swap. We checked the analytical quasi-periodic development by performing a numerical simulation and find both results in good agreement. To complete this study, the results regarding the short-period librations are studied with the help of an adiabatic-like approach.

One major objective of MESSENGER and BepiColombo spatial missions is to accurately measure Mercury's rotation and its obliquity in order to obtain constraints on internal structure of the planet. Which is the obliquity's dynamical behavior deriving from a complete spin-orbit motion of Mercury simultaneously integrated with planetary interactions? We have used our SONYR model integrating the spin-orbit N-body problem applied to the solar System (Sun and planets). For lack of current accurate observations or ephemerides of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical integration, we have built an original method for finding the libration center of the spin-orbit system and, as a consequence, for avoiding arbitrary amplitudes in librations of the spin-orbit motion as well as in Mercury's obliquity. The method has been carried out in two cases: (1) the spin-orbit motion of Mercury in the 2-body problem case (Sun-Mercury) where an uniform precession of the Keplerian orbital plane is kinematically added at a fixed inclination (S2K case), (2) the spin-orbit motion of Mercury in the N-body problem case (Sun and planets) (Sn case). We find that the remaining amplitude of the oscillations in the Sn case is one order of magnitude larger than in the S2K case, namely 4 versus 0.4 arcseconds (peak-to-peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for a difference of 10.8 arcseconds. These theoretical results are in a good agreement with recent radar observations but it is not excluded that it should be possible to push farther the convergence process by drawing nearer still more precisely to the libration center.

The upcoming space missions, MESSENGER and BepiColombo with onboard instrumentation capable of measuring the rotational parameters stimulate the objective to reach an accurate theory of the rotational motion of Mercury. For obtaining the real motion of Mercury, we have used our relativistic BJV model of solar system integration including the coupled spin-orbit motion of the Moon. We have extended this model by generalizing the spin-orbit couplings to the terrestrial planets, notably Mercury. The updated model is called SONYR (acronym of Spin-Orbit N-BodY Relativistic model). The SONYR model giving an accurate solution of the spin-orbit motion of Mercury permits to analyze the different families of the Hermean rotational librations. The spin-orbit motion of Mercury is characterized by two proper frequencies (namely 15.847 and 1066 years) and its 3:2 resonance presents a second synchronism which can be understood as a spin-orbit secular resonance (278 898 years). By using the SONYR model, we find a new determination of the mean obliquity, namely 1.6 arcminutes. Besides, we identify in the Hermean librations the impact of the uncertainty of the greatest principal moment of inertia (\cmr2) on the obliquity and on the libration in longitude (2.3 mas and 0.45 as respectively for an increase of 1% on the \cmr2 value). These determinations prove to be suitable for providing constraints on the internal structure of Mercury.

The orbits in the HD 160691 planetary system at first appeared highly unstable, but using the MEGNO and FLI techniques of global dynamics analysis in the orbital parameter space we have found a stabilizing mechanism that could be the key to its existence. In order to be dynamically stable, the HD 160691 planetary system has to satisfy the following conditions: (1) a 2:1 mean motion resonance, combined with (2) an apsidal secular resonance in (3) a configuration $P_{c}(ap) - S - P_{b}(ap)$ where the two apsidal lines are anti-aligned, and (4) specific conditions on the respective sizes of the eccentricities (high eccentricity for the outer orbit is in particular the most probable necessary condition). More generally, in this original orbital topology, where the resonance variables $\theta_{1}$ and $\theta_{3}$ librate about $180^{\circ}$ while $\theta_{2}$ librates about $0^{\circ}$, the HD 160691 system and its mechanism have revealed aspects of the 2:1 orbital resonances that have not been observed nor analyzed before. The present topology combined with the 2:1 resonance is indeed more wide-ranging than the particular case of the HD 160691 planetary system. It is a new theoretical possibility suitable for a stable regime despite relatively small semi-major axes with respect to the important masses in interactions.

The 3:2 spin-orbit resonance between the rotational and orbital motions of Mercury results from a functional dependance on a tidal friction adding to a non-zero eccentricity with a permanent asymmetry in the equatorial plane of the planet. The upcoming space missions, MESSENGER and BepiColombo with onboard instrumentation capable of measuring the Mercury's rotational parameters, stimulate the objective to attempt to an accurate theory of the planet's rotation. We have used our BJV relativistic model of solar system integration including the spin-orbit motion of the Moon. This model had been previously built in accordance with the requirements of the Lunar Laser Ranging observational accuracy. We extended this model to the spin-orbit couplings of the terrestrial planets including Mercury; the updated model is called SONYR (acronym of Spin-Orbit N-BodY Relativistic model). An accurate rotation of Mercury has been then obtained. Moreover, the conception of the SONYR model is suitable for analyzing the different families of hermean rotational librations. We accurately identify the non-linear relations between the rotation of Mercury and its dynamical figure (\cmr2, $C_{20}$, and $C_{22}$). Notably, for a variation of 1% on the \cmr2 value, signatures in the $\phi$ hermean libration in longitude as well as in the $\eta$ obliquity of the planet are respectively 0.45 arcseconds (as) and 2.4 milliarcseconds (mas). These determinations provide new constraints on the internal structure of Mercury to be discussed with the expected accuracy forecasted in the BepiColombo mission (respectively 3.2 and 3.7 as according to Milani et al 2001).

In order to distinguish between regular and chaotic planetary orbits we apply a new technique called MEGNO in a wide neighbourhood of orbital parameters determined using standard two-body Keplerian fits for HD 12661, HD 38529, HD 37124 and HD 160691 planetary systems. We show that the currently announced orbital parameters place these systems in very different situations from the point of view of dynamical stability. While HD 38529 and HD 37124 are located within large stability zones in the phase space around their determined orbits, the preliminary orbits in HD 160691 are highly unstable. The orbital parameters of the HD 12661 planets are located in a border region between stable and unstable dynamical regimes, so while its currently determined orbital parameters produce stable regular orbits, a minor change within the margin of error of just one parameter may result in a chaotic dynamical system.