We present a model to study secularly and tidally evolving three-body systems composed by two low-mass planets orbiting a star, in the case where the bodies rotation axes are always perpendicular to the orbital plane. The tidal theory allows us to study the spin and orbit evolution of both stiff Earth-like planets and predominantly gaseous Neptune-like planets. The model is applied to study two recently-discovered exoplanetary systems containing potentially habitable exoplanets (PHE): LHS-1140 b-c and K2-18 b-c. For the former system, we show that both LHS-1140 b and c must be in nearly-circular orbits. For K2-18 b-c, the combined analysis of orbital evolution timescales with the current eccentricity estimation of K2-18 b allows us to conclude that the inner planet (K2-18 c) must be a Neptune-like gaseous body. Only this would allow for the eccentricity of K2-18 b to be in the range of values estimated in recent works ($e=0.20 \pm 0.08$), provided that the uniform viscosity coefficient of K2-18 b is greater than $2.4 \times 10^{19} \ \textrm{Pa s}$ (which is a value characteristic of stiff bodies) and supposing that such system has an age of some Gyr.

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.