Relationship of comets and planets

We analyze our earlier data on the numerical integration of the equations of motion for 274 short-period comets (with the period P<200 yr) on a time interval of 6000 yr. As many as 54 comets had no close approaches to planets, 13 comets passed through the Saturnian sphere of action, and one comet passed through the Uranian sphere of action. The orbital elements of these 68 comets changed by no more than ±3 percent in a space of 6000 yr. As many as 206 comets passed close to Jupiter. We confirm Everhart’s conclusion that Jupiter can capture long-period comets with q = 4–6 AU and i < 9° into short-period orbits. We show that nearly parabolic comets cross the solar system mainly in the zone of terrestrial planets. No relationship of nearly parabolic comets and terrestrial planets was found for the epoch of the latest apparition of comets. Guliev’s conjecture about two trans-Plutonian planets is based on the illusory excess of cometary nodes at large heliocentric distances. The existence of cometary nodes at the solar system periphery turns out to be a solely geometrical effect.     

Long-period comet flux in the planetary region: Dynamical evolution from the Oort cloud

This study analyzes the evolution of 2 × 10⁵ orbits with initial parameters corresponding to the orbits of comets of the Oort cloud under the action of planetary, galactic, and stellar perturbations over 2 × 10⁹ years. The dynamical evolution of comets of the outer (orbital semimajor axes a > 10⁴ AU) and inner (5 × 10³ < a (AU) < 10⁴) parts of the comet cloud is analyzed separately. The estimates of the flux of “new” and long-period comets for all perihelion distances q in the planetary region are reported. The flux of comets with a > 10⁴ AU in the interval 15 AU < q < 31 AU is several times higher than the flux of comets in the region q < 15 AU. We point out the increased concentration of the perihelia of orbits of comets from the outer cloud, which have passed several times through the planetary system, in the Saturn-Uranus region. The maxima in the distribution of the perihelia of the orbits of comets of the inner Oort cloud are located in the Uranus-Neptune region. “New” comets moving in orbits with a < 2 × 10⁴ AU and arriving at the outside of the planetary system (q > 25 AU) subsequently have a greater number of returns to the region q < 35 AU. The perihelia of the orbits of these comets gradually drift toward the interior of the Solar System and accumulate beyond the orbit of Saturn. The distribution of the perihelia of long-period comets beyond the orbit of Saturn exhibits a peak. We discuss the problem of replenishing the outer Oort cloud by comets from the inner part and their subsequent dynamical evolution. The annual rate of passages of comets of the inner cloud, which replenish the outer cloud, in the region q < 1 AU in orbits with a > 10⁴ AU (～ 5.0 × 10^?14 yr^?1) is one order of magnitude lower than the rate of passage of comets from the outer Oort cloud (～ 9.1 × 10^?13 yr^?1).     

Are There Many Inactive Jupiter-Family Comets among the Near-Earth Asteroid Population?

We analyze the dynamical evolution of Jupiter-family (JF) comets and near-Earth asteroids (NEAs) with aphelion distances Q>3.5 AU, paying special attention to the problem of mixing of both populations, such that inactive comets may be disguised as NEAs. From numerical integrations for 2×10⁶ years we find that the half lifetime (where the lifetime is defined against hyperbolic ejection or collision with the Sun or the planets) of near-Earth JF comets (perihelion distances q<1.3 AU) is about 1.5×10⁵ years but that they spend only a small fraction of this time (∼ a few 10³ years) with q<1.3 AU. From numerical integrations for 5×10⁶ years we find that the half lifetime of NEAs in “cometary” orbits (defined as those with aphelion distances Q>4.5 AU, i.e., that approach or cross Jupiter's orbit) is 4.2×10⁵ years, i.e., about three times longer than that for near-Earth JF comets. We also analyze the problem of decoupling JF comets from Jupiter to produce Encke-type comets. To this end we simulate the dynamical evolution of the sample of observed JF comets with the inclusion of nongravitational forces. While decoupling occurs very seldom when a purely gravitational motion is considered, the action of nongravitational forces (as strong as or greater than those acting on Encke) can produce a few Enckes. Furthermore, a few JF comets are transferred to low-eccentricity orbits entirely within the main asteroid belt (Q<4 AU and q>2 AU). The population of NEAs in cometary orbits is found to be adequately replenished with NEAs of smaller Q's diffusing outward, from which we can set an upper limit of ∼20% for the putative component of deactivated JF comets needed to maintain such a population in steady state. From this analysis, the upper limit for the average time that a JF comet in near-Earth orbit can spend as a dormant, asteroid-looking body can be estimated to be about 40% of the time spent as an active comet. More likely, JF comets in near-Earth orbits will disintegrate once (or shortly after) they end their active phases.     

Capture of comets from the Oort cloud into Halley-type and Jupiter-family orbits

This paper analyzes the capture of comets into Halley-type and Jupiter-family orbits from the nearparabolic flux of the Oort cloud. Two types of capture into Halley-type orbits are found. The first type is the evolution of near-parabolic orbits into short-period orbits (with heliocentric orbital periods P < 200 years) as a result of close encounters with giant planets. This process is followed by a very slow drift of cometary orbits into the inner part of the Solar System. Only those comets may pass from short-period orbits into Halley-type and Jupiter-family orbits, which move in orbits with perihelion distances q < 13 au. In the second type of capture, the perihelion distances of cometary orbits become rather small (< 1.5 au) during the first stage of dynamic evolution under the action of perturbations from the Galaxy, and then their semimajor axes decrease as a result of diffusion. The capture takes place, on average, in 500 revolutions of the comet about the Sun, whereas in the first case, the comet is captured, on average, after 12500 revolutions. The region of initial orbital perihelion distances q > 4 au is found to be at least as important a source of Halley-type comets as the region of perihelion distances q < 4 au. More than half of the Halley-type comets are captured from the nearly parabolic flux with q > 4 au. The analysis of the dynamic evolution of objects moving in short-period orbits shows that the distribution of Centaurs orbits agrees well with the observed distribution corrected for observational selection effects. Hence, the hypothesis associating the origin of Centaurs with the Edgeworth-Kuiper belt and the trans-Neptunian region exclusively should be rejected.     

Vaporization of comet nuclei: Light curves and life times

We have examined the effects of vaporization from the nucleus of a comet and show that a latitude dependence of vaporization can, in some cases, explain asymmetries in cometary light curves. We also find that a non-uniform distribution of solar radiation over a comet can considerably shorten the vaporization lifetime compared to the results normally obtained by assuming that the nuclear surface is isothermal.Independent of any latitude effects, comets with CO₂-dominated nuclei and with perihelion distances less than 0.5 AU have vaporization lifetimes less than or comparable to their dynamical ejection times. This may explain the observed deficit of comets with small perihelion distances. Similarly comets with CO₂-dominated nuclei and perihelia near Jupiter's orbit have vaporization lifetimes that are shorter than the time for capture into short-period orbits. We suggest, therefore, that at least some new comets are composed in large part of CO₂, while only H₂O-dominated comets, with lower vaporization rates, can survive to be captured into short-period orbits.     

Evolution of comet orbits under the perturbing influence of the giant planets and nearby stars

The orbital evolution of 500 hypothetical comets during 10⁹ years is studied numerically. It is assumed that the birthplace of such comets was the region of Uranus and Neptune from where they were deflected into very elongated orbits by perturbations of these planets. Then, we adopted the following initial orbital elements: perihelion distances between 20 and 30 AU, inclinations to the ecliptic plane smaller than 20°, and semimajor axes from 5 × 10³ to 5 × 10⁴ AU. Gravitational perturbations by the four giant planets and by hypothetical stars passing at distances from the Sun smaller than 5 × 10⁵ AU are considered. During the simulation, somewhat more than 50% of the comets were lost from the solar system due to planetary or stellar perturbations. The survivors were removed from the planetary region and left as members of what is generally known as the cometary cloud. At the end of the studied period, the semimajor axes of the surviving comets tend to be concentrated in the interval 2 × 10⁴ < a < 3 × 10⁴ AU. The orbital planes of the comets with initial $\text{a ≧ 3 × 10}{}^4\text{AU}$ acquired a complete randomization while the others still maintain a slight predominance of direct orbits. In addition, comet orbits with final a < 6 × 10⁴AU preserve high eccentricities with an average value greater than 0.8 Most “new” comets from the sample entering the region interior to Jupiter's orbit had already registered earlier passages through the planetary region. By scaling up the rate of paritions of hypothetical new comets with the observed one, the number of members of the cometary cloud is estimated to be about 7 × 10¹⁰ and the conclusion is drawn that Uranus and Neptune had to remove a number of comets ten times greater.     

Invisible comets on evolutionary track of short-period comets

We systematically surveyed the orbits of short-period (SP) comets that show a large change of perihelion distance (q) between 1–2 AU (visible comets) and 4–5 AU (invisible comets) during 4400 years. The data are taken from Cosmo-DICE (Nakamura and Yoshikawa 1991a), which is a long-term orbital evolution project for SP comets. Recognizing that q is the most critical element for observability of comets, an invisibility factor (f), defined as the ratio of unobservable time span to observable span during 4400 years, is calculated for each of the large-q-change comets. A detection limit for each comet is obtained from the heliocentric distance at discovery and/or the absolute magnitude at recent apparitions. A mean f value for 35 SP comets with 2.9 J (J is the Tisserand's invariant) is found to be 19.8. This implies that for each visible SP comet of this J-range, at every epoch of time, there exist about 20 invisible comets near the capture orbits by Jupiter, under the assumptions of steady-state flux and ergodicity for the SP-comet population.     

Close encounters of nearly parabolic comets and planets

An overview is given of close encounters of nearly parabolic comets (NPCs; with periods of P > 200 years and perihelion distances of q > 0.1 AU; the number of the comets is N = 1041) with planets. The minimum distances Δ_min between the cometary and planetary orbits are calculated to select comets whose Δ_min are less than the radius of the planet’s sphere of influence. Close encounters of these comets with planets are identified by numerical integration of the comets’ equations of motion over an interval of ±50 years from the time of passing the perihelion. Close encounters of NPCs with Jupiter in 1663–2011 are reported for seven comets. An encounter with Saturn is reported for comet 2004 F2 (in 2001).     

End-states of short-period comets and their role in maintaining the zodiacal dust cloud

Physical lifetimes and end-states of short-period comets are analysed in connection with the problem of the maintainance of the zodiacal dust cloud. In particular, the problem of the comet-asteroid relationship is addressed. Recent studies of the physical properties of Apollo-Amor asteroids and short-period comets (e.g., Hartmann et al., 1987) show significant differences between them, suggesting that they are distinct classes of objects. A few percent of the active SP comets might become asteroidal-like bodies in comet-type orbits due to the buildup of dust mantles. The remainder probably disintegrate as they consume their volatile content so their debris can only be observed as fireballs when they meet the Earth. Unobservable faint SP comets — i.e., comets so small (m 10¹⁴ g) that quickly disintegrate before being detected, might be a complementary source of dust material. They might be completely sublimated even at rather large heliocentric distances (r - 3 AU). Yet the released dust grains can reach the vicinity of the Sun by Poynting-Robertson drag. The mass associated with unobservable SP comets with perihelion distances q 3 AU might be comparable to that computed for the sample of observed SP co-mets with q 1.5 AU. It is concluded that SP comets (from the large to the unobservable small ones) may supply an average of several tons/sec of meteoric matter to the zodiacal dust cloud.     

On the aging of comets

This study is based primarily on the calculations of comet orbits over ~ 10⁶ years for 160 short-period comets by Harold F. Levison and Martin J. Duncan from which there are calculated “ablation AGES”. There are positive statistical correlations (having many deviations) with radial nongravitational forces, comet activity measures, and dust-to-gas ratios in the spectra, in the sense that comets of greater “AGES” tend to be less active and to show less dust in their spectra than comets of lesser “AGES”.     

Orbital evolution of long-period comets I. Trivariate Monte Carlo simulation

The mean orbital evolution of long-period comets for 16 representative initial orbits to short-period comets is calculated by a Monte Carlo method. First, trivariate perturbation distributions of barycentric Kepler energy, total angular momentum, and its z component in single encounters of comets with Jupiter are obtained numerically. Their characteristics are examined in detail and the distributions are found to be simple, symmetric, and easy to handle. Second, utilizing these distributions, we have done trivariate Monte Carlo simulations of the orbital evolution of long-period comets, with special emphasis on high-inclination orbits. About half of the 16 initial orbits are traced up to 5000 returns. For each of these orbits, the mean values of semimajor axis, perihelion distance, and inclination; their standard deviations, survival, and capture rates; as well as time scales of orbital evolution are calculated as functions of return number. Survival rates of the initial orbits with high inclination (～90°) and small perihelion distance (～1–2 AU) have been found to be only two or three times smaller than those of the main-source orbits of short-period comets established quantitatively by Everhart. The time scales of orbitsl evolution of the former, however, are nearly 10 times longer than the latter. There is a general trend that, for smaller perihelion distance, the survival efficiency becomes higher. The results of this paper should be considered a basis for a succeeding paper (Paper II) in which the physical lifetime of comets will be determined, and a comparison with the orbital data will be done.     

Searching for the source of short-period comet nuclei: The direction of the spatial migration of comets

An attempt is made to determine the spatial location of the main source of short-period comet nuclei. Numerical calculations for the orbital evolution of Jupiter family comets, medium-period comets, and Centaurs are used to show that the orbits of small solar system bodies tend to evolve in the direction of increasing semimajor axes. This relates to bodies that can experience encounters with planets and whose orbital evolution is shaped by gravitational perturbations. It is concluded that there is good reason to search for the main source of the nuclei of Jupiter family comets at distances of 6 AU or less from the sun.     

Orbits of short period comets captured from the Oort cloud

Oort cloud comets occasionally obtain orbits which take them through the planetary region. The perturbations by the planets are likely to change the orbit of the comet. We model this process by using a Monte Carlo method and cross sections for orbital changes, i.e. changes in energy, inclination and perihelion distance, in a single planet-comet encounter. The influence of all major planets is considered. We study the distributions of orbital parameters of observable comets, i.e. those which have perihelion distance smaller than a given value. We find that enough comets are captured from the Oort cloud in order to explain the present populations of short period comets. The median value of cos i for the Jupiter family is 0.985 while it is 0.27 for the Halley types. The results may explain the orbital features of short period comets, assuming that the active lifetime of a comet is not much greater than 400 orbital revolutions.     

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.     

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.     

Near-parabolic cometary flux in the outer solar system

The flux of near-parabolic comets in the outer-planetary region is estimated on the presumption that the major planets and the galactic tide control the dynamics of comets. It is found that the flux of the Oort cloud comets (semi-major axis > 20000 AU) is similar to the case of a strong comet shower derived on the presumption that the galactic tidal force were not operative. On the other hand, the flux of comets with semi-major axes <- 20000 AU is found to be an increasing function of q (perihelion distance) until q reaches 20 AU, while for a 45000 AU it is a rapidly increasing function for q 12 AU. In other words, for comets of the inner extension of the Oort cloud the planetary perturbation acts as a strong barrier for them to penetrate into the inner planetary region.     

Injection of Oort Cloud comets: the fundamental role of stellar perturbations

We present Monte Carlo simulations of the dynamical evolution of the Oort cloud over the age of the Solar System, using an initial sample of one million test comets without any cloning. Our model includes perturbations due to the Galactic tide (radial and vertical) and passing stars. We present the first detailed analysis of the injection mechanism into observable orbits by comparing the complete model with separate models for tidal and stellar perturbations alone. We find that a fundamental role for injecting comets from the region outside the loss cone (perihelion distance q > 15 AU) into observable orbits (q < 5 AU) is played by stellar perturbations. These act in synergy with the tide such that the total injection rate is significantly larger than the sum of the two separate rates. This synergy is as important during comet showers as during quiescent periods and concerns comets with both small and large semi-major axes. We propose different dynamical mechanisms to explain the synergies in the inner and outer parts of the Oort Cloud. We find that the filling of the observable part of the loss cone under normal conditions in the present-day Solar System rises from <1% for a < 20 000 AU to about 100% for a ? 100 000 AU.     

The role of collisions with interplanetary particles in the physical evolution of comets

Effects of collisions with interplanetary particles are investigated. To this purpose, collision probabilities for comets with different orbital elements are computed. It is found that collisions may have a non-negligible effect on the physical evolution of comets. In this connection, it is shown that under certain conditions collisional lifetimes may be shorter than dynamical or vaporization lifetimes. In particular, collisional lifetimes are on average shorter for comets in retrograde orbits than those for direct ones. It is further suggested that catastrophic collisions may contribute to prevent long-period comets in retrograde orbits from reaching short-period orbits by orbital diffusion. Collisions may also produce irregularities of the nucleus brightness by leaving exposed regions of fresh volatile material and may in this way lead to a rejuvenation of old dusty short-period comets. Catastrophic collision probabilities are too low to account for the observed comet splittings, so other trigger mechanisms should be at work. However, it is shown that collisional mini-bursts (increases in brightness of one magnitude or so) caused by decimeter-sized bodies may occur rather frequently on short-period comets when they pass through the asteroid belt. The burst observed in comet Tempel-2 at 3 AU in December, 1978 could be an example of such collisional mini-bursts. The systematic observation of periodic comets when they pass through the asteroid belt could give valuable information about the spatial density of decimeter and meter-sized bodies. In particular, collisional effects for comet Halley, for which a continuous surveillance is planned, are evaluated.     

The frequency and intensity of comet showers from the Oort cloud

We simulate the Oort comet cloud to study the rate and properties of new comets and the intensity and frequency of comet showers. An ensemble of ∼10⁶ comets is perturbed at random times by a population of main sequence stars and white dwarfs that is described by the Bahcall-Soneira Galaxy model. A cloning procedure allows us to model a large ensemble of comets efficiently, without wasting computer time following a large number of low eccentricity orbits. For comets at semimajor axis a = 20,000 AU, about every 100 myr a star with mass in the range 1M_⊙−2M_⊙ passes within ∼10,000 AU of the Sun and triggers a shower that enhances the flux of new comets by more than a factor of 10. The time-integrated flux is dominated by the showers for comets with semimajor axes less than ∼30,000 AU. For semimajor axes greater than ∼30,000 AU the comet loss rate is roughly constant and strong showers do not occur. In some of our simulations, comets are also perturbed by the Galactic tidal field. The inclusion of tidal effects increases the loss rate of comets with semimajor axes between 10,000 and 20,000 AU by about a factor of 4. Thus the Galactic tide, rather than individual stellar perturbations, is the dominant mechanism which drives the evolution of the Oort cloud.     

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.