Canonical solution of the two critical argument problem

The solution by Sessin and Ferraz-Mello (Celes. Mech. 32, 307–332) of the Hori auxiliary system for the motion of two planets with periods nearly commensurate in the ratio 21 is considerably simplified by the introduction of canonical variables. An analogous canonical transformation simplifies the elliptic restricted problem.     

Capture of dust grains in exterior resonances with planets

The temporary capture of the dust grains in the exterior resonances with planets is studied in the frames of the planar circular three-body problem with Poynting-Robertson (PR) drag. For the Earth and particles ~ 10 Μm the resonances 4/5, 5/6, 6/7, 7/8 are shown to be most effective. The capture is only temporary (of order 10⁵ years) and the position of resonance may be calculated from semi-analytical model using averaged disturbing function. These semi-analytical results are confirmed by numerical integration. For various planet this picture changes as with increasing planetary mass the more exterior resonances become more important. We showed that for Jupiter (at least in the space between Jupiter and Saturn) the resonance 1/2 plays the dominant role. The capture time is here several myr but again eccentricity is evolving to eccentricity e ₀ ~ 0.48 of libration point for this resonance.     

Families of periodic orbits in the planar three-body problem

A periodic orbit of the restricted circular three-body problem, selected arbitrarily, is used to generate a family of periodic motions in the general three-body problem in a rotating frame of reference, by varying the massm ₃ of the third body. This family is continued numerically up to a maximum value of the mass of the originally small body, which corresponds to a mass ratiom ₁:m ₂:m ₃?5:5:3. From that point on the family continues for decreasing massesm ₃ until this mass becomes again equal to zero. It turns out that this final orbit of the family is a periodic orbit of the elliptic restricted three body problem. These results indicate clearly that families of periodic motions of the three-body problem exist for fixed values of the three masses, since this continuation can be applied to all members of a family of periodic orbits of the restricted three-body problem. It is also indicated that the periodic orbits of the circular restricted problem can be linked with the periodic orbits of the elliptic three-body problem through periodic orbits of the general three-body problem.     

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.

     

Recent progress in the theory of Trojan asteroids

In previous publications the author has constructed a long-periodic solution of the problem of the motion of the Trojan asteroids, treated as the case of 1:1 resonance in the restricted problem of three bodies. The recent progress reported here is summarized under three headings:

1. The nature on the long-periodic family of orbits is re-examined in the light of the results of the numerical integrations carried out by Deprit and Henrard (1970). In the vicinity of the critical divisor $$D_k \equiv \omega _1 - k\omega _2 ,$$ not accessible to our solution, the family is interrupted by bifurcations and shortperiodic bridges. Parametrized by the normalized Jacobi constant α², our family may, accordingly, be defined as the intersection of admissible intervals, in the form $$L = \mathop \cap \limits_j \left\{ {\left| {\alpha - \alpha _j } \right| > \varepsilon _j } \right\};j = k,k + 1, \ldots \infty .$$ Here, {α_j(m)} is the sequence of the critical α_j corresponding to the exactj: 1 commensurability between the characteristic frequencies ω₁ and ω₂ for a given value of the mass parameterm. Inasmuch as the ‘critical’ intervals |α?α_j|<ε_j can be shown to be disjoint, it follows that, despite the clustering of the sequence {α_j} at α=1, asj→∞, the family extends into the vicinity of the separatrix α=1, which terminates the ‘tadpole’ branch of the family.
2. Our analysis of the epicyclic terms of the solution, carrying the critical divisorD _k , supports the Deprit and Henrard refutation of the E. W. Brown conjecture (1911) regarding the termination of the tadpole branch at the Lagrangian pointL ₃. However, the conjecture may be revived in a refined form. “The separatrix α=1 of the tadpole branch spirals asymptotically toward a limit cycle centered onL ₃.”
3. The periodT(α,m) of the libration in the mean synodic longitude λ in the range $$\lambda _1 \leqslant \lambda \leqslant \lambda _2$$ is given by a hyperelliptic integral. This integral is formally expanded in a power series inm and α² or \(\beta \equiv \sqrt {1 - \alpha ^2 }\) .

The large amplitude of the libration, peculiar to our solution, is made possible by the mode of the expansion of the disturbing functionR. Rather than expanding about Lagrangian pointL ₄, with the coordinatesr=1, θ=π/3, we have expandedR about the circler=1. This procedure is equivalent to analytic continuation, for it replaces the circle of convergence centered atL ₄ by an annulus |r?1|<ε with 0≤θ<2π.     

Parameter distribution of small periodic librations about the equilateral points of the elliptic restricted problem

As previously shown (Rabe, 1970), two classes of small periodic librations exist, in the plane, elliptic restricted problem, for an infinite sequence of easily specified oscillation frequenciesZ _j . The present paper considers the dependence ofZ on the eccentricitye of the primary motion, in addition to its dependence on the mass parameter , and determines the resulting relations between ande, for any given periodic frequencyZ _j . These relationships are obtained from the conditionD (Z _j ,, e)=0, where the basic determinantD has been expanded up to terms of orderZ ²⁰, ⁵, ande ⁴.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 27–September 2, 1972.     

The General Solution of the Planar Laplace Problem

The famous Laplace problem is the three-body, secular, planetary problem. Its plane version has the great theoretical advantage of being integrable (Ferraz-Mello, private correspondence, 2001) and Ferraz-Mello et al. (Chaotic World: from order to Disorder in Gravitational N-Body Systems, Kluwer Academic Publisher, 2004)). Nevertheless it remains a very complex problem with many singularities and many possibilities of collisions. Large eccentricities lead generally to large perturbations especially if the two planets have the same direction of revolution about their star.     

Periodic motions around a collinear equilibrium solution of the general three-body problem

The collinear equilibrium position of the circular restricted problem with the two primaries at unit distance and the massless body at the pointL ₃ is extended to the planar three-body problem with respect to the massm ₃ of the third body; the mass ratio μ of the two primaries is considered constant and the constant angular velocity of the straight line on which the three masses stay at rest is taken equal to 1. As regards periodic motions ‘around’ the equilibrium pointL ₃, four possible extensions from the restricted to the general problem are presented each of them starting with a simple or a doubly periodic orbit of the family α of the Copenhagen category (μ=0.50). Form ₃=0.10, μ=0.50 (i.e. for fixed masses of all three bodies) the characteristic curve of the extended family α is found. The qualitative differences of the families corresponding tom ₃=0 andm ₃=0.10 are discussed.     

Capture of dust grains in exterior resonances with planets

The temporary capture of the dust grains in the exterior resonances with planets is studied in the frames of the planar circular three-body problem with Poynting-Robertson (PR) drag. For the Earth and particles ~ 10 m the resonances 4/5, 5/6, 6/7, 7/8 are shown to be most effective. The capture is only temporary (of order 10⁵ years) and the position of resonance may be calculated from semi-analytical model using averaged disturbing function. These semi-analytical results are confirmed by numerical integration. For various planet this picture changes as with increasing planetary mass the more exterior resonances become more important. We showed that for Jupiter (at least in the space between Jupiter and Saturn) the resonance 1/2 plays the dominant role. The capture time is here several myr but again eccentricity is evolving to eccentricity e ₀ ~ 0.48 of libration point for this resonance.     

A third-order solution of Vinti's problem with explicit expressions for the poisson brackets

We have derived an explicit third-order solution of Vinti's problem including J₃. Poisson brackets for the elements a, e, s, M_s, ψ_s and $\text{Ω}{}_{\mathrm{s}}$ are given. These may be used in the construction of a third-order theory of artificial satellites.     

Chirikov diffusion in the asteroidal three-body resonance (5, −2, −2)

The theory of diffusion in many-dimensional Hamiltonian system is applied to asteroidal dynamics. The general formulation developed by Chirikov is applied to the Nesvorny-Morbidelli analytic model of three-body (three-orbit) mean-motion resonances (Jupiter-Saturn-asteroid). In particular, we investigate the diffusion along and across the separatrices of the (5, −2, −2) resonance of the (490) Veritas asteroidal family and their relationship to diffusion in semi-major axis and eccentricity. The estimations of diffusion were obtained using the Melnikov integral, a Hadjidemetriou-type sympletic map and numerical integrations for times up to 10⁸ years.     

Dust acoustic shock waves in a suprathermal dusty plasma with dust charge fluctuation

The properties of small but finite amplitude dust acoustic (DA) shock waves are studied in a charge varying dusty plasma with ions and electrons having kappa velocity distribution. We obtain the global Debye length including the influence of suprathermality effects and dust charge fluctuations. It is shown that the effects of suprathermality of ions/electrons and dust charge fluctuation significantly modify the basic properties of DA shock wave. We observe that only negative DA shock waves will be excited in this model. The amplitude of DA shock wave increases with deviation of electrons or ions from Maxwellian distribution via decrease of spectral index, κ _j (j=i,e denotes, ions and electrons, respectively). Also, it is indicated that the amplitude and steepness of the shock front decreases with an increase in the ion temperature.     

Expansion of the perturbing function of a satellite restricted four-body problem

A satellite four-body problem is the problem of motion of an artificial satellite of a planet in a region of the space where perturbations due to the gravitational field of the planet are of the same order as perturbations due to influences of two perturbing bodies. In this paper an expansion of the perturbing function into a Fourier series in terms of angular Keplerian elements ( _j , _j ,M _j :j=0,1,2) (designations are standard) is obtained taking into account a sharp commensurability of the typen/ ₀=(p+q)/p (n is the mean motion of the artificial satellite and ₀ is the angular velocity of rotation of the planet,p andq are integers).The coefficients of the Fourier series are the functions of the positional Keplerian elements (a _j ,e _j ,i _j ;j=0, 1, 2) (designations are standard) and, in particular, are series in terms ofe _j that, generally speaking, can be written out to an accuracy ofe _j ¹⁹ .The expansion obtained can be used for the construction of a semianalytical theory of motion of resonant satellites on the basis of conditionally periodic solutions of the restricted four-body problem.     

Expansion theory for the elliptic motion of arbitrary eccentricity and semi-major axis

In this paper of the series, the time transform and the explicit exact forms of the time will be established in terms of the sectorial variables _j ⁽ⁱ⁾ introduced in Paper IV (Sharaf, 1982) to regularize the highly oscillating perturbation force. Simple recurrence formulae are given to facilitate the computations. The formulations are general in the sense that they are valid whatever the types and the number of sectors forming the divisions situation of the elliptic orbit may be. Moreover, the constants of integration for the explicit forms of the time are determined in a way that it gives for these forms its generality during any revolution of the body in its Keplerian orbit.     

The elliptic restricted problem at the 3 : 1 resonance

Four 3 : 1 resonant families of periodic orbits of the planar elliptic restricted three-body problem, in the Sun-Jupiter-asteroid system, have been computed. These families bifurcate from known families of the circular problem, which are also presented. Two of them, I _c , II _c bifurcate from the unstable region of the family of periodic orbits of the first kind (circular orbits of the asteroid) and are unstable and the other two, I _e , II _e , from the stable resonant 3 : 1 family of periodic orbits of the second kind (elliptic orbits of the asteroid). One of them is stable and the other is unstable. All the families of periodic orbits of the circular and the elliptic problem are compared with the corresponding fixed points of the averaged model used by several authors. The coincidence is good for the fixed points of the circular averaged model and the two families of the fixed points of the elliptic model corresponding to the families I _c , II _c , but is poor for the families I _e , II _e . A simple correction term to the averaged Hamiltonian of the elliptic model is proposed in this latter case, which makes the coincidence good. This, in fact, is equivalent to the construction of a new dynamical system, very close to the original one, which is simple and whose phase space has all the basic features of the elliptic restricted three-body problem.     

Approximation functions to the Emden and associated emden functions near the first zero

Analytic expressions are presented which approximate the Emden function θ, and the associated Emden functions,ψ _o andψ ₂, and their derivatives, near the first zero of θ,ξ _o. The range of accurate representation extends toξ= 1.5ξ ₁. This range is sufficient to encompass the boundary of a critically rotating polytrope.     

Stability of L4 for radiated rigid primaries in resonant cases

The non-linear stability of the triangular libration point L₄ of the restricted three-body problem is studied under the presence of third- and fourth-order resonances, when the more massive primary is a triaxial rigid body and source of radiation. In this study, Markeev's theorems are applied with the help of Moser's theorem. It is found that the stability of the triangular libration point is unstable in the third-order resonance case and in the fourth-order resonance case, this is stable or unstable depending on A₁ and A₂, and a source of radiation parameter α, where A₁, A₂ depend upon the lengths of the semi-axes of the triaxial rigid body.     

Behaviour of Magnetoacoustic Waves in the Neighbourhood of a Two-Dimensional Null Point: Initially Cylindrically Symmetric Perturbations

The propagation of magnetoacoustic waves in the neighbourhood of a 2D null point is investigated for both β=0 and β ≠ 0 plasmas. Previous work has shown that the Alfvén speed, here v _A ∝r, plays a vital role in such systems and so a natural choice is to switch to polar coordinates. For β=0 plasma, we derive an analytical solution for the behaviour of the fast magnetoacoustic wave in terms of the Klein–Gordon equation. We also solve the system with a semi-analytical WKB approximation which shows that the β=0 wave focuses on the null and contracts around it but, due to exponential decay, never reaches the null in a finite time. For the β ≠ 0 plasma, we solve the system numerically and find the behaviour to be similar to that of the β=0 system at large radii, but completely different close to the null. We show that for an initially cylindrically-symmetric fast magnetoacoustic wave perturbation, there is a decrease in wave speed along the separatrices and so the perturbation starts to take on a quasi-diamond shape; with the corners located along the separatrices. This is due to the growth in pressure gradients that reach a maximum along the separatrices, which in turn reduces the acceleration of the fast wave along the separatrices leading to a deformation of the wave morphology.