The close triple approach

A new method is described for representing the motion in the planar problem of three bodies when all three point masses simultaneously come close to each other. The main results are (1) that the motion during the critical phase of closest approach is intimately connected with triple parabolic escape and (2) that a sufficiently close triple approach generally leads to the escape of one body witharbitrarily high asymptotic velocity.     

Dynamical evolution of triple systems with big differences in the masses of bodies. A criterion for stability

The dynamical evolution of about 1.5 million planar hierarchical triple systems with a negative total energy and different-mass bodies is investigated by computer simulations. We considered both cases — prograde and retrograde motions of bodies. For every system, calculations were carried out either till a time when the Marchal'set al. (1984) criterion of escape of a body from a triple system was satisfied (the unstable triple systems) or during 1000 rotations of a total system (the stable triple systems). Computations were carried out on three computers-Sunstations in the Physical Research Laboratory, Ahmedabad, India during several months continuously. We changed smoothly the initial value of the coefficient of hierarchy of triples $$q = r_{3 - 12} /r_{12} $$ Wherer ₁₂ is a distance between close bodiesM1,M2 andr _3–12 is a distance between their center of masses and a distant bodyM3. We define critical (minimum) values of the coefficientq of hierarchy of stable triple systems with a relative accuracy δq=1%. Ratios of masses of bodies belong to the interval [0.13, 244.00]. A possibility of extention of these results for hierarchical subsystems with different multiplicities inside clusters is discussed.     

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.     

Strong triple interactions in the general three-body problem

Strong three-body interactions play a decisive role in the course of the dynamical evolution of triple systems having positive as well as negative total energies. These interactions may produce qualitative changes in the relative motions of the components. In triple systems with positive or zero total energies the processes of formation, disruption or exchange of the components of binaries take place as the result of close approaches of the three single bodies or as the result of the passages of single bodies past wide or hard binaries. In the triple systems with negative energies, the strong triple interactions may result in an escape from the system as well as a formation of a hard final binary. This paper summarizes the general results of the studies of the strong three-body interactions in the triple systems with positive and negative energies. These studies were conducted at the Leningrad University Observatory by computer simulations during 1968–1989.     

Close triple approaches and escape in the three-body problem

We have studied a total of 5000 close triple approaches resulting in escape, for equal-mass systems with zero initial velocities. Escape is shown to take place in the majority of the cases after a fly-by close triple approach when the escaper passes near the centre of mass along an almost straight-line orbit. A number of configurational and kinematical parameters are introduced in order to characterize the triple approach. The distributions of these parameters are investigated. A comparison with 831 examples in the vicinity of the so-called Pythagorean problem is carried out. We find that the general features of close triple approaches which result in escape are the same for both types of systems.     

Global chaoticity in the Pythagorean three-body problem

The effects of small changes in the initial conditions of the Pythagorean three-body problem are investigated by computer simulations. This problem consists of three interacting bodies with masses 3, 4 and 5 placed with zero velocities at the apices of a triangle with sides 3, 4 and 5. The final outcome of this motion is that two bodies form a binary and the third body escapes. We attempt to establish regions of the initial positions which give regular and chaotic motions. The vicinity of a small neighbourhood around the standard initial position of each body defines a regular region. Other regular regions also exist. Inside these regions the parameters of the triple systems describing the final outcome change continuously with the initial positions. Outside the regular regions the variations of the parameters are abrupt when the initial conditions change smoothly. Escape takes place after a close triple approach which is very sensitive to the initial conditions. Time-reversed solutions are employed to ensure reliable numerical results and distinguish between predictable and non-predictable motions. Close triple approaches often result in non-predictability, even when using regularization; this introduces fundamental difficulties in establishing chaotic regions.     

The rectilinear restricted problem of three bodies

In this paper we discuss some aspects of the isosceles case of the rectilinear restricted problem of three bodies, where two primaries of equal mass move on rectilinear ellipses, and the particle is confined to the symmetry axis of the system. In particular, the behaviour near a collision of the primaries and also near a collision of all three bodies is investigated. It is shown that this latter singularity is a triple collision in the sense of Siegel's theory. Furthermore, asymptotic expansions for the particle's motion during a parabolic and a hyperbolic escape are derived.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 27–September 2, 1972.     

Main features of dynamical escape from three-dimensional triple systems

The dynamical evolution of triple systems with equal and unequal-mass components and different initial velocities is studied. It is shown that, in general, the statistical results for the planar and three-dimensional triple systems do not differ significantly. Most (about 85%) of the systems disrupt; the escape of one component occurs after a triple approach of the components. In a system with unequal masses, the escaping body usually has the smallest mass. A small fraction (about 15%) of stable or long-lived systems is formed if the angular momentum is non-zero. Averages, distributions and coefficients of correlations of evolutionary characteristics are presented: the life-time, angular momentum, numbers of wide and close triple approaches of bodies, relative energy of escapers, minimum perimeter during the last triple approach resulting in escape, elements of orbits of the final binary and escaper.     

Classification of orbits in the plane isosceles three-body problem

The general plane isosceles three-body problem is considered for different ratios of the central body mass to the masses of other bodies. The central body goes through the middle of the segment connecting the other bodies along the perpendicular to this segment. The initial conditions are chosen by two parameters: the virial ratio k and the parameter     , where r˙ is the relative velocity of the 'outer' bodies, and R˙ is the velocity of the 'central' body with respect to the mass centre of the 'outer' bodies. The equations of motion are numerically integrated until one of three times: the time of escape of the central body, its time of ejection with   R >100 d   , or 1000 τ (here d is the mean size, and τ is the mean crossing time of the triple system). The regions corresponding to escapes of the central body after different numbers of triple approaches are found at the plane of parameters   k ∈(0,1)  and   μ ∈(-1,1)  . The regions of stable motions are revealed. The zones of regular and stochastic orbits are outlined. The fraction of stochastic trajectories increases with the central mass. The fraction of stable orbits is highest for equal masses of the bodies.     

Mass effects in the problem of three bodies

This is a study of the dynamical behavior of three point masses moving under their mutual gravitational attraction in a plane. The initial positions and velocities are identical for all cases studied and only the masses of the participating bodies change in the series of numerical experiments. In this way the effect of the coupling terms in the differential equations of motion are investigated. The motion in all 125 cases begins with an interplay between the three bodies, followed by temporary ejections or by an eventual escape. The total mass of the system is kept constant while the massratios change from 1 to 5. The initial velocities being zero, the total energy is negative in all cases.Approximately 74% of the cases disintegrated (i.e. two bodies formed a binary and the third body escaped) in less than 140 time units, 47% in less than 50 time units and 10% ended in escape in less than 10 time units. Considering three stars with total mass 12M , initially placed at 3, 4 and 5 parsec distances (or three galaxies with mass 2.4×10¹² M , initially placed 30, 40 and 50 kpc apart), the unit of time (approximately the crossing time) becomes 1.5×10⁷ y (3.2×10⁷ y). The average time of disintegration was found to be of the order of 10⁹ y. The average semi-major axis of the binaries left behind after disintegration was 0.7 parsec and the average value of the eccentricity was 0.76. The effect of the masses on the escapes was established and it was found that the bodynot with the smallest mass escaped in 13% of the disintegrated cases. The cases which did not disintegrate in 150 time units were analyzed in detail and the time of their eventual escape was estimated.The numerical results are tabulated regarding escape time, ejection period, total energy, escape energy, terminal velocity, semi-major axis, and eccentricity.The evolution of triple systems is followed from interplays through ejections to escapes and the orbital parameters for the separation of these classes are estimated.     

Numerical integration methods in dynamical astronomy

We review what kinds of numerical integrators are used by astronomers in the field of dynamical astronomy and to what problems they are applied. This review is based on the questionaires distributed mainly to the members of IAU Commission 7 (Celestial Mechanics). Because of the restriction to the Commission 7 members, the answers are mainly from astronomers in the solar system dynamics and problems mentioned in the answers are also related to celestial bodies in the solar system. Other than above, two questions, how to check the precision or accuracy of numerically integrated results and how to treat a close approach, are also surveyed. The problem of the suitable choice of a numerical integrator from various numerical integrators is out of the scope of this review, and it depends strongly on the dynamical nature of a particular dynamical system and the required accuracy.     

The three-body problem near triple collision

The behaviour of three gravitationally interacting particles in a plane, which approach each other almost on a central configuration, is studied. Linearization near a Lagrangean solution and matching methods lead to the following results: (i) After a close triple encounter in the planar problem of three bodies, one particle generally escapes with an arbitrarily large asymptotic velocity. (ii) Particular cases of actual triple collisions may be extended by the method of Easton.     

Dynamics of rotating triple systems

We investigate the dynamical evolution of 100 000 rotating triple systems with equal-mass components. The system rotation is specified by the parameter ω=?c²E, where c and E are the angular momentum and total energy of the triple system, respectively. We consider ω=0.1,1, 2, 4, 6 and study 20 000 triple systems with randomly specified coordinates and velocities of the bodies for each ω. We consider two methods for specifying initial conditions: with and without a hierarchical structure at the beginning of the evolution. The evolution of each system is traced until the escape of one of the bodies or until the critical time equal to 1000 mean system crossing times. For each set of initial conditions, we computed parameters of the final motions: orbital parameters for the final binary and the escaping body. We analyze variations in the statistical characteristics of the distributions of these parameters with ω. The mean disruption time of triple systems and the fraction of the systems that have not been disrupted in 1000 mean crossing times increase with ω. The final binaries become, on average, wider at larger angular momenta. The distribution of their eccentricities does not depend on ω and generally agrees with the theoretical law f(e)=2e. The velocities of the escaping bodies, on average, decrease with increasing angular momentum of the triple system. The fraction of the angles between the escaping-body velocity vector and the triple-system angular momentum close to 90° increases with ω. Escapes in the directions opposite to rotation and prograde motions dominate at small and large angular momenta, respectively. For slowly rotating systems, the angular momentum during their disruption is, on average, evenly divided between the escaping body and the final binary, whereas in rapidly rotating systems, about 80% of the angular momentum is carried away by the escaping component. We compare our numerical simulations with the statistical theory of triple-system disruption.     

Treatment of close approaches in the numerical integration of the gravitational problem ofN bodies

One of the main difficulties encountered in the numerical integration of the gravitationaln-body problem is associated with close approaches. The singularities of the differential equations of motion result in losses of accuracy and in considerable increase in computer time when any of the distances between the participating bodies decreases below a certain value. This value is larger than the distance when tidal effects become important, consequently,numerical problems are encounteredbefore the physical picture is changed. Elimination of these singularities by transformations is known as the process of regularization. This paper discusses such transformations and describes in considerable detail the numerical approaches to more accurate and faster integration. The basic ideas of smoothing and regularization are explained and applications are given.     

Dynamical evolution of triple systems

This article reviews numerical experiments on the three-body problem carried out at the Leningrad University Astronomical Observatory during the past 20 years. Systematic studies of triple systems with negative total energy have yielded the following main results. Most (95%) of the systems decay; the decay always occurs after a close triple approach of the components. In a system with unequal masses, the escaping body usually has the smallest mass. A small fraction (5%) of quasi-stable systems is formed if the angular momentum is non-zero. The qualitative evolution in three-dimensional cases is the same as for planar systems. Small changes in initial conditions sometimes lead to substantial differences in the final outcome. The decay of triple systems is a stochastic process similar to radioactive decay. The estimated mean lifetime is 100 crossing times for equal-mass components and decreases for increasing mass dispersion.A classification of the close triple approaches which lead to immediate escape is given for equal-mass systems as well as for selected sets of unequal components. Detailed studies of close triple approaches by computer simulations reveal that the early evolutions is determined by the initial ratio of the interaction forces. The review concludes by discussing applications of the results to observational problems of stellar and extragalactic systems.     

Numerical explorations of the rectilinear problem of three bodies

We present some results of a numerical exploration of the rectilinear problem of three bodies, with the two outer masses equal. The equations of motion are first given in relative coordinates and in regularized variables, removing both binary collision singularities in a single coordinate transformation. Among our most important results are seven periodic solutions and three symmetric triple collision solutions. Two of these periodic solutions have been continued into families, the outer massm ₃ being the family parameter. One of these families exists for all masses while the second family is a branch of the first at a second-kind critical orbit. This last family ends in a triple collision orbit.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.     

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.     

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.     

Trojan-type orbits in the system of two gravitational centers with variable masses and separation

Trojan type orbits in the system of two gravitational centers with variable separation are studied within the framework of the restricted problem of three bodies. The backward numerical integration of the equations of motion of the bodies starting in the triangular libration pointsL ₄ andL ₅ (reverse problem) finds the breakdown of librations as the separation decreases because of the mass gain of the smaller component and an approach of the body of negligible, mass to the latter up to the distance below its sphere of action with a relative velocity approximately equal to the escape one on this sphere. The breakdown of librations aboutL ₅ occurs under the mass gain of the smaller component considerably larger than in the case ofL ₄ and implications are made for the asymmetry of the number of librators aboutL ₄ andL ₅ in the solar system (Greeks and Trojans). Other parameters of the libration motion near 1/1 commensurability are obtained, namely, the asymmetry of the libration amplitudes about the triangular points as well as the values of periods and amplitudes within the limits of those for real Trojan asteroids. Trojans could be supposedly, formed inside the Proto-jupiter and escape during its intensive mass loss.