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.     

The Expansion of the Hamiltonian of the Planetary Problem into the Poisson Series in All Keplerian Elements (Theory)

The study of the evolution of planetary systems, primarily of the Solar System, is one of the basic problems of celestial mechanics. The stability of motion of giant planets on cosmogonic time scales was established by numerical and analytical methods, but the question about the evolution of orbits of terrestrial planets and arbitrary solar-type planetary systems remained open. This work initiates a series of papers allowing one to advance in solving the problem of the evolution of the solar-type planetary systems on cosmogonic time scales by using powerful analytical tools. In the first paper of this series, we choose the optimum reference system and obtain the Poisson series expansion of the Hamiltonian of the problem in all Keplerian elements. We propose to use the integral representation of the corresponding coefficients or the Poisson processor means instead of conventionally addressing any possible special functions. This approach extremely simplifies the algorithm. The next paper of this series deals with the calculation of the expansion coefficients.     

Unstable modes of Keplerian discs

We derive the softened analogue of the Laplace–Lagrange secular theory of planetary motion, and use it to show that a small fraction of counter-rotating stars is all it takes for a hot Keplerian disc to grow unstable lopsided modes.     

The stability of our solar system

Most investigations of the stability of the solar system have been concerned with the question as to whether the very long term effect of the gravitational attractions of the planets on each other will be to alter the nearly coplanar, nearly circular nature of the orbits in which they move. Analytical investigations in the traditions of Laplace, Lagrange, Poisson and Poincaré strongly indicate stability, though rely on asymptotic expansions with difficult analytical properties. The question is related to the existence of invariant tori, which have been proved to exist in certain motions. Numerical integration experiments have thrown considerable light on possible types of motions, especially in fictitious solar systems in which the planetary masses have been increased to enhance the perturbations, and in testing how critical are stability boundary estimates given by Hill surface type methods.     

Gyldén systems: Rotation of pericenters

A canonical transformation in phase space and a rescaling of time are proposed to reduce a Keplerian system with a time-dependent Gaussian parameter to a perturbed Keplerian system with a constant Gaussian parameter. When the time variation is slow, the perturbation through second order in the reduced problem is conservative, and, as a result, the orbital space of the averaged system is a sphere on which the phase flow causes a differential rotation representing a circulation of the line of apsides. The flow presents two isolated singularities corresponding to circular orbits travelled respectively in the direct and in the retrograde sense, and a degenerate manifold of fixed points corresponding to the collision orbits. Normalization beyond order two does not break the degeneracy. Adiabatic invariants, which are conservative functions, may be computed from the normalized Hamiltonian evaluated here to the fourth order. Nonetheless so high an approximation gives little information because the normalizing Lie transformation depends explicitly on the time through mixed secular-periodic terms. As an application, an estimate is offered for the apsidal rotation that a second order time derivative in the mass of the sun would induce on planetary orbits. This suggests an observational method for determining the latter parameter for the solar wind, but the predicted motions are too slow for the current level of observational precision.     

Gauss's method for secular dynamics, softened

We show that the algorithm proposed by Gauss to compute the secular evolution of gravitationally interacting Keplerian rings extends naturally to softened gravitational interactions. The resulting tool is ideal for the study of the secular dynamical evolution of nearly Keplerian systems such as stellar clusters surrounding black holes in galactic nuclei, cometary clouds or planetesimal discs. We illustrate its accuracy, efficiency and versatility on a variety of configurations. In particular, we examine a secularly unstable system of counterrotating discs, and follow the unfolding and saturation of the instability into a global, uniformly precessing, lopsided  ( m = 1)  mode.     

The present status of periodic orbits

Recent results on periodic orbits are presented and it is shown that the periodic orbits can be used in the study of planetary systems and triple or multiple stellar systems. Triple stellar systems are stable even for close approaches of the three components. Also stable triple systems exist with nearly zero angular momentum. For the planetary systems a global view is obtained from which it is clear which configurations are stable or unstable and also what factors affect the stability. Also, the relation between resonance and instability is studied by making use of periodic orbits.     

The 1/1 resonance in extrasolar systems

We present families of symmetric and asymmetric periodic orbits at the 1/1 resonance, for a planetary system consisting of a star and two small bodies, in comparison to the star, moving in the same plane under their mutual gravitational attraction. The stable 1/1 resonant periodic orbits belong to a family which has a planetary branch, with the two planets moving in nearly Keplerian orbits with non zero eccentricities and a satellite branch, where the gravitational interaction between the two planets dominates the attraction from the star and the two planets form a close binary which revolves around the star. The stability regions around periodic orbits along the family are studied. Next, we study the dynamical evolution in time of a planetary system with two planets which is initially trapped in a stable 1/1 resonant periodic motion, when a drag force is included in the system. We prove that if we start with a 1/1 resonant planetary system with large eccentricities, the system migrates, due to the drag force, along the family of periodic orbits and is finally trapped in a satellite orbit. This, in principle, provides a mechanism for the generation of a satellite system: we start with a planetary system and the final stage is a system where the two small bodies form a close binary whose center of mass revolves around the star.     

Noncanonical Parametrization of Poisson Brackets in Celestial Mechanics

The formulas for the Poisson bracket of a perturbed two-body problem and a perturbed planetary problem are found in different systems of Keplerian elements. As with canonical parametrization, the Poisson bracket is equal to a linear combination of partial brackets, but it contains coefficients depending on semimajor axis, eccentricity, and inclination. A simple relation between the Poisson brackets and matrices of coefficients of Lagrange-type equations determining the variations of osculating elements is derived. The Poisson bracket of D'Alembertian functions is proved to be a D'Alembertian one by itself.     

Canonical planetary perturbation equations for velocity-dependent forces,and the Lense-Thirring precession

A form of planetary perturbation theory based on canonical equations of motion, rather than on the use of osculating orbital elements, is developed and applied to several problems of interest. It is proved that, with appropriately selected initial conditions on the orbital elements, the two forms of perturbation theory give rise to identical predictions for the observable coordinates and velocities, while the orbital elements themselves may be strikingly different. Differences between the canonical form of perturbation theory and the classical Lagrange planetary perturbation equations are discussed. The canonical form of perturbation theory in some cases has advantages when the perturbing forces are velocity-dependent, but the two forms of perturbation theory are equivalent if the perturbing forces are dependent only on position and not on velocity. The canonical form of the planetary perturbation equations are derived and applied to the Lense Thirring precession of a test body in a Keplerian orbit around a rotating mass source.     

The mid-term and long-term solar quasi-periodic cycles and the possible relationship with planetary motions

This work investigates the solar quasi-periodic cycles with multi-timescales and the possible relationships with planetary motions. The solar cycles are derived from long-term observations of the relative sunspot number and microwave emission at frequency of 2.80 GHz. A series of solar quasi-periodic cycles with multi-timescales are registered. These cycles can be classified into three classes: (1) the strong PLC (PLC is defined as the solar cycle with a period very close to the ones of some planetary motions, named as planetary-like cycle) which is related strongly with planetary motions, including nine periodic modes with relatively short period (P<12 yr), and related to the motions of the inner planets and of Jupiter; (2) the weak PLC, which is related weakly to planetary motions, including two periodic modes with relatively long period (P>12 yr), and possibly related to the motions of outer planets; (3) the non-PLC, for which so far there has been found no clear evidence to show the relationship with any planetary motions. Among the planets, Jupiter plays a key role in most periodic modes due to its sidereal motion or spring tidal motions associated with other planets. Among planetary motions, the spring tidal motion of the inner planets and of Jupiter dominates the formation of most PLCs. The relationships between multi-timescale solar periodic modes and the planetary motions will help us to understand the essential nature and prediction of solar activities.     

Regularizing Time Transformations in Symplectic and Composite Integration

We consider numerical integration of nearly integrable Hamiltonian systems. The emphasis is on perturbed Keplerian motion, such as certain cases of the problem of two fixed centres and the restricted three-body problem. We show that the presently known methods have useful generalizations which are explicit and have a variable physical timestep which adjusts to both the central and perturbing potentials. These methods make it possible to compute accurately fairly close encounters. In some cases we suggest the use of composite (instead of symplectic) alternatives which typically seem to have equally good energy conservation properties.This revised version was published online in October 2005 with corrections to the Cover Date.     

GJ 581 update: Additional evidence for a Super‐Earth in the habitable zone

We present an analysis of the significantly expanded HARPS 2011 radial velocity data set for GJ 581 that was presented by Forveille et al. (2011). Our analysis reaches substantially different conclusions regarding the evidence for a Super‐Earth‐mass planet in the star's Habitable Zone. We were able to reproduce their reported χ²_ν and RMS values only after removing some outliers from their models and refitting the trimmed down RV set. A suite of 4000 N‐body simulations of their Keplerian model all resulted in unstable systems and revealed that their reported 3.6σ detection of e = 0.32 for the eccentricity of GJ 581e is manifestly incompatible with the system's dynamical stability. Furthermore, their Keplerian model, when integrated only over the time baseline of the observations, significantly increases the χ²_ν and demonstrates the need for including non‐Keplerian orbital precession when modeling this system. We find that a four‐planet model with all of the planets on circular or nearly circular orbits provides both an excellent self‐consistent fit to their RV data and also results in a very stable configuration. The periodogram of the residuals to a 4‐planet all‐circular‐orbit model reveals significant peaks that suggest one or more additional planets in this system. We conclude that the present 240‐point HARPS data set, when analyzed in its entirety, and modeled with fully self‐consistent stable orbits, by and of itself does offer significant support for a fifth signal in the data with a period near 32 days. This signal has a false alarm probability of <4% and is consistent with a planet of minimum mass 2.2 M_⊕, orbiting squarely in the star's habitable zone at 0.13 AU, where liquid water on planetary surfaces is a distinct possibility (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the relation between resonance and instability in planetary systems

The factors which affect the linear stability of a periodic planetary orbit in the plane are studied. It is proved that planetary systems with two planets describing nearly circular orbits in the same direction are linearly stable and no perturbation exists which destroys the stability, unless a resonance of the form 1/3, 3/5, 5/7, ... among the orbits of the planets occurs. This latter resonant case is always unstable. Retrograde motion is always linearly stable. Planetary systems with three or more planets in nearly circular orbits in the same direction are proved to be unstable, in the sense that a Hamiltonian perturbation always exists which destroys the stability. The generation of instability in the case of three or more planets is not only due to the existence of resonances, as in the case of two planets, but also to the nonexistence of integrals of motion, apart from the energy and angular momentum integrals. It is also proved that planetary systems with nearly elliptic orbits of the planets are unstable.     

系外行星系统掩星周期变化研究进展

在掩星法发现的系外行星系统中,如果存在其他未知的伴星绕同一颗恒星运动,掩星行星由于受到伴星引力的影响,运动轨道将发生变化,轨道周期不再是常数,而是变化的。利用这种变化探测掩星系统中的其他行星,已成为一种新的方法。主要介绍了未知行星与掩星行星之间的引力作用引起的掩星周期变化效应,以及掩星周期变化法探测系外行星的理论和研究进展状况,最后简要讨论了几种影响掩星周期变化的其他因素:共轨行星、卫星、潮汐效应、相对论效应及恒星的引力四极矩等。     

The 1/1 resonance in extrasolar planetary systems

We study orbits of planetary systems with two planets, for planar motion, at the 1/1 resonance. This means that the semimajor axes of the two planets are almost equal, but the eccentricities and the position of each planet on its orbit, at a certain epoch, take different values. We consider the general case of different planetary masses and, as a special case, we consider equal planetary masses. We start with the exact resonance, which we define as the 1/1 resonant periodic motion, in a rotating frame, and study the topology of the phase space and the long term evolution of the system in the vicinity of the exact resonance, by rotating the orbit of the outer planet, which implies that the resonance and the eccentricities are not affected, but the symmetry is destroyed. There exist, for each mass ratio of the planets, two families of symmetric periodic orbits, which differ in phase only. One is stable and the other is unstable. In the stable family the planetary orbits are in antialignment and in the unstable family the planetary orbits are in alignment. Along the stable resonant family there is a smooth transition from planetary orbits of the two planets, revolving around the Sun in eccentric orbits, to a close binary of the two planets, whose center of mass revolves around the Sun. Along the unstable family we start with a collinear Euler–Moulton central configuration solution and end to a planetary system where one planet has a circular orbit and the other a Keplerian rectilinear orbit, with unit eccentricity. It is conjectured that due to a migration process it could be possible to start with a 1/1 resonant periodic orbit of the planetary type and end up to a satellite-type orbit, or vice versa, moving along the stable family of periodic orbits.     

On the gravitational instability in thin gaseous Kepler disks

The idea that the Titius-Bode law reflects an unstable mode of a self-gravitationl instability in very thin Keplerian disks makes a careful discussion of the Poisson equation especially necessary. Due to the planetary distances in the solar system (δr/r ⋍ 0·5) the well-known short-wave approximation is not appropriate for definite asertions. We will here use a simple series expansion of the relation between the radial and vertical wave numbers of the disturbances which is additionally valid for medium-scale and non-zonal modes. The numerical solution of the dispersion relation reveals an extra unstable branch for wave-lengths of rings and spirals two orders of magnitudes larger than those already known. Though we are not yet able to consider modes long enough for application to the planetary system, we feel the existence of the medium-wave instability (δr/r ⋍ 0·1) to be a serious challenge for a better, i.e. non-local theory.     

Generation of a secular system in the analytical theory for the motion of the moon

In order to generate an analytical theory of the motion of the Moon by considering planetary perturbations, a procedure of general planetary theory (GPT) is used. In this case, the Moon is considered as an addition planet to the eight principal planets. Therefore, according to the GPT procedure, the theory of the Moon’s orbital motion can be presented in the form of series with respect to the evolution of eccentric and oblique variables with quasi-periodic coefficients, which are the functions of mean longitudes for principal planets and the Moon. The relationship between evolution variables and the time is determined by a trigonometric solution for the independent secular system that describes the secular motion of a perigee and the Moon node by considering secular planetary inequalities. Principal planetary coordinates required for generating the theory of the motion of the Moon includes only Keplerian terms, the intermediate orbit, and the linear theory with respect to eccentricities and inclinations in the first order relative to the masses. All analytical calculations are performed by means of the specialized echeloned Poisson Series Processor EPSP.     

On the differential geometry of trajectory deviation with particular emphasis on Keplerian motion

The Pontryagin/Lawden scalar product of the deviation phase vector and the adjoint phase vector may be identified with Lagrange's reciprocal formula for two variant motions if the acceleration field is conservative. Hence for two slightly different trajectories with common end-points, the terminal velocity differences must have equal scalar product with Lawden's primer vectors. The final velocity difference is orthogonal to the constant time of flight locus for isoenergetic motions from a common initial point. A Keplerian trajectory behaves like a rigid curve as far as radial deviation is concerned in the case of a small change of direction of the initial velocity vector.