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).     

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.     

Perihelion evolution of observed new comets implies the dominance of the galactic tide in making Oort cloud comets discernable

For an Oort cloud comet to be seen as a new comet, its perihelion must be moved from a point exterior to the loss cylinder boundary to a point interior to observable limits in a single orbit. The galactic tide can do this continuously, in a non-impulsive manner. Near-parabolic comets, with specific angular momentum , will most easily be made observable. Therefore, to reduce the perihelion distance H must decrease. Since weakly perturbed comets are, in general, more numerous than strongly perturbed comets, we can anticipate that new comets made observable by a weak tidal torque will more likely be first observed when their slowly changing perihelion distances are approaching their minimum osculating values under the action of the tide, rather than receding from their minimum values. That is, defining ΔH_tide as the vector change due to the galactic tidal torque during the prior orbit, and H_obs as the observed vector, the sign S≡Sign(H_obs·ΔH_tide) will more likely be −1 than +1 if a weak galactic tidal perturbation indeed dominates in making comets observable. Using comet data of the highest quality class (1A) for new comets (a>10,000 AU), we find that 49 comets have S=−1 and 22 have S=+1. The binomial probability that as many or more would exhibit this characteristic if in fact S=?1 were equally likely is only 0.0009. This characteristic also persists in other long-period comet populations, lending support to the notion that they are dominated by comets recently arrived from the outer Oort cloud. The preponderance of S=−1 also correlates with weakly perturbed (i.e., smaller semimajor axis) new comets in a statistically significant manner. This is strong evidence that the data are of sufficiently high quality and sufficiently free of observational selection effects to detect this unique imprint of the tide.     

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.     

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.     

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.     

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.     

Galactic perturbations on nearly-parabolic cometary orbits

The effect of galactic perturbations on long-period comet orbits is examined via numerical and analytical means. Relations are found between a comet's initial perihelion position and the positions of succeeding perihelia. It was found that the galactic effects were strongest on the comets initially at galactic latitudes close to 40°. In such cases the galactic perturbations caused the orbit to become almost circular before becoming nearly parabolic again. This effect allows comets with semimajor axes of about 25 000 AU to make only a few passages through the inner solar system in a time interval of 10⁹yr. Thus the galactic field is an important factor in the evolution of long-period comet orbits. The observed distribution of perihelia of long-period comets indicates that galactic effects have been active.     

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.     

A model of the galactic tidal interaction with the oort comet cloud

We consider a model of the in situ Oort cloud which is isotropic with a random distrihution of perihelia directions and angular momenta. The energy distribution adopted has a continuous range of values appropriate for long-period (>200 yr) comets. Only the tidal torque of the Galaxy is included as a perturbation of comet orbits and it is approximated to be that due to a quasi-steady state distribution of matter with disk-like symmetry. The time evolution of all orbital elements can be analytically obtained for this case. In particular, the change in the perihelion distance per orbit and its dependence on other orbital elements is readily found. We further make the assumption that a comet whose perihelion distance was beyond 15 AU during its last passage through the Solar System would have orbit parameters that are essentially unchanged by planetary perturbations. Conversely, if the prior passage was inside 15 AU we assume that planetary perturbations would have removed the comet from the in situ energy distribution accessible by the galactic tide. Comets which had their perihelia changed from beyond 15 AU to within 5 AU in a single orbit are taken to be observable. We are able to track the evolution of 10⁶ comets as they are made observable by the galactic tidal touque. Detailed results are obtained for the predicted distribution of new (0 < 1/ < 10^–4 AU^–1) comets. Further, correlations between orbital elements can be studied. We present predictions of observed distributions and compare them with the random in situ results as well as with the actual observed distributions of class I comets. The predictions are in reasonable agreement with actual observations and, in many cases, are significantly different from random when perihelia directions are separated into galactic northern and southern hemispheres. However the well-known asymmetry in the north-south populations of perihelia remains to be explained. Such an asymmetry is consistent with the dominance of tidal torques today if a major stochastic event produced it in the past since tidal torques are unable to cause the migration of perihelia across the latitude barriers ±26°.6 in the disk model.     

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.     

A probability of encounter with interstellar comets and the likelihood of their existence

A theory of the probability of encounter of the Sun with an interstellar comet at a distance comparable to the Earth-Sun distance is formulated, and a general expression is derived establishing the relationship among the influx rate of interstellar comets, the perihelion distance, the space density of the comets, the Maxwellian distribution of comet velocities in the interstellar cloud, and the cloud's systematic velocity relative to the Sun. The fact that no comet with a strongly hyperbolic orbit has so far been observed is used to determine an upper limit of 6 × 10^?4 solar masses per cubic parsec (4 × 10^?26 gcm^?3) for the space density of interstellar comets. The theoretical distribution of semimajor axes of interstellar comets is derived to show that a strong hyperbolic excess must be present in the orbits of a majority of interstellar comets regardless of the dynamical characteristics of the comet cloud, except when the cloud is moving along with the Sun and the distribution of individual velocities has a very low dispersion. This case, however, implies a possibility of capture by the Sun and thus becomes a problem of an Oort-type cloud.     

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.     

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).     

From the Oort cloud to observable short-period comets – I. The initial stage of cometary capture

We investigate the first stage of the dynamical evolution of Oort cloud comets entering the planetary region for the first time. To this purpose, we integrate numerically the motions of a large number of fictitious comets pertaining to two samples, both with perihelion distances up to 5.7 au and random inclinations; the first sample is composed of comets whose orbits have at least one node close to 5.2 au, while the second is not subject to this constraint. We examine the orbits when the comets come to aphelion after their first perihelion passage within the planetary region, and find that there is a clear statistical dependence of the energy perturbations on the Tisserand parameter. There appear to be two main processes, of comparable importance, governing the shortening of semimajor axes to values of less than 1000 au, i.e. planetary close encounters, especially with Jupiter, and indirect perturbations due to the shifting of the motion from barycentric to heliocentric and back; the former process mostly affects comets crossing the ecliptic at about 5.2 au, or on low-inclination orbits, while the latter mostly affects comets of small perihelion distance. This last result may help to understand the relative paucity of Halley-type comets with perihelion distances larger than about 1.5 au.     

Tidal gravitational forces: The infall of “new” comets and comet showers

The tidal gravitational field of the Galaxy directed into the galactic plane changes the angular momentum of comets in the Oort cloud. For comet orbits with semimajor axis greater than 2 × 10⁴ AU, the change of angular momentum in one orbit is sufficient to bring comets from the Oort cloud into the visible region, causing the infall of “new” comets. The limiting size orbit is weakly dependent on the angle between the major axis of the comet orbit and the galactic plane. The flux of comets into the inner Solar System caused by the galactic tidal field will be continuous and nearly isotropic. This effect appears to exclude any determination of the trajectories of passing stars by analysis of the angular distribution of new comets. The production of intense comet showers by the tidal field of a solar companion or of an interstellar cloud is considered. We show that the direction of a solar companion cannot be found from the present distribution of observable comets. The frequency of comet showers induced by encounters with interstellar clouds is found to be much lower than that from passing stars, and the tidal fields of interstellar clouds are not strong enough to cause comet showers of sufficient intensity to result in Earth impacts.     

Orientation-dependent effects in Oort's theory of comet origin II. Anisotropies in the distribution of long-period comet orbits

The statistical significance of anisotropies in the distribution of orbital orientations among the long-period and nearly parabolic comets is evaluated. It is suggested that these anisotropies are not the result of observational selection effects. A numerical model for the distribution of orbital orientations is constructed, based on Oort's theory of comet origin and the assumption that the observed anisotropies are caused by multiple planetary perturbations over the course of many perihelion passages. The model, which is restricted to comets with peri-helion distances less than 0.3AU, does not predict any significant anisotropies.     

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.     

Characteristics and Frequency of Weak Stellar Impulses of the Oort Cloud

We have developed a model of the response of the outer Oort cloud of comets to simultaneous tidal perturbations of the adiabatic galactic force and a stellar impulse. The six-dimensional phase space of near-parabolic comet orbital elements has been subdivided into cells. A mapping of the evolution of these elements from beyond the loss cylinder boundary into the inner planetary region over the course of a single orbit is possible. This is done by treating each perturbation separately, and in combination, during a time interval of 5 Myr. We then obtain the time dependence of a wide range of observable comet flux characteristics, which provides a fingerprint of the dynamics. These include the flux distributions of energy, perihelion distance, major axis orientation, and angular momentum orientation. Correlations between these variables are also determined. We show that substantive errors occur if one superposes the separately obtained flux results of the galactic tide and the stellar impulse rather than superposing the tidal and impulsive perturbations in a single analysis. Detailed illustrations are given for an example case where the stellar mass and relative velocity have the ratio M∗/V_rel=0.043 M⊙/km s⁻¹ and the solar impact parameter is 45,000 AU. This case has features similar to the impending Gliese 710 impulse with the impact parameter selected to be close to the low end of the predicted range. We find that the peak in the observable comet flux exceeds that due to the galactic tide alone by ≈41%. We also present results for the time dependence of the flux enhancements and for the mean encounter frequency of weak stellar impulse events as functions of M∗/V_rel and solar impact parameter.     

Cometary impacts with the Sun: Physical and dynamical considerations

D. J. Michels, N. R. Sheeley, Jr., R. A. Howard, and M. J. Koomen (Science215, 1097–1102, 1982) observed a comet which appears to have impacted the Sun. Z. Sekanina (Astron. J..87, 1059–1072, 1982) showed that the comet, 1979XI, was probably a member of the Kreutz group of sungrazing comets. The sungrazers typically have perihelia of 1.2–1.9 solar radii but Sekanina found q = 0.35 R_⊙ for 1979XI. It is interesting to speculate how the perihelion may have been reduced to this small value. The change in perihelion can not be explained by planetary, stellar, or nongravitational perturbations. Tidal splitting of the nucleus on a previous perihelion passage is also ruled out, through a random splitting event near aphelion of the comet's orbit is a remote possibility. The most plausible explanation is collision with another body, most likely a comet, at large heliocentric distance. However, the expected probability of such an event is exceedingly small. Another aspect of the problem is whether the nucleus of 1979XI sublimated completely before impacting the Sun. Assuming a water ice nucleus, it is shown that a surface layer of only 5–15 m thickness would be sublimated prior to impact. Although it is likely that the nucleus tidally disrupted after crossing the solar Roche limit, the ultimate destruction of the nucleus probably resulted from the shock of hitting the denser regions of the solar atmosphere, just above the photosphere.