In Gayon & Bois (2008) and Gayon etal (2009), (i) we studied the theoretical feasibility and efficiency of retrograde mean motion resonances (i.e. two planets are both in orbital resonance and in counter-revolving configuration), (ii) we showed that retrograde resonances can generate interesting mechanisms of stability, and (iii) we obtained a dynamical fit involving a counter-revolving configuration that is consistent with the observations of the HD73526 planetary system. In the present paper, we present and analyze data reductions assuming counter-revolving configurations for eight compact multi-planetary systems detected through the radial velocity method. In each case, we select the best fit leading to a dynamically stable solution. The resulting data reductions obtained in rms and chi values for counter-revolving configurations are of the same order, and sometimes slightly better, than for prograde configurations. In the end, these fits tend to show that, over the eight studied multi-planetary systems, six of them could be regulated by a mechanism involving a counter-revolving configuration.

In a previous paper (Gayon & Bois 2008a), we have shown the general efficiency of retrograde resonances for stabilizing compact planetary systems. Such retrograde resonances can be found when two-planets of a three-body planetary system are both in mean motion resonance and revolve in opposite directions. For a particular two-planet system, we have also obtained a new orbital fit involving such a counter-revolving configuration and consistent with the observational data. In the present paper, we analytically investigate the three-body problem in this particular case of retrograde resonances. We therefore define a new set of canonical variables allowing to express correctly the resonance angles and obtain the Hamiltonian of a system harboring planets revolving in opposite directions. The acquiring of an analytical "rail" may notably contribute to a deeper understanding of our numerical investigation and provides the major structures related to the stability properties. A comparison between our analytical and numerical results is also carried out.

Most multi-planetary systems are characterized by hot-Jupiters close to their central star, moving on eccentric orbits. From a dynamical point of view, compact multi-planetary systems form a specific class of the general N-body problem (where N>=3). Moreover, extrasolar planets are found in prograde orbits about their host star, and often in mean motion resonances (MMR). In a first step, we study theoretically a new stabilizing mechanism suitable for compact two-planet systems. This mechanism involves counter-revolving orbits forming a retrograde MMR. In a second step, we investigate the feasibility of planetary systems hosting counter-revolving planets. Dynamical stability, observations, and formation processes of these systems are analyzed and discussed. To characterize the dynamical behavior of multi-dimensional planetary systems, we apply our technique of global dynamics analysis based on the MEGNO indicator (Mean Exponential Growth factor of Nearby Orbits) that provides the fine structure of the phase space. In a few cases of possible counter-revolving configurations, we carry out new fits to the observations using the Pikaia genetic algorithm. A statistical study of the stability in the neighborhood of different observed, planetary-systems is completed using a Monte-Carlo method. We analyse the observational data for the HD73526 planetary system and find that counter-revolving configurations may be consistent with the observational data. We highlight the fine and characteristic structure of retrograde MMRs. We demonstrate that retrograde resonances open a family of stabilizing mechanisms involving new apsidal precession behaviors. Considering two possible formation mechanisms (free-floating planet and Slingshot model), we conclude that counter-revolving configurations are feasible.

Multi-planet systems detected until now are in most cases characterized by hot-Jupiters close to their central star as well as high eccentricities. As a consequence, from a dynamical point of view, compact multi-planetary systems form a variety of the general N-body problem (with N >= 3), whose solutions are not necessarily known. Extrasolar planets are up to now found in prograde (i.e. direct) orbital motions about their host star and often in mean-motion resonances (MMR). In the present paper, we investigate a theoretical alternative suitable for the stability of compact multi-planetary systems. When the outer planet moves on a retrograde orbit in MMR with respect to the inner planet, we find that the so-called retrograde resonances present fine and characteristic structures particularly relevant for dynamical stability. We show that retrograde resonances and their resources open a family of stabilizing mechanisms involving specific behaviors of apsidal precessions. We also point up that for particular orbital data, retrograde MMRs may provide more robust stability compared to the corresponding prograde MMRs.

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

We discuss the usual account of causal structure that relies on the temporal precedence constraint between cause-effect pairs. In particular, we consider the subtle interplay between local and global characters of time and causality encoded in General Relativity. We find that the temporal precedence constraint must be relaxed. We argue that these results could lead to an enriched structure of the physical causality.

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.

We consider dynamical effects of additional perturbative forces due to the non-point mass nature of stars and planets: effects such as quadrupolar distortion and tidal friction in the systems of exo-planets. It is shown that these forces should not be neglected while modelling the dynamics of planetary systems, especially taking into account the undefined real masses of the planets due to unknown orbital inclinations and the unsatisfactory application of Keplerian fits to the radial velocity data in multiple planetary systems.