Third and fourth order resonances in Hamiltonian systems

The stability of the origin of an autonomous Hamiltonian system is investigated when the system possesses a third or fourth-order resonance.H ₂, the quadratic part ofH isH ₂=ⁿ _{i=1 i} J _i and the resonance condition is ⁿ _i=1 k _{i i} where thek 0,i = 1, 2, ...,n are the natural or fundamental frequencies. It is shown that the only case in which the origin can be unstable is ifk _i0,i=1,2,...,n. The condition for instability is then given in terms of the coefficients of the higher order terms in the Hamiltonian. The transfer of energy between modes is also investigated when a near-resonant condition exists.     

On the stability of strongly non-homogeneous self-gravitating equilibria

The stability analysis of several stronglynon-homogeneous, self-gravitating, one-dimensionalunstable equilibrium systems is carried out with the help of numerical techniques. The evolution of the perturbed unstable equilibria is studied by following the motion of the boundary curves of water bag configurations defining the systems.It is found that initial perturbations drive the unstable equilibrium states out of equilibrium at rates depending on the typical scale length of the perturbations : the instability rates increase with .     

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.     

On the chaotic diffusion in multidimensional Hamiltonian systems

We present numerical evidence that diffusion in the herein studied multidimensional near-integrable Hamiltonian systems departs from a normal process, at least for realistic timescales. Therefore, the derivation of a diffusion coefficient from a linear fit on the variance evolution of the unperturbed integrals fails. We review some topics on diffusion in the Arnold Hamiltonian and yield numerical and theoretical arguments to show that in the examples we considered, a standard coefficient would not provide a good estimation of the speed of diffusion. However, numerical experiments concerning diffusion would provide reliable information about the stability of the motion within chaotic regions of the phase space. In this direction, we present an extension of previous results concerning the dynamical structure of the Laplace resonance in Gliese-876 planetary system considering variations of the orbital parameters accordingly to the error introduced by the radial velocity determination. We found that a slight variation of the eccentricity of planet c would destabilize the inner region of the resonance that, though chaotic, shows stable when adopting the best fit values for the parameters.     

A proof of Nekhoroshev's theorem for the stability times in nearly integrable Hamiltonian systems

In the present paper we give a proof of Nekhoroshev's theorem, which is concerned with an exponential estimate for the stability times in nearly integrable Hamiltonian systems. At variance with the already published proof, which refers to the case of an unperturbed Hamiltonian having the generic property of steepness, we consider here the particular case of a convex unperturbed Hamiltonian. The corresponding simplification in the proof might be convenient for an introduction to the subject.     

On the stability of hierarchical four-body systems

Generalized Jacobian coordinates can be used to decompose anN-body dynamical system intoN-1 2-body systems coupled by perturbations. Hierarchical stability is defined as the property of preserving the hierarchical arrangement of these 2-body subsystems in such a way that orbit crossing is avoided. ForN=3 hierarchical stability can be ensured for an arbitrary span of time depending on the integralz=c ² h (angular momentum squared times energy): if it is smaller than a critical value, defined by theL ₂ collinear equilibrium configuration, then the three possible hierarchical arrangements correspond to three disconnected subsets of the invariant manifold in the phase space (and in the configuration space as well; see Milani and Nobili, 1983a). The same definitions can be extended, with the Jacobian formalism, to an arbitrary hierarchical arrangement ofN≥4 bodies, and the main confinement condition, the Easton inequality, can also be extended but it no longer provides separate regions of trapped motion, whatever is the value ofz for the wholeN-body system,N≥4. However, thez criterion of hierarchical stability applies to every 3-body subsystem, whosez ‘integral’ will of course vary in time because of the perturbations from the other bodies. In theN=4 case we decompose the system into two 3-body subsystems whosec ² h ‘integrals’,z ₂₃ andz ₃₄, att=0 are assumed to be smaller than the corresponding critical values \(\tilde z_{23} \) and \(\tilde z_{34} \) , so that both the subsystems are initially hierarchically stable. Then the hierarchical arrangement of the 4 bodies cannot be broken until eitherz ₂₃ orz ₃₄ is changed by an amount \(\tilde z_{ij} - z_{ij} \left( 0 \right)\) ; that is the whole system is hierarchically stable for a time spain not shorter than the minimum between \(\Delta t_{23} = {{\left( {\tilde z_{23} - z_{23} \left( 0 \right)} \right)} \mathord{\left/ {\vphantom {{\left( {\tilde z_{23} - z_{23} \left( 0 \right)} \right)} {\dot z_{23} }}} \right. \kern-0em} {\dot z_{23} }}\) and \(\Delta t_{34} = {{\left( {\tilde z_{34} - z_{34} \left( 0 \right)} \right)} \mathord{\left/ {\vphantom {{\left( {\tilde z_{34} - z_{34} \left( 0 \right)} \right)} {\dot z_{34} }}} \right. \kern-0em} {\dot z_{34} }}\) . To estimate how long is this stability time, two main steps are required. First the perturbing potentials have to be developed in series; the relevant small parameters are some combinations of mass ratios and length ratios, the? _ij of Roy and Walker. When an appropriate perturbation theory is based on the? _ij , the asymptotic expansions are much more rapidly decreasing than the usual expansions in powers of the mass ratios (as in the classical Lagrange perturbation theory) and can be extended also to cases such as lunar theory or double binaries. The second step is the computation of the time derivatives \(\dot z_{ij} \) (we limit ourselves to the planar case). To assess the long term behaviour of the system, we can neglect the short-periodic perturbations and discuss only the long-periodic and the secular perturbations. By using a Poisson bracket formalism, a generalization of Lagrange theorem for semimajor axes and a generalization of the classical first order theories for eccentricities and pericenters, we prove that thez _ij do not undergo any secular perturbation, because of the interaction with the other subsystem, at the first order in the? _ik . After the long-periodic perturbations have been accounted for, and apart from the small divisors problems that could arise both from ordinary and secular resonances, only the second order terms have to be considered in the computation of Δt ₂₃, Δt ₃₄. A full second order perturbative theory is beyond the scope of this paper; however an order-of-magnitude lower estimate of the Δt _ij can be obtained with the very pessimistic assumption that essentially all the second order terms affect in a secular way thez _ij . The same method could be applied also toN≥5 body systems. Since almost everyN-body system existing in nature is strongly hierarchical, the product of two? _ij is very small for almost all the real astronomical problems. As an example, the hierarchical stability of the 4-body system Sun, Mercury, Venus, and Jupiter is investigated; this system turns out to be stable for at least 110 million years. Although this hierarchical stability time is ～10 times less than the real age of the Solar System, taking into account that many pessimistic assumptions have been done we can conclude that the stability of the Solar System is no more a forbidden problem for Celestial Mechanics.     

On the stability of circumbinary planetary systems

The dynamics of circumbinary planetary systems (the systems in which the planets orbit a central binary) with a small binary mass ratio discovered to date is considered. The domains of chaotic motion have been revealed in the “pericentric distance–eccentricity” plane of initial conditions for the planetary orbits through numerical experiments. Based on an analytical criterion for the chaoticity of planetary orbits in binary star systems, we have constructed theoretical curves that describe the global boundary of the chaotic zone around the central binary for each of the systems. In addition, based on Mardling’s theory describing the separate resonance “teeth” (corresponding to integer resonances between the orbital periods of a planet and the binary), we have constructed the local boundaries of chaos. Both theoretical models are shown to describe adequately the boundaries of chaos on the numerically constructed stability diagrams, suggesting that these theories are efficient in providing analytical criteria for the chaoticity of planetary orbits.     

Spectral stability of relative equilibria

The spectral stability of synchronous circular orbits in a rotating conservative force field is treated using a recently developed Hamiltonian method. A complete set of necessary and sufficient conditions for spectral stability is derived in spherical geometry. The resulting theory provides a general unified framework that encompasses a wide class of relative equilibria, including the circular restricted three-body problem and synchronous satellite motion about an aspherical planet. In the latter case we find an interesting class of stable nonequatorial circular orbits. A new and simplified treatment of the stability of the Lagrange points is given for the restricted three-body problem.     

On the dynamics in the asteroids belt. Part II: Detailed study of the main resonances

The general theory exposed in the first part of this paper is applied to the following resonances with Jupiter's motion : 3/2, 2/1, 5/2, 3/1, 7/2, 4/1; these are the most relevant resonances for the asteroids. The whole analysis is performed in the framework of the spatial problem of three bodies, both in the circular and in the elliptic case. The results are also compared with the observed distribution of the asteroids.     

Effects of resonances on the stability of retrograde satellites

The stability evolution of family f of the planar circular restricted three-body problem in the Earth–Moon case is explored numerically using the Poincaré surface of section. It is shown that third order resonances are the main cause of the reduction of the stability region of retrograde satellites. Several branches of family f are also computed and these are seen by the configuration of their family characteristics to roughly determine the stability region. Previous results on smaller mass ratios of primaries are thus extended to the Earth–Moon system.     

On the stability of dust orbits in mean-motion resonances perturbed by from an interstellar wind

Circumstellar dust particles can be captured in a mean-motion resonance (MMR) with a planet and simultaneously be affected by non-gravitational effects. It is possible to describe the secular variations of a particle orbit in the MMR analytically using averaged resonant equations. We derive the averaged resonant equations from the equations of motion in near-canonical form. The secular variations of the particle orbit depending on the orientation of the orbit in space are taken into account. The averaged resonant equations can be derived/confirmed also from Lagrange’s planetary equations. We apply the derived theory to the case when the non-gravitational effects are the Poynting–Robertson effect, the radial stellar wind, and an interstellar wind. The analytical and numerical results obtained are in excellent agreement. We found that the types of orbits correspond to libration centers of the conservative problem. The averaged resonant equations can lead to a system of equations which holds for stationary points in a subset of resonant variables. Using this system we show analytically that for the considered non-gravitational effects, all stationary points should correspond to orbits which are stationary in interplanetary space after an averaging over a synodic period. In an exact resonance, the stationary orbits are stable. The stability is achieved by a periodic repetition of the evolution during the synodic period. Numerical solutions of this system show that there are no stationary orbits for either the exact or non-exact resonances.     

Normalization algorithms of Hamiltonian near an equilibrium point

Algorithms of linear and nonlinear normalization of a Hamiltonian system near an equilibrium point are described. They were applied to determination of stability libration point in photogravitational restricted problem and to refinement of periodic orbits in the restricted three body problem.     

On the stability of Earth-like planets in multi-planet systems

We present a continuation of our numerical study on planetary systems with similar characteristics to the Solar System. This time we examine the influence of three giant planets on the motion of terrestrial-like planets in the habitable zone (HZ). Using the Jupiter–Saturn–Uranus configuration we create similar fictitious systems by varying Saturn’s semi-major axis from 8 to 11 AU and increasing its mass by factors of 2–30. The analysis of the different systems shows the following interesting results: (i) Using the masses of the Solar System for the three giant planets, our study indicates a maximum eccentricity (max-e) of nearly 0.3 for a test-planet placed at the position of Venus. Such a high eccentricity was already found in our previous study of Jupiter–Saturn systems. Perturbations associated with the secular frequency g ₅ are again responsible for this high eccentricity. (ii) An increase of the Saturn-mass causes stronger perturbations around the position of the Earth and in the outer HZ. The latter is certainly due to gravitational interaction between Saturn and Uranus. (iii) The Saturn-mass increased by a factor 5 or higher indicates high eccentricities for a test-planet placed at the position of Mars. So that a crossing of the Earth’ orbit might occur in some cases. Furthermore, we present the maximum eccentricity of a test-planet placed in the Earth’ orbit for all positions (from 8 to 11 AU) and masses (increased up to a factor of 30) of Saturn. It can be seen that already a double-mass Saturn moving in its actual orbit causes an increase of the eccentricity up to 0.2 of a test-planet placed at Earth’s position. A more massive Saturn orbiting the Sun outside the 5:2 mean motion resonance (a _S  ≥9.7 AU) increases the eccentricity of a test-planet up to 0.4.     

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.     

On the global stability of single‐planet systems

In this paper, we estimate the global stability properties of single‐planet systems by using a catalogue of stability maps. The data of the catalogue were used to generate probability values on the mass parameter–eccentricity plane for the occurrence of stable orbits. We showed that the probability data can be well approximated by a second order surface. Using the resulted formula the likelihood of finding Earth‐like planets in single‐planet systems can be easily estimated. As an example, we derived estimations for four known exoplanetary systems. Our formula can be useful in selecting target stars for future space missions. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the stability of close binaries in hierarchical three-body systems

The Hill-type stability (cf. closure of the zero-velocity curves in the circular restricted three-body problem) of general hierarchical three-body systems is examined analytically in the case where the total mass of the binary is small in comparison to the mass of the external body (e.g. systems of the type Planet-Satellite-Sun, Planet-Planet-Star, etc.). This is compared with results derived by Szebehely, Markellos and Roy in the Planet-Satellite-Sun case of the circular restricted three-body problem. It is demonstrated how the Hill-type stability is affected by the sense of revolution of the binary, i.e. corotational or contrarotational, and the mass ratio within the binary. The effect of the difference in longitudes of the bodies in their orbits is also examined.     

On the Schumann resonances on Mars

A conductive ionosphere and a totally non-conductive layer of the atmosphere close to the surface of the planet form a quasispherical concentric resonator. This provides in principle for the possibility of the existence of global resonances of an electromagnetic field generated by thunderstorm activity or by hydromagnetic waves excited in an upper ionosphere and transformed into ordinary electromagnetic waves while penetrating the resonator. We have obtained an estimate of resonance frequencies of a Martian resonator: ƒ₁ 13–14 Hz,ƒ₂ 24–26 Hz, ƒ₃ 35–38 Hz, etc. for two essentially different models of electron density distribution in the low ionosphere of Mars. The corresponding estimated quality values are low: Q_n 2–4. A relatively wide range of the quality variation depending on a model of averaged altitudinal electron density distribution in the low ionosphere of Mars yields the criterion for an adequate model.     

On the stability of small quasiperiodic motions in the hamiltonian systems

Orbital stability of quasiperiodic motions in the many dimensional autonomic hamiltonian systems is considered. Studied motions are supposed to be not far from equilibrium, the number of their basic frequencies may be not equal to the number of degrees of freedom, and the procedure of their construction is supposed to be converged. The stability problem is solved in the strict nonlinear mode.Obtained results are used in the stability investigation of small plane motions near the lagrangian solutions of the three-dimensional circular restricted three-body problem. The values of parameters for which the plane motions are unstable have been found.

---

. , , . . , . , .

     

Secondary resonances and the boundary of effective stability of Trojan motions

One of the most interesting features in the libration domain of co-orbital motions is the existence of secondary resonances. For some combinations of physical parameters, these resonances occupy a large fraction of the domain of stability and rule the dynamics within the stable tadpole region. In this work, we present an application of a recently introduced ‘basic Hamiltonian model’ \(H_\mathrm{b}\) for Trojan dynamics (Páez and Efthymiopoulos in Celest Mech Dyn Astron 121(2):139, 2015; Páez et al. in Celest Mech Dyn Astron 126:519, 2016): we show that the inner border of the secondary resonance of lowermost order, as defined by \(H_\mathrm{b}\), provides a good estimation of the region in phase space for which the orbits remain regular regardless of the orbital parameters of the system. The computation of this boundary is straightforward by combining a resonant normal form calculation in conjunction with an ‘asymmetric expansion’ of the Hamiltonian around the libration points, which speeds up convergence. Applications to the determination of the effective stability domain for exoplanetary Trojans (planet-sized objects or asteroids) which may accompany giant exoplanets are discussed.