The Spin-Orbit Resonant Rotation of Mercury: A Two Degree of Freedom Hamiltonian Model

The paper develops a hamiltonian formulation describing the coupled orbital and spin motions of a rigid Mercury rotation about its axis of maximum moment of inertia in the frame of a 3:2 spin orbit resonance; the (ecliptic) obliquity is not constant, the gravitational potential of mercury is developed up to the second degree terms (the only ones for which an approximate numerical value can be given) and is reduced to a two degree of freedom model in the absence of planetary perturbations. Four equilibria can be calculated, corresponding to four different values of the (ecliptic) obliquity. The present situation of Mercury corresponds to one of them, which is proved to be stable. We introduce action-angle variables in the neighborhood of this stable equilibrium, by several successive canonical transformations, so to get two constant frequencies, the first one for the free spin-orbit libration, the other one for the 1:1 resonant precession of both nodes (orbital and rotational) on the ecliptic plane. The numerical values obtained by this simplified model are in perfect agreement with those obtained by Rambaux and Bois [Astron. Astrophys. 413, 381–393]. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the oscillations in Mercury's obliquity

Mercury's capture into the 3:2 spin–orbit resonance can be explained as a result of its chaotic orbital dynamics. 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. Analytical approaches at the first-order level using the Cassini state assumptions give the obliquity constant or quasi-constant. 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 (acronym of Spin–Orbit N-bodY Relativistic 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 (S_2K case), (2) the spin–orbit motion of Mercury in the N-body problem case (Sun and planets) (S_n case). We find that the remaining amplitude of the oscillations in the S_n case is one order of magnitude larger than in the S_2K 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. We note that the dynamically driven spin precession, which occurs when the planetary interactions are included, is more complex than the purely kinematic case. Nevertheless, in such a N-body problem, we find that the 3:2 spin–orbit resonance is really combined to a synchronism where the spin and orbit poles on average precess at the same rate while the orbit inclination and the spin axis orientation on average decrease at the same rate. As a consequence and whether it would turn out that there exists an irreducible minimum of the oscillation amplitude, quasi-periodic oscillations found in Mercury's obliquity should be to geometrically understood as librations related to these synchronisms that both follow a Cassini state. Whatever the open question on the minimal amplitude in the obliquity's oscillations and in spite of the planetary interactions indirectly acting by the solar torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium state that proceeds from a 2-body Cassini state.     

Mercury's capture into the 3/2 spin-orbit resonance including the effect of core-mantle friction

The rotation of Mercury is presently captured in a 3/2 spin-orbit resonance with the orbital mean motion. The capture mechanism is well understood as the result of tidal interactions with the Sun combined with planetary perturbations [Goldreich, P., Peale, S., 1966. Astron. J. 71, 425-438; Correia, A.C.M., Laskar, J., 2004. Nature 429, 848-850]. However, it is now almost certain that Mercury has a liquid core [Margot, J.L., Peale, S.J., Jurgens, R.F., Slade, M.A., Holin, I.V., 2007. Science 316, 710-714] which should induce a contribution of viscous friction at the core-mantle boundary to the spin evolution. According to Peale and Boss [Peale, S.J., Boss, A.P., 1977. J. Geophys. Res. 82, 743-749] this last effect greatly increases the chances of capture in all spin-orbit resonances, being 100% for the 2/1 resonance, and thus preventing the planet from evolving to the presently observed configuration. Here we show that for a given resonance, as the chaotic evolution of Mercury's orbit can drive its eccentricity to very low values during the planet's history, any previous capture can be destabilized whenever the eccentricity becomes lower than a critical value. In our numerical integrations of 1000 orbits of Mercury over 4 Gyr, the spin ends 99.8% of the time captured in a spin-orbit resonance, in particular in one of the following three configurations: 5/2 (22%), 2/1 (32%) and 3/2 (26%). Although the present 3/2 spin-orbit resonance is not the most probable outcome, we also show that the capture probability in this resonance can be increased up to 55% or 73%, if the eccentricity of Mercury in the past has descended below the critical values 0.025 or 0.005, respectively.     

Librations and obliquity of Mercury from the BepiColombo radio-science and camera experiments

A major goal of the BepiColombo mission to Mercury is the determination of the structure and state of Mercury's interior. Here the BepiColombo rotation experiment has been simulated in order to assess the ability to attain the mission goals and to help lay out a series of constraints on the experiment's possible progress. In the rotation experiment pairs of images of identical surface regions taken at different epochs are used to retrieve information on Mercury's rotation and orientation. The idea is that from observations of the same patch of Mercury's surface at two different solar longitudes of Mercury the orientation of Mercury can be determined, and therefore also the obliquity and rotation variations with respect to the uniform rotation.The estimation of the libration amplitude and obliquity through pattern matching of observed surface landmarks is challenging. The main problem arises from the difficulty to observe the same landmark on the planetary surface repeatedly over the MPO mission lifetime, due to the combination of Mercury's 3:2 spin-orbit resonance, the absence of a drift of the MPO polar orbital plane and the need to combine data from different instruments with their own measurement restrictions.By assuming that Mercury occupies a Cassini state and that the spacecraft operates nominally we show that under worst case assumptions the annual libration amplitude and obliquity can be measured with a precision of, respectively, 1.4 arcseconds (as) and 1.0 as over the nominal BepiColombo MPO lifetime with about 25 landmarks for rather stringent illumination restrictions. The outcome of the experiment cannot be easily improved by simply relaxing the observational constraints, or increasing the data volume.     

The proximity of Mercury's spin to Cassini state 1 from adiabatic invariance

In determining Mercury's core structure from its rotational properties, the value of the normalized moment of inertia, C/MR², from the location of Cassini 1 is crucial. If Mercury's spin axis occupies Cassini state 1, its position defines the location of the state, where the axis is fixed in the frame precessing with the orbit. Although tidal and core-mantle dissipation drive the spin to the Cassini state with a time scale O(¹⁰⁵) years, the spin might still be displaced from the Cassini state if the variations in the orbital elements induced by planetary perturbations, which change the position of the Cassini state, cause the spin to lag behind as it attempts to follow the state. After being brought to the state by dissipative processes, the spin axis is expected to follow the Cassini state for orbit variations with time scales long compared to the 1000 year precession period of the spin about the Cassini state because the solid angle swept out by the spin axis as it precesses is an adiabatic invariant. Short period variations in the orbital elements of small amplitude should cause displacements that are commensurate with the amplitudes of the short period terms. The exception would be if there are forcing terms in the perturbations that are nearly resonant with the 1000 year precession period. The precision of the radar and eventual spacecraft measurements of the position of Mercury's spin axis warrants a check on the likely proximity of the spin axis to the Cassini state. How confident should we be that the spin axis position defines the Cassini state sufficiently well for a precise determination of C/MR²? By following simultaneously the spin position and the Cassini state position during long time scale orbital variations over past 3 million years [Quinn, T.R., Tremaine, S., Duncan, M., 1991. Astron. J. 101, 2287-2305] and short time scale variations for 20,000 years [JPL Ephemeris DE 408; Standish, E.M., private communication, 2005], we show that the spin axis will remain within one arcsec of the Cassini state after it is brought there by dissipative torques. In this process the spin is located in the orbit frame of reference, which in turn is referenced to the inertial ecliptic plane of J2000. There are no perturbations with periods resonant with the precession period that could cause large separations. We thus expect Mercury's spin to occupy Cassini state 1 well within the uncertainties for both radar and spacecraft measurements, with correspondingly tight constraints on C/MR² and the extent of Mercury's molten core. Two unlikely caveats for this conclusion are: (1) an excitation of a free spin precession by an unknown mechanism or (2) a displacement by a dissipative core mantle interaction that exceeds the measurement uncertainties.     

Possibility of self-sustaining bombardment of inner planets

The great strengthening the material undergoes under high confining pressure, and jet pattern of matter outflowing from large impact craters make possible the ejection of asteroid-size bodies from the Earth into space. The ejected bodies, after gaining energy in planetary perturbations, may fall back with a velocity higher than that of their ejection. This solves, in particular, the problem of shower bombardments with ~ 25 Myr interval (Drobyshevski, Sov. Astron. Let. 16(3), 193, 1990), and a question arises whether this process could become self-sustained, like a chain reaction, when secondary impacts release an energy higher than that of primary impact. Estimates show that such a possibility could have been realized for Mercury (Drobyshevski, Lunar Planet. Sci. Conf. Abstr. 23(1), 317, 1992) due to its low escape and high orbital velocities. Self-sustained bombardment can account for the loss of the silicate mantle from Mercury. The energy and angular momentum conservation laws imply that its orbit contracted toward the Sun in the course of ejection of the mantle fragments by Mercury's perturbations beyond its orbit. Straight-forward calculations show the initial orbit to have practically coincided with the Venusian orbit. This puts the old hypothesis of Mercury being a lost satellite of Venus on a solid ground and provides an explanation for many facts from the origin of the Imbrium bombardment to the observed locks in the axial and orbital rotation of Mercury, Venus, and the Earth.     

Global MHD modeling of Mercury's magnetosphere with applications to the MESSENGER mission and dynamo theory

We use a global magnetohydrodynamic (MHD) model to simulate Mercury's space environment for several solar wind and interplanetary magnetic field (IMF) conditions in anticipation of the magnetic field measurements by the MESSENGER spacecraft. The main goal of our study is to assess what characteristics of the internally generated field of Mercury can be inferred from the MESSENGER observations, and to what extent they will be able to constrain various models of Mercury's magnetic field generation. Based on the results of our simulations, we argue that it should be possible to infer not only the dipole component, but also the quadrupole and possibly even higher harmonics of the Mercury's planetary magnetic field. We furthermore expect that some of the crucial measurements for specifying the Hermean internal field will be acquired during the initial fly-bys of the planet, before MESSENGER goes into orbit around Mercury.     

Resonant forcing of Mercury's libration in longitude

The period of free libration of Mercury's longitude about the position it would have had if it were rotating uniformly at 1.5 times its orbital mean motion is close to resonance with Jupiter's orbital period. The Jupiter perturbations of Mercury's orbit thereby lead to amplitudes of libration at the 11.86 year period that may exceed the amplitude of the 88 day forced libration determined by radar. Mercury's libration in longitude may be thus dominated by only two periods of 88 days and 11.86 years, where other periods from the planetary perturbations of the orbit have much smaller amplitudes.     

Obliquity variations of a moonless Earth

We numerically explore the obliquity (axial tilt) variations of a hypothetical moonless Earth. Previous work has shown that the Earth’s Moon stabilizes Earth’s obliquity such that it remains within a narrow range, between 22.1° and 24.5°. Without lunar influence, a frequency map analysis by Laskar et al. (Laskar, J., Joutel, F., Robutel, P. [1993]. Nature 361, 615–617) showed that the obliquity could vary between 0° and 85°. This has left an impression in the astrobiology community that a big moon is necessary to maintain a habitable climate on an Earth-like planet. Using a modified version of the orbital integrator mercury, we calculate the obliquity evolution for moonless Earths with various initial conditions for up to 4 Gyr. We find that while obliquity varies significantly more than that of the actual Earth over 100,000 year timescales, the obliquity remains within a constrained range, typically 20–25° in extent, for timescales of hundreds of millions of years. None of our Solar System integrations in which planetary orbits behave in a typical manner show obliquity accessing more than 65% of the full range allowed by frequency-map analysis. The obliquities of moonless Earths that rotate in the retrograde direction are more stable than those of prograde rotators. The total obliquity range explored for moonless Earths with rotation periods less than 12 h is much less than that for slower-rotating moonless Earths. A large moon thus does not seem to be needed to stabilize the obliquity of an Earth-like planet on timescales relevant to the development of advanced life.     

Moonlet induced wakes in planetary rings: Analytical model including eccentric orbits of moon and ring particles

Saturn’s rings host two known moons, Pan and Daphnis, which are massive enough to clear circumferential gaps in the ring around their orbits. Both moons create wake patterns at the gap edges by gravitational deflection of the ring material (Cuzzi, J.N., Scargle, J.D. [1985]. Astrophys. J. 292, 276-290; Showalter, M.R., Cuzzi, J.N., Marouf, E.A., Esposito, L.W. [1986]. Icarus 66, 297-323). New Cassini observations revealed that these wavy edges deviate from the sinusoidal waveform, which one would expect from a theory that assumes a circular orbit of the perturbing moon and neglects particle interactions. Resonant perturbations of the edges by moons outside the ring system, as well as an eccentric orbit of the embedded moon, may partly explain this behavior (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S., Burns, J.A., Dones, L. [2005]. Bull. Am. Astron. Soc. 37, 767; Weiss, J.W., Porco, C.C., Tiscareno, M.S. [2009]. Astron. J. 138, 272-286). Here we present an extended non-collisional streamline model which accounts for both effects. We describe the resulting variations of the density structure and the modification of the nonlinearity parameter q. Furthermore, an estimate is given for the applicability of the model. We use the streamwire model introduced by Stewart (Stewart, G.R. [1991]. Icarus 94, 436-450) to plot the perturbed ring density at the gap edges.We apply our model to the Keeler gap edges undulated by Daphnis and to a faint ringlet in the Encke gap close to the orbit of Pan. The modulations of the latter ringlet, induced by the perturbations of Pan (Burns, J.A., Hedman, M.M., Tiscareno, M.S., Nicholson, P.D., Streetman, B.J., Colwell, J.E., Showalter, M.R., Murray, C.D., Cuzzi, J.N., Porco, C.C., and the Cassini ISS team [2005]. Bull. Am. Astron. Soc. 37, 766), can be well described by our analytical model. Our analysis yields a Hill radius of Pan of 17.5 km, which is 9% smaller than the value presented by Porco (Porco, C.C., and 34 colleagues [2005]. Science 307, 1226-1236), but fits well to the radial semi-axis of Pan of 17.4 km. This supports the idea that Pan has filled its Hill sphere with accreted material (Porco, C.C., Thomas, P.C., Weiss, J.W., Richardson, D.C. [2007]. Science 318, 1602-1607). A numerical solution of a streamline is used to estimate the parameters of the Daphnis-Keeler gap system, since the close proximity of the gap edge to the moon induces strong perturbations, not allowing an application of the analytic streamline model. We obtain a Hill radius of 5.1 km for Daphnis, an inner edge variation of 8 km, and an eccentricity for Daphnis of 1.5 × 10⁻⁵. The latter two quantities deviate by a factor of two from values gained by direct observations (Jacobson, R.A., Spitale, J., Porco, C.C., Beurle, K., Cooper, N.J., Evans, M.W., Murray, C.D. [2008]. Astron. J. 135, 261-263; Tiscareno, M.S., Burns, J.A., Hedman, M.M., Spitale, J.N., Porco, C.C., Murray, C.D., and the Cassini Imaging team [2005]. Bull. Am. Astron. Soc. 37, 767), which might be attributed to the neglect of particle interactions and vertical motion in our model.     

Nonsingular vectorial reformulation of the short-period corrections in Kozai’s oblateness solution

We present fully three-dimensional equations to describe the rotations of a body made of a deformable mantle and a fluid core. The model in its essence is similar to that used by INPOP19a (Integration Planétaire de l’Observatoire de Paris) Fienga et al. (INPOP19a planetary ephemerides. Notes Scientifiques et Techniques de l’Institut de Mécanique Céleste, vol 109, 2019), and by JPL (Jet Propulsion Laboratory) (Park et al. The JPL Planetary and Lunar Ephemerides DE440 and DE441. Astron J 161(3):105, 2021. doi:10.3847/1538-3881/abd414), to represent the Moon. The intended advantages of our model are: straightforward use of any linear-viscoelastic model for the rheology of the mantle; easy numerical implementation in time-domain (no time lags are necessary); all parameters, including those related to the “permanent deformation”, have a physical interpretation. The paper also contains: (1) A physical model to explain the usual lack of hydrostaticity of the mantle (permanent deformation). (2) Formulas for free librations of bodies in and out-of spin-orbit resonance that are valid for any linear viscoelastic rheology of the mantle. (3) Formulas for the offset between the mantle and the idealised rigid-body motion (Peale’s Cassini states). (4) Applications to the librations of Moon, Earth, and Mercury that are used for model validation.

     

Analytical model of the long-period forced longitude librations of Mercury

The shaking of Mercury’s orbit by the planets forces librations in longitude in addition to those at harmonics of the orbital period that have been used to detect Mercury’s molten core. We extend the analytical formulation of Peale et al. (Peale, S.J., Margot, J.L., Yseboodt, M. [2009]. Icarus 199, 1-8) in order to provide a convenient means of determining the amplitudes and phases of the forced librations without resorting to numerical calculations. We derive an explicit relation between the amplitude of each forced libration and the moment of inertia parameter (B-A)/C_m. Far from resonance with the free libration period, the libration amplitudes are directly proportional to (B-A)/C_m. Librations with periods close to the free libration period of ∼12 years may have measurable (∼arcsec) amplitudes. If the free libration period is sufficiently close to Jupiter’s orbital period of 11.86 years, the amplitude of the forced libration at Jupiter’s period could exceed the 35 arcsec amplitude of the 88-day forced libration. We also show that the planetary perturbations of the mean anomaly and the longitude of pericenter of Mercury’s orbit completely determine the libration amplitudes.While these signatures do not affect spin rate at a detectable level (as currently measured by Earth-based radar), they have a much larger impact on rotational phase (affecting imaging, altimetry, and gravity sensors). Therefore, it may be important to consider planetary perturbations when interpreting future spacecraft observations of the librations.     

Numerical integration of dynamical systems with Lie series

The integration of the equations of motion in gravitational dynamical systems—either in our Solar System or for extra-solar planetary systems—being non integrable in the global case, is usually performed by means of numerical integration. Among the different numerical techniques available for solving ordinary differential equations, the numerical integration using Lie series has shown some advantages. In its original form (Hanslmeier and Dvorak, Astron Astrophys 132, 203 1984), it was limited to the N-body problem where only gravitational interactions are taken into account. We present in this paper a generalisation of the method by deriving an expression of the Lie terms when other major forces are considered. As a matter of fact, previous studies have been done but only for objects moving under gravitational attraction. If other perturbations are added, the Lie integrator has to be re-built. In the present work we consider two cases involving position and position-velocity dependent perturbations: relativistic acceleration in the framework of General Relativity and a simplified force for the Yarkovsky effect. A general iteration procedure is applied to derive the Lie series to any order and precision. We then give an application to the integration of the equation of motions for typical Near-Earth objects and planet Mercury.     

The free precession and libration of Mercury

An analysis based on the direct torque equations including tidal dissipation and a viscous core-mantle coupling is used to determine the damping time scales of O(¹⁰⁵) years for free precession of the spin about the Cassini state and free libration in longitude for Mercury. The core-mantle coupling dominates the damping over the tides by one to two orders of magnitude for the plausible parameters chosen. The short damping times compared with the age of the Solar System means we must find recent or on-going excitation mechanisms if such free motions are found by the current radar experiments or the future measurement by the MESSENGER and BepiColombo spacecraft that will orbit Mercury. We also show that the average precession rate is increased by about 30% over that obtained from the traditional precession constant because of a spin-orbit resonance induced contribution by the C₂₂ term in the expansion of the gravitational field. The C₂₂ contribution also causes the path of the spin during the precession to be slightly elliptical with a variation in the precession rate that is a maximum when the obliquity is a minimum. An observable free precession will compromise the determination of obliquity of the Cassini state and hence of C/MMR² for Mercury, but a detected free libration will not compromise the determination of the forced libration amplitude and thus the verification of a liquid core.     

The effects of ablation on the cross section of planetary envelopes at capturing planetesimals

We explore the cross section of giant planet envelopes at capturing planetesimals of different sizes. For this purpose we employ two sets of realistic planetary envelope models (computed assuming for the protoplanetary nebula masses of 10 and 5 times the mass of the minimum mass solar nebula), account for drag and ablation effects and study the trajectories along which planetesimals move. The core accretion of these models has been computed in the oligarchic growth regime [Fortier, A., Benvenuto, O.G., Brunini, A., 2007. Astron. Astrophys. 473, 311-322], which has also been considered for the velocities of the incoming planetesimals. This regime predicts velocities larger that those used in previous studies of this problem. As the rate of ablation is dependent on the third power of velocity, ablation is more important in the oligarchic growth regime. We compute energy and mass deposition, fractional ablated masses and the total cross section of planets for a wide range of values of the critical parameter of ablation. In computing the total cross section of the planet we have included the contributions due to mass deposited by planetesimals moving along unbound orbits. Our results indicate that, for the case of small planetary cores and low velocities for the incoming planetesimals, ablation has a negligible impact on the capture cross section in agreement with the results presented in Inaba and Ikoma [Inaba, S., Ikoma, M., 2003. Astron. Astrophys. 410, 711-723]. However for the case of larger cores and high velocities of the incoming planetesimals as predicted by the oligarchic growth regime, we find that ablation is important in determining the planetary cross section, being several times larger than the value corresponding ignoring ablation. This is so regardless of the size of the incoming planetesimals.     

Possibility of self-sustaining bombardment of inner planets

The great strengthening the material undergoes under high confining pressure, and jet pattern of matter outflowing from large impact craters make possible the ejection of asteroid-size bodies from the Earth into space. The ejected bodies, after gaining energy in planetary perturbations, may fall back with a velocity higher than that of their ejection. This solves, in particular, the problem of shower bombardments with ~ 25 Myr interval (Drobyshevski, Sov. Astron. Let. 16(3), 193, 1990), and a question arises whether this process could become self-sustained, like a chain reaction, when secondary impacts release an energy higher than that of primary impact. Estimates show that such a possibility could have been realized for Mercury (Drobyshevski, Lunar Planet. Sci. Conf. Abstr. 23(1), 317, 1992) due to its low escape and high orbital velocities. Self-sustained bombardment can account for the loss of the silicate mantle from Mercury. The energy and angular momentum conservation laws imply that its orbit contracted toward the Sun in the course of ejection of the mantle fragments by Mercury's perturbations beyond its orbit. Straight-forward calculations show the initial orbit to have practically coincided with the Venusian orbit. This puts the old hypothesis of Mercury being a lost satellite of Venus on a solid ground and provides an explanation for many facts from the origin of the Imbrium bombardment to the observed locks in the axial and orbital rotation of Mercury, Venus, and the Earth.     

The sodium exosphere of Mercury: Comparison between observations during Mercury's transit and model results

In this study we compare the sodium exosphere observations made by Schleicher et al. [Schleicher, H., and 4 colleagues, 2004. Astron. Astrophys. 425, 1119-1124] with the result of a detailed numerical simulation. The observations, made during the transit of Mercury across the solar disk on 7 May 2003, show a maximum of sodium emission near the polar regions, with north prevalence, and the presence of a dawn-dusk asymmetry. We interpret this distribution as the resulting effect of two combined processes: the solar wind proton precipitation causing chemical alteration of the surface, freeing the sodium atoms from their bounds in the crystalline structure on the surface, and the subsequent photon-stimulated and thermal desorption of the sodium atoms. While we find that the velocity distribution of photon desorbed sodium can explain the observed exosphere population, thermal desorption seems to play a minor role only causing a smearing at the locations where Na atoms are released on the dayside. The observed and simulated distributions agree very well with this hypothesis and indicate that the combination of the proposed processes is able to explain the observed features.     

Planetary long periodic terms in Mercury’s rotation: a two dimensional adiabatic approach

The averaged spin-orbit resonant motion of Mercury is considered, with e the orbital eccentricity, and i _o the orbital inclination introduced as very slow functions of time, given by any secular planetary theory. The basis is our Hamiltonian approach (D’Hoedt, S., Lemaître, A.: Celest. Mech. Dyn. Astron. 89:267–283, 2004) in which Mercury is considered as a rigid body. The model is based on two degrees of freedom; the first one is linked to the 3:2 resonant spin-orbit motion, and the second one to the commensurability of the rotational and orbital nodes. Mercury is assumed to be very close to the Cassini equilibrium of the model. To follow the motion of rotation close to this equilibrium, which varies with respect to time through e and i _o , we use the adiabatic invariant theory, extended to two degrees of freedom. We calculate the corrections (remaining functions) introduced by the time dependence of e and i _o in the three steps necessary to characterize the frequencies at the equilibrium. The conclusion is that Mercury follows the Cassini equilibrium (stays in the Cassini forced state), in an adiabatic behavior: the area around the equilibrium does not change by more than ${\varepsilon}$ for times smaller than ${\frac{1}{\varepsilon}}$ . The role of the inclination and the eccentricity can be dissociated and measured in each step of the canonical transformation.     

Models of the collisional damping scenario for ice-giant planets and Kuiper belt formation

Chiang et al. [Chiang, E., Lithwick, Y., Murray-Clay, R., Buie, M., Grundy, W., Holman, M., 2007. In: Protostars and Planets V, pp. 895-911] have recently proposed that the observed structure of the Kuiper belt could be the result of a dynamical instability of a system of ∼5 primordial ice-giant planets in the outer Solar System. According to this scenario, before the instability occurred, these giants were growing in a highly collisionally damped environment according to the arguments in Goldreich et al. [Goldreich, P., Lithwick, Y., Sari, R., 2004. Astrophys. J. 614, 497-507; Annu. Rev. Astron. Astrophys. 42, 549-601]. Here we test this hypothesis with a series of numerical simulations using a new code designed to incorporate the dynamical effects of collisions. We find that we cannot reproduce the observed Solar System. In particular, Goldreich et al. [Goldreich, P., Lithwick, Y., Sari, R., 2004. Astrophys. J. 614, 497-507; Annu. Rev. Astron. Astrophys. 42, 549-601] and Chiang et al. [Chiang, E., Lithwick, Y., Murray-Clay, R., Buie, M., Grundy, W., Holman, M., 2007. In: Protostars and Planets V, pp. 895-911] argue that during the instability, all but two of the ice giants would be ejected from the Solar System by Jupiter and Saturn, leaving Uranus and Neptune behind. We find that ejections are actually rare and that instead the systems spread outward. This always leads to a configuration with too many planets that are too far from the Sun. Thus, we conclude that both Goldreich et al.'s scheme for the formation of Uranus and Neptune and Chiang et al.'s Kuiper belt formation scenario are not viable in their current forms.     

Planetary perturbations on Mercury’s libration in longitude

Two space missions dedicated to Mercury (MESSENGER and BepiColombo) aim at understanding its rotation and confirming the existence of a liquid core. This double challenge requires much more accurate models for the spin-orbit resonant rotation of Mercury. The purpose of this paper is to introduce planetary perturbations on Mercury’s rotation using an analytical method and to analyse the influence of the perturbations on the libration in longitude. Applying a perturbation theory based on the Lie triangle, we were able to re-introduce short periodic terms into the averaged Hamiltonian and to compute the evolution of the rotational variables. The perturbations on Mercury’s forced libration in longitude mainly come from the orbital motion of Mercury (with an amplitude around 41 arcsec that depends on the momenta of inertia). It is completed by various effects from Jupiter (11.86 and 5.93 year-periods), Venus (with a 5.66 year-period), Saturn (14.73 year-period), and the Earth (6.58 year-period). The amplitudes of the oscillations due to Jupiter and Venus are approximately 33% and 10% of those from the orbital motion of Mercury and the amplitudes of the oscillations due to Saturn and the Earth are approximately 3% and 2%. We compare the analytical results with the solution obtained from the spin-orbit numerical model SONYR.