Paper on the creep tide theory and its applications to satellites and planets with emphasis on a new set of differential equations allowing easier numerical studies. The creep tide theory is a new paradigm that does not fix a priori the tidal deformation of the body, but considers the deformation as a low-Reynolds-number flow. The evolution under tidal forces is ruled by an approximate solution of the Navier-Stokes equation depending on the body's viscosity with no ad hoc assumptions on its shape and orientation. It reproduces closely the results of Darwinian theories in the case of gaseous planets and stars, but the results are completely different in the case of stiff satellites and planets. It explains the tidal dissipations of Enceladus and Mimas. The extension of the theory to nonhomogeneous icy satellites with a subsurface ocean allows the amplitude of the forced oscillations around synchronization (librations) to be better determined.

This paper deals with the rotation and figure evolution of a planet near the 3/2 spin-orbit resonance and the exploration of a new formulation of the creep tide theory (Folonier et al. 2018). This new formulation is composed by a system of differential equations for the figure and the rotation of the body simultaneously (which is the same system of equations used in Folonier et al. 2018), different from the original one (Ferraz-Mello, 2013, 2015a) in which rotation and figure were considered separately. The time evolution of the figure of the body is studied for both the 3/2 and 2/1 spin-orbit resonances. Moreover, we provide a method to determine the relaxation factor gamma of non-rigid homogeneous bodies whose endpoint of rotational evolution from tidal interactions is the 3/2 spin-orbit resonance, provided that (i) an initially faster rotation is assumed and (ii) no permanent components of the flattenings of the body existed at the time of the capture in the 3/2 spin-orbit resonance. The method is applied to Mercury, since it is currently trapped in a 3/2 spin-orbit resonance with its orbital motion and we obtain 4.8 times 10 -8 s -1 lower than gamma lower than 4.8 times 10 -9 s -1 . The equatorial prolateness and polar oblateness coefficients obtained for Mercury's figure with such range of values of gamma are the same as the ones given by the Darwin-Kaula model (Matsuyama and Nimmo 2009). However, comparing the values of the flattenings obtained for such range of gamma with those obtained from MESSENGER's measurements (Perry et al. 2015), we see that the current values for Mercury's equatorial prolateness and polar oblateness are 2-3 orders of magnitude larger than the values given by the tidal theories.

We present a self-consistent model for the tidal evolution of circumbinary planets. Based on the weak-friction model, we derive expressions of the resulting forces and torques considering complete tidal interactions between all the bodies of the system. Although the tidal deformation suffered by each extended mass must take into account the combined gravitational effects of the other two bodies, the only tidal forces that have a net effect on the dynamic are those that are applied on the same body that exerts the deformation, as long as no mean-motion resonance exists between the masses. We apply the model to the Kepler-38 binary system. The evolution of the spin equations shows that the planet reaches a stationary solution much faster than the stars, and the equilibrium spin frequency is sub-synchronous. The binary components evolve on a longer timescale, reaching a super-synchronous solution very close to that derived for the 2-body problem. After reaching spin stationarity, the eccentricity is damped in all bodies and for all the parameters analyzed here. A similar effect is noted for the binary separation. The semimajor axis of the planet, on the other hand, may migrate inwards or outwards, depending on the masses and orbital parameters. In some cases the secular evolution of the system may also exhibit an alignment of the pericenters, requiring to include additional terms in the tidal model. Finally, we derived analytical expressions for the variational equations of the orbital evolution and spin rates based on low-order elliptical expansions in the semimajor axis ratio and the eccentricities. These are found to reduce to the 2-body case when one of the masses is taken equal to zero. This model allow us to find a close and simple analytical expression for the stationary spin rates of all the bodies, as well as predicting the direction and magnitude of the orbital migration.

This paper deals with a new formulation of the creep tide theory (Ferraz-Mello, Cel. Mech. Dyn. Astron. {\bf 116}, 109, 2013 $-$ Paper I) and with the tidal dissipation predicted by the theory in the case of stiff bodies whose rotation is not synchronous but is oscillating around the synchronous state with a period equal to the orbital period. We show that the tidally forced libration influences the amount of energy dissipated in the body and the average perturbation of the orbital elements. This influence depends on the libration amplitude and is generally neglected in the study of planetary satellites. However, they may be responsible for a 27 percent increase in the dissipation of Enceladus. The relaxation factor necessary to explain the observed dissipation of Enceladus ($\gamma=1.2-3.8\times 10^{-7}\ {\rm s}^{-1}$) has the expected order of magnitude for planetary satellites and corresponds to the viscosity $0.6-1.9 \times 10^{14}$ Pa s, which is in reasonable agreement with the value recently estimated by Efroimsky (2018) ($0.24 \times 10^{14}$ Pa s) and with the value adopted by Roberts and Nimmo (2008) for the viscosity of the ice shell ($10^{13}-10^{14}$ Pa s). For comparison purposes, the results are extended also to the case of Mimas and are consistent with the negligible dissipation and the absence of observed tectonic activity. The corrections of some mistakes and typos of paper II (Ferraz-Mello, Cel. Mech. Dyn. Astron. {\bf 122}, 359, 2015) are included at the end of the paper.

This paper presents one analytical tidal theory for a viscoelastic multi-layered body with an arbitrary number of homogeneous layers. Starting with the static equilibrium figure, modified to include tide and differential rotation, and using the Newtonian creep approach, we find the dynamical equilibrium figure of the deformed body, which allows us to calculate the tidal potential and the forces acting on the tide generating body, as well as the rotation and orbital elements variations. In the particular case of the two-layer model, we study the tidal synchronization when the gravitational coupling and the friction in the interface between the layers is added. For high relaxation factors (low viscosity), the stationary solution of each layer is synchronous with the orbital mean motion (n) when the orbit is circular, but the spin rates increase if the orbital eccentricity increases. For low relaxation factors (high viscosity), as in planetary satellites, if friction remains low, each layer can be trapped in different spin-orbit resonances with frequencies n/2,n,3n/2,... . We apply the theory to Titan. The main results are: i) the rotational constraint does not allow us confirm or reject the existence of a subsurface ocean in Titan; and ii) the crust-atmosphere exchange of angular momentum can be neglected. Using the rotation estimate based on Cassini's observation, we limit the possible value of the shell relaxation factor, when a subsurface ocean is assumed, to 10^-9 Hz, which correspond to a shell's viscosity 10^18 Pa s, depending on the ocean's thickness and viscosity values. In the case in which the ocean does not exist, the maximum shell relaxation factor is one order of magnitude smaller and the corresponding minimum shell's viscosity is one order higher.

This paper deals with the application of the creep tide theory (Ferraz-Mello, Cel. Mech. Dyn. Astron. vol. 116, 109, 2013) to the study of the rotation of stars hosting massive close-in planets. The stars have nearly the same tidal relaxation factors as gaseous planets and the evolution of their rotation is similar to that of close-in hot Jupiters: they tidally evolve towards a stationary solution. However, stellar rotation may also be affected by stellar wind braking. Thus, while the rotation of a quiet host star evolves towards a stationary attractor with a frequency ($1+6e^2$) times the orbital mean-motion of the companion, the continuous loss of angular momentum in an active star displaces the stationary solution towards slower values: Active host stars with big close-in companions tend to have rotational periods larger than the orbital periods of their companions. The study of some hypothetical examples shows that because of tidal evolution, the rules of gyrochronology cannot be used to estimate the age of one system with a large close-in companion, no matter if the star is quiet or active, if the current semi-major axis of the companion is smaller than 0.03--0.04 AU. Details on the evolution of the systems: CoRoT LRc06E21637, CoRoT-27, Kepler-75, CoRoT-2, CoRoT-18, CoRoT-14 and on hypothetical systems with planets of mass 1--4 M_Jup in orbit around a star similar to the Sun are given.

We consider the Clairaut theory of the equilibrium ellipsoidal figures for differentiated non-homogeneous bodies in non-synchronous rotation adding to it a tidal deformation due to the presence of an external gravitational force. We assume that the body is a fluid formed by $n$ homogeneous layers of ellipsoidal shape and we calculate the external polar flattenings and the mean radius of each layer, or, equivalently, their semiaxes. To first order in the flattenings, the general solution can be written as $\epsilon_k={\cal H}_k*\epsilon_h$ and $\mu_k={\cal H}_k*\mu_h$, where $\cal{H}_k$ is a characteristic coefficient for each layer which only depends on the internal structure of the body and $\epsilon_h, \mu_h$ are the flattenings of the equivalent homogeneous problem. For the continuous case, we study the Clairaut differential equation for the flattening profile, using the Radau transformation to find the boundary conditions when the tidal potential is added. Finally, the theory is applied to several examples: i) a body composed of two homogeneous layers; ii) bodies with simple polynomial density distribution laws and iii) bodies following a polytropic pressure-density law.

We study the capture and crossing probabilities into the 3:1 mean motion resonance with Jupiter for a small asteroid that migrates from the inner to the middle Main Belt under the action of the Yarkovsky effect. We use an algebraic mapping of the averaged planar restricted three-body problem based on the symplectic mapping of Hadjidemetriou (1993), adding the secular variations of the orbit of Jupiter and non-symplectic terms to simulate the migration. We found that, for fast migration rates, the captures occur at discrete windows of initial eccentricities whose specific locations depend on the initial resonant angles, indicating that the capture phenomenon is not probabilistic. For slow migration rates, these windows become narrower and start to accumulate at low eccentricities, generating a region of mutual overlap where the capture probability tends to 100%, in agreement with the theoretical predictions for the adiabatic regime. Our simulations allow to predict the capture probabilities in both the adiabatic and non-adiabatic cases, in good agreement with results of Gomes (1995) and Quillen (2006). We apply our model to the case of the Vesta asteroid family in the same context as Roig et al. (2008), and found results indicating that the high capture probability of Vesta family members into the 3:1 mean motion resonance is basically governed by the eccentricity of Jupiter and its secular variations.