In about 6 Giga years our Sun will evolve into a red giant and finally end its life as a white dwarf. This stellar metamorphosis will occur to virtually all known host stars of exo-planetary systems and is therefore crucial for their final fate. It is clear that the innermost planets will be engulfed and evaporated during the giant phase and that planets located farther out will survive. However, the destiny of planets in-between, at ~ 1 - 10 au, has not yet been investigated with a multi-planet tidal treatment. We here combine for the first time multi-planet interactions, stellar evolution, and tidal effects in an N-body code to study the evolution of a Neptune-Jupiter planetary system. We report that the fate of the Neptune-mass planet, located closer to the star than the Jupiter-mass planet, can be very different from the fate of a single Neptune. The simultaneous effects of gravitational interactions, mass loss and tides can drive the planetary system towards mean motion resonances. Crossing these resonances affects particularly the eccentricity of the Neptune and thereby also its fate, which can be engulfment, collision with the Jupiter-mass planet, ejection from the system, or survival at a larger separation.

Close-in giant planets represent the most significant evidence of planetary migration. If large exomoons form around migrating giant planets which are more stable (e.g. those in the Solar System), what happens to these moons after migration is still under intense research. This paper explores the scenario where large regular exomoons escape after tidal-interchange of angular momentum with its parent planet, becoming small planets by themselves. We name this hypothetical type of object a \textit{ploonet}. By performing semi-analytical simulations of tidal interactions between a large moon with a close-in giant, and integrating numerically their orbits for several Myr, we found that in $\sim$50 per cent of the cases a young ploonet may survive ejection from the planetary system, or collision with its parent planet and host star, being in principle detectable. Volatile-rich ploonets are dramatically affected by stellar radiation during both planetocentric and siderocentric orbital evolution, and their radius and mass change significantly due to the sublimation of most of their material during time-scales of hundred of Myr. We estimate the photometric signatures that ploonets may produce if they transit the star during the phase of evaporation, and compare them with noisy lightcurves of known objects ("Kronian" stars and non-periodical dips in dusty lightcurves). Additionally, the typical transit timing variations (TTV) induced by the interaction of a ploonet with its planet are computed. We find that present and future photometric surveys' capabilities can detect these effects and distinguish them from those produced by other nearby planetary encounters.

Context. Dynamical studies suggest that most circumbinary discs (CBDs) should be coplanar. However, some theoretical works show that under certain conditions a CBD can become polar, which means that its rotation vector is orthogonal with respect to the binary orbital plane. Interestingly, very recent observations show that polar CBDs exist in nature (e.g. HD 98800). Aims. We test the predictions of CBD alignment around eccentric binaries based on linear theory. In particular, we compare prograde and retrograde CBD configurations. Then, assuming planets form in these systems, we thoroughly characterise the orbital behaviour and stability of misaligned (P-type) particles. This is done for massless and massive particles. Methods. The evolution of the CBD alignment for various configurations was modelled through three-dimensional hydrodynamical simulations. For the orbital characterisation and the analysis stability, we relied on long-term N-body integrations and structure and chaos indicators, such as $\Delta e$ and MEGNO. Results. We confirm previous analytical predictions on CBD alignment, but find an unexpected symmetry breaking between prograde and retrograde configurations. Therefore, the likelihood of becoming polar for a highly misaligned CBD is higher than previously thought. Regarding the stability of circumbinary P-type planets (also know as Tatooines), polar orbits are stable over a wide range of binary parameters. In particular, for binary eccentricities below 0.4 the orbits are stable for any value of the binary mass ratio. In the absence of gas, planets with masses below $10^{-5}\,M_{\odot}$ have negligible effects on the binary orbit. Finally, we suggest that mildly eccentric equal-mass binaries should be searched for polar Tatooines.

Dynamical studies suggest that most of the circumbinary discs (CBDs) should be coplanar. However, under certain initial conditions, the CBD can evolve toward polar orientation. Here we extend the parametric study of polar configurations around detached close-in binaries through $N$-body simulations. For polar configurations around binaries with mass ratios $q$ below $0.7$, the nominal location of the mean motion resonance (MMR) $1~:~4$ predicts the limit of stability for $e_{B} > 0.1$. Alternatively, for $e_{B} < 0.1$ or $q \sim 1$, the nominal location of the MMR $1~:~3$ is the closest stable region. The presence of a} giant planet increases the region of forbidden polar configurations around low mass ratio binaries with eccentricities $e_B\sim0.4$ with respect to rocky earth-like planets. For equal mass stars, the eccentricity excitation $\Delta e$ of polar orbits smoothly increases with decreasing distance to the binary. For $q<1$, $\Delta e$ can reach values as high as $0.4$. Finally, we studied polar configurations around $HD~98800BaBb$ and show that the region of stability is strongly affected by the relative positions of the nodes. The most stable configurations in the system correspond to polar particles, which are not expected to survive on longer time-scales due to the presence of the external perturber HD~$98800AaAb$.

Dwarf stars with debris discs and planets appear to be excellent laboratories to study the core accretion theory of planets formation. These systems are however, insufficiently studied. In this paper we present the main metallicity and lithium abundance properties of these stars together with stars with only debris discs and stars with only planets. Stars without detected planets nor discs are also considered. The analysed sample is formed by main-sequence FGK field single stars. Apart from the basic stellar parameters, we include the use of dusty discs masses. The main results show for the first time that the dust mass of debris disc stars with planets correlate with metallicity. We confirm that these disc dust masses are related to their central stellar masses. Separately, the masses of stars and those of planets also correlate with metallicity. We conclude that two conditions are necessary to form giant planets: to have a sufficient metallicity and also a sufficient protoplanetary mass of gas and dust. The debris discs masses of stars without giant planets do not correlate with metallicity, because they do not fulfil these two conditions. Concerning lithium, by adopting a stellar model for lithium depletion based on a strong interaction between the star and a protoplanetary disc, we found that in agreement with the model predictions, observations indicate that the main lithium depletion occurs during this initial protoplanetary evolution stage. We show that the ultimately lithium depletion is independent of the presence or absence of planets and appears to be only age dependent.

The stability of satellites in the solar system is affected by the so-called evection resonance. The moons of Saturn, in particular, exhibit a complex dynamical architecture in which co-orbital configurations occur, especially close to the planet where this resonance is present. We address the dynamics of the evection resonance, with particular focus on the Saturn system, and compare the known behavior of the resonance for a single moon to that of a pair of moons in co-orbital trojan configuration. We developed an analytic expansion of the averaged Hamiltonian of a trojan pair of bodies, including the perturbation from a distant massive body. {The analysis of the corresponding equilibrium points was restricted to the asymmetric apsidal corotation solution of the co-orbital dynamics.} We also performed numerical N-body simulations to construct dynamical maps of the stability of the evection resonance in the Saturn system, and to study the effects of this resonance under the migration of trojan moons caused by tidal dissipation. The structure of the phase space of the evection resonance for trojan satellites is similar to that of a single satellite, differing in that the libration centers are displaced from their standard positions by an angle that depends on the periastron difference $\varpi_2-\varpi_1$ and on the mass ratio $m_2/m_1$ of the trojan pair. The interaction with the inner evection resonance may have been relevant during the early evolution of the Saturn moons Tethys, Dione, and Rhea. In particular, Rhea may have had trojan companions in the past that were lost when it crossed the evection resonance, while Tethys and Dione may either have retained their trojans or have never crossed the evection. This may help to constrain the dynamical processes that led to the migration of these satellites and to the evection itself.

The stability of planets in the alpha-Centauri AB stellar system has been studied extensively. However, most studies either focus on the orbital plane of the binary or consider inclined circular orbits. Here, we numerically investigate the stability of a possible planet in the alpha-Centauri AB binary system for S-type orbits in an arbitrary spatial configuration. In particular, we focus on inclined orbits and explore the stability for different eccentricities and orientation angles. We show that large stable and regular regions are present for very eccentric and inclined orbits, corresponding to libration in the Lidov-Kozai resonance. We additionally show that these extreme orbits can survive over the age of the system, despite the effect of tides. Our results remain qualitatively the same for any compact binary system.

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 test a crossing orbit stability criterion for eccentric planetary systems, based on Wisdom's criterion of first order mean motion resonance overlap (Wisdom, 1980). We show that this criterion fits the stability regions in real exoplanet systems quite well. In addition, we show that elliptical orbits can remain stable even for regions where the apocenter distance of the inner orbit is larger than the pericenter distance of the outer orbit, as long as the initial orbits are aligned. The analytical expressions provided here can be used to put rapid constraints on the stability zones of multi-planetary systems. As a byproduct of this research, we further show that the amplitude variations of the eccentricity can be used as a fast-computing stability indicator.

In this paper, we study the behavior of a pair of co-orbital planets, both orbiting a central star on the same plane and undergoing tidal interactions. Our goal is to investigate final orbital configurations of the planets, initially involved in the 1/1 mean-motion resonance (MMR), after long-lasting tidal evolution. The study is done in the form of purely numerical simulations of the exact equations of motions accounting for gravitational and tidal forces. The results obtained show that, at least for equal mass planets, the combined effects of the resonant and tidal interactions provoke the orbital instability of the system, often resulting in collision between the planets. We first discuss the case of two hot-super-Earth planets, whose orbital dynamics can be easily understood in the frame of our semi-analytical model of the 1/1 MMR. Systems consisting of two hot-Saturn planets are also briefly discussed.

The existence of multiple planetary systems involved in mean motion conmensurabilities has increased significantly since the Kepler mission. Although most correspond to 2-planet resonances, multiple resonances have also been found. The Laplace resonance is a particular case of a three-body resonance where the period ratio between consecutive pairs is n_1/n_2 near to n_2/n_3 near to 2/1. It is not clear how this triple resonance can act in order to stabilize (or not) the systems. The most reliable extrasolar system located in a Laplace resonance is GJ876 because it has two independent confirmations. However best-fit parameters were obtained without previous knowledge of resonance structure and no exploration of all the possible stable solutions for the system where done. In the present work we explored the different configurations allowed by the Laplace resonance in the GJ876 system by varying the planetary parameters of the third outer planet. We find that in this case the Laplace resonance is a stabilization mechanism in itself, defined by a tiny island of regular motion surrounded by (unstable) highly chaotic orbits. Low eccentric orbits and mutual inclinations from -20 to 20 degrees are compatible with the observations. A definite range of mass ratio must be assumed to maintain orbital stability. Finally we give constrains for argument of pericenters and mean anomalies in order to assure stability for this kind of systems.

We investigate the stability of prograde versus retrograde planets in circular binary systems using numerical simulations. We show that retrograde planets are stable up to distances closer to the perturber than prograde planets. We develop an analytical model to compute the prograde and retrograde mean motion resonances' locations and separatrices. We show that instability is due to single resonance forcing, or caused by nearby resonances' overlap. We validate our results regarding the role of single resonances and resonances' overlap on orbit stability, by computing surfaces of section of the CR3BP. We conclude that the observed enhanced stability of retrograde planets with respect to prograde planets is due to essential differences between the phase-space topology of retrograde versus prograde resonances (at p/q mean motion ratio, prograde resonance is of order p - q while retrograde resonance is of order p + q).

The HD\,196885 system is composed of a binary star and a planet orbiting the primary. The orbit of the binary is fully constrained by astrometry, but for the planet the inclination with respect to the plane of the sky and the longitude of the node are unknown. Here we perform a full analysis of the HD\,196885 system by exploring the two free parameters of the planet and choosing different sets of angular variables. We find that the most likely configurations for the planet is either nearly coplanar orbits (prograde and retrograde), or highly inclined orbits near the Lidov-Kozai equilibrium points, i = 44^{\circ} or i = 137^{\circ} . Among coplanar orbits, the retrograde ones appear to be less chaotic, while for the orbits near the Lidov-Kozai equilibria, those around \omega= 270^{\circ} are more reliable, where \omega_k is the argument of pericenter of the planet's orbit with respect to the binary's orbit. From the observer's point of view (plane of the sky) stable areas are restricted to (I1, \Omega_1) \sim (65^{\circ}, 80^{\circ}), (65^{\circ},260^{\circ}), (115^{\circ},80^{\circ}), and (115^{\circ},260^{\circ}), where I1 is the inclination of the planet and \Omega_1 is the longitude of ascending node.

We analyze the possibilities of detection of hypothetical exoplanets in coorbital motion from synthetic radial velocity (RV) signals, taking into account different types of stable planar configurations, orbital eccentricities and mass ratios. For each nominal solution corresponding to small-amplitude oscillations around the periodic solution, we generate a series of synthetic RV curves mimicking the stellar motion around the barycenter of the system. We then fit the data sets obtained assuming three possible different orbital architectures: (a) two planets in coorbital motion, (b) two planets in a 2/1 mean-motion resonance, and (c) a single planet. We compare the resulting residuals and the estimated orbital parameters. For synthetic data sets covering only a few orbital periods, we find that the discrete radial velocity signal generated by a coorbital configuration could be easily confused with other configurations/systems, and in many cases the best orbital fit corresponds to either a single planet or two bodies in a 2/1 resonance. However, most of the incorrect identifications are associated to dynamically unstable solutions. We also compare the orbital parameters obtained with two different fitting strategies: a simultaneous fit of two planets and a nested multi-Keplerian model. We find that the nested models can yield incorrect orbital configurations (sometimes close to fictitious mean-motion resonances) that are nevertheless dynamically stable and with orbital eccentricities lower than the correct nominal solutions. Finally, we discuss plausible mechanisms for the formation of coorbital configurations, by the interaction between two giant planets and an inner cavity in the gas disk. For equal mass planets, both Lagrangian and anti-Lagrangian configurations can be obtained from same initial condition depending on final time of integration.

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 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.