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 analyze the distribution of extrasolar planets (both confirmed and Kepler candidates) according to their orbital periods P and planetary radii R. Among confirmed planets, we find compelling evidence for a paucity of bodies with 3 < R < 10 R_\oplus, where R_\oplus in the Earth's radius, and P < 2-3 days. We have christened this region a "sub-Jovian Pampas". The same trend is detected in multiplanet Kepler candidates. Although approximately 16 Kepler single-planet candidates inhabit this Pampas, at least 7 are probable false positives (FP). This last number could be significantly higher if the ratio of FP is higher than 10%, as suggested by recent studies. In a second part of the paper we analyze the distribution of planets in the (P,R) plane according to stellar metallicities. We find two interesting trends: (i) a lack of small planets (R < 4 R_\oplus) with orbital periods P < 5 days in metal-poor stars, and (ii) a paucity of sub-Jovian planets (4 R_\oplus < R < 8 R_\oplus) with P < 100 days, also around metal-poor stars. Although all these trends are preliminary, they appear statistically significant and deserve further scrutiny. If confirmed, they could represent important constraints on theories of planetary formation and dynamical evolution.

In the last years several exoplanets have been discovered that orbit one component of a compact binary system (separation < 50 astronomical units), the probably best-known case is gamma-Cephei. So far, all attempts to explain the in-situ formation of these planets has been unsuccessful, in part because of the strong gravitational perturbations of the secondary star on any initial planetesimal swarm. Here we test whether planetary bodies in compact binaries, in particular gamma-Cephei, could have originated from a close encounter with a passing star, assuming initial configurations for the stellar system suitable for planetary formation. In other words, we analyze whether the orbital configuration of the current binary system might have been generated after the formation of the planet, and as a consequence of a close encounter with a third star in hyperbolic orbit. We performed a series of time-reverse N-body simulations of stellar scattering events in which the present-day configuration of gamma-Cephei was used as the initial condition plus a hypothetical third star as an impactor. We analyzed which configurations and system parameters could have given birth to the current system. Depending on the maximum impact velocity allowed for accretional collisions, we find that between 1% and 5% of stellar encounters correspond to an "original" system in which planetary formation around the primary star is not inhibited by the secondary, but is acceptable within the classical core-accretion scenario. Thus, although not highly probable, it is plausible that stellar encounters may have played a significant role in shaping these types of exoplanetary systems.

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.

We present the dynamical structure of the phase space of the planar planetary 2/1 mean-motion resonance (MMR). Inside the resonant domain, there exist two families of periodic orbits, one associated to the librational motion of the critical angle ($\sigma$-family) and the other related to the circulatory motion of the angle between the pericentres ($\Delta\varpi$-family). The well-known apsidal corotation resonances (ACR) appear at the intersections of these families. A complex web of secondary resonances exists also for low eccentricities, whose strengths and positions are dependent on the individual masses and spatial scale of the system. Depending on initial conditions, a resonant system is found in one of the two topologically different states, referred to as \textit{internal} and \textit{external} resonances. The internal resonance is characterized by symmetric ACR and its resonant angle is $2\,\lambda_2-\lambda_1-\varpi_1$, where $\lambda_i$ and $\varpi_i$ stand for the planetary mean longitudes and longitudes of pericentre, respectively. In contrast, the external resonance is characterized by asymmetric ACR and the resonant angle is $2\,\lambda_2-\lambda_1-\varpi_2$. We show that systems with more massive outer planets always envolve inside internal resonances. The limit case is the well-known asteroidal resonances with Jupiter. At variance, systems with more massive inner planets may evolve in either internal or external resonances; the internal resonances are typical for low-to-moderate eccentricity configurations, whereas the external ones for high eccentricity configurations of the systems. In the limit case, analogous to Kuiper belt objects in resonances with Neptune, the systems are always in the external resonances characterized by asymmetric equilibria.

Exoplanets show a pile-up of Jupiter-size planets in orbits with a 3-day period. A fraction of these hot Jupiters have retrograde orbits with respect to the parent star's rotation. To explain these observations we performed a series of numerical integrations of planet scattering followed by the tidal circularization. We considered planetary systems having 3 and 4 planets initially. We found that the standard Kozai migration is an inefficient mechanism for the formation of hot Jupiters. Our results show the formation of two distinct populations of hot Jupiters. The inner population of hot Jupiters with semimajor axis a < 0.03 AU formed in the systems where no planetary ejections occurred. This group contained a significant fraction of highly inclined and retrograde orbits, with distributions largely independent of the initial setup. However, our follow-up integrations showed that this populations was transient with most planets falling inside the Roche radius of the star in <1 Gyr. The outer population of hot Jupiters formed in systems where at least one planet was ejected. This population survived the effects of tides over >1 Gyr. The semimajor axis distribution of Population II fits nicely the observed 3-day pile-up. The inclination distribution of the outer hot planets depends on the number of planets in the initial systems and the 4-planet case showed a larger proportion (up to 10%), and a wider spread in inclination values. As the later results roughly agrees with observations, this may suggest that the planetary systems with observed hot Jupiters were originally rich in the number of planets, some of which were ejected. In a broad perspective, our work therefore hints on an unexpected link between the hot Jupiters and recently discovered free floating planets.

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.

This investigation uses the excellent HARPS radial velocity measurements of CoRoT-7 to re-determine the planet masses and to explore techniques able to determine mass and elements of planets discovered around active stars when the relative variation of the radial velocity due to the star activity cannot be considered as just noise and can exceed the variation due to the planets. The main technique used here is a self-consistent version of the high-pass filter used by Queloz et al. (2009) in the first mass determination of CoRoT-7b and CoRoT-7c. The results are compared to those given by two alternative techniques: (1) The approach proposed by Hatzes et al. (2010) using only those nights in which 2 or 3 observations were done; (2) A pure Fourier analysis. In all cases, the eccentricities are taken equal to zero as indicated by the study of the tidal evolution of the system; the periods are also kept fixed at the values given by Queloz et al. Only the observations done in the time interval BJD 2,454,847 - 873 are used because they include many nights with multiple observations; otherwise it is not possible to separate the effects of the rotation fourth harmonic (5.91d = Prot/4) from the alias of the orbital period of CoRoT-7b (0.853585 d). The results of the various approaches are combined to give for the planet masses the values 8.0 \pm 1.2 MEarth for CoRoT-7b and 13.6 \pm 1.4 MEarth for CoRoT 7c. An estimation of the variation of the radial velocity of the star due to its activity is also given.The results obtained with 3 different approaches agree to give masses larger than those in previous determinations. From the existing internal structure models they indicate that CoRoT-7b is a much denser super-Earth. The bulk density is 11 \pm 3.5 g.cm-3 . CoRoT-7b may be rocky with a large iron core.

The lower limit to the distribution of orbital periods P for the current population of close-in exoplanets shows a distinctive discontinuity located at approximately one Jovian mass. Most smaller planets have orbital periods longer than P~2.5 days, while higher masses are found down to P~1 day. We analyze whether this observed mass-period distribution could be explained in terms of the combined effects of stellar tides and the interactions of planets with an inner cavity in the gaseous disk. We performed a series of hydrodynamical simulations of the evolution of single-planet systems in a gaseous disk with an inner cavity mimicking the inner boundary of the disk. The subsequent tidal evolution is analyzed assuming that orbital eccentricities are small and stellar tides are dominant. We find that most of the close-in exoplanet population is consistent with an inner edge of the protoplanetary disk being located at approximately P>2 days for solar-type stars, in addition to orbital decay having been caused by stellar tides with a specific tidal parameter on the order of Q'*=10^7. The data is broadly consistent with planets more massive than one Jupiter mass undergoing type II migration, crossing the gap, and finally halting at the interior 2/1 mean-motion resonance with the disk edge. Smaller planets do not open a gap in the disk and remain trapped in the cavity edge. CoRoT-7b appears detached from the remaining exoplanet population, apparently requiring additional evolutionary effects to explain its current mass and semimajor axis.

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.

The Transit Timing Variation (TTV) method relies on monitoring changes in timing of transits of known exoplanets. Non-transiting planets in the system can be inferred from TTVs by their gravitational interaction with the transiting planet. The TTV method is sensitive to low-mass planets that cannot be detected by other means. Here we describe a fast algorithm that can be used to determine the mass and orbit of the non-transiting planets from the TTV data. We apply our code, ttvim.f, to a wide variety of planetary systems to test the uniqueness of the TTV inversion problem and its dependence on the precision of TTV observations. We find that planetary parameters, including the mass and mutual orbital inclination of planets, can be determined from the TTV datasets that should become available in near future. Unlike the radial velocity technique, the TTV method can therefore be used to characterize the inclination distribution of multi-planet systems.

The three-dimensional secular behavior of a system composed of a central star and two massive planets is modeled semi-analytically in the frame of the general three-body problem. The main dynamical features of the system are presented in geometrical pictures allowing us to investigate a large domain of the phase space of this problem without time-expensive numerical integrations of the equations of motion and without any restriction on the magnitude of the planetary eccentricities, inclinations and mutual distance. Several regimes of motion of the system are observed. With respect to the secular angle $\Delta\varpi$, possible motions are circulations, oscillations (around 0 and $180^\circ$), and high eccentricity/inclination librations in secular resonances. With respect to the arguments of pericenter, $\omega_1$ and $\omega_2$, possible motions are direct circulation and high-inclination libration around $\pm 90^\circ$ in the Lidov-Kozai resonance. The regions of transition between domains of different regimes of motion are characterized by chaotic behavior. We apply the analysis to the case of the two outer planets of the $\upsilon$ Andromedae system, observed edge-on. The topology of the 3--D phase space of this system is investigated in detail by means of surfaces of section, periodic orbits and dynamical spectra, mapping techniques and numerical simulations. We obtain the general structure of the phase space, and the boundaries of the spatial secular stability. We find that this system is secularly stable in a large domain of eccentricities and inclinations.

This paper is a review of the dynamics of a system of planets. It includes the study of averaged equations in both non-resonant and resonant systems and shows the great deal of situations in which the angle between the two semi-major axes oscillates around a constant value. It introduces the Hamiltonian equations of the $N$-planet problem and Poincar\'e's reduction of them to $3N$ degrees of freedom with a detailed discussion of the non-osculating ``canonical'' heliocentric Keplerian elements that should be used with Poincar\'e relative canonical variables. It also includes Beaug\'e's approximation to expand the disturbing function in the exoplanetary case where masses and eccentricities are large. The paper is concluded with a discussion of systems captured into resonance and their evolution to symmetric and asymmetric stationary solutions with apsidal corotation.

We analyze the possible relationship between the current orbital elements fits of known exoplanets in the 2/1 mean-motion resonance and the expected orbital configuration due to migration. It is found that, as long as the orbital decay was sufficiently slow to be approximated by an adiabatic process, all captured planets should be in apsidal corotations. In other words, they should show a simultaneous libration of both the resonant angle and the difference in longitudes of pericenter. We present a complete set of corotational solutions for the 2/1 commensurability, including previously known solutions and new results. Comparisons with observed exoplanets show that current orbital fits of three known planetary systems in this resonance are either consistent with apsidal corotations (GJ876 and HD82943) or correspond to bodies with uncertain orbits (HD160691). Finally, we discuss the applicability of these results as a test for the planetary migration hypothesis itself. If all future systems in this commensurability are found to be consistent with corotational solutions, then resonance capture of these bodies through planetary migration is a working hypothesis. Conversely, If any planetary pair is found in a different configuration, then either migration did not occur for those bodies, or it took a different form than currently believed.

In this paper we present numerical simulations of the evolution of planets or massive satellites captured in the 2/1 and 3/1 resonances, under the action of an anti-dissipative tidal force. The evolution of resonant trapped bodies show a richness of solutions: librations around stationary symmetric solutions with aligned periapses (w1-w2=0) or anti-aligned periapses (w1-w2=180 deg), and librations around stationary solutions in which the periapses configuration is fixed, but with w1-w2 taking values in a wide range of angles. Many of these solutions exist for large values of the eccentricities and, during the semimajor axis drift, the solutions show turnabouts from one configuration to another. These results are valid for other non-conservative forces leading to adiabatic covergent migration and capture into one of these resonances.

In recent years several pairs of extrasolar planets have been discovered in the vicinity of mean-motion commensurabilities. In some cases, such as the Gliese 876 system, the planets seem to be trapped in a stationary solution, the system exhibiting a simultaneous libration of the resonant angle theta_1 = 2 lambda_2 - lambda_1 - varpi_1 and of the relative position of the pericenters. In this paper we analyze the existence and location of these stable solutions, for the 2/1 and 3/1 resonances, as function of the masses and orbital elements of both planets. This is undertaken via an analytical model for the resonant Hamiltonian function. The results are compared with those of numerical simulations of the exact equations. In the 2/1 commensurability, we show the existence of three principal families of stationary solutions: (i) aligned orbits, in which theta_1 and varpi_1 - varpi_2 both librate around zero, (ii) anti-aligned orbits, in which theta_1=0 and the difference in pericenter is 180 degrees, and (iii) asymmetric stationary solutions, where both the resonant angle and varpi_1 - varpi_2 are constants with values different of 0 or 180 degrees. Each family exists in a different domain of values of the mass ratio and eccentricities of both planets. Similar results are also found in the 3/1 resonance. We discuss the application of these results to the extrasolar planetary systems and develop a chart of possible planetary orbits with apsidal corotation. We estimate, also, the maximum planetary masses in order that the stationary solutions are dynamically stable.