Application of Hori's technique in general planetary theory

In this part we present the complete solution of the planetary canonical equations of motion by the method of G. Hori through successive changes of canonical variables using the Lie series. Thus, we can eliminate the long or critical terms of the planetary perturbing function, in our general planetary theory. In our formulas, we neglect perturbation terms of order higher than the third with respect to planetary masses.     

Progress in general planetary theory

Celestial Mechanics and Dynamical Astronomy - When on searches for a planetary theory valid over 1 million years, one can leave in the solution the short period terms whose amplitude are small, and...     

Secular perturbations in general planetary theory

Based on a general planetary theory, the secular perturbations in the motion of the eight major planets (excluding Pluto) have been derived in polynomial form. The results are presented in the tables. The linear terms of second order with respect to the planetary masses and the nonlinear terms of first order up to the fifth (and partly seventh) degree with respect to eccentricities and inclinations were taken into account in the right-hand members of the secular system. Calculations were carried out by computer with the use of a system that performed analytic operations on power series with complex coefficients.

---

qA ( ). . ( ) . .

     

Determination of secular terms in general planetary theory

A method to calculate secular terms of the two parts of the planetary disturbing function— when it is expressed in terms of the true anomalies or the eccentric anomalies instead of the mean anomalies - is described. Also an alternative method is outlined.     

The construction of a general third-order planetary theory

We review in this part the outline of a third-order general planetary theory established through Von Zeipel's method and in terms of Poincaré's canonical variables We consider our system to consist of the Sun as the primary body, one disturbed planet, and one disturbing planet.     

Elimination of the critical terms in a first-order general planetary theory

We construct a first order canonical general planetary theory, assuming the solar system to be composed of 8 planets excluding Pluto, referring to common fixed plane and applying the Jacobi-Radau set of origins. We eliminated by von Zeipel's method the 2:5 and 1:2 critical terms of Jupiter-Saturn and Uranus-Neptune inequalities. Our variables are those of Poincaré, and we expanded up to power three in the eccentricities and sines of the inclinations.     

Construction of a general planetary theory of the first order

All the necessary formulae for constructing a general solution for the motion of a planet, in rectangular coordinates, at the first order of the disturbing masses, in purely literal form in eccentricities and inclinations, are given. The authors present the transformation formulae in the two-body problem which give the correspondence between the constants of integration introduced in the theory and the classical keplerian elements. The practical elaboration of the algorithm and some partial results for the couple of planets Jupiter and Saturn are described.     

On the expansions of intermediate orbit for general planetary theory

Intermediate orbit for general planetary theory is constructed in the form of multivariate Fourier series with numerical coefficients. The structure and efficiency of the derived series are illustrated by giving various statistical properties of the coefficients.The ability of the recently proposed elliptic function approach to compress the Fourier series representing the intermediate orbit is investigated. Our results confirm that when mutual perturbations of a pair of planets are considered the elliptic function approach is quite efficient and allows one to compress the series substantially. However, when perturbations of three or more planets are under study the elliptic function approach does not give any advantages.     

Precession/nutation solution consistent with the general planetary theory

In the present paper the equations of the translatory motion of the major planets and the Moon and the Poisson equations of the Earth’s rotation in Euler parameters are reduced to the secular system describing the evolution of the planetary and lunar orbits (independent of the Earth’s rotation) and the evolution of the Earth’s rotation (depending on the planetary and lunar evolution). Hence, the theory of the Earth’s rotation is presented by means of the series in powers of the evolutionary variables with quasi-periodic coefficients.     

On the elimination of short-period terms in second-order general planetary theory investigated by Hori's method

1. The short-period terms of a second-order general planetary theory are removed through the Hori's method based on a development of the HamiltonianF in a Lie series which involves a determining functionS not depending upon mixed canonical variables as in the Von Zeipel's method but upon all the canonical variables resulting from the elimination of the short period terms ofF. Canonical variables adopted are the slow Delaunay variables. Eccentricitiese _j and sines γ_j of the semi inclinations are respectively replaced by the Jacques Henrard variablesE _j ,J _j which lead to formulas remarkably simple.F is reduced to the sumF ₀+F ₁ of its terms of degrees 0,1 in small parameter ε of the order of the masses. Only one disturbing planet is considered.F ₁ is not calculated beyond its terms of degree 3 inE _j ,E _j ,J _j , the determining functionS ₂ of degree 2 in ε not being therefore calculated beyond its terms of degree 2 inE′ _j ,E _j ,J _j and the expressions of slow Delaunay canonical variables of the disturbed planetP ₁ and the disturbing planetP ₂ in terms of the new slow Delaunay canonical variables ofP ₁ andP ₂ which result from the elimination of the short period terms ofF ₁ being therefore reduced to their terms of degree <1 in theE′ _j ,E′ _j ,J′ _j . Calculation of the principal partF _1m ofF ₁ is carried out through Laplace coefficients and operatorD=α(d/dα) applied to Laplace coefficients, α ratio of the semi major axis ofP ₁ andP ₂. Eccentricitye ₂ of the disturbed planetP ₂ is assumed to be zero, such an assumption not restricting our aim which is to investigate the mechanism of the elimination of short period terms in a second order general planetary theory carried out through the Hori's method, not to perform the elimination of those terms for a complete second order general planetary theory. Expressions of the slow Delaunay canonical variables in terms of the new ones resulting from the elimination of the short period terms ofF ₁ are written down only for the disturbed planetP ₁.
2. Small divisors in 1/E′ ₁ and 1/E′ ₁ ² appear in the longitude ?₁ of perihelia ofP ₁. No small divisors appear in the other five slow Delaunay variables ofP ₁. The only Jacques Henrard variables which appear in the longitude Ω₁ of the ascending node ofP ₁ are the J _j′ j=1, 2 and no Jacques Henrard variables appear in the slow Delaunay canonical variablesX ₁,Y ₁,Z ₁, λ₁. The solving of the ten canonical equations ofP ₁ andP ₂ in the slow Delaunay canonical variablesX′ _j ,Y′ ₁,Z′ _j ,λ′ _j ,ω′ _j ,Ω′ _j resulting from the elimination of the short period terms ofF ₁ reduces to that of four canonical equations inZ′ _j ,©′ _j and to six quadratures three of them expressing theX′ _j ,Y′ ₁ are constants and the three others expressingλ′ _j ,?′ _j as functions of timet. Solving of the four canonical equations inZ′ _j ,Ω′ _j reduces to that of a first order non linear differential equation and to two quadratures. Sinceγ′ ₁ is then constant, so is the Jacques Henrard variableE′ ₁. If the eccentricitye ₂ ofP ₂ is no more assumed to be zero, additive small divisors inE′ ₂/E′ ² ₁ appear in longitude ?′₁ of perihelia ofP ₁ and the solving of the twelve canonical equations ofP ₁ andP ₂ inX′ _j ,Y′ _j ,Z′ _j ,λ′ _j ,?′ _j ,Ω′ _j is reduced to that of eight canonical equations inY′ _j ,?′ _j ,Z′ _j ,Ω′ _j and to four quadratures expressingX′ _j are constants andλ′ _j as functions oft. Those eight canonical equations split into two systems of four canonical equations, one of them inY′ _j ,?′ _j and the other one inZ′ _j ,Ω′ _j . Each of those two systems is identical to the system inZ′ _j ,Ω′ _j corresponding toe ₂=0 and its solving reduces to that of a first order non linear differential equation and to two quadratures identical to those of the casee ₂=0.
3. Expressions ofX ₁,Y ₁,Z ₁,λ ₁,? ₁,Ω ₁ as functions ofX′ _j ,Y′ ₁,Z′ _j ,λ′ _j ,?′ ₁,Ω′ _j ;j=1, 2 are sums of sines and cosines of the multiples ofλ′ _j ,?′ ₁,Ω′ _j for the terms arising from the indirect partF _1j ofF ₁, Fourier series in those sines and cosines or products of two such Fourier series for the terms arising from the principal partF _1m ofF ₁, coefficients of those sums and Fourier series having one of the eight forms: $$A,{\text{ }}\frac{B}{{E'}},{\text{ }}\frac{C}{{E'^2 }},{\text{ }}D\frac{{j'^{2_1 } }}{{E'^{2_1 } }},{\text{ }}E\frac{{j'^{2_2 } }}{{E'^{2_1 } }},{\text{ }}F\frac{{j'^{_1 } j'^2 }}{{E'^{2_1 } }},{\text{ }}G\frac{{j'^2 }}{{j'^{_1 } }},{\text{ }}H\frac{{j'^{22} }}{{j'^{2_1 } }}{\text{.}}$$ A,..., H being constants which depend upon ratio α. Numerical calculation of the constantsA,..., H arising from the terms ofF _1j is easily carried out; that of theA,..., H arising from the terms ofF _1m require more manipulations, Fourier series in sines and cosines of the multiples ofλ′ _j ,?′ _j ,Ω _ij and products of two such Fourier series having then to be reduced to sums of a finite number of terms and treated through the methods of harmonic analysis. Divisors inp+qα^3/2;p, q relative integers, or products of such divisors appear inA,..., H.
4. the method extends to the case whenF ₁ is calculated beyond its terms of degree 3 in the Jacques Henrard variables.F ₁ being calculated up to its terms of degree 8 in the Jacques Henrard variables which is the precision required to eliminate the short period terms of a complete second order general planetary theory,S ₂ has to be calculated up to its terms of degree 7 and the expression of the slow Delaunay canonical variables ofP ₁ andP ₂ in terms of the slow Delaunay canonical variables ofP ₁ andP ₂ resulting from the elimination of the short period terms ofF ₁ have, therefore, to be calculated up to their terms of degree 5 in the Jacques Henrard variables.

     

The elimination of short-periodic terms in a Uranus-Neptune general planetary theory by von Zeipel's method

We eliminate by the method of von Zeipel the short-period terms in a first order-with respect to planetary masses—general planetary Uranus-Neptune theory. We exclude in the expansion terms of eccentricities and sines of inclinations higher than the third power.Our variables are the Poincaré canonical variables. We use the Jacobi-Radau set of origins, and we refer the planes of the osculating ellipses to a common fixed plane, the longitudes to a common origin. The short-periodic terms arising from the indirect and principal parts of the disturbing functions, are eliminated separately. The Fourier series of the principal part of the disturbing function, is reduced to the sum of only the first three terms.     

Third-order Uranus-Neptune planetary theory

In this part we calculate the secular and critical terms arising from the indirect part of the classical planetary Hamiltonian for Uranus and Neptune. We neglect in our expansions powers higher than the second in the eccentricity-inclination. Our required results, are expressed in terms of Poincaré variables.     

General planetary theory in elliptic functions

It is shown that the first-order general planetary theory, i.e. the theory without secular terms, developed in (Brumberg and Chapront, 1973) may be re-constructed and presented by the series in powers of the eccentricity and inclination variables with the closed form coefficients expressed in terms of elliptic functions. The intermediate solution of the zero degree in eccentricities and inclinations has been given explicitly with the aid of elliptic functions and the Hansen type quadratures with trigonometric function kernels. In determining the first and higher degree terms in eccentricities and inclinations one meets the Hansen type quadratures with elliptic function kernels. The secular evolution is described by the autonomous polynomial differential system.     

Third order uranus-neptune planetary theory

We calculate in this paper the secular and critical terms arising from the principal part of the classical planetary Hamiltonian. This is the first step to establish a third order canonical planetary theory of Uranus-Neptune through the Hori-Lie technique. We truncate our expansions at the second degree of eccentricity-inclination. Our planetary theory is expressed in terms of the canonical variables of H. Poincaré.     

The floccule theory for planetary formation

The accumulation of floccules into protoplanets is discussed, and it is pointed out that the simplifications which have been introduced into recent numerical models may result in the incorrect conclusion being reached.     

A second order Jupiter-Saturn planetary theory

We complete by this part II the establishment of a second order secular Jupiter-Saturn theory. This is achieved by taking into consideration the influence of the indirect part of the planetary disturbing function, and expressing the second order secular Hamiltonian in terms of Poincaré's canonical variables.     

Third-order Jupiter — Saturn planetary theory

The construction of a third order J-S theory is presented. The Hori theory of planetary perturbations is employed. No Critical J-S terms due to the 2:5 commensurabilities and its multiples exist, when we take into account the periodic terms of order 0, 1, 2 with respect to the eccentricity- inclination. In this case the Lie series transformation degenerates and is meaningless. The J-S equations of motion for secular perturbations are solved when we neglect in our treatment, the Poisson terms of degree > 2 in the Poincaré canonical variables H _u , K _u , P _u Q _u (u = 1, 2). The Jacobi-Radau referential is adopted, and the theory is expressed in terms of the canonical variables of H. Poincaré.Now at the Jet Propulsion Laboratory, California Institute of Technology, Pasadena, California, U.S.A.     

A second order Jupiter-Saturn planetary theory

We present a second order secular Jupiter-Saturn planetary theory through Poincaré canonical variables, von Zeipel's method and Jacobi-Radau referential. We neglect in our expansions terms of power higher than the fourth with respect to eccentricities and sines of inclinations. We assume that the disturbing function is composed of secular and critical terms only. We shall deriveF _2si and writeF _2s in terms of Poincaré canonical variables in Part II of this problem.     

An analytical theory of planetary rotation rates

An approximate analytical theory is derived for the rate of rotation acquired by a planet as it grows from the solar nebula. This theory was motivated by a numerical study by Giuli, and yields fair agreement with his results. The periods of planetary rotation obtained are proportional to planetesimal encounter velocity, and appear to suggest lower values of this velocity than are commonly assumed to have existed during planetary formation.     

A determination of secular terms in first-order planetary theory

The secular terms of the first-order planetary Hamiltonian is determined, by two methods, in terms of the variables of H. Poincaré, neglecting powers higher than the second in the eccentricity-inclination.