Periodic orbits and escapes in dynamical systems

We study the periodic orbits and the escapes in two different dynamical systems, namely (1) a classical system of two coupled oscillators, and (2) the Manko-Novikov metric which is a perturbation of the Kerr metric (a general relativistic system). We find their simple periodic orbits, their characteristics and their stability. Then we find their ordered and chaotic domains. As the energy goes beyond the escape energy, most chaotic orbits escape. In the first case we consider escapes to infinity, while in the second case we emphasize escapes to the central ??bumpy?? black hole. When the energy reaches its escape value, a particular family of periodic orbits reaches an infinite period and then the family disappears (the orbit escapes). As this family approaches termination it undergoes an infinity of equal period and double period bifurcations at transitions from stability to instability and vice versa. The bifurcating families continue to exist beyond the escape energy. We study the forms of the phase space for various energies, and the statistics of the chaotic and escaping orbits. The proportion of these orbits increases abruptly as the energy goes beyond the escape energy.     

Bifurcations of triple-periodic orbits

We consider families of periodic orbits in potentials symmetric with respect to thex-axis. The characteristics of triple-periodic orbits (i.e. orbits intersecting thex-axis three times) that bifurcate from the central characteristic do not have their maximum or minimum energy (or perturbation) at the point of intersection. We explain theoretically that this happens only for triple-periodic orbits and not for any other type of resonant periodic orbits and verify this fact by numerical calculations.     

Evolution of the velocity field in model potentials

We consider the evolution of the structure of the velocity field formed by the loops of the trajectory in a generalized Hénon-Heiles model potential. Box-shaped orbits alternating with periodic trajectories are shown to dominate at low values of the energy integral. Signatures of shell-and tube-shaped structures appear as the energy increases. Tube-shaped orbits are associated with stable periodic trajectories with small resonance ratios. Zones of stochastic orbits appear at values of the energy integral close to its critical value, which corresponds to the opening of the zero-velocity contour.     

Boundary-value problem formulations for computing invariant manifolds and connecting orbits in the circular restricted three body problem

We demonstrate the remarkable effectiveness of boundary value formulations coupled to numerical continuation for the computation of stable and unstable manifolds in systems of ordinary differential equations. Specifically, we consider the circular restricted three-body problem (CR3BP), which models the motion of a satellite in an Earth–Moon-like system. The CR3BP has many well-known families of periodic orbits, such as the planar Lyapunov orbits and the non-planar vertical and halo orbits. We compute the unstable manifolds of selected vertical and halo orbits, which in several cases leads to the detection of heteroclinic connections from such a periodic orbit to invariant tori. Subsequent continuation of these connecting orbits with a suitable end point condition and allowing the energy level to vary leads to the further detection of apparent homoclinic connections from the base periodic orbit to itself, or the detection of heteroclinic connections from the base periodic orbit to other periodic orbits. Some of these connecting orbits are of potential interest in space mission design.     

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.     

Energy and stability in the Full Two Body Problem

The conditions for relative equilibria and their stability in the Full Two Body Problem are derived for an ellipsoid–sphere system. Under constant angular momentum it is found that at most two solutions exist for the long-axis solutions with the closer solution being unstable while the other one is stable. As the non-equilibrium problem is more common in nature, we look at periodic orbits in the F2BP close to the relative equilibrium conditions. Families of periodic orbits can be computed where the minimum energy state of one family is the relative equilibrium state. We give results on the relative equilibria, periodic orbits and dynamics that may allow transition from the unstable configuration to a stable one via energy dissipation.      

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.     

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.     

Numerical periodic orbits of charged grains around magnetic planets

We apply a numerical searching method to investigate three-dimensional periodic orbits of charged dust particles in planetary magnetospheres. A classic generalized Stormer model of magnetic planets along with the parameters of Saturn is employed. More periodic orbits are found, besides the already known circular periodic orbits in or parallel to the equatorial plane. We divide all these orbits into six categories based on their appearances. By calculating the characteristic multipliers of the orbits, we investigate the stabilities of these periodic orbits.     

Simple three-dimensional periodic orbits in a galactic-type potential

We study the families of simple periodic orbits in a three-dimensional system that represents the inner parts of a perturbed triaxial galaxy. The perturbations depend on two control parameters. We find the regions where each family is stable, simply unstable, doubly unstable, or complex unstable. the stable and simply unstable families produce other families by bifurcation. Several families reach a maximum (or minimum) perturbation and then are continued by other families. The bifurcations are direct or inverse. The transition from one type of bifurcation to the other is theoretically explained. Another important phenomenon is the splitting of one family into two, or the joining of two families into one. We do not have any complex instability in the limiting cases of two-dimensional motions (when one control parameter is zero).The two main families of periodic orbits are in most cases stable when the energy is smaller than the escape energy. Most high energy orbits are unstable. However, we found stable orbits even for energies about four times larger than the escape energy.     

Is the third integral a function of the Hamiltonian?

Numerical evidence is presented which indicates that, although the third integral is tangent to the Hamiltonian (energy integral) along some periodic orbits (as has been shown by Goudas), it is not tangent to it along non-periodic orbits; therefore it is not a function of the Hamiltonian. The set of periodic orbits is probably dense in general, but a given form of the third integral is valid in the neighbourhood of a limited number of them; no form of the third integral is valid for all periodic orbits, except in integrable cases.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

Periodic Orbits in Rotating Second Degree and Order Gravity Fields

Periodic orbits in an arbitrary 2nd degree and order uniformly rotating gravity field are studied. We investigate the four equilibrium points in this gravity field. We see that close relation exists between the stability of these equilibria and the existence and stability of their nearby periodic orbits. We check the periodic orbits with non-zero periods. In our searching procedure for these periodic orbits, we remove the two unity eigenvalues from the state transition matrix to find a robust, non-singular linear map to solve for the periodic orbits. The algorithm converges well, especially for stable periodic orbits. Using the searching procedure, which is relatively automatic, we find five basic families of periodic orbits in the rotating second degree and order gravity field for planar motion, and discuss their existence and stability at different central body rotation rates.     

The rectilinear three-body problem using symbol sequence II: role of the periodic orbits

We study the change of phase space structure of the rectilinear three-body problem when the mass combination is changed. Generally, periodic orbits bifurcate from the stable Schubart periodic orbit and move radially outward. Among these periodic orbits there are dominant periodic orbits having rotation number (n − 2)/n with n ≥ 3. We find that the number of dominant periodic orbits is two when n is odd and four when n is even. Dominant periodic orbits have large stable regions in and out of the stability region of the Schubart orbit (Schubart region), and so they determine the size of the Schubart region and influence the structure of the Poincaré section out of the Schubart region. Indeed, with the movement of the dominant periodic orbits, part of complicated structure of the Poincaré section follows these orbits. We find stable periodic orbits which do not bifurcate from the Schubart orbit.     

The evolution of periodic orbits close to homoclinic points

For conservative dynamical systems having two degrees of freedom Birkhoff has established the existence of two classes of periodic orbits. The first consists of stable-unstable pairs close to periodic orbits of the stable type, and the second of orbits having fixed points (in a suitable surface of section) close to homoclinic points. In this paper orbits of the latter type are listed, and their evolution followed as a function of the energy. For the energy at which they were first computed, all were unstable; but they evolved, with diminishing energy, into one orbit of the stable type which appears to be a member of the first class of orbits mentioned above.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 27–September 2, 1972.     

2/1 resonant periodic orbits in three dimensional planetary systems

We consider the general spatial three body problem and study the dynamics of planetary systems consisting of a star and two planets which evolve into 2/1 mean motion resonance and into inclined orbits. Our study is focused on the periodic orbits of the system given in a suitable rotating frame. The stability of periodic orbits characterize the evolution of any planetary system with initial conditions in their vicinity. Stable periodic orbits are associated with long term regular evolution, while unstable periodic orbits are surrounded by regions of chaotic motion. We compute many families of symmetric periodic orbits by applying two schemes of analytical continuation. In the first scheme, we start from the 2/1 (or 1/2) resonant periodic orbits of the restricted problem and in the second scheme, we start from vertical critical periodic orbits of the general planar problem. Most of the periodic orbits are unstable, but many stable periodic orbits have been, also, found with mutual inclination up to 50^?–60^?, which may be related with the existence of real planetary systems.     

The evolution of periodic orbits close to heteroclinic points

Several families of periodic orbits close to heteroclinic points were computed, and their evolution, as a function of the energy, was followed. All appeared spontaneously at the occasion of the formation of the heteroclinic points, and had no genealogical links with other families. This behavior, which confirms a prediction by Dr. Contopoulos, is in contrast to the behavior of periodic orbits close to homoclinic points, some of which have genealogical links with periodic orbits existing in the ‘quasi-integrable’ phase of the system. The orbits found here appeared in regions where stochasticity was well established; so their appearance does not seem to be connected with the onset of stochasticity.     

Asymmetric Periodic Solutions in the Restricted Problem of Three Bodies

A systematic numerical exploration of the families of asymmetric periodic orbits of the restricted three-body problem when a) the primary bodies are equal and b) for the Earth-Moon mass ratio, is presented. Decades families of asymmetric periodic solutions were found and three of the simplest ones, in the first case, and ten of the second one are illustrated. All of these families consist of periodic orbits which are asymmetric with respect to x-axis while are simple symmetric periodic orbits with respect to y-axis (i.e. the orbit has only one perpendicular intersection at half period with y-axis). Many asymmetric periodic orbits, members of these families, are calculated and plotted. We studied the stability of all the asymmetric periodic orbits we found. These families consist, mainly, of unstable periodic solutions but there exist very small, with respect to x, intervals where these families have stable periodic orbits. We also found, using appropriate Poincaré surface of sections, that a relatively large region of phase space extended around all these stable asymmetric periodic orbits shows chaotic motion.     

A new family of periodic orbits for the restricted problem

A new family of periodic orbits of the three-dimensional restricted three-body problem which continue off from a consecutive collision orbit are numerically studied. Their behavior for varying energy is unexpected. In particular, associated with our system is a countable set of resonant energy values and each time the energy passes through one of them the periodic orbit forms a loop by self-intersection. Any number of loops can form by this process and the resulting orbits take on an interesting appearance.     

Symmetric and Nonsymmetric Periodic Orbits in the Exterior Mean Motion Resonances with Neptune

We present families of periodic orbits and their stability for the exterior mean motion resonances 1:2, 1:3 and 1:4 with Neptune in the framework of the planar circular restricted three-body problem. We found that in each resonance there exist two branches of symmetric elliptic periodic orbits with stable and unstable segments. Asymmetric periodic orbits bifurcate from the corresponding symmetric ones. Asymmetric periodic orbits are stable and the motion in their neighbourhood is a libration with respect to the resonant angle variable. In all the families of asymmetric periodic orbits the eccentricity extends to high values. Poincaré sections reveal the changes of the topology in phase space.