This study presents a study of equilibrium points, periodic orbits, stabilities, and manifolds in a rotating plane symmetric potential field. It has been found that the dynamical behaviour near equilibrium points is completely determined by the structure of the submanifolds and subspaces. The non-degenerate equilibrium points are classified into twelve cases. The necessary and sufficient conditions for linearly stable, non resonant unstable and resonant equilibrium points are established. Furthermore, the results show that a resonant equilibrium point is a Hopf bifurcation point. In addition, if the rotating speed changes, two non degenerate equilibria may collide and annihilate each other. The theory developed here is lastly applied to two particular cases, motions around a rotating, homogeneous cube and the asteroid 1620 Geographos. We found that the mutual annihilation of equilibrium points occurs as the rotating speed increases, and then the first surface shedding begins near the intersection point of the x axis and the surface. The results can be applied to planetary science, including the birth and evolution of the minor bodies in the Solar system, the rotational breakup and surface mass shedding of asteroids, etc.

The rotational speeds of irregular shaped minor celestial bodies can be changed by the YORP effect. This variation in speed can make the numbers, positions, stabilities, and topological cases of the minor body s relative equilibrium points vary. The numbers of relative equilibrium points can be reduced through the collision and annihilation of relative equilibrium points, or increase through the creation and separation of relative equilibrium points. Here we develop a classification system of multiple annihilation behaviors of the equilibrium points for irregular shaped minor celestial bodies. Most minor bodies have five equilibrium points; there are twice the number of annihilations per equilibrium point when the number of equilibrium points is between one and five. We present the detailed annihilation classifications for equilibria of objects which have five equilibrium points. Additionally, the annihilation classification for the seven equilibria of Kleopatra shaped objects and the nine equilibria of Bennu shaped objects are also discussed. Equilibria of different objects fall into different annihilation classifications. By letting the rotational speeds vary, we studied the annihilations and creations of relative equilibria in the gravitational field of ten minor bodies.

We used binary octahedrons to investigate the dynamical behaviors of binary asteroid systems. The mutual potential of the binary polyhedron method is derived from the fourth order to the sixth order. The irregular shapes, relative orbits, attitude angles, as well as the angular velocities of the binary asteroid system are included in the model. We investigated the relative trajectory of the secondary relative to the primary, the total angular momentum and total energy of the system, the three-axis attitude angular velocity of the binary system, as well as the angular momentum of the two components. The relative errors of total angular momentum and total energy indicate the calculation has a high precision. It is found that the influence of the orbital and attitude motion of the primary from the gravitational force of the secondary is obvious. This study is useful for understanding the complicated dynamical behaviors of the binary asteroid systems discovered in our Solar system.

We investigate the orbital stability close to the unique L4-point Jupiter binary Trojan asteroid 624 Hektor. The gravitational potential of 624 Hektor is calculated using the polyhedron model with observational data of 2038 faces and 1021 vertexes. Previous studies have presented three different density values for 624 Hektor. The equilibrium points in the gravitational potential of 624 Hektor with different density values have been studied in detail. There are five equilibrium points in the gravitational potential of 624 Hektor no matter the density value. The positions, Jacobian, eigenvalues, topological cases, stability, as well as the Hessian matrix of the equilibrium points are investigated. For the three different density values the number, topological cases, and the stability of the equilibrium points with different density values are the same. However, the positions of the equilibrium points vary with the density value of the asteroid 624 Hektor. The outer equilibrium points move away from the asteroid s mass center when the density increases, and the inner equilibrium point moves close to the asteroid s mass center when the density increases. There exist unstable periodic orbits near the surface of 624 Hektor. We calculated an orbit near the primary s equatorial plane of this binary Trojan asteroid, the results indicate that the orbit remains stable after 28.8375 d.

We present the analysis and computational results for the inclination relative effect of moonlets of triple asteroidal systems. Perturbations on moonlets due to the primarys non-sphericity gravity, the solar gravity, and moonlets relative gravity are discussed. The inclination vector for each moonlet follows a periodic elliptical motion; the motion period depends on the moonlets semi-major axis and the primarys J2 perturbations. Perturbation on moonlets from the Solar gravity and moonlets relative gravity makes the motion of the x component of the inclination vector of moonlet 1 and the y component of the inclination vector of moonlet 2 to be periodic.The mean motion of x component and the y component of the inclination vector of each moonlet forms an ellipse. However, the instantaneous motion of x component and the y component of the inclination vector may be an elliptical disc due to the coupling effect of perturbation forces. Furthermore, the x component of the inclination vector of moonlet 1 and the y component of the inclination vector of moonlet 2 form a quasi-periodic motion. Numerical calculation of dynamical configurations of two triple asteroidal systems (216) Kleopatra and (153591) 2001 SN263 validates the conclusion.

In this paper we analyze the dynamical behavior of large dust grains in the vicinity of a cometary nucleus. To this end we consider the gravitational field of the irregularly shaped body, as well as its electric and magnetic fields. Without considering the effect of gas friction and solar radiation, we find that there exist grains which are static relative to the cometary nucleus; the positions of these grains are the stable equibria. There also exist grains in the stable periodic orbits close to the cometary nucleus. The grains in the stable equibria or the stable periodic orbits will not escape or impact on the surface of the cometary nucleus. The results are applicable for large charge dusts with small area-mass ratio which are near the cometary nucleus and far from the Solar. It is found that the resonant periodic orbit can be stable, and there exist stable non-resonant periodic orbits, stable resonant periodic orbits and unstable resonant periodic orbits in the potential field of cometary nuclei. The comet gravity force, solar gravity force, electric froce, magnetic force, solar radiation pressure, as well as the gas drag force are all considered to analyze the order of magnitude of these forces acting on the grains with different parameters. Let the distance of the dust grain relative to the mass centre of the cometary nucleus, the charge and the mass of the dust grain vary, respectively, fix other parameters, we clacluated the strengths of different forces.

The discovery of binary and triple asteroids in addition to the execution of space missions to minor celestial bodies in the past several years have focused increasing attention on periodic orbits around irregular-shaped celestial bodies. In the present work, we adopt a polyhedron shape model for providing an accurate representation of irregular-shaped bodies, and employ the model to calculate their corresponding gravitational and effective potentials. We also investigate the characteristics of periodic orbit families and the continuation of periodic orbits. We prove a fact, which provides a conserved quantity that permits restricting the number of periodic orbits in a fixed energy curved surface about an irregular-shaped body. The collisions of Floquet multipliers are maintained during the continuation of periodic orbits around the comet 1P/Halley. Multiple bifurcations in the periodic orbit families about irregular-shaped bodies are also discussed. Three bifurcations in the periodic orbit family have been found around the asteroid 216 Kleopatra, which include two real saddle bifurcations and one period-doubling bifurcation.

This manuscript considers the main features of the nonlinear dynamics of multiple irregular celestial body systems. The gravitational potential, static electric potential, and magnetic potential are considered. Based on the three established potentials, we show that three conservative values exist for this system, including a Jacobi integral. The equilibrium conditions for the system are derived and their stability analyzed. The equilibrium conditions of a celestial system comprised of n irregular bodies are reduced to 12n minus 9 equations. The dynamical results are applied to simulate the motion of multiple-asteroid systems. The simulation is useful for the study of the stability of multiple irregular celestial body systems and for the design of spacecraft orbits to triple asteroid systems discovered in the solar system. The dynamical configurations of the five triple-asteroid systems 45 Eugenia, 87 Sylvia, 93 Minerva, 216 Kleopatra, and 136617 1994CC, and the six-body system 134340 Pluto are calculated and analyzed.

We study the motion and equilibria of the grains on the surface of the irregular celestial body (hereafter called irregular bodies). Motions for the grains on the smooth and unsmooth surfaces are discussed, respectively. The linearized equations of motion relative to a surface equilibrium point and its characteristic equations are presented. Considering the stick-slip effect, the damping forces and the spring forces for the grain are calculated, then the linearized equations of motion and the characteristic equations relative to the surface equilibrium points are derived. The number of non-degenerate surface equilibria is an even number. We compute the motion of a grain released above three different regions relative to the irregular asteroid 6489 Golevka, including the flat surface, the concave region, and the convex region. Following the grain release and initial bounce, three kinds of motions exist, the orbital motion, the impact motion and the surface motion. We find that the maximum height of the next hop may be bigger than the maximum height of the current hopping. We also used Monte Carlo simulations to calculate 100 grains hopping motions, the results shows that the stable surface equilibria are on the concave region and flat surface of the asteroid.

We investigate the multiple bifurcations in periodic orbit families in the potential field of a highly irregular-shaped celestial body. Topological cases of periodic orbits and four kinds of basic bifurcations in periodic orbit families are studied. Multiple bifurcations in periodic orbit families consist of four kinds of basic bifurcations. We found both binary period-doubling bifurcations and binary tangent bifurcations in periodic orbit families around asteroid 433 Eros. The periodic orbit family with binary period-doubling bifurcations is nearly circular, with almost zero inclination, and is reversed relative to the body of the asteroid 433 Eros. This implies that there are two stable regions separated by one unstable region for the motion around this asteroid. In addition, we found triple bifurcations which consist of two real saddle bifurcations and one period-doubling bifurcation. A periodic orbit family generated from an equilibrium point of asteroid 433 Eros has five bifurcations, which are one real saddle bifurcation, two tangent bifurcations, and two period-doubling bifurcations.

Some asteroids contain water ice, and a space mission landing on an asteroid may take liquid to the surface of the asteroid. Gas pressure is very weak on the surface of asteroids. Here we consider the capillary action in a crack on the surface of irregular asteroids. The crack is modelled as a capillary which has a fixed radius. An asteroid s irregular gravitational potential influences the height of the liquid in the capillary. The height of the liquid in the capillary on the surface of such asteroids is derived from the asteroid s irregular gravitational potential. Capillary mechanisms are expected to produce an inhomogeneaous distribution of emergent liquid on the surface. This result is applied to asteroid 433 Eros, which has an irregular, elongated, and concave shape. Two cases are considered 1) we calculate the height of the liquid in the capillary when the direction of the capillary is perpendicular to the local surface of the asteroid; 2) we calculate the height of the liquid in the capillary when the direction of the capillary is parallel to the vector from the center of mass to the surface position. The projected height in the capillary on the local surface of the asteroid seems to depend on the assumed direction of the capillary.

We present a design concept for a space engine that can continuously remove the orbit debris by using the debris as a propellant. Space robotic cleaner is adopted to capture the targeting debris and to transfer them into the engine. Debris with larger size is first disintegrated into small pieces by using a mechanical method. The planetary ball mill is then adopted to grind the pieces into micrometer or smaller powder. The energy needed in this process is get from the nuclear and solar power. By the effect of gamma-ray photoelectric or the behavior of tangently rub of tungsten needles, the debris powered is charged. This behavior can be used to speed up the movement of powder in a tandem electrostatic particle accelerator. By ejecting the high-temperture and high-pressure charged powered from the nozzle of the engine,the continuously thrust is obtained. This thrust can be used to perform orbital maneuver and debris rendezvous for the spacecraft and robotic cleaner. The ejected charged particle will be blown away from the circumterrestrial orbit by the solar wind. By digesting the space debris, we obtain not only the previous thrust but also the clean space. In the near future, start trek will not just a dream, human exploration will extend to deep universe. The analysis shown, the magnitude of the specific impulse for debris engine is determined by the accelerating electrostatic potential and the charge-to-mass ratio of the powder.

We are interested in the periodic motion and bifurcations near the surface of an asteroid. The gravity field of an irregular asteroid and the equation of motion of a particle near the surface of an asteroid are studied. The periodic motions around the major body of triple asteroid 216 Kleopatra and the OSIRIS REx mission target asteroid 101955 Bennu are discussed. We find that motion near the surface of an irregular asteroid is quite different from the motion near the surface of a homoplastically spheroidal celestial body. The periodic motions around the asteroid 101955 Bennu and 216 Kleopatra indicate that the geometrical shapes of the orbits are probably very sophisticated. There exist both stable periodic motions and unstable periodic motions near the surface of the same irregular asteroid. This periodic motion which is unstable can be resonant or non resonant. The period doubling bifurcation and pseudo period doubling bifurcation of periodic orbits coexist in the same gravity field of the primary of the triple asteroid 216 Kleopatra. It is found that both of the period doubling bifurcations of periodic orbits and pseudo period-doubling bifurcation of periodic orbits have four different paths. The pseudo period doubling bifurcation found in the potential field of primary of triple asteroid 216 Kleopatra shows that there exist stable periodic orbits near the primary s equatorial plane, which gives an explanation for the motion stability of the triple asteroid 216 Kleopatra s two moonlets, Alexhelios and Cleoselene.

The order and chaos of the motion near equilibrium points in the potential of a rotating highly irregular-shaped celestial body are investigated from point of view of the dynamical system theory. The positions of the non-degenerate equilibrium points vary continuously when the parameter changes. The topological structures in the vicinity of equilibrium points are classified into several different cases. Bifurcations at equilibrium points and the topological transfers between different cases for equilibrium points are also discussed. The conclusions can be applied to all kinds of rotating celestial bodies, simple-shaped or highly irregular-shaped, including asteroids, comets, planets and satellites of planets to help one to understand the dynamical behaviors around them. Applications to asteroids 216 Kleopatra, 2063 Bacchus, and 25143 Itokawa are significant and interesting: eigenvalues affiliated to the equilibrium points for the asteroid 216 Kleopatra move and always belong to the same topological cases; while eigenvalues affiliated to two different equilibrium points for the asteroid 2063 Bacchus and 25143 Itokawa move through the resonant cases of equilibrium points, and the collision of eigenvalues in the complex plane occurs. Poincar\'e sections in the potential of the asteroid 216 Kleopatra show the chaos behaviors of the orbits in large scale.

In this work, we investigate the bifurcations of relative equilibria in the gravitational potential of asteroids. A theorem concerning a conserved quantity, which is about the eigenvalues and number of relative equilibria, is presented and proved. The conserved quantity can restrict the number of non-degenerate equilibria in the gravitational potential of an asteroid. It is concluded that the number of non-degenerate equilibria in the gravitational field of an asteroid varies in pairs and is an odd number. In addition, the conserved quantity can also restrict the kinds of bifurcations of relative equilibria in the gravitational potential of an asteroid when the parameter varies. Furthermore, studies have shown that there exist transcritical bifurcations, quasi-transcritical bifurcations, saddle-node bifurcations, saddle-saddle bifurcations, binary saddle-node bifurcations, supercritical pitchfork bifurcations, and subcritical pitchfork bifurcations for the relative equilibria in the gravitational potential of asteroids. It is found that for the asteroid 216 Kleopatra, when the rotation period varies as a parameter, the number of relative equilibria changes from 7 to 5 to 3 to 1, and the bifurcations for the relative equilibria are saddle-node bifurcations and saddle-saddle bifurcations.

This paper studies the distribution of characteristic multipliers, the structure of submanifolds, the phase diagram, bifurcations and chaotic motions in the potential field of rotating highly irregular-shaped celestial bodies (hereafter called irregular bodies). The topological structure of the submanifolds for the orbits in the potential field of an irregular body is shown to be classified into 34 different cases, including 6 ordinary cases, 3 collisional cases, 3 degenerate real saddle cases, 7 periodic cases, 7 period-doubling cases, 1 periodic and collisional case, 1 periodic and degenerate real saddle case, 1 period-doubling and collisional case, 1 period-doubling and degenerate real saddle case, and 4 periodic and period-doubling cases. The different distribution of the characteristic multipliers has been shown to fix the structure of the submanifolds, the type of orbits, the dynamical behaviour and the phase diagram of the motion. Classifications and properties for each case are presented. Moreover, tangent bifurcations, period-doubling bifurcations, Neimark-Sacker bifurcations and the real saddle bifurcations of periodic orbits in the potential field of an irregular body are discovered. Submanifolds appear to be Mobius strips and Klein bottles when the period-doubling bifurcation occurs.

In this study, the potential locations of asteroidal small satellites (also called moonlets) with quasi-circular mutual orbit are analyzed. For the motion of the moonlets, only the solar gravity perturbation and the primary's 2nd degree-and-order gravity field are considered. By eliminating of short periodic terms, the dynamical behavior of the Hamiltonian for the moonlets is investigated. The observational data of some high size ratio binary asteroids show that the orbits of the moonlets lie close to the classical Laplace equilibria, which reach global minimum values of the Hamiltonian. It is found that tides or Yarkovsky effects alone cannot account for the reason why the orbits of asteroidal moonlets are not exactly at the classical Laplace equilibria. The analysis in this study is expected to provide useful information for the potential locations of asteroidal moonlets, and contribute to principles to relate predictions to observations.

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 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 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 four, and we find at least 1 family of periodic orbits near the resonant equilibrium points. Besides, this theory is applied to asteroids 216 Kleopatra, 1620 Geographos, 4769 Castalia, and 6489 Golevka.

This study investigates the different novel forms of the dynamical equations of a particle orbiting a rotating asteroid and the effective potential, the Jacobi integral, etc. on different manifolds. Nine new forms of the dynamical equations of a particle orbiting a rotating asteroid are presented, and the classical form of the dynamical equations has also been found. The dynamical equations with the potential and the effective potential in scalar form in the arbitrary body-fixed frame and the special body-fixed frame are presented and discussed. Moreover, the simplified forms of the effective potential and the Jacobi integral have been derived. The dynamical equation in coefficient-matrix form has been derived. Other forms of the dynamical equations near the asteroid are presented and discussed, including the Lagrange form, the Hamilton form, the symplectic form, the Poisson form, the Poisson-bracket form, the cohomology form, and the dynamical equations on the K\"ahler manifold and another complex manifold. Novel forms of the effective potential and the Jacobi integral are also presented. The dynamical equations in scalar form and coefficient-matrix form can aid in the study of the dynamical system, the bifurcation, and the chaotic motion of the orbital dynamics of a particle near a rotating asteroid. The dynamical equations of a particle near a rotating asteroid are presented on several manifolds, including the symplectic manifold, the Poisson manifold, and complex manifolds, which may lead to novel methods of studying the motion of a particle in the potential field of a rotating asteroid.

The abstract is additional with repect to the paper published in JGCD. Ordinary Earth satellites are usually placed into five categories of special orbits: sun-synchronous orbits, orbits at the critical inclination, frozen orbits, repeating ground track orbits, and geostationary orbits. This paper investigates their counterparts around Mars and examines the basic nature of these orbits, which are of special interest for missions conducted around Mars, including Mars reconnaissance. Mars' gravity field is much more complicated, with relatively smaller J2, compared to Earth's, which makes the behaviors of these Martian orbits different from those of Earth. Analytical formulations and numerical simulations are used to analyze these Martian orbits and compare them with their Earth counterparts. First, mean element theory is employed to describe variations of orbital elements and give the constraint conditions for achieving these special orbits. Then, numerical verifications based on the PSODE algorithm (particle swarm optimization combined with differential evolution) are adopted to provide more accurate conditions for achieving these orbits when considering an Mars gravity field. Using the numerical method can significantly improve the design in the full gravity field, and it is therefore possible to select these usable orbits for Mars that can reduce or eliminate the need for stationkeeping.

The triple asteroids and triple Kuiper belt objects (collectively called the triple minor planets) in the Solar system are of particular interest to the scientific community since the discovery of the first triple asteroid system in 2004. In this paper, the Hill stability of the nine known triple minor planets in the Solar system is investigated. Seven of the systems are of large size ratio, i.e. they consist of a larger primary and two moonlets, while the other two systems have components of comparable size. Each case is treated separately. For the triple minor planets that have large size ratio, the sufficient condition for Hill stability is expressed in closed form. This is not possible for the systems with comparable size components, for which the Hill stability is assessed by a combination of analytical and numerical means. It is found that all the known triple minor planets are Hill stable, except 3749 Balam, for which the incomplete orbital parameters make the Hill stability of the system uncertain. This suggests that there might be more such stable triple minor planets in the Solar system yet to be observed. It is also shown that the Hill stability regions increase as the mutual inclination between the inner orbit and outer orbit decreases, the semimajor axis ratio of the inner orbit with respect to the outer orbit decreases, and the mass ratio of the outer satellite with respect to the inner satellite increases. This study therefore provides useful information about dynamical properties of the triple minor planets in the Solar system.

This study investigates the problem of areostationary orbits around Mars in the three-dimensional space. Areostationary orbits are expected to be used to establish a future telecommunication network for the exploration of Mars. However, no artificial satellites have been placed in these orbits thus far. In this paper, the characteristics of the Martian gravity field are presented, and areostationary points and their linear stability are calculated. By taking linearized solutions in the planar case as the initial guesses and utilizing the Levenberg-Marquardt method, families of periodic orbits around areostationary points are shown to exist. Short-period orbits and long-period orbits are found around linearly stable areostationary points, and only short-period orbits are found around unstable areostationary points. Vertical periodic orbits around both linearly stable and unstable areostationary points are also examined. Satellites in these periodic orbits could depart from areostationary points by a few degrees in longitude, which would facilitate observation of the Martian topography. Based on the eigenvalues of the monodromy matrix, the evolution of the stability index of periodic orbits is determined. Finally, heteroclinic orbits connecting the two unstable areostationary points are found, providing the possibility for orbital transfer with minimal energy consumption.

In the current study, a double-averaged analytical model including the action of the perturbing body's inclination is developed to study third-body perturbations. The disturbing function is expanded in the form of Legendre polynomials truncated up to the second-order term, and then is averaged over the periods of the spacecraft and the perturbing body. The efficiency of the double-averaged algorithm is verified with the full elliptic restricted three-body model. Comparisons with the previous study for a lunar satellite perturbed by Earth are presented to measure the effect of the perturbing body's inclination, and illustrate that the lunar obliquity with the value 6.68\degree is important for the mean motion of a lunar satellite. The application to the Mars-Sun system is shown to prove the validity of the double-averaged model. It can be seen that the algorithm is effective to predict the long-term behavior of a high-altitude Martian spacecraft perturbed by Sun. The double-averaged model presented in this paper is also applicable to other celestial systems.

Recently, Near Earth Objects (NEOs) have been attracting great attention, and thousands of NEOs have been found to date. This paper examines the NEOs' orbital dynamics using the framework of an accurate solar system model and a Sun-Earth-NEO three-body system when the NEOs are close to Earth to search for NEOs with low-energy orbits. It is possible for such an NEO to be temporarily captured by Earth; its orbit would thereby be changed and it would become an Earth-orbiting object after a small increase in its velocity. From the point of view of the Sun-Earth-NEO restricted three-body system, it is possible for an NEO whose Jacobian constant is slightly lower than C1 and higher than C3 to be temporarily captured by Earth. When such an NEO approaches Earth, it is possible to change its orbit energy to close up the zero velocity surface of the three-body system at point L1 and make the NEO become a small satellite of the Earth. Some such NEOs were found; the best example only required a 410m/s increase in velocity.

This paper investigates the dynamics of a particle orbiting around a rotating homogeneous cube, and shows fruitful results that have implications for examining the dynamics of orbits around non-spherical celestial bodies. This study can be considered as an extension of previous research work on the dynamics of orbits around simple shaped bodies, including a straight segment, a circular ring, an annulus disk, and simple planar plates with backgrounds in celestial mechanics. In the synodic reference frame, the model of a rotating cube is established, the equilibria are calculated, and their linear stabilities are determined. Periodic orbits around the equilibria are computed using the traditional differential correction method, and their stabilities are determined by the eigenvalues of the monodromy matrix. The existence of homoclinic and heteroclinic orbits connecting periodic orbits around the equilibria is examined and proved numerically in order to understand the global orbit structure of the system. This study contributes to the investigation of irregular shaped celestial bodies that can be divided into a set of cubes.

The critical inclination is of special interest in artificial satellite theory. The critical inclination can maintain minimal deviations of eccentricity and argument of pericentre from the initial values, and orbits at this inclination have been applied to some space missions. Most previous researches about the critical inclination were made under the assumption that the oblateness term J2 is dominant among the harmonic coefficients. This paper investigates the extension of the critical inclination where the concept of the critical inclination is different from that of the traditional sense. First, the study takes the case of Venus for instance, and provides some preliminary results. Then for general cases, given the values of argument of pericentre and eccentricity, the relationship between the multiplicity of the solutions for the critical inclination and the values of J2 and J4 is analyzed. Besides, when given certain values of J2 and J4, the relationship between the multiplicity of the solutions for the critical inclination and the values of semimajor axis and eccentricity is studied. The results show that for some cases, the value of the critical inclination is far away from that of the traditional sense or even has multiple solutions. The analysis in this paper could be used as starters of correction methods in the full gravity field of celestial bodies

In the current study, the existence of periodic orbits around a fixed homogeneous cube is investigated, and the results have powerful implications for examining periodic orbits around non-spherical celestial bodies. In the two different types of symmetry planes of the fixed cube, periodic orbits are obtained using the method of the Poincar\'e surface of section. While in general positions, periodic orbits are found by the homotopy method. The results show that periodic orbits exist extensively in symmetry planes of the fixed cube, and also exist near asymmetry planes that contain the regular Hex cross section. The stability of these periodic orbits is determined on the basis of the eigenvalues of the monodromy matrix. This paper proves that the homotopy method is effective to find periodic orbits in the gravity field of the cube, which provides a new thought of searching for periodic orbits around non-spherical celestial bodies. The investigation of orbits around the cube could be considered as the first step of the complicated cases, and helps to understand the dynamics of orbits around bodies with complicated shapes. The work is an extension of the previous research work about the dynamics of orbits around some simple shaped bodies, including a straight segment, a circular ring, an annulus disk, and simple planar plates.

Frozen orbits are always important foci of orbit design because of their valuable characteristics that their eccentricity and argument of pericentre remain constant on average. This study investigates quasi-circular frozen orbits and examines their basic nature analytically using two different methods. First, an analytical method based on Lagrangian formulations is applied to obtain constraint conditions for Martian frozen orbits. Second, Lie transforms are employed to locate these orbits accurately, and draw the contours of the Hamiltonian to show evolutions of the equilibria. Both methods are verified by numerical integrations in an 80\times80 Mars gravity field. The simulations demonstrate that these two analytical methods can provide accurate enough results. By comparison, the two methods are found well consistent with each other, and both discover four families of Martian frozen orbits: three families with small eccentricities and one family near the critical inclination. The results also show some valuable conclusions: for the majority of Martian frozen orbits, argument of pericentre are kept at 270 degrees because J3 has the same sign with J2; while for a minority of ones with low altitude and low inclination, argument of pericentre are possible to be kept at 90 degrees because of the effect of the higher degree odd zonals; for the critical inclinations cases, argument of pericentre can also be kept at 90 degrees. It is worthwhile to note that there exist some special frozen orbits with extremely small eccentricity, which could provide much convenience for reconnaissance. Finally, the stability of Martian frozen orbits is estimated based on the trace of the monodromy matrix. The analytical investigations can provide good initial conditions for numerical correction methods in the more complex models.

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.