Kuiper带天体的轨道动力学

主要评述太阳系动力学研究的一个新方向——Kuiper带的轨道动力学。早期的研究是为了探讨短周期彗星的起源。在发现第一颗Kuiper带小天体之后，人们开始将注意力转到Kuiper带共振区的相空间结构上，Morbidelli和Malhotra分别采用不同的模型研究了这些共振区的大小。其中主要研究对象是3:2共振区。冥王星也处在这一共振区中。从冥王星的轨道特性来看，冥王星应是一颗较大的Kuiper带天体，它还拥有另外两种共振——Kozai共振和1：1超级共振。正是由于这些共振的存在，冥王星的运动才得以长期保持稳定。观测表明许多Kuiper带天体也处的海王星的平运动共振中。早期的理论认为这些平运动共振起源于灾难性事件，如碰撞。然而这都是一些小概率事件，无法对共振的形成进行合理的解释。Malhotra通过行星迁移成功地解释了冥王星被共振俘获的机制。这一机制的概率非常大，同样可以用来解释Kuiper带天体共振的形成。     

两个可能存在大质量类冥王星的区域

1992年以来，在海王星外的太阳系发现了近千个小天体，称为Kuiper带天体(KBO)或Edgeworth—Kuiper带天体，其中有一部分偏心率和倾角较大的小天体与海王星之间存在3：2平运动共振，轨道特征类似冥王星，命名为类冥王星，自KBO发现以来，天文学家们进行了多次小天区的搜索，发现了几个质量较大的KBO，通过数值计算，在轨道参数空间发现了两个和冥王星一样同时具有3种共振的区域，在这两个区域里的小天体既避免了海王星的强摄动又不会与冥王星密切交会，轨道非常稳定，因此有可能在其中发现质量较大的类冥王星。     

A discussion of the solution for the motion of Pluto

In a previous paper, a semi-analytical solution for the long-term motion of Pluto was presented. The present paper contains: (1) a comparison of the present solution with the solution by Williams and Benson; (2) a discussion of the effect of the near resonance between Pluto and Uranus; and, (3) a calculation of the librational period of the eccentricity, inclination and perihelion.The semi-analytical solution is shown to agree very closely with the long-term solution for Pluto obtained by Williams and Benson using numerical integration of the averaged equations of motion. A small difference between the two solutions is attributed to neglecting the eccentricity and inclination of Neptune in the semi-analytical solution.     

A Symplectic Mapping Model for the Study of 2:3 Resonant Trans-Neptunian Motion

A symplectic mapping is constructed for the study of the dynamical evolution of Edgeworth-Kuiper belt objects near the 2:3 mean motion resonance with Neptune. The mapping is six-dimensional and is a good model for the Poincaré map of the real system, that is, the spatial elliptic restricted three-body problem at the 2:3 resonance, with the Sun and Neptune as primaries. The mapping model is based on the averaged Hamiltonian, corrected by a semianalytic method so that it has the basic topological properties of the phase space of the real system both qualitatively and quantitatively. We start with two dimensional motion and then we extend it to three dimensions. Both chaotic and regular motion is observed, depending on the objects' initial inclination and phase. For zero inclination, objects that are phase-protected from close encounters with Neptune show ordered motion even at eccentricities as large as 0.4 and despite being Neptune-crossers. On the other hand, not-phase-protected objects with eccentricities greater than 0.15 follow chaotic motion that leads to sudden jumps in their eccentricity and are removed from the 2:3 resonance, thus becoming short period comets. As inclination increases, chaotic motion becomes more widespread, but phase-protection still exists and, as a result, stable motion appears for eccentricities up to e = 0.3 and inclinations as high as i = 15°, a region where plutinos exist.     

Long-term behavior of the motion of Pluto over 5.5 billion years

The motion of Pluto is said to be chaotic in the sense that the maximum Lyapunov exponent is positive: the Lyapunov time (the inverse of the Lyapunov exponent) is about 20 million years. So far the longest integration up to now, over 845 million years (42 Lyapunov times), does not show any indication of a gross instability in the motion of Pluto. We carried out the numerical integration of Pluto over the age of the solar system (5.5 billion years ≈ 280 Lyapunov times). This integration also did not give any indication of chaotic evolution of Pluto. The divergences of Keplerian elements of a nearby trajectory at first grow linearly with the time and then start to increase exponentially. The exponential divergences stop at about 420 million years. The divergences in the semi-major axis and the mean anomaly ( equivalently the longitude and the distance) saturate. The divergences of the other four elements, the eccentricity, the inclination, the argument of perihelion, and the longitude of node still grow slowly after the stop of the exponential increase and finally saturate.     

Stability of higher order resonances in the restricted three-body problem

Third and fourth order mean motion resonances are studied in the model of the restricted three-body problem by numerical methods for mass parameters corresponding approximately to the Sun?CJupiter and Sun?CNeptune systems. In the case of inner resonances, it is shown that there are two regions of libration in the 8:5 and 7:4 resonances, one at low, the other at high eccentricities. In the 9:5 and 7:3 resonances libration can exist only in one region at high eccentricities. The 5:2 and 4:1 resonances are very regular, with one librational zone existing for all eccentricities. There is no visible region of libration at any eccentricities in the 5:1 resonance, the transition between the regions of direct and retrograde circulation is very sharp. In the case of outer resonances, the 8:5 and 7:4 resonances have also two regions of libration, but the 9:5 resonance has three, the 7:3 resonance two librational zones. The 5:2 resonance is again very regular, but it is parted for two regions of libration at high eccentricities. Libration is possible in the 4:1 resonance only at high eccentricities. The 5:1 resonance is very symmetric. In the case of outer resonances, a comparison is made with trans-Neptunian objects (TNO) in higher order mean motion resonances. Several new librating TNOs are identified.     

Motions of the perihelions of Neptune and Pluto

Five outer planets are numerically integrated over five million years in the Newtonian frame. The argument of Pluto's perihelion librates about 90 degrees with an amplitude of about 23 degrees. The period of the libration depends on the mass of Pluto: 4.0×10⁶ years forM _pluto=2.78×10^–6 M _sun and 3.8×10⁶ years forM _pluto=7.69×10^–9 M _sun, which is the newly determined mass. The motion of Neptune's perihelion is more sensitive to the mass of Pluto. ForM _pluto=7.69×10^–9 M _sun, the perihelion of Neptune does circulate counter-clockwise and forM _pluto=2.78×10^–6 M _sun, it does not circulate and the Neptune's eccentricity does not have a minimum. With the initial conditions which do not lie in the resonance region between Neptune and Pluto, a close approach between them takes place frequently and the orbit of Pluto becomes unstable and irregular.     

A dust cloud around Pluto and Charon

We suggest that Pluto and Charon are immersed in a tenuous dust cloud. The cloud consists of ejecta from Pluto and—especially—Charon, released from their surfaces by impacts of micrometeoroids originating from Edgeworth-Kuiper belt objects. The motion of the ejected grains is dominated by the gravity of Pluto and Charon, which determines a pear-shape of the densest part of the cloud. While the production rates of escaping particles from both sides are comparable, the lifetimes of the Charon particles inside the Hill sphere of Pluto-Charon with respect to the Sun are much longer than of the Pluto ejecta, so that the cloud is composed predominantly of Charon grains. The dust cloud is dense enough to be detected with an in situ dust detector onboard a future space mission to Pluto. The cloud's maximum optical depth of τ≈3×10⁻¹¹ is, however, too low to allow remote sensing observations.     

Long-term behavior of the motion of Pluto over 5.5 billion years

The motion of Pluto is said to be chaotic in the sense that the maximum Lyapunov exponent is positive: the Lyapunov time (the inverse of the Lyapunov exponent) is about 20 million years. So far the longest integration up to now, over 845 million years (42 Lyapunov times), does not show any indication of a gross instability in the motion of Pluto. We carried out the numerical integration of Pluto over the age of the solar system (5.5 billion years 280 Lyapunov times). This integration also did not give any indication of chaotic evolution of Pluto. The divergences of Keplerian elements of a nearby trajectory at first grow linearly with the time and then start to increase exponentially. The exponential divergences stop at about 420 million years. The divergences in the semi-major axis and the mean anomaly ( equivalently the longitude and the distance) saturate. The divergences of the other four elements, the eccentricity, the inclination, the argument of perihelion, and the longitude of node still grow slowly after the stop of the exponential increase and finally saturate.     

On the long-term motion of Pluto

A modified periodic orbit of the third kind is introduced that is closely related to periodic orbits of the third kind as defined by Poincaré. It is shown that Pluto librates about the periodic orbit with apparent stability. This further explains the librational motion of the resonant argument of Pluto and the avoidance of a Pluto-Neptune close approach as found by Cohen and Hubbard and the long-term motion of Pluto and the librational motion of the perihelion as found by Williams and Benson. With libration about a periodic orbit, the numerical solution of Williams and Benson can be extrapolated to longer times in the past and future.     

The mass ratio of Charon to Pluto from Hubble Space Telescope astrometry with the fine guidance sensors

The mass ratio of Charon to Pluto is a basic parameter describing the binary system and is necessary for determining the individual masses and densities of these two bodies. Previous measurements of the mass ratio have been made, but the solutions differ significantly (Null et al., 1993; Young et al., 1994; Null and Owen, 1996; Foust et al., 1997; Tholen and Buie, 1997). We present the first observations of Pluto and Charon with a well-calibrated astrometric instrument—the fine guidance sensors on the Hubble Space Telescope. We observed the motion of Pluto and Charon about the system barycenter over 4.4 days (69% of an orbital period) and determined the mass ratio to be 0.122±0.008 which implies a density of 1.8 to 2.1 g cm⁻³ for Pluto and 1.6 to 1.8 g cm⁻³ for Charon. The resulting rock-mass fractions for Pluto and Charon are higher than expected for bodies formed in the outer solar nebula, possibly indicating significant postaccretion loss of volatiles.     

Meteorite transport—Revisited II

In Wisdom (2017), I presented new simulations of meteorite transport from the chaotic zones associated with major resonances in the asteroid belt: the ν₆ secular resonance, the 3:1 mean motion resonance with Jupiter, and the 5:2 mean motion resonance with Jupiter. I found that the observed afternoon excess (the fact that approximately twice as many meteorites fall in the afternoon as in the morning) of the ordinary chondrites is consistent with chaotic transport from the 3:1 resonance, contradicting prior reports. Here I report an additional study of the transport of meteorites from ν₆ secular resonance and the 3:1 mean motion resonance. I use an improved integration algorithm, and study the evolution of more particles. I confirm that the afternoon excess of the ordinary chondrites is consistent with transport from the 3:1 resonance.     

主带小行星的动力学模拟

采用接近真实太阳系的动力学模型,对主带小行星的动力学演化进行了数值模拟。计算的起始时间是儒略日JD=2.4540005×10~6,计算的时间长度为100万年。力学模型采用n+m体模型,计算程序基于小行星轨道演化的软件Orbit9。对演化结果进行分析可以发现测试粒子与木星的平运动共振对测试粒子稳定性的不同作用,以及在2:3、3:4共振处不同初始ω值对测试粒子演化结果的影响。     

Origin of chaos in the Prometheus-Pandora system

We demonstrate that the chaotic orbits of Prometheus and Pandora are due to interactions associated with the 121:118 mean motion resonance. Differential precession splits this resonance into a quartet of components equally spaced in frequency. Libration widths of the individual components exceed the splitting, resulting in resonance overlap which causes the chaos. Mean motions of Prometheus and Pandora wander chaotically in zones of width 1.8 and 3.1 deg yr⁻¹, respectively. A model with 1.5 degrees of freedom captures the essential features of the chaotic dynamics. We use it to show that the Lyapunov exponent of 0.3 yr⁻¹ arises because the critical argument of the dominant member of the resonant quartet makes approximately two separatrix crossings every 6.2 year precessional cycle.     

Meteorite transport—Revisited

Meteorites are delivered from the asteroid belt by way of chaotic zones (Wisdom 1985a). The dominant sources are believed to be the chaotic zones associated with the ν₆ secular resonance, the 3:1 mean motion resonance, and the 5:2 mean motion resonance. Though the meteorite transport process has been previously studied, those studies have limitations. Here I reassess the meteorite transport process with fewer limitations. Prior studies have not been able to reproduce the afternoon excess (the fact that approximately twice as many meteorites fall in the afternoon as in the morning) and suggested that the afternoon excess is an observational artifact; here it is shown that the afternoon excess is in fact consistent with the transport of meteorites by way of chaotic zones in the asteroid belt. By studying models with and without the inner planets it is found that the inner planets significantly speed up the transport of meteorites.     

Speckle interferometric observations of Pluto and Charon

We report speckle interferometric observations of Pluto and its moon (1978 P1) Charon obtained on 5 June 1980 with a single 1.8-m mirror of the Multiple Mirror Telescope. Our observations yield a separation of 0″.31 (±0″.05) between Pluto and Charon at position angle 285° (±7°) for JD 2444395.75. This result and other direct observations indicate an adjustment of +4.0 hr to the orbital epoch of R. S. Harrington and J. W. Christy [Astron.J.86, 442–443 (1981)]. Our observation, which represents the first resolution of the system near minimum separation, also suggests that the inclination of the orbit to the plane of the sky should be increased by 3°; this will delay the onset of the predicted eclipsee season by one apparition to 1984 or 1985. Our data are consistent with Pluto diameter 0″.14 (±0″.02) = 3000 (±400) km and Charon diameter 0″..05 (±0″.03) = 1100 (±600) km.     

On the orbits and masses of the satellites of the Pluto-Charon system

Two small satellites of Pluto, S/2005 P1 (hereafter P1) and S/2005 P2 (hereafter P2), have recently been discovered outside the orbit of Charon, and their orbits are nearly circular and nearly coplanar with that of Charon. Because the mass ratio of Charon-Pluto is ∼0.1, the orbits of P2 and P1 are significantly non-Keplerian even if P2 and P1 have negligible masses. We present an analytic theory, with P2 and P1 treated as test particles, which shows that the motion can be represented by the superposition of the circular motion of a guiding center, the forced oscillations due to the non-axisymmetric components of the potential rotating at the mean motion of Pluto-Charon, the epicyclic motion, and the vertical motion. The analytic theory shows that the azimuthal periods of P2 and P1 are shorter than the Keplerian orbital periods, and this deviation from Kepler's third law is already detected in the unperturbed Keplerian fit of Buie and coworkers. In this analytic theory, the periapse and ascending node of each of the small satellites precess at nearly equal rates in opposite directions. From direct numerical orbit integrations, we show the increasing influence of the proximity of P2 and P1 to the 3:2 mean-motion commensurability on their orbital motion as their masses increase within the ranges allowed by the albedo uncertainties. If the geometric albedos of P2 and P1 are high and of order of that of Charon, the masses of P2 and P1 are sufficiently low that their orbits are well described by the analytic theory. The variation in the orbital radius of P2 due to the forced oscillations is comparable in magnitude to that due to the best-fit Keplerian eccentricity, and there is at present no evidence that P2 has any significant epicyclic eccentricity. However, the orbit of P1 has a significant epicyclic eccentricity, and the prograde precession of its longitude of periapse with a period of 5300 days should be easily detectable. If the albedos of P2 and P1 are as low as that of comets, the large inferred masses induce significant short-term variations in the epicyclic eccentricities and/or periapse longitudes on the 400-500-day timescales due to the proximity to the 3:2 commensurability. In fact, for the maximum inferred masses, P2 and P1 may be in the 3:2 mean-motion resonance, with the resonance variable involving the periapse longitude of P1 librating. Observations that sample the orbits of P2 and P1 well on the 400-500-day timescales should provide strong constraints on the masses of P2 and P1 in the near future.     

Atlas of the mean motion resonances in the Solar System

The aim of this work is to present a systematic survey of the strength of the mean motion resonances (MMRs) in the Solar System. We know by applying simple formulas where the resonances with the planets are located but there is no indication of the strength that these resonances have. We propose a numerical method for the calculation of this strength and we present an atlas of the MMRs constructed with this method. We found there exist several resonances unexpectedly strong and we look and find in the small bodies population several bodies captured in these resonances. In particular in the inner Solar System we find one asteroid in the resonance 6:5 with Venus, five asteroids in resonance 1:2 with Venus, three asteroids in resonance 1:2 with Earth and six asteroids in resonance 2:5 with Earth. We find some new possible co-orbitals of Earth, Mars, Saturn, Uranus and Neptune. We also present a discussion about the behavior of the resonant disturbing function and where the stable equilibrium points can be found at low and high inclination resonant orbits.     

The orbital evolution of NEA 30825 1900 TG1

The orbital evolution of the near-Earth asteroid (NEA) 30825 1990 TG1 has been studied by numerical integration of the equations of its motion over the 100 000-year time interval with allowance for perturbations from eight major planets and Pluto, and the variations in its osculating orbit over this time interval were determined. The numerical integrations were performed using two methods: the Bulirsch-Stoer method and the Everhart method. The comparative analysis of the two resulting orbital evolutions of motion is presented for the time interval examined. The evolution of the asteroid motion is qualitatively the same for both variants, but the rate of evolution of the orbital elements is different. Our research confirms the known fact that the application of different integrators to the study of the long-term evolution of the NEA orbit may lead to different evolution tracks.     

Dynamical capture in the Pluto–Charon system

This paper explores the possibility that the progenitors of the small satellites of Pluto got captured in the Pluto?CCharon system from the massive heliocentric planetesimal disk in which Pluto was originally embedded into. We find that, if the dynamical excitation of the disk is small, temporary capture in the Pluto?CCharon system can occur with non- negligible probability, due to the dynamical perturbations exerted by the binary nature of the Pluto?CCharon pair. However, the captured objects remain on very elliptic orbits and the typical capture time is only ~ 100?years. In order to explain the origin of the small satellites of Pluto, we conjecture that some of these objects got disrupted during their Pluto-bound phase by a collision with a planetesimal of the disk. This could have generated a debris disk, which damped under internal collisional evolution, until turning itself into an accretional disk that could form small satellites on circular orbits, co-planar with Charon. Unfortunately, we find that objects large enough to carry a sufficient amount of mass to generate the small satellites of Pluto have collisional lifetimes orders of magnitude longer than the capture time. Thus, this scenario cannot explain the origin of the small satellites of Pluto, which remains elusive.