Computing secular motion under slowly rotating quadratic perturbation

We consider secular perturbations of nearly Keplerian two-body motion under a perturbing potential that can be approximated to sufficient accuracy by expanding it to second order in the coordinates. After averaging over time to obtain the secular Hamiltonian, we use angular momentum and eccentricity vectors as elements. The method of variation of constants then leads to a set of equations of motion that are simple and regular, thus allowing efficient numerical integration. Some possible applications are briefly described.     

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.     

Dynamics of close binary systems with mass exchange

The effects of non-isotropic ejection of mass from either component of a binary system on the orbital elements are studied, for the case of a small initial eccentricity of the relative orbit, when all the ejected mass falls on the other component. The problem is transformed to an equivalent two-body problem with isotropic variation of mass, plus a perturbing force which is a function of the intial conditions of ejection of the particles and their final, positions and velocities when they fall on the surface of the other star. The variation of the orbital elements are derived. It is shown that, to first-order terms in the eccentricity, the secular change of the semimajor axis is equal to the one corresponding to the case of zero initial eccentricity. On the contrary, the secular change of the eccentricity is smaller and it depends on the variations of mass ejection due to the finite eccentricity.     

On the 3/1 planar and circular resonant problem

The simplest model of a resonant problem of second order is the planar and circular case. Simplification like this is very old and for 3/1 resonance, several authors have studied this problem with different purposes. In this work, we test this model for the available asteroids, by applying Hori's perturbation method. Explicit solutions of the intermediate orbit are obtained. In the plane of two constants of the problem, all types of motion are described. By testing the model, it is shown that, in general, one can confirm results of numerical integrations indicating libration for a few number of asteroids and circulation for most of them. However, agreement in numerical values for amplitude and period of librations seems to be not possible mainly if Jupiter's eccentricity is neglected. On the other hand, even though there might be some physical reasons determining that only asteroids with high eccentricity may librate, it is shown that, from mathematical point of view, libration may occur even in the case of small eccentricities provided that some relations are satisfied.     

Gravitational Radiation Damping and Evolution of the Orbit of Compact Binary Stars (Solution by the Second Perturbation Method)

The influence of the gravitational radiation damping on the evolution of the orbital elements of compact binary stars is examined by using the method of perturbation. The perturbation equations with the true anomaly as an independent variable are given. This effect results in both the secular and periodic variation of the semi-major axis, the eccentricity, the mean longitude at the epoch and the mean longitude. However, the longitude of periastron exhibits no secular variation, but only periodic variation. The effect of secular variation of the orbit would lead to collapse of the system of binary stars. The deduced formulae are applied to the calculation of secular variation of the orbital elements for three compact binary stars: PSR19 13 + 16, PSR J0737-3039 and M33X-7. The results obtained are discussed.     

Tidal Evolution of the Kepler Candidate Two-planet Systems

The tidal evolution of ten Kepler candidate two-planet systems is investigated by using the general secular perturbation theory, and then a general picture of tidal evolution for these systems is described. Taking the KOI 1239 system as an example, the tidal effect of the system is studied in detail, the results indicate that the dissipative term of the tidal effect causes the attenuation of planets’ orbital eccentricities, and it plays a dominant role in the process of orbital evolution, however, the conservative term of the tidal effect and the relativistic effect may damp the excitation of the eccentricity of the inner planet under the secular perturbation of the outer planet. In addition, the process of tidal evolution is also affected by both the initial eccentricity of the outer planet and the planet's tidal dissipation coeffcient. At the same time, the numerical simulation on the tidal evolution of the KOI 1239 system is also made, and the numerical results are consistent with those of general secular perturbation theory.     

在电离中电感应阻力对带电卫星轨道根数的摄动影响

研究了在高空电离层中运动的带电荷的卫星受电感应阻力后对轨道根数产生的摄动影响。研究结果表明 ,电感应阻力对带电卫星的轨道半长轴、轨道偏心率、近地点赤经、历元平赤经均有周期摄动影响 ,但除对半长轴有长期摄动效应外对其它轨道根数均无长期摄动。轨道倾角和升交点赤经不受摄动影响。文中以飞行在高度 1 50 0km的电离层中的导体卫星作为算例。计算结果显示 :带电导体卫星在高空电离层中带有一定电量时电感应阻力对轨道半长轴的缩短产生显著效应     

太阳多方模型对行星轨道要素变化的长期摄动影响

本文利用天体力学中的摄动理论和天体物理学中的气体星多方模型理论研究了太阳多方模型对行星轨道要素变化的长期摄动影响。文中给出了太阳日多方指数ｎ＝３的模型由于自转、扁度和内部密度分布等因素对行星轨道要素变化的长期摄动效应的理论结果。研究结果表明：行星轨道要素除长轴、偏心率和轨道倾角不受长期摄动外升交点经度、近日点经度以及平近点角均受长期摄动的影响。最后利用理论结果对行星轨道要素的长期投动效应做了数值计算，数值结果在表１中给出。     

Secular dynamics of the three-body problem: application to the υ Andromedae planetary system

The discovery of extra-solar planetary systems with multiple planets in highly eccentric orbits (∼0.1-0.6), in contrast with our own Solar System, makes classical secular perturbation analysis very limited. In this paper, we use a semi-numerical approach to study the secular behavior of a system composed of a central star and two massive planets in co-planar orbits. We show that the secular dynamics of this system can be described using only two parameters, the ratios of the semi-major axes and the planetary masses. The main dynamical features of the system are presented in geometrical pictures that allows us to investigate a large domain of the phase space of this three-body problem without time-expensive numerical integrations of the equations of motion, and without any restriction on the magnitude of the planetary eccentricities. The topology of the phase space is also investigated in detail by means of spectral map techniques, which allow us to detect the separatrix of a non-linear secular apsidal resonance. Finally, the qualitative study is supplemented by direct numerical integrations. Three different regimes of secular motion with respect to the secular angle Δ? are possible: they are circulation, oscillation (around 0° and 180°), and high eccentricity libration in a non-linear secular resonance. The first two regimes are a continuous extension of the classical linear secular perturbation theory; the last is a new feature, hitherto unknown, in the secular dynamics of the three-body problem. We apply the analysis to the case of the two outer planets in the υ Andromedae system, and obtain its periodic and ordinary orbits, the general structure of its secular phase space, and the boundaries of its secular stability; we find that this system is secularly stable over a large domain of eccentricities. Applying this analysis to a wide range of planetary mass and semi-major axis ratios (centered about the υ Andromedae parameters), we find that apsidal oscillation dominates the secular phase space of the three-body coplanar system, and that the non-linear secular resonance is also a common feature.     

Lunisolar Apsidal Resonances at low Satellite Orbits

Eccentricity resonances between the secular motion of an Earth satellite's orbit and the longitudes of the Sun and the Moon are studied within a Hamiltonian framework. The problem is approximated in a traditional manner, with the Earth's potential including only the second zonal harmonic, and a Hill‐type approximation for perturbing bodies. For a family of 10 resonances, stable and unstable points are identified and libration widths are estimated. Numerical values are given for the maximum variation of eccentricity available at each resonance. The respective amplitudes of the perigee heights' librations range from 2 to 750 km. The resonances of the solar origin are generally stronger than their lunar counterparts. This revised version was published online in July 2006 with corrections to the Cover Date.     

Numerical Instruments for the Analysis of Secular Dynamics of Exoplanetary Systems

A review is given of modern numerical methods for the analysis of resonant and chaotic dynamics: calculation of the Lyapunov characteristic exponents, the MEGNO method, and the maximum eccentricity method. These methods are used to construct stability diagrams for the planetary systems γ Cep, HD 196885, and HD 41004. The diagrams are analyzed to determine the most probable values taken by the orbital parameters of the exoplanets and obtain estimates for the Lyapunov time of their orbital dynamics. The stability diagrams constructed using the different methods are compared to analyze their effectiveness in the study of secular dynamics of exoplanetary systems.     

On the artificial satellite theory

Some of the basic ideas of an analytical orbiter theory which is being developed by Hubert Claes in Namur are presented.The theory is based on the Lie transform technique and will be expressed in a closed form up to second order. The inclusion of additional terms of the third order (expanded in power series of the eccentricity) will be considered.Special attention is being given to the choice of the elements and to the final form of the theory. Three main criteria are used. The removal of the virtual singularities of small inclination and eccentricity. The simplicity of the final form of the theory once the elements have been given their numerical values. The numerical stability of the evaluation of the theory.     

Improved equations for eccentricity generation in hierarchical triple systems

In a series of papers, we developed a technique for estimating the inner eccentricity in hierarchical triple systems, with the inner orbit being initially circular. However, for certain combinations of the masses and the orbital elements, the secular part of the solution failed. In this paper, we derive a new solution for the secular part of the inner eccentricity, which corrects the previous weakness. The derivation applies to hierarchical triple systems with coplanar and initially circular orbits. The new formula is tested numerically by integrating the full equations of motion for systems with mass ratios from 10⁻³ to 10³. We also present more numerical results for short-term eccentricity evolution, in order to get a better picture of the behaviour of the inner eccentricity.     

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.     

A numerical investigation of secular terms of the planetary disturbing function

The secular terms of the planetary disturbing function are given, after elimination of short period terms by von Zeipel's transformation. The adequacy of this expansion up to terms of eighth order in the inclination and eccentricity is investigated by numerical processes, as a function of the Keplerian elementsa, e andi. The eccentricityé of the outer planet, is taken equal to zero. It is concluded that for values ofi which are not small the inclusion of additional terms in the expression for the disturbing function, results to drastic changes in its values, while larger values ofe do not have an equaly large effect on the disturbing function.     

Long period variations of the motion of a satellite due to non-resonant tesseral harmonics of a gravity potential

The main effects of tesseral harmonics of a gravity potential expansion on the motion of a satellite, are short period variations as well as long period variations due to resonances. However, other smaller long period and secular variations can arise from interactions between tesseral terms of the same order. The analytical integration of these effects is developed, using numerical evaluation of Kaula eccentricity and inclination functions. Examples for some Earth's geodetic satellites show that secular effects can reach a few decameters per year. The secular variations can even reach several hundred of meters per year for the Mars natural satellite Phobos.     

A second order analytical atmospheric drag theory based on the TD88 thermospheric density model

A second order atmospheric drag theory based on the usage of TD88 model is constructed. It is developed to the second order in terms of TD88 small parameters K _n,j . The short periodic perturbations, of all orbital elements, are evaluated. The secular perturbations of the semi-major axis and of the eccentricity are obtained. The theory is applied to determine the lifetime of the satellites ROHINI (1980 62A), and to predict the lifetime of the microsatellite MIMOSA. The secular perturbations of the nodal longitude and of the argument of perigee due to the Earth’s gravity are taken into account up to the second order in Earth’s oblateness.     

A simple analytical model for the secular resonance ν6 in the asteroidal belt

When asteroids are in the secular resonance ₆, the variation of the eccentricity becomes very large. In this paper, the dynamics of this secular resonance ₆ is investigated by a simple analytical model, in which the third degree terms of the eccentricity and inclination are taken into account. The eccentricity variations of asteroids located near this resonance are represented clearly by the diagrams of equi-Hamiltonian curves on the plane of versuse ( the longitude of perihelion of asteroids and Saturn,e: the eccentricity of asteroids). These diagrams predict that the eccentricity of these asteroids suffers a large increase or decrease, and that the secular resonance argument librates about 0° and 180°. In order to confirm these predictions, numerical integrations are carried out over one million years. By these integrations, it is found that the eccentricity of secular resonant asteroids becomes more than 0.8, and that the libration about 0° also exists, as well as the libration about 180°. The strongly depopulated region in the asteroidal belt, which corresponds to the position of the secular resonance ₆, is also explained well by this analytical model.     

Parameterized post-Newtonian orbital effects in extrasolar planets

Perturbative post-Newtonian variations of the standard osculating orbital elements are obtained by using the two-body equations of motion in the parameterized post-Newtonian theoretical framework. The results obtained are applied to the Einstein and Brans–Dicke theories. As a results, the semi-major axis and eccentricity exhibit periodic variation, but no secular changes. The longitude of periastron and mean longitude at epoch experience both secular and periodic shifts. The post-Newtonian effects are calculated and discussed for six extrasolar planets.