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.     

Sensitivity study of high eccentricity orbits for Mars gravity recovery

By linear perturbation theory, a sensitivity study is presented to calculate the contribution of the Mars gravity field to the orbital perturbations in velocity for spacecrafts in both low eccentricity Mars orbits and high eccentricity orbits(HEOs). In order to improve the solution of some low degree/order gravity coefficients, a method of choosing an appropriate semimajor axis is often used to calculate an expected orbital resonance, which will significantly amplify the magnitude of the position and velocity perturbations produced by certain gravity coefficients. We can then assess to what degree/order gravity coefficients can be recovered from the tracking data of the spacecraft. However, this existing method can only be applied to a low eccentricity orbit, and is not valid for an HEO. A new approach to choosing an appropriate semimajor axis is proposed here to analyze an orbital resonance. This approach can be applied to both low eccentricity orbits and HEOs. This small adjustment in the semimajor axis can improve the precision of gravity field coefficients and does not affect other scientific objectives.     

On the long-period evolution of the sun-synchronous orbits

The dynamic evolution of sun-synchronous orbits at a time interval of 20 years is considered. The numerical motion simulation has been carried out using the Celestial Mechanics software package developed at the Institute of Astronomy of the University of Bern. The dependence of the dynamic evolution on the initial value of the ascending node longitude is examined for two families of sun-synchronous orbits with altitudes of 751 and 1191 km. Variations of the semimajor axis and orbit inclination are obtained depending on the initial value of the ascending node longitude. Recommendations on the selection of orbits, in which spent sun-synchronous satellites can be moved, are formulated. Minimal changes of elements over a time interval of 20 years have been observed for orbits in which at the initial time the angle between the orbit ascending node and the direction of the Sun measured along the equator have been close to 90° or 270°. In this case, the semimajor axis of the orbit is not experiencing secular perturbations arising from the satellite’s passage through the Earth’s shadow.     

(42355) Typhon-Echidna: Scheduling observations for binary orbit determination

We describe a strategy for scheduling astrometric observations to minimize the number required to determine the mutual orbits of binary transneptunian systems. The method is illustrated by application to Hubble Space Telescope observations of (42355) Typhon-Echidna, revealing that Typhon and Echidna orbit one another with a period of 18.971±0.006 days and a semimajor axis of 1628±29 km, implying a system mass of (9.49±0.52)×¹⁰¹⁷ kg. The eccentricity of the orbit is 0.526±0.015. Combined with a radiometric size determined from Spitzer Space Telescope data and the assumption that Typhon and Echidna both have the same albedo, we estimate that their radii are and , respectively. These numbers give an average bulk density of only , consistent with very low bulk densities recently reported for two other small transneptunian binaries.     

Optimal two-impulse rendezvous using constrained multiple-revolution Lambert solutions

A solution to the fixed-time minimum-fuel two-impulse rendezvous problem for the general non-coplanar elliptical orbits is provided. The optimal transfer orbit is obtained using the constrained multiple-revolution Lambert solution. Constraints consist of lower bound for perigee altitude and upper bound for apogee altitude. The optimal time-free two-impulse transfer problem between two fixed endpoints implies finding the roots of an eighth order polynomial, which is done using a numerical iterative technique. The set of feasible solutions is determined by using the constraints conditions to solve for the short-path and long-path orbits semimajor axis ranges. Then, by comparing the optimal time-free solution with the feasible solutions, the optimal semimajor axis for the two fixed-endpoints transfer is identified. Based on the proposed solution procedure for the optimal two fixed-endpoints transfer, a contour of the minimum cost for different initial and final coasting parameters is obtained. Finally, a numerical optimization algorithm (e.g., evolutionary algorithm) can be used to solve this global minimization problem. A numerical example is provided to show how to apply the proposed technique.     

Bounded orbits in the elliptic restricted three-body problem

Special solutions of the planar rectilinear elliptic restricted 3-body problem are investigated for the limiting case e=1. Numerical integration is performed for primaries of equal masses. Starting values which define circular orbit solutions lead to bounded solutions if the initial radius a₀ is larger than 3.74 in units of the primaries' semimajor axis a. A comparison with the Eulerian two-fixedcentre problem is presented in order to understand qualitatively the characteristic features of bounded orbits and the transition to escape orbits.     

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.     

Averaged model to study long-term dynamics of a probe about Mercury

This paper provides a method for finding initial conditions of frozen orbits for a probe around Mercury. Frozen orbits are those whose orbital elements remain constant on average. Thus, at the same point in each orbit, the satellite always passes at the same altitude. This is very interesting for scientific missions that require close inspection of any celestial body. The orbital dynamics of an artificial satellite about Mercury is governed by the potential attraction of the main body. Besides the Keplerian attraction, we consider the inhomogeneities of the potential of the central body. We include secondary terms of Mercury gravity field from \(J_2\) up to \(J_6\), and the tesseral harmonics \(\overline{C}_{22}\) that is of the same magnitude than zonal \(J_2\). In the case of science missions about Mercury, it is also important to consider third-body perturbation (Sun). Circular restricted three body problem can not be applied to Mercury–Sun system due to its non-negligible orbital eccentricity. Besides the harmonics coefficients of Mercury’s gravitational potential, and the Sun gravitational perturbation, our average model also includes Solar acceleration pressure. This simplified model captures the majority of the dynamics of low and high orbits about Mercury. In order to capture the dominant characteristics of the dynamics, short-period terms of the system are removed applying a double-averaging technique. This algorithm is a two-fold process which firstly averages over the period of the satellite, and secondly averages with respect to the period of the third body. This simplified Hamiltonian model is introduced in the Lagrange Planetary equations. Thus, frozen orbits are characterized by a surface depending on three variables: the orbital semimajor axis, eccentricity and inclination. We find frozen orbits for an average altitude of 400 and 1000 km, which are the predicted values for the BepiColombo mission. Finally, the paper delves into the orbital stability of frozen orbits and the temporal evolution of the eccentricity of these orbits.     

A Mechanism for Producing Short-Period Binaries

When three point masses form a hierarchical triple system, the short-period orbit can be severely modified by the long-period orbit if the two orbits are inclined to each other by more than about 39^deg (). Such an inclination can induce ‘Kozai cycles’ (Kozai, 1962), in which the eccentricity of the inner orbit cycles by a large amount while its period and therefore semimajor axis remains roughly constant. During those periastra when the eccentricity is largest, tidal friction may become important, and this can result in a secular shrinkage of the orbit, until it becomes circularised at a period of a few days.However, apsidal motion due to either GR or to the quadrupolar distortion of the components (if they are no longer treated as point masses) can reduce the range of eccentricity. We explore the limits on outer and inner orbital period that these perturbations imply.If the components are F/G/K/M dwarfs, then rotationally-driven dynamo activity can become important at the short periods that can occur in the right circumstances. It can cause the period to shorten further. The result may be a contact binary, and/or a merger in which the two stars of the inner pair coalesce to form a single rapidly rotating star. We suggest that this may be the origin of AB Dor, a very rapidly rotating K dwarf that is probably about 50 Myr old.     

Sun-synchronous orbits near critical inclination

The long period dynamics of Sun-synchronous orbits near the critical inclination 116.6° are investigated. It is known that, at the critical inclination, the average perigee location is unchanged by Earth oblateness. For certain values of semimajor axis and eccentricity, orbit plane precession caused by Earth oblateness is synchronous with the mean orbital motion of the apparent Sun (a Sun-synchronism). Sun-synchronous orbits have been used extensively in meteorological and remote sensing satellite missions. Gravitational perturbations arising from an aspherical Earth, the Moon, and the Sun cause long period fluctuations in the mean argument of perigee, eccentricity, inclination, and ascending node. Double resonance occurs because slow oscillations in the perigee and Sun-referenced ascending node are coupled through the solar gravity gradient. It is shown that the total number and infinitesimal stability of equilibrium solutions can change abruptly over the Sun-synchronous range of semimajor axis values (1.54 to 1.70 Earth radii). The effect of direct solar radiation pressure upon certain stable equilibria is investigated.     

Ephemeris of a highly eccentric orbit: Explorer 28

An ephemeris has been obtained for Explorer 28 (IMP 3) which agrees well with 2 years of radio observations and with SAO observations a year later. This ephemeris is generated over the 3 year lifetime by a numerical integration method utilizing a set of initial conditions, at launch and without requiring further differential correction. Because highly eccentric orbits are difficult to compute with acceptable accuracy and because a long continuous arc has been obtained which compares with actual data to a known precision, this ephemeris may be used as a standard for computing highly eccentric orbits in the Earth-Moon system.Orbit improvement was used to obtain the initial conditions which generated the ephemeris. This improvement was based on correcting the energy by adjusting the semimajor axis to match computed times of perigee passage with the observed. This procedure may generate errors in semimajor axis to compensate for model errors in the energy; however this compensation error is also implicit in orbit determination itself.     

The Trojan asteroid belt: Proper elements,stability, chaos and families

I have computed proper elements for 174 asteroids in the 1 : 1 resonance with Jupiter, that is for all the reliable orbits available (numbered and multi-opposition). The procedure requires numerical integration, under the perturbations by the four major planets, for 1,000,000 years; the output is digitally filtered and compressed into a synthetic theory (as defined within theLONGSTOP project). The proper modes of oscillation of the variables related to eccentricity, perihelion, inclination and node define proper elements. A third proper element is defined as the amplitude of the oscillation of the semimajor axis associated with the libration period; because of the strong nonlinearity of the problem, this component cannot be determined by a simple Fourier transform to the frequency domain. I therefore give another definition, which results in very good stability with time. For 87% of the computed orbits, the stability of the proper elements-at least over 1M yr-is within the following bounds: 0.001AU in semimajor axis, 0.0025 in eccentricity and sine of inclination. Half of the cases with degraded stability of the proper elements are found to be chaotic, with e-folding times between 16,000 and 660,000yr; in some other cases, chaotic behaviour does not result in a significantly decreased stability of the proper elements (stable chaos). The accuracy and stability of these proper elements is good enough to allow a search for asteroid families; however, the dynamical structure of the Trojan belt is very different from the one of the main belt, and collisional events among Trojans can result in a distribution of fragments difficult to identify. The occurrence of couples of Trojans with very close proper elements is proven not to be statistically significant in almost all cases. As the only exception, the couple 1583 Antilochus — 3801 Thrasimedes is significant; however, it is not easy to account for it by a conventional collisional theory. The Menelaus group is confirmed as a strong candidate collisional family; Teucer and Sarpedon could be considered as significant clusters. A number of other clumps are detected (by the same automated clustering method used for the main belt by Zappalà et al., 1990, 1992), but the total number of Trojans with reliable orbits is not large enough to detect many significant candidate families.     

The Poynting-Robertson effect in the Newtonian potential with a Yukawa correction

We consider a Yukawa-type gravitational potential combined with the Poynting-Robertson effect. Dust particles originating within the asteroid belt and moving on circular and elliptic trajectories are studied and expressions for the time rate of change of their orbital radii and semimajor axes, respectively, are obtained. These expressions are written in terms of basic particle parameters, namely their density and diameter. Then, they are applied to produce expressions for the time required by the dust particles to reach the orbit of Earth. For the Yukawa gravitational potential, dust particles of diameter \(10^{ - 3}\) m in circular orbits require times of the order of \(8.557 \times 10^{6}\) yr and for elliptic orbits of eccentricities \(e =0.1, 0.5\) require times of \(9.396 \times 10^{6}\) and \(2.129 \times 10^{6}\) yr respectively to reach Earth’s orbit. Finally, various cases of the Yukawa potential are studied and the corresponding particle times to reach Earth’s are derived per case along with numerical results for circular and various elliptical orbits.     

Theory of satellite orbit-orbit resonance

On the basis of the strong mathematical and physical parallels between orbit-orbit and spin-orbit resonances, the dynamics of mutual orbit perturbations between two satellites about a massive planet are examined, exploiting an approach previously adopted in the study of spin-orbit coupling. The satellites are assumed to have arbitrary mass ratio and to move in non-intersecting orbits of arbitrary size and eccentricity. Resonances are found to exist when the mean orbital periods are commensurable with respect to some rotating axis, which condition also involves the apsidal and nodal motions of both satellites. In any resonant state the satellites are effectively trapped in separate potential wells, and a single variable is found to describe the simultaneous librations of both satellites. The librations in longitude are 180° out-of-phase, with fixed amplitude ratio that depends only on their relative masses and semimajor axes. At the same time the stroboscopic longitude of conjunction also librates about the commensurate axis with the same period. The theory is applicable to Saturn's resonant pairs Titan-Hyperion and Mimas-Tethys, and in these cases our calculated libration periods are in reasonably good agreement with the observed periods.This research supported under a grant from the California Institute of Technology President's Fund and NASA Contract NAS 7-100.     

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.     

Where are the Uranus Trojans?

The area of stable motion for fictitious Trojan asteroids around Uranus’ equilateral equilibrium points is investigated with respect to the inclination of the asteroid’s orbit to determine the size of the regions and their shape. For this task we used the results of extensive numerical integrations of orbits for a grid of initial conditions around the points L ₄ and L ₅, and analyzed the stability of the individual orbits. Our basic dynamical model was the Outer Solar System (Jupiter, Saturn, Uranus and Neptune). We integrated the equations of motion of fictitious Trojans in the vicinity of the stable equilibrium points for selected orbits up to the age of the Solar system of 5 × 10⁹ years. One experiment has been undertaken for cuts through the Lagrange points for fixed values of the inclinations, while the semimajor axes were varied. The extension of the stable region with respect to the initial semimajor axis lies between 19.05 ≤ a ≤ 19.3 AU but depends on the initial inclination. In another run the inclination of the asteroids’ orbit was varied in the range 0° < i < 60° and the semimajor axes were fixed. It turned out that only four ‘windows’ of stable orbits survive: these are the orbits for the initial inclinations 0° < i < 7°, 9° < i < 13°, 31° < i < 36° and 38° < i < 50°. We postulate the existence of at least some Trojans around the Uranus Lagrange points for the stability window at small and also high inclinations.     

Tidal evolution of close binary asteroid systems

We provide a generalized discussion of tidal evolution to arbitrary order in the expansion of the gravitational potential between two spherical bodies of any mass ratio. To accurately reproduce the tidal evolution of a system at separations less than 5 times the radius of the larger primary component, the tidal potential due to the presence of a smaller secondary component is expanded in terms of Legendre polynomials to arbitrary order rather than truncated at leading order as is typically done in studies of well-separated system like the Earth and Moon. The equations of tidal evolution including tidal torques, the changes in spin rates of the components, and the change in semimajor axis (orbital separation) are then derived for binary asteroid systems with circular and equatorial mutual orbits. Accounting for higher-order terms in the tidal potential serves to speed up the tidal evolution of the system leading to underestimates in the time rates of change of the spin rates, semimajor axis, and mean motion in the mutual orbit if such corrections are ignored. Special attention is given to the effect of close orbits on the calculation of material properties of the components, in terms of the rigidity and tidal dissipation function, based on the tidal evolution of the system. It is found that accurate determinations of the physical parameters of the system, e.g., densities, sizes, and current separation, are typically more important than accounting for higher-order terms in the potential when calculating material properties. In the scope of the long-term tidal evolution of the semimajor axis and the component spin rates, correcting for close orbits is a small effect, but for an instantaneous rate of change in spin rate, semimajor axis, or mean motion, the close-orbit correction can be on the order of tens of percent. This work has possible implications for the determination of the Roche limit and for spin-state alteration during close flybys.     

There are ten,not eight planets

In 1946, E. Sevin postulated the global vibrations of the Sun with a period P ₀ = 1/9 day and a “wavelength” L ₀ = c × P ₀ = 19.24 AU and predicted the tenth planet at a mean distance of 4.0 × L ₀ ≈ 77.0 AU from the Sun (c is the speed of light). The global vibrations of the Sun, precisely with the period of 1/9 day, were actually detected in 1974. Recently, the largest Kuiper Bell object 2003 UB₃₁₃, or Eris, with an orbital semimajor axis ≈ 3.5 × L ₀ ≈ 67.5 AU was discovered. We adduce arguments for the status of Eris as our tenth planet: (i) the object is larger and farther from the Sun than Pluto and (ii) the semimajor axis of Eris agrees well with the sequence of planetary distances that follows from the resonance spectrum of the Solar system dimensions (with the scale L ₀ and for all 11 orbits, including those of Pluto, Eris, and the asteroid belt). We point to a mistake of the Prague (2006) IAU Assembly, which excluded Pluto from the family of planets by introducing a new, highly controversial class of objects—“dwarf planets.”     

The dynamics of tadpole and horseshoe orbits: II. The coorbital satellites of saturn

The recently discovered coorbital satellites of Saturn, 1980S1 and 1980S3, are shown to be librating in horseshoe orbits. By considering the effects of tangential forces on the semimajor axes of the satellite orbits, we derive an accurate relation between the sum of the satellite masses and (a) their minimum angular separation, (b) the variation of their angular separation with time and (c) the libration period. Observations of (b) and (c) are the most practical methods of determining the satellite masses. The orbits of the coorbital satellites of Dione and Tethys are discussed. We demonstrate the possibility of calculating a new value for the mass of Dione and we show that one of the coorbital satellites of Tethys could be moving in a horseshoe orbit even though another satellite is librating in a tadpole orbit about the leading Lagrangian equilibrium point L₄. The origin of coorbital satellites and the stability of their orbits are discussed.     

On the Instability of Jupiter's Trojans

We have numerically explored the mechanisms that destabilize Jupiter's Trojan orbits outside the stability region defined by Levison et al. (1997, Nature385, 42-44). Different models have been exploited to test various possible sources of instability on timescales on the order of ∼10⁸ years.In the restricted three-body model, only a few Trojan orbits become unstable within 10⁸ years. This intrinsic instability contributes only marginally to the overall instability found by Levison et al.In a model where the orbital parameters of both Jupiter and Saturn are fixed, we have investigated the role of Saturn and its gravitational influence. We find that a large fraction of Trojan orbits become unstable because of the direct nonresonant perturbations by Saturn. By shifting its semimajor axis at constant intervals around its present value we find that the near 5:2 mean motion resonance between the two giant planets (the Great Inequality) is not responsible for the gross instability of Jupiter's Trojans since short-term perturbations by Saturn destabilize Trojans, even when the two planets are far out of the resonance.Secular resonances are an additional source of instability. In the full six-body model with the four major planets included in the numerical integration, we have analyzed the effects of secular resonances with the node of the planets. Trojan asteroids have relevant inclinations, and nodal secular resonances play an important role. When a Trojan orbit becomes unstable, in most cases the libration amplitude of the critical argument of the 1:1 mean motion resonance grows until the asteroid encounters the planet. Libration amplitude, eccentricity, and nodal rate are linked for Trojan orbits by an algebraic relation so that when one of the three parameters is perturbed, the other two are affected as well. There are numerous secular resonances with the nodal rate of Jupiter that fall inside the region of instability and contribute to destabilize Trojans, in particular the ν₁₆. Indeed, in the full model the escape rate over 50 Myr is higher compared to the fixed model.Some secular resonances even cross the stability region delimited by Levison et al. and cause instability. This is the case of the 3:2 and 1:2 nodal resonances with Jupiter. In particular the 1:2 is responsible for the instability of some clones of the L4 Trojan (3540) Protesilaos.