Stationary Orbits in Resonant Extrasolar Planetary Systems

We present a catalog of stable and unstable apsidal corotation resonance (ACR) for the resonant planar planetary three-body problem, including both symmetric and asymmetric solutions. Calculations are performed with a new approach based on a numerical determination of the averaged Hamiltonian function. It has the advantage of being very simple to use and, with the exception of the immediate vicinity of the collision curve, yields precise results for any values of the eccentricities and semimajor axes. The present catalog includes results for the 3/2, 3/1, and 4/1 mean-motion resonances. The 5/1 and 5/2 commensurabilities are also discussed briefly. These results complement our previous results for the 2/1 (Beaugé et al. 2006, MNRAS, 365, 1160–1170), and give a broad picture of the structure of many important planetary resonances.     

On the 2/1 resonant planetary dynamics – periodic orbits and dynamical stability

The 2/1 resonant dynamics of a two-planet planar system is studied within the framework of the three-body problem by computing families of periodic orbits and their linear stability. The continuation of resonant periodic orbits from the restricted to the general problem is studied in a systematic way. Starting from the Keplerian unperturbed system, we obtain the resonant families of the circular restricted problem. Then, we find all the families of the resonant elliptic restricted three-body problem, which bifurcate from the circular model. All these families are continued to the general three-body problem, and in this way we can obtain a global picture of all the families of periodic orbits of a two-planet resonant system. The parametric continuation, within the framework of the general problem, takes place by varying the planetary mass ratio ρ. We obtain bifurcations which are caused either due to collisions of the families in the space of initial conditions or due to the vanishing of bifurcation points. Our study refers to the whole range of planetary mass ratio values  [ρ∈ (0, ∞)]  and, therefore we include the passage from external to internal resonances. Thus, we can obtain all possible stable configurations in a systematic way. As an application, we consider the dynamics of four known planetary systems at the 2/1 resonance and we examine if they are associated with a stable periodic orbit.     

Secular dynamics of the three-body problem: application to the υ Andromedae planetary system

The discovery of extra-solar planetary systems with multiple planets in highly eccentric orbits (∼0.1-0.6), in contrast with our own Solar System, makes classical secular perturbation analysis very limited. In this paper, we use a semi-numerical approach to study the secular behavior of a system composed of a central star and two massive planets in co-planar orbits. We show that the secular dynamics of this system can be described using only two parameters, the ratios of the semi-major 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 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. The topology of the phase space is also investigated in detail by means of spectral map techniques, which allow us to detect the separatrix of a non-linear secular apsidal resonance. Finally, the qualitative study is supplemented by direct numerical integrations. Three different regimes of secular motion with respect to the secular angle Δ? are possible: they are circulation, oscillation (around 0° and 180°), and high eccentricity libration in a non-linear 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 apply the analysis to the case of the two outer planets in the υ Andromedae system, and obtain its periodic and ordinary orbits, the general structure of its secular phase space, and the boundaries of its secular stability; we find that this system is secularly stable over a large domain of eccentricities. Applying this analysis to a wide range of planetary mass and semi-major axis ratios (centered about the υ Andromedae parameters), we find that apsidal oscillation dominates the secular phase space of the three-body coplanar system, and that the non-linear secular resonance is also a common feature.     

Confined chaotic motion in three-body resonances: trapping of trans-Neptunian material induced by gas-drag

Jiang & Yeh proposed gas-drag-induced resonant capture as a mechanism able to explain the dominant 3:2 resonance observed in the trans-Neptunian belt. Using a model of a disc–star–planet system they concluded that gaseous drag in a protoplanetary disc can trap trans-Neptunian object (TNO) embryos into the 3:2 resonance rather easily although it could not trap objects into the 2:1 resonance. Here we further investigate this scenario using numerical simulations within the context of the planar restricted four-body problem by including both present-day Uranus and Neptune. Our results show that mean motion and corotation resonances are possible and trapping into both the 3:2 and 2:1 resonances as well as other resonances is observed. The associated corotation centres may easily form larger planetesimals from smaller ones. Corotation resonances evolve into pure Lindblad resonances in a time-scale of 0.5 Myr. The non-linear corotation and mean motion resonances produced are very size selective. The 3:2 resonance is dominant for submetric particles but for larger particles the 2:1 resonance is stronger. In summary, our calculations show that confined chaotic motion around the resonances not only increases trapping efficiency but also the orbital eccentricities of the trapped material, modifying the relative abundance of trapped particles in different resonances. If we assume a more compact planetary system, instead of using the present-day values of the orbital elements of Uranus and Neptune, our results remain largely unchanged.     

The long-term stability of extrasolar system HD 37124. Numerical study of resonance effects

We describe numerical tools for the stability analysis of extrasolar planetary systems. In particular, we consider the relative Poincaré variables and symplectic integration of the equations of motion. We apply the tangent map to derive a numerically efficient algorithm of the fast indicator Mean Exponential Growth factor of Nearby Orbits (MEGNO), a measure of the maximal Lyapunov exponent, that helps to distinguish chaotic and regular configurations. The results concerning the three-planet extrasolar system HD 37124 are presented and discussed. The best-fitting solutions found in earlier works are studied more closely. The system involves Jovian planets with similar masses. The orbits have moderate eccentricities, nevertheless the best-fitting solutions are found in dynamically active region of the phase space. The long-term stability of the system is determined by a net of low-order two-body and three-body mean motion resonances. In particular, the three-body resonances may induce strong chaos that leads to self-destruction of the system after Myr of apparently stable and bounded evolution. In such a case, numerically efficient dynamical maps are useful to resolve the fine structure of the phase space and to identify the sources of unstable behaviour.     

Resonances in Multiple Planetary Systems

Since the first extrasolar planet was discovered about 10 years ago, a major point of dynamical investigations was the determination of stable regions in extrasolar planetary systems where additional planets may exist. Using numerical methods, we investigate the dynamical stability in known multiple planetary systems (HD74156, HD38529, HD168443, HD169830) with special interest on the region between the two known planets and on the mean motion resonances inside this region. As a dynamical model we take the restricted 4-body problem containing the host star, the two planets and massless test-planets. For our numerical integrations, we used the Lie-integrator and additionally the Fast Lyapunov Indicators as a tool for detecting chaotic motion. We also investigated the inner resonances with the outer planet and the outer resonances with the inner planet with test-planets located inside the resonances.     

Partial averaging near a resonance in planetary dynamics

An analytical treatment of the evolutionary dynamics of a three-body planetary system subject to dynamical friction with an interplanetary medium is presented. The analysis presented here is in connection with the results of numerical integrations of such systems recently published by Haghighipour. Using the method of partial averaging near a resonance, the dynamics of a restricted, circular, planar three-body system, with the inner body more massive, is studied and the time variation of quantities such as the orbital angular momentum and the eccentricity of the outer planet, which were previously obtained from numerical integrations, is analytically verified.     

Chirikov diffusion in the asteroidal three-body resonance (5, −2, −2)

The theory of diffusion in many-dimensional Hamiltonian system is applied to asteroidal dynamics. The general formulation developed by Chirikov is applied to the Nesvorny-Morbidelli analytic model of three-body (three-orbit) mean-motion resonances (Jupiter-Saturn-asteroid). In particular, we investigate the diffusion along and across the separatrices of the (5, −2, −2) resonance of the (490) Veritas asteroidal family and their relationship to diffusion in semi-major axis and eccentricity. The estimations of diffusion were obtained using the Melnikov integral, a Hadjidemetriou-type sympletic map and numerical integrations for times up to 10⁸ years.     

A numerical investigation of coorbital stability and libration in three dimensions

Motivated by the dynamics of resonance capture, we study numerically the coorbital resonance for inclination \(0\le I\le 180^\circ \) in the circular restricted three-body problem. We examine the similarities and differences between planar and three dimensional coorbital resonance capture and seek their origin in the stability of coorbital motion at arbitrary inclination. After we present stability maps of the planar prograde and retrograde coorbital resonances, we characterize the new coorbital modes in three dimensions. We see that retrograde mode I (R1) and mode II (R2) persist as we change the relative inclination, while retrograde mode III (R3) seems to exist only in the planar problem. A new coorbital mode (R4) appears in 3D which is a retrograde analogue to an horseshoe-orbit. The Kozai–Lidov resonance is active for retrograde orbits as well as prograde orbits and plays a key role in coorbital resonance capture. Stable coorbital modes exist at all inclinations, including retrograde and polar obits. This result confirms the robustness the coorbital resonance at large inclination and encourages the search for retrograde coorbital companions of the solar system’s planets.     

Analytical and numerical results on the stability of a planetary precessional model

A model for planetary precession is investigated using analytical and numerical techniques. A Hamiltonian function governing the model is derived in terms of action-angle Andoyer-Déprit variables under the assumption of equatorial symmetry. As a first approximation a simplified Hamiltonian with zero-eccentricity is considered and stability estimates are derived using KAM theory. A validation of the analytical results is performed computing Poincaré surfaces of section for the circular and elliptical model. We also investigate the role of the eccentricity and its connection with the appearance of resonances. Special attention is devoted to the particular case of the Earth–Moon system. This revised version was published online in July 2006 with corrections to the Cover Date.     

Motion of dust in mean motion resonances with planets

Effect of stellar electromagnetic radiation on the motion of spherical dust particle in mean motion orbital resonances with a planet is investigated. Planar circular restricted three-body problem with the Poynting–Robertson (P–R) effect yields monotonic secular evolution of eccentricity when the particle is trapped in the resonance. Planar elliptic restricted three-body problem with the P–R effect enables nonmonotonous secular evolution of eccentricity and the evolution of eccentricity is qualitatively consistent with the published results for the complicated case of interaction of electromagnetic radiation with nonspherical dust grain. Thus, it is sufficient to allow either nonzero eccentricity of the planet or nonsphericity of the grain and the orbital evolutions in the resonances are qualitatively equal for the two cases. This holds both for exterior and interior mean motion orbital resonances. Evolutions of argument of perihelion in the planar circular and elliptical restricted three-body problems are shown. Numerical integrations show that an analytic expression for the secular time derivative of the particle’s argument of perihelion does not exist, if only dependence on semimajor axis, eccentricity and argument of perihelion is admitted. Connection between the shift of perihelion and oscillations in secular eccentricity is presented for the planar elliptic restricted three-body problem with the P–R effect. Period of the oscillations corresponds to the period of one revolution of perihelion. Change of optical properties of the spherical grain with the heliocentric distance is also considered. The change of the optical properties: (i) does not have any significant influence on the secular evolution of eccentricity, (ii) causes that the shift of perihelion is mainly in the same direction/orientation as the particle motion around the Sun. The statements hold both for circular and noncircular planetary orbits.     

Stability of higher order resonances in the restricted three-body problem

Third and fourth order mean motion resonances are studied in the model of the restricted three-body problem by numerical methods for mass parameters corresponding approximately to the Sun?CJupiter and Sun?CNeptune systems. In the case of inner resonances, it is shown that there are two regions of libration in the 8:5 and 7:4 resonances, one at low, the other at high eccentricities. In the 9:5 and 7:3 resonances libration can exist only in one region at high eccentricities. The 5:2 and 4:1 resonances are very regular, with one librational zone existing for all eccentricities. There is no visible region of libration at any eccentricities in the 5:1 resonance, the transition between the regions of direct and retrograde circulation is very sharp. In the case of outer resonances, the 8:5 and 7:4 resonances have also two regions of libration, but the 9:5 resonance has three, the 7:3 resonance two librational zones. The 5:2 resonance is again very regular, but it is parted for two regions of libration at high eccentricities. Libration is possible in the 4:1 resonance only at high eccentricities. The 5:1 resonance is very symmetric. In the case of outer resonances, a comparison is made with trans-Neptunian objects (TNO) in higher order mean motion resonances. Several new librating TNOs are identified.     

A Review on Co-orbital Motion in Restricted and Planetary Three-body Problems

The 1:1 mean motion resonance may be referred to as the lowest order mean motion resonance in restricted or planetary three-body problems. The five well-known libration points of the circular restricted three-body problem are five equilibriums of the 1:1 resonance. Coorbital motion may take different shapes of trajectory. In case of small orbital eccentricities and inclinations, tadpole-shape and horseshoe-shape orbits are well-known. Other 1:1 libration modes different from the elementary ones can exist at moderate or large eccentricities and inclinations. Coorbital objects are not rare in our solar system, for example the Trojans asteroids and the coorbital satellite systems of Saturn. Recently, dozens of coorbital bodies have been identified among the near-Earth asteroids. These coorbital asteroids are believed to transit recurrently between different 1:1 libration modes mainly due to orbital precessions, planetary perturbations, and other possible effects. The Hamiltonian system and the Hill’s three-body problem are two effective approaches to study coorbital motions. To apply the perturbation theory to the Hamiltonian system, standard procedures involve the development of the disturbing function, averaging and normalization, theory of ideal resonance model, secular perturbation theory, etc. Global dynamics of coorbital motion can be revealed by the Hamiltonian approach with a suitable expansion. The Hill’s problem is particularly suitable for the studies on the relative motion of two coorbital bodies during their close encounter. The Hill’s equation derived from the circular restricted three-body problem is well known. However, the general Hill’s problem whose equation of motion takes exactly the same form applies to the non-restricted case where the mass of each body is non-negligible, namely the planetary case. The Hill’s problem can be transformed into a “canonical shape” so that the averaging principle can be applied to construct a secular perturbation theory. Besides the two analytical theories, numerical methods may be consulted, for example the approach of periodic orbit, the surface of section, and the computation of invariant manifolds carried by equilibriums or periodic orbits.     

Exoplanetary systems: The role of an equilibrium at high mutual inclination in shaping the global behavior of the 3-D secular planetary three-body problem

On the basis of a high-order (order 12) expansion of the perturbative potential in powers of the eccentricities and the inclinations, we analyze the secular interactions of two non-coplanar planets which are not in mean-motion resonance. The model is based on the planetary three-body problem which can be reduced to two degrees of freedom by the well-known elimination of the nodes [Jacobi, C.G.J., 1842. Astron. Nachr. XX, 81-102]. We introduce non-singular canonical variables which bring forward the symmetries of the problem. The main dynamical features depend on the location and stability of the equilibria which are easily found with our analytical model. We find that there exists an equilibrium when both eccentricities are zero. When the mutual inclination is small, this equilibrium is stable, but for larger mutual inclination it becomes unstable, generating a large chaotic zone and, by bifurcation, two regular regions, the so-called Kozai resonances. This analytical study which depends on only two parameters (the ratio of the semi-major axes and the mass ratio of the planets) makes possible a large survey of the problem and enables us to identify and quantify its main dynamical features, periodic orbits, regular and chaotic zones, etc. The results of our analytical model are illustrated and confirmed by numerical integrations.     

Resonant periodic orbits in the exoplanetary systems

The planetary dynamics of 4/3, 3/2, 5/2, 3/1 and 4/1 mean motion resonances is studied by using the model of the general three body problem in a rotating frame and by determining families of periodic orbits for each resonance. Both planar and spatial cases are examined. In the spatial problem, families of periodic orbits are obtained after analytical continuation of vertical critical orbits. The linear stability of orbits is also examined. Concerning initial conditions nearby stable periodic orbits, we obtain long-term planetary stability, while unstable orbits are associated with chaotic evolution that destabilizes the planetary system. Stable periodic orbits are of particular importance in planetary dynamics, since they can host real planetary systems. We found stable orbits up to 60^° of mutual planetary inclination, but in most families, the stability does not exceed 20^°–30^°, depending on the planetary mass ratio. Most of these orbits are very eccentric. Stable inclined circular orbits or orbits of low eccentricity were found in the 4/3 and 5/2 resonance, respectively.     

The resonant structure of Jupiter's Trojan asteroids – II. What happens for different configurations of the planetary system

In a previous paper, we have found that the resonance structure of the present Jupiter Trojan swarms could be split up into four different families of resonances. Here, in a first step, we generalize these families in order to describe the resonances occurring in Trojan swarms embedded in a generic planetary system. The location of these families changes under a modification of the fundamental frequencies of the planets and we show how the resonant structure would evolve during a planetary migration. We present a general method, based on the knowledge of the fundamental frequencies of the planets and on those that can be reached by the Trojans, which makes it possible to predict and localize the main events arising in the swarms during migration. In particular, we show how the size and stability of the Trojan swarms are affected by the modification of the frequencies of the planets. Finally, we use this method to study the global dynamics of the Jovian Trojan swarms when Saturn migrates outwards. Besides the two resonances found by Morbidelli et al. which could have led to the capture of the current population just after the crossing of the 2:1 orbital resonance, we also point out several sequences of chaotic events that can influence the Trojan population.     

On the stability of circumbinary planetary systems

The dynamics of circumbinary planetary systems (the systems in which the planets orbit a central binary) with a small binary mass ratio discovered to date is considered. The domains of chaotic motion have been revealed in the “pericentric distance–eccentricity” plane of initial conditions for the planetary orbits through numerical experiments. Based on an analytical criterion for the chaoticity of planetary orbits in binary star systems, we have constructed theoretical curves that describe the global boundary of the chaotic zone around the central binary for each of the systems. In addition, based on Mardling’s theory describing the separate resonance “teeth” (corresponding to integer resonances between the orbital periods of a planet and the binary), we have constructed the local boundaries of chaos. Both theoretical models are shown to describe adequately the boundaries of chaos on the numerically constructed stability diagrams, suggesting that these theories are efficient in providing analytical criteria for the chaoticity of planetary orbits.     

Capture of dust grains in exterior resonances with planets

The temporary capture of the dust grains in the exterior resonances with planets is studied in the frames of the planar circular three-body problem with Poynting-Robertson (PR) drag. For the Earth and particles ~ 10 Μm the resonances 4/5, 5/6, 6/7, 7/8 are shown to be most effective. The capture is only temporary (of order 10⁵ years) and the position of resonance may be calculated from semi-analytical model using averaged disturbing function. These semi-analytical results are confirmed by numerical integration. For various planet this picture changes as with increasing planetary mass the more exterior resonances become more important. We showed that for Jupiter (at least in the space between Jupiter and Saturn) the resonance 1/2 plays the dominant role. The capture time is here several myr but again eccentricity is evolving to eccentricity e ₀ ~ 0.48 of libration point for this resonance.     

Explicit Semianalytical Theory of Asteroid Motion

A new semianalytical theory of asteroid motion is presented. The theory is developed on the basis of Kaula's expansion of the disturbing function including terms up to the second order with respect to the masses of disturbing bodies. The theory is constructed in explicit form that gives the possibility to study separately the influence of different perturbations in the dynamics of minor planets. The mean-motion resonances with major planets as well as mixed three-body resonances can also be taken into account. For the non-resonant case the formulas obtained can be used for deriving the second transformation to calculate the proper elements of an asteroid orbit in closed form with respect to inclinations and eccentricities. This revised version was published online in July 2006 with corrections to the Cover Date.     

Reducing the probability of capture into resonance

A migrating planet can capture planetesimals into mean motion resonances. However, resonant trapping can be prevented when the drift or migration rate is sufficiently high. Using a simple Hamiltonian system for first- and second-order resonances, we explore how the capture probability depends on the order of the resonance, drift rate and initial particle eccentricity. We present scaling factors as a function of the planet mass and resonance strength to estimate the planetary migration rate above which the capture probability drops to less than half. Applying our framework to multiple extrasolar planetary systems that have two planets locked in resonance, we estimate lower limits for the outer planet's migration rate, allowing resonance capture of the inner planet.
Mean motion resonances are comprised of multiple resonant subterms. We find that the corotation subterm can reduce the probability of capture when the planet eccentricity is above a critical value. We present factors that can be used to estimate this critical planet eccentricity. Applying our framework to the migration of Neptune, we find that Neptune's eccentricity is near the critical value that would make its 2 : 1 resonance fail to capture twotinos. The capture probability is affected by the separation between resonant subterms and so is also a function of the precession rates of the longitudes of periapse of both planet and particle near resonance.