A Review on Co-orbital Motion in Restricted and Planetary Three-body Problems

The 1:1 mean motion resonance may be referred to as the lowest order mean motion resonance in restricted or planetary three-body problems. The five well-known libration points of the circular restricted three-body problem are five equilibriums of the 1:1 resonance. Coorbital motion may take different shapes of trajectory. In case of small orbital eccentricities and inclinations, tadpole-shape and horseshoe-shape orbits are well-known. Other 1:1 libration modes different from the elementary ones can exist at moderate or large eccentricities and inclinations. Coorbital objects are not rare in our solar system, for example the Trojans asteroids and the coorbital satellite systems of Saturn. Recently, dozens of coorbital bodies have been identified among the near-Earth asteroids. These coorbital asteroids are believed to transit recurrently between different 1:1 libration modes mainly due to orbital precessions, planetary perturbations, and other possible effects. The Hamiltonian system and the Hill’s three-body problem are two effective approaches to study coorbital motions. To apply the perturbation theory to the Hamiltonian system, standard procedures involve the development of the disturbing function, averaging and normalization, theory of ideal resonance model, secular perturbation theory, etc. Global dynamics of coorbital motion can be revealed by the Hamiltonian approach with a suitable expansion. The Hill’s problem is particularly suitable for the studies on the relative motion of two coorbital bodies during their close encounter. The Hill’s equation derived from the circular restricted three-body problem is well known. However, the general Hill’s problem whose equation of motion takes exactly the same form applies to the non-restricted case where the mass of each body is non-negligible, namely the planetary case. The Hill’s problem can be transformed into a “canonical shape” so that the averaging principle can be applied to construct a secular perturbation theory. Besides the two analytical theories, numerical methods may be consulted, for example the approach of periodic orbit, the surface of section, and the computation of invariant manifolds carried by equilibriums or periodic orbits.     

Precession of the orbital nodes of Jupiter and Saturn triggered by the mutual perturbation: A model of two rings

The problem of the precession of the orbital planes of Jupiter and Saturn under the influence of mutual gravitational perturbations was formulated and solved using a simple dynamical model. Using the Gauss method, the planetary orbits are modeled by material circular rings, intersecting along the diameter at a small angle α. The planet masses, semimajor axes and inclination angles of orbits correspond to the rings. What is new is that each ring has an angular momentum equal to the orbital angular momentum of the planet. Contrary to popular belief, it was proved that the orbital resonance 5: 2 does not preclude the use of the ring model. Moreover, the period of averaging of the disturbing force (T ≈ 1332 yr) proves to be appreciably greater than a conventionally used period (≈900 yr). The mutual potential energy of rings and the torque of gravitational forces between the rings were calculated. We compiled and solved the system of differential equations for the spatial motion of rings. It was established that a perturbing torque causes the precession and simultaneous rotation of the orbital planes of Jupiter and Saturn. Moreover, the opposite orbit nodes on the Laplace plane coincide and perform a secular movement in retrograde direction with the same velocity of 25.6″/yr and the period T _J = T _S ≈ 50687 yr. These results are close to those obtained in the general theory (25.93″/yr), which confirms the adequacy of the developed model. It was found that the vectors of the angular velocity of orbital rings move counterclockwise over circular cones and describe circles on the celestial sphere with radii β₁ ≈ 0.8403504° (Saturn) and β₂ ≈ 0.3409296° (Jupiter) around the point which is located at an angular distance of 1.647607° from the ecliptic pole.     

Periodic orbits of the planetary type and their stability

A review is presented of periodic orbits of the planetary type in the general three-body problem and fourbody problem and the restricted circular and elliptic tnreebody problem. These correspond to planetary systems with one Sun and two or three planets (or a planet and its satellites), the motion of asteoids and also planetary systems with two Suns. The factors which affect the stability of the above configurations are studied in connection with resonance or additional perturbations. Finally, the correspondence of the periodic orbits in the restricted three-body problem with the fixed points obtained by the method of averaging or the method of surface of section is indicated.     

Intensity of gravitational wave emitted by an oscillating Keplerian binary

This paper attempts to formulate a way for calculating the intensity of gravitational wave from two point masses in Keplerian circular and elliptic orbits. The intensity is calculated with the assumption that the orbital plane of the binary undergoes small oscillation about the equilibrium x-y plane. This problem is simplification of a physically possible orbit where one of the point masses is spinning whereby the spin-orbit force drives the orbital plane to wobble in a complicated manner. It is shown that the total energy of gravitational wave emitted by the binary in this case is dominated by the parameters which take into account the oscillation of the plane. The results presented are in fact a generalization of the classic results of Landau and Lifshitz.     

An approximate integral and solution for eccentric planetary type orbits in the restricted problem of three bodies

The circular restricted problem of three bodies is investigated analytically with respect to the problem of deriving a second integral of motion besides the well known Jacobian Integral. The second integral is searched for as a correction the angular momentum integral valid in the two body case. A partial differential equation equivalent to the problem is derived and solved approximately by an asymptotic Fourier method assuming either sufficiently small values for the dimensionless mass parameter or sufficiently large distances from the barycentre. The solution of the partial equation then leads to a function of the coordinates, velocities and time being nearly constant, which means that its variation with time is about 40–300 times less than that of the pure angular momentum. By averaging over the remaining fluctuating part of the quasi-integral we are able to integrate the first order equations using a renormalization transformation. This leads to an explicit expression for the approximate solution of the circular problem which describes the motion of the third body orbiting both primaries with nonvanishing initial eccentricity (eccentric planetary type orbits). One of the main results is an explicit formula for the frequency of the perihelion motion of the third body which depends on the mass parameter, the initial distance of the third body from the barycentre and the initial eccentricity. Finally we study orbits of the P-Type, being defined as solutions of the restricted problem with circular initial conditions (vanishing initial eccentricity).     

Precession of the Orbit of a Planet around Stars Revolving along: Low-Eccentricity Orbits in a Binary System: Analytical Solution

In our previous paper (hereafter, paper I) we presented analytical results on the non-planar motion of a planet around a binary star for the cases of the circular orbits of the components of the binary. We found that the orbital plane of the planet (the plane containing the “unperturbed” elliptical orbit of the planet), in addition to precessing about the angular momentum of the binary, undergoes simultaneously the precession within the orbital plane. We demonstrated that the analytically calculated frequency of this additional precession is not the same as the frequency of the precession of the orbital plane about the angular momentum of the binary, though the frequencies of both precessions are of the same order of magnitude. In the present paper we extend the analytical results from paper I by relaxing the assumption that the binary is circular – by allowing for a relatively small eccentricity ε of the stars orbits in the binary. We obtain an additional, ε-dependent term in the effective potential for the motion of the planet. By analytical calculations we demonstrate that in the particular case of the planar geometry (where the planetary orbit is in the plane of the stars orbits), it leads to an additional contribution to the frequency of the precession of the planetary orbit. We show that this additional, ε-dependent contribution to the precession frequency of the planetary orbit can reach the same order of magnitude as the primary, ε-independent contribution to the precession frequency. Besides, we also obtain analytical results for another type of the non-planar configuration corresponding to the linear oscillatory motion of the planet along the axis of the symmetry of the circular orbits of the stars. We show that as the absolute value of the energy increases, the period of the oscillations decreases.     

A physical interpretation of the Titius-Bode rule and its connection to the closed orbits of Bertrand's theorem

We consider the geometric Titius-Bode rule for the semimajor axes of planetary orbits. We derive an equivalent rule for the midpoints of the segments between consecutive orbits along the radial direction and we interpret it physically in terms of the work done in the gravitational field of the Sun by particles whose orbits are perturbed around each planetary orbit. On such energetic grounds, it is not surprising that some exoplanets in multiple-planet extrasolar systems obey the same relation. However,it is surprising that this simple interpretation of the Titius-Bode rule also reveals new properties of the bound closed orbits predicted by Bertrand's theorem, which has been known since 1873.     

Origin and early evolution of the planetary system

It is shown by linear stability analysis that a preplanetary (presatellite) disk of dust and gas with Keplerian velocity field can become unstable due to the collective self-gravity of the disk. The radial distribution of rings, which may result from this instability, is derived. These rings later on can be the formation sites for planets around the Sun and for satellites around the planets. The derived orbits are shown to be in good agreement with that of the planets and the satellites (of Jupiter, Saturn, and Uranus). Predictions and conclusions seem to be possible for the existence of three yet unknown Uranian satellites, the origin of the early Moon and the possible radial extension of the planetary system.     

Peculiarities of Uranus’s satellite system

The absence of Uranus’s equatorial satellites in the region of approximately equal influence of its oblateness and solar perturbations is explained in terms of an improved physical model. This model is more complete than the previously studied case of an integrable averaged problem. The model improvement stems from the fact that the inclination of Uranus’s equator to the ecliptic differs by 90° and that the orbital evolution of Uranus due to secular planetary perturbations is taken into account. The lifetime of Uranus’s hypothetical satellites in orbits with semimajor axes 1.3–7 million km can be estimated by numerically integrating the evolution equations to be ～10⁴ yr. This is the time scale on which the evolution of the orbits leads to their intersection with the orbits of inner satellites.     

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.     

Extension of the Kepler problem towards minimization of energy and gravity softening

The classical Kepler Problem consists in the determination of the relative orbital motion of a secondary body (planet) with respect to the primary body (Sun), for a given time. However, any natural system tends to have minimum energy and is subjected to differential gravitational or tidal forces (called into play mainly due to the finite size and deformability of the secondary body). We formulate the Kepler Problem taking into account the finite size of the secondary body and consider an approximation which tends towards minimum energy orbits, by increasing the dimensionality of the problem. This formulation leads to a conceivable natural explanation of the fact that the planetary orbits are characterized by small eccentricities.     

Analytical method for the effects of the asteroid belt on planetary orbits

Analytic expressions are derived for the perturbation of planetary orbits due to a thick constant density asteroid belt. The derivations include extensions and adaptations of Plakhov's analytic expressions for the perturbations in five of the orbital elements for closed orbits around Saturn's rings. The equations of Plakhov are modified to include the effect of ring thickness and additional equations are derived for the perturbations in the sixth orbital element, the mean anomaly. The gravitational potential and orbital perturbations are derived for the asteroid belt with and without thickness, and for a hoop approximation to the belt. The procedures are also applicable to Saturn's rings and the newly discovered rings of Uranus.The effects of the asteroid belt thickness on the gravitational potential coefficients and the orbital motions are demonstrated. Comparisons between the Mars orbital perturbations obtained using the analytic expressions and those obtained using numerical integration are discussed. The effects of the asteroid belt on the Earth based ranging to Mars are also demonstrated.     

Resonant periodic orbits in the exoplanetary systems

The planetary dynamics of 4/3, 3/2, 5/2, 3/1 and 4/1 mean motion resonances is studied by using the model of the general three body problem in a rotating frame and by determining families of periodic orbits for each resonance. Both planar and spatial cases are examined. In the spatial problem, families of periodic orbits are obtained after analytical continuation of vertical critical orbits. The linear stability of orbits is also examined. Concerning initial conditions nearby stable periodic orbits, we obtain long-term planetary stability, while unstable orbits are associated with chaotic evolution that destabilizes the planetary system. Stable periodic orbits are of particular importance in planetary dynamics, since they can host real planetary systems. We found stable orbits up to 60^° of mutual planetary inclination, but in most families, the stability does not exceed 20^°–30^°, depending on the planetary mass ratio. Most of these orbits are very eccentric. Stable inclined circular orbits or orbits of low eccentricity were found in the 4/3 and 5/2 resonance, respectively.     

On the habitability of the OGLE-2006-BLG-109L planetary system

We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analogue of the Solar system. Our work is inspired by the work of Malhotra & Minton. Using the linear Laplace–Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the subsystem of Jupiters may be found in a dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semimajor axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.     

Interplanetary dust particles and solar radiation

The problem of the action of the solar radiation on the motion of interplanetary dust particle is discussed. Differences between the action of electromagnetic solar radiation and that of the solar wind are explained not only from the point of view of the physical nature of these phenomena but also from the point of view of dust particle's orbital evolution. As for the electromagnetic solar radiation, general equation of motion for the particle is written and the most important consequences are: (i) the process of inspiralling toward the Sun is not the only possible motion - even spiralling from the Sun is also possible, and, (ii) the orbital plane of the particle (its inclination) may change in time. As for the solar wind, the effect corresponding to the fact that solar wind particles spread out from the Sun in nonradial direction causes that the process of inspiralling toward the Sun is in more than 50% less effective than for radial spread out; in the region of the asteroid belt (long period orbits) the process of inspiralling is changed into offspiralling. Also shift in the perihelion of dust particle's orbit exists.     

Slow modes in stellar systems with nearly harmonic potentials: II. Gravitational loss-cone instability

Using a consistent perturbation theory for collisionless disk-like and spherical star clusters, we construct a theory of slow modes for systems having an extended central region with a nearly harmonic potential due to the presence of a fairly homogeneous (on the scales of the stellar system) heavy, dynamically passive halo. In such systems, the stellar orbits are slowly precessing, centrally symmetric ellipses (2: 1 orbits). We consider star clusters with monoenergetic distribution functions that monotonically increase with angular momentum in the entire range of angular momenta (from purely radial orbits to circular ones) or have a growing region only at low angular momenta. In these cases, there are orbits with a retrograde precession, i.e., in a direction opposite to the orbital rotation of the star. The presence of a gravitational loss-cone instability, which is also observed in systems of 1: 1 orbits in near-Keplerian potentials, is associated with such orbits. In contrast to 1: 1 systems, the loss-cone instability takes place even for distribution functions monotonically increasing with angular momentum, including those for systems with circular orbits. The regions of phase space with retrograde orbits do not disappear when the distribution function is smeared in energy. We investigate the influence of a weak inhomogeneity of a heavy halo with a density that decreases with distance from the center.     

Three-dimensional time optimal double angular momentum reversal trajectory using solar sails

A new concept of three dimensional non-Keplerian trajectories with double angular momentum reversal is investigated with high performance solar sails. The main discussion of this paper is about such 3D solar inverse orbits with inner constraints. The problem is addressed in a time optimal control framework solved by an indirect method. Two typical solar inverse orbits have been achieved and presented in a 3D non-dimensional dynamic model in the Heliocentric Inertial Frame. Starting from the Earth orbit ecliptic plane, a sailcraft in the inverse orbit exhibits a butterfly shape trajectory. As such, the new orbits are symmetrical with respect to a plane which contains the Sun-perihelion line. The relation of the sail attitude angles between the two symmetrical parts of the orbits are used to reduce the simulation effort. The quasi-heliostationary property at its aphelia is demonstrated with variation of the orbital radius. Evolutions of the orbital velocity and optimal sail orientations are also outlined and discussed to benefit future design work. As is suited for space observation guaranteed by its butterfly shape, the inverse orbits are thoroughly studied in terms of the concerned parameters. The discussion of the parametric influence is ranked in order as perihelion distance r _E , required maximum position z _max, perihelion position z _f and the sail lightness number β. Suitable ranges of each parameter are adopted to illustrate the orbital variation trend. Through numerical simulations the features of such inverse orbits are further emphasized to provide an initial reference for future researchers.     

Retrograde precession of stellar orbits and gravitational loss-cone instability in galactic centers

We consider disk and spherical subsystems of stars with nearly radial orbits under conditions when the well-known radial orbit instability is not possible. This requires that the precession of stellar orbits be retrograde, i.e., in the direction opposite to the orbital rotation of stars. We show that an instability that is an analogue of the loss-cone instability known in plasma physics can then develop in the presence of a “loss cone” in the angular momentum distribution of stars, which ensures a deficit or even absence of stars with low angular momenta. Examples of systems with a loss cone are the centers of galaxies or star clusters with massive black holes. The instability can produce a flux of stars onto the galactic center, i.e., it can serve as a mechanism of fueling the nuclear activity of galaxies. Mathematically, the problem is reduced to analyzing simple characteristic equations that describe small perturbations in a disk and a sphere of radially highly elongated stellar orbits. In turn, these characteristics equations are derived through a number of successive simplifications of the general linearized Vlasov equations (i.e., the system that includes the collisionless Boltzmann kinetic equation and the Poisson equation) in action—angle variables.     

On the 2/1 resonant planetary dynamics – periodic orbits and dynamical stability

The 2/1 resonant dynamics of a two-planet planar system is studied within the framework of the three-body problem by computing families of periodic orbits and their linear stability. The continuation of resonant periodic orbits from the restricted to the general problem is studied in a systematic way. Starting from the Keplerian unperturbed system, we obtain the resonant families of the circular restricted problem. Then, we find all the families of the resonant elliptic restricted three-body problem, which bifurcate from the circular model. All these families are continued to the general three-body problem, and in this way we can obtain a global picture of all the families of periodic orbits of a two-planet resonant system. The parametric continuation, within the framework of the general problem, takes place by varying the planetary mass ratio ρ. We obtain bifurcations which are caused either due to collisions of the families in the space of initial conditions or due to the vanishing of bifurcation points. Our study refers to the whole range of planetary mass ratio values  [ρ∈ (0, ∞)]  and, therefore we include the passage from external to internal resonances. Thus, we can obtain all possible stable configurations in a systematic way. As an application, we consider the dynamics of four known planetary systems at the 2/1 resonance and we examine if they are associated with a stable periodic orbit.     

超短周期系外行星研究进展

超短周期(ultra-short-period,USP)行星是指轨道周期小于1 d的系外行星,是近年来系外行星研究领域中一个新的前沿目标。USP行星的搜寻与确认需要借助傅里叶变换(Fourier transform,FT)和盒最小二乘法(the box least,BLS)等光变曲线分析算法,以筛选和确认精准的周期信号。利用统计方法可得到目前USP行星的轨道周期、行星半径、宿主恒星类型等分布特征。大部分USP行星半径小于2R⊕,受行星质量限制,大多数USP行星无法通过视向速度信号测得精确的行星质量。根据已有的观测结果可算出,部分USP行星的质量小于10M⊕,由此推测这些USP的组成更接近金属与岩石混合的类地行星。由于密近轨道可能发生光致蒸发等物质损失过程,USP行星大气的存在情况尚不明确。目前,USP行星被认为起源于热木星(hot-Jupiters)或亚海王星(sub-Neptunes),但USP行星与热木星的主星金属丰度的分布存在较大差异,亚海王星的光致蒸发起源理论可能性更高。USP行星轨道演化机制包括低偏心率轨道迁移和潮汐耗散的原位起源模型等。