Mapping the three-body system – decay time and reversibility

In this paper we carry out a quantitative analysis of the three-body systems and map them as a function of decaying time and initial configuration, look at this problem as an example of a simple deterministic system and ask to what extent the orbits are really predictable. We have investigated the behaviour of about 200 000 general Newtonian three-body systems using the simplest initial conditions. Within our resolution these cover all the possible states where the objects are initially at rest and have no angular momentum. We have determined the decay time-scales of the triple systems and show that the distribution of this parameter is fractal in appearance. Some areas that appear stable on large scales exhibit very narrow strips of instability and the overall pattern, dominated by resonances, reminds us of a traditional Maasai warrior shield. Also an attempt is made to recover the original starting configuration of the three bodies by backward integration. We find there are instances where the evolution to the future and to the past lead to different orbits, in spite of time symmetric initial conditions. This implies that even in simple deterministic systems there exists an arrow of time.     

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.     

High-inclination planets and asteroids in multistellar systems

The Kozai mechanism often destabilizes high-inclination orbits. It couples changes in the eccentricity and inclination, and drives high inclination, circular orbits to low inclination, eccentric orbits. In a recent study of the dynamics of planetesimals in the quadruple star system HD 98800, there were significant numbers of stable particles in circumbinary polar orbits about the inner binary pair which are apparently able to evade the Kozai instability.
Here, we isolate this feature and investigate the dynamics through numerical and analytical models. The results show that the Kozai mechanism of the outer star is disrupted by a nodal libration induced by the inner binary pair on a shorter time-scale. By empirically modelling the period of the libration, a criteria for determining the high-inclination stability limits in general triple systems is derived. The nodal libration feature is interesting and, although affecting inclination and node only, shows many parallels to the Kozai mechanism. This raises the possibility that high-inclination planets and asteroids may be able to survive in multistellar systems.     

New methods for large dynamic range problems in planetary formation

Modern N -body techniques for planetary dynamics are generally based on symplectic algorithms specially adapted to the Kepler problem. These methods have proven very useful in studying planet formation, but typically require the time-step for all objects to be set to a small fraction of the orbital period of the innermost body. This computational expense can be prohibitive for even moderate particle number for many physically interesting scenarios, such as recent models of the formation of hot exoplanets, in which the semimajor axis of possible progenitors can vary by orders of magnitude. We present new methods which retain most of the benefits of the standard symplectic integrators but allow for radial zones with distinct time-steps. These approaches should make simulations of planetary accretion with large dynamic range tractable. As proof-of-concept, we present preliminary science results from an implementation of the algorithm as applied to an oligarchic migration scenario for forming hot Neptunes.     

Some properties of a two-body system under the influence of the Galactic tidal field

In this paper the effect of the Galactic tidal field on a Sun–comet pair will be considered when the comet is situated in the Oort cloud and planetary perturbations can be neglected. First, two averaged models were created, one of which can be solved analytically in terms of Jacobi elliptic functions. In the latter system, switching between libration and circulation of the argument of perihelion is prohibited. The non-averaged equations of motion are integrated numerically in order to determine the regions of the ( e ,  i ) phase space in which chaotic orbits can be found, and an effort is made to explain why the chaotic orbits manifest in these regions only. It is evident that for moderate values of semimajor axis, a ∼50 000 au , chaotic orbits are found for ( e <0.15 , 40°≤ i ≤140°) as determined by integrating the evolution of the comet over a period of 10⁴ orbits. These regions of chaos increase in size with increasing semimajor axis. The typical e-folding times for these orbits range from around 600 Myr to 1 Gyr and thus are of little practical interest, as the time-scales for chaos arising from passing stars are much shorter. As a result of Galactic rotation, the chaotic regions in ( e ,  i ) phase space are not symmetric for prograde and retrograde orbits.     

The resonant structure of Jupiter's Trojan asteroids – I. Long-term stability and diffusion

We study the global dynamics of the jovian Trojan asteroids by means of the frequency map analysis. We find and classify the main resonant structures that serve as skeleton of the phase space near the Lagrangian points. These resonances organize and control the long-term dynamics of the Trojans. Besides the secondary and secular resonances, that have already been found in other asteroid sets in mean motion resonance (e.g. main belt, Kuiper belt), we identify a new type of resonance that involves secular frequencies and the frequency of the great inequality, but not the libration frequency. Moreover, this new family of resonances plays an important role in the slow transport mechanism that drives Trojans from the inner stable region to eventual ejections. Finally, we relate this global view of the dynamics with the observed Trojans, identify the asteroids that are close to these resonances and study their long-term behaviour.     

On the habitability of the OGLE-2006-BLG-109L planetary system

We investigate the dynamics of putative Earth-mass planets in the habitable zone (HZ) of the extrasolar planetary system OGLE-2006-BLG-109L, a close analogue of the Solar system. Our work is inspired by the work of Malhotra & Minton. Using the linear Laplace–Lagrange theory, they identified a strong secular resonance that may excite large eccentricity of orbits in the HZ. However, due to uncertain or unconstrained orbital parameters, the subsystem of Jupiters may be found in a dynamically active region of the phase space spanned by low-order mean-motion resonances. To generalize this secular model, we construct a semi-analytical averaging method in terms of the restricted problem. The secular orbits of large planets are approximated by numerically averaged osculating elements. They are used to calculate the mean orbits of terrestrial planets by means of a high-order analytic secular theory developed in our previous works. We found regions in the parameter space of the problem in which stable, quasi-circular orbits in the HZ are permitted. The excitation of eccentricity in the HZ strongly depends on the apsidal angle of jovian orbits. For some combinations of that angle, eccentricities and semimajor axes consistent with the observations, a terrestrial planet may survive in low eccentric orbits. We also study the effect of post-Newtonian gravity correction on the innermost secular resonance.     

Periodic variation of Oort Cloud flux and cometary impacts on the Earth and Jupiter

Time variation in impact probability is studied by assuming that the periodic flux of the Oort Cloud comets within 15 au arises from the motion of the Sun with respect to the Galactic mid-plane. The periodic flux clearly shows up in the impact rate of the captured Oort Cloud cometary population, with a phase shift caused by the orbital evolution. Depending on the assumed flux of comets and the size distribution of comets, the impact rate of the Oort Cloud comets of 1 km in diameter or greater is from 5 to 700 impacts Myr⁻¹ on the Earth and from 0.5 to 70 impacts per 1000 yr on Jupiter. The relative fractions of impacts are 0.09, 0.11, 0.26 and 0.54 for long-period comets, Halley type comets, Jupiter family comets and near-Earth objects, respectively. For Jupiter, the corresponding fractions in the first three categories are 0.18, 0.31 and 0.51. If we consider physical fading of comet activity that is compatible with the observations, then the impact rates of active comets are two orders of magnitude smaller than the total impact rates by all kinds of comets and cometary asteroids of size 1 km or greater.     

An extension of Newton's apsidal precession theorem

Newton's apsidal precession theorem in Proposition 45 of Book I of the 'Principia' has great mathematical, physical, astronomical and historical interest. The lunar theory and the precession of the perihelion of the planet Mercury are but two examples of the applications of this theorem. We have examined the precession of orbits under varying force laws as measured by the apsidal angle θ( N , e ), where N is the index for the centripetal force law, for varying eccentricity e . The paper derives a general function for the apsidal angle, dependent only on e and N as the potential is spherically symmetric. Further, we explore approximate ways of the solution of this equation, in the neighbourhood of   N = 2  which happens to be the case of greatest historical interest. Exact solutions are derived where they are possible. The first derivatives  ∂θ/∂ N   and  ∂θ/∂ h   [where h ( N , e ) is the angular momentum] are analytically expressed in the neighbourhood of   N = 2  (case of the inverse square law). The value of  ∂θ/∂ N   is computed numerically as well for  1 ≤ N < 3  . The resulting integrals are interesting improper integrals with singularities at both limits. Some of the integrals, especially for   N = 2  , can be given in closed form in terms of generalized hypergeometric functions which are reducible in terms of algebraic and logarithmic functions. No evidence was found for isolated cases of zero precession as e was increased. The   N = 1  case of the logarithmic potential is also briefly discussed in view of its interest for the dynamics of eccentric orbits and its relevance to realistic galaxy models. The possibility of apsidal precession was also examined for a few cases of high-eccentricity asteroids and extrasolar planets. We find that these systems may provide interesting new laboratories for studies of gravity.