A preliminary analysis of the orbit of the Mars Trojan asteroid (5261) Eureka

Observations and results of orbit determination of the first known Mars Trojan asteroid (5261) Eureka are presented. We have numerically calculated the evolution of the orbital elements, and have analyzed the behavior of the motion during the next 2 Myr. Strong perturbations by planets other than Mars seem to stabilize the eccentricity of the asteroid by stirring the high order resonances present in the elliptic restricted problem. As a result, the orbit appears stable at least on megayear timescales. The difference of the mean longitudes of Mars and Eureka and the semimajor axis of the asteroid form a pair of variables that essentially behave in an adiabatic manner, while the evolution of the other orbital elements is largely determined by the perturbations due to other planets.     

Planetary theories with the aid of the expansions of elliptic functions

For coplanar circular orbits, the mutual perturbations between two bodies can be expressed in term of the argument of Jacobian elliptic functions instead of the difference of the mean longitudes. For a given pair of planets, such a change of time variable improves the convergence of the developments. At the first order of planetary masses an integration of Lagrange's equations for the osculating elements is performed. When compared to classical developments the results are reduced by an important factor. The method is then extended to the mutual perturbations of Jupiter and Saturn, at any order of planetary masses, either with Fourier series with two arguments, or with one argument solely, taking advantage of the close commensurability of the mean motions.     

Numerical Computation of Hansen-like Expansions

Fourier expansions of elliptic motion functions in multiples of the true, eccentric, elliptic and mean anomalies are computed numerically by means of the fast Fourier transform. Both Hansen-like coefficients and their derivatives with respect to eccentricity of the orbit are considered. General behavior of the coefficients and the efficiency (compactness) of the expansions are investigated for various values of eccentricity of the orbit. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

Generation of a secular system in the analytical theory for the motion of the moon

In order to generate an analytical theory of the motion of the Moon by considering planetary perturbations, a procedure of general planetary theory (GPT) is used. In this case, the Moon is considered as an addition planet to the eight principal planets. Therefore, according to the GPT procedure, the theory of the Moon’s orbital motion can be presented in the form of series with respect to the evolution of eccentric and oblique variables with quasi-periodic coefficients, which are the functions of mean longitudes for principal planets and the Moon. The relationship between evolution variables and the time is determined by a trigonometric solution for the independent secular system that describes the secular motion of a perigee and the Moon node by considering secular planetary inequalities. Principal planetary coordinates required for generating the theory of the motion of the Moon includes only Keplerian terms, the intermediate orbit, and the linear theory with respect to eccentricities and inclinations in the first order relative to the masses. All analytical calculations are performed by means of the specialized echeloned Poisson Series Processor EPSP.     

On the attitude dynamics of perturbed triaxial rigid bodies

Attitude dynamics of perturbed triaxial rigid bodies is a rather involved problem, due to the presence of elliptic functions even in the Euler equations for the free rotation of a triaxial rigid body. With the solution of the Euler–Poinsot problem, that will be taken as the unperturbed part, we expand the perturbation in Fourier series, which coefficients are rational functions of the Jacobian nome. These series converge very fast, and thus, with only few terms a good approximation is obtained. Once the expansion is performed, it is possible to apply to it a Lie-transformation. An application to a tri-axial rigid body moving in a Keplerian orbit is made.     

Artificial frozen orbits around Mercury

Orbits around Mercury are influenced by the strong elliptic third-body perturbation, especially for high eccentricity orbits, the periapsis altitude changes dramatically. Frozen orbits whose mean eccentricity and argument of perigee remain constants are obviously a good choice for space missions, but the forming conditions are too harsh to meet practical needs. To deal with this problem, a continuous control method that combines analytical theory and parameter optimization is proposed to build an artificial frozen orbit. The artificial frozen orbits are investigated on the basis of double averaged Hamiltonian, of which the second and third zonal harmonics and the perturbation of elliptic third-body gravity are considered. In this paper, coefficients of perturbations which satisfy the conditions of frozen orbits are involved as control parameters, and the relevant artificial perturbations are compensated by the control strategy. So probes around Mercury can be kept on frozen orbit under the influence of continuous control force. Then complex method of optimization is used to search for the energy optimized artificial frozen orbits. The choosing of optimal parameters, the objective function setting and other issues are also discussed in the study. Evolution of optimal control parameters are given in large ranges of semi-major axis and eccentricity, through the variation of these curves, the fuel efficiency is discussed. The result shows that the control method proposed in this paper can effectively maintain the eccentricity and argument of perigee frozen.     

A dynamical investigation of the conjecture that Mercury is an escaped satellite of Venus

The possibility that Mercury might once have been satellite of a Venus, suggested by a number of anomalies, is investigated by a series of numerical computer experiments. Tidal interaction between Mercury and Venus would result in the escape of Mercury into a solar orbit. Only two escape orbits are possible, one exterior and one interior to the Venus orbit. For the interior orbit, subsequent encounters are sufficiently distant to avoid recapture or large perturbations. The perihelion distance of Mercury tends to decrease, while the orientation of perihelion librates for the first few thousand revolutions. If dynamical evolution or nonconservative forces were large enough in the early solar system, the present semimajor axes could have resulted. The theoretical minimum quadrupole moment of the inclined rotating Sun would rotate the orbital planes out of coplanarity. Secular perturbations by the other planets would evolve the eccentricity and inclination of Mercury's orbit through a range of possible configurations, including the present orbit. Thus the conjecture that Mercury is an escaped satellite of Venus remains viable, and is rendered more attractive by our failure to disprove it dynamically.     

Analytical direct planetary perturbations on the orbital motion of satellites

An analytical expansion of the disturbing function arising from direct planetary perturbations on the motion of satellites is derived. As a Fourier series, it allows the investigation of the secular effects of these direct perturbations, as well as of every argument present in the perturbation. In particular, we construct an analytical model describing the evection resonance between the longitude of pericenter of the satellite orbit and the longitude of a planet, and study briefly its dynamic. The expansion developed in this paper is valid in the case of planar and circular planetary orbits, but not limited in eccentricity or inclination of the satellite orbit.     

Un formulaire pour le calcul des perturbations d'ordres élevés dans les problèmes planétaires

The authors present formulas in compact form for constructing high order planetary perturbations with respect to the disturbing masses. They have been built by an iterative process and give the variations of osculating elements. Singularities due to vanishing eccentricities and inclinations are not present in the differential equations. All elementary operations are manipulations of Fourier series with numerical coefficients, and great care has been taken to economize algebraic operations. Results are presented in three forms:

1. vectorial form, with real components which may be useful in numerical integrations;
2. complex form, to put in evidence the symmetries of the system of variables;
3. scalar form, which is the most elaborate. This last form has been used for constructing the first order perturbations for any pair of planets. Two illustrations are given (Jupiter and Saturn, Venus and Earth). Further remarks are made about the practical manipulation of Fourier series, resolution of Kepler's equation in complex form and construction by iteration of the inverse of the distance between two bodies.

     

A resonance problem of two degrees of freedom

An analytical theory is presented for determining the motion described by a Hamiltonian of two degrees of freedom. Hamiltonians of this type are representative of the problem of an artificial Earth satellite in a near-circular orbit or a near-equatorial orbit and in resonance with a longitudinal dependent part of the geopotential. Using the classical Bohlin-von Zeipel procedure the variation of the elements is developed through a generating function expressed as a trigonometrical series. The coefficients of this series, determined in ascending powers of an auxiliary parameter, are the solutions of paired sets of ordinary differential equations and involve elliptic functions and quadrature. The first order solution accounts for the full variation of the resonance terms with the second coordinate.     

Stability of Co-Orbital Motion in Exoplanetary Systems

The stability of co-orbital motions is investigated in such exoplanetary systems, where the only known giant planet either moves fully in the habitable zone, or leaves it for some part of its orbit. If the regions around the triangular Lagrangian points are stable, they are possible places for smaller Trojan-like planets. We have determined the nonlinear stability regions around the Lagrangian point L₄ of nine exoplanetary systems in the model of the elliptic restricted three-body problem by using the method of the relative Lyapunov indicators. According to our results, all systems could possess small Trojan-like planets. Several features of the stability regions are also discussed. Finally, the size of the stability region around L₄ in the elliptic restricted three-body problem is determined as a function of the mass parameter and eccentricity.     

Numerical models for the study of motion of lunar satellites

The present study deals with numerical modeling of the elliptic restricted three-body problem as well as of the perturbed elliptic restricted three-body (Earth-Moon-Satellite) problem by a fourth body (Sun). Two numerical algorithms are established and investigated. The first is based on the method of the series solution of the differential equations and the second is based on a 5th-order Runge-Kutta method. The applications concern the solution of the equations and integrals of motion of the circular and elliptical restricted three-body problem as well as the search for periodic orbits of the natural satellites of the Moon in the Earth-Moon system in both cases in which the Moon describes circular or elliptical orbit around the Earth before the perturbations induced by the Sun. After the introduction of the perturbations in the Earth-Moon-Satellite system the motions of the Moon and the Satellite are studied with the same initial conditions which give periodic orbits for the unperturbed elliptic problem.     

Small perturbations in flat galaxies

The aim of this series of papers is to develop straightforward methods of computing the response of flat galaxies to small perturbations. This Paper I considers steady state problems; Paper II considers time varying perturbations and the effects of resonances; and Paper III applies the methods developed in Papers I and II to a numerical study of the stability of flat galaxies.The general approach is to study the dynamics of each individual orbit. The orbits are described by their apocentric and pericentric radii,r _a andr _p , and the distribution function of an equilibrium model is a function ofr _a andr _p . The mass density and potential corresponding to a distribution function is found by means of an expansion in Hankel-Laguerre functions; the coefficients of the expansion being found by taking moments of the mass density of the individual orbits. This leads to a simple method of constructing equilibrium models.The response to a small perturbation is found by seeking the response of each orbit. When the perturbations are axisymmetric and slowly varying, the response can be easily found using adiabatic invariants. The potential is expanded in a series of Hankel-Laguerre functions, and the response operator becomes a discrete matrix. The condition that the model is stable against adiabatic radial perturbations is that the largest eigenvalue of the response matrix should be less than one.An analytic approximation to the response matrix is derived, and applied to estimate the eccentricity needed for stability against local perturbations.     

Model orbits of asteroid 3040 Kozai

The orbital evolutions of the asteroid 3040 Kozai and model asteroids with similar orbits have been investigated. Their osculating orbits for an epoch 1991 December 10 were numerically integrated forward within the interval of 20,000 years, using a dynamical model of the solar system consisting of all inner planets, Jupiter, and Saturn.The orbit of the asteroid Kozai is stable. Its motion is affected only by long-period perturbations of planets. With change of the argument of perihelion of the asteroid Kozai, the evolution of the model asteroid orbits changes essentially, too. The model orbits with the argument of perihelion changed by the order of 10% show that asteroids with such orbital parameters may approach the Earth orbit, while asteroids with larger changes may even cross it, at least after 10,000 years. Long-term orbital evolution of asteroids with these orbital parameters is very sensitive on their angular elements.     

The influence of the great inequality on the secular disturbing function of the planetary system

This paper derives the contributionF ₂ ^* by the great inequality to the secular disturbing function of the principal planets. Andoyer's expansion of the planetary disturbing function and von Zeipel's method of eliminating the periodic terms is employed; thereby, the corrected secular disturbing function for the planetary system is derived. An earlier solution suggested by Hill is based on Leverrier's equations for the variation of elements of Jupiter and Saturn and on the semi-empirical adjustment of the coefficients in the secular disturbing function. Nowadays there are several modern methods of eliminating periodic terms from the Hamiltonian and deriving a purely secular disturbing function. Von Zeipel's method is especially suitable. The conclusion is drawn that the canonicity of the equations for the secular variation of the heliocentric elements can be preserved if there be retained, in the secular disturbing function, terms only of the second and fourth order relative to the eccentricity and inclinations.The Krylov-Bogolubov method is suggested for eliminating periodic terms, if it is desired to include the secular perturbations of the fifth and higher order in the heliocentric elements. The additional part of the secular disturbing functionF ₂ ^* derived in this paper can be included in existing theories of the secular effects of principal planets. A better approach would be to preserve the homogeneity of the theory and rederive all the secular perturbations of principal planets using Andoyer's symbolism, including the part produced by the great inequality.     

Explicit Semianalytical Theory of Asteroid Motion

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

Calculation of the intermediate perturbed orbit from three or more positions of a small body on the celestial sphere

Two new methods are described for finding the orbit of a small celestial body from three or more pairs of angular measurements and the corresponding time points. The methods are based on, first, the approach that has been developed previously by the author to the determination, from a minimum number of observations, of intermediate orbit considering most of the perturbations in the bodies’ motion and, second, Herget’s algorithmic procedure enabling the introduction of additional observations. The errors of orbital parameters calculated by the proposed methods are two orders of magnitude smaller than the corresponding errors of the traditional approach based on the construction of an unperturbed Keplerian orbit. The thus-calculated orbits of the minor planets 1566 Icarus, 2002 EC1, and 2010 TO48 are used to compare the results of Herget’s multiposition procedure and the new methods. The comparison shows that the new methods are highly effective in the study of perturbed motion. They are particularly beneficial if high-precision observational data covering short orbital arcs are available.     

On the perturbations of a close-Earth satellite due to lunar inequalities

The lunar disturbing function for a close-Earth satellite is expressed as a sum of products of harmonics of the satellite's position and harmonics of the Moon's position, and the latter are expanded about a rotating and precessing elliptic orbit inclined to the ecliptic. The deviations of the Moon from this approximate orbit are computed from Brown's lunar theory andthe perturbations in satellite orbital elements due to these inequalities are derived. Numerical calculations indicate that several perturbations in the position of the satellite's node and perigee have magnitudes on the order of one meter.The author is supported in part by a National Science Foundation Graduate Fellowship.     

Orbital structure of the Kuiper belt and its contribution in the near-Earth asteroid and meteor population

Numerical calculations are given to describe evolution of orbits of simulated and real Kuiper belt objects for large intervals of time. Gravitational perturbations caused by all major planets have been taken into account, and, when considering small particles, Poynting-Robertson nongravitational effect has also been incorporated. Large orbit scattering of the Kuiper belt objects regarding the semimajor axes and eccentricities is shown to be due to their evolution over millions of years. Relative contribution of great objects and meteor particles from the Kuiper belt into the near-Earth population is believed to be extremely small.