Secular and tidal evolution of circumbinary systems

We investigate the secular dynamics of three-body circumbinary systems under the effect of tides. We use the octupolar non-restricted approximation for the orbital interactions, general relativity corrections, the quadrupolar approximation for the spins, and the viscous linear model for tides. We derive the averaged equations of motion in a simplified vectorial formalism, which is suitable to model the long-term evolution of a wide variety of circumbinary systems in very eccentric and inclined orbits. In particular, this vectorial approach can be used to derive constraints for tidal migration, capture in Cassini states, and stellar spin–orbit misalignment. We show that circumbinary planets with initial arbitrary orbital inclination can become coplanar through a secular resonance between the precession of the orbit and the precession of the spin of one of the stars. We also show that circumbinary systems for which the pericenter of the inner orbit is initially in libration present chaotic motion for the spins and for the eccentricity of the outer orbit. Because our model is valid for the non-restricted problem, it can also be applied to any three-body hierarchical system such as star–planet–satellite systems and triple stellar systems.     

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.     

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.     

On stationary solutions of some canonical systems of differential equations

In our article (Zhuravlev, 1979) a formal method of constructing conditionally periodic solutions of canonical systems of differential equations with a quick-rotating phase in the case of sharp commensurability was presented. The existence of stationary (or periodic) solutions of an averaged system of differential equations corresponding to the initial system of differential equations is necessary for an effective application of the method for different problems.Evidently, the stationary solutions do not always exist but in numerous papers on stationary solutions (oscillations or motions), the conditions of existence of such solutions are very often not considered at all. Usually a simple assumption is used that the stationary solutions do exist.Otherwise it is well known that Poincaré's theory of periodic solutions (Poincaré, 1892) let one set up conditions of existence of periodic solutions in different systems of differential equations. Particularly, in papers,Mah (1949, 1956), see alsoexmah (1971), the necessary and sufficient conditions of the existence of periodic solutions of (non-canonical) systems of differential equations which are close to arbitrary non-linear systems are given. For canonical autonomous systems of differential equations the conditions of existence of periodic solutions and a method of calculation are presented in the paperMepmah (1952).In our paper another approach is given and the conditions of existence of stationary solutions of canonical systems of differential equations with a quick-rotating phase are proved. For this purpose Delaunay-Zeipel's transformation and Poincaré's small parameter method are used.     

The positions of secular resonance surfaces

The surfaces for the three strongest secular resonances have been located as a function of proper semimajor axis, eccentricity, and inclination for semimajor axes between 1.25 and 3.5 AU. The results are presented graphically. The ν₅ resonance only occurs at high inclinations (?23°). The ν₆ resonance passes through both the main belt and Mars-crossing space. The ν₁₆ resonance starts near the inner edge of the belt and, at low inclinations at least, folds around a portion of the Mars-crossing space until it runs nearly parallel with the Earth-crossing boundary.     

Secular dynamics of S-type planetary orbits in binary star systems: applicability domains of first- and second-order theories

We analyse the secular dynamics of planets on S-type coplanar orbits in tight binary systems, based on first- and second-order analytical models, and compare their predictions with full N-body simulations. The perturbation parameter adopted for the development of these models depends on the masses of the stars and on the semimajor axis ratio between the planet and the binary. We show that each model has both advantages and limitations. While the first-order analytical model is algebraically simple and easy to implement, it is only applicable in regions of the parameter space where the perturbations are sufficiently small. The second-order model, although more complex, has a larger range of validity and must be taken into account for dynamical studies of some real exoplanetary systems such as \(\gamma \) Cephei and HD 41004A. However, in some extreme cases, neither of these analytical models yields quantitatively correct results, requiring either higher-order theories or direct numerical simulations. Finally, we determine the limits of applicability of each analytical model in the parameter space of the system, giving an important visual aid to decode which secular theory should be adopted for any given planetary system in a close binary.     

Conditionally-periodic solutions of the Hamiltonian systems with multiple resonance

A method is developed to construct conditionally periodic solutions of Hamiltonian systems in the case of deep resonance and in the presence of more than one resonant vector (multiple resonance). The method depends on modification of a standard perturbation method. It is applicable to systems with more than one degree of freedom.     

Secular spin dynamics of inner main-belt asteroids

Understanding the evolution of asteroid spin states is challenging work, in part because asteroids have a variety of orbits, shapes, spin states, and collisional histories but also because they are strongly influenced by gravitational and non-gravitational (YORP) torques. Using efficient numerical models designed to investigate asteroid orbit and spin dynamics, we study here how several individual asteroids have had their spin states modified over time in response to these torques (i.e., 951 Gaspra, 60 Echo, 32 Pomona, 230 Athamantis, 105 Artemis). These test cases which sample semimajor axis and inclination space in the inner main belt, were chosen as probes into the large parameter space described above. The ultimate goal is to use these data to statistically characterize how all asteroids in the main belt population have reached their present-day spin states. We found that the spin dynamics of prograde-rotating asteroids in the inner main belt is generally less regular than that of the retrograde-rotating ones because of numerous overlapping secular spin-orbit resonances. These resonances strongly affect the spin histories of all bodies, while those of small asteroids (?40 km) are additionally influenced by YORP torques. In most cases, gravitational and non-gravitational torques cause asteroid spin axis orientations to vary widely over short (?1 My) timescales. Our results show that (951) Gaspra has a highly chaotic rotation state induced by an overlap of the s and s₆ spin-orbit resonances. This hinders our ability to investigate its past evolution and infer whether thermal torques have acted on Gaspra's spin axis since its origin.     

Geostationary secular dynamics revisited: application to high area-to-mass ratio objects

The long-term dynamics of the geostationary Earth orbits (GEO) is revisited through the application of canonical perturbation theory. We consider a Hamiltonian model accounting for all major perturbations: geopotential at order and degree two, lunisolar perturbations with a realistic model for the Sun and Moon orbits, and solar radiation pressure. The long-term dynamics of the GEO region has been studied both numerically and analytically, in view of the relevance of such studies to the issue of space debris or to the disposal of GEO satellites. Past studies focused on the orbital evolution of objects around a nominal solution, hereafter called the forced equilibrium solution, which shows a particularly strong dependence on the area-to-mass ratio. Here, we (i) give theoretical estimates for the long-term behavior of such orbits, and (ii) we examine the nature of the forced equilibrium itself. In the lowest approximation, the forced equilibrium implies motion with a constant non-zero average ‘forced eccentricity’, as well as a constant non-zero average inclination, otherwise known in satellite dynamics as the inclination of the invariant ‘Laplace plane’. Using a higher order normal form, we demonstrate that this equilibrium actually represents not a point in phase space, but a trajectory taking place on a lower-dimensional torus. We give analytical expressions for this special trajectory, and we compare our results to those found by numerical orbit propagation. We finally discuss the use of proper elements, i.e., approximate integrals of motion for the GEO orbits.     

Study and application of the resonant secular dynamics beyond Neptune

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (Celest Mech Dyn Astron, doi: 10.1007/s10569-016-9700-5, 2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or \(2012\text {VP}_{113}\), which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at \(a\in [100;300]\) AU, \(q\in [50;70]\) AU and \(I\in [30;50]^{\circ }\). In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.     

Secular frequencies of 3-D exoplanetary systems

Using a 12th order expansion of the perturbative potential in powers of the eccentricities and the inclinations, we study the secular effects of two non-coplanar planets which are not in mean–motion resonance. By means of Lie transformations (which introduce an action–angle formulation of the Hamiltonian), we find the four fundamental frequencies of the 3-D secular three-body problem and compute the long-term time evolutions of the Keplerian elements. To find the relations between these elements, the main combinations of the fundamental frequencies common to these evolutions are identified by frequency analysis. This study is performed for two different reference frames: a general one and the Laplace plane. We underline the known limitations of the linear Laplace–Lagrange theory and point out the great sensitivity of the 3-D secular three-body problem to its initial values. This analytical approach is applied to the exoplanetary system Andromedae in order to search whether the eccentricities evolutions and the apsidal configuration (libration of ) observed in the coplanar case are maintained for increasing initial values of the mutual inclination of the two orbital planes. Anne-Sophie Libert is FNRS Research Fellow.     

Secular apsidal configuration of non-resonant exoplanetary systems

Using a high-order (order 12) expansion of the perturbative potential in powers of eccentricities [Libert, A.-S., Henrard, J., 2005. Celest. Mech. Dynam. Astron. 93, 187-200], we study the secular effects of two coplanar planets which are not in mean motion resonances. The main results concern eccentricity variations, oscillation amplitude of the angular difference of the apsidal lines (Δ?) and frequency of such an oscillation. We show that this analytical approach describes correctly the behaviour of most of the exosystems and underlines the known limitations of the linear Laplace-Lagrange theory. Apsidal configuration of υ Andromedae, HD 168443, HD 169830, HD 38529, HD 74156 and HD 12661 are examined. We also point out the great sensitivity of the υ Andromedae system to the initial values (e₁(0), e₂(0) or Δ?(0)).     

Onset of secular chaos in planetary systems: period doubling and strange attractors

As a result of resonance overlap, planetary systems can exhibit chaotic motion. Planetary chaos has been studied extensively in the Hamiltonian framework, however, the presence of chaotic motion in systems where dissipative effects are important, has not been thoroughly investigated. Here, we study the onset of stochastic motion in presence of dissipation, in the context of classical perturbation theory, and show that planetary systems approach chaos via a period-doubling route as dissipation is gradually reduced. Furthermore, we demonstrate that chaotic strange attractors can exist in mildly damped systems. The results presented here are of interest for understanding the early dynamical evolution of chaotic planetary systems, as they may have transitioned to chaos from a quasi-periodic state, dominated by dissipative interactions with the birth nebula.     

Construction of conditionally-periodic solutions of canonical systems of differential equations with multiple resonance

A procedure of building of conditionally-periodic solutions of canonical systems of differential equations in a case of multiple resonances is given. A modification of Delaunay-Zeipel's method is used as a base of the procedure.     

Angular momentum exchange during secular migration of two-planet systems

We investigate the secular dynamics of two-planet coplanar systems evolving under mutual gravitational interactions and dissipative forces. We consider two mechanisms responsible for the planetary migration: star-planet (or planet-satellite) tidal interactions and interactions of a planet with a gaseous disc. We show that each migration mechanism is characterized by a specific law of orbital angular momentum exchange. Calculating stationary solutions of the conservative secular problem and taking into account the orbital angular momentum leakage, we trace the evolutionary routes followed by the planet pairs during the migration process. This procedure allows us to recover the dynamical history of two-planet systems and constrain parameters of the involved physical processes.     

Secular resonance,solar spin down,and the orbit of Mercury

A mechanism capable of accounting for the large mean eccentricity (0.175) and inclination (7°.2) of Mercury is discussed. Provided the gravitational field of the rapidly rotating primordial Sun had a sufficiently large second degree harmonic (i.e., J₂ ? order 10^?3), subsequent solar spin down would drive the orbit of Mercury through two secular resonances with Venus, one involving the precession of the line of apsides, the other one involving the regression of the nodal line. Resonance passage generates contributions to the eccentricity and inclination that are proportional to the square root of the characteristic solar spin down time. We find that an initial solar rotation l period of $\text{P ? 5}\text{1}\text{2}\text{hr}$ guarantees passage through resonance and that a spin down time of $\text{τ = Ω|dΩ/dt|}{}^{\mathrm{?1}}$ of order 10⁶ years could have produced the observed eccentricity and inclination. Such a primordial rotation rate is comparable to the measured rotations of very young stars and the spin down time appears consistent with the time scale derived for magnetic braking of the Sun's rotation by an intense solar wind during a T-Tauri stage of solar evolution.     

Analytical development of the lunisolar disturbing function and the critical inclination secular resonance

We provide a detailed derivation of the analytical expansion of the lunar and solar disturbing functions. Although there exist several papers on this topic, many derivations contain mistakes in the final expansion or rather (just) in the proof, thereby necessitating a recasting and correction of the original derivation. In this work, we provide a self-consistent and definite form of the lunisolar expansion. We start with Kaula’s expansion of the disturbing function in terms of the equatorial elements of both the perturbed and perturbing bodies. Then we give a detailed proof of Lane’s expansion, in which the elements of the Moon are referred to the ecliptic plane. Using this approach the inclination of the Moon becomes nearly constant, while the argument of perihelion, the longitude of the ascending node, and the mean anomaly vary linearly with time. We make a comparison between the different expansions and we profit from such discussion to point out some mistakes in the existing literature, which might compromise the correctness of the results. As an application, we analyze the long-term motion of the highly elliptical and critically-inclined Molniya orbits subject to quadrupolar gravitational interactions. The analytical expansions presented herein are very powerful with respect to dynamical studies based on Cartesian equations, because they quickly allow for a more holistic and intuitively understandable picture of the dynamics.