Orbital stability of systems of closely-spaced planets,II: configurations with coorbital planets

We numerically investigate the stability of systems of 1 \({{\rm M}_{\oplus}}\) planets orbiting a solar-mass star. The systems studied have either 2 or 42 planets per occupied semimajor axis, for a total of 6, 10, 126, or 210 planets, and the planets were started on coplanar, circular orbits with the semimajor axes of the innermost planets at 1 AU. For systems with two planets per occupied orbit, the longitudinal initial locations of planets on a given orbit were separated by either 60° (Trojan planets) or 180°. With 42 planets per semimajor axis, initial longitudes were uniformly spaced. The ratio of the semimajor axes of consecutive coorbital groups in each system was approximately uniform. The instability time for a system was taken to be the first time at which the orbits of two planets with different initial orbital distances crossed. Simulations spanned virtual times of up to 1 × 10⁸, 5 × 10⁵, and 2 × 10⁵ years for the 6- and 10-planet, 126-planet, and 210-planet systems, respectively. Our results show that, for a given class of system (e.g., five pairs of Trojan planets orbiting in the same direction), the relationship between orbit crossing times and planetary spacing is well fit by the functional form log(t _c /t ₀) = b β + c, where t _c is the crossing time, t ₀ = 1 year, β is the separation in initial orbital semimajor axis (in terms of the mutual Hill radii of the planets), and b and c are fitting constants. The same functional form was observed in the previous studies of single planets on nested orbits (Smith and Lissauer 2009). Pairs of Trojan planets are more stable than pairs initially separated by 180°. Systems with retrograde planets (i.e., some planets orbiting in the opposite sense from others) can be packed substantially more closely than can systems with all planets orbiting in the same sense. To have the same characteristic lifetime, systems with 2 or 42 planets per orbit typically need to have about 1.5 or 2 times the orbital separation as orbits occupied by single planets, respectively.     

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 NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

The Apsidal Antialignment of the HD 82943 System

We perform numerical simulations to explore the dynamical evolution of the HD 82943 planetary system. By simulating diverse planetary configurations, we find two mechanisms of stabilizing the system: the 2:1 mean motion resonance (MMR) between the two planets can act as the first mechanism for all stable orbits. The second mechanism is a dynamical antialignment of the apsidal lines of the orbiting planets, which implies that the difference of the periastron longitudes ₃ librates about 180° in the simulations. We also use a semi-analytical model to explain the numerical results for the system under study.     

On the stability of circular ‘asteroid’ orbits in anN-planetary system

The restricted problem of the motion of a point of negligible mass (asteroid) in anN-planetary system is considered. It is assumed that all the planets move about the central body (Sun) along circular orbits in the same plane and the mean motions of the asteroid and the planets are incommensurable. The asteroid orbit evolution is described as a first approximation by secular equations with the perturbing function averaged by the mean longitudes of the asteroid and the planets. For small values of the asteroid orbit eccentricity an expression for the secular part of the perturbing function has been obtained. This expression holds for the arbitrary values of the asteroid orbit semiaxis which are different from those of the planet orbit radii. The stability of the asteroid circular orbits in a linear approximation with respect to the eccentricity is studied. The critical inclinations for a Solar system model are calculated.     

Encounter frequency of Halley-like comets with the planets

We present the results of a numerical study on encounter frequencies of fictitious Halley-like comets with the planets in a dynamical model of the solar system, in which we take into account the gravitational forces of the Sun and the planets Venus through Neptune. The change of the orbital elements during a close approach with a planet was carefully monitored with the aid of a thoroughly tested numerical integration method with automatic step size control. We computed the encounter frequencies of the comets' orbits using two different spheres of influence and compared the results. In both cases, it turned out that the encounter frequency of the fictitious Halley-like comets with Jupiter and Saturn is about a factor 10 to 100 higher than for the other planets. Concerning the changes of the semi-major axes and inclinations our results show that an increase and decrease of these elements is equally probable after an encounter.     

External perturbations on distant planetary orbits and objects in the Kuiper disk

We investigate the potential importance of molecular cloud and stellar perturbations on the orbits of Pluto and more distant (hypothetical) planets up to 500 AU from the Sun. It is found that stellar and molecular cloud-core perturbations are of roughly equal importance. It also is found that the likelihood of substantial perturbations on Pluto is not insignificant, and that numerous substantial stellar and molecular cloud perturbations are likely to have influenced the orbits of any planets beyond 200 AU. These perturbations may contribute to a prevalence of moderate eccentricities and inclinations for planets beyond the orbit of Neptune, and may be a characteristic of distant planetary orbits in other solar systems. Given the recent discovery of chaotic behavior in Pluto's orbit (Sussman and Wisdom 1988), the effects of external perturbations on the long-term stability of Pluto's orbit warrant continued study.     

On the solving of the Hill's differential equation through the method of operational calculus

Using a method previously applied to the treatment of the Mathieu differential equation, we solve the Hill's differential equation of lunar theory through the way of operational calculus, which avoids the cumbersome infinite determinants of the classical procedure. The one-sided Laplace transformation changes it into a difference equation with an infinite number of terms and variable coefficients. When its first member is divided by a suitable factor, this difference equation is the image of an integral equation of the Volterra type which is equivalent to the initial Hill's differential equation. Solution of this Volterra integral equation is unique and it is the general solution of the Hill's differential equation. This solution is a series in the powers of a small dimensionless parameter ² which appears as a factor in the second member of the Hill's differential equation. We reduce it to the sum of its terms of degree 12 with respect to which is the precision usually required in a lunar theory and we write down effectively the coefficients of the terms in ², ( ²)² and the coefficient of the term in ( ²)³ which depends upon the initial valuey(0) of the Hill's differential equation.     

Fractal Dimensions as Chaos Indicators

The correlation dimension D ₂ is used to develop a method of classification for phase space orbits. D ₂ depends only on the mutual distances of the orbits points, therefore the time development of the orbit is reflected in the time development of the correlation dimension approximants. It is shown, that this technique allows to investigate the dynamical properties of a phase space orbit, a classification of chaotic and regular orbits and a detection of sticky orbits.     

Long term evolution of comet Halleys orbit

The aim of the paper is to study the long term evolution of comet Halleys orbit taking into account small errors in the initial conditions. Recent papers deal with mapping methods to model cometary dynamics; (e.g. Petrosky and Broucke, 1987 and Chirikov and Vecheslavov, 1986). They will be discussed critically and compared with our own results. We then tested the model using numerical integration methods. For the moment we limited our calculation to 2.10⁵ years, but a 10⁶ year integration is still in progress. We show the expected dynamical evolution of Hallyes orbit taking into account also smaller and larger errors of the initial conditions (nongravitational effects are only roughly estimated). Finally we discuss alsothe controversal opinions concerning the role of the planets (especially the earth).     

Families of periodic orbits in the general three-body problem for the Sun-Jupiter-Saturn mass-ratio and their stability

Several families of planar planetary-type periodic orbits in the general three-body problem, in a rotating frame of reference, for the Sun-Jupiter-Saturn mass-ratio are found and their stability is studied. It is found that the configuration in which the orbit of the smaller planet is inside the orbit of the larger planet is, in general, more stable.We also develop a method to study the stability of a planar periodic motion with respect to vertical perturbations. Planetary periodic orbits with the orbits of the two planets not close to each other are found to be vertically stable. There are several periodic orbits that are stable in the plane but vertically unstable and vice versa. It is also shown that a vertical critical orbit in the plane can generate a monoparametric family of three-dimensional periodic orbits.     

Numerical analysis of the symmetric methods

Aimed at the initial value problem of the particular second-order ordinary differential equations,y =f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.     

Solar structure in terms of Gauss' hypergeometric function

Hydrostatic equilibrium and energy conservation determine the conditions in the gravitationally stabilized solar fusion reactor. We assume a matter density distribution varying non-linearly through the central region of the Sun. The analytic solutions of the differential equations of mass conservation, hydrostatic equilibrium, and energy conservation, together with the equation of state of the perfect gas and a nuclear energy generation rate = ₀ ⁿ T ⁿT ^m ,are given in terms of Gauss' hypergeometric function. This model for the structure of the Sun gives the run of density, mass, pressure, temperature, and nuclear energy generation through the central region of the Sun. Because of the assumption of a matter density distribution, the conditions of hydrostatic equilibrium and energy conservation are separated from the mode of energy transport in the Sun.     

Synthetic Proper Elements for Outer Main Belt Asteroids

For the orbits with low to moderate inclination and eccentricity, in the asteroid main belt, the analytically computed proper elements are accurate to a level very close to the best result achievable by any analytical theory. This fundamental limitation results from the infinite web of resonances and because of the occurrence of chaotic motions. Still, there are some regions of the belt in which these proper elements are of degraded accuracy, thus preventing a reliable definition of asteroid families and detailed studies of the dynamical structure. We have used a different method to compute asteroid proper elements, following the approach introduced in the LONGSTOP project to describe the secular dynamics of the major outer planets. By applying purely numerical techniques, we produced so-called synthetic proper elements for a catalog of 10,256 asteroids with osculating semimajor axes between 2.5 and 4.0AU.The procedure consisted of simultaneous integration of asteroid and planetary orbits for 2Myr, with online filtering of the short-periodic perturbations. The output of the integration was spectrally resolved, and the principal harmonics (proper values) extracted from the time series. For each asteroid we have also tested the accuracy and stability in time of the proper elements, and estimated the maximum Lyapunov Characteristic Exponent to monitor the chaotic behaviors. This provided information on the reliability of the data for each orbit, in particular allowing to select 1,852 cases for an extended integration (10Myr) of the orbits showing instability. The results indicate that for more than half of the cases the proper elements have a time stability improved by more than a factor 3 with respect to the elements computed by the previous analytical theory. But of course there are also unstable cases for which the proper elements are less accurate and reliable, the extreme examples being 23 orbits exhibiting hyperbolic escape from the solar system. This form of escape from the asteroid belt could be responsible for a significant mass loss over the age of the solar system.     

Minor planet observations and the fundamental reference system

A 15 year project to establish a dynamical reference system utilizing ground-based and Space Telescope observations of 34 minor planets is being undertaken. The orbits of these minor planets will be knit into a common system through the use of crossing point observations. The system of orbits thus established can be used to measure long arcs in the sky (similar to the function of a transit circle) and can be used to detect individual star errors as well as residual periodic effects in the fundamental reference system. The minor planet dynamical reference system will also provide an independent method to establish the zero point and the solid-body rotation of the HIPPARCOS reference system.     

Three-Dimensional p–q Resonant Orbits Close to Second Species Solutions

The purpose of this paper is to study, for small values of , the three-dimensional p–q resonant orbits that are close to periodic second species solutions (SSS) of the restricted three-body problem. The work is based on an analytic study of the in- and out-maps. These maps are associated to follow, under the flow of the problem, initial conditions on a sphere of radius around the small primary, and consider the images of those initial points on the same sphere. The out-map is associated to follow the flow forward in time and the in-map backwards. For both mappings we give analytical expressions in powers of the mass parameter. Once these expressions are obtained, we proceed to the study of the matching equations between both, obtaining initial conditions of orbits that will be 'periodic' with an error of the order ^1–, for some (1/3,1/2). Since, as 0, the inner solution and the outer solution will collide with the small primary, these orbits will be close to SSS.     

Why do Trojan ASCs (not) Escape?

The orbits of 12 Trojan asteroids, which have Lyapunov times T _L10⁵ years and were previously classified as ASCs(=asteroids in stable chaos), are integrated for 50 Myrs, along with a group of neighbouring initial conditions for each nominal orbit. About 40% of the orbits present strong instabilities in the inclination, which may be attributed primarily to the action of the ₁₆ secular resonance; two escapes are also recorded. Higher-order secular resonances, involving the nodes of the outer planets, are also found to be responsible for chaotic motion. Orbital stability depends critically on the choice of initial conditions and, thus, these objects can be regarded as being on the edge of strong chaos.     

Relativistic effects and solar oblateness from radar observations of planets and spacecraft

We used more than 250 000 high-precision American and Russian radar observations of the inner planets and spacecraft obtained in the period 1961–2003 to test the relativistic parameters and to estimate the solar oblateness. Our analysis of the observations was based on the EPM ephemerides of the Institute of Applied Astronomy, Russian Academy of Sciences, constructed by the simultaneous numerical integration of the equations of motion for the nine major planets, the Sun, and the Moon in the post-Newtonian approximation. The gravitational noise introduced by asteroids into the orbits of the inner planets was reduced significantly by including 301 large asteroids and the perturbations from the massive ring of small asteroids in the simultaneous integration of the equations of motion. Since the post-Newtonian parameters and the solar oblateness produce various secular and periodic effects in the orbital elements of all planets, these were estimated from the simultaneous solution: the post-Newtonian parameters are β = 1.0000 ± 0.0001 and γ = 0.9999 ± 0.0002, the gravitational quadrupole moment of the Sun is J₂ = (1.9 ± 0.3) × 10^?7, and the variation of the gravitational constant is ?/G = (?2 ± 5) × 10^?14 yr^?1. The results obtained show a remarkable correspondence of the planetary motions and the propagation of light to General Relativity and narrow significantly the range of possible values for alternative theories of gravitation.     

Searches for transplutonian planets using long-period comets

We discuss the dynamical connection of long-period and nearly parabolic comets with hypothetical transplutonian planets. The statistics includes 792 comets with periods P > 200 years. The orbital plane of the parent planet can be determined from the observed distribution of the perihelia and poles of cometary orbits. The radius of a planetary orbit can be calculated using the Radzievsky-Tisseran criterion. We calculated the minimum distance of each of the 792 orbits to 11 hypothetical planetary orbits. Testing for the kinematic connection of comets with transplutonian planets yielded a negative result. The presence of the nodes of cometary orbits in the transplutonian region is shown to be the result of a geometric effect. We found a high concentration of the nodes and perihelia of cometary orbits in the zone of the terrestrial planets.     

Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some “fast-track” dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.