Search for and study of nearly periodic orbits in the plane problem of three equal-mass bodies

We analyze nearly periodic solutions in the plane problem of three equal-mass bodies by numerically simulating the dynamics of triple systems. We identify families of orbits in which all three points are on one straight line (syzygy) at the initial time. In this case, at fixed total energy of a triple system, the set of initial conditions is a bounded region in four-dimensional parameter space. We scan this region and identify sets of trajectories in which the coordinates and velocities of all bodies are close to their initial values at certain times (which are approximately multiples of the period). We classify the nearly periodic orbits by the structure of trajectory loops over one period. We have found the families of orbits generated by von Schubart’s stable periodic orbit revealed in the rectilinear three-body problem. We have also found families of hierarchical, nearly periodic trajectories with prograde and retrograde motions. In the orbits with prograde motions, the trajectory loops of two close bodies form looplike structures. The trajectories with retrograde motions are characterized by leafed structures. Orbits with central and axial symmetries are identified among the families found.     

The rectilinear three-body problem

The rectilinear equal-mass and unequal-mass three-body problems are considered. The first part of the paper is a review that covers the following items: regularization of the equations of motion, integrable cases, triple collisions and their vicinities, escapes, periodic orbits and their stability, chaos and regularity of motions. The second part contains the results of our numerical simulations in this problem. A classification of orbits in correspondence with the following evolution scenarios is suggested: ejections, escapes, conditional escapes (long ejections), periodic orbits, quasi-stable long-lived systems in the vicinity of stable periodic orbits, and triple collisions. Homothetic solutions ending by triple collisions and their dependence on initial parameters are found. We study how the ejection length changes in response to the variation of the triple approach parameters. Regions of initial conditions are outlined in which escapes occur after a definite number of triple approaches or a definite time. In the vicinity of a stable Schubart periodic orbit, we reveal a region of initial parameters that corresponds to trajectories with finite motions. The regular and chaotic structure of the manifold of orbits is mostly defined by this periodic orbit. We have studied the phase space structure via Poincaré sections. Using these sections and symbolic dynamics, we study the fine structure of the region of initial conditions, in particular the chaotic scattering region.     

A search for collision orbits in the free-fall three-body problem I. Numerical procedure

A numerical procedure is devised to find binary collision orbits in the free-fall three-body problem. Applying this procedure, families of binary collision orbits are found and a sequence of triple collision orbits are positioned. A property of sets of binary collision orbits which is convenient to search triple collision orbits is found. Important numerical results are formulated and summarized in the final section.     

Qualitative study of the planar isosceles three-body problem

We consider the particular case of the planar three body problem obtained when the masses form an isosceles triangle for all time. Various authors [1, 2, 12, 8, 9, 13, 10] have contributed in the knowledge of the triple collision and of several families of periodic orbits in this problem. We study the flow on a fixed level of negative energy. First we obtain a topological representation of the energy manifold including the triple collision and infinity as boundaries of that manifold. The existence of orbits connecting the triple collision and infinity gives some homoclinic and heteroclinic orbits. Using these orbits and the homothetic solutions of the problem we can characterize orbits which pass near triple collision and near infinity by pairs of sequences. One of the sequences describes the regions visited by the orbit, the other refers to the behaviour of the orbit between two consecutive passages by a suitable surface of section. This symbolic dynamics which has a topological character is given in an abstract form and after it is applied to the isosceles problem. We try to keep globality as far as possible. This strongly relies on the fact that the intersection of some invariant manifolds with an equatorial plane (v=0) have nice spiraling properties. This can be proved by analytical means in some local cases. Numerical simulations given in Appendix A make clear that these properties hold globally.     

A Search for Collision Orbits in the Free-Fall Three-Body Problem II

A numerical procedure to systematically find collision orbits in the planar three-body problem has been developed in the preceding paper (Tanikawa et al., 1995). Using this procedure, a search for binary and triple collision orbits has been carried out in the free-fall three-body problem. Some detailed structures of a part of the initial value space are discussed. Various interesting orbits have been found. Examples are oscillatory orbits in which ejected particles change from ejection to ejection, and orbits which are not isosceles initially but nearly isosceles after escape. Some results of isosceles problems (Simó and Martínez, 1988) are extended to non-isosceles problems. This revised version was published online in July 2006 with corrections to the Cover Date.     

Multiple Periodic Orbits in the Hill Problem with Oblate Secondary

We study the multiple periodic orbits of Hill’s problem with oblate secondary. In particular, the network of families of double and triple symmetric periodic orbits is determined numerically for an arbitrary value of the oblateness coefficient of the secondary. The stability of the families is computed and critical orbits are determined. Attention is paid to the critical orbits at which families of non-symmetric periodic orbits bifurcate from the families of symmetric periodic orbits. Six such bifurcations are found, one for double-periodic and five for triple-periodic orbits. Critical orbits at which families of sub-multiple symmetric periodic orbits bifurcate are also discussed. Finally, we present the full network of families of multiple periodic orbits (up to multiplicity 12) together with the parts of the space of initial conditions corresponding to escape and collision orbits, obtaining a global view of the orbital behavior of this model problem.     

Non-schubart periodic orbits in the rectilinear three-body problem

In the present paper, in the rectilinear three-body problem, we qualitatively follow the positions of non-Schubart periodic orbits as the mass parameter changes. This is done by constructing their characteristic curves. In order to construct characteristic curves, we assume a set of properties on the shape of areas corresponding to symbol sequences. These properties are assured by our preceding numerical calculations. The main result is that characteristic curves always start at triple collision and end at triple collision. This may give us some insight into the nature of periodic orbits in the N-body problem.     

Families of periodic collision orbits in the general three-body problem

In the general three-body problem, in a rotating frame of reference, a symmetric periodic solution with a binary collision is determined by the abscissa of one body and the energy of the system. For different values of the masses of the three bodies, the symmetric periodic collision orbits form a two-parametric family. In the case of equal masses of the two bodies and small mass of the third body, we found several symmetric periodic collision orbits similar to the corresponding orbits in the restricted three-body problem. Starting with one symmetric periodic collision orbit we obtained two families of such orbits. Also starting with one collision orbit in the Sun-Jupiter-Saturn system we obtained, for a constant value of the mass ratio of two bodies, a family of symmetric periodic collision orbits.     

Orbits of galactic globular clusters

New absolute proper motions of the two globular clusters NGC 4147 and NGC 6218 and a specific model of galactic mass distribution are used to integrate their orbits numerically. The resulting values of orbital parameters and their variation due to the uncertainty of the initial values are discussed. Furthermore, the deviations from the time average of the virial theorem are determined.     

Periodic orbits in the general three-body problem and the relationship between them

We study the regions of finite motions in the vicinity of three simple stable periodic orbits in the general problem of three equal-mass bodies with a zero angular momentum. Their distinctive feature is that one of the moving bodies periodically passes through the center of mass of the triple system. We consider the dynamical evolution of plane nonrotating triple systems for which the initial conditions are specified in such a way that one of the bodies is located at the center of mass of the triple system. The initial conditions can then be specified by three parameters: the virial coefficient k and the two angles, φ₁ and φ₂, that characterize the orientation of the velocity vectors for the bodies. We scanned the region of variation in these parameters k∈(0, 1); φ₁, φ₂∈(0, π) at steps of δk=0.01; δφ₁=δφ₂=1° and identified the regions of finite motions surrounding the periodic orbits. These regions are isolated from one another in the space of parameters (k, φ₁, φ₂). There are bridges that correspond to unstable orbits with long lifetimes between the regions. During the evolution of these metastable systems, the phase trajectory can “stick” to the vicinity of one of the periodic orbits or move from one vicinity to another. The evolution of metastable systems ends with their breakup.     

A method for calculating probability of collision between space objects

A method is developed to calculate probability of collision. Based on geometric features of space objects during the encounter, it is reasonable to separate the radial orbital motions from those in the cross section for most encounter events that occur in a near-circular orbit. Therefore, the probability of collision caused by differences in both altitude of the orbit in the radial direction and the probability of collision caused by differences in arrival time in the cross section are calculated. The net probability of collision is expressed as an explicit expression by multiplying the above two components. Numerical cases are applied to test this method by comparing the results with the general method. The results indicate that this method is valid for most encounter events that occur in near-circular orbits.     

Superosculating intermediate orbits: Fourth-and fifth-order tangencies to trajectories of perturbed motion

The theory of superosculating intermediate orbits previously suggested by the author is developed. A new class of orbits with a fourth-order tangency to the actual trajectory of a celestial body at the initial time is constructed. Orbits with a fifth-order tangency have been constructed for the first time. The motion in the constructed orbits is represented as a combination of two motions: the motion of a fictitious attracting center with a variable mass and the motion relative to this center. The first motion is generally parabolic, while the second motion is described by the equations of the Gylden—Mestschersky problem. The variation in the mass of the fictitious center obeys Mestschersky’s first and combined laws. The new orbits represent more accurately the actual motion in the initial segment of the trajectory than an osculating Keplerian orbit and other existing analogues. Encke’s generalized methods of special perturbations in which the constructed intermediate orbits are used as reference orbits are presented. Numerical simulations using the approximations of the motions of Asteroid Toutatis and Comet P/Honda—Mrkos—Pajdu?áková as examples confirm that the constructed orbits are highly efficient. Their application is particularly beneficial in investigating strongly perturbed motion.     

A survey of near-mean-motion resonances between Venus and Earth

It is known since the seminal study of Laskar (1989) that the inner planetary system is chaotic with respect to its orbits and even escapes are not impossible, although in time scales of billions of years. The aim of this investigation is to locate the orbits of Venus and Earth in phase space, respectively, to see how close their orbits are to chaotic motion which would lead to unstable orbits for the inner planets on much shorter time scales. Therefore, we did numerical experiments in different dynamical models with different initial conditions—on one hand the couple Venus–Earth was set close to different mean motion resonances (MMR), and on the other hand Venus’ orbital eccentricity (or inclination) was set to values as large as e = 0.36 (i = 40°). The couple Venus–Earth is almost exactly in the 13:8 mean motion resonance. The stronger acting 8:5 MMR inside, and the 5:3 MMR outside the 13:8 resonance are within a small shift in the Earth’s semimajor axis (only 1.5 percent). Especially Mercury is strongly affected by relatively small changes in initial eccentricity and/or inclination of Venus, and even escapes for the innermost planet are possible which may happen quite rapidly.     

Orbital Distribution Arbitrarily Close to the Homothetic Equilateral Triple Collision in the Free‐Fall Three‐Body Problem with Equal Masses

The existence of escape and nonescape orbits arbitrarily close to the homothetic equilateral triplecollision orbit is considered analytically in the threebody problem with zero initial velocities and equal masses. It is proved that escape orbits in the initial condition space are distributed around three kinds of isosceles orbits. It is also proved that nonescape orbits are distributed in between the escape orbits where different particles escape. In order to show this, it is proved that the homotheticequilateral orbit is isolated from other triplecollision orbits as far as the collision at the first triple encounter is concerned. Moreover, the escape criterion is formulated in the planarisosceles problem and translated into the words of regularizing variables. The result obtained by us explains the orbital structure numerically.     

Stability of interplay motions

A family of rectilinear periodic solutions of the three-body problem, in which the central body collides alternately with each of the two other bodies, is investigated numerically for all values of the three masses. It is found that for every mass combination there exists just one solution of this kind. The linear stability of the orbits with respect to arbitrary three-dimensional perturbations is also investigated. Domains of stability and instability are displayed in a triangular mass diagram. Their boundaries form one-parameter families of critical orbits, which are tabulated. Limiting cases where one or two masses vanish are studied in detail. The domains of stability cover nearly one half of the total area in the mass diagram: this reinforces the conclusion that real triple stars might have motions of a kind entirely different from the usual hierarchical arrangement.     

Öpik-type collision probability for high-inclination orbits

The classical Öpik theory provides an estimate of the collision probability between two bodies on bound, heliocentric or planetocentric orbits under restrictive assumptions of: (i) constant eccentricity and inclination, and (ii) uniform circulation of the longitude of node and argument of pericenter. These assumptions are violated whenever either of the orbits has a large inclination with respect to the local Laplace plane or large eccentricity, and their motion is perturbed by an exterior (tidal) gravitational field of a planet or the Sun. In this situation, known as the Lidov–Kozai regime, the eccentricity and inclination values exhibit large and correlated oscillations. At the same time, the longitude of node and the argument of pericenter may have strongly nonlinear time evolution, with the latter being even bound to a small interval of values. Here we develop a new Öpik-type collision probability theory which is valid even for highly inclined and/or eccentric orbits of the projectile. We assume that the orbit of the target is circular and in the local Laplace plane. Such a generalized setting is necessary, as an example, to correctly estimate the terrestrial impact fluxes of sporadic micrometeoroids on high-inclination orbits (notably those from the toroidal source and the associated helion and anti-helion arcs).     

Collisions of two solid bodies in Keplerian orbits; new orbits and systematic effects

It has been shown in various papers dealing with systems of colloding bodies in a Keplerian field that the dynamical evolution does not depend only on the initial orbital conditions. This is a consequence of the wide range of orbits generated by the collision process. From the study of a few pairs of orbits we examine what factors which produce that variety of orbits, and search for systematic effects. The role of the positions along the orbits, of inelasticity, of size, of mass and of relative inclination is emphasized.     

The possible motions of a satellite about an oblate planet

The regions of motions of a satellite for given values of energy and angular momentum about polar axes are shown. Special attention is paid to the circular equatorial orbits which have been shown to be Hill stable. The anomalistic and the nodal period for the motions near to the circular equatorial orbits have been found.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation in the present configuration of the two planets

The dynamics of the two Jupiter triangular libration points perturbed by Saturn is studied in this paper. Unlike some previous works that studied the same problem via the pure numerical approach, this study is done in a semianalytic way. Using a literal solution, we are able to explain the asymmetry of two orbits around the two libration points with symmetric initial conditions. The literal solution consists of many frequencies. The amplitudes of each frequency are the same for both libration points, but the initial phase angles are different. This difference causes a temporary spatial asymmetry in the motions around the two points, but this asymmetry gradually disappears when the time goes to infinity. The results show that the two Jupiter triangular libration points should have symmetric spatial stable regions in the present status of Jupiter and Saturn. As a test of the literal solution, we study the resonances that have been extensively studied in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006). The resonance structures predicted by our analytic theory agree well with those found in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006) via a numerical approach. Two kinds of chaotic orbits are discussed. They have different behaviors in the frequency map. The first kind of chaotic orbits (inner chaotic orbits) is of small to moderate amplitudes, while the second kind of chaotic orbits (outer chaotic orbits) is of relatively larger amplitudes. Using analytical theory, we qualitatively explain the transition process from the inner chaotic orbits to the outer chaotic orbits with increasing amplitudes. A critical value of the diffusion rate is given to separate them in the frequency map. In a forthcoming paper, we will study the same problem but keep the planets in migration. The time asymmetry, which is unimportant in this paper, may cause an observable difference in the two Jupiter Trojan groups during a very fast planet migration process.