Dynamical evolution of a cometary swarm in the outer planetary region

The dynamical evolution of bodies under the gravitational influence of the accreting proto-Uranus and proto-Neptune is investigated. The main aim of this study is to analyze the interrelations between the accretion of Uranus and Neptune with other processes of cosmological importance as, for example, the formation of a cometary reservoir from bodies placed into near-parabolic orbits by planetary perturbations and the scattering of bodies to the region of the terrestrial planets. Starting with a mass ratio (initial mass/present mass) of 0.1, Uranus and Neptune acquire masses close to their present ones in a time scale of 10⁸ years. Neptune is found to be the most important contributor of comets to the cometary reservoir. The time scale of bodies scattered by Neptune to reach near-parabolic orbits (semimajor axes a > 10⁴ AU)is about 10⁹ years. The contribution of Uranus was partially inhibited because a large part of the residual bodies of its accretion zone fell under the strong gravitational influence of Jupiter and Saturn. A significant fraction of the bodies dispersed by Uranus and Neptune reached the region of the terrestrial planets in a time scale of some 10⁸ years.     

Apsidal precession of the outer solar planetary orbits due to the pioneer anomaly

Despite the now common position that the Pioneer anomaly is not a real gravitational effect but an effect due to the on-board thermal recoil forces – for curiosity’s sake, we here take the suggestion of Nyambuya (2015) where it has been assumed that the Pioneer anomaly – can, in-principle, be attributed to a gravitational effect due to these spacecrafts accreting some material from a rarefied Interplanetary Medium (IPM) in the domain where the Pioneer anomaly has manifested [20 AU ≲ r ≲ 70 AU]. If this assumption is correct, then, the expected Pioneer acceleration of these spacecrafts maybe much smaller than the Pioneer acceleration to cause as noticeable apsidal precession of the outer Solar planets Uranus, Neptune and Pluto, thus making it difficult to rule-out a gravitational origin of the Pioneer anomaly.     

The 1/1 resonance in extrasolar systems

We present families of symmetric and asymmetric periodic orbits at the 1/1 resonance, for a planetary system consisting of a star and two small bodies, in comparison to the star, moving in the same plane under their mutual gravitational attraction. The stable 1/1 resonant periodic orbits belong to a family which has a planetary branch, with the two planets moving in nearly Keplerian orbits with non zero eccentricities and a satellite branch, where the gravitational interaction between the two planets dominates the attraction from the star and the two planets form a close binary which revolves around the star. The stability regions around periodic orbits along the family are studied. Next, we study the dynamical evolution in time of a planetary system with two planets which is initially trapped in a stable 1/1 resonant periodic motion, when a drag force is included in the system. We prove that if we start with a 1/1 resonant planetary system with large eccentricities, the system migrates, due to the drag force, along the family of periodic orbits and is finally trapped in a satellite orbit. This, in principle, provides a mechanism for the generation of a satellite system: we start with a planetary system and the final stage is a system where the two small bodies form a close binary whose center of mass revolves around the star.     

Mass of the Kuiper belt

The Kuiper belt includes tens of thousand of large bodies and millions of smaller objects. The main part of the belt objects is located in the annular zone between 39.4 and 47.8 au from the Sun; the boundaries correspond to the average distances for orbital resonances 3:2 and 2:1 with the motion of Neptune. One-dimensional, two-dimensional, and discrete rings to model the total gravitational attraction of numerous belt objects are considered. The discrete rotating model most correctly reflects the real interaction of bodies in the Solar system. The masses of the model rings were determined within EPM2017—the new version of ephemerides of planets and the Moon at IAA RAS—by fitting spacecraft ranging observations. The total mass of the Kuiper belt was calculated as the sum of the masses of the 31 largest trans-Neptunian objects directly included in the simultaneous integration and the estimated mass of the model of the discrete ring of TNO. The total mass is \((1.97 \pm 0.35)\times 10^{-2} \ m_{\oplus }\). The gravitational influence of the Kuiper belt on Jupiter, Saturn, Uranus, and Neptune exceeds at times the attraction of the hypothetical 9th planet with a mass of \(\sim 10 \ m_{\oplus }\) at the distances assumed for it. It is necessary to take into account the gravitational influence of the Kuiper belt when processing observations and only then to investigate residual discrepancies to discover a possible influence of a distant large planet.     

Possible relation between galactic flat rotational curves and the Pioneers’ anomalous acceleration

We consider a generic minimal modification of the Newtonian potential, that is a modification that introduces only one additional dimensional parameter. The modified potential depends on a function whose behavior for large and small distances can be fixed in order to obtain: (i) galactic flat rotational curves and (ii) a universal constant acceleration independent of the masses of the interacting bodies (Pioneer anomaly). Then using a dimensional argument we show that the Tully–Fisher relation for the maximal rotational velocity of spiral galaxies follows without any further assumptions. This result suggests that the Pioneer anomalous acceleration and the flat rotational curves of galaxies could have a common origin in a modified gravitational theory. The relation of these results with the Modified Newtonian Dynamics (MOND) is discussed.     

A Model for the Common Origin of Jupiter Family and Halley Type Comets

A numerical simulation of the Oort cloud is used to explain the observed orbital distributions and numbers of Jupiter-family (JF) and Halley-type (HT) short-period (SP) comets. Comets are given initial orbits with perihelion distances between 5 and 36 au, and evolve under planetary, stellar and Galactic perturbations for 4.5 Gyr. This process leads to the formation of an Oort cloud (which we define as the region of semimajor axes a > 1,000 au), and to a flux of cometary bodies from the Oort cloud returning to the planetary region at the present epoch. The results are consistent with the dynamical characteristics of SP comets and other observed cometary populations: the near-parabolic flux, Centaurs, and high-eccentricity trans-Neptunian objects. To achieve this consistency with observations, the model requires that the number of comets versus initial perihelion distance is concentrated towards the outer planetary region. Moreover, the mean physical lifetime of observable comets in the inner planetary region (q < 2.5 au) at the present epoch should be an increasing function of the comets’ initial perihelion distances. Virtually all observed HT comets and nearly half of observed JF comets come from the Oort cloud, and initially (4.5 Gyr ago) from orbits concentrated near the outer planetary region. Comets that have been in the Oort cloud also return to the Centaur (5 < q < 28 au, a < 1,000 au) and near-Neptune high-eccentricity regions. Such objects with perihelia near Neptune are hard to discover, but Centaurs with characteristics predicted by the model (e.g. large semimajor axes, above 60 au, or high inclinations, above 40°) are increasingly being found by observers. The model provides a unified picture for the origin of JF and HT comets. It predicts that the mean physical lifetime of all comets in the region q < 1.5 au is less than ～200 revolutions.     

Embedded star clusters and the formation of the Oort cloud: II. The effect of the primordial solar nebula

This paper deals with Oort cloud formation while the Sun was in an embedded cluster and surrounded by its primordial nebula. This work is a continuation of Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82], building on the model presented therein, and adding the aerodynamic drag and gravitational potential of the primordial solar nebula. Results are presented of numerical simulations of comets subject to the gravitational influence of the Sun, Jupiter, Saturn, star cluster and primordial solar nebula; some of the simulations included the gravitational influence of Uranus and Neptune as well. The primordial solar nebula was approximated by the minimum-mass Hayashi model [Hayashi, C., Nakozawa, K., Nakagawa, Y., 1985. In: Black, D.C., Matthews, M.S. (Eds.). Protostars and Planets II. Univ. of Arizona Press, Tucson, AZ] whose inner and outer radii have been truncated at various distances from the Sun. A comet size of 1.7 km was used for most of our simulations. In all of our simulations, the density of the primordial solar nebula decayed exponentially with an e-folding time of 2 Myr. It turns out that when the primordial solar nebula extends much beyond Saturn or Neptune, virtually no material will end up in the Oort cloud (OC) during this phase. Instead, the majority of the material will be on circular orbits inside of Jupiter if the inner edge of the disk is well inside Jupiter's orbit. If the disk's inner edge is beyond Jupiter's orbit, most comets end up on orbits in exterior mean-motion resonances with Saturn when Uranus and Neptune are not present. In those cases where the outer edge of the disk is close to Saturn or Neptune, the fraction of material that ends up in the subsequently formed OC is much less than that found in Brasser et al. [Brasser, R., Duncan, M., Levison, H., 2006. Icarus 184, 59-82] for the same cluster densities. This implies that for comets of roughly 2 km in size, the presence of the primordial solar nebula hinders OC formation. A byproduct of some of our simulations are endresults with a substantial fraction of the comets in the Uranus-Neptune scattered disk. A subsequent followup of this material is planned for the near future. In order to determine the effect of the size of the comets on OC formation efficiency, a set of runs with the same initial conditions but different cometary radii have been performed as well, from which it is determined that the threshold comet size to begin producing significant Oort clouds is roughly 20 km. This implies that the presence of the primordial solar nebula acts as a size-sorting mechanism, with large bodies unaffected by the gas drag and ending up in the OC while small bodies remain trapped in the planetary region, in the models studied.     

Satellite aided capture

A two body, patched conic analysis is presented for a planetary capture mode in which a gravity assist by an existing natural satellite of the planet aids in the capture. An analytical condition sufficient for capture is developed and applied for the following planet/satellite systems: Earth/Moon, Jupiter/Ganymede, Jupiter/Callisto, Saturn/Titan and Neptune/Triton. Co-planar, circular planetary orbits are assumed. Three sources of bodies to be captured are considered: spacecraft launched from Earth, bodies entering the solar system from interstellar space, and bodies already in orbit around the Sun. Results show that the Neptune/Triton system has the most capability for satellite aided capture of those studied. It can easily capture bodies entering the Solar System from interstellar space. Its ability to capture spacecraft launched from Earth is marginal and can only be decided with better definition of physical properties. None of the other systems can capture bodies from these two sources, but all can capture bodies already in orbit around the Sun under appropriate conditions.     

What do the orbital motions of the outer planets of the Solar System tell us about the Pioneer anomaly?

In this paper, we investigate the effects that an anomalous acceleration as that experienced by the Pioneer spacecraft after they passed the 20 AU threshold would induce on the orbital motions of the Solar System planets placed at heliocentric distances of 20 AU or larger as Uranus, Neptune and Pluto. It turns out that such an acceleration, with a magnitude of 8.74 × 10⁻¹⁰ m s⁻², would affect their orbits with secular and short-period signals large enough to be detected according to the latest published results by E.V. Pitjeva, even by considering errors up to 30 times larger than those released. The absence of such anomalous signatures in the latest data rules out the possibility that in the region 20–40 AU of the Solar System an anomalous force field inducing a constant and radial acceleration with those characteristics affects the motion of the major planets.     

A model for accretion of the terrestrial planets

One possible origin of the terrestrial planets involves their formation by gravitational accretion of particles originally in Keplerian orbits about the sun. Some implications of this theory are considered. A formal expression for the rate of mass accretion by a planet is developed. The formal singularity of the gravitational collision cross-section for low relative velocities is shown to be without physical significance when the accreting bodies are in heliocentric orbits. The distribution of particle velocities relative to an accreting planet is considered; the mean velocity increases with time. The internal temperature of an accreting planet is shown to depend simply on the accretion rate. A simple and physically reasonable approximate expression for a planetary accretion rate is proposed.     

Secular Effects on Orbits of Binary Stars Induced by Temporal Variation of Gravitational Constant (Case for Elliptical Orbit)

In this paper we investigate the influence of a varying gravitation constant on the orbits of celestial bodies. Regarding the eccentric anomaly as an independent variable, we find the solutions to the perturbed equations of motion. In the first order solutions, we find the secular and periodic variations in semi-major axis. For the other orbital elements only periodic variations exhibit. However in the second order solutions, the longitude of periastron and the mean longitude have secular terms. Applying the calculations to six selected binaries, we give the numerical estimations of the variations of orbits. These results are then carefully compared and discussed.     

GLISSE: A GPU-optimized planetary system integrator with application to orbital stability calculations.

We present a GPU-accelerated numerical integrator specifically optimized for stability calculations of small bodies in planetary systems. Specifically, the integrator is designed for cases when large numbers of test particles (tens or hundreds of thousands) need to be followed for long durations (millions of orbits) to assess the orbital stability of their initially “close-encounter free” orbits. The GLISSE (Gpu Long-term Integrator for Solar System Evolution) code implements several optimizations to achieve a roughly factor of 100 speed increase over running the same code on a CPU. We explain how various hardware speed bottlenecks can be avoided by the careful code design, although some of the choices restrict the usage to specific types of application.As a first application, we study the long-term stability of small bodies initially on orbits between Uranus and Neptune. We map out in detail the small portion of the phase space in which small bodies can survive for 4.5 billion years of evolution; the ability to integrate large numbers of particles allow us to identify for the first time how instability-inducing mean-motion resonance overlaps sharply define the stable regions. As a second application, we map the boundaries of 4 Gyr stability for transneptunian objects in the 5:2 and 3:1 mean-motion resonances, demonstrating that long-term perturbations remove the initially stable Neptune-crossing members.     

Orbital resonances in the solar nebula: Implications for planetary accretion

The influence of gas drag and gravitational perturbations by a planetary embryo on the orbit of a planetesimal in the solar nebula was examined. Non-Keplerian rotation of the gas causes secular decay of the orbit. If the planetesimal's orbit is exterior to the perturber's, resonant perturbations oppose this drag and can cause it to be trapped in a stable orbit at a commensurability of order j/(j + 1), where j is an integer. Numerical and analytical demonstrations show that resonant trapping occurs for wide ranges of perturbing mass, planetesimal size, and j. Induced eccentricities are large, causing overlap of orbits for bodies in different resonances with j > 2. Collisions between planetesimals in different resonances, or between resonant and nonresonant bodies, result in their disruption. Fragments smaller than a critical size can pass through resonances under the influence of drag and be accreted by the embryo. This effect speeds accretion and tends to prevent dynamical isolation of planetary embryos, making gas-rich scenarios for planetary formation more plausible.     

Multiple Periodic Orbits in the Ring Problem: Families of Triple Periodic Orbits

The ring problem deals with the motion of a small body which is subjected to the combined gravitational attraction of N massive bodies arranged in an annular configuration. In this paper we study the distribution of the triple periodic orbits in the phase space of the initial conditions and we discuss their evolution and their principal features. This revised version was published online in July 2006 with corrections to the Cover Date.     

Secular influence of change in the heliocentric gravitation constant <Emphasis Type="Italic">GM</Emphasis><Subscript>⊙</Subscript> on evolution of orbits of Meteor Streams

The Secular influence of the change in the heliocentric gravitational constant on the evolution of orbits of Meteor Streams is examined by using the method of celestial mechanics with variable mass and variable gravitational constant. The change in the heliocentric gravitational constant includes the combined changes in the sun’s mass and gravitational constant obtained from the modern observation of planets and spacecraft. The perturbation equations are solved by expanding series with mean anomaly. The solutions of the secular and periodic variation of orbital elements are derived. The theoretical results for the secular variables of the semi-major axes, solar distances at perihelion and orbital periods are given for three Meteor Streams: Dracorids, Quadrantids, and Ursids. The numerical results are shown in Table 2. The discussion and conclusion are drawn.     

On the Dynamical Foundations of the Lidov–Kozai Theory

The Lidov–Kozai theory developed by each of the authors independently in 1961–1962 is based on qualitative methods of studying the evolution of orbits for the satellite version of the restricted three-body problem (Hill’s problem). At present, this theory is in demand in various fields of science: in the field of planetary research within the Solar system, the field of exoplanetary systems, and the field of high-energy physics in interstellar and intergalactic space. This has prompted me to popularize the ideas that underlie the Lidov–Kozai theory based on the experience of using this theory as an efficient tool for solving various problems related to the study of the secular evolution of the orbits of artificial planetary satellites under the influence of external gravitational perturbations with allowance made for the perturbations due to the polar planetary oblateness.     

Planetary close encounters: geometry of approach and post-encounter orbital parameters

Öpik's assumptions on the geometry of particle trajectories leading to and through planetary close encounters are used to compute the distribution of changes in heliocentric orbital elements that result from such encounters for a range of initial heliocentric orbits. Behaviour at encounter is assumed to follow two-body (particle—planet) gravitational scattering, while before and after encounter particle motion is only governed by the force of the Sun. Derivation of these distributions allows precise analysis of the probability of various outcomes in terms of the physical characteristics of the bodies involved. For example, they allow an explanation and prediction of the asymmetry of the extreme energy perturbations for different initial orbits. The formulae derived here may be applied to problems including the original accumulation of planets and satellites, and the continuing evolution of populations of small bodies, such as asteroids and comets.     

The provenance and evolution of comets

The process of comet formation through the hierarchical aggregation of originally submicron-sized interstellar grains to form micron-sized particles and then larger bodies in the protoplanetary disc, culminating in the formation of planetesimals in the disc extending from Jupiter to beyond Neptune, is briefly reviewed. The ‘planetesimal’ theory for the origin of comets implies the existence of distinct cometary reservoirs, with implications for the immediate provenance of observed comets (both long-period and short-period) and their evolution as a result of planetary perturbations and physical decay, for example splitting and sublimation. The principal mode of cometary decay and collisional interaction with the terrestrial planets is through the formation and evolution of streams of cometary debris and hitherto undiscovered ‘families’ of cometary asteroids. Recent dynamical results, in particular the sungrazing and sun-colliding end-state for short-period comet and asteroid orbits, are briefly discussed.     

Long-Period Comets and the Oort Cloud

Long-period (LP) comets, Halley-type (HT) comets, and even some comets of the Jupiter family, probably come from the Oort cloud, a huge reservoir of icy bodies that surrounds the solar system. Therefore, these comets become important probes to learn about the distant Oort cloud population. We review the fundamental dynamical properties of LP comets, and what is our current understanding of the dynamical mechanisms that bring these bodies from the distant Oort cloud region to the inner planetary region. Most new comets have original reciprocal semimajor axes in the range2 × 10^-5 < 1/a_orig < 5 × 10^-5AU^-1. Yet, this cannot be taken to represent the actual space distribution of Oort cloud comets, but only the region in the energy space in which external perturbers have the greatest efficiency in bringing comets to the inner planetary region. The flux of Oort cloud comets in the outer planetary region is found to be at least several tens times greater than the flux in the inner planetary region. The sharp decrease closer to the Sun is due to the powerful gravitational fields of Jupiter and Saturn that prevent most Oort cloud comets from reaching the Earth’s neighborhood (they act as a dynamical barrier). A small fraction of ～10^-2 Oort cloud comets become Halley type (orbital periods P < 200 yr), and some of them can reach short-period orbits with P < 20 yr. We analyze whether we can distinguish the latter, very ‘old” LP comets, from comets of the Jupier family coming from the Edgeworth-Kuiper belt.