Source author record

Li-Yong Zhou

Li-Yong Zhou appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2020arXiv

Asteroid migration due to the Yarkovsky effect and the distribution of the Eos family

Based on a linearized model of the Yarkovsky effect, we investigate in this paper the dependence of the semimajor axis drift $Δa$ of a celestial body on its size, spinning obliquity, initial orbit and thermal parameters on its surface. With appropriate simplification and approximation, we obtain the analytical solutions to the perturbation equations for the motion of asteroids influenced by the Yarkovsky effect, and they are then verified by numerical simulations of the full equations of motion. These solutions present explicitly the dependencies of $Δa$ on the thermal and dynamical parameters of the asteroid. With these analytical formulae for $Δa$, we investigate the combined seasonal and diurnal Yarkovsky effects. The critical points where the migration direction reverses are calculated and the consequent selective effects according to the size and rotation state of asteroids are discussed. %Solely the Yarkovsky effect is found to be able to produce some ring structure in the aged circumstellar debris disk. Finally, we apply the analytical formulae to calculate the migration of Eos family members. The space distribution of asteroids is well reproduced. Our calculations suggest that statistically the orientations of spin axes of family members satisfy a random-obliquity distribution, and the rotation rate $ω_{\rm rot}$ of asteroid depends on its size $R$ by $ω_{\rm rot}\propto R^{-1}$.

preprint2020arXiv

Dynamics of charged dust in the orbit of Venus

We study the dynamics of co-orbital dust in the inner solar system, i.e. the role of the solar radiation pressure, Poynting-Robertson effect, solar wind, and the interplanetary magnetic field on the location, width and stability of resonant motion of charged, and micron sized dust grains situated in the 1:1 mean motion resonance with planet Venus. We find deviations and asymmetry between $L_4$ and $L_5$ in locations of libration centers and libration width under the influence of non-gravitational effects via both analytical and numerical methods. The triangular Lagrangian points become unstable once we take into consideration solar radiation pressure, the Poynting-Robertson effect and solar wind drag. The Lorentz force could further destabilize the orbits, especially for small dust particles. We also make a comparison between the circular, elliptic restricted three-body model and a more complete model including all planets.

preprint2020arXiv

Transit of asteroids across the 7/3 Kirkwood gap under the Yarkovsky effect

Many asteroids in the main belt are continuously pushed by Yarkovsky effect into regions of different mean motion resonances (MMRs) and then ejected out. They are considered as the principal source of near-Earth objects. We investigate in this paper the effects of the 7/3 MMR with Jupiter (J7/3 MMR) on the transportation of asteroids from Koronis and Eos families that reside respectively on the inner and outer side of the resonance. The fraction of asteroids that make successful crossing through the resonance and the escaping rate from the resonance are found to depend on the Yarkovsky drifting rate, the initial inclination and the migrating direction. The excitation of eccentricity and inclination due to the combined influence from both the resonance and Yarkovsky effect is discussed. Only the eccentricity can be pumped up considerably, and it is attributed mainly to the resonance. In the observational data, family members are also found in the resonance and on the opposite side of the resonance with respect to the corresponding family centre. The existence of these family members is explained using our results of numerical simulations. Finally, the replenishment of asteroids in the J7/3 MMR and the transportation of asteroids by it are discussed.

preprint2019arXiv

A study of the high-inclination population in the Kuiper belt -- III. The 4:7 mean motion resonance

The high-inclination population in the 4:7 mean motion resonance (MMR) with Neptune has also substantial eccentricities ($e\gtrsim0.1$), with more inclined objects tending to occupy more eccentric orbits. For this high-order resonance, there are two different resonant modes. The principal one is the eccentricity-type mode, and we find that libration is permissible for orbits with $e\ge e_c^0$, where the critical eccentricity $e_c^0$ increases as a function of increasing inclination $i$. Correspondingly, we introduce a limiting curve $e_c^0(i)$, which puts constraints on the $(e, i)$ distribution of possible 4:7 resonators. We then perform numerical simulations on the sweep-up capture and long-term stability of the 4:7 MMR, and the results show that the simulated resonators are well-constrained by this theoretical limiting curve. The other 4:7 resonant mode is the mixed-$(e, i)$-type, and we show that stable resonators should exist at $i\gtrsim20^{\circ}$. We predict that the intrinsic number of these mixed-$(e, i)$-type resonators may provide a new clue into the Solar system's evolution, but, so far, only one real object has been observed resonating in this mode.

preprint2019arXiv

A systematic survey of the dynamics of Uranus Trojans

We aim to locate the stability region for Uranus Trojans (UT hereafter) and find out the dynamical mechanisms responsible for the structures in the phase space. Using the spectral number as the stability indicator, we construct the dynamical maps on the (a0, i0) plane. The proper frequencies of UTs are determined precisely so that we can depict the resonance web via a semi-analytical method. Two main stability regions are found, one each for the low-inclination (0-14deg) and high-inclination regime (32-59deg). There is also an instability strip in each of them, at 9deg and 51deg respectively. All stability regions are in the tadpole regime and no stable horseshoe orbits exist for UTs. The lack of moderate-inclined UTs is caused by the nu5 and nu7 secular resonances. The fine structures in the dynamical maps are shaped by high-degree secular resonances and secondary resonances. During the planetary migration, about 36.3% and 0.4% of the pre-formed orbits survive the fast and slow migrations (with migrating time scales of 1 and 10Myr) respectively, most of which are in high inclination. Since the low-inclined UTs are more likely to survive the age of the solar system, they make up 77% of all such long-life orbits by the end of the migration, making a total fraction up to 4.06E-3 and 9.07E-5 of the original population for the fast and slow migrations, respectively. About 3.81% UTs are able to survive the age of the solar system, among which 95.5% are on low-inclined orbits with i0<7.5deg. However, the depletion of the planetary migration seems to prevent a large fraction of such orbits, especially for the slow migration model.

preprint2014arXiv

A study of the high-inclination population in the Kuiper belt - II. The Twotinos

As the second part of our study, in this paper we proceed to explore the dynamics of the high-inclination Twotinos in the 1:2 Neptune mean motion resonance (NMMR). Depending on the inclination $i$, we show the existence of two critical eccentricities $e_a(i)$ and $e_c(i)$, which are lower limits of the eccentricity $e$ for the resonant angle $σ$ to exhibit libration and asymmetric libration, respectively. Accordingly, we have determined the libration centres $σ_0$ for inclined orbits, which are strongly dependent on $i$. With initial $σ=σ_0$ on a fine grid of $(e, i)$, the stability of orbits in the 1:2 NMMR is probed by 4-Gyr integrations. It is shown that symmetric librators are totally unstable for $i\ge30^{\circ}$; while stable asymmetric librators exist for $i$ up to $90^{\circ}$. We further investigate the 1:2 NMMR capture and retention of planetesimals with initial inclinations $i_0\le90^{\circ}$ in the planet migration model using a time-scale of $2\times10^7$ yr. We find that: (1) the capture efficiency of the 1:2 NMMR decreases drastically with the increase of $i_0$, and it goes to 0 when $i_0\ge60^{\circ}$; (2) the probability of discovering Twotinos with $i>25^{\circ}$, beyond observed values, is roughly estimated to be $\le0.1$ per cent; (3) more particles are captured into the leading rather than the trailing asymmetric resonance for $i_0\le10^{\circ}$, but this number difference appears to be the opposite at $i_0=20^{\circ}$ and is continuously varying for even larger $i_0$; (4) captured Twotinos residing in the trailing resonance or having $i>15^{\circ}$ are practically outside the Kozai mechanism, like currently observed samples.

preprint2013arXiv

A study of the high-inclination population in the Kuiper belt -- 1. The Plutinos

The dynamics of the high-inclination Plutinos is systematically studied. We first present the peculiar features of the 2:3 Neptune mean motion resonance (NMMR) for inclined orbits, especially for the correlation of resonant amplitude A_σ with inclination i. Using the numerical integrations for the age of the Solar system, the dynamical structure of the 2:3 NMMR is mapped out on the plane of semi-major axis versus i for different eccentricities. We have shown that i of stable resonant orbits could be as high as 90 deg; and the stable region is roughly surrounded by the contours of A_σ = 120 deg. These new findings allow us to further explore the 2:3 NMMR capture and retention of planetesimals with initial inclinations i0 =< 90 deg in the frame of the planet migration model. We find that the outward transportation of Plutinos is possible for any inclined or even perpendicular orbits. The role of i0 in the formation of Plutinos during Neptune's migration is highlighted and interesting results are obtained: (1) The capture efficiency of the 2:3 NMMR decreases drastically first with the increase of i0, but it then raises instead when i0 exceeds ~ 50 deg; (2) The magnitude of i-variation is limited to less than 5 deg for any i0, and moreover, for Plutinos with i > 48 deg, their i are forced to decrease throughout the outward migration; (3) Plutinos with i > 48 deg are certainly outside the Kozai mechanism, since an inclination increase is prohibited by the migrating 2:3 NMMR; (4) The 7:11 inclination-type NMMR could be responsible for nearly-circular Plutinos, and a minimum i0 ~ 15 deg is required to intrigue this mechanism.

preprint2013arXiv

Modeling the Migration of Neptune and the Corresponding Resonant Captures

Due to the angular momentum exchange with planetesimals, Neptune might have migrated outward to the current position, and captured many Kuiper belt objects (KBOs) into resonances. We set up a semi-analytic model to simulate the outward migration of Neptune, and the processes of resonant captures. Our model can naturally explain Neptune's currently observed semi-major axis and eccentricity. The results show that the current population ratio between 3:2 and 2:1 is mainly due to the original density distribution of KBOs, which might be related to drag-induced inward migrations of proto-KBOs.

preprint2012arXiv

Hyperbolic Structure and Stickiness Effect: A case of a 2D Area-Preserving Twist Mapping

The stickiness effect suffered by chaotic orbits diffusing in the phase space of a dynamical system is studied in this paper. Previous works have shown that the hyperbolic structures in the phase space play an essential role in causing the stickiness effect. We present in this paper the relationship between the stickiness effect and the geometric property of hyperbolic structures. Using a two-dimensional area-preserving twist mapping as the model, we develop the numerical algorithms for computing the positions of the hyperbolic periodic orbits and for calculating the angle between the stable and unstable manifolds of the hyperbolic periodic orbit. We show how the stickiness effect and the orbital diffusion speed are related to the angle.

preprint2012arXiv

Trailing (L5) Neptune Trojans: 2004 KV18 and 2008 LC18

The population of Neptune Trojans is believed to be bigger than that of Jupiter Trojans and that of asteroids in the main belt, although only eight members of this far distant asteroid swarm have been observed up to now. Six leading Neptune Trojans around the Lagrange point L4 discovered earlier have been studied in detail, but two trailing ones found recently around the L5 point, 2004 KV18 and 2008 LC18, have not been investigated yet. In this paper, we report our investigations on the dynamical behaviors of these two new Neptune Trojans. Our calculations show that the asteroid 2004 KV18 is a temporary Neptune Trojan. Most probably, it was captured into the trailing Trojan cloud no earlier than 203kyr ago, and it will not keep this identity no later than 165kyr in future. Based on the statistics on our orbital simulations, we argue that this object is more like a scattered Kuiper belt object. On the contrary, the orbit of asteroid 2008 LC18 is much more stable. Among the clone orbits spread within the orbital uncertainties, a considerable portion of clones may survive on the L5 tadpole orbits for 4Gyr. The strong dependence of the stability on the semimajor axis and resonant angle suggests that further observations are badly needed to confine the orbit in the stable region. We also discuss the implications of the existence and dynamics of these two trailing Trojans on the Solar system history.

preprint2010arXiv

Dynamics of Neptune's Trojans: II. Eccentric orbits and observed ones

In a previous paper, we have presented a global view of the stability of Neptune Trojan (NT hereafter) on inclined orbit. We discuss in this paper the dependence of stability of NT orbits on the eccentricity. High-resolution dynamical maps are constructed using the results of extensive numerical integrations of orbits initialized on the fine grids of initial semimajor axis (a0) versus eccentricity (e0). The extensions of regions of stable orbits on the (a0, e0) plane at different inclinations are shown. The maximum eccentricities of stable orbits in three most stable regions at low (0, 12deg.), medium (22,36deg.) and high (51, 59deg.) inclination, are found to be 0.10, 0.12 and 0.04, respectively. The fine structures in the dynamical maps are described. Via the frequency analysis method, the mechanisms that portray the dynamical maps are revealed. The secondary resonances, concerning the frequency of the librating resonant angle and the frequency of the quasi 2:1 mean motion resonance between Neptune and Uranus, are found deeply involved in the motion of NTs. Secular resonances are detected and they also contribute significantly to the triggering of chaos in the motion. Particularly, the effects of the secular resonance v8, v18 are clarified. We also investigate the orbital stabilities of six observed NTs by checking the orbits of hundreds clones of them generated within the observing error bars. We conclude that four of them, except 2001 QR322 and 2005 TO74, are deeply inside the stable region. The 2001 QR322 is in the close vicinity of the most significant secondary resonance. The 2005 TO74 locates close to the boundary separating stable orbits from unstable ones, and it may be influenced by a secular resonance.

preprint2010arXiv

Where are the Uranus Trojans?

The area of stable motion for fictitious Trojan asteroids around Uranus' equilateral equilibrium points is investigated with respect to the inclination of the asteroid's orbit to determine the size of the regions and their shape. For this task we used the results of extensive numerical integrations of orbits for a grid of initial conditions around the points L4 and L5, and analyzed the stability of the individual orbits. Our basic dynamical model was the Outer Solar System (Jupiter, Saturn, Uranus and Neptune). We integrated the equations of motion of fictitious Trojans in the vicinity of the stable equilibrium points for selected orbits up to the age of the Solar system of 5 billion years. One experiment has been undertaken for cuts through the Lagrange points for fixed values of the inclinations, while the semimajor axes were varied. The extension of the stable region with respect to the initial semimajor axis lies between 19.05 < a < 19.3 AU but depends on the initial inclination. In another run the inclination of the asteroids' orbit was varied in the range 0 < i < 60 and the semimajor axes were fixed. It turned out that only four 'windows' of stable orbits survive: these are the orbits for the initial inclinations 0 < i < 7, 9 < i < 13, 31 < i < 36 and 38 < i < 50. We postulate the existence of at least some Trojans around the Uranus Lagrange points for the stability window at small and also high inclinations.

preprint2009arXiv

The Dynamics of Neptune Trojan: I. the Inclined Orbits

The stability of Trojan type orbits around Neptune is studied. As the first part of our investigation, we present in this paper a global view of the stability of Trojans on inclined orbits. Using the frequency analysis method based on the FFT technique, we construct high resolution dynamical maps on the plane of initial semimajor axis $a_0$ versus inclination $i_0$. These maps show three most stable regions, with $i_0$ in the range of $(0^\circ,12^\circ), (22^\circ,36^\circ)$ and $(51^\circ,59^\circ)$ respectively, where the Trojans are most probably expected to be found. The similarity between the maps for the leading and trailing triangular Lagrange points $L_4$ and $L_5$ confirms the dynamical symmetry between these two points. By computing the power spectrum and the proper frequencies of the Trojan motion, we figure out the mechanisms that trigger chaos in the motion. The Kozai resonance found at high inclination varies the eccentricity and inclination of orbits, while the $ν_8$ secular resonance around $i_0\sim44^\circ$ pumps up the eccentricity. Both mechanisms lead to eccentric orbits and encounters with Uranus that introduce strong perturbation and drive the objects away from the Trojan like orbits. This explains the clearance of Trojan at high inclination ($>60^\circ$) and an unstable gap around $44^\circ$ on the dynamical map. An empirical theory is derived from the numerical results, with which the main secular resonances are located on the initial plane of $(a_0,i_0)$. The fine structures in the dynamical maps can be explained by these secular resonances.