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 fast Fourier transform technique, we construct high-resolution dynamical maps on the plane of initial semimajor axis a ₀ versus inclination i ₀. These maps show three most stable regions, with i ₀ in the range of  (0°, 12°), (22°, 36°)  and  (51°, 59°),  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 ₄ and L ₅ 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  ν₈  secular resonance around   i ₀∼ 44°  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°)  and an unstable gap around  44°  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 ₀, i ₀)  . The fine structures in the dynamical maps can be explained by these secular resonances.     

冥族小天体与海王星特洛伊的近密交会

海王星外天体中的冥族小天体与海王星成2:3的平运动轨道共振,且具有较大的轨道偏心率,因此它们能与海王星特洛伊的轨道发生重叠,导致近密交会和碰撞,从而深刻地影响两者的动力学演化。利用数值模拟的方法,有效地获得了这两群小天体间近密交会的信息,讨论了可能影响两者近密交会频率的因素,包括小天体质量、轨道倾角和轨道偏心率等。在合理近似条件下,建立了估算两群小天体近密交会和碰撞次数的理论公式。结合已有的数值模拟结果,以及对冥族小天体观测数据的分析,对实际情况下冥族小天体群与典型特洛伊小天体之间的近密交会和碰撞次数进行了估算,证明近密交会较为频繁地发生,而碰撞则极其罕见,并且各尺寸范围的小天体对近密交会和碰撞次数的贡献各有不同。这一套分析和估算的方法可以直接应用在其他类似小天体间交会过程的估算上。     

Neptune Trojans: a Revisit and a New Membertwo

Up to now, 17 Neptune Trojan asteroids have been detected with their orbits being well determined by continuous observations. This paper analyzes systematically their orbital dynamics. Our results show that except for two temporary members with relatively short lifespans on Trojan orbits, the vast majority of Neptune Trojans located within their orbital uncertainties may survive in the solar system age. The escaping probability of Neptune Trojans, through slow diffusion in the orbital element space in 4.5 billion years, is estimated to be ～50%. The asteroid 2012 UW177 classified as a Centaur asteroid by the IAU Minor Planet Center currently is in fact a Neptune Trojan. Numerical simulations indicate that it is librating on the tadpole-shaped orbit around the Neptune's L₄ point. It was captured into the current orbit approximately 0.23 million years ago, and will stay there for at least another 1.3 million years in the future. Its high inclination of i ≈ 54^° not only makes it the most inclined Neptune Trojan, but also makes it exhibit the complicated and interesting co-orbital transitions between the leading and trailing Trojans via the quasi-satellite orbit phase.     

The dynamics of two massive planets on inclined orbits

The significant orbital eccentricities of most giant extrasolar planets may have their origin in the gravitational dynamics of initially unstable multiple planet systems. In this work, we explore the dynamics of two close planets on inclined orbits through both analytical techniques and extensive numerical scattering experiments. We derive a criterion for two equal mass planets on circular inclined orbits to achieve Hill stability, and conclude that significant radial migration and eccentricity pumping of both planets occurs predominantly by 2:1 and 5:3 mean motion resonant interactions. Using Laplace-Lagrange secular theory, we obtain analytical secular solutions for the orbital inclinations and longitudes of ascending nodes, and use those solutions to distinguish between the secular and resonant dynamics which arise in numerical simulations. We also illustrate how encounter maps, typically used to trace the motion of massless particles, may be modified to reproduce the gross instability seen by the numerical integrations. Such a correlation suggests promising future use of such maps to model the dynamics of more coplanar massive planet systems.     

Study and application of the resonant secular dynamics beyond Neptune

We use a secular representation to describe the long-term dynamics of transneptunian objects in mean-motion resonance with Neptune. The model applied is thoroughly described in Saillenfest et al. (Celest Mech Dyn Astron, doi: 10.1007/s10569-016-9700-5, 2016). The parameter space is systematically explored, showing that the secular trajectories depend little on the resonance order. High-amplitude oscillations of the perihelion distance are reported and localised in the space of the orbital parameters. In particular, we show that a large perihelion distance is not a sufficient criterion to declare that an object is detached from the planets. Such a mechanism, though, is found unable to explain the orbits of Sedna or \(2012\text {VP}_{113}\), which are insufficiently inclined (considering their high perihelion distance) to be possibly driven by such a resonant dynamics. The secular representation highlights the existence of a high-perihelion accumulation zone due to resonances of type 1:k with Neptune. That region is found to be located roughly at \(a\in [100;300]\) AU, \(q\in [50;70]\) AU and \(I\in [30;50]^{\circ }\). In addition to the flux of objects directly coming from the Scattered Disc, numerical simulations show that the Oort Cloud is also a substantial source for such objects. Naturally, as that mechanism relies on fragile captures in high-order resonances, our conclusions break down in the case of a significant external perturber. The detection of such a reservoir could thus be an observational constraint to probe the external Solar System.     

The formation of Uranus and Neptune in solid-rich feeding zones: Connecting chemistry and dynamics

The core accretion theory of planet formation has at least two fundamental problems explaining the origins of Uranus and Neptune: (1) dynamical times in the trans-saturnian solar nebula are so long that core growth can take >15 Myr and (2) the onset of runaway gas accretion that begins when cores reach ∼10M_⊕ necessitates a sudden gas accretion cutoff just as Uranus and Neptune’s cores reach critical mass. Both problems may be resolved by allowing the ice giants to migrate outward after their formation in solid-rich feeding zones with planetesimal surface densities well above the minimum-mass solar nebula. We present new simulations of the formation of Uranus and Neptune in the solid-rich disk of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) using the initial semimajor axis distribution of the Nice model (Gomes, R., Levison, H.F., Tsiganis, K., Morbidelli, A. [2005]. Nature 435, 466-469; Morbidelli, A., Levison, H.F., Tsiganis, K., Gomes, R. [2005]. Nature 435, 462-465; Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F. [2005]. Nature 435, 459-461), with one ice giant forming at 12 AU and the other at 15 AU. The innermost ice giant reaches its present mass after 3.8-4.0 Myr and the outermost after 5.3-6 Myr, a considerable time decrease from previous one-dimensional simulations (e.g. Pollack, J.B., Hubickyj, O., Bodenheimer, P., Lissauer, J.J., Podolak, M., Greenzweig, Y. [1996]. Icarus 124, 62-85). The core masses stay subcritical, eliminating the need for a sudden gas accretion cutoff.Our calculated carbon mass fractions of 22% are in excellent agreement with the ice giant interior models of Podolak et al. (Podolak, M., Weizman, A., Marley, M. [1995]. Planet. Space Sci. 43, 1517-1522) and Marley et al. (Marley, M.S., Gómez, P., Podolak, M. [1995]. J. Geophys. Res. 100, 23349-23354). Based on the requirement that the ice giant-forming planetesimals contain >10% mass fractions of methane ice, we can reject any Solar System formation model that initially places Uranus and Neptune inside of Saturn’s orbit. We also demonstrate that a large population of planetesimals must be present in both ice giant feeding zones throughout the lifetime of the gaseous nebula. This research marks a substantial step forward in connecting both the dynamical and chemical aspects of planet formation. Although we cannot say that the solid-rich solar nebula model of Dodson-Robinson et al. (Dodson-Robinson, S.E., Willacy, K., Bodenheimer, P., Turner, N.J., Beichman, C.A. [2009]. Icarus 200, 672-693) gives exactly the appropriate initial conditions for planet formation, rigorous chemical and dynamical tests have at least revealed it to be a viable model of the early Solar System.     

The fate of primordial lunar Trojans

Multiple large impact basins on the lunar nearside formed in a relatively-short interval around 3.8-3.9 Gyr ago, in what is known as the Lunar Cataclysm (LC; also known as Late Heavy Bombardment). It is widely thought that this impact bombardment has affected the whole Solar System or at least all the inner planets. But with non-lunar evidence for the cataclysm being relatively weak, a geocentric cause of the Lunar Cataclysm cannot yet be completely ruled out [Ryder, G., 1990. Eos 71, 313, 322-323]. In principle, late destabilization of an additional Earth satellite could result in its tidal disruption during a close lunar encounter (cf. [Asphaug, E., Agnor, C.B., Williams, Q., 2006. Nature 439, 155-160]). If the lost satellite had D>500 km, the resulting debris can form multiple impact basins in a relatively short time, possibly explaining the LC. Canup et al. [Canup, R.M., Levison, H.F., Stewart, G.R., 1999. Astron. J. 117, 603-620] have shown that any additional satellites of Earth formed together with (and external to) the Moon would be unable to survive the rapid initial tidally-driven expansion of lunar orbit. Here we explore the fate of objects trapped in the lunar Trojan points, and find that small lunar Trojans can survive the Moon's orbital evolution until they and the Moon reach 38 Earth radii, at which point they are destabilized by a strong solar resonance. However, the dynamics of Trojans containing enough mass to cause the LC (diameters >150 km) is more complex; we find that such objects do not survive the passage through a weaker solar resonance at 27 Earth radii. This distance was very likely reached by the Moon long before the LC, which seems to rule out the disruption of lunar Trojans as a cause of the LC.     

Dynamics of Mars Trojans

In this paper we explore the dynamical stability of the Mars Trojan region applying mainly Laskar's Frequency Map Analysis. This method yields the chaotic diffusion rate of orbits and allows to determine the most stable regions. It also gives the frequencies which are responsible for the instability of orbits. The most stable regions are found for inclinations between about 15° and 30°. For inclinations smaller than 15°, we confirm, by applying a synthetic secular theory, that the secular resonances ν₃, ν₄, ν₁₃, ν₁₄ rapidly excite asteroid orbits within a few Myrs, or even faster. The asteroids are removed from the Trojan region after a close encounter with Mars. For large inclinations, the secular resonance ν₅ clears a small region around 30° while the Kozai resonance rapidly removes bodies for inclinations larger than 35°. The dynamical lifetimes of the three L5 Trojans, (5261) Eureka, 1998 VF31, 2001 DH47, and the only L4 Trojan 1999 UJ7 are determined by numerically integrating clouds of corresponding clones over the age of the Solar System. All four Trojans reside in the most stable region with smallest diffusion coefficients. Their dynamical half-lifetime is of the order of the age of the Solar System. The Yarkovsky force has little effect on the known Trojans but for bodies smaller than about 1-5 m the drag is strong enough to destabilize Trojans on a timescale shorter than 4.5 Gyr.     

A peculiar family of Jupiter Trojans: The Eurybates

The Eurybates family is a compact core inside the Menelaus clan, located in the L₄ swarm of Jupiter Trojans. Fornasier et al. (Fornasier, S., Dotto, E., Hainaut, O., Marzari, F., Boehnhardt, H., De Luise, F., Barucci, M.A. [2007]. Icarus 190, 622-642) found that this family exhibits a peculiar abundance of spectrally flat objects, similar to Chiron-like Centaurs and C-type main belt asteroids. On the basis of the visible spectra available in literature, Eurybates family’s members seemed to be good candidates for having on their surfaces water/water ice or aqueous altered materials.To improve our knowledge of the surface composition of this peculiar family, we carried out an observational campaign at the Telescopio Nazionale Galileo (TNG), obtaining near-infrared spectra of 7 members. Our data show a surprisingly absence of any spectral feature referable to the presence of water, ices or aqueous altered materials on the surface of the observed objects. Models of the surface composition are attempted, evidencing that amorphous carbon seems to dominate the surface composition of the observed bodies and some amount of silicates (olivine) could be present.     

The rotational periods of Uranus and Neptune

Observations of tilts of spectral lines in the spectrum of Uranus and Neptune yield the following rotational periods: “Uranus,” 24 ± 3 hr; “Neptune,” 22 ± 4 hr. Neptune is confirmed to rotate in a direct sense. The position angle of the pole of Uranus, projected onto the plane of the sky, is found to be 283 ± 4°. The value for Neptune is 32 ± 11°. These results agree with the direction of the pole of Uranus inferred from the common plane of its four brightest satellites and with the direction of the pole of Neptune as inferred from the precession of Triton's orbit. The rotational period of Uranus is found to be consistent with modern values of its optical and dynamical oblateness and the theory of solid-body rotation with hydrostatic equilibrium. This is barely the case for the period derived for Neptune and we suspect that future observations made under better seeing conditions may lead to a shorter rotation period between 15 and 18 hr. Because of a substantial difference between our results and those of earlier spectroscopic and photometric investigations we include an assessment of several previously published photometric studies and a new reduction of the original Lowell and Slipher spectroscopic plates of Uranus [Lowell Obs. Bull. 2, 17–18, 19–20 (1912)]. The early visual photometry of Campbell (Uranus) and Hall (Neptune) is found to be more satisfactorily accounted for by periods of 21.6 and 23.1 hr, respectively, than by the periods originally suggested by the observers. Our reduction of the Lowell and Slipher Uranus plates yields a period near 33 hr uncorrected for seeing. This value is consistent with the results based on the 4-m echelle date.     

Variability of Neptune

Earth-based observers of Neptune have found that the planet varies in brightness at various wavelengths in ways that suggest that changes occur in the planet's atmosphere on several different time scales. Global inhomogeneities in high-altitude haze distribution that are stable for several days permit measurements of the planet's rotation period (about 18 hr), but this stability sometimes breaks down, obscuring the diurnal lightcurve. In addition, there is an apparent long-term variability of the brightness of Neptune in anticorrelation with the cycle of solar activity. This slow variability of low amplitude may be punctuated by outburst of high-altitude condensation of particles in the atmosphere whose decay time is several months.     

Spectroscopy and photometry of Mars Trojans

Mars is the only terrestrial planet known to have co-orbiting “Trojan” asteroids. We have obtained visible and near-IR reflectance spectra of three of these objects: 5261 Eureka and 1998 VF31 in the L5 region and 1999 UJ7 in the L4 region. We also obtained JHK spectrophotometry and a visible lightcurve for 5261 Eureka. The asteroid 5261 Eureka has a visible spectrum that is classified as Sr in the Bus taxonomy, and has infrared colors consistent with the A-class asteroids. The data for 1998 VF31 have a restricted wavelength range, but are most consistent with the Sr or Sa class, though we note a marginal consistency with the D class. We can rule out a C-class classification. 1999 UJ7 has an X-class or T-class spectrum, which is unlike that of the other two Mars Trojans. The photometric data for Eureka are limited, but we can constrain the period to longer than 5 hours (likely 5.5-6 hours) and lightcurve amplitude of at least 0.15 magnitude at this viewing geometry. The spectral differences among the Mars Trojans suggests that either they did not all form at their present solar distances or that they have not always been at their present sizes.     

On the prediction of Neptune

An investigation into the problems and queries surrounding the theoretical procedures involved in the theory of John Couch Adams, in order to determine the effect of variations in the time and mean distance on the accuracy of the predictions.This work forms part of a thesis approved by the University of London for the award of the degree of Doctor of Philosophy.     

Long term stability of Earth Trojans

We explore the long-term stability of Earth Trojans by using a chaos indicator, the Frequency Map Analysis. We find that there is an extended stability region at low eccentricity and for inclinations lower than about $50^{\circ }$ even if the most stable orbits are found at $i \le 40^{\circ }$ . This region is not limited in libration amplitude, contrary to what found for Trojan orbits around outer planets. We also investigate how the stability properties are affected by the tidal force of the Earth–Moon system and by the Yarkovsky force. The tidal field of the Earth–Moon system reduces the stability of the Earth Trojans at high inclinations while the Yarkovsky force, at least for bodies larger than 10 m in diameter, does not seem to strongly influence the long-term stability. Earth Trojan orbits with the lowest diffusion rate survive on timescales of the order of $10^9$  years but their evolution is chaotic. Their behaviour is similar to that of Mars Trojans even if Earth Trojans appear to have shorter lifetimes.     

On the Instability of Jupiter's Trojans

We have numerically explored the mechanisms that destabilize Jupiter's Trojan orbits outside the stability region defined by Levison et al. (1997, Nature385, 42-44). Different models have been exploited to test various possible sources of instability on timescales on the order of ∼10⁸ years.In the restricted three-body model, only a few Trojan orbits become unstable within 10⁸ years. This intrinsic instability contributes only marginally to the overall instability found by Levison et al.In a model where the orbital parameters of both Jupiter and Saturn are fixed, we have investigated the role of Saturn and its gravitational influence. We find that a large fraction of Trojan orbits become unstable because of the direct nonresonant perturbations by Saturn. By shifting its semimajor axis at constant intervals around its present value we find that the near 5:2 mean motion resonance between the two giant planets (the Great Inequality) is not responsible for the gross instability of Jupiter's Trojans since short-term perturbations by Saturn destabilize Trojans, even when the two planets are far out of the resonance.Secular resonances are an additional source of instability. In the full six-body model with the four major planets included in the numerical integration, we have analyzed the effects of secular resonances with the node of the planets. Trojan asteroids have relevant inclinations, and nodal secular resonances play an important role. When a Trojan orbit becomes unstable, in most cases the libration amplitude of the critical argument of the 1:1 mean motion resonance grows until the asteroid encounters the planet. Libration amplitude, eccentricity, and nodal rate are linked for Trojan orbits by an algebraic relation so that when one of the three parameters is perturbed, the other two are affected as well. There are numerous secular resonances with the nodal rate of Jupiter that fall inside the region of instability and contribute to destabilize Trojans, in particular the ν₁₆. Indeed, in the full model the escape rate over 50 Myr is higher compared to the fixed model.Some secular resonances even cross the stability region delimited by Levison et al. and cause instability. This is the case of the 3:2 and 1:2 nodal resonances with Jupiter. In particular the 1:2 is responsible for the instability of some clones of the L4 Trojan (3540) Protesilaos.     

Trojans in Stable Chaotic Motion

The orbits of 13 Trojan asteroids have been calculated numerically in the model of the outer solar system for a time interval of 100 million years. For these asteroids Milani et al. (1997) determined Lyapunov times less than 100 000 years and introduced the notion "asteroids in stable chaotic motion". We studied the dynamical behavior of these Trojan asteroids (except the asteroid Thersites which escaped after 26 million years) within 11 time intervals - i.e. subintervals of the whole time - by means of: (1) a numerical frequency analysis (2) the root mean square (r.m.s.) of the orbital elements and (3) the proper elements. For each time interval we compared the root mean squares of the orbital elements (a, e and i) with the corresponding proper element. It turned out that the variations of the proper elements e_p in the different time intervals are correlated with the corresponding r.m.s.(e); this is not the case for sin I_p with r.m.s.(i). This revised version was published online in July 2006 with corrections to the Cover Date.     

The satellites of Neptune and the origin of Pluto

Pluto and the chaotic satellite system of Neptune may have originated from a single encounter of Neptune with a massive solar system body. A series of numerical experiments has been carried out to try to set limits on the circumstances of such an encounter. These experiments show that orbits very much like those of Pluto, Triton, and Nereid can result from a single close encounter of such a body with Neptune. The implied mass range and encounter velocities limit the source of the encountering body to a former trans-Neptunian planet in the 2- to 5-Earth-mass range.