We find an interesting fact that fictitious retrograde co-orbitals of Saturn, or small bodies inside the retrograde 1:1 resonance with Saturn, are highly unstable in our numerical simulations. It is shown that in the presence of Jupiter, the retrograde co-orbitals will get ejected from Saturn's co-orbital space within a timescale of 10 Myr. This scenario reminds us of the instability of Saturn Trojans caused by both the Great Inequality and the secular resonances. Therefore, we carry out in-depth inspections on both mechanisms and prove that the retrograde resonance overlap, raised by Great Inequality, cannot serve as an explanation for the instability of retrograde co-orbitals, due to the weakness of the retrograde 2:5 resonance with Jupiter at a low eccentricity. However, we discover that both $\nu_5$ and $\nu_6$ secular resonances contribute to the slow growth of the eccentricity, therefore, are possibly the primary causes of the instability inside Saturn's retrograde co-orbital space.

In the preliminary trajectory design of the multi-target rendezvous problem, a model that can quickly estimate the cost of the orbital transfer is essential. The estimation of the transfer time using solar sail between two arbitrary orbits is difficult and usually requires to solve an optimal control problem. Inspired by the successful applications of the deep neural network in nonlinear regression, this work explores the possibility and effectiveness of mapping the transfer time for solar sail from the orbital characteristics using the deep neural network. Furthermore, the Monte Carlo Tree Search method is investigated and used to search the optimal sequence considering a multi-asteroid exploration problem. The obtained sequences from preliminary design will be solved and verified by sequentially solving the optimal control problem. Two examples of different application backgrounds validate the effectiveness of the proposed approach.

As the discoveries of more minor bodies in retrograde resonances with giant planets, such as 2015 BZ509 and 2006 RJ2, our curiosity about the Kozai-Lidov dynamics inside the retrograde resonance has been sparked. In this study, we focus on the 3D retrograde resonance problem and investigate how the resonant dynamics of a minor body impacts on its own Kozai-Lidov cycle. Firstly we deduce the action-angle variables and canonical transformations that deal with the retrograde orbit specifically. After obtaining the dominant Hamiltonian of this problem, we then carry out the numerical averaging process in closed form to generate phase-space portraits on a $e-\omega$ space. The retrograde 1:1 resonance is particularly scrutinized in detail, and numerical results from a CRTBP model shows a great agreement with the our semi-analytical portraits. On this basis, we inspect two real minor bodies currently trapped in retrograde 1:1 mean motion resonance. It is shown that they have different Kozai-Lidov states, which can be used to analyze the stability of their unique resonances. In the end, we further inspect the Kozai-Lidov dynamics inside the 2:1 and 2:5 retrograde resonance, and find distinct dynamical bifurcations of equilibrium points on phase-space portraits.

Aims. 2015 BZ509 is the first asteroid confirmed to be in retrograde co-orbit resonance (or 1/-1 resonance) with the giant planets in the solar system. While Saturn is the only giant planet whose Trojans are not discovered until now, we identify some small bodies among Centaurs and Damocloids that are potentially in 1/-1 resonance with Saturn in the present study. Methods. We integrate numerically the motion of the 1000 clones (include the nominal orbit) of each Centaur whose orbit has a semi-major axis between 9.3 au and 9.8 au and an inclination i > 90 deg. To confirm and evaluate the 1/-1 resonant configurations mentioned above, we introduce a useful one-degree integrable approximation for planar 1/-1 resonance. Results. We identify four candidates potentially in 1/-1 resonance with Saturn. The capture in this particular resonant state during the 40000 yr integration timespan is very common for 2006 RJ2 (906/1000 clones), 2006 BZ8 (878/1000 clones), and 2017 SV13 (998/1000 clones), and it is less likely for 2012 YE8 (426/1000 clones). According to our statistical results, 2006 RJ2 is the best candidate to be currently in a 1/-1 mean motion resonance with Saturn, and 2017 SV13 is another important potential candidate. Moreover, 2012 YE8 and 2006 BZ8 are also Centaurs of interest but their current and long-term 1/-1 resonant state with Saturn is less likely. The proportions of the clones captured in the relative long-term stable co-orbit resonance (over 10000 yr) are also given. Conclusions. Small bodies in retrograde co-orbit resonance with giant planets may be more common than previously expected. Identification of these potential mysterious minor bodies encourages the search for such objects on a larger scale in our solar system. The findings of this paper are also useful for understanding the origin and dynamical evolution of the Centaurs and Damocloids on retrograde orbits.

Asteroids in mean motion resonances with giant planets are common in the solar system, but it was not until recently that several asteroids in retrograde mean motion resonances with Jupiter and Saturn were discovered. A retrograde co-orbital asteroid of Jupiter, 2015 BZ509 is confirmed to be in a long-term stable retrograde 1:1 mean motion resonance with Jupiter, which gives rise to our interests in its unique resonant dynamics. In this paper, we investigate the phase-space structure of the retrograde 1:1 resonance in detail within the framework of the circular restricted three-body problem. We construct a simple integrable approximation for the planar retrograde resonance using canonical contact transformation and numerically employ the averaging procedure in closed form. The phase portrait of the retrograde 1:1 resonance is depicted with the level curves of the averaged Hamiltonian. We thoroughly analyze all possible librations in the co-orbital region and uncover a new apocentric libration for the retrograde 1:1 resonance inside the planet's orbit. We also observe the significant jumps in orbital elements for outer and inner apocentric librations, which are caused by close encounters with the perturber.

This paper studies the binary disruption problem and asteroid capture mechanism in a sun-planet-binary four-body system. Firstly, the binary disruption condition is studied and the result shows that the binary is always disrupted at the perigee of their orbit instantaneously. Secondly, an analytic expression to describe the energy exchange between the binary is derived based on the instantaneous disruption hypothesis. The analytic result is validated through numerical integration. We obtain the energy exchange in encounters simultaneously by the analytic expression and numerical integration. The maximum deviation of this two results is always less than 25% and the mean deviation is about 8.69%. The analytic expression can give us an intuitive description of the energy exchange between the binary. It indicates that the energy change depends on the hyperbolic shape of the binary orbit with respect to the planet, the masses of planet and the primary member of the binary, the binary phase at perigee. We can illustrate the capture/escape processes and give the capture/escape region of the binary clearly by numerical simulation. We analyse the influence of some critical factors to the capture region finally.

The stability and topological structure of equilibrium points in the potential field of the asteroid 101955 Bennu have been investigated with a variable density and rotation period. A dimensionless quantity is introduced for the nondimensionalization of the equations of motion, and this quantity can indicate the effect of both the rotation period and bulk density of the asteroid. Using the polyhedral model of the asteroid Bennu, the number and position of the equilibrium points are calculated and illustrated by a contour plot of the gravitational effective potential field. The topological structure and the stability of the equilibrium points are also investigated using the linearized method. The results show that there are nine equilibrium points in the potential field of the asteroid Bennu, eight in the exterior of the body and one in the interior of the body. Moreover, bifurcation will occur with a variation of the density and rotation period. Different equilibrium points will encounter each other and mix together. Thus, the number of equilibrium points will change. The stability and topological structure of the equilibrium points will also change because of the variation of the density and rotation period of the asteroid. When considering the error of the density of Bennu, the range of the dimensionless quantity covers the critical values that will lead to bifurcation. This means that the stability of the equilibrium points is uncertain, making the dynamical environment of Bennu much more complicated. These bifurcations can help better understand the dynamic environment of an irregular-shaped asteroid.

A method for classifying orbits near asteroids under a polyhedral gravitational field is presented, and may serve as a valuable reference for spacecraft orbit design for asteroid exploration. The orbital dynamics near asteroids are very complex. According to the variation in orbit characteristics after being affected by gravitational perturbation during the periapsis passage, orbits near an asteroid can be classified into 9 categories: (1) surroundingto-surrounding, (2) surrounding-to-surface, (3) surroundingto-infinity, (4) infinity-to-infinity, (5) infinity-to-surface, (6) infinity-to-surrounding, (7) surface-to-surface, (8) surfaceto-surrounding, and (9) surface-to- infinity. Assume that the orbital elements are constant near the periapsis, the gravitation potential is expanded into a harmonic series. Then, the influence of the gravitational perturbation on the orbit is studied analytically. The styles of orbits are dependent on the argument of periapsis, the periapsis radius, and the periapsis velocity. Given the argument of periapsis, the orbital energy before and after perturbation can be derived according to the periapsis radius and the periapsis velocity. Simulations have been performed for orbits in the gravitational field of 216 Kleopatra. The numerical results are well consistent with analytic predictions.

The equilibrium points of the gravitational potential field of minor celestial bodies, including asteroids, comets, and irregular satellites of planets, are studied. In order to understand better the orbital dynamics of massless particles moving near celestial minor bodies and their internal structure, both internal and external equilibrium points of the potential field of the body are analyzed. In this paper, the location and stability of the equilibrium points of 23 minor celestial bodies are presented. In addition, the contour plots of the gravitational effective potential of these minor bodies are used to point out the differences between them. Furthermore, stability and topological classifications of equilibrium points are discussed, which clearly illustrate the topological structure near the equilibrium points and help to have an insight into the orbital dynamics around the irregular-shaped minor celestial bodies. The results obtained here show that there is at least one equilibrium point in the potential field of a minor celestial body, and the number of equilibrium points could be one, five, seven, and nine, which are all odd integers. It is found that for some irregular-shaped celestial bodies, there are more than four equilibrium points outside the bodies while for some others there are no external equilibrium points. If a celestial body has one equilibrium point inside the body, this one is more likely linearly stable.

Near Earth Asteroids have a possibility of impacting with the Earth and always have a thread on the Earth. This paper proposes a way of changing the trajectory of the asteroid to avoid the impaction. Solar sail evolving in a H-reversal trajectory is utilized for asteroid deflection. Firstly, the dynamics of solar sail and the characteristics of the H-reversal trajectory are analyzed. Then, the attitude of the solar sail is optimized to guide the sail to impact with the object asteroid along a H-reversal trajectory. The impact velocity depends on two important parameters: the minimum solar distance along the trajectory and lightness number. A larger lightness number and a smaller solar distance lead to a higher impact velocity. Finally, the deflection capability of a solar sail impacting with the asteroid along the H-reversal is discussed. The results show that a 10 kg solar sail with a lead-time of one year can move Apophis out of a 600-m keyhole area in 2029 to eliminate the possibility of its resonant return in 2036.

Advanced solar sailing has been an increasingly attractive propulsion system for highly non-Keplerian orbits. Three new applications of the orbital angular momentum reversal (H-reversal) trajectories using solar sails are presented in this paper: space observation, heliocentric orbit transfer, and collision orbits with asteroids. A theoretical proof for the existence of double H-reversal trajectories (referred to as 'H2RTs') is given, and the characteristics of the H2RTs are introduced before the discussion of the mission applications. A new family of H2RTs was obtained using a 3D dynamic model of the two-body frame. In a time-optimal control model, the minimum period H2RTs both inside and outside the ecliptic plane were examined using an ideal solar sail. Due to the quasi-heliostationary property at its two symmetrical aphelia, the H2RTs were deemed suitable for space observation. For the second application, the heliocentric transfer orbit was able to function as the time-optimal H-reversal trajectory, as its perihelion velocity is circular or elliptic velocity. Such a transfer orbit can place the sailcraft into a clockwise orbit in the ecliptic plane, with a high inclination or displacement above or below the Sun. The third application of the H-reversal trajectory was simulated impacting an asteroid passing near Earth in a head-on collision. The collision point can be designed through selecting different perihelia or different launch windows. Sample orbits of each application were presented through numerical simulation. The results can serve as a reference for theoretical research and engineering design.