Long-term evolution of the spin of Mercury: I. Effect of the obliquity and core-mantle friction

The present obliquity of Mercury is very low (less than 0.1°), which led previous studies to always adopt a nearly zero obliquity during the planet’s past evolution. However, the initial orientation of Mercury’s rotation axis is unknown and probably much different than today. As a consequence, we believe that the obliquity could have been significant when the rotation rate of the planet first encountered spin-orbit resonances. In order to compute the capture probabilities in resonance for any evolutionary scenario, we present in full detail the dynamical equations governing the long-term evolution of the spin, including the obliquity contribution.The secular spin evolution of Mercury results from tidal interactions with the Sun, but also from viscous friction at the core-mantle boundary. Here, this effect is also regarded with particular attention. Previous studies show that a liquid core enhances drastically the chances of capture in spin-orbit resonances. We confirm these results for null obliquity, but we find that the capture probability generally decreases as the obliquity increases. We finally show that, when core-mantle friction is combined with obliquity evolution, the spin can evolve into some unexpected configurations as the synchronous or the 1/2 spin-orbit resonance.     

Observable consequences of planet formation models in systems with close-in terrestrial planets

To date, two planetary systems have been discovered with close-in, terrestrial-mass planets     . Many more such discoveries are anticipated in the coming years with radial velocity and transit searches. Here we investigate the different mechanisms that could form 'hot Earths' and their observable predictions. Models include: (1) in situ accretion; (2) formation at larger orbital distance followed by inward 'type 1' migration; (3) formation from material being 'shepherded' inward by a migrating gas giant planet; (4) formation from material being shepherded by moving secular resonances during dispersal of the protoplanetary disc; (5) tidal circularization of eccentric terrestrial planets with close-in perihelion distances and (6) photoevaporative mass-loss of a close-in giant planet. Models 1–4 have been validated in previous work. We show that tidal circularization can form hot Earths, but only for relatively massive planets     with very close-in perihelion distances (≲0.025 au), and even then the net inward movement in orbital distance is at most only 0.1–0.15 au. For planets of less than     , photoevaporation can remove the planet's envelope and leave behind the solid core on a Gyr time-scale, but only for planets inside 0.025–0.05 au. Using two quantities that are observable by current and upcoming missions, we show that these models each produce unique signatures, and can be observationally distinguished. These observables are the planetary system architecture (detectable with radial velocities, transits and transit timing) and the bulk composition of transiting close-in terrestrial planets (measured by transits via the planet's radius).     

Secular and tidal evolution of circumbinary systems

We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin–orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star–planet–satellite systems and triple stellar systems.     

Empirical evidence for tidal evolution in transiting planetary systems

Most transiting planets orbit very close to their parent star, causing strong tidal forces between the two bodies. Tidal interaction can modify the dynamics of the system through orbital alignment, circularization, synchronization and orbital decay by exchange of angular moment. Evidence for tidal circularization in close-in giant planet is well known. Here, we review the evidence for excess rotation of the parent stars due to the pull of tidal forces towards spin-orbit synchronization. We find suggestive empirical evidence for such a process in the present sample of transiting planetary systems. The corresponding angular momentum exchange would imply that some planets have spiralled towards their star by substantial amounts since the dissipation of the protoplanetary disc. We suggest that this could quantitatively account for the observed mass–period relation of close-in gas giants. We discuss how this scenario can be further tested and point out some consequences for theoretical studies of tidal interactions and for the detection and confirmation of transiting planets from radial velocity and photometric surveys.     

Tidal effects in differentiated viscoelastic bodies: a numerical approach

The majority of confirmed terrestrial exoplanets orbits close to their host stars and their evolution was likely altered by tidal interaction. Nevertheless, due to their viscoelastic properties on the tidal frequencies, their response cannot be described exactly by standardly employed constant-lag models. We therefore introduce a tidal model based on the numerical evaluation of a continuum mechanics problem describing the deformation of viscoelastic (Maxwell or Andrade) planetary mantles subjected to external force. We apply the method on a model Earth-size planet orbiting a low-mass star and study the effect of the orbital eccentricity, the mantle viscosity and the chosen rheology on the tidal dissipation, the complex Love numbers and the tidal torque. The number of stable spin states (i.e., zero tidal torque) grows with increasing mantle viscosity, similarly to the analytical model of Correia et al. (Astron Astrophys 571:A50, 2014) for homogeneous bodies. This behavior is only slightly influenced by the rheology used. Similarly, the Love numbers do not distinctly depend on the considered rheological model. The increase in viscosity affects the amplitude of their variations. The tidal heating described by the Maxwell rheology attains local minima associated with low spin-orbit resonances, with depth and shape depending on both the eccentricity and the viscosity. For the Andrade rheology, the minima at low resonances are very shallow and the tidal heating for all viscosities resembles a “fluid limit.” The tidal heating is the quantity influenced the most by the rheology, having thus possible impact on the internal thermal evolution.     

Dynamics of rotation of super-Earths

We numerically investigate the dynamics of rotation of several close-in terrestrial exoplanet candidates. In our model, the rotation of the planet is disturbed by the torque of the central star due to the asymmetric equilibrium figure of the planet. We model the shape of the planet by a Jeans spheroid. We use surfaces of section and spectral analysis to explore numerically the rotation phase space of the systems adopting different sets of parameters and initial conditions close to the main spin–orbit resonant states. One of the parameters, the orbital eccentricity, is critically discussed here within the domain of validity of orbital circularization timescales given by tidal models. We show that, depending on some parameters of the system like the radius and mass of the planet, eccentricity etc., the rotation can be strongly perturbed and a chaotic layer around the synchronous state may occupy a significant region of the phase space. 55 Cnc e is an example.     

Tidal evolution of hierarchical and inclined systems

We investigate the dynamical evolution of hierarchical three-body systems under the effect of tides, when the ratio of the orbital semi-major axes is small and the mutual inclination is relatively large (greater than 20°). Using the quadrupolar non-restricted approximation for the gravitational interactions and the viscous linear model for tides, we derive the averaged equations of motion in a vectorial formalism which is suitable to model the long-term evolution of a large variety of exoplanetary systems in very eccentric and inclined orbits. In particular, it can be used to derive constraints for stellar spin-orbit misalignment, capture in Cassini states, tidal-Kozai migration, or damping of the mutual inclination. Because our model is valid for the non-restricted problem, it can be used to study systems of identical mass or for the outer restricted problem, such as the evolution of a planet around a binary of stars. Here, we apply our model to various situations in the HD 11964, HD 80606, and HD 98800 systems.     

Effect of the rotation and tidal dissipation history of stars on the evolution of close-in planets

Since 20 years, a large population of close-in planets orbiting various classes of low-mass stars (from M-type to A-type stars) has been discovered. In such systems, the dissipation of the kinetic energy of tidal flows in the host star may modify its rotational evolution and shape the orbital architecture of the surrounding planetary system. In this context, recent observational and theoretical works demonstrated that the amplitude of this dissipation can vary over several orders of magnitude as a function of stellar mass, age and rotation. In addition, stellar spin-up occurring during the Pre-Main-Sequence (PMS) phase because of the contraction of stars and their spin-down because of the torque applied by magnetized stellar winds strongly impact angular momentum exchanges within star–planet systems. Therefore, it is now necessary to take into account the structural and rotational evolution of stars when studying the orbital evolution of close-in planets. At the same time, the presence of planets may modify the rotational dynamics of the host stars and as a consequence their evolution, magnetic activity and mixing. In this work, we present the first study of the dynamics of close-in planets of various masses orbiting low-mass stars (from \(0.6~M_\odot \) to \(1.2~M_\odot \)) where we compute the simultaneous evolution of the star’s structure, rotation and tidal dissipation in its external convective envelope. We demonstrate that tidal friction due to the stellar dynamical tide, i.e. tidal inertial waves excited in the convection zone, can be larger by several orders of magnitude than the one of the equilibrium tide currently used in Celestial Mechanics, especially during the PMS phase. Moreover, because of this stronger tidal friction in the star, the orbital migration of the planet is now more pronounced and depends more on the stellar mass, rotation and age. This would very weakly affect the planets in the habitable zone because they are located at orbital distances such that stellar tide-induced migration happens on very long timescales. We also demonstrate that the rotational evolution of host stars is only weakly affected by the presence of planets except for massive companions.     

Motion of dust in mean motion resonances with planets

Effect of stellar electromagnetic radiation on the motion of spherical dust particle in mean motion orbital resonances with a planet is investigated. Planar circular restricted three-body problem with the Poynting–Robertson (P–R) effect yields monotonic secular evolution of eccentricity when the particle is trapped in the resonance. Planar elliptic restricted three-body problem with the P–R effect enables nonmonotonous secular evolution of eccentricity and the evolution of eccentricity is qualitatively consistent with the published results for the complicated case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity of the grain and the orbital evolutions in the resonances are qualitatively equal for the two cases. This holds both for exterior and interior mean motion orbital resonances. Evolutions of argument of perihelion in the planar circular and elliptical restricted three-body problems are shown. Numerical integrations show that an analytic expression for the secular time derivative of the particle’s argument of perihelion does not exist, if only dependence on semimajor axis, eccentricity and argument of perihelion is admitted. Connection between the shift of perihelion and oscillations in secular eccentricity is presented for the planar elliptic restricted three-body problem with the P–R effect. Period of the oscillations corresponds to the period of one revolution of perihelion. Change of optical properties of the spherical grain with the heliocentric distance is also considered. The change of the optical properties: (i) does not have any significant influence on the secular evolution of eccentricity, (ii) causes that the shift of perihelion is mainly in the same direction/orientation as the particle motion around the Sun. The statements hold both for circular and noncircular planetary orbits.     

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.     

On the stability of dust orbits in mean-motion resonances perturbed by from an interstellar wind

Circumstellar dust particles can be captured in a mean-motion resonance (MMR) with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the MMR analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange’s planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting–Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which holds for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.     

On the early history of the asteroids

Among the major features of the asteroids as a group is the fact that they are small and numerous rather than being a single planet, and that they have unusually high eccentricities and inclinations. Regarding the first, this paper presents two lines of argument concerning the concept that the asteroids formed with a mass-distribution similar to what we observe today. Considerations of planet accretion are used to clarify the role of Jupiter in preventing coalescence of asteroids into a single planet. Regarding the eccentricities and inclinations, it is proposed that they were excited by secular resonances associated with the presence of Jupiter within a dissipating solar nebula.     

Secular spin dynamics of inner main-belt asteroids

Understanding the evolution of asteroid spin states is challenging work, in part because asteroids have a variety of orbits, shapes, spin states, and collisional histories but also because they are strongly influenced by gravitational and non-gravitational (YORP) torques. Using efficient numerical models designed to investigate asteroid orbit and spin dynamics, we study here how several individual asteroids have had their spin states modified over time in response to these torques (i.e., 951 Gaspra, 60 Echo, 32 Pomona, 230 Athamantis, 105 Artemis). These test cases which sample semimajor axis and inclination space in the inner main belt, were chosen as probes into the large parameter space described above. The ultimate goal is to use these data to statistically characterize how all asteroids in the main belt population have reached their present-day spin states. We found that the spin dynamics of prograde-rotating asteroids in the inner main belt is generally less regular than that of the retrograde-rotating ones because of numerous overlapping secular spin-orbit resonances. These resonances strongly affect the spin histories of all bodies, while those of small asteroids (?40 km) are additionally influenced by YORP torques. In most cases, gravitational and non-gravitational torques cause asteroid spin axis orientations to vary widely over short (?1 My) timescales. Our results show that (951) Gaspra has a highly chaotic rotation state induced by an overlap of the s and s₆ spin-orbit resonances. This hinders our ability to investigate its past evolution and infer whether thermal torques have acted on Gaspra's spin axis since its origin.     

A note on the tidal expansion of the orbit of Mercury

Under the assumption that the secular decrease in the sidereal mean motion of Mercury is due to the Sun-Mercury tidal torque, the expansion of the orbit of Mercury and the tidal phase-lag angle of the tidal bulge are estimated.     

Ringlet–satellite interactions

We study the interaction of a satellite and a nearby ringlet on eccentric and inclined orbits. Secular torques originate from mean motion resonances and the secular interaction potential which represents the m  = 1 global modes of the ring. The torques act on the relative eccentricity and inclination. The resonances damp the relative eccentricity. The inclination instability owing to the resonances is turned off by a finite differential eccentricity of the order of 0.27 for nearly coplanar systems. The secular potential torque damps the eccentricity and inclination and does not affect the relative semi-major axis; also, it suppresses the inclination instability that persists at small differential eccentricities. The damping of the relative eccentricity and inclination forces an initially circular and planar small mass ringlet to reach the eccentricity and inclination of the satellite. When the planet is oblate, the interaction of the satellite damps the proper precession of a small mass ringlet so that it precesses at the satellite's rate independently of their relative distance. The oblateness of the primary modifies the long-term eccentricity and inclination magnitudes and introduces a constant shift in the apsidal and nodal lines of the ringlet with respect to those of the satellite. These results are applied to Saturn's F-ring, which orbits between the moons Prometheus and Pandora.     

Normal Form Invariants Around Spin-orbit Periodic Orbits

We consider a model of spin-orbit interaction, describing the motion of an oblate satellite rotating about an internal spin-axis and orbiting about a central planet. The resulting second order differential equation depends upon the parameters provided by the equatorial oblateness of the satellite and its orbital eccentricity. Normal form transformations around the main spin-orbit resonances are carried out explicitly. As an outcome, one can compute some invariants; the fact that these quantities are not identically zero is a necessary condition to prove the existence of nearby periodic orbits (Birkhoff fixed point theorem). Moreover, the nonvanishing of the invariants provides also the stability of the spin-orbit resonances, since it guarantees the existence of invariant curves surrounding the periodic orbit.     

Mean-Motion Resonances as a Source for Infalling Comets toward β Pictoris

Repeated time-variable redshifted absorption features in the spectrum of β Pictoris (β Pic) have been attributed to comet-like bodies falling toward the star, when evaporating in its immediate vicinity. This model explains now a large number of observational characteristics, but the exact mechanism that could generate these numerous star-grazers is still controversial, even if planetary perturbations are thought to be the basic process. The different models proposed up to now are here reviewed, and we discuss in particular a recent one, involving the effect of secular resonances in the β Pic system. We stress that it seems highly improbable that such a mechanism could apply to the β Pic case, because the extremely strong power of secular resonances is connected to the very specific structure of the Solar System. Therefore, the secular resonance mechanism is highly non-generic. Conversely, we propose a model involving the eccentricity-pumping effect of mean-motion resonances with a massive planet on a moderately eccentric orbit. We show in particular that the 4:1 mean-motion resonance is a very active source of star-grazers as soon as the eccentricity of the perturbing planet is ?0.05, while the 3:1 mean-motion resonance is less efficient. We stress that this mechanism is very generic. These theoretical predictions are confirmed by numerical integrations using the Extended Schubart Integrator. The time-scale of the process is discussed, and we show that if the eccentricity of the perturbing planet fluctuates, due to secular perturbations, this time-scale is compatible with the age of β Pic's system.     

The secular variation of pulsar magnetic dipole moments

The time dependences of the inertia tensor and of a dissipative torque caused by the nonleptonic weak interaction have been investigated for a certain class of pulsars with no solid core. Early in the life of the pulsar, the angular velocity vector is predicted to move with respect to fixed body axes in such a way that it becomes perpendicular to the magnetic dipole moment. During this motion, the solid outer shell suffers plastic deformation so that the dipole moment becomes approximately collinear with a principal axis. After 10⁴ or 10⁵ yr, the dissipative torque is negligibly small compared with the electromagnetic torque, the Euler equations are those for a simple rigid body, and alignment of spin and dipole moment occurs. If the dipole moment discussed by Lyneet al. (1975) is interpreted as being equal to the component perpendicular to the spin, its secular decay is a natural property of this model and is not a consequence of field decay through electrical resistivity.     

Obliquity variations of terrestrial planets in habitable zones

We have investigated obliquity variations of possible terrestrial planets in habitable zones (HZs) perturbed by a giant planet(s) in extrasolar planetary systems. All the extrasolar planets so far discovered are inferred to be jovian-type gas giants. However, terrestrial planets could also exist in extrasolar planetary systems. In order for life, in particular for land-based life, to evolve and survive on a possible terrestrial planet in an HZ, small obliquity variations of the planet may be required in addition to its orbital stability, because large obliquity variations would cause significant climate change. It is known that large obliquity variations are caused by spin-orbit resonances where the precession frequency of the planet's spin nearly coincides with one of the precession frequencies of the ascending node of the planet's orbit. Using analytical expressions, we evaluated the obliquity variations of terrestrial planets with prograde spins in HZs. We found that the obliquity of terrestrial planets suffers large variations when the giant planet's orbit is separated by several Hill radii from an edge of the HZ, in which the orbits of the terrestrial planets in the HZ are marginally stable. Applying these results to the known extrasolar planetary systems, we found that about half of these systems can have terrestrial planets with small obliquity variations (smaller than 10°) over their entire HZs. However, the systems with both small obliquity variations and stable orbits in their HZs are only 1/5 of known systems. Most such systems are comprised of short-period giant planets. If additional planets are found in the known planetary systems, they generally tend to enhance the obliquity variations. On the other hand, if a large/close satellite exists, it significantly enhances the precession rate of the spin axis of a terrestrial planet and is likely to reduce the obliquity variations of the planet. Moreover, if a terrestrial planet is in a retrograde spin state, the spin-orbit resonance does not occur. Retrograde spin, or a large/close satellite might be essential for land-based life to survive on a terrestrial planet in an HZ.