A Hill Problem with Oblate Primaries and Effect of Oblateness on Hill Stability of Orbits

We introduce a new version of Hill's problem to include the effect of oblateness of the primaries, and briefly discuss its equilibrium points and zero velocity curves. As a first application we use this to study Hill stability of direct orbits around the small primary. This can be employed to study the stability of a planet's moon perturbed by an oblate Sun, or of a star's planet perturbed by a distant disk-shaped galaxy. Oblateness of the `Sun' is found to decrese the maximum distance of Hill stable direct `moon' orbits. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Photogravitational Hill Problem: Numerical Exploration

We consider the recently introduced version of the classical Lunar Hill problem, the photogravitational Hill problem, and study it's equilibrium points and zero-velocity curves. The full network of families of periodic orbits is numerically explored, their stability is computed and critical orbits are determined. Non-periodic orbits are also computed as points on a surface of section, providing an outlook of the stability regions, chaotic motions and escape.     

The planar Hill problem with oblate primary

The regularized equations of motion of the planar Hill problem which includes the effect of the oblateness of the larger primary body, is presented. Using the Levi-Civita coordinate transformation as well as the corresponding time transformation, we obtain a simple regularized polynomial Hamiltonian of the dynamical system that corresponds to that of two uncoupled harmonic oscillators perturbed by polynomial terms. The relations between the synodic and regularized variables are also given. The convenient numerical computations of the regularized equations of motion, allow derivation of a map of the group of families of simple-periodic orbits, free of collision cases, of both the classical and the Hill problem with oblateness. The horizontal stability of the families is calculated and we determine series of horizontally critical symmetric periodic orbits of the basic families g and g'. This revised version was published online in July 2006 with corrections to the Cover Date.     

Hill stability of satellite orbits

Szebehely's criterion for Hill stability of satellites is derived from Hill's problem and a more exact result is obtained. Direct, Hill stable, circular satellites can exist almost twice as far from the planet as retrograde satellites. For direct satellites the new result agrees with Kuiper's empirical estimate that such satellites are stable up to a distance of half the radius of action of the planet. Comparison with the results of numerical experiments shows that Hill 'stability is valid for direct satellites but meaningless for retrograde satellites. Further accuracy for the maximum distance of Hill stable orbits is obtained from the restricted problem formulation. This provides estimates for planetary distances in double star systems.     

The dynamics of the elliptic Hill problem: periodic orbits and stability regions

The motion of a satellite around a planet can be studied by the Hill model, which is a modification of the restricted three body problem pertaining to motion of a satellite around a planet. Although the dynamics of the circular Hill model has been extensively studied in the literature, only few results about the dynamics of the elliptic model were known up to now, namely the equations of motion and few unstable families of periodic orbits. In the present study we extend these results by computing a large set of families of periodic orbits and their linear stability and classify them according to their resonance condition. Although most of them are unstable, we were able to find a considerable number of stable ones. By computing appropriate maps of dynamical stability, we study the effect of the planetary eccentricity on the stability of satellite orbits. We see that, even for large values of the planetary eccentricity, regular orbits can be found in the vicinity of stable periodic orbits. The majority of irregular orbits are escape orbits.     

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.     

Hill stability of the general problem of three bodies

We discuss Hill stability in the general three-body problem. The Hill curves in the general problem are the same as in the planar problem. We show that the bifurcation points correspond to the five equilibrium solutions, and derive the criterion for stability in the general case. Application of this criterion to 19 natural satellites of the Solar system leads to the result that, apart from Neptune 1, all the other 18 satellites are unstable in the sense of Hill. The dominant factor in producing this result is the finite eccentricity of the planetary orbits around the Sun.     

On the stability of particular solutions of the singly averaged Hill problem

We consider the particular solutions of the evolutionary system of equations in elements that correspond to planar and spatial circular orbits of the singly averaged Hill problem. We analyze the stability of planar and spatial circular orbits to inclination and eccentricity, respectively. We construct the instability regions of both particular solutions in the plane of parameters of the problem.     

Hill stability of the satellite in the elliptic restricted four-body problem

We consider an elliptic restricted four-body system including three primaries and a massless particle. The orbits of the primaries are elliptic, and the massless particle moves under the mutual gravitational attraction. From the dynamic equations, a quasi-integral is obtained, which is similar to the Jacobi integral in the circular restricted three-body problem (CRTBP). The energy constant \(C\) determines the topology of zero velocity surfaces, which bifurcate at the equilibrium point. We define the concept of Hill stability in this problem, and a criterion for stability is deduced. If the actual energy constant \(C_{\mathrm{ac}}\ ( {>} 0 ) \) is bigger than or equal to the critical energy constant \(C_{\mathrm{cr}}\), the particle will be Hill stable. The critical energy constant is determined by the mass and orbits of the primaries. The criterion provides a way to capture an asteroid into the Earth–Moon system.     

On the periodically evolving orbits in the singly averaged Hill problem

We continue to analyze the periodic solutions of the singly averaged Hill problem. We have numerically constructed the families of solutions that correspond to periodically evolving satellite orbits for arbitrary initial values of their eccentricities and inclinations to the plane of motion of the perturbing body. The solutions obtained are compared with the numerical solutions of the rigorous (nonaveraged) equations of the restricted circular three-body problem. In particular, we have constructed a periodically evolving orbit for which the well-known Lidov-Kozai mechanism manifests itself, just as in the doubly averaged problem.     

The Hill stability of inclined small mass binary systems in three-body systems with special application to triple star systems, extrasolar planetary systems and Binary Kuiper Belt systems

The dynamical stability of a bound triple system composed of a small binary or minor planetary system moving on a orbit inclined to a central third body is discussed in terms of Hill stability for the full three-body problem. The situation arises in the determination of stability of triple star systems against disruption and component exchange and the determination of stability of extrasolar planetary systems and minor planetary systems against disruption, component exchange or capture. The Hill stability criterion is applied to triple star systems and extrasolar planetary systems, the Sun-Earth-Moon system and Kuiper Belt binary systems to determine the critical distances for stable orbits. It is found that increasing the inclination of the third body decreases the Hill regions of stability. Increasing the eccentricity of the binary also produces similar effects.These type of changes make exchange or disruption of the component masses more likely. Increasing the eccentricity of the binary orbit relative to the third body substantially decreases stability regions as the eccentricity reaches higher values. The Kuiper Belt binaries were found to be stable if they move on circular orbits. Taking into account the eccentricity, it is less clear that all the systems are stable.     

Stability of periodic solutions for Hill’s averaged problem with allowance for planetary oblateness

We analyze the stability of periodic solutions for Hill’s double-averaged problem by taking into account a central planet’s oblateness. They are generated by steady-state solutions that are stable in the linear approximation. By numerically calculating the monodromy matrix of variational equations, we plot its trace against the integral of the problem—an averaged perturbing function, for two model systems, [(Sun + Moon)-Earth-satellite] and (Sun-Uranus-satellite). We roughly estimate the ranges of values for the parameters of satellite orbits corresponding to periodic solutions of the evolutionary system that are stable in the linear approximation.     

THE HILL PROBLEM WITH OBLATE SECONDARY: NUMERICAL EXPLORATION

We introduce a three-dimensional version of Hill’s problem with oblate secondary, determine its equilibrium points and their stability and explore numerically its network of families of simple periodic orbits in the plane, paying special attention to the evolution of this network for increasing oblateness of the secondary. We obtain some interesting results that differentiate this from the classical problem. Among these is the eventual disappearance of the basic family g′ of the classical Hill problem and the existence of out-of-plane equilibrium points and a family of simple-periodic plane orbits non-symmetric with respect to the x-axis.     

Energy accumulating accretion disk orbiting a Kerr black hole

A geometrically thin, energy accumulating $aL-disk is suggested which orbits a Kerr black hole. With increasing internal forces, the “standard” disks develop into energy accumulating disks. These accumulating disks are geometrically thin as long as their internal forces remain below a certain bound, allowing nearly geodesic orbits.     

On the Coplanar Integrable Case of the Twice-Averaged Hill Problem with Central Body Oblateness

The twice-averaged Hill problem with the oblateness of the central planet is considered in the case where its equatorial plane coincides with the plane of its orbital motion relative to the perturbing body. A qualitative study of this so-called coplanar integrable case was begun by Y. Kozai in 1963 and continued by M.L. Lidov and M.V. Yarskaya in 1974. However, no rigorous analytical solution of the problem can be obtained due to the complexity of the integrals. In this paper we obtain some quantitative evolution characteristics and propose an approximate constructive-analytical solution of the evolution system in the form of explicit time dependences of satellite orbit elements. The methodical accuracy has been estimated for several orbits of artificial lunar satellites by comparison with the numerical solution of the evolution system.     

Dynamical Effects of Mars on Asteroidal Dust Particles

Asteroidal dust particles resulting from family-forming events migrate from their source locations in the asteroid belt inwards towards the Sun under the effect of Poynting-Robertson (PR) drag. Understanding the distribution of these dust particle orbits in the inner solar system is of great importance to determining the asteroidal contribution to the zodiacal cloud, the accretion rate by the Earth, and the threat that these particles pose to spacecraft and satellites in near-Earth space. In order to correctly describe this distribution of orbits in the inner solar system, we must track the dynamical perturbations that the dust particle orbits experience as they migrate inwards. In a seminal paper Öpik (1951) determines that very few of the μm-cm sized dust particles suffer a collision with the planet face as they decay inwards past Mars. Here we re-analyze this problem, considering additionally the likelihood that the dust particle orbits pass through the Hill sphere of Mars (to various depths) and experience potentially significant perturbations to their orbits. We find that a considerable fraction of dust particle orbits will enter the Hill sphere of Mars. Furthermore, we find that there is a bias with inclination, particle size, and eccentricity of the particle orbits that enter the Martian Hill sphere. In particular the bias with inclination may create a bias towards higher-inclination sources in the proportions of asteroid family particles that reach near-Earth space.     

Hill stability in the many-body problem

We established a criterion for the Hill stability of motions in the problem of many spherical bodies with a spherical density distribution. The region of Hill stability was determined. The sizes of this region are comparable to the total volume of all of the bodies in the system, which sharply increases the probability of mutual collisions. This result may be considered as a confirmation that a supermassive core can be formed at the center of a globular star cluster. The motions in the n-body problem are shown to be unstable according to Hill.     

The photogravitational Hill problem with oblateness: equilibrium points and Lyapunov families

We introduce a new version of Hill’s problem that incorporates the effects of radiation of the primary and oblateness of the secondary and study the basic dynamical features of this new model-problem. This formulation is more appropriate for some astronomical applications as an approximation to the corresponding restricted three-body problem. We use iterative methods for deriving approximate expressions of the equilibrium point locations and study their stability properties by using a linear stability analysis. All equilibrium points are unstable. We also employ singular perturbations methods for obtaining approximate expressions of the Lyapunov families emanating from equilibrium points, in both coplanar and spatial case, and numerical techniques for their continuation.     

On the Dynamical Foundations of the Lidov–Kozai Theory

The Lidov–Kozai theory developed by each of the authors independently in 1961–1962 is based on qualitative methods of studying the evolution of orbits for the satellite version of the restricted three-body problem (Hill’s problem). At present, this theory is in demand in various fields of science: in the field of planetary research within the Solar system, the field of exoplanetary systems, and the field of high-energy physics in interstellar and intergalactic space. This has prompted me to popularize the ideas that underlie the Lidov–Kozai theory based on the experience of using this theory as an efficient tool for solving various problems related to the study of the secular evolution of the orbits of artificial planetary satellites under the influence of external gravitational perturbations with allowance made for the perturbations due to the polar planetary oblateness.