Motion near the 3/1 resonance of the planar elliptic restricted three body problem

The global semi-numerical perturbation method proposed by Henrard and Lemaître (1986) for the 2/1 resonance of the planar elliptic restricted three body problem is applied to the 3/1 resonance and is compared with Wisdom's perturbative treatment (1985) of the same problem. It appears that the two methods are comparable in their ability to reproduce the results of numerical integration especially in what concerns the shape and area of chaotic domains. As the global semi-numerical perturbation method is easily adapted to more general types of perturbations, it is hoped that it can serve as the basis for the analysis of more refined models of asteroidal motion. We point out in our analysis that Wisdom's uncertainty zone mechanism for generating chaotic domains (also analysed by Escande 1985 under the name of slow Hamiltonian chaotic layer) is not the only one at work in this problem. The secondary resonance _p = 0 plays also its role which is qualitatively (if not quantitatively) important as it is closely associated with the random jumps between a high eccentricity mode and a low eccentricity mode.     

Non-linear stability of L 4 in the restricted problem when the primaries are finite straight segments under resonances

The non-linear stability of L ₄ in the restricted three-body problem when both primaries are finite straight segments in the presence of third and fourth order resonances has been investigated. Markeev’s theorem (Markeev in Libration Points in Celestial Mechanics and Astrodynamics, 1978) is used to examine the non-linear stability for the resonance cases 2:1 and 3:1. It is found that the non-linear stability of L ₄ depends on the lengths of the segments in both resonance cases. It is also found that the range of stability increases when compared with the classical restricted problem. The results have been applied in the following asteroids systems: (i) 216 Kleopatra–951 Gaspara, (ii) 9 Metis–433 Eros, (iii) 22 Kalliope–243 Ida.     

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.     

Sur de nouvelles séries pour le problème de masses critiques de Routh dans le problème restreint plan des trois corps

The planar restricted 3-body problem, linearized in the neighborhood of Lagrangian equilibriaL ₄ andL ₅, has in general two distinct eigenvalues and their opposites. When they are pure imaginary and not multiples of each other, they generate two families of periodic solutions called long and short periodic families. This is essentially a consequence of the famous theorem of Liapunov (Siegel, 1956). We showed (Roels, 1971b) how to solve the problem when the eigenvalues are multiples of each other in building series with negative exponents instead of the integer expansions of Siegel (Roels and Lauterman, 1970). When the eigenvalues are equal, which is the case for the mass ratio of Routh, the problem was solved by Deprit and Henrard (1968) using formal series in ordinary unnormalized variables. That leads to very complicated series because of the use of variables that are not well adapted to the problem. The convergence of the series was proven by Meyer and Schmidt (1971). In this paper we solve the problem by using normalized variables. This brings us to build expansions with fractional exponents. So in summary, normalized variables generate integer series in the non-resonant cases, series with negative exponents in the case of resonancek≥3, and series with fractional exponents when the resonance is 1.     

Effect of oblateness on the non-linear stability of the triangular liberation points of the restricted three-body problem in the presence of resonances

The non-linear stability of the triangular libration points of the restricted three-body problem is studied under the presence of third and fourth order resonance's, when the more massive primary is an oblate spheroid. In this study Markeev's theorem are utilised with the help of KAM theorem. It is found that the stability of the triangular libration points are unstable in the third order resonance case and stable in the fourth order resonance case, for all the values of oblateness factor A₁. This revised version was published online in August 2006 with corrections to the Cover Date.     

A problem of libration with two degrees of freedom

Moore (1983) presented a theory of resonance with two degrees of freedom based on the Bohlin-von Zeipel procedure. This procedure is now applied to librational motion with all constants re-evaluated in terms of values of the momenta given either by the initial conditions, or, in the case of the momentumy ₁ conjugate to the critical argument x₁, by its value at the libration centre. Numerical results are presented for a resonant satellite in a near 12 hr orbit and for a geosynchronous satellite. The theory is further developed to include near-circular orbits by recasting the problem in terms of the Poincaré eccentric variables.     

On the problem of the flare-like emission in sodium lines

Unusual and intense emission in the D ₁, D ₂ lines has been registered photoelectrically in an undisturbed region near the center of the solar disk on 28 July, 1966. Various aspects of the possibility that the emission might originate in either artificial or natural comets by processes of resonance scattering or flu orescence have been analyzed. It is shown that the observed emission can only be located on the solar surface. Some aspects of the problem are discussed briefly.     

On the oscillations in Mercury's obliquity

Mercury's capture into the 3:2 spin–orbit resonance can be explained as a result of its chaotic orbital dynamics. One major objective of MESSENGER and BepiColombo spatial missions is to accurately measure Mercury's rotation and its obliquity in order to obtain constraints on internal structure of the planet. Analytical approaches at the first-order level using the Cassini state assumptions give the obliquity constant or quasi-constant. Which is the obliquity's dynamical behavior deriving from a complete spin–orbit motion of Mercury simultaneously integrated with planetary interactions? We have used our SONYR model (acronym of Spin–Orbit N-bodY Relativistic model) integrating the spin–orbit N-body problem applied to the Solar System (Sun and planets). For lack of current accurate observations or ephemerides of Mercury's rotation, and therefore for lack of valid initial conditions for a numerical integration, we have built an original method for finding the libration center of the spin–orbit system and, as a consequence, for avoiding arbitrary amplitudes in librations of the spin–orbit motion as well as in Mercury's obliquity. The method has been carried out in two cases: (1) the spin–orbit motion of Mercury in the 2-body problem case (Sun–Mercury) where an uniform precession of the Keplerian orbital plane is kinematically added at a fixed inclination (S_2K case), (2) the spin–orbit motion of Mercury in the N-body problem case (Sun and planets) (S_n case). We find that the remaining amplitude of the oscillations in the S_n case is one order of magnitude larger than in the S_2K case, namely 4 versus 0.4 arcseconds (peak-to-peak). The mean obliquity is also larger, namely 1.98 versus 1.80 arcminutes, for a difference of 10.8 arcseconds. These theoretical results are in a good agreement with recent radar observations but it is not excluded that it should be possible to push farther the convergence process by drawing nearer still more precisely to the libration center. We note that the dynamically driven spin precession, which occurs when the planetary interactions are included, is more complex than the purely kinematic case. Nevertheless, in such a N-body problem, we find that the 3:2 spin–orbit resonance is really combined to a synchronism where the spin and orbit poles on average precess at the same rate while the orbit inclination and the spin axis orientation on average decrease at the same rate. As a consequence and whether it would turn out that there exists an irreducible minimum of the oscillation amplitude, quasi-periodic oscillations found in Mercury's obliquity should be to geometrically understood as librations related to these synchronisms that both follow a Cassini state. Whatever the open question on the minimal amplitude in the obliquity's oscillations and in spite of the planetary interactions indirectly acting by the solar torque on Mercury's rotation, Mercury remains therefore in a stable equilibrium state that proceeds from a 2-body Cassini state.     

Analytical investigation of non-linear stability of the Lagrangian pointL 4 around the commensurability 1:2

The problem of stability of the Lagrangian pointL ₄ in the circular restricted problem of three bodies is investigated close to the 1 : 2 commensurability of the long and short period libration. By stability we define boundedness of the solution for a given initial finite displacement from the equilibrium point as function of the mass parameter close to the commensurability. A rigorous treatment close to the resonance condition is possible using a transformation that diagonalizes the matrix related to the linear part of the equations of motion. The so obtained equations are further transformed to action angle type variables. Then using an isolated resonance approach, only the slowly varying terms are kept in the equations and two independent isolating first integrals can be found. These integrals finally enable us to solve the stability problem in an exact way. The so obtained results are compared to numeric integration of the equations of motion and are found to be in perfect agreement.     

Resonance stability of triangular equilibrium points in elliptical restricted three body problem under the radiating primaries

The main aim of this paper is to study the existence of resonance and linear stability of the triangular equilibrium points of the planar elliptical restricted three body problem considering the photo gravitational effect of both the primaries in circular and elliptical case. A practical application of this case could be the study of the dynamical system around the binary systems. For this the Hamiltonian function, convergent in nature and describing the motion of the infinitesimal body in the neighborhood of the triangular equilibrium solutions is derived. Also, the Hamiltonian for the system is expanded in powers of the generalized components of momenta. Further, canonical transformation has also been used to study the stability of the triangular equilibrium points. The study primarily focuses on establishing the relation for determining the range of stability at and near the resonance frequency ω ₂=1/2 around the binary systems using simulation technique. It is observed that the parametric resonance is only possible at the resonance frequency ω ₂=1/2 in both circular and elliptical cases.     

On the global solution in the resonance problem of Poincaré

Poincaré formulated a general problem of resonance in the case of a dynamical system which is reducible to one degree of freedom. He introduced the concept of a global solution; in essence, this means that the domain of the solution(s) covers the entire phase plane, comprising regions of libration and circulation. It is the author's opinion that the technique proposed by Poincaré for the construction of a global solution is impractical. Indeed, in §§201 and 211 ofLes méthodes nouvelles de la méchanique céleste, where he describes the passage from shallow resonance to deep resonance, Poincaré asserts an erroneous conclusion. An alternative procedure, which admits secular terms into the determining function and introduces a regularizing function, is outlined. The latter method has been successfully applied to the Ideal Resonance Problem, which is a special case of the more general problem considered by Poincaré, (Garfinkelet al. (1971); Garfinkel (1972).     

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.     

Galilean satellites: Evolutionary paths in deep resonance

The Laplace resonance among the inner three Galilean satellites (mean motions n₁ − 3n₂ + 2n₃ = 0) has stable configurations in “deep resonance,” i.e., where mean motions taken by pairs are in ratios very close to 2:1. The present satellite configuration, with the resonance variable φ ≡ λ₁ − 3λ₂ + 2λ₃ stable at 180°, is unstable near this exact commensurability. But there is a continuous path of stable conditions branching from φ = 180° to higher and lower values of φ and toward very deep resonance, according to a theory extended to third order in orbital eccentricity. This path provides a track for tidal evolution of the system. Thus, scenarios involving evolution (probably episodic) from deep resonance are viable, and eliminate the requirement by the alternative equilibrium hypothesis for rapid tidal dissipation in Jupiter. Evolution out from deep resonance is consistent with the free eccentricity of Ganymede, the free libration of φ, and observational constraints on Io's secular acceleration. Also, the relatively large forced eccentricities in deep resonance may have controlled geophysical processes in the satellites by much greater tidal heating and global stress than at present.     

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.     

Stability of L4 and L5 against radiation pressure

The restricted three-body problem is generalized with the inclusion of solar radiation pressure. For small particles (typically 1 m to 1 mm) the familiar equilibrium triangular points L₄ and L₅ no longer exist. However libration orbits are not completely destroyed, although an effect of resonance causes their amplitude to be very large, for a particle initially at rest at either of the triangular point. Finally the results of a study of the linearized equations of motion, supplemented by a numerical integration, rule out the possibility of an accumulation of dust at the Earth-Moon lagrangian triangular points.     

On The Ideal Resonance Problem

The Ideal Resonance Problem in its normal form is defined by the Hamiltonian (1) $$F = A (y) + 2B (y) sin^2 x$$ with (2) $$A = 0(1),B = 0(\varepsilon )$$ where ? is a small parameter, andx andy a pair of canonically conjugate variables. A solution to 0(?^1/2) has been obtained by Garfinkel (1966) and Jupp (1969). An extension of the solution to 0(?) is now in progress in two papers ([Garfinkel and Williams] and [Hori and Garfinkel]), using the von Zeipel and the Hori-Lie perturbation methods, respectively. In the latter method, the unperturbed motion is that of a simple pendulum. The character of the motion depends on the value of theresonance parameter α, defined by (3) $$\alpha = - A\prime /|4A\prime \prime B\prime |^{1/2} $$ forx=0. We are concerned here withdeep resonance, (4) $$\alpha< \varepsilon ^{ - 1/4} ,$$ where the classical solution with a critical divisor is not admissible. The solution of the perturbed problem would provide a theoretical framework for an attack on a problem of resonance in celestial mechanics, if the latter is reducible to the Ideal form: The process of reduction involves the following steps: (1) the ration ₁/n₂ of the natural frequencies of the motion generates a sequence. (5) $$n_1 /n_2 \sim \left\{ {Pi/qi} \right\},i = 1, 2 ...$$ of theconvergents of the correspondingcontinued fraction, (2) for a giveni, the class ofresonant terms is defined, and all non-resonant periodic terms are eliminated from the Hamiltonian by a canonical transformation, (3) thedominant resonant term and itscritical argument are calculated, (4) the number of degrees of freedom is reduced by unity by means of a canonical transformation that converts the critical argument into an angular variable of the new Hamiltonian, (5) the resonance parameter α (i) corresponding to the dominant term is then calculated, (6) a search for deep resonant terms is carried out by testing the condition (4) for the function α(i), (7) if there is only one deep resonant term, and if it strongly dominates the remaining periodic terms of the Hamiltonian, the problem is reducible to the Ideal form.     

On the spiral structure of our Galaxy

It is shown in the present paper that properties of the spiral wave in the Galaxy are determined by the mass distribution of its flat subsystem rather than by the full mass distribution. Then it turns out that better agreement with the observed spiral pattern furnish the ‘long’ waves in contrast to the ‘short’ waves in the Linet al. (1969) theory. With the surface density σI=40M _⊙ /ps ₂ which is taken in the first approximation as independent on the galacto-centric distance, and the pattern velocityΩ _p=23 km/s kps, the evaluated spiral pattern fits surprisingly well with the Weaver (1970) map of the HI-distribution in the Galaxy, and is in good agreement with the Kerr (1969) map. The inner Lindblad resonance occurs at 2 kps from the Galaxy center, where Weaver has placed the ring condensation of the gas, and the outer resonance lies close by 14 kps. At the outer resonance the nonlinear phenomena are expected, which lead to chaotization of the regular structure. This seems to be consistent with the Weaver (1970) and Kerr (1969) maps. The hypothesis is suggested which associates the generating mechanism of spiral waves with the rotating bar of old stars in the center of the Galaxy. Depending on the velocity of the bar rotation and the bar length, different combinations of the normal wave pattern and bar-like structure may occur, which possibly explains the great variety of transition forms between normal and barred spirals. In the proposed theory the packet of spiral waves moves from the inner Lindblad resonance outwards and could be permanently maintained by the ‘generator’ in the center of the Galaxy. Therefore, the difficulty associated with the rapid obliteration of the packet (Toomre, 1969) does not arise.     

Normality condition in the ideal resonance problem

The Ideal Resonance Problem, defined by the Hamiltonian $$F = B(y) + 2\mu ^2 A(y)\sin ^2 x,\mu \ll 1,$$ has been solved in Garfinkelet al. (1971). As a perturbed simple pendulum, this solution furnishes a convenient and accurate reference orbit for the study of resonance. In order to preserve the penduloid character of the motion, the solution is subject to thenormality condition, which boundsAB" andB' away from zero indeep and inshallow resonance, respectively. For a first-order solution, the paper derives the normality condition in the form $$pi \leqslant max(|\alpha /\alpha _1 |,|\alpha /\alpha _1 |^{2i} ),i = 1,2.$$ Herep _i are known functions of the constant ‘mean element’y', α is the resonance parameter defined by $$\alpha \equiv - {\rm B}'/|4AB\prime \prime |^{1/2} \mu ,$$ and $$\alpha _1 \equiv \mu ^{ - 1/2}$$ defines the conventionaldemarcation point separating the deep and the shallow resonance regions. The results are applied to the problem of the critical inclination of a satellite of an oblate planet. There the normality condition takes the form $$\Lambda _1 (\lambda ) \leqslant e \leqslant \Lambda _2 (\lambda )if|i - tan^{ - 1} 2| \leqslant \lambda e/2(1 + e)$$ withΛ ₁, andΛ ₂ known functions of λ, defined by $$\begin{gathered} \lambda \equiv |\tfrac{1}{5}(J_2 + J_4 /J_2 )|^{1/4} /q, \hfill \\ q \equiv a(1 - e). \hfill \\ \end{gathered}$$      

Nonlinear dynamical analysis for displaced orbits above a planet

Nonlinear dynamical analysis and the control problem for a displaced orbit above a planet are discussed. It is indicated that there are two equilibria for the system, one hyperbolic (saddle) and one elliptic (center), except for the degenerate h _z ^max, a saddle-node bifurcation point. Motions near the equilibria for the nonresonance case are investigated by means of the Birkhoff normal form and dynamical system techniques. The Kolmogorov–Arnold–Moser (KAM) torus filled with quasiperiodic trajectories is measured in the τ ₁ and τ ₂ directions, and a rough algorithm for calculating τ ₁ and τ ₂ is proposed. A general iterative algorithm to generate periodic Lyapunov orbits is also presented. Transitions in the neck region are demonstrated, respectively, in the nonresonance, resonance, and degradation cases. One of the important contributions of the paper is to derive necessary and sufficiency conditions for stability of the motion near the equilibria. Another contribution is to demonstrate numerically that the critical KAM torus of nontransition is filled with the (1,1)-homoclinic orbits of the Lyapunov orbit.