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.

We present an extended version of the Constant Time Lag analytical approach for the tidal evolution of circumbinary planets introduced in our previous work. The model is self-consistent, in the sense that all tidal interactions between pairs are computed, regardless of their size. We derive analytical expressions for the variational equations governing the spin and orbital evolution, which are expressed as high-order elliptical expansions in the semimajor axis ratio but retain closed form in terms of the binary and planetary eccentricities. These are found to reproduce the results of the numerical simulations with arbitrary eccentricities very well, as well as reducing to our previous results in the low-eccentric case. Our model is then applied to the well-characterised Kepler circumbinary systems by analysing the tidal timescales and unveiling the tidal flow around each different system. In all cases we find that the spins reach stationary values much faster than the characteristic timescale of the orbital evolution, indicating that all Kepler circumbinary planets are expected to be in a sub-synchronous state. On the other hand, all systems are located in a tidal flow leading to outward migration; thus the proximity of the planets to the orbital instability limit may have been even greater in the past. Additionally, Kepler systems may have suffered a significant tidally induced eccentricity damping, which may be related to their proximity to the capture eccentricity. To help understand the predictions of our model, we also offer a simple geometrical interpretation of our results.

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.

Based on the model described in Ramos et al., 2017, we present an analytical+numerical study of the resonance capture under Type-I migration for the Kepler-25 (Marcy et al., 2014) and K2-24 (Petigura et al., 2016) Kepler systems, both close to a 2/1 mean-motion resonance. We find that, depending on the flare index and the proximity to the central star, the average value of the period-ratio between two consecutive planets show a significant deviation with respect to the resonant nominal value, up to values well in agreement with the observations.

We present a series of dynamical maps for fictitious 3-planets systems in initially circular coplanar orbits. These maps have unveiled a rich resonant structure involving two or three planets, as well as indicating possible migration routes from secular to double resonances or pure 3-planet commensurabilities. These structures are then compared to the present-day orbital architecture of observed resonant chains. In a second part of the paper we describe N-body simulations of type-I migration. Depending on the orbital decay timescale, we show that 3-planet systems may be trapped in different combinations of independent commensurabilities: (i) double resonances, (ii) intersection between a 2-planet and a first-order 3-planet resonance, and (iii) simultaneous libration in two first-order 3-planet resonances. These latter outcomes are found for slow migrations, while double resonances are almost always the final outcome in high-density disks. Finally, we discuss an application to the TRAPPIST-1 system. We find that, for low migration rates and planetary masses of the order of the estimated values, most 3-planet sub-systems are able to reach the observed double resonances after following evolutionary routes defined by pure 3-planet resonances. The final orbital configuration shows resonance offsets comparable with present-day values without the need of tidal dissipation. For the 8/5 resonance proposed to dominate the dynamics of the two inner planets, we find little evidence of its dynamical significance; instead, we propose that this relation between mean motions could be a consequence of the interaction between a pure 3-planet resonance and a 2-planet commensurability between planets c and d.

We present an analytical and numerical study of the orbital migration and resonance capture of fictitious two-planet systems with masses in the super-Earth range undergoing Type-I migration. We find that, depending on the flare index and proximity to the central star, the average value of the period ratio, $P_2/P_1$, between both planets may show a significant deviation with respect to the nominal value. For planets trapped in the 2:1 commensurability, offsets may reach values on the order of $0.1$ for orbital periods on the order of $1$ day, while systems in the 3:2 mean-motion resonance (MMR) show much smaller offsets for all values of the semimajor axis. These properties are in good agreement with the observed distribution of near-resonant exoplanets, independent of their detection method. We show that 2:1-resonant systems far from the star, such as HD82943 and HR8799, are characterized by very small resonant offsets, while higher values are typical of systems discovered by Kepler with orbital periods approximately a few days. Conversely, planetary systems in the vicinity of the 3:2 MMR show little offset with no significant dependence on the orbital distance. In conclusion, our results indicate that the distribution of Kepler planetary systems around the 2:1 and 3:2 MMR are consistent with resonant configurations obtained as a consequence of a smooth migration in a laminar flared disk, and no external forces are required to induce the observed offset or its dependence with the commensurability or orbital distance from the star.

We present a new method that allows long-term and large-scale hydrodynamical simulations of migrating planets over a grid-based Eulerian code. This technique, which consists in a remapping of the disk by tracking the planetary migration, enables runs of migrating planets over a time comparable to the age of protoplanetary disks. This method also has the potential to address efficiently problems related with migration of multi-planet systems in gaseous disks, and to improve current results of migration of massive planets by including global viscous evolution as well as detailed studies of the co-orbital region during migration. We perform different tests using the public code FARGO3D to validate this method and compare its results with those obtained using a classical fixed grid.

Chaotic diffusion is supposed to be responsible for orbital instabilities in planetary systems after the dissipation of the protoplanetary disk, and a natural consequence of irregular motion. In this paper we show that resonant multi-planetary systems, despite being highly chaotic, not necessarily exhibit significant diffusion in phase space, and may still survive virtually unchanged over timescales comparable to their age.Using the GJ-876 system as an example, we analyze the chaotic diffusion of the outermost (and less massive) planet. We construct a set of stability maps in the surrounding regions of the Laplace resonance. We numerically integrate ensembles of close initial conditions, compute Poincar\'e maps and estimate the chaotic diffusion present in this system. Our results show that, the Laplace resonance contains two different regions: an inner domain characterized by low chaoticity and slow diffusion, and an outer one displaying larger values of dynamical indicators. In the outer resonant domain, the stochastic borders of the Laplace resonance seem to prevent the complete destruction of the system. We characterize the diffusion for small ensembles along the parameters of the outermost planet. Finally, we perform a stability analysis of the inherent chaotic, albeit stable Laplace resonance, by linking the behavior of the resonant variables of the configurations to the different sub-structures inside the three-body resonance.

We present a series of numerical integrations of observed and fictitious Jupiter Trojan asteroids, under the gravitational effects of the four outer planets, for time-spans comparable with the age of the Solar System. From these results we calculate the escape rate from each Lagrange point, and construct dynamical maps of "permanence" time in different regions of the phase space. Fictitious asteroids in L4 and L5 show no significant difference, showing almost identical dynamical maps and escape rates. For real Trojans, however, we found that approximately 23% o f the members of the leading swarm escaped after 4.5 Gyrs, while this number increased to 28.3% for L5. This implies that the asymmetry between the two populations increases with time, indicating that it may have been smaller at the time of formation/capture of these asteroids. Nevertheless, the difference in chaotic diffusion cannot, in itself, account for the current observed asymmetry (~40%), and must be primarily primordial and characteristic of the capture mechanism of the Trojans. Finally, we calculate new proper elements for all the numbered Trojans using the semi-analytical approach of Beaug\'e and Roig (2001), and compare the results with the numerical estimations by Bro\v{z} and Rosehnal (2011). For asteroids that were already numbered in 2011, both methods yield very similar results, while significant differences were found for those bodies that became numbered after 2011.

We analyse the secular dynamics of planets on S-type coplanar orbits in tight binary systems, based on first- and second-order analytical models, and compare their predictions with full N-body simulations. The perturbation parameter adopted for the development of these models depends on the masses of the stars and on the semimajor axis ratio between the planet and the binary. We show that each model has both advantages and limitations. While the first-order analytical model is algebraically simple and easy to implement, it is only applicable in regions of the parameter space where the perturbations are sufficiently small. The second-order model, although more complex, has a larger range of validity and must be taken into account for dynamical studies of some real exoplanetary systems such as $\gamma$-Cephei and HD 41004A. However, in some extreme cases, neither of these analytical models yields quantitatively correct results, requiring either higher-order theories or direct numerical simulations. Finally, we determine the limits of applicability of each analytical model in the parameter space of the system, giving an important visual aid to decode which secular theory should be adopted for any given planetary system in a close binary.

We review the orbital stability of the planar circular restricted three-body problem, in the case of massless particles initially located between both massive bodies. We present new estimates of the resonance overlap criterion and the Hill stability limit, and compare their predictions with detailed dynamical maps constructed with N-body simulations. We show that the boundary between (Hill) stable and unstable orbits is not smooth but characterized by a rich structure generated by the superposition of different mean-motion resonances which does not allow for a simple global expression for stability. We propose that, for a given perturbing mass $m_1$ and initial eccentricity $e$, there are actually two critical values of the semimajor axis. All values $a < a_{\rm Hill}$ are Hill-stable, while all values $a > a_{\rm unstable}$ are unstable in the Hill sense. The first limit is given by the Hill-stability criterion and is a function of the eccentricity. The second limit is virtually insensitive to the initial eccentricity, and closely resembles a new resonance overlap condition (for circular orbits) developed in terms of the intersection between first and second-order mean-motion resonances.

HD 82943 hosts a mysterious multi-planet system in the 2:1 mean-motion resonance that puzzles astronomers for more than a decade. We describe our new analysis of all radial velocity data currently available for this star, including both the most recent Keck data and the older but more numerous CORALIE measurements. Here we pay a major attention to the task of optimal scheduling of the future observation of this system. Applying several optimality criteria, we demonstrate that in the forthcoming observational season of HD 82943 (the winter 2014/2015) rather promising time ranges can be found. Observations of the near future may give rather remarkable improvement of the orbital fit, but only if we choose their time carefully.

Hot Jupiters (HJs) are usually defined as giant Jovian-size planets with orbital periods $P \le 10$ days. Although they lie close to the star, several have finite eccentricities and significant misalignment angle with respect to the stellar equator. Two mechanisms have been proposed to explain the excited and misaligned sub-population of HJs: Lidov-Kozai migration and planet-planet scattering. Although both are based on completely different dynamical phenomena, they appear to be equally effective in generating hot planets. Nevertheless, there has been no detailed analysis comparing the predictions of both mechanisms. In this paper we present numerical simulations of Lidov-Kozai trapping of single planets in compact binary systems. Both the planet and the binary are initially placed in coplanar orbits, although the inclination of the impactor is assumed random. After the passage of the third star, we follow the orbital and spin evolution of the planet using analytical models based on the octupole expansion of the secular Hamiltonian. The present work aims at the comparison of the two mechanisms, as an explanation for the excited and inclined HJs in binary systems. We compare the results obtained through this paper with results in Beaug\'e & Nesvorn\'y 2012, where the authors analyze how the planet-planet scattering mechanisms works. Several of the orbital characteristics of the simulated HJs are caused by tidal trapping from quasi-parabolic orbits, independent of the driving mechanism. These include both the 3-day pile-up and the distribution in the eccentricity vs semimajor axis plane. However, the distribution of the inclinations shows significant differences. While Lidov-Kozai trapping favors a more random distribution, planet-planet scattering shows a large portion of bodies nearly aligned with the equator of the central star.

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.

We carry out a new analysis of the published radial velocity data for the planet-hosting star HD82943. We include the recent Keck/HIRES measurements as well as the aged but much more numerous CORALIE data. We find that the CORALIE radial velocity measurements are polluted by a systematic annual variation which affected the robustness of many previous results. We show that after purging this variation, the residuals still contain a clear signature of an additional $\sim 1100$ days periodicity. The latter variation leaves significant hints in all three independent radial velocity subsets that we analysed: the CORALIE data, the Keck data acquired prior to a hardware upgrade, and the Keck data taken after the upgrade. We mainly treat this variation as a signature of a third planet in the system, although we cannot rule out other interpretations, such as long-term stellar activity. We find it easy to naturally obtain a stable three-planet radial-velocity fit close to the three-planet mean-motion resonance 1:2:5, with the two main planets (those in the 1:2 resonance) in an aligned apsidal corotation. The dynamical status of the third planet is still uncertain: it may reside in as well as slightly out of the 5:2 resonance. We obtain the value of $\sim 1075$ days for its orbital period and of $\sim 0.3 M_{\rm Jup}$ for its minimum mass, while the eccentric parameters are uncertain.

In this paper we revise the model proposed by Rein et al. (2010) for the origin of the HD45364 exoplanetary system, currently known to host two planets close to the 3/2 mean-motion commensurability (MMR). We show that, due to high surface density of the protoplanetary disk needed for Type III migration, this model could only lead to planets in a quasi-resonant regime of motion, and thus not consistent with the resonant configuration obtained by Correia et al. (2009). Although both resonant and quasi-resonant solutions are statistically indistinguishable with respect to radial velocity measurements, their distinct dynamical behaviour is intriguing. We use the semi-analytical model to confirm the quantitative difference between two configurations. In order to form a system evolving inside the 3/2 resonance, we develop a different model. Our scenario includes an interaction between different (but slower) planetary migration types, planet growth, and gap formation in the protoplanetary disk. The choice of the described evolutionary path was due to a detailed analysis of the structure of the phase space in the vicinity of the 3/2 MMR employing dynamical mapping techniques. The outcomes of our simulations are able to reproduce very closely the 3/2 resonant dynamics obtained from the best-fit presented by Correia et al. (2009). In addition, varying the strength of the eccentricity damping, we can also simulate the quasi-resonant configuration similar to that in Rein et al. (2010). We furthermore show that our scenario is robust with respect to the physical parameters involved in the resonance trapping process. However, the confirmation of our scenario will be possible only with additional radial velocities measurements.

We present an algebraic map (MAMA) for the dynamical and collisional evolution of a planetesimal swarm orbiting the main star of a tight binary system (TBS). The orbital evolution of each planetesimal is dictated by the secular perturbations of the secondary star and gas drag due to interactions with a protoplanetary disk. The gas disk is assumed eccentric with a constant precession rate. Gravitational interactions between the planetesimals are ignored. All bodies are assumed coplanar. A comparison with full N-body simulations shows that the map is of the order of 100 times faster, while preserving all the main characteristics of the full system. In a second part of the work, we apply MAMA to the \gamma-Cephei, searching for friendly scenarios that may explain the formation of the giant planet detected in this system. For low-mass protoplanetary disks, we find that a low-eccentricity static disk aligned with the binary yields impact velocities between planetesimals below the disruption threshold. All other scenarios appear hostile to planetary formation.

The motion of two planets around a Sun-like star under the combined effects of mutual interaction and tidal dissipation is investigated. The secular behaviour of the system is analyzed using two different approaches. First, we solve the exact equations of motion through the numerical simulation of the system evolution. In addition to the orbital decay and circularization, we show that the final configuration of the system is affected by the shrink of the inner orbit. Our second approach consist in the analysis of the stationary solutions of mean equations of motion based on a Hamiltonian formalism. We consider the case of a hot super-Earth planet with a more massive outer companion. As a real example, the CoRoT-7 system is analyzed solving the exact and mean equations of motion. The star-planet tidal interaction produces orbital decay and circularization of the orbit of CoRoT-7b. In addition, the long-term tidal evolution is such that the eccentricity of CoRoT-7c is also circularized and a pair of final circular orbits is obtained. A curve in the space of eccentricities can be constructed through the computation of stationary solutions of mean equations including dissipation. The application to CoRoT-7 system shows that the stationary curve agrees with the result of numerical simulations of exact equations. A similar investigation performed in a super-Earth-Jupiter two-planet system shows that the doubly circular state is accelerated when there is a significant orbital migration of the inner planet, in comparison with previous results were migration is neglected.

We study the stability regions and families of periodic orbits of two planets locked in a co-orbital configuration. We consider different ratios of planetary masses and orbital eccentricities, also we assume that both planets share the same orbital plane. Initially we perform numerical simulations over a grid of osculating initial conditions to map the regions of stable/chaotic motion and identify equilibrium solutions. These results are later analyzed in more detail using a semi-analytical model. Apart from the well known quasi-satellite (QS) orbits and the classical equilibrium Lagrangian points L4 and L5, we also find a new regime of asymmetric periodic solutions. For low eccentricities these are located at $(\sigma,\Delta\omega) = (\pm 60\deg, \mp 120\deg)$, where \sigma is the difference in mean longitudes and \Delta\omega is the difference in longitudes of pericenter. The position of these Anti-Lagrangian solutions changes with the mass ratio and the orbital eccentricities, and are found for eccentricities as high as ~ 0.7. Finally, we also applied a slow mass variation to one of the planets, and analyzed its effect on an initially asymmetric periodic orbit. We found that the resonant solution is preserved as long as the mass variation is adiabatic, with practically no change in the equilibrium values of the angles.

Constructing dynamical maps from the filtered output of numerical integrations, we analyze the structure of the $\nu_\odot$ secular resonance for fictitious irregular satellites in retrograde orbits. This commensurability is associated to the secular angle $\theta = \varpi - \varpi_\odot$, where $\varpi$ is the longitude of pericenter of the satellite and $\varpi_\odot$ corresponds to the (fixed) planetocentric orbit of the Sun. Our study is performed in the restricted three-body problem, where the satellites are considered as massless particles around a massive planet and perturbed by the Sun. Depending on the initial conditions, the resonance presents a diversity of possible resonant modes, including librations of $\theta$ around zero (as found for Sinope and Pasiphae) or 180 degrees, as well as asymmetric librations (e.g. Narvi). Symmetric modes are present in all giant planets, although each regime appears restricted to certain values of the satellite inclination. Asymmetric solutions, on the other hand, seem absent around Neptune due to its almost circular heliocentric orbit. Simulating the effects of a smooth orbital migration on the satellite, we find that the resonance lock is preserved as long as the induced change in semimajor axis is much slower compared to the period of the resonant angle (adiabatic limit). However, the librational mode may vary during the process, switching between symmetric and asymmetric oscillations. Finally, we present a simple scaling transformation that allows to estimate the resonant structure around any giant planet from the results calculated around a single primary mass.

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than $\sim 4$. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are presently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically over estimated, although the inner planet's eccentricity suffers a much smaller effect. If the initial conditions correspond to small amplitude oscillations around stable apsidal corotation resonances (ACR), the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

The published orbits of the planets HD 82943b and HD 82943c correspond to a system bound to a catastrophic event in less than 100,000 years. Alternative sets of elements and masses, which fit the available observational data and correspond to regular motions, are presented in this paper. The planets HD 82943 c,b are in a 2/1 mean-motion resonance and are trapped in the neighborhood of a (0,0)-apsidal corotation.