Dynamics of the 3/1 planetary mean-motion resonance: an application to the HD60532 b-c planetary system

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. The case where the outer planet is more massive than its inner companion is considered. We show that the resonant dynamics can be described using two fundamental parameters, the total angular momentum and the spacing parameter. The topology of the Hamiltonian function describing the resonant behaviour is investigated on a large domain of the phase space 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 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 the 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.     

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.     

On the dynamics of extrasolar planetary systems under dissipation: Migration of planets

We study the dynamics of planetary systems with two planets moving in the same plane, when frictional forces act on the two planets, in addition to the gravitational forces. The model of the general three-body problem is used. Different laws of friction are considered. The topology of the phase space is essential in understanding the evolution of the system. The topology is determined by the families of stable and unstable periodic orbits, both symmetric and non symmetric. It is along the stable families, or close to them, that the planets migrate when dissipative forces act. At the critical points where the stability along the family changes, there is a bifurcation of a new family of stable periodic orbits and the migration process changes route and follows the new stable family up to large eccentricities or to a chaotic region. We consider both resonant and non resonant planetary systems. The 2/1, 3/1 and 3/2 resonances are studied. The migration to larger or smaller eccentricities depends on the particular law of friction. Also, in some cases the semimajor axes increase and in other cases they are stabilized. For particular laws of friction and for special values of the parameters of the frictional forces, it is possible to have partially stationary solutions, where the eccentricities and the semimajor axes are fixed.     

Dynamics Of 2/1 Resonant Extrasolar Systems Application to HD82943 and GLIESE876

A complete study is made of the resonant motion of two planets revolving around a star, in the model of the general planar three body problem. The resonant motion corresponds to periodic motion of the two planets, in a rotating frame, and the position and stability properties of the periodic orbits determine the topology of the phase space and consequently play an important role in the evolution of the system. Several families of symmetric periodic orbits are computed numerically, for the 2/1 resonance, and for the masses of some observed extrasolar planetary systems. In this way we obtain a global view of all the possible stable configurations of a system of two planets. These define the regions of the phase space where a resonant extrasolar system could be trapped, if it had followed in the past a migration process.The factors that affect the stability of a resonant system are studied. For the same resonance and the same planetary masses, a large value of the eccentricities may stabilize the system, even in the case where the two planetary orbits intersect. The phase of the two planets (position at perihelion or aphelion when the star and the two planets are aligned) plays an important role, and the change of the phase, other things being the same, may destabilize the system. Also, the ratio of the planetary masses, for the same total mass of the two planets, plays an important role and the system, at some resonances and some phases, is destabilized when this ratio changes.The above results are applied to the observed extrasolar planetary systems HD 82943, Gliese 876 and also to some preliminary results of HD 160691. It is shown that the observed configurations are close to stable periodic motion.     

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.     

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.     

Mean Motion Resonances, Gas Drag, and Supersonic Planetesimals in the Solar Nebula

We examine the orbital evolution of planetesimals under the influence of Jupiter's perturbations and nebular gas drag, under the assumption that gas persisted in the asteroid region for some time after Jupiter attained its final mass. Two distinct mechanisms, associated with the 2 : 1 and 3 : 2 mean motion resonances, can excite eccentricities to high values, despite the damping effect of drag. If Jupiter's eccentricity was comparable to its present value, planetesimals can be temporarily trapped in the 2 : 1 resonance. Bodies crossing the 3 : 2 resonance can enter a region of phase space with overlapping high-order resonances. Both mechanisms can produce eccentricities greater than 0.5 for asteroid-sized planetesimals. The combination of resonant perturbations and drag causes secular decay of semimajor axes, resulting in migration of bodies from the outer to inner belt. Inclinations remain low, implying significant collisional evolution during this migration. Velocities of resonant bodies relative to the gas are highly supersonic; these would have been a source of shock waves in the solar nebula.This revised version was published online in October 2005 with corrections to the Cover Date.     

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.     

The 1/1 resonance in extrasolar planetary systems

We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.     

Capture of grains into resonances through Poynting-Robertson drag

We review here some relevant problems connected to the evolution of circumstellar dust grains, subjected to Poynting-Robertson (PR) drag, and perturbed by first-order resonances with a planet on a circular orbit. We show that only outer mean motion resonances are able to counteract the damping effect of PR drag. However, the high orbital eccentricities reached by the particle lead to orbit crossings with the planet. This is a serious difficulty for a permanent trapping to be achieved. In any case, we show that the time spent in the resonance is long enough for statistical effects (accumulation at the resonant radius) to be significant. We underline some difficulties associated with this problem, namely, the non-adiabaticity of motion in the resonance phase space and the existence of close encounters with the planet at high eccentricities.     

Symmetric and asymmetric 3:1 resonant periodic orbits with an application to the 55Cnc extra-solar system

We study the dynamics of 3:1 resonant motion for planetary systems with two planets, based on the model of the general planar three body problem. The exact mean motion resonance corresponds to periodic motion (in a rotating frame) and the basic families of symmetric and asymmetric periodic orbits are computed. Four symmetric families bifurcate from the family of circular orbits of the two planets. Asymmetric families bifurcate from the symmetric families, at the critical points, where the stability character changes. There exist also asymmetric families that are independent of the above mentioned families. Bounded librations exist close to the stable periodic orbits. Therefore, such periodic orbits (symmetric or asymmetric) determine the possible stable configurations of a 3:1 resonant planetary system, even if the orbits of the two planets intersect. For the masses of the system 55Cnc most of the periodic orbits are unstable and they are associated with chaotic motion. There exist however stable symmetric and asymmetric orbits, corresponding to regular trajectories along which the critical angles librate. The 55Cnc extra-solar system is located in a stable domain of the phase space, centered at an asymmetric periodic orbit.     

Symmetric and asymmetric librations in extrasolar planetary systems: a global view

We present a global view of the resonant structure of the phase space of a planetary system with two planets, moving in the same plane, as obtained from the set of the families of periodic orbits. An important tool to understand the topology of the phase space is to determine the position and the stability character of the families of periodic orbits. The region of the phase space close to a stable periodic orbit corresponds to stable, quasi periodic librations. In these regions it is possible for an extrasolar planetary system to exist, or to be trapped following a migration process due to dissipative forces. The mean motion resonances are associated with periodic orbits in a rotating frame, which means that the relative configuration is repeated in space. We start the study with the family of symmetric periodic orbits with nearly circular orbits of the two planets. Along this family the ratio of the periods of the two planets varies, and passes through rational values, which correspond to resonances. At these resonant points we have bifurcations of families of resonant elliptic periodic orbits. There are three topologically different resonances: (1) the resonances (n + 1):n, (2:1, 3:2, ...), (2) the resonances (2n + 1):(2n-1), (3:1, 5:3, ...) and (3) all other resonances. The topology at each one of the above three types of resonances is studied, for different values of the sum and of the ratio of the planetary masses. Both symmetric and asymmetric resonant elliptic periodic orbits exist. In general, the symmetric elliptic families bifurcate from the circular family, and the asymmetric elliptic families bifurcate from the symmetric elliptic families. The results are compared with the position of some observed extrasolar planetary systems. In some cases (e.g., Gliese 876) the observed system lies, with a very good accuracy, on the stable part of a family of resonant periodic orbits.     

The 1/1 resonance in extrasolar systems

We present families of symmetric and asymmetric periodic orbits at the 1/1 resonance, for a planetary system consisting of a star and two small bodies, in comparison to the star, moving in the same plane under their mutual gravitational attraction. The stable 1/1 resonant periodic orbits belong to a family which has a planetary branch, with the two planets moving in nearly Keplerian orbits with non zero eccentricities and a satellite branch, where the gravitational interaction between the two planets dominates the attraction from the star and the two planets form a close binary which revolves around the star. The stability regions around periodic orbits along the family are studied. Next, we study the dynamical evolution in time of a planetary system with two planets which is initially trapped in a stable 1/1 resonant periodic motion, when a drag force is included in the system. We prove that if we start with a 1/1 resonant planetary system with large eccentricities, the system migrates, due to the drag force, along the family of periodic orbits and is finally trapped in a satellite orbit. This, in principle, provides a mechanism for the generation of a satellite system: we start with a planetary system and the final stage is a system where the two small bodies form a close binary whose center of mass revolves around the star.     

Modeling the 3-D secular planetary three-body problem: Discussion on the outer υ Andromedae planetary system

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 Δ?, possible motions are circulations, oscillations (around 0° and 180°), and high-eccentricity/inclination librations in secular resonances. With respect to the arguments of pericenter, ω₁ and ω₂, possible motions are direct circulation and high-inclination libration around ±90° 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 υ 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.     

Inclined asymmetric librations in exterior resonances

Librational motion in Celestial Mechanics is generally associated with the existence of stable resonant configurations and signified by the existence of stable periodic solutions and oscillation of critical (resonant) angles. When such an oscillation takes place around a value different than 0 or \(\pi \), the libration is called asymmetric. In the context of the planar circular restricted three-body problem, asymmetric librations have been identified for the exterior mean motion resonances (MMRs) 1:2, 1:3, etc., as well as for co-orbital motion (1:1). In exterior MMRs the massless body is the outer one. In this paper, we study asymmetric librations in the three-dimensional space. We employ the computational approach of Markellos (Mon Not R Astron Soc 184:273–281,  https://doi.org/10.1093/mnras/184.2.273, 1978) and compute families of asymmetric periodic orbits and their stability. Stable asymmetric periodic orbits are surrounded in phase space by domains of initial conditions which correspond to stable evolution and librating resonant angles. Our computations were focused on the spatial circular restricted three-body model of the Sun–Neptune–TNO system (TNO = trans-Neptunian object). We compare our results with numerical integrations of observed TNOs, which reveal that some of them perform 1:2 resonant, inclined asymmetric librations. For the stable 1:2 TNO librators, we find that their libration seems to be related to the vertically stable planar asymmetric orbits of our model, rather than the three-dimensional ones found in the present study.     

Gravitational radiation damping and the three-body problem

A model of three-body motion is developed which includes the effects of gravitational radiation reaction. The radiation reaction due to the emission of gravitational waves is the only post-Newtonian effect that is included here. For simplicity, all of the motion is taken to be planar. Two of the masses are viewed as a binary system, and the third mass, whose motion will be a fixed orbit around the centre-of-mass of the binary system, is viewed as a perturbation. This model aims to describe the motion of a relativistic binary pulsar that is perturbed by a third mass. Numerical integration of this simplified model reveals that, given the right initial conditions and parameters, one can see resonances. These ( m , n ) resonances are defined by the resonance condition,   mω =2 n Ω  , where m and n are relatively prime integers, and ω and Ω are the angular frequencies of the binary orbit and third mass orbit (around the centre-of-mass of the binary), respectively. The resonance condition consequently fixes a value for the semimajor axis of the binary orbit for the duration of the resonance; therefore the binary energy remains constant on average, while its angular momentum changes during the resonance.     

A symplectic mapping model for the 3/4 exterior

We present a symplectic mapping model to study the evolution of a small body at the 3/4 exterior resonance with Neptune, for planar and for three dimensional motion. The mapping is based on the averaged Hamiltonian close to this resonance and is constructed in such a way that the topology of its phase space is similar to that of the Poincaré map of the elliptic restricted three-body problem. Using this model we study the evolution of a small object near the 3/4 resonance. Both chaotic and regular motions are found, and it is shown that the initial phase of the object plays an important role on the appearance of chaos. In the planar case, objects that are phase-protected from close encounters with Neptune have regular orbits even at eccentricities up to 0.44. On the other hand objects that are not phase protected show chaotic behaviour even at low eccentricities. The introduction of the inclination to our model affects the stable areas around the 3/4 mean motion resonance, which now become thinner and thinner and finally at is=10° the whole resonant region becomes chaotic. This may justify the absence of a large population of objects at this resonance.     

Capture into first-order resonances and long-term stability of pairs of equal-mass planets

Massive planets form within the lifetime of protoplanetary disks, and therefore, they are subject to orbital migration due to planet–disk interactions. When the first planet reaches the inner edge of the disk, its migration stops and consequently the second planet ends up locked in resonance with the first one. We detail how the resonant trapping works comparing semi-analytical formulae and numerical simulations. We restrict to the case of two equal-mass coplanar planets trapped in first-order resonances, but the method can be easily generalized. We first describe the family of resonant stable equilibrium points (zero-amplitude libration orbits) using series expansions up to different orders in eccentricity as well as a non-expanded Hamiltonian. Then we show that during convergent migration the planets evolve along these families of equilibrium points. Eccentricity damping from the disk leads to a final equilibrium configuration that we predict precisely analytically. The fact that observed multi-exoplanetary systems are rarely seen in resonances suggests that in most cases the resonant configurations achieved by migration become unstable after the removal of the protoplanetary disk. Here we probe the stability of the resonances as a function of planetary mass. For this purpose, we fictitiously increase the masses of resonant planets, adiabatically maintaining the low-amplitude libration regime until instability occurs. We discuss two hypotheses for the instability, that of a low-order secondary resonance of the libration frequency with a fast synodic frequency of the system, and that of minimal approach distance between planets. We show that secondary resonances do not seem to impact resonant systems at low amplitude of libration. Resonant systems are more stable than non-resonant ones for a given minimal distance at close encounters, but we show that the latter nevertheless play the decisive role in the destabilization of resonant pairs. We show evidence that as the planetary mass increases and the minimal distance between planets gets smaller in terms of mutual Hill radius, the region of stability around the resonance center shrinks, until the equilibrium point itself becomes unstable.     

Chaotic transitions in resonant asteroidal dynamics

The utilization of chaotic dynamics approaches allowed the identification of many modes of motion in resonant asteroidal dynamics. As these dynamical systems are not integrable, the motion modes are not separated and one orbit may transit from one mode to another. In some cases, as in the \31 resonance, these transitions may lead, in a relatively short time scale, to eccentricities so high that the asteroid may approach the Sun and be destroyed. In the \21 and \32 resonances these transitions are much slower and only indirect estimations of the time which is needed for a generic asteroid to leave the resonance are possible. It may reach hundreds of million years in the more robust regions of the \21 resonance and a time of the order of billions of years in those of the \32 resonance. These values are consistent with the observed depletion of the \21 resonance (only a few asteroids known while almost 60 asteroids are known in the \32 resonance).     

Symmetric and Nonsymmetric Periodic Orbits in the Exterior Mean Motion Resonances with Neptune

We present families of periodic orbits and their stability for the exterior mean motion resonances 1:2, 1:3 and 1:4 with Neptune in the framework of the planar circular restricted three-body problem. We found that in each resonance there exist two branches of symmetric elliptic periodic orbits with stable and unstable segments. Asymmetric periodic orbits bifurcate from the corresponding symmetric ones. Asymmetric periodic orbits are stable and the motion in their neighbourhood is a libration with respect to the resonant angle variable. In all the families of asymmetric periodic orbits the eccentricity extends to high values. Poincaré sections reveal the changes of the topology in phase space.