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.     

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.     

On the elimination of the critical term in a general first-order Uranus-Neptune theory

A solution of the Uranus-Neptune planetary canonical equations of motion through the Von Zeipel technique is presented. A unique determinging function which depends upon mixed canonical variables, reduces the 12 critical terms of the Hamiltonian to the set of its secular terms. The Poincaré canonical variables are used. We refer to a common fixed plane, and apply the Jacobi-Radau set of origins. In our expansion we neglected terms of power higher than the fourth with respect to the eccentricities and sines of the inclinations.     

An elimination of the critical terms of a first-order Uranus-Neptune theory by Hori's method

We eliminate the 1:2 critical terms — after a previous elimination of the short period terms — in the Hamiltonian of a first order U-N theory. We take into account terms of degree 0, 1, 2, 3, 4 in the eccentricity-inclination. We apply for this elimination the Hori-Lie technique through the Poincaré canonical variables and the Jacobi coordinates. The purely principal first order secular U-N Hamiltonian admits a complete solution. We obtained the U-N equations of motion generated by the principal first order long period U-N Hamiltonian which will be solved later. This part III is closely related to the two previous papers (Kamel, 1982, 1983).     

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.     

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.

     

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 ( ). . ( ) . .

     

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.     

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.     

Application of Hori technique in general planetary theory

We explain how the first step of Hori-Lie procedure is applied in general planetary theory to eliminate short-period terms. We extend the investigation to the third-order planetary theory. We solved the canonical equations of motion for secular and periodic perturbations by this method, and obtained the first integrals of the system of canonical equations. Also we showed the relation between the determining function in the sense of Hori and the determining function in the sense of Von Zeipel.     

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.     

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.     

Solution of canonical equations which result from elimination of short-periodic terms in second-order planetary theory

We solve the first order non-linear differential equation and we calculate the two quadratures to which are reduced the canonical differential equations resulting from the elimination of the short period terms in a second order planetary theory carried out through Hori's method and slow Delaunay canonical variables when powers of eccentricities and the sines of semi-inclinations which are >3 are neglected and the eccentricity of the disturbing planet is identically equal to zero. The procedure can be extended to the case when the eccentricity of the disturbing planet is not identically equal to zero. In this latter general case, we calculatedthe two quadratures expressing angular slow Delaunay canonical variable ₁ of the disturbed planet and angular slow Delaunay canonical variable ₂ of the disturbing planet in terms of timet.     

The elimination of the critical terms of a first order Uranus-Neptune theory through Hori's method

An outline for the elimination of the critical terms of a first order Uranus-Neptune theory is presented with a stress on the application of Hori's procedure to the problem.     

The elimination of the critical terms of a first order Uranus-Neptune theory by Hori's method

In this part we determine the value ofS ₁, and in terms of the canonical variables of H. Poincaré. A complete solution of the auxiliary system of equations generated by the Hamiltonian is presented.     

Expansion of the first-order planetary disturbing function by Smart's method

We expand the principal part of the planetary disturbing function, by Smart's method, using Taylor's theorem. In our expansion we neglect terms of degree higher than the fourth in the eccentricities and tangents of the inclinations.Now at the JPL Pasadena, California.     

Evidence for convection in planetary interiors from first-order topography

Center of mass-center of figure offsets are known for the Earth, Moon, Mars and Venus. Such an offset requires a density distribution asymmetric about the center of mass. Observational evidence indicates that the terrestrial, lunar and Martian offsets result from crusts of variable thickness rather than lateral density inhomogeneities and that the thickness variations are more likely caused by internal convection than impact.Paper dedicated to Professor Harold C. Urey on the occasion of his 80th birthday on 29 April, 1973.     

Elimination of the gravitational singularities within the framework of Einstein's general theory of relativity

In this paper two examples are given, showing that the existence of black holes in the Universe violates, in its consequences, the principle of causality. The solution presented is based on the idea that the primordial black holes have zero-mass-energy and consequently zero-radius of the event horizon. Despite the existence of the surface of last influence, gravitational collapse does not produce black holes during a finite time interval as measured by an external observer. The only singularity, possible to accept (if any), is the initial and final cosmological singularity.