Lie transforms and the Hamiltonization of non-Hamiltonian systems

To develop the perturbation solution of the non-Hamiltonian system of differential equationsy=g(y, t; ), it is sufficient to obtain the perturbation solution of a Hamiltonian system represented by the HamiltonianK=Y·g(y, t; ) which is linear in the adjoint vectorY. This Hamiltonization allows the direct use of the perturbation methods already established for Hamiltonian systems. To demonstrate this fact, a Hamiltonian algorithm developed by this author and based on the Lie-Deprit transform is applied to the Hamiltonized system and is shown to be equivalent to the application of the non-Hamiltonian form of this same algorithm to the original non-Hamiltonian system.     

The constrained normal form algorithm

This paper describes an algorithm which brings a regularizable polynomial perturbation of a three degree of freedom Kepler problem into a normal form which Poisson commutes with the Kepler Hamiltonian. We illustrate the alogrithm with an example: the quadratic Zeeman effect. For other applications of this algorithm see [1],[4], and [5]. The authors have written a program in MAPLE which implements the constrained normal form.     

Perturbations des systèmes différentiels et résultats de géométrie différentielle

We use classical definitions and results of differential geometry in studying properties of transformations depending on a small parameter, acting on differential systems. Hori's and Deprit's algorithms can be defined for these systems. A lemma is given to show these algorithms are equivalent. The so-called property of covariance is merely established. The canonical systems are then considered as associated with Hamiltonian vectorfields on symplectic manifolds. The property that the infinitesimal generator of a canonical transformation is an Hamiltonian vectorfield permits to establish separately the generality of Hori's and Deprit's algorithms. We suggest that the Hamiltonian vectorfield property can be extended to the generators of transformations depending on several parameters.     

Expansion of the Hamiltonian of the Two-Planetary Problem into the Poisson Series in All Elements: Application of the Poisson Series Processor

This paper is the third in a series of articles devoted to one of the basic problems of celestial mechanics: the study of the evolution of solar-type planetary systems. In the previous papers a brief review of the history and current state of the problem was given; the plan of the study was outlined; the Jacobi coordinates and the related osculating elements were introduced; the form of the Poisson expansion of the Hamiltonian in all elements was given; and the expansion coefficients for the Hamiltonian of the two-planetary Sun–Jupiter–Saturn problem were obtained (though with impure accuracy) by a simple algorithm that is reduced to the calculation of multiple integrals of elementary functions. In the present paper the expansion of the Hamiltonian of the two-planetary Sun–Jupiter–Saturn problem into the Poisson series in all elements is constructed with the help of the PSP Poisson series processor, which is capable of required accuracy.     

A note on the algorithm of symplectic integrators

For a Hamiltonian that can be separated into N+1(N\geq 2) integrable parts, four algorithms can be built for a symplectic integrator. This research compares these algorithms for the first and second order integrators. We found that they have similar local truncation errors represented by error Hamiltonian but rather different numerical stability. When the computation of the main part of the Hamiltonian, H ₀, is not expensive, we recommend to use S ^* type algorithm, which cuts the calculation of the H ₀ system into several small time steps as Malhotra(1991) did. As to the order of the N+1 parts in one step calculation, we found that from the large to small would get a slower error accumulation. This revised version was published online in July 2006 with corrections to the Cover Date.     

Note on Mercury’s rotation: the four equilibria of the Hamiltonian model

Mercury is observed in a stable Cassini’s state, close to a 3:2 spin-orbit resonance, and a 1:1 node resonance. This present situation is not the only possible mathematical stable state, as it is shown here through a simple model limited to the second-order in harmonics and where Mercury is considered as a rigid body. In this framework, using a Hamiltonian formalism, four different sets of resonant angles are computed from the differential Hamiltonian equations, and each of them corresponds to four values of the obliquity; thanks to the calculation of the corresponding eigenvalues, their linear stability is analyzed. In this simplified model, two equilibria (one of which corresponding to the present state of Mercury) are stable, one is unstable, and the fourth one is degenerate. This degenerate status disappears with the introduction of the orbit (node and pericenter) precessions. The influence of these precession rates on the proper frequencies of the rotation is also analyzed and quantified, for different planetary models.     

Attitude stability of artificial satellites subject to gravity gradient torque

The stability of the rotational motion of artificial satellites is analyzed considering perturbations due to the gravity gradient torque, using a canonical formulation, and Andoyer’s variables to describe the rotational motion. The stability criteria employed requires the reduction of the Hamiltonian to a normal form around the stable equilibrium points. These points are determined through a numerical study of the Hamilton’s equations of motion and linear study of their stability. Subsequently a canonical linear transformation is used to diagonalize the matrix associated to the linear part of the system resulting in a normalized quadratic Hamiltonian. A semi-analytic process of normalization based on Lie–Hori algorithm is applied to obtain the Hamiltonian normalized up to the fourth order. Lyapunov stability of the equilibrium point is performed using Kovalev and Savchenko’s theorem. This semi-analytical approach was applied considering some data sets of hypothetical satellites, and only a few cases of stable motion were observed. This work can directly be useful for the satellite maintenance under the attitude stability requirements scenario.     

Ptolemaic Transformation in Keplerian Problem

We propose the Ptolemaic transformation: a canonical change of variables reducing the Keplerian motion to the form of a perturbed Hamiltonian problem. As a solution of the unperturbed case, the Ptolemaic variables define an intermediary orbit, accurate up to the first power of eccentricity, like in the kinematic model of Claudius Ptolemy. In order to normalize the perturbed Hamiltonian we modify the recurrent Lie series algorithm of HoriuuMersman. The modified algorithm accounts for the loss of a term's order during the evaluation of a Poisson bracket, and thus can be also applied in resonance problems. The normalized Hamiltonian consists of a single Keplerian term; the mean Ptolemaic variables occur to be trivial, linear functions of the Delaunay actions and angles. The generator of the transformation may serve to expand various functions in Poisson series of eccentricity and mean anomaly.     

A secular theory of coplanar, non-resonant planetary system

We present the secular theory of coplanar N -planet system, in the absence of mean motion resonances between the planets. This theory relies on the averaging of a perturbation to the two-body problem over the mean longitudes. We expand the perturbing Hamiltonian in Taylor series with respect to the ratios of semimajor axes which are considered as small parameters, without direct restrictions on the eccentricities. Next, we average out the resulting series term by term. This is possible thanks to a particular but in fact quite elementary choice of the integration variables. It makes it possible to avoid Fourier expansions of the perturbing Hamiltonian. We derive high-order expansions of the averaged secular Hamiltonian (here, up to the order of 24) with respect to the semimajor axes ratio. The resulting secular theory is a generalization of the octupole theory. The analytical results are compared with the results of numerical (i.e. practically exact) averaging. We estimate the convergence radius of the derived expansions, and we propose a further improvement of the algorithm. As a particular application of the method, we consider the secular dynamics of three-planet coplanar system. We focus on stationary solutions in the HD 37124 planetary system.     

Non-integrability of the dumbbell and point mass problem

This paper discusses a constrained gravitational three-body problem with two of the point masses separated by a massless inflexible rod to form a dumbbell. This problem is a simplification of a problem of a symmetric rigid body and a point mass, and has numerous applications in Celestial Mechanics and Astrodynamics. The non-integrability of this system is proven. This was achieved thanks to an analysis of variational equations along a certain particular solution and an investigation of their differential Galois group. Nowadays this approach is the most effective tool for study integrability of Hamiltonian and non-Hamiltonian systems.     

A method of symplectic integrations with adaptive time-steps for individual Hamiltonians in the planetary <Emphasis Type="Italic">N</Emphasis>-body problem

A new algorithm is developed for long-term integrations of the N-body problem. The method uses symplectic integrations of the Hamiltonian equations of motion for each body. This allows one to employ individual adaptive time-steps in computations. The efficiency of this technique is demonstrated by several tests performed for typical problems of Solar System dynamics.     

On the convergence of an algorithm constructing the normal form for elliptic lower dimensional tori in planetary systems

We give a constructive proof of the existence of elliptic lower dimensional tori in nearly integrable Hamiltonian systems. In particular we adapt the classical Kolmogorov normalization algorithm to the case of planetary systems, for which elliptic tori may be used as replacements of elliptic Keplerian orbits in Lagrange-Laplace theory. With this paper we support with rigorous convergence estimates the semi-analytic work in our previous article (Sansottera et al., Celest Mech Dyn Astron 111:337–361, 2011), where an explicit calculation of an invariant torus for a planar model of the Sun-Jupiter-Saturn-Uranus system has been made. With respect to previous works on the same subject we exploit the characteristic of Lie series giving a precise control of all terms generated by our algorithm. This allows us to slightly relax the non-resonance conditions on the frequencies.     

Extension of Cassini's Laws

We investigate the Cassini's laws which describe the rotational motion in a 1:1 spin-orbit resonance. When this rotational motion follows the conventional Cassini's laws, the figure axis coincides with the angular momentum axis. In this case we underline the differences between the rotational Hamiltonian for a 'slow rotating' body like the Moon and for a 'fast rotating' body like Phobos. Then, we study a more realistic rotational Hamiltonian where the angle J between the figure axis and the angular momentum axis could be different from zero. This Hamiltonian has not been studied before. We have found a new particular solution for this Hamiltonian which could be seen as an extension of the Cassini's laws. In this new solution the angle J is constant, which is not zero, and the precession of the angular momentum plane is equal to the mean motion of the argument of pericenter of the rotating body. This type of rotational motion is only possible when the orbital eccentricity of the rotating body is not zero. This new law enables describing in particular, the Moon mean rotational motion for which the mean value of the angle J is found to be equal to 103.9±0.7 s of arc.     

On one special case of parametric resonance in problems of celestial mechanics

We consider a periodic (in time) linear Hamiltonian system that depends on a small parameter. At a zero value of this parameter, the matrix of the system is constant, has two identical pairs of purely imaginary roots, and is not reducible to diagonal form. Therefore, the unperturbed system is unstable. We propose an algorithm for determining the boundaries of the instability regions for the system at nonzero values of the small parameter. This algorithm was used to analyze the stability of triangular libration points in the elliptical restricted three-body problem and in the stability problem in one special case of stationary rotation of a satellite relative to the center of mass.     

The Hamiltonian transformation in quadratic lie transforms

The algorithm for Hamiltonian transformation in the quadratic perturbation technique of one of the authors admits of various equivalent forms. Using as a criterion the number of inter-term multiplications required for transformation, however, the amount of effort required to obtain the transformed Hamiltonian is not equivalent among these forms. Each is considered in some detail, and general guidelines for the choice of most efficient algorithm to be used in a given problem are provided. Their utility is demonstrated by application to Duffing's equation.     

几类辛方法的数值稳定性研究

主要对一阶隐式Euler辛方法M1、二阶隐式Euler中点辛方法M2、一阶显辛Euler方法M3和二阶leapfrog显辛积分器M4共4种辛方法及一些组合算法进行了通常意义下的线性稳定性分析．针对线性哈密顿系统,理论上找到每个数值方法的稳定区,然后用数值方法检验其正确性．对于哈密顿函数为实对称二次型的情况,为了理论推导便利,特推荐采用相似变换将二次型的矩阵对角化来研究辛方法的线性稳定性．当哈密顿分解为一个主要部分和一个小摄动次要部分且二者皆可积时,无论是线性系统还是非线性系统,这种主次分解与哈密顿具有动势能分解相比,明显扩大了辛方法的稳定步长范围．     

Non-integrability of first order resonances in Hamiltonian systems in three degrees of freedom

The normal forms of the Hamiltonian 1:2:ω resonances to degree three for ω = 1, 3, 4 are studied for integrability. We prove that these systems are non-integrable except for the discrete values of the parameters which are well known. We use the Ziglin–Morales–Ramis method based on the differential Galois theory.     

Periodic solutions with alternating singularities in the collinear four-body problem

This paper gives an analytic proof of the existence of Schubart-like orbit, a periodic orbit with singularities in the symmetric collinear four-body problem. In each period of the Schubart-like orbit, there is a binary collision (BC) between the inner two bodies and a simultaneous binary collision (SBC) of the two clusters on both sides of the origin. The system is regularized and the existence is proved by using a “turning point” technique and a continuity argument on differential equations of the regularized Hamiltonian.     

On the use of the autonomous Birkhoff equations in Lie series perturbation theory

In this article, we present the Lie transformation algorithm for autonomous Birkhoff systems. Here, we are referring to Hamiltonian systems that obey a symplectic structure of the general form. The Birkhoff equations are derived from the linear first-order Pfaff–Birkhoff variational principle, which is more general than the Hamilton principle. The use of 1-form in formulating the equations of motion in dynamics makes the Birkhoff method more universal and flexible. Birkhoff’s equations have a tensorial character, so their form is independent of the coordinate system used. Two examples of normalization in the restricted three-body problem are given to illustrate the application of the algorithm in perturbation theory. The efficiency of this algorithm for problems of asymptotic integration in dynamics is discussed for the case where there is a need to use non-canonical variables in phase space.     

Geometry of Poincaré's Variables and the Secular Planetary Problem

A method for the expansion of the perturbative Hamiltonian in the planetary problem is presented, which allows one to immediately detect the terms vanishing under the averaging process. The method bases itself on a geometrical analysis, through the groups SO(3) and SU(2), of the Poincaré canonical variables or of the similar Laplace variables. As an outcome, one obtains a MAPLE program, which calculates the first averaged terms of the perturbative Hamiltonian. This revised version was published online in July 2006 with corrections to the Cover Date.