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.

In this work, we construct and test an analytical and a semianalytical secular models for two planets locked in a coorbital non-coplanar motion, comparing some results with the case of restricted three body problem. The analytical average model replicates the numerical N-body integrations, even for moderate eccentricities ($\lesssim$ 0.3) and inclinations ($\lesssim10^\circ$), except for the regions corresponding to quasi-satellite and Lidov-Kozai configurations. Furthermore, this model is also useful in the restricted three body problem, assuming very low mass ratio between the planets. We also describe a four-degree-of-freedom semianalytical model valid for any type of coorbital configuration in a wide range of eccentricities and inclinations. {Using a N-body integrator, we have found that the phase space of the General Three Body Problem is different to the restricted case for inclined systems, and establish the location of the Lidov-Kozai equilibrium configurations depending on mass ratio. We study the stability of periodic orbits in the inclined systems, and find that apart from the robust configurations $L_4$, $AL_4$, and $QS$ is possible to harbour two Earth-like planets in orbits previously identified as unstable $U$ and also in Euler $L_3$ configurations, with bounded chaos.

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 secular dynamics of small planetesimals in tight binary systems play a fundamental role in establishing the possibility of accretional collisions in such extreme cases. The most important secular parameters are the forced eccentricity and secular frequency, which depend on the initial conditions of the particles, as well as on the mass and orbital parameters of the secondary star. We construct a second-order theory (with respect to the masses) for the planar secular motion of small planetasimals and deduce new expressions for the forced eccentricity and secular frequency. We also reanalyze the radial velocity data available for Gamma-Cephei and present a series of orbital solutions leading to residuals compatible with the best fits. Finally, we discuss how different orbital configurations for Gamma-Cephei may affect the dynamics of small bodies in circunmstellar motion. For Gamma-Cephei, we find that the classical first-order expressions for the secular frequency and forced eccentricity lead to large inaccuracies around 50 % for semimajor axes larger than one tenth the orbital separation between the stellar components. Low eccentricities and/or masses reduce the importance of the second-order terms. The dynamics of small planetesimals only show a weak dependence with the orbital fits of the stellar components, and the same result is found including the effects of a nonlinear gas drag. Thus, the possibility of planetary formation in this binary system largely appears insensitive to the orbital fits adopted for the stellar components, and any future alterations in the system parameters (due to new observations) should not change this picture. Finally, we show that planetesimals migrating because of gas drag may be trapped in mean-motion resonances with the binary, even though the migration is divergent.

We analyze the dynamics of individual kilometer-size planetesimals in circumstellar orbits of a tight binary system. We include both the gravitational perturbations of the secondary star and a non-linear gas drag stemming from an eccentric gas disk with a finite precession rate. We consider several precession rates and eccentricities for the gas, and compare the results with a static disk in circular orbit. The disk precession introduces three main differences with respect to the classical static case: (i) The equilibrium secular solutions generated by the gas drag are no longer fixed points in the averaged system, but limit cycles with frequency equal to the precession rate of the gas. The amplitude of the cycle is inversely dependent on the body size, reaching negligible values for $\sim 50$ km size planetesimals. (ii) The maximum final eccentricity attainable by small bodies is restricted to the interval between the gas eccentricity and the forced eccentricity, and apsidal alignment is no longer guaranteed for planetesimals strongly coupled with the gas. (iii) The characteristic timescales of orbital decay and secular evolution decrease significantly with increasing precession rates, with values up to two orders of magnitude smaller than for static disks. Finally, we apply this analysis to the $\gamma$-Cephei system and estimate impact velocities for different size bodies and values of the gas eccentricity. For high disk eccentricities, we find that the disk precession decreases the velocity dispersion between different size planetesimals, thus contributing to accretional collisions in the outer parts of the disk. The opposite occurs for almost circular gas disks, where precession generates an increase in the relative velocities.

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.

In the framework of the planar CR3BP for mass parameter mu=0.0121505, corresponding to the Earth-Moon system, we identify and describe 80 families of periodic orbits encircling both the Earth and the Moon ("transfer" orbits). All the orbits in these families have very low energies, most of them corresponding to values of the Jacobi constant C for which the Hill surface is closed at the Lagrangian point L2. All of these orbits have also short period T, generally under six months. Most of the families are composed of orbits that are asymmetric with respect to the Earth-Moon axis. The main results presented for each family are: (i) the characteristic curves T(h), y(h), v_y(h), and v_x(h) on the Poincare section Sigma_1={x=0.836915310,y,v_x>0,v_y} normal to the Earth-Moon axis at the Lagrangian point L1, parameterized by their energy h=-C/2 in the synodic coordinate system; (ii) the stability parameter along each family; (iii) the intersections x_i(h) of the orbits with the Earth-Moon axis, on the Poincare section Sigma_2={x,y=0,v_x},v_y>0}; (iv) plots of some selected orbits and details of their circumlunar region; and (v) numerical data for the intersection of an orbit with Sigma_1 at a reference value of h. Some possible extensions and applications of this work are also discussed.