The present paper is the culminating one of a series aimed to contribute to the understanding of the kinematic structures of the solar neighbourhood (SN), explaining the origin of the Local Arm and relating the moving groups with the spiral-arms resonances in the disk. With a model for the Galactic potential, with the Sun inside the spiral corotation resonance (CR), we integrate the 2D orbits of test particles distributed in birthplaces along the main spiral arms, the Local Arm, and in the axisymmetric disk. A comparison of the resulting U-V plane of the SN with that provided by Gaia DR2 confirms our previous conclusion that the moving groups of Coma Berenices, Pleiades, and Hyades are associated with the CR, and that the Hercules stream is formed by the bulk of high-order inner Lindblad resonances. The kinematic structures result from stellar orbits trapped by the spiral resonances in a timespan of ~ 1 Gyr, indicating the long-living nature of the spiral structure and challenging recent arguments in favor of short-lived structures originated from incomplete phase mixing in the Galactic disk. As a byproduct, our simulations give some insight into the birthplaces of the stars presently located in the SN; the majority of stars of the main moving groups and the Hercules stream were likely born in the Local Arm, while stars of the Sirius group possibly originated from the outer segment of the Sagittarius-Carina arm. We also propose the spiral resonances as the dynamical origin for the diagonal ridges in the Galactic distribution of rotation velocities.

In this work, we adapt a module for planetary formation within the hydrodynamic code FARGO3D. Planetary formation is modeled by a solid core accretion scenario, with the core growing in oligarchic regime. The initial superficial density of planetesimals is proportional to the initial superficial density of gas in the disc. We include a numerical approach to describe the evolution of the eccentricity and the inclination of planetesimals during the formation. This approach impacts directly on the accretion rate of solids. When the core reaches a critical mass, gas accretion begins, following the original FARGO scheme adapted to the FARGO3D code. To exemplify how the module for planetary formation can be used, we investigate the migration of a planet in a two-dimensional, locally isothermal gas disc with a prescribed accretion rate, analyzing the timescale involved in the planetary migration process along with the timescale for planetary formation. The analysis reveals that the mass of the nucleus must be close to its critical value when crossing the ice line to avoid the planet's fall into the stellar envelope. This will allow enough time for the planet to initiate runaway gas accretion, leading to a rapid mass increase and entering type II planetary migration.

FARGO3D has been proposed to investigate numerically the gravitational interaction of a planet with the gas disc, providing original outcomes. However, FARGO3D does not consider the gas accretion of the planet which may affect the planetary migration process. Thus, the aim of this paper is to include a model for the gas accretion in the FARGO3D code. For this task, we choose Russell's scheme, which is an adaptation of Kley's model for the gas accretion on to migrating planets. Initially, we study the influence of the gas accretion on to the planet on type II migration of giant planets. For this purpose, we follow the evolution of planets in a 2D locally isothermal disc with a specific accretion rate and different values of viscosity and planetary mass considering two scenarios. In the first, the gas mass is withdrawn from the disc and is not added to the planet. In the second, the planets migrate, while their masses grow due to the gas accretion. Then, we extend our study to estimate the time needed for a low-mass planet to open a gap in the gas disc and compare its value with the characteristic time of type I migration.

We analyze the orbital motion of two natural satellites initially in co-orbital configuration with a third (guiding) satellite embedded into a circumplanetary gas disc and undergoing tidal interactions with the central planet. By solving the exact equations of motion, including the dissipative effects and the mutual gravitational perturbations, we investigate the configuration of the system soon after the guiding satellite is disrupted when crossing the Roche limit with the central planet during its orbital decay. The application is done for the system composed of Saturn and the two small co-orbital satellites Janus and Epimetheus. We perform a large number of numerical simulations, varying the mass of the hypothetical guiding satellite, Saturn's tidal quality factor and the initial configuration and the masses of the small satellites. Analyzing the orbital configurations of Janus and Epimetheus soon after the disruption of the guiding satellite, we obtain that the mutual horseshoe co-orbital orbits appear as a natural outcome of the numerical simulations.

Using the Gaia data release 2 (DR2), we analyzed the distribution of stars in the close vicinity of the Sun in the full 3D position-velocity space. We have found no evidence of incomplete phase mixing in the vertical direction of the disk, which could be originated by some external events. We show that the vertical phase space spiral $Z$-$V_z$ is produced by the well-known moving groups (MGs), mainly by Coma-Berenices, Pleiades-Hyades and Sirius, when the statistical characteristics (mean, median, or mode) of the azimuthal velocity $V_\varphi$ are used to analyze the distribution in the vertical position-velocity plane. This result does not invoke external perturbations and is independent on the internal dynamical mechanisms that originate the MGs. Our conclusions counterbalance current arguments in favor of short-lived (between 300 and 900 Myr) structures in the solar neighborhood. Contrarily, they support the hypothesis of a longer formation time scale (around a few Gyr) for the MGs.

The aim of this work is to contribute to the understanding of the stellar velocity distribution in the solar neighborhood (SN). We propose that the structures on the $U-V$ planes, known as the moving groups, can be mainly explained by the spiral arms perturbations. The applied model of the Galactic disk and spiral arms, with the parameters defined by observational data and with pattern speed $\Omega_p=$28.0 km s$^{-1}$ kpc$^{-1}$, is the same that allowed us to explain the origin of the Local Arm and the Sun's orbit trapped inside the corotation resonance (CR). We show that the $U-V$ picture of the SN consists of the main component, associated with the CR, and the inner and outer structures, which we could associate with the Hercules and Sirius streams, respectively. The Coma-Berenices and Hyades-Pleiades groups and the Sun itself belong to the main part. The substructures of Hercules are formed mainly by the nearby 8/1, 12/1, and even 6/1 inner Lindblad resonances, while Sirius is shaped by the bulk of overlapping outer Lindblad resonances, -8/1, -12/1, -16/1, which are stuck to the CR. This richness in resonances only exists near corotation, which should be of the spiral arms, not of the Galactic bar, whose stable corotation zone is far away from the Sun. The model's predictions of the velocity distribution match qualitatively and quantitatively the distribution provided by Gaia DR2.

Observational data indicate that the Milky Way is a barred spiral galaxy. Computation facilities and availability of data from Galactic surveys stimulate the appearance of models of the Galactic structure. More efforts to build dynamical models containing both spiral arms and the central bar/bulge are needed. We expand the study of the stellar dynamics in the Galaxy by adding the bar/bulge component to a model with spiral arms introduced in our previous paper. The model is tested by applying it to the solar neighborhood, where observational data are more precise. We model analytically the potential of the Galaxy to derive the force field in its equatorial plane. The model comprises an axisymmetric disc derived from the observable rotation curve, four arms with Gaussian-shaped groove profiles, and a classical elongated/oblate ellipsoidal bar/bulge structure. The parameters describing the bar/bulge are constrained by observations and the stellar dynamics, and their possible limits are determined. A basic model results in a bar of 2.9 kpc in length, with a mass of the order of a few 10$^9M_\odot$. The size and orientation of the bar are also restricted by the position of masers with VLBI distances. The bar's rotation speed is constrained to $\Omega_{\rm bar}<50$ km s$^{-1}$ kpc$^{-1}$ taking into account the allowed mass range. We conclude that our basic model is compatible with observations and with the dynamical constraints. The model explains simultaneously the bulk of the main moving groups, associated here with the spiral corotation resonance, and the Hercules stream, associated with several inner high-order spiral resonances; in particular, with the 8/1 resonance. From the dynamical constraints on the bar's angular speed, it is unlikely that the bar's OLR lies near the solar circle; moreover, its proximity would compromise the stability of the Local Arm structure.

The Local arm of the Milky Way, a short spiral feature near the Sun whose existence is known for decades, was recently observed in detail with different tracers. Many efforts have been dedicated to elaborate plausible hypotheses concerning the origin of the main spiral arms of the Galaxy; however, up to now, no specific mechanism for the origin of the Local arm was proposed. Here we explain, for the first time, the Local arm as an outcome of the spiral corotation resonance, which traps arm tracers and the Sun inside it. We show that the majority of maser sources belonging to the Local arm, together with the Sun, evolve inside the corotation resonance, never crossing the main spiral arms but instead oscillating in the region between them. This peculiar behavior of the Sun could have numerous consequences to our understanding of the local kinematics of stars, the Galactic Habitable Zone, and the Solar System evolution.

Context: Resonances in the stellar orbital motion under perturbations from spiral arms structure play an important role in the evolution of the disks of spiral galaxies. The epicyclic approximation allows the determination of the corresponding resonant radii on the equatorial plane (for nearly circular orbits), but is not suitable in general. Aims: We expand the study of resonant orbits by analysing stellar motions perturbed by spiral arms with Gaussian-shaped profiles without any restriction on the stellar orbital configurations, and we expand the concept of Lindblad (epicyclic) resonances for orbits with large radial excursions. Methods: We define a representative plane of initial conditions, which covers the whole phase space of the system. Dynamical maps on representative planes are constructed numerically, in order to characterize the phase-space structure and identify the precise location of resonances. The study is complemented by the construction of dynamical power spectra, which provide the identification of fundamental oscillatory patterns in the stellar motion. Results: Our approach allows a precise description of the resonance chains in the whole phase space, giving a broader view of the dynamics of the system when compared to the classical epicyclic approach, even for objects in retrograde motion. The analysis of the solar neighbourhood shows that, depending on the current azimuthal phase of the Sun with respect to the spiral arms, a star with solar kinematic parameters may evolve either inside the stable co-rotation resonance or in a chaotic zone. Conclusions: Our approach contributes to quantifying the domains of resonant orbits and the degree of chaos in the whole Galactic phase-space structure. It may serve as a starting point to apply these techniques to the investigation of clumps in the distribution of stars in the Galaxy, such as kinematic moving groups.

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.

In this paper, we use a semi-analytical approach to analyze the global structure of the phase space of the planar planetary 3/1 mean-motion resonance, in cases where the outer planet is more massive than its inner companion. We show that the resonant dynamics can be described using only two fundamental parameters, the total angular momentum and the scaling parameter. The topology of the Hamiltonian function describing the resonant behaviour is studied on the representative planes that allows us to investigate a large domain of the phase space of the three-body problem without time-expensive numerical integrations of the equations of motion, and without any restriction on the magnitude of the planetary eccentricities. The families of the well known Apsidal Corotation Resonances (ACR) parameterized by the planetary mass ratio are obtained and their stability is analyzed. The main dynamical features in the domains around ACR are also investigated in detail by means of spectral analysis techniques, which allow us to detect the regions of different regimes of motion of resonant systems. The construction of dynamical maps for various values of the total angular momentum shows the evolution of domains of stable motion with the eccentricities, identifying possible configurations suitable for exoplanetary systems.

It has been suggested that HD 200964 is the first exoplanetary system with two Jovian planets evolving in the 4/3 mean- motion resonance. Previous scenarios to simulate the formation of two giant planets in the stable 4/3 resonance configuration have failed. Moreover, the orbital parameters available in the literature point out an unstable configuration of the planetary pair. The purpose of this paper is i) to determine the orbits of the planets from the RV measurements and update the value of the stellar mass (1.57 M), ii) to analyse the stability of the planetary evolution in the vicinity and inside the 4/3 MMR, and iii) to elaborate a possible scenario for the formation of systems in the 4/3 MMR. The results of the formation simulations are able to very closely reproduce the 4/3 resonant dynamics of the best-fit config- uration obtained in this paper. Moreover, the confidence interval of the fit matches well with the very narrow stable region of the 4/3 mean-motion resonance. The formation process of the HD 200964 system is very sensitive to the planetary masses and protoplanetary disk parameters. Only a thin, flat disk allows the embryo-sized planets to reach the 4/3 resonant configuration. The stable evolution of the resonant planets is also sensitive to the mass of the central star, because of overlapping high-order resonances inside the 4/3 resonance. Regardless of the very narrow domain of stable motion, the confidence interval of our fit closely matches the stability area.

In this paper, we investigate whether hypothetical Earth-like planets have high probability of remaining on stable orbits inside the habitable zones around the stars A and B of {\alpha} Centauri, for lengths of time compatible with the evolution of life. We introduce a stability criterion based on the solution of the restricted three-body problem and apply it to the {\alpha} Centauri system. In this way, we determine the regions of the short-term stability of the satellite-type (S-type) planetary orbits, in both planar and three-dimensional cases. We also study the long-term stability of hypothetical planets through the dynamical mapping of the habitable zones of the stars. The topology of the maps is analyzed using the semi-analytical secular Hamiltonian model and possible processes responsible for long-lasting instabilities are identified. We verify that the planetary motion inside the habitable zones is regular, regardless of high eccentricities, for inclinations smaller than 40{\deg}. We show that the variation of the orbital distance of the planet located in the habitable zones of the binary is comparable to that of Earth, if the planet is close to the Mode I stationary solution. This result brings positive expectations for finding habitable planets in binary stars.

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.

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.

We investigate the spin behavior of close-in rocky planets and the implications for their orbital evolution. Considering that the planet rotation evolves under simultaneous actions of the torque due to the equatorial deformation and the tidal torque, both raised by the central star, we analyze the possibility of temporary captures in spin-orbit resonances. The results of the numerical simulations of the exact equations of motions indicate that, whenever the planet rotation is trapped in a resonant motion, the orbital decay and the eccentricity damping are faster than the ones in which the rotation follows the so-called pseudo-synchronization. Analytical results obtained through the averaged equations of the spin-orbit problem show a good agreement with the numerical simulations. We apply the analysis to the cases of the recently discovered hot super-Earths Kepler-10 b, GJ 3634 b and 55 Cnc e. The simulated dynamical history of these systems indicates the possibility of capture in several spin-orbit resonances; particularly, GJ 3634 b and 55 Cnc e can currently evolve under a non-synchronous resonant motion for suitable values of the parameters. Moreover, 55 Cnc e may avoid a chaotic rotation behavior by evolving towards synchronization through successive temporary resonant trappings.

We investigate the secular dynamics of two-planet coplanar systems evolving under mutual gravitational interactions and dissipative forces. We consider two mechanisms responsible for the planetary migration: star-planet (or planet-satellite) tidal interactions and interactions of a planet with a gaseous disc. We show that each migration mechanism is characterized by a specific law of orbital angular momentum exchange. Calculating stationary solutions of the conservative secular problem and taking into account the orbital angular momentum leakage, we trace the evolutionary routes followed by the planet pairs during the migration process. This procedure allows us to recover the dynamical history of two-planet systems and constrain parameters of the involved physical processes.

We have done a new analysis of the available observations for the GJ581 exoplanetary system. Today this system is controversial due to choices that can be done in the orbital determination. The main ones are the ocurrence of aliases and the additional bodies - the planets f and g - announced in Vogt et al. 2010. Any dynamical study of exoplanets requires the good knowledge of the orbital elements and the investigations involving the planet g are particularly interesting, since this body would lie in the Habitable Zone (HZ) of the star GJ581. This region,for this system, is very attractive of the dynamical point of view due to several resonances of two and three bodies present there. In this work, we investigate the conditions under which the planet g may exist. We stress the fact that the planet g is intimately related with the orbital elements of the planet d; more precisely, we conclude that it is not possible to disconnect its existence from the determination of the eccentricity of the planet d. Concerning the planet f, we have found one solution with period $\approx 450$ days, but we are judicious about any affirmation concernig this body because its signal is in the threshold of detection and the high period is in a spectral region where the ocorruence of aliases is very common. Besides, we outline some dynamical features of the habitable zone with the dynamical map and point out the role played by some resonances laying there.

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.

Bottke et al. (2007) linked the catastrophic formation of Baptistina Asteroid Family (BAF) to the K/T impact event. This linkage was based on dynamical and compositional evidence, which suggested the impactor had a composition similar to CM2 carbonaceous chondrites. However, our recent study (Reddy et al. 2009) suggests that the composition of (298) Baptistina is similar to LL-type ordinary chondrites rather than CM2 carbonaceous chondrites. This rules out any possibility of it being related to the source of the K/T impactor, if the impactor was of CM-type composition. Mineralogical study of asteroids in the vicinity of BAF has revealed a plethora of compositional types suggesting a complex formation and evolution environment. A detailed compositional analysis of 16 asteroids suggests several distinct surface assemblages including ordinary chondrites (Gaffey SIV subtype), primitive achondrites (Gaffey SIII subtype), basaltic achondrites (Gaffey SVII subtype and V-type), and a carbonaceous chondrite. Based on our mineralogical analysis we conclude that (298) Baptistina is similar to ordinary chondrites (LL-type) based on olivine and pyroxene mineralogy and moderate albedo. S-type and V-type in and around the vicinity of BAF we characterized show mineralogical affinity to (8) Flora and (4) Vesta and could be part of their families. Smaller BAF asteroids with lower SNR spectra showing only a 'single' band are compositionally similar to (298) Baptistina and L/LL chondrites. It is unclear at this point why the silicate absorption bands in spectra of asteroids with formal family definition seem suppressed relative to background population, despite having similar mineralogy.

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.

The subject of this paper is the secular behaviour of a pair of planets evolving under dissipative forces. In particular, we investigate the case when dissipative forces affect the planetary semi-major axes and the planets move inward/outward the central star, in a process known as planet migration. To perform this investigation, we introduce fundamental concepts of conservative and dissipative dynamics of the three-body problem. Based on these concepts, we develop a qualitative model of the secular evolution of the migrating planetary pair. Our approach is based on analysis of the energy and the orbital angular momentum exchange between the two-planet system and an external medium; thus no specific kind of dissipative forces is invoked. We show that, under assumption that dissipation is weak and slow, the evolutionary routes of the migrating planets are traced by the Mode I and Mode II stationary solutions of the conservative secular problem. The ultimate convergence and the evolution of the system along one of these secular modes of motion is determined uniquely by the condition that the dissipation rate is sufficiently smaller than the proper secular frequency of the system. We show that it is possible to reassemble the starting configurations and migration history of the systems on the basis of their final states and consequently to constrain the parameters of the physical processes involved.

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.

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.

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

Among the largest objects in the main belt, asteroid 4 Vesta is unique in showing a basaltic crust. It is also the biggest member of the Vesta family, which is supposed to originate from a large cratering event about 1 Gyr ago (Marzari et al. 1996). Most of the members of the Vesta family for which a spectral classification is available show a V-type spectra. Before the discovery of 1459 Magnya (Lazzaro et al. 2000) and of several V-type NEA (Xu 1995), all the known V-type asteroids were members of the Vesta family. Recently two V-type asteroids, 809 Lundia and 956 Elisa, (Florczak et al. 2002) have been discovered well outside the limits of the family, near the Flora family. We currently know 22 V-type asteroids outside the family, in the inner asteroid belt. In this work we investigate the possibility that these objects are former family members that migrated to their current positions via the interplay of Yarkovsky effect and nonlinear secular resonances. The main dynamical feature of 956 Elisa and 809 Lundia is that they are currently inside the 2(g-g6)+s-s6 (z2 by Milani and Knezevic, 1993) secular resonance. Our investigations show that members of the Vesta dynamical family may drift in three-body and weak secular resonances until they are captured in the strong z2 secular resonance. Only asteroids with diameters larger than 16 km can remain in one of the three-body or secular resonances long enough to reach the region of the z2 resonance. This two-step mechanism of capture into the z2 resonance could explain: i) the current resonant orbits of 956 Elisa and 809 Lundia, ii) why their size is significantly larger than that of the typical member of the Vesta family, and iii) provide a lower limit on the Vesta family age.

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.

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.

(abbreviated) We use a semi-numerical approach to study the secular behavior of a system composed of a central star and two massive planets in eccentric co-planar orbits. We show that the secular dynamics of this system can be described using only two parameters, the ratios of the semimajor axes and the planetary masses. The main dynamical features of the system are presented in geometrical pictures that allows us to investigate a large domain of the secular phase space of this three-body problem without time-expensive numerical integrations of the equations of motion, and without any restriction on the magnitude of the planetary eccentricities. We show that three different regimes of secular motion with respect to the secular angle $\Delta\varpi$ are possible: they are circulation, oscillation (around 0 and $180^\circ$), and high eccentricity libration in a nonlinear secular resonance. The first two regimes are a continuous extension of the classical linear secular perturbation theory; the last is a new feature, hitherto unknown, in the secular dynamics of the three-body problem.

We investigate the survivability of Trojan-type companions of Neptune during primordial radial migration of the giant planets Jupiter, Saturn, Uranus, and Neptune. Loss of Neptune Trojans during planetary migration is not a random diffusion process. Rather, losses occur almost exclusively during discrete episodes when Trojan particles are swept by secondary resonances associated with mean-motion commensurabilities of Uranus with Neptune. The single greatest episode of loss ejects nearly 75% of existing Neptune Trojans and occurs just prior to Neptune reaching its final orbit.

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.