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.     

Binary and triple collisions causing instability in the free-fall three-body problem

Dominant factors for escape after the first triple-encounter are searched for in the three-body problem with zero initial velocities and equal masses. By a global numerical survey on the whole initial-value space, it is found that not only a triple-collision orbit but also a particular family of binary-collision orbits exist in the set of escape orbits. This observation is justified from various viewpoints. Binary-collision orbits experiencing close triple-encounter turn out to be close to isosceles orbits after the encounter and hence lead to escape. Except for a few cases, binary-collision orbits of near-isosceles slingshot also escape. This revised version was published online in July 2006 with corrections to the Cover Date.     

Keplerian periodic orbits in the isosceles problem

The planar isosceles three-body problem where the two symmetric bodies have small masses is considered as a perturbation of the Kepler problem. We prove that the circular orbits can be continued to saddle orbits of the Isosceles problem. This continuation is not possible in the elliptic case. Their perturbed orbits tend to a continued circular one or approach a triple collision. The basic tool used is the study of the Poincaré maps associated with the periodic solutions. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

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.     

Reduction and relative equilibria for the two-body problem on spaces of constant curvature

We consider the two-body problem on surfaces of constant nonzero curvature and classify the relative equilibria and their stability. On the hyperbolic plane, for each \(q>0\) we show there are two relative equilibria where the masses are separated by a distance q. One of these is geometrically of elliptic type and the other of hyperbolic type. The hyperbolic ones are always unstable, while the elliptic ones are stable when sufficiently close, but unstable when far apart. On the sphere of positive curvature, if the masses are different, there is a unique relative equilibrium (RE) for every angular separation except \(\pi /2\). When the angle is acute, the RE is elliptic, and when it is obtuse the RE can be either elliptic or linearly unstable. We show using a KAM argument that the acute ones are almost always nonlinearly stable. If the masses are equal, there are two families of relative equilibria: one where the masses are at equal angles with the axis of rotation (‘isosceles RE’) and the other when the two masses subtend a right angle at the centre of the sphere. The isosceles RE are elliptic if the angle subtended by the particles is acute and is unstable if it is obtuse. At \(\pi /2\), the two families meet and a pitchfork bifurcation takes place. Right-angled RE are elliptic away from the bifurcation point. In each of the two geometric settings, we use a global reduction to eliminate the group of symmetries and analyse the resulting reduced equations which live on a five-dimensional phase space and possess one Casimir function.     

Analysis of some degenerate quadruple collisions

We consider the trapezoidal problem of four bodies. This is a special problem where only three degrees of freedom are involved. The blow up method of McGehee can be used to deal with the quadruple collision. Two degenerate cases are studied in this paper: the rectangular and the collinear problems. They have only two degrees of freedom and the analysis of total collapse can be done in a way similar to the one used for the collinear and isosceles problems of three bodies. We fully analyze the flow on the total collision manifold, reducing the problem of finding heteroclinic connections to the study of a single ordinary differential equation. For the collinear case, from which arises a one parameter family of equations, the analysis for extreme values of the parameter is done and numerical computations fill up the gap for the intermediate values. Dynamical consequences for possible motions near total collision as well as for regularization are obtained.Paper presented at the 1981 Oberwolfach Conference on Mathematical Methods in Celestial Mechanics.Dedicated to Prof. Szebehely on the occasion of his sixtieth birthday.     

Short Time Lyapunov Indicators in the Restricted Three-Body Problem

We applied the method of the short time Lyapunov indicators to the planar circular and to the planar elliptic restricted three-body problem in order to study the structure of the phase space in some selected regions. In the circular case we computed the short-time averages of the stretching numbers to distinguish between regular and chaotic domains of the phase space. The results obtained in this way are in good agreement with the corresponding Poincaré's surface of sections. Using the short time Lyapunov indicators we determined the detailed structure of the phase space in the semi-major axis-eccentricity plane of the test particle both in the circular and in the elliptic restricted problem (in the latter case for some values of the eccentricities of the primaries) and we studied the main features of the phase space.This revised version was published online in October 2005 with corrections to the Cover Date.     

Relative Equilibrium in the Three-Body Problem with a Rigid Body

In the three-body problem, where two bodies are punctual and the third is rigid, we prove the existence of some relative equilibrium configurations where the rigid body is either an homogeneous ball, an oblate or an elongated ball. In particular, we found conditions of relative equilibrium of Euler and Lagrange type and several families of relative equilibrium configurations, where the triangle of the two punctual bodies and the mass center of the rigid body is isosceles or having unequal sides. This revised version was published online in July 2006 with corrections to the Cover Date.     

Motion near total collapse in the planar isosceles three-body problem

We describe in detail the qualitative behavior of solutions of the planar isosceles problem which come close to triple collision. Using a method of McGehee, one may describe a neighborhood of triple collision for any mass ratios. For sufficiently small mass ratios, we exhibit infinitely many distinct periodic and collision/ejection solutions of this problem.Partially supported by N.S.F. Grant MCS 81-01855.Paper presented at the 1981 Oberwolfach Conference on Mathematical Methods in Celestial Mechanics.     

Triple approaches in the plane isosceles equal-mass three-body problem

We analyze flyby-type triple approaches in the plane isosceles equal-mass three-body problem and in its vicinity. At the initial time, the central body lies on a straight line between the other two bodies. Triple approaches are described by two parameters: virial coefficient k and parameter $\mu = \dot r/\sqrt {\dot r^2 + \dot R^2 }$ , where $\dot r$ is the relative velocity of the extreme bodies and $\dot R$ is the velocity of the central body relative to the center of mass of the extreme bodies. The evolution of the triple system is traceable until the first turn or escape of the central body. The ejection length increases with closeness of the triple approach (parameter k). The longest ejections and escapes occur when the extreme bodies move apart with a low velocity at the time of triple approach. We determined the domain of escapes; it corresponds to close triple approaches (k>0.8) and to μ in the range ?0.2<μ<0.7. For small deviations from the isosceles problem, the evolution does not differ qualitatively from the isosceles case. The domain of escapes decreases with increasing deviations. In general, the ejection length increases for wide approaches and decreases for close approaches.     

High-order expansion of the solution of preliminary orbit determination problem

A method for high-order treatment of uncertainties in preliminary orbit determination is presented. The observations consist in three couples of topocentric right ascensions and declinations at three observation epochs. The goal of preliminary orbit determination is to compute a trajectory that fits with the observations in two-body dynamics. The uncertainties of the observations are usually mapped to the phase space only when additional observations are available and a least squares fitting problem is set up. A method based on Taylor differential algebra for the analytical treatment of observation uncertainties is implemented. Taylor differential algebra allows for the efficient computation of the arbitrary order Taylor expansion of a sufficiently continuous multivariate function. This enables the mapping of the uncertainties from the observation space to the phase space as high-order multivariate Taylor polynomials. These maps can then be propagated forward in time to predict the observable set at successive epochs. This method can be suitably used to recover newly discovered objects when a scarce number of measurements is available. Simulated topocentric observations of asteroids on realistic orbits are used to assess the performances of the method.     

Chaotic Asteroidal Dynamics and Maximum Lyapunov Exponents

This paper describes the results of studies of dynamical chaos in the problem of the orbital dynamics of asteroids near the 3 : 1 mean-motion resonance with Jupiter. Maximum Lyapunov characteristic exponents (MLCEs) are used as an indicator and a measure of the chaoticity of motion. MLCE values are determined for trajectories calculated by the numerical integration of equations of motion in the planar elliptical restricted three-body problem. The dependence of the MLCE on the problem parameters and on the initial data is analyzed. The inference is made that the domain of chaos in the phase space of the problem considered consists of two components of different nature. The values of the MLCEs observed for one of the components (namely, for the component corresponding to low-eccentricity asteroidal orbits) are compared to the theoretical estimates obtained within the framework of model of the resonance as a perturbed nonlinear pendulum.     

Families of symmetric relative periodic orbits originating from the circular Euler solution in the isosceles three-body problem

We study symmetric relative periodic orbits in the isosceles three-body problem using theoretical and numerical approaches. We first prove that another family of symmetric relative periodic orbits is born from the circular Euler solution besides the elliptic Euler solutions. Previous studies also showed that there exist infinitely many families of symmetric relative periodic orbits which are born from heteroclinic connections between triple collisions as well as planar periodic orbits with binary collisions. We carry out numerical continuation analyses of symmetric relative periodic orbits, and observe abundant families of symmetric relative periodic orbits bifurcating from the two families born from the circular Euler solution. As the angular momentum tends to zero, many of the numerically observed families converge to heteroclinic connections between triple collisions or planar periodic orbits with binary collisions described in the previous results, while some of them converge to “previously unknown” periodic orbits in the planar problem.     

Low-energy Lunar Trajectories with Lunar Flybystwo

The low-energy lunar trajectories with lunar flybys are investigated based on the Sun-Earth-Moon bicircular problem (BCP). The characteristics of the distribution of trajectories in the phase space are summarized. Using the invariant manifolds in the BCP system, the low-energy lunar trajectories with lunar flybys are sought. Then, take time as an augmented dimension in the phase space of a nonautonomous system, we present the state space map and reveal the distribution of these lunar trajectories in the phase space. Consequently, we find that the low-energy lunar trajectories exist as families, and that the every moment in the Sun-Earth-Moon synodic period can be the departure date. Finally, we analyse the velocity increment, transfer duration, and system energy for the different trajectory families, and obtain the velocity-impulse optimal family and the transfer-duration optimal family, respectively.     

Some results on the global dynamics of the regularized restricted three-body problem with dissipation

We perform an analysis of the dynamics of the circular, restricted, planar three-body problem under the effect of different kinds of dissipation (linear, Stokes and Poynting–Robertson drags). Since the problem is singular, we implement a regularization technique in the style of Levi–Civita. The effect of the dissipation is often to decrease the semi-major axis; as a consequence the minor body collides with one of the primaries. In general, it is quite difficult to find non-collision orbits using random initial conditions. However, by means of the computation of the Fast Lyapunov Indicators (FLI), we obtain a global view of the dynamics. Precisely, we detect the regions of the phase space potentially belonging to basins of attraction. This investigation provides information on the different regions of the phase space, showing both collision and non-collision trajectories. Moreover, we find periodic orbit attractors for the case of linear and Stokes drags, while in the case of the Poynting–Robertson effect no other attractors are found beside the primaries, unless a fourth body is added to counterbalance the dissipative effect.     

Tame and chaotic behavior in the planar isosceles 3-body problem

We show that every planar isosceles solution of the three-body problem encounters a collision of the symmetric particles, either forwards or backwards in time. Regularizing analytically this collision, the solution has at least a syzygy configuration and/or leads to a total collapse. Some further simple results support the intuitive image on the tame local behavior of the motion as long as it does not lead to a triple collision. As a main result we prove that total collapse singularities, can be regularized in aC ¹-fashion with respect to time, for all values of the masses. Using symbolic dynamics, the chaotic character of theC ¹-regularized solutions is pointed out.     

Dynamics of Two Planets in the 2/1 Mean-Motion Resonance

The dynamics of two planets near a first-order mean-motion resonance is modeled in the domain of the general three-body planar problem. The system studied is the pair Uranus-Neptune (2/1 resonance). The phase space of the resonance and near-resonance regions is studied by means of surfaces of section and spectral analysis techniques. After a thorough investigation of the topology of the phase space, we find that several regimes of motion are possible for the Uranus-Neptune system, and the regions of transition between the regimes of motion are the seats of chaotic motion. This revised version was published online in July 2006 with corrections to the Cover Date.