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.     

Third-order Uranus-Neptune planetary theory

In this part we calculate the secular and critical terms arising from the indirect part of the classical planetary Hamiltonian for Uranus and Neptune. We neglect in our expansions powers higher than the second in the eccentricity-inclination. Our required results, are expressed in terms of Poincaré variables.     

On the expansion of the secular part of the perturbing function of mutual attraction in the satellite systems of planets

We propose a special representation for the secular part of the perturbing function describing the mutual attraction of satellites. In contrast to the known representations, it has a single analytical form for any ratio between the semimajor axes of the perturbed and perturbing satellites. The resulting expression is a partial sum of a power series with respect to the small eccentricities and planet-equatorial inclinations of the satellites’ orbits. This sum includes terms up to and including the fourth degree with respect to these small parameters. The proposed expansion is compared with one of the known expansions for the secular part of the perturbing function.     

A semi-analytical model of the Galilean satellites’ dynamics

The Galilean satellites’ dynamics has been studied extensively during the last century. In the past it was common to use analytical expansions in order to get simple models to integrate, but with the new generation of computers it became prevalent the numerical integration of very sophisticated and almost complete equations of motion. In this article we aim to describe the resonant and secular motion of the Galilean satellites through a Hamiltonian, depending on the slow angles only, obtained with an analytical expansion of the perturbing functions and an averaging operation. In order to have a model as near as possible to the actual dynamics, we added perturbations and we considered terms that in similar studies of the past were neglected, such as the terms involving the inclinations and the Sun’s perturbation. Moreover, we added the tidal dissipation into the equations, in order to investigate how well the model captures the evolution of the system.     

Third order uranus-neptune planetary theory

We calculate in this paper the secular and critical terms arising from the principal part of the classical planetary Hamiltonian. This is the first step to establish a third order canonical planetary theory of Uranus-Neptune through the Hori-Lie technique. We truncate our expansions at the second degree of eccentricity-inclination. Our planetary theory is expressed in terms of the canonical variables of H. Poincaré.     

A second order Jupiter-Saturn planetary theory

We present a second order secular Jupiter-Saturn planetary theory through Poincaré canonical variables, von Zeipel's method and Jacobi-Radau referential. We neglect in our expansions terms of power higher than the fourth with respect to eccentricities and sines of inclinations. We assume that the disturbing function is composed of secular and critical terms only. We shall deriveF _2si and writeF _2s in terms of Poincaré canonical variables in Part II of this problem.     

Nodal regression of the quadrantid meteor stream: An analytic approach

The mean orbit of the Quadrantid meteor stream has a high eccentricity and inclination with an aphelion close to the orbit of Jupiter. The nodal regression rate, a quantity which has been well determined from observations, cannot be calculated with sufficient accuracy using standard low-order expansions of the disturbing function. By using a high-order expansion of the disturbing function we show how the behavior of the longitude of ascending node of the Quadrantid stream is a result of both secular and resonant effects. Our analysis illustrates how the proximity of the stream's orbit to the 2: 1 commensurability with Jupiter dominates the short-term variations in orbital elements.     

The secular variations in the restricted problem's sixth-order theory and stability: An explicit numerical approach

In this paper, we develop and implement an explicit numerical technique for studying equilibrium solutions concerning the secular variations in the restricted problem of three bodies, and their stability. In this implementation, we employ a sixth-order theory for the secular terms.     

Partial Reduction in the N-Body Planetary Problem using the Angular Momentum Integral

We present a new set of variables for the reduction of the planetary n-body problem, associated to the angular momentum integral, which can be of any use for perturbation theory. The construction of these variables is performed in two steps. A first reduction, called partial is based only on the fixed direction of the angular momentum. The reduction can then be completed using the norm of the angular momentum. In fact, the partial reduction presents many advantages. In particular, we keep some symmetries in the equations of motion (d'Alembert relations). Moreover, in the reduced secular system, we can construct a Birkhoff normal form at any order. Finally, the topology of this problem remains the same as for the non-reduced system, contrarily to Jacobi's reduction where a singularity is present for zero inclinations. For three bodies, these reductions can be done in a very simple way in Poincaré's rectangular variables. In the general n-body case, the reduction can be performed up to a fixed degree in eccentricities and inclinations, using computer algebra expansions. As an example, we provide the truncated expressions for the change of variable in the 4-body case, obtained using the computer algebra system TRIP.     

The problem of small divisors in planetary motion

Small divisors caused by certain linear combinations of frequencies appear in all analytical planetary theories. With the exception of the deep resonance between Neptune and Pluto, they can be removed at the expense of introducing secular and mixed secular terms, limiting the domain in which the solution is valid. Because of them classical solutions are known not to converge uniformly; Poincaré referred to them as asymptotic. The KAM theory shows that if one is far enough from exact commensurability and has small enough planetary masses, expansions exist which will converge to quasi-periodic orbits. Solutions showing very small divisors are excluded from this region of convergence. The question of whether they are intrinsic to the problem or are just manifestations of the method of solution is not settled. Problems with a single commensurabily that can be isolated from the rest of the Hamiltonian may have solutions with no small divisors. The problem of two or more commensurabilities remains unsolved.     

Secular dynamics of the three-body problem: application to the υ Andromedae planetary system

The discovery of extra-solar planetary systems with multiple planets in highly eccentric orbits (∼0.1-0.6), in contrast with our own Solar System, makes classical secular perturbation analysis very limited. In this paper, we use a semi-numerical approach to study the secular behavior of a system composed of a central star and two massive planets in co-planar orbits. We show that the secular dynamics of this system can be described using only two parameters, the ratios of the semi-major axes and the planetary masses. The main dynamical features of the system are presented in geometrical pictures that allows us to investigate a large domain of the phase space of this three-body problem without time-expensive numerical integrations of the equations of motion, and without any restriction on the magnitude of the planetary eccentricities. The topology of the phase space is also investigated in detail by means of spectral map techniques, which allow us to detect the separatrix of a non-linear secular apsidal resonance. Finally, the qualitative study is supplemented by direct numerical integrations. Three different regimes of secular motion with respect to the secular angle Δ? are possible: they are circulation, oscillation (around 0° and 180°), and high eccentricity libration in a non-linear secular resonance. The first two regimes are a continuous extension of the classical linear secular perturbation theory; the last is a new feature, hitherto unknown, in the secular dynamics of the three-body problem. We apply the analysis to the case of the two outer planets in the υ Andromedae system, and obtain its periodic and ordinary orbits, the general structure of its secular phase space, and the boundaries of its secular stability; we find that this system is secularly stable over a large domain of eccentricities. Applying this analysis to a wide range of planetary mass and semi-major axis ratios (centered about the υ Andromedae parameters), we find that apsidal oscillation dominates the secular phase space of the three-body coplanar system, and that the non-linear secular resonance is also a common feature.     

Polynomial expansions of single-mode motions around equilibrium points in the circular restricted three-body problem

In this work, the single-mode motions around the collinear and triangular libration points in the circular restricted three-body problem are studied. To describe these motions, we adopt an invariant manifold approach, which states that a suitable pair of independent variables are taken as modal coordinates and the remaining state variables are expressed as polynomial series of them. Based on the invariant manifold approach, the general procedure on constructing polynomial expansions up to a certain order is outlined. Taking the Earth–Moon system as the example dynamical model, we construct the polynomial expansions up to the tenth order for the single-mode motions around collinear libration points, and up to order eight and six for the planar and vertical-periodic motions around triangular libration point, respectively. The application of the polynomial expansions constructed lies in that they can be used to determine the initial states for the single-mode motions around equilibrium points. To check the validity, the accuracy of initial states determined by the polynomial expansions is evaluated.     

Secular perturbations of asteroids in secular resonance

When the precessional rate of the orbital plane of an asteroid is nearly equal to that of Jupiter, the orbital inclination of the asteroid changes quite largely due to this near equality of their precessional rates, which is called a secular resonance. In the vicinity of the exact resonance the difference of their longitudes of nodes librates with quite a long period of order of 1×10⁶ yr. In this paper we treat this secular resonance by a method of semianalytical secular perturbations with use of numerical averaging for both non-resonant and resonant asteroids and show that the results by the semi-analytical treatment agrees qualitatively with those obtained by direct numerical integrations of asteroid's orbits.     

Secular Effects on Orbits of Binary Stars Induced by Temporal Variation of Gravitational Constant (Case for Elliptical Orbit)

In this paper we investigate the influence of a varying gravitation constant on the orbits of celestial bodies. Regarding the eccentric anomaly as an independent variable, we find the solutions to the perturbed equations of motion. In the first order solutions, we find the secular and periodic variations in semi-major axis. For the other orbital elements only periodic variations exhibit. However in the second order solutions, the longitude of periastron and the mean longitude have secular terms. Applying the calculations to six selected binaries, we give the numerical estimations of the variations of orbits. These results are then carefully compared and discussed.     

Improved equations for eccentricity generation in hierarchical triple systems

In a series of papers, we developed a technique for estimating the inner eccentricity in hierarchical triple systems, with the inner orbit being initially circular. However, for certain combinations of the masses and the orbital elements, the secular part of the solution failed. In this paper, we derive a new solution for the secular part of the inner eccentricity, which corrects the previous weakness. The derivation applies to hierarchical triple systems with coplanar and initially circular orbits. The new formula is tested numerically by integrating the full equations of motion for systems with mass ratios from 10⁻³ to 10³. We also present more numerical results for short-term eccentricity evolution, in order to get a better picture of the behaviour of the inner eccentricity.     

The dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

Planetary elements for 10 000 000 years

In 1950 Brouwer and van Woerkom published a secular theory of the variations of the planetary elements in analytical form. In the present paper we provide a graphical representation of this theory in the form of element plots for a time span of ten million years.     

The construction of a third order secular

In this part we find out the 24 equations of secular perturbation equations for the subsystem J-S-U-N. The solution of these equations by the Lagrange-Laplace procedure and the Eigen value Eigen vector is analysed. Also we refer to Hurwitz theorem to test stability.     

Secular Post-NEWTONian Effects for the Restricted Three-Body Problem in TREDER'S Dynamics with Induced Inertia

In the framework of TREDER'S dynamics, in which according to the MACH -EINSTEIN -doctrin the inertial masses are induced by the gravitational interaction with the particles of the cosmos, we calculate the secular post-NEWTON ian effects for the restricted three-body problem. The dominant secular post-NEWTON ian variation of nodes and apsides is shown to be the same as in EINSTEIN'S theory of gravitation. The formulation of the gravodynamics in the HERTZ ian configuration space allows – as in General Relativity – the explanation of the effect as an “geodetic precession” of the lines of nodes and apsides.