Families of periodic collision orbits in the general three-body problem

In the general three-body problem, in a rotating frame of reference, a symmetric periodic solution with a binary collision is determined by the abscissa of one body and the energy of the system. For different values of the masses of the three bodies, the symmetric periodic collision orbits form a two-parametric family. In the case of equal masses of the two bodies and small mass of the third body, we found several symmetric periodic collision orbits similar to the corresponding orbits in the restricted three-body problem. Starting with one symmetric periodic collision orbit we obtained two families of such orbits. Also starting with one collision orbit in the Sun-Jupiter-Saturn system we obtained, for a constant value of the mass ratio of two bodies, a family of symmetric periodic collision orbits.     

Planets in the asteroid belt

Abstract— The main asteroid belt has lost >99.9% of its solid mass since the time at which the planets were forming, according to models for the protoplanetary nebula. Here we show that the primordial asteroid belt could have been cleared efficiently if much of the original mass accreted to form planetsized bodies, which were capable of perturbing one another into unstable orbits. We provide results from 25 N‐body integrations of up to 200 planets in the asteroid belt, with individual masses in the range 0.017–0.33 Earth masses. In the simulations, these bodies undergo repeated close encounters which scatter one another into unstable resonances with the giant planets, leading to collision with the Sun or ejection from the solar system. In response, the giant planets' orbits migrate radially and become more circular. This reduces the size of the main‐belt resonances and the clearing rate, although clearing continues. If ～3 Earth masses of material was removed from the belt this way, Jupiter and Saturn would initially have had orbital eccentricities almost twice their current values. Such orbits would have made Jupiter and Saturn 10–100x more effective at clearing material from the belt than they are on their current orbits. The time required to remove 90% of the initial mass from the belt depends sensitively on the giant planets' orbits, and weakly on the masses of the asteroidal planets. 18 of the 25 simulations end with no planets left in the belt, and the clearing takes up to several hundred million years. Typically, the last one or two asteroidal planets are removed by interactions with planets in the terrestrial region     

Searching for the source of short-period comet nuclei: The direction of the spatial migration of comets

An attempt is made to determine the spatial location of the main source of short-period comet nuclei. Numerical calculations for the orbital evolution of Jupiter family comets, medium-period comets, and Centaurs are used to show that the orbits of small solar system bodies tend to evolve in the direction of increasing semimajor axes. This relates to bodies that can experience encounters with planets and whose orbital evolution is shaped by gravitational perturbations. It is concluded that there is good reason to search for the main source of the nuclei of Jupiter family comets at distances of 6 AU or less from the sun.     

Temporary satellite capture of comets by Jupiter

This paper studies the dynamical evolution of 97 Jupiter-family comets over an 800-year time period. More than two hundred encounters with Jupiter are investigated, with the observed comets moving during a certain period of time in an elliptic jovicentric orbit. In most cases this is an ordinary temporary satellite capture of a comet in Everhart??s sense, not associated with a transition of the small body into Jupiter??s family of satellites. The phenomenon occurs outside the Hill sphere with comets with a high Tisserand constant relative to Jupiter; the comets?? orbits have a small inclination to the ecliptic plane. An analysis of 236 encounters has allowed the determination within the planar pair two-body problem of a region of orbits in the plane (a, e) whose semimajor axes and eccentricities contribute to the phenomenon under study. Comets with orbits belonging to this region experience a temporary satellite capture during some of their encounters; the jovicentric distance function has several minima; and the encounters are characterized by reversions of the line of apsides and some others features of their combination that are intrinsic to comets in this region. Therefore, this region is called a region of comets with specific features in their encounters with Jupiter. Twenty encounters (out of 236), whereby the comet enters an elliptic jovicentric orbit in the Hill sphere, are identified and investigated. The size and shape of the elliptic heliocentric orbits enabling this transition are determined. It is found that in 11 encounters the motion of small bodies in the Hill sphere has features the most important of which is multiple minima of the jovicentric distance function. The study of these 20 encounters has allowed the introduction of the concept of temporary gravitational capture of a small body into the Hill sphere. An analysis of variations in the Tisserand constant in these (20) encounters of the observable comets shows that their motion is unstable in Hill??s sense.     

Periodic orbits of the elliptic restricted problem for the Sun-Jupiter-Saturn system

A systematic approach to generate periodic orbits in the elliptic restricted problem of three bodies in introduced. The approach is based on (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to periodic orbits in the elliptic restricted problem. Two families of periodic orbits of the elliptic restricted problem are found by this approach. The mass ratio of the primaries of these orbits is equal to that of the Sun-Jupiter system. The sidereal mean motions between the infinitesimal body and the smaller primary are in a 2:5 resonance, so as to approximate the Sun-Jupiter-Saturn system. The linear stability of these periodic orbits are studied as functions of the eccentricities of the primaries and of the infinitesimal body. The results show that both stable and unstable periodic orbits exist in the elliptic restricted problem that are close to the actual Sun-Jupiter-Saturn system. However, the periodic orbit closest to the actual Sun-Jupiter-Saturn system is (linearly) stable.     

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.     

Dynamical evolution of weakly disturbed two-planetary system on cosmogonic time scales: The Sun-Jupiter-Saturn system

In the present paper, we used the Hori-Deprit method to construct the averaged Hamiltonian of the two-planetary problem accurate to the second order of a small parameter, the generating function of the transform, the change of variables formulas, and the right-hand sides of the equations in average elements. The evolution of the two-planet Sun-Jupiter-Saturn system was investigated by numerical integration over 10 billion years. The motion of the planets has an almost periodic character. The eccentricities and inclinations of Jupiter’s and Saturn’s orbits remain small but different from zero. The short-term disturbances remain small over the entire period considered in the study.     

Periodic orbits of the general three-body problem for the Sun-Jupiter-Saturn system

Two families of symmetric periodic orbits of the planar, general, three-body problem are presented. The masses of the three bodies include ratios equal to the Sun-Jupiter-Saturn system and the periods of the orbits of Jupiter and Saturn are in a 25 resonance. The (linear) stability of the orbits are studied in relation to eccentricity and mass variations. The generation of the two families of periodic orbits follows a systematic approach and employs (numerical) continuation from periodic orbits of the first and second kind in the circular restricted problem to the elliptic restricted problem and from the circular and elliptic problems to the general problem through bifurcation phenomena relating the three dynamical systems. The approach also provides insight into the evolutionary process of periodic orbits continued from the restricted problems to the general problem.     

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.     

Computer studies of the evolution of planetary and satellite systems

This paper describes two computer experiments carried out with a CDC-Cyber 74 program (Barricelliet al., 1979) for computer simulation of a large number of objects in orbit about a central body or primary. The first experiment was started with 125 planets of which the two largest ones had coplanar orbits and masses comparable to those of Jupiter and Saturn, respectively. Their semimajor axes and eccentricities were, however, much larger. The smaller planets had a distribution promoting the formation of an axial meeting area. The experiment gives information relevant to the question of focusing of planetary orbits into a common plane and to the question of the formation and stability of an axial meeting area. Together with the next experiment, it also gives information about the development of commensurabilities (or resonances) with the largest planets.The second experiment started with 55 planets none of them with a mass greater than about 20% of Jupiter's but several of them with orbits close to a common plane. The aim of the experiment was to investigate whether successive captures followed by planetary fusion (Barricelli, 1972a) could lead to the formation of major planets comparable to Jupiter and Saturn, and in similar orbits. Also this experiment gives information relevant to the commensurability problem.     

Planet mass stability margin of two-planetary systems

The Lagrange stability of the Sun-Jupiter-Saturn and 47 UMa two-planetary systems at a time scale of 10⁶ yr was studied using the method of averaging. When the masses of Jupiter and Saturn increase by 19 times, these planets can closely converge. The study of Lagrange stability in the case of successive mass increase allows for the obtainment of upper estimates of possible masses of extrasolar planets. Close approachs in the 47 UMa system are possible when minimal masses increase by 38 times. approachs are revealed when analyzing osculating elements; in averaged elements, approachs are absent. Resonant properties of six extrasolar two-planetary systems where the outer planet is less massive than the inner one are studied. The values of semi-major axes of planet orbits in HD 82943 and HD 73526 systems lie in a narrow resonant zone; in 47 UMa, μ Ara and HD 108874 systems lie in a wide resonant zone. In the HD 12661, the system resonances of a lower order were not revealed.     

Eccentric Extrasolar Planets: The Jumping Jupiter Model

Most extrasolar planets discovered to date are more massive than Jupiter, in surprisingly small orbits (semimajor axes less than 3 AU). Many of these have significant orbital eccentricities. Such orbits may be the product of dynamical interactions in multiplanet systems. We examine outcomes of such evolution in systems of three Jupiter-mass planets around a solar-mass star by integration of their orbits in three dimensions. Such systems are unstable for a broad range of initial conditions, with mutual perturbations leading to crossing orbits and close encounters. The time scale for instability to develop depends on the initial orbital spacing; some configurations become chaotic after delays exceeding 10⁸ y. The most common outcome of gravitational scattering by close encounters is hyperbolic ejection of one planet. Of the two survivors, one is moved closer to the star and the other is left in a distant orbit; for systems with equal-mass planets, there is no correlation between initial and final orbital positions. Both survivors may have significant eccentricities, and the mutual inclination of their orbits can be large. The inner survivor's semimajor axis is usually about half that of the innermost starting orbit. Gravitational scattering alone cannot produce the observed excess of “hot Jupiters” in close circular orbits. However, those scattered planets with large eccentricities and small periastron distances may become circularized if tidal dissipation is effective. Most stars with a massive planet in an eccentric orbit should have at least one additional planet of comparable mass in a more distant orbit.     

Dynamical evolution of a cometary swarm in the outer planetary region

The dynamical evolution of bodies under the gravitational influence of the accreting proto-Uranus and proto-Neptune is investigated. The main aim of this study is to analyze the interrelations between the accretion of Uranus and Neptune with other processes of cosmological importance as, for example, the formation of a cometary reservoir from bodies placed into near-parabolic orbits by planetary perturbations and the scattering of bodies to the region of the terrestrial planets. Starting with a mass ratio (initial mass/present mass) of 0.1, Uranus and Neptune acquire masses close to their present ones in a time scale of 10⁸ years. Neptune is found to be the most important contributor of comets to the cometary reservoir. The time scale of bodies scattered by Neptune to reach near-parabolic orbits (semimajor axes a > 10⁴ AU)is about 10⁹ years. The contribution of Uranus was partially inhibited because a large part of the residual bodies of its accretion zone fell under the strong gravitational influence of Jupiter and Saturn. A significant fraction of the bodies dispersed by Uranus and Neptune reached the region of the terrestrial planets in a time scale of some 10⁸ years.     

High-resolution simulations of the final assembly of Earth-like planets I. Terrestrial accretion and dynamics

The final stage in the formation of terrestrial planets consists of the accumulation of ∼1000-km “planetary embryos” and a swarm of billions of 1-10 km “planetesimals.” During this process, water-rich material is accreted by the terrestrial planets via impacts of water-rich bodies from beyond roughly 2.5 AU. We present results from five high-resolution dynamical simulations. These start from 1000-2000 embryos and planetesimals, roughly 5-10 times more particles than in previous simulations. Each simulation formed 2-4 terrestrial planets with masses between 0.4 and 2.6 Earth masses. The eccentricities of most planets were ∼0.05, lower than in previous simulations, but still higher than for Venus, Earth and Mars. Each planet accreted at least the Earth's current water budget. We demonstrate several new aspects of the accretion process: (1) The feeding zones of terrestrial planets change in time, widening and moving outward. Even in the presence of Jupiter, water-rich material from beyond 2.5 AU is not accreted for several millions of years. (2) Even in the absence of secular resonances, the asteroid belt is cleared of >99% of its original mass by self-scattering of bodies into resonances with Jupiter. (3) If planetary embryos form relatively slowly, then the formation of embryos in the asteroid belt may have been stunted by the presence of Jupiter. (4) Self-interacting planetesimals feel dynamical friction from other small bodies, which has important effects on the eccentricity evolution and outcome of a simulation.     

Orbital Evolution of the Sun–Jupiter–Saturn–Uranus–Neptune Four-Planet System on Long-Time Scales

The four-planet problem is solved by constructing an averaged semi-analytical theory of secondorder motion by planetary masses. A discussion is given of the results obtained by numerical integration of the averaged equations of motion for the Sun–Jupiter–Saturn–Uranus–Neptune system over a time interval of 10 Gyr. The integration is based on high-order Runge–Kutta and Everhart methods. The motion of the planets is almost periodic in nature. The eccentricities and inclinations of the planetary orbits remain small. Short-period perturbations remain small over the entire interval of integration. Conclusions are drawn about the resonant properties of the motion. Estimates are given for the accuracy of the numerical integration.     

Chaotic dynamics of planet‐encountering bodies

The dynamics of two families of minor inner solar system bodies that suffer frequent close encounters with the planets is analyzed. These families are: Jupiter family comets (JF comets) and Near Earth Asteroids (NEAs). The motion of these objects has been considered to be chaotic in a short time scale,and the close encounters are supposed to be the cause of the fast chaos. For a better understanding of the chaotic behavior we have computed Lyapunov Characteristic Exponents (LCEs) for all the observed members of both populations. LCEs are a quantitative measure of the exponential divergence of initially close orbits. We have observed that most members of the two families show a concentration of Lyapunov times (inverse of LCE) around 50–100yr. The concentration is more pronounced for JF comets than for NEAs, among which a lesser spread is observed for those that actually cross the Earth's orbit (mean perihelion distance q < 1.05 AU). It is also observed that a general correspondence exists between Lyapunov times and the time between consecutive encounters. A simple model is introduced to describe the basic characteristics of the dynamical evolution. This model considers an impulsive approach, where the particles evolve unperturbedly between encounters and suffer ‘kicks’ in semimajor axis at the encounters. It also reproduces successfully the short Lyapunov times observed in the numerical integrations and is able to estimate the dynamical lifetimes of comets during a stay in the Jupiter family in correspondence with previous estimates. It has been demonstrated with the model that the encounters with the largest effect on the exponential growth of the distance between initially nearby orbits are neither the infrequent deep encounters, nor the frequent and far ones; instead, the intermediate approaches have the most relevant contribution to the error growth. Such encounters are at a distance a few times the radius of the Hill's sphere of the planet (e.g. 3). An even simpler model allows us to get analytical estimates of the Lyapunov times in good agreement with the values coming from the model above and the numerical integrations. The predictability of the medium‐term evolution and the hazard posed to the Earth by those objects are analysed in the Discussion section. This revised version was published online in July 2006 with corrections to the Cover Date.     

The provenance and evolution of comets

The process of comet formation through the hierarchical aggregation of originally submicron-sized interstellar grains to form micron-sized particles and then larger bodies in the protoplanetary disc, culminating in the formation of planetesimals in the disc extending from Jupiter to beyond Neptune, is briefly reviewed. The ‘planetesimal’ theory for the origin of comets implies the existence of distinct cometary reservoirs, with implications for the immediate provenance of observed comets (both long-period and short-period) and their evolution as a result of planetary perturbations and physical decay, for example splitting and sublimation. The principal mode of cometary decay and collisional interaction with the terrestrial planets is through the formation and evolution of streams of cometary debris and hitherto undiscovered ‘families’ of cometary asteroids. Recent dynamical results, in particular the sungrazing and sun-colliding end-state for short-period comet and asteroid orbits, are briefly discussed.     

Building the terrestrial planets: Constrained accretion in the inner Solar System

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include: (1) the orbits, in particular the small eccentricities, and (2) the masses of the terrestrial planets - Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; (3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; (4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and (5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N = 1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the “Nice model”. We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars’ mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. Our simulations leave us with a problem: we can reasonably satisfy the observed constraints (except for Earth’s water) with a configuration of Jupiter and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which are consistent with the Nice model fail to reproduce Mars’ small size.     

A two-parameter survey of periodic orbits in the restricted problem of three bodies

Within the context of the restricted problem of three bodies, we wish to show the effects, caused by varying the mass ratio of the primaries and the eccentricity of their orbits, upon periodic orbits of the infinitesimal mass that are numerical continuations of circular orbits in the ordinary problem of two bodies. A recursive-power-series technique is used to integrate numerically the equations of motion as well as the first variational equations to generate a two-parameter family of periodic orbits and to identify the linear stability characteristics thereof. Seven such families (comprised of a total of more than 2000 orbits) with equally spaced mass ratios from 0.0 to 1.0 and eccentricities of the orbits of the primaries in a range 0.0 to 0.6 are investigated. Stable orbits are associated with large distances of the infinitesimal mass from the perturbing primary, with nearly circular motion of the primaries, and, to a slightly lesser extent, with small mass ratios of the primaries.Conversely, unstable orbits for the infinitesimal mass are associated with small distances from the perturbing primary, with highly elliptic orbits of the primaries, and with large mass ratios.