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.     

Metrics in Keplerian orbits quotient spaces

Quotient spaces of Keplerian orbits are important instruments for the modelling of orbit samples of celestial bodies on a large time span. We suppose that variations of the orbital eccentricities, inclinations and semi-major axes remain sufficiently small, while arbitrary perturbations are allowed for the arguments of pericentres or longitudes of the nodes, or both. The distance between orbits or their images in quotient spaces serves as a numerical criterion for such problems of Celestial Mechanics as search for common origin of meteoroid streams, comets, and asteroids, asteroid families identification, and others. In this paper, we consider quotient sets of the non-rectilinear Keplerian orbits space \(\mathbb H\). Their elements are identified irrespective of the values of pericentre arguments or node longitudes. We prove that distance functions on the quotient sets, introduced in Kholshevnikov et al. (Mon Not R Astron Soc 462:2275–2283, 2016), satisfy metric space axioms and discuss theoretical and practical importance of this result. Isometric embeddings of the quotient spaces into \(\mathbb R^n\), and a space of compact subsets of \(\mathbb H\) with Hausdorff metric are constructed. The Euclidean representations of the orbits spaces find its applications in a problem of orbit averaging and computational algorithms specific to Euclidean space. We also explore completions of \(\mathbb H\) and its quotient spaces with respect to corresponding metrics and establish a relation between elements of the extended spaces and rectilinear trajectories. Distance between an orbit and subsets of elliptic and hyperbolic orbits is calculated. This quantity provides an upper bound for the metric value in a problem of close orbits identification. Finally the invariance of the equivalence relations in \(\mathbb H\) under coordinates change is discussed.     

Distance in the space of energetically bounded Keplerian orbits

In this paper we derive an explicit, analytic formula for the geodesic distance between two points in the space of bounded Keplerian orbits in a particular topology. The specific topology we use is that of a cone passing through the direct product of two spheres. The two spheres constitute the configuration manifold for the space of bounded orbits of constant energy. We scale these spheres by a factor equal to the semi-major axis of the orbit, forming a linear cone. This five-dimensional manifold inherits a Riemannian metric, which is induced from the Euclidean metric on \mathbbR⁶{\mathbb{R}^6}, the space in which it is embedded. We derive an explicit formula for the geodesic distance between any two points in this space, each point representing a physical, gravitationally bound Keplerian orbit. Finally we derive an expression for the Riemannian metric that we used in terms of classical orbital elements, which may be thought of as local coordinates on our configuration manifold.     

Symplectic Mapping for Satellites and Space Debris Including Nongravitational Forces

The paper presents an efficient algorithm for the study of satellite and space debris orbits on long time intervals. The averaged equations of motion are integrated by means of the implicit midpoint method. This approach is known as a symplectic mapping technique. The perturbing forces included in the mapping are: the geopotential, the atmospheric drag, lunisolar perturbations and the direct radiation pressure (without shadow effects). The influence of the atmosphere is approximated by simple methods for the estimation of integrals. The described mapping is valid for the wide range of orbits including the resonant and the eccentric ones; it can be helpful in practical and theoretical problems. The lifetime of GPS transfer orbits is discussed as an exemplary application.     

Orbits and manifolds near the equilibrium points around a rotating asteroid

We study the orbits and manifolds near the equilibrium points of a rotating asteroid. The linearised equations of motion relative to the equilibrium points in the gravitational field of a rotating asteroid, the characteristic equation and the stable conditions of the equilibrium points are derived and discussed. First, a new metric is presented to link the orbit and the geodesic of the smooth manifold. Then, using the eigenvalues of the characteristic equation, the equilibrium points are classified into 8 cases. A theorem is presented and proved to describe the structure of the submanifold as well as the stable and unstable behaviours of a massless test particle near the equilibrium points. The linearly stable, the non-resonant unstable, and the resonant equilibrium points are discussed. There are three families of periodic orbits and four families of quasi-periodic orbits near the linearly stable equilibrium point. For the non-resonant unstable equilibrium points, there are four relevant cases; for the periodic orbit and the quasi-periodic orbit, the structures of the submanifold and the subspace near the equilibrium points are studied for each case. For the resonant equilibrium points, the dimension of the resonant manifold is greater than 4, and we find at least one family of periodic orbits near the resonant equilibrium points. As an application of the theory developed here, we study relevant orbits for the asteroids 216 Kleopatra, 1620 Geographos, 4769 Castalia and 6489 Golevka.     

Natural Metrics in the Spaces of Elliptic Orbits

Different natural metrizations by Hölder type on the five dimensional space of Keplerian elliptic orbits are introduced. Certain applications of topological and metrical properties of the space of Keplerian elliptic orbits to several problems of Celestial Mechanics are discussed.     

银河系中球状星团的空间分布和运动轨道

选用银河系中29个累积光谱型为F型的球状星团样本。根据它们的视向速度，绝对自行等参数，归算处理后得出了各样本星团的空间分布和运动速度。并以此作为初始条件，在给定的3种银河系引力势模型中，采用数值积分方法计算出各样本星团的运动轨道。计算结果表明：（1）大部分样本星团都位于银心距5kpc-10kpc的范围内，相对于银心呈球对称分布，它们的速度也呈椭球分布；（2）29个样本星团按其金属度大小和基本性发类，可分属HB和MP两个次系，且样本星团数随金属度[Fe/H]而变化，在[Fe/H]＝－1．6处出现一个峰值；（3）所有样本星团的轨道运动都呈周期性，大都在一个有界而不封闭的周期轨道上运动，其最大银心距大都在40kpc以内。不同的引力势模型对球状星团轨道的具体形态影响不大，在给定的引力势模型下，当某些星团的运动轨道穿越距银心1kpc附近的区域时会出现“混沌”行为。而样本星团的金属度与其轨道形态之间的相关性并不明显；（4）29个样本星团的轨道半长轴、远银心距和方位周期随金属度的变化规律基本相似。轨道偏心率与金属度有关，对于所选的晕族样本星团而言，大约有24％的样本星团的轨道偏心率低于0．4，不同的引力势模型对近银心距、偏心率和参数的不确定度等量影响较小，但是对远银心距、径向周期和方位周期等参数影响较为明显。     

Families of double symmetric “Schubart-like” periodic orbits

We consider a class of Hamiltonian systems with two degrees of freedom with singularities. This class includes several symmetric subproblems of the $n$ -body problem where the singularities are due to collisions involving two or more bodies. “Schubart-like” periodic orbits having two collisions in one period, are present in most of these subproblems. The purpose of this paper is to study the existence of families of such a periodic orbits in a general setting. The blow up techniques of total collision and infinity are applied to our class of Hamiltonian system. This allows us to derive sufficient conditions to ensure the existence of families of double symmetric “Schubart-like” periodic orbits having many singularities. The orbits in the family can be parametrized by the number of singularities in one period. The results are applied to some subproblems of the gravitational $n$ -body problem.     

Predicting the maximum sunspot number and the associated geomagnetic activity indices a a $aa$ and A p $Ap$ for solar cycle 25

We are considering the spacetime described by the metric proposed by Mannheim and Kazanas. The effective potential and the circular orbits are discussed. The rotational velocity derived from the geodesics equation agrees with the observed flat galactic rotation curves. Finally, solutions to the Gordon equation for massless bosons evolving in this spacetime are obtained in terms of Heun general functions.

     

Structures in the phase space of a four dimensional symplectic map

We numerically investigate the projections of non periodic orbits in a 4-dimensional (4-D) symplectic map composed of two coupled 2-dimensional (2-D) maps. We describe in detail the structures that are produced in different planes of projection and we find how the morphology of the 4-D orbits is influenced by the features of the 2-D maps as the coupling parameter increases. We give an empirical law that describes this influence.     

Nodal and periastron precession of inclined orbits in the field of a rapidly rotating neutron star

We derive a formula for the nodal precession frequency and the Keplerian period of a particle at an arbitrary orbital inclination (with a minimum latitudinal angle reached at the orbit) in the post-Newtonian approximation in the external field of an oblate rotating neutron star (NS). We also derive formulas for the nodal precession and periastron rotation frequencies of slightly inclined low-eccentricity orbits in the field of a rapidly rotating NS in the form of asymptotic expansions whose first terms are given by the Okazaki-Kato formulas. The NS gravitational field is described by the exact solution of the Einstein equation that includes the NS quadrupole moment induced by rapid rotation. Convenient asymptotic formulas are given for the metric coefficients of the corresponding space-time in the form of Kerr metric perturbations in Boyer-Lindquist coordinates.     

关于共线平动点的特征及其在深空探测中的应用

系统阐述了小天体运动对应的圆型限制性三体问题共线平动点的强不稳定性特征,以及其附近的条件周期轨道——晕轨道(Halo Orbit)的存在、相应解的构造。这种特殊的轨道形式和共线平动点附近的弱稳定走廊,可分别用于在深空特殊位置附近定点有各种科学探测目标的探测器和向节能轨道过渡的通道。     

Fractal Dimensions as Chaos Indicators

The correlation dimension D ₂ is used to develop a method of classification for phase space orbits. D ₂ depends only on the mutual distances of the orbits points, therefore the time development of the orbit is reflected in the time development of the correlation dimension approximants. It is shown, that this technique allows to investigate the dynamical properties of a phase space orbit, a classification of chaotic and regular orbits and a detection of sticky orbits.     

Perihelion precession in binary systems: higher order corrections

Higher order corrections (up to n-th order) are obtained for the perihelion precession in binary systems like OJ287 using the Schwarzschild metric and complex integration. The corrections are performed considering the third root of the motion equation and developing the expansion in terms of \(r_{s}/ (a(1-e^{2}) )\).The results are compared with other expansions that appear in the literature giving corrections to second and third order. Finally, we simulate the shape of relativistic orbits for binary systems with different masses.     

Newtonian analogue of force and motion of a free particle in the gravitational field of Kerr-de Sitter space-time

In the first part of the paper the Newtonian analogue of force for the Kerr-de Sitter metric has been investigated. To the first-order of approximation, a component of the force vector corresponds to the Newtonian gravitational force and a cosmic force arising due to the cosmological constant A. Int the higher order of approximation, the relativistic correction terms due to rotation and the presence of A are obtained. In the second part of the paper the motion of a freely-falling body has been investigated. It is found that plane orbits are not possible. Also a radial fall is not possible and there is a rotational drag on the particle which has no Newtonian analogue.     

The warping of the galactic plane and changes of ellipticity

Observational evidence on the widespread occurrence of warping of the outer part of the galactic plane in many galaxies is presented and various hypotheses for its explanation are reviewed. None is found to be able to account for all the cases reported. Other phenomena considered are: (1) deviations from the galactic plane in the innermost nuclear region; (2) differences in orientation of the nucleus and disk in spiral galaxies; (3) changing ellipticity with distance from the center.After discussing the observations available, a theory has been developed which explains the phenomena mentioned as natural consequences of the non-steady nature of galactic systems due to time-dependent metric. This manifests itself in the appearance of tangential acceleration which leads necessarily to variability of the orbital plane and orbital eccentricity in dependence on the radius vector of the orbits.     

Hysteresis effects and diagnostics of the shock formation in low angular momentum axisymmetric accretion in the Kerr metric

The secular evolution of the purely general relativistic low angular momentum accretion flow around a spinning black hole is shown to exhibit hysteresis effects. This confirms that a stationary shock is an integral part of such an accretion disc in the Kerr metric. The equations describing the space gradient of the dynamical flow velocity of the accreting matter have been shown to be equivalent to a first order autonomous dynamical systems. Fixed point analysis ensures that such flow must be multi-transonic for certain astrophysically relevant initial boundary conditions. Contrary to the existing consensus in the literature, the critical points and the sonic points are proved not to be isomorphic in general, they can form in a completely different length scales. Physically acceptable global transonic solutions must produce odd number of critical points. Homoclinic orbits for the flow possessing multiple critical points select the critical point with the higher entropy accretion rate, confirming that the entropy accretion rate is the degeneracy removing agent in the system. However, heteroclinic orbits are also observed for some special situation, where both the saddle type critical points of the flow configuration possesses identical entropy accretion rate. Topologies with heteroclinic orbits are thus the only allowed non-removable degenerate solutions for accretion flow with multiple critical points, and are shown to be structurally unstable. Depending on suitable initial boundary conditions, a homoclinic trajectory can be combined with a standard non-homoclinic orbit through an energy preserving Rankine-Hugoniot type of stationary shock, and multi-critical accretion flow then becomes truly multi-transonic. An effective Lyapunov index has been proposed to analytically confirm why certain class of transonic flow cannot accommodate shock solutions even if it produces multiple critical points.     

Interplanetary radar time delays in different theories of gravitation

The first order (in the massm of the central body) time delay is deduced for all theories of gravitation that satisfy the 3-dimensional principle of equivalence in the one-body problem. It is shown that the curvature of the light ray causes no first order time delay but that the first order time delay is sensitive to the non-linear term of a general metric when distancesr are measured by the periods of circular planetary orbits.     

Natural and sail-displaced doubly-symmetric Lagrange point orbits for polar coverage

This paper proposes the use of doubly-symmetric, eight-shaped orbits in the circular restricted three-body problem for continuous coverage of the high-latitude regions of the Earth. These orbits, for a range of amplitudes, spend a large fraction of their period above either pole of the Earth. It is shown that they complement Sun-synchronous polar and highly eccentric Molniya orbits, and present a possible alternative to low thrust pole-sitter orbits. Both natural and solar-sail displaced orbits are considered. Continuation methods are described and used to generate families of these orbits. Starting from ballistic orbits, other families are created either by increasing the sail lightness number, varying the period or changing the sail attitude. Some representative orbits are then chosen to demonstrate the visibility of high-latitude regions throughout the year. A stability analysis is also performed, revealing that the orbits are unstable: it is found that for particular orbits, a solar sail can reduce their instability. A preliminary design of a linear quadratic regulator is presented as a solution to stabilize the system by using the solar sail only. Finally, invariant manifolds are exploited to identify orbits that present the opportunity of a ballistic transfer directly from low Earth orbit.