The integrable cases of the planetary three-body problem at first-order resonance

This paper investigates the stability of the motion in the averaged planar general three-body problem in the case of first-order resonance. The equations of the averaged motion of bodies near the resonance surface is obtained and is analytically integrated by quadratures. The stability of the averaged motion is analytically investigated in relation to the semi-major axes, the eccentricities and the resonance phases. An autonomous second-order equation is obtained for the deviation of semiaxes from the resonance surface. This equation has an energy integral and is analytically integrated by quadratures. The quasi-periodic dependence on time with two-frequency basis of the averaged motion of bodies is found. The basic frequencies are analytically calculated. With the help of the mean functionals calculated along integral curves of the averaged problem the new analytic first integrals are constructed with coefficients periodic in time. The analytic conditions of librations of resonance phases are obtained.     

Disappearance of integrals in systems of more than two degrees of freedom

The disappearance of some integrals of motion when two or more resonance conditions are approached at the same time is explained. As an example a Hamiltonian of three degrees of freedom is considered in action-angle variables which in zero order represents three harmonic oscillators, while the perturbation contains two trigonometric terms. One integral disappears if two appropriate resonant conditions are approached sufficiently closely.     

About the first integrals of the generalized problem of translatory-rotary motion of rigid bodies

The existence of ten first integrals for the classical problem of the motion of a system of material points, mutually attracting according to Newtonian law, is well known.The existence of the analogous ten first integrals for the more complicated problem of the motion of a system of absolutely rigid bodies, whose elementary particles mutually attract according to the Newtonian law, was established by the author (Duboshin, 1958, 1963, 1968).In his later papers (Duboshin, 1969, 1970), the problem of the motion of a system of material points, attracting each other according to a more general law, was considered and, in particular, it was shown under what conditions the ten first integrals, analogous to the classical integrals, may exist for this problem.In the present paper, the generalized problem of translatory-rotatory motion of rigid bodies, whose elementary particles acting upon each other according to arbitrary laws of forces along the straight line joining them, is discussed.The author has shown that the first integrals for this general problem, analogous to the integrals of the problem of the translatory-rotatory motion of rigid bodies, whose elementary particles acting according to the Newtonian law, exist under certain well known conditions.That is, it has been established that if the third axiom of dynamics (action = reaction) is satisfied, then the integrals of the motion of centre of inertia and the integrals of the moment of momentum exist for this generalized problem.If the third axiom is not satisfied, then the above mentioned integrals do not exist.The third axiom is a necessary but not a sufficient condition for the existence of the tenth integral-the energy integral. The tenth integral always exists if the elementary particles of the bodies acting with a force, depend only on the mutual distances between them. In this case the force function exists for the problem and the energy integral can be expressed in a well known form.The tenth integral may exist for some more general case, without expressing the principle of conservation of energy, but permitting calculation of the kinetic energy, if the configuration of a system is given.The problem, in which the elementary particles acting according to the generalized Veber's law (Tisserand, 1896) has been cited as an example of this more general case.     

Non—integrability of Hill's lunar problem

We prove that Hill's lunar problem does not possess a second analytic integral of motion, independent of the Hamiltonian. In order to obtain this result, we avoid the usual normalization in which the angular velocity of the rotating reference frame is put equal to unit. We construct an artificial Hamiltonian that includes an arbitrary parameter b and show that this Hamiltonian does not possess an analytic integral of motion for in an open interval around zero. Then, by selecting suitable values of , b and using the invariance of the Hamiltonian under scaling in the units of length and time, we show that the Hamiltonian of Hill's problem does not possess an integral of motion, analytically continued from the integrable two–body problem in a rotating frame.     

About the third integral of charged particle motion in strongly curved magnetic fields

Charged particle motion is studied in magnetic fields with an inversion of field direction and a strong curvature in the reversal region, that means field structures which are expected to play a crucial role in energetics and dynamics of space plasma. Investigations are performed in a typical field model. Due to its symmetry and stationarity two of the integrals of motion of the resulting Hamiltonian immediately arise. The question is analyzed whether the system is integrable (this means regular solutions exist) and of what kind the third integral is. Two different third integrals are found, which are valid in different parts of the phase space. Their validity is estimated analytically for particles, crossing the neutral plane (z = 0) -- in addition the action integral I_z previously unusual in plasma applications is verified numerically. Further numerical research establishes that there are no more third integrals -- all remaining parts of the phase space are filled with chaotic solutions. Because in strongly curved field reversals the conservation of the magnetic moment as an integral of motion is shown to be restricted to very small energies or very large pitch angles, the unusual I_z-integral will be an important tool for solution of plasma problems of cosmical current sheets and plasma boundaries.     

On the number of effective integrals in galactic models

Three different numerical techniques are tested to determine the number of integrals of motion in dynamical systems with three degrees of freedom.

1. The computation of the whole set of Lyapunov Characteristic Exponents (LCE).
2. The triple sections in the configurations space.
3. The Stine-Noid box-counting technique.

These methods are applied to a triple oscillator with coupling terms of the third order. Cases are found for which one effective integral besides the Hamiltonian subsists during a very long time. Such orbits display simultaneously chaotic and quasi-periodic motion, according to which coordinates are considered. As an application, the LCE procedure is applied to a triaxial elliptical galaxy model. Contrary to similar 2-dimensional systems, this 3-dimensional one presents noticeable zones in the phase space without any non-classical integral.     

Averaging the Equations of a Planetary Problem in an Astrocentric Reference Frame

A system of averaged equations of planetary motion around a central star is constructed. An astrocentric coordinate system is used. The two-planet problem is considered, but all constructions are easily generalized to an arbitrary number N of planets. The motion is investigated in modified (complex) Poincarécanonical elements. The averaging is performed by the Hori–Deprit method over the fast mean longitudes to the second order relative to the planetary masses. An expansion of the disturbing function is constructed using the Laplace coefficients. Some terms of the expansion of the disturbing function and the first terms of the expansion of the averaged Hamiltonian are given. The results of this paper can be used to investigate the evolution of orbits with moderate eccentricities and inclinations in various planetary systems.     

A Map for a Group of Resonant Cases in a quartic Galactic Hamiltonian

We present a map for the study of resonant motion in a potential made up of two harmonic oscillators with quartic perturbing terms. This potential can be considered to describe motion in the central parts of non-rotating elliptical galaxies. The map is based on the averaged Hamiltonian. Adding on a semi-empirical basis suitable terms in the unperturbed averaged Hamiltonian, corresponding to the 1:1 resonant case, we are able to construct a map describing motion in several resonant cases. The map is used in order to find thex − p _x Poincare phase plane for each resonance. Comparing the results of the map, with those obtained by numerical integration of the equation of motion, we observe, that the map describes satisfactorily the broad features of orbits in all studied cases for regular motion. There are cases where the map describes satisfactorily the properties of the chaotic orbits as well.     

Influence of the resonance in a gravity-gradient stabilized satellite

Using Hamiltonian formalism the translational-rotational motion of a satellite is studied near a resonance considering the orbital and rotational motions. A first order perturbation theory is derived by Hori's transformation in order to eliminate short and long periodic terms, preserving in the new Hamiltonian secular and resonant terms. This theory is again applied to study the resonant system whose analysis lead us to a system of equations equivalent to the equations of a simple pendulum which is integrable in terms of elliptical integrals.     

On the double averaged three-body problem

The Hamiltonian of three point masses is averaged over fast variablel and ll (mean anomalies) The problem is non-planar and it is assumed that two of the bodies form a close pair (stellar three-body problem). Only terms up to the order of (a/á)⁴ are taken into account in the Hamiltonian, wherea andá are the corresponding semi-major axes. Employing the method of elimination of the nodes, the problem may be reduced to one degree of freedom. Assuming in addition that the angular momentum of the close binary is much smaller than the angular momentum of the motion of the binary around a third body, we were able to solve the equation for the eccentricity changes in terms of the Jacobian elliptic functions.     

On an analytical solution of the optimum trajectory problem in a gravitational field

Mayer's variational problem for a point with a limited mass flow rate is described by differential equations of the fourteenth order, allowing for a few first integrals. By reducing the equations to closed canonical form, these integrals are analyzed from the viewpoint of finding a possible solution to the problem via quadratures on zero, intermediate, and maximum thrust sections. In addition to confirming well-known cases of total integrability, this approach enabled us to establish that the essential difficulty of the solution of the space problem with intermediate thrust is reduced to finding one integral, and the solution of the problem with maximum thrust requires two integrals in involution. It is shown that these integrals can be applied to find particular solutions.     

Averaging on the motion of a fast revolving body. Application to the stability study of a planetary system

Exploring the global dynamics of a planetary system involves computing integrations for an entire subset of its parameter space. This becomes time-consuming in presence of a planet close to the central star, and in practice this planet will be very often omitted. We derive for this problem an averaged Hamiltonian and the associated equations of motion that allow us to include the average interaction of the fast planet. We demonstrate the application of these equations in the case of the μ Arae system where the ratio of the two fastest periods exceeds 30. In this case, the effect of the inner planet is limited because the planet’s mass is one order of magnitude below the other planetary masses. When the inner planet is massive, considering its averaged interaction with the rest of the system becomes even more crucial.     

Integral surfaces with space-time coordinates,in the gravitational field of a rotating system

A new integration theory is formulated for dynamical systems with two degrees of freedom, in the gravitational field of a rotating system. Four integrals of motion may be determined from complete solutions of a system of three first-order, partial differential equations in three independent variables. The solutions of this system define two integral surfaces with space-time coordinates. These surfaces represent two independent solutions of a second-order kinematic system to which the original fourth-order system has been reduced. An integral curve may be represented as the locus of intersection points of the integral surfaces. The new theory is the theoretical basis for a method of analytic continuation of periodic orbits of the circular restricted problem.     

Zero velocity surfaces for the general planar three-body problem

The angular momentum and the energy integral of the planar three-body problem are used to establish regions of the physical space where motion is allowed to take place. Although forbidden regions exist for both negative and positive values of the energy of the system, the known integrals of the motion always allow for at least one of the three bodies to escape.     

Exact solution of a triaxial gyrostat with one rotor

The problem of the attitude dynamics of a triaxial gyrostat under no external torques and one constant internal rotor, is a three degrees-of-freedom system, although thanks to the existence of integrals of motion it can be reduced to only one degree-of-freedom problem. We introduce coordinates to represent the orbits of constant angular momentum as a flow on a sphere. This representation shows that the problem is equivalent to a quadratic Hamiltonian depending on two parameters. We find the exact solution of the orbits in terms of elliptic functions. By making use of properties of elliptic functions we find the solution at each region of the parametric partition from the solution of one region. We also prove that heteroclinic orbits are planar curves.     

Dynamical evolution of weakly disturbed two-planetary system on cosmogonic time scales: The Sun-Jupiter-Saturn system

In the present paper, we used the Hori-Deprit method to construct the averaged Hamiltonian of the two-planetary problem accurate to the second order of a small parameter, the generating function of the transform, the change of variables formulas, and the right-hand sides of the equations in average elements. The evolution of the two-planet Sun-Jupiter-Saturn system was investigated by numerical integration over 10 billion years. The motion of the planets has an almost periodic character. The eccentricities and inclinations of Jupiter’s and Saturn’s orbits remain small but different from zero. The short-term disturbances remain small over the entire period considered in the study.     

Quasiintegrals of the photogravitational eccentric restricted three-body problem with Poynting–Robertson drag

A restricted three-body problem for a dust particle, in presence of a spherical cometary nucleus in an eccentric (elliptic, parabolic or hyperbolic) orbit about the Sun, is considered. The force of radiation pressure and the Poynting– Robertson effect are taken into account. The differential equations of the particle’s non-inertial spatial motion are investigated both analytically and numerically. With the help of a complex representation, a new single equation of the motion is obtained. Conversion of the equations of motion system into a single equation allows the derivation of simple expressions similar to the integral of energy and integrals of areas. The derived expressions are named quasiintegrals. Relative values of terms of the energy quasiintegral for a smallest, largest, and a mean comet are calculated. We have found that in a number of cases the quasiintegrals are related to the regular integrals of motion, and discuss how the quasiintegrals may be applied to find some significant constraints on the motion of a body of infinitesimal mass.     

New formulation of the two body problem using a continued fractional potential

In the present work, the two body problem using a potential of a continued fractions procedure is reformulated. The equations of motion for two bodies moving under their mutual gravity is constructed. The integrals of motion, angular momentum integral, center of mass integral, total mechanical energy integral are obtained. New orbit equation is obtained. Some special cases are followed directly. Some graphical illustrations are shown. The only included constant of the continued fraction procedure is adjusted so as to represent the so called J ₂ perturbation term of the Earth’s potential.     

Motion near the triangular points in the elliptic restricted problem of three bodies

In this paper the first variational equations of motion about the triangular points in the elliptic restricted problem are investigated by the perturbation theories of Hori and Deprit, which are based on Lie transforms, and by taking the mean equations used by Grebenikov as our upperturbed Hamiltonian system instead of the first variational equations in the circular restricted problem. We are able to remove the explicit dependence of transformed Hamiltonian on the true anomaly by a canonical transformation. The general solution of the equations of motion which are derived from the transformed Hamiltonian including all the constant terms of any order in eccentricity and up to the periodic terms of second order in eccentricity of the primaries is given.