On The Concept Of Periapsis In Hill’s Problem

The concept of closest approach is analyzed in Hill’s problem, resulting in a partitioning of the position space. The different behavior between the direct and retrograde motion is explained analytically, resulting in a simple estimate of the variation of Hill’s periodic and quasi-circular orbits as a function of the Jacobi constant. The local behavior of the orbits on the zero velocity surfaces and an analytical definition of local escape and capture in Hill’s problem are also given.     

Radar observations of Itokawa in 2004 and improved shape estimation

Abstract— We present June 2004 radar images of asteroid 25143 Itokawa (1998 SF36) that improve upon the longitude‐latitude coverage of images obtained in 2001 by Ostro et al. (2004) and use the 2001–2004 data to refine that paper's constraints on Itokawa's shape. The 2004 images, the first of the asteroid's southern side, look distinctly different from the 2001 images, revealing leading edges that are much more curved and rugged than the nearly convex leading edges seen at northern latitudes in 2001. Itokawa is shaped like a slightly asymmetrical, bent, lumpy ellipsoid with dimensions along the principal axes within 10% of 594 times 320 times 288 m. To illustrate the uncertainty space associated with shape reconstruction from images with suboptimal orientational coverage, we present two alternative three‐dimensional models of the object.     

The Restricted Hill Four-Body Problem with Applications to the Earth–Moon–Sun System

Starting from the four-body problem a generalization of both the restricted three-body problem and the Hill three-body problem is derived. The model is time periodic and contains two parameters: the mass ratio ν of the restricted three-body problem and the period parameter m of the Hill Variation orbit. In the proper coordinate frames the restricted three-body problem is recovered as m → 0 and the classical Hill three-body problem is recovered as ν → 0. This model also predicts motions described by earlier researchers using specific models of the Earth–Moon–Sun system. An application of the current model to the motion of a spacecraft in the Sun perturbed Earth–Moon system is made using Hill's Variation orbit for the motion of the Earth–Moon system. The model is general enough to apply to the motion of an infinitesimal mass under the influence of any two primaries which orbit a larger mass. Using the model, numerical investigations of the structure of motions around the geometric position of the triangular Lagrange points are performed. Values of the parameter ν range in the neighborhood of the Earth–Moon value as the parameter m increases from 0 to 0.195 at which point the Hill Variation orbit becomes unstable. Two families of planar periodic orbits are studied in detail as the parameters m and ν vary. These families contain stable and unstable members in the plane and all have the out-of-plane stability. The stable and unstable manifolds of the unstable periodic orbits are computed and found to be trapped in a geometric area of phase space over long periods of time for ranges of the parameter values including the Earth–Moon–Sun system. This model is derived from the general four-body problem by rigorous application of the Hill and restricted approximations. The validity of the Hill approximation is discussed in light of the actual geometry of the Earth–Moon–Sun system. This revised version was published online in July 2006 with corrections to the Cover Date.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation when the two planets are in migration

In a previous paper (Hou et al. in Celest Mech Dyn Astron 119:119–142, 2014a), the problem of dynamical symmetry between two Jupiter triangular libration points (TLPs) with Saturn’s perturbation in the present configuration of the two planets was studied. A small short-time scale spatial asymmetry exists but gradually disappears with the time going, so the planar stable regions around the two Jupiter TLPs should be dynamically symmetric from a longtime perspective. In this paper, the symmetry problem is studied when the two planets are in migration. Several mechanisms that can cause asymmetries are discussed. Studies show that three important ones are the large short-time scale spatial asymmetry when Jupiter and Saturn are in resonance, the changing orbits of Jupiter and Saturn in the planet migration process, and the chaotic nature of Trojan orbits during the planet migration process. Their joint effects can cause an observable difference to the two Jupiter Trojan swarms. The thermal Yarkovsky effect is also found to be able to cause dynamical differences to the two TLPs, but generally they are too small to be practically observed.     

The use of invariant manifolds for transfers between unstable periodic orbits of different energies

Techniques from dynamical systems theory have been applied to the construction of transfers between unstable periodic orbits that have different energies. Invariant manifolds, trajectories that asymptotically depart or approach unstable periodic orbits, are used to connect the initial and final orbits. The transfer asymptotically departs the initial orbit on a trajectory contained within the initial orbit’s unstable manifold and later asymptotically approaches the final orbit on a trajectory contained within the stable manifold of the final orbit. The manifold trajectories are connected by the execution of impulsive maneuvers. Two-body parameters dictate the selection of the individual manifold trajectories used to construct efficient transfers. A bounding sphere centered on the secondary, with a radius less than the sphere of influence of the secondary, is used to study the manifold trajectories. A two-body parameter, κ, is computed within the bounding sphere, where the gravitational effects of the secondary dominate. The parameter κ is defined as the sum of two quantities: the difference in the normalized angular momentum vectors and eccentricity vectors between a point on the unstable manifold and a point on the stable manifold. It is numerically demonstrated that as the κ parameter decreases, the total cost to complete the transfer decreases. Preliminary results indicate that this method of constructing transfers produces a significant cost savings over methods that do not employ the use of invariant manifolds.     

Mutual Potential of Homogeneous Polyhedra

The mutual gravitational potential between a pair of homogeneous polyhedra is expressed using an infinite series. The nested volume integrals are evaluated analytically and result in simple tensor expressions containing no special functions. However, complexity increases as O(6 ⁿ ), where n is the term degree. An alternate formulation due to Liebenthal is also presented.     

The OSIRIS‐REx target asteroid (101955) Bennu: Constraints on its physical,geological, and dynamical nature from astronomical observations

We review the results of an extensive campaign to determine the physical, geological, and dynamical properties of asteroid (101955) Bennu. This investigation provides information on the orbit, shape, mass, rotation state, radar response, photometric, spectroscopic, thermal, regolith, and environmental properties of Bennu. We combine these data with cosmochemical and dynamical models to develop a hypothetical timeline for Bennu's formation and evolution. We infer that Bennu is an ancient object that has witnessed over 4.5 Gyr of solar system history. Its chemistry and mineralogy were established within the first 10 Myr of the solar system. It likely originated as a discrete asteroid in the inner Main Belt approximately 0.7–2 Gyr ago as a fragment from the catastrophic disruption of a large (approximately 100‐km), carbonaceous asteroid. It was delivered to near‐Earth space via a combination of Yarkovsky‐induced drift and interaction with giant‐planet resonances. During its journey, YORP processes and planetary close encounters modified Bennu's spin state, potentially reshaping and resurfacing the asteroid. We also review work on Bennu's future dynamical evolution and constrain its ultimate fate. It is one of the most Potentially Hazardous Asteroids with an approximately 1‐in‐2700 chance of impacting the Earth in the late 22nd century. It will most likely end its dynamical life by falling into the Sun. The highest probability for a planetary impact is with Venus, followed by the Earth. There is a chance that Bennu will be ejected from the inner solar system after a close encounter with Jupiter. OSIRIS‐REx will return samples from the surface of this intriguing asteroid in September 2023.     

Dynamics of the Jupiter Trojans with Saturn’s perturbation in the present configuration of the two planets

The dynamics of the two Jupiter triangular libration points perturbed by Saturn is studied in this paper. Unlike some previous works that studied the same problem via the pure numerical approach, this study is done in a semianalytic way. Using a literal solution, we are able to explain the asymmetry of two orbits around the two libration points with symmetric initial conditions. The literal solution consists of many frequencies. The amplitudes of each frequency are the same for both libration points, but the initial phase angles are different. This difference causes a temporary spatial asymmetry in the motions around the two points, but this asymmetry gradually disappears when the time goes to infinity. The results show that the two Jupiter triangular libration points should have symmetric spatial stable regions in the present status of Jupiter and Saturn. As a test of the literal solution, we study the resonances that have been extensively studied in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006). The resonance structures predicted by our analytic theory agree well with those found in Robutel and Gabern (Mon Not R Astron Soc 372:1463–1482, 2006) via a numerical approach. Two kinds of chaotic orbits are discussed. They have different behaviors in the frequency map. The first kind of chaotic orbits (inner chaotic orbits) is of small to moderate amplitudes, while the second kind of chaotic orbits (outer chaotic orbits) is of relatively larger amplitudes. Using analytical theory, we qualitatively explain the transition process from the inner chaotic orbits to the outer chaotic orbits with increasing amplitudes. A critical value of the diffusion rate is given to separate them in the frequency map. In a forthcoming paper, we will study the same problem but keep the planets in migration. The time asymmetry, which is unimportant in this paper, may cause an observable difference in the two Jupiter Trojan groups during a very fast planet migration process.     

Small body surface gravity fields via spherical harmonic expansions

Conventional gravity field expressions are derived from Laplace’s equation, the result being the spherical harmonic gravity field. This gravity field is said to be the exterior spherical harmonic gravity field, as its convergence region is outside the Brillouin (i.e., circumscribing) sphere of the body. In contrast, there exists its counterpart called the interior spherical harmonic gravity field for which the convergence region lies within the interior Brillouin sphere that is not the same as the exterior Brillouin sphere. Thus, the exterior spherical harmonic gravity field cannot model the gravitation within the exterior Brillouin sphere except in some special cases, and the interior spherical harmonic gravity field cannot model the gravitation outside the interior Brillouin sphere. In this paper, we will discuss two types of other spherical harmonic gravity fields that bridge the null space of the exterior/interior gravity field expressions by solving Poisson’s equation. These two gravity fields are obtained by assuming the form of Helmholtz’s equation to Poisson’s equation. This method renders the gravitational potentials as functions of spherical Bessel functions and spherical harmonic coefficients. We refer to these gravity fields as the interior/exterior spherical Bessel gravity fields and study their characteristics. The interior spherical Bessel gravity field is investigated in detail for proximity operation purposes around small primitive bodies. Particularly, we apply the theory to asteroids Bennu (formerly 1999 RQ36) and Castalia to quantify its performance around both nearly spheroidal and contact-binary asteroids, respectively. Furthermore, comparisons between the exterior gravity field, interior gravity field, interior spherical Bessel gravity field, and polyhedral gravity field are made and recommendations are given in order to aid planning of proximity operations for future small body missions.     

The Influence of Reactive Torques on Comet Nucleus Rotation

Reactive torques, due to anisotropic sublimation on a comet nucleus surface, produce slow variations of its rotation. In this paper the secular effects of this sublimation are studied. The general rotational equations of motion are averaged over unperturbed fast rotation around the mass center (Euler-Poinsot motion) and over the orbital comet motion. We discuss the parameters that define typical properties of the rotational evolution and discover different classifications of the rotational evolution. As an example we discuss some possible scenarios of rotational evolution for the nuclei of the comets Halley and Borrelly.