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.     

Impact cratering and the Oort Cloud

We calculate the expected flux profile of comets into the planetary system from the Oort Cloud arising from Galactic tides and encounters with molecular clouds. We find that both periodic and sporadic bombardment episodes, with amplitudes an order of magnitude above background, occur on characteristic time-scales ∼25–35 Myr. Bombardment episodes occurring preferentially during spiral arm crossings may be responsible both for mass extinctions of life and the transfer of viable microorganisms from the bombarded Earth into the disturbing nebulae. Good agreement is found between the theoretical expectations and the age distribution of large, well-dated terrestrial impact craters of the past 250 Myr. A weak periodicity of ∼36 Myr in the cratering record is consistent with the Sun's recent passage through the Galactic plane, and implies a central plane density  ∼0.15 M_⊙ pc⁻³  . This leaves little room for a significant dark matter component in the disc.     

The Oort Cloud and long‐period comets

This review starts with a brief historical overview of the subject, after which some recent papers attempting to improve the understanding of comet injection from the Oort Cloud and the origin of new comets are discussed. Special attention is paid to the importance of nongravitational effects in comet orbit determination, the synergy between stellar encounters and the galactic tides for the injection dynamics, and the role of planetary perturbations. The field is thus shown to be advancing rapidly, and brief comments on possible implications for studying the origin of the cloud are made.     

Planet X and the Extended Scattered Disk

The effects that a hypothetical trans-Plutonian planet would produce on theorbital distribution of the Classical Edgeworth-Kuiper-Belt, has beensurveyed for different physical and orbital parameters of the hypotheticalbody in Melita et al. (2003a). The best fits were obtained by a moderatelyeccentric and inclined Earth-sized object with a semimajor axis of ～ 70AU. However the history of some objects in the `Extended Scattered disk’still represent a puzzle. One possibility is that they can be `extracted’from the Scattered disk by the planetoid. In this work we confirm that such anhypothesis would not explain the present orbit of 2000 CR105, given theconditions for a gap as observed to be formed in the Classical EKB.     

Embedded star clusters and the formation of the Oort Cloud

Observations suggest most stars originate in clusters embedded in giant molecular clouds [Lada, C.J., Lada, E.A., 2003. Annu. Rev. Astron. Astrophys. 41, 57-115]. Our Solar System likely spent 1-5 Myrs in such regions just after it formed. Thus the Oort Cloud (OC) possibly retains evidence of the Sun's early dynamical history and of the stellar and tidal influence of the cluster. Indeed, the newly found objects (90377) Sedna and 2000 CR₁₀₅ may have been put on their present orbits by such processes [Morbidelli, A., Levison, H.F., 2004. Astron. J. 128, 2564-2576]. Results are presented here of numerical simulations of the orbital evolution of comets subject to the influence of the Sun, Jupiter and Saturn (with their current masses on orbits appropriate to the period before the Late Heavy Bombardment (LHB) [Tsiganis, K., Gomes, R., Morbidelli, A., Levison, H.F., 2005. Nature 435, 459-461]), passing stars and tidal force associated with the gas and stars of an embedded star cluster. The cluster was taken to be a Plummer model with 200-400 stars, with a range of initial central densities. The Sun's orbit was integrated in the cluster potential together with Jupiter and Saturn and the test particles. Stellar encounters were incorporated by directly integrating the effects of stars passing within a sphere centred on the Sun of radius equal to the Plummer radius for low-density clusters and half a Plummer radius for high-density clusters. The gravitational influence of the gas was modeled using the tidal force of the cluster potential. For a given solar orbit, the mean density, 〈ρ〉, was computed by orbit-averaging the density of material encountered. This parameter proved to be a good measure for predicting the properties of the OC. On average 2-18% of our initial sample of comets end up in the OC after 1-3 Myr. A comet is defined to be part of the OC if it is bound and has q>35 AU. Our models show that the median distance of an object in the OC scales approximately as 〈ρ〉^−1/2 when . Our models easily produce objects on orbits like that of (90377) Sedna [Brown, M.E., Trujillo, C., Rabinowitz, D., 2004. Astrophys. J. 617, 645-649] within ∼1 Myr in cases where the mean density is or higher; one needs mean densities of order to create objects like 2000 CR₁₀₅ by this mechanism, which are reasonable (see, e.g., Guthermuth, R.A., Megeath, S.T., Pipher, J.L., Williams, J.P., Allen, L.E., Myers, P.C., Raines, S.N., 2005. Astrophys. J. 632, 397-420). Thus the latter object may also be part of the OC. Close stellar passages can stir the primordial Kuiper Belt to sufficiently high eccentricities (e?0.05; Kenyon, S.J., Bromley, B.C., 2002. Astron. J. 123, 1757-1775) that collisions become destructive. From the simulations performed it is determined that there is a 50% or better chance to stir the primordial Kuiper Belt to eccentricities e?0.05 at 50 AU when . The orbit of the new object 2003 UB₃₁₃ [Brown, M.E., Trujillo, C.A., Rabinowitz, D.L., 2005. Astrophys. J. 635, L97-L100] is only reproduced for mean cluster densities of the order of , but in the simulations it could not come to be on its current orbit by this mechanism without disrupting the formation of bodies in the primordial Kuiper Belt down to 20 AU. It is therefore improbable that the latter object is created by this mechanism.     

Sedna and the Oort Cloud around a migrating Sun

Recent numerical simulations have demonstrated that the Sun’s dynamical history within the Milky Way may be much more complex than that suggested by its current low peculiar velocity (Sellwood, J.A., Binney, J.J. [2002]. Mon. Not. R. Astron. Soc. 336, 785-796; Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79-L82). In particular, the Sun may have radially migrated through the galactic disk by up to 5-6 kpc (Roškar, R., Debattista, V.P., Quinn, T.R., Stinson, G.S., Wadsley, J. [2008]. Astrophys. J. 684, L79-L82). This has important ramifications for the structure of the Oort Cloud, as it means that the Solar System may have experienced tidal and stellar perturbations that were significantly different from its current local galactic environment. To characterize the effects of solar migration within the Milky Way, we use direct numerical simulations to model the formation of an Oort Cloud around stars that end up on solar-type orbits in a galactic-scale simulation of a Milky Way-like disk formation. Surprisingly, our simulations indicate that Sedna’s orbit may belong to the classical Oort Cloud. Contrary to previous understanding, we show that field star encounters play a pivotal role in setting the Oort Cloud’s extreme inner edge, and due to their stochastic nature this inner edge sometimes extends to Sedna’s orbit. The Sun’s galactic migration heightens the chance of powerful stellar passages, and Sedna production occurs around ∼20-30% of the solar-like stars we study. Considering the entire Oort Cloud, we find its median distance depends on the minimum galactocentric distance attained during the Sun’s orbital history. The inner edge also shows a similar dependence but with increased scatter due to the effects of powerful stellar encounters. Both of these Oort Cloud parameters can vary by an order of magnitude and are usually overestimated by an Oort Cloud formation model that assumes a fixed galactic environment. In addition, the amount of material trapped in outer Oort Cloud orbits (a > 20,000 AU) can be extremely low and may present difficulties for traditional models of Oort Cloud formation and long-period comet production.     

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.     

Evidence for an Extended Scattered Disk

By telescopic tracking, we have established that the transneptunian object (TNO) 2000 CR₁₀₅ has a semimajor axis of 220±1 AU and perihelion distance of 44.14±0.02 AU, beyond the domain which has heretofore been associated with the “scattered disk” of Kuiper Belt objects interacting via gravitational encounters with Neptune. We have also firmly established that the TNO 1995 TL₈ has a high perihelion (of 40.08±0.02 AU). These objects, and two other recent discoveries which appear to have perihelia outside 40 AU, have probably been placed on these orbits by a gravitational interaction which is not strong gravitational scattering off of any of the giant planets on their current orbits. Their existence may thus have profound cosmogonic implications for our understanding of the formation of the outer Solar System. We discuss some viable scenarios which could have produced these objects, including long-term diffusive chaos and scattering off of other massive bodies in the outer Solar System. This discovery implies that there must be a large population of TNOs in an “extended scattered disk” with perihelia above the previously suggested 38 AU boundary. The total population is difficult to estimate due to the ease with which such objects would have been lost. This illustrates the great value of frequent and well time-sampled recovery observations of trans-neptunian objects within their discovery opposition.     

Monte Carlo Simulations of Comet Capture from the Oort Cloud

There is a very large number of small bodies in the Solar System and their orbits are varied and complicated. Some types of orbits and events are so rare that they occur in numerical simulations only when millions or billions of orbits have been calculated. In order to study these orbits or events an efficient Monte Carlo simulation is useful. Here we describe a new Monte Carlo simulation method and test it against some existing simulations of orbits of small bodies which have been obtained by different methods. We find good agreement with many earlier calculations, and study briefly the possibility of the Oort Cloud capture origin of short period comets. This revised version was published online in July 2006 with corrections to the Cover Date.     

The Role of Giant Molecular Clouds in the Evolution of the Oort Comet Cloud

We estimated the gravitational influence of giant molecular clouds passing near the Solar system on the orbital evolution of Oort cloud comets. We performed a comparative analysis of the accuracies of the following two methods of allowance for the perturbations from giant molecular clouds: the impulse approximation and numerical integration. The impulse approximation yields fairly accurate estimates of the change in the energy of Oort cloud comets and the probability of their ejection under the influence of a molecular cloud if the path of the Solar system does not cross its boundary and if the molecular cloud may be treated as a point perturbing mass. The comet survival probability in the Oort cloud depends significantly on the internal structure of the perturbing molecular cloud and the impact parameter of the encounter. The most massive injection of comets into the planetary region and their ejection from the Oort cloud take place if the Solar system passes through a giant molecular cloud composed of several high-mass condensations. In this case, most of the comets injected into the planetary region were initially comets of the inner Oort cloud (a 10^–4 AU) with high orbital eccentricities.     

Long-term evolution of Oort Cloud comets: methods and comparisons

Two long-term simulation methods for cometary orbits, a Monte Carlo method and a direct integration method, are compared with each other. The comparison is done in seven inclination and perihelion distance intervals, and shows differences in dynamical lifetime and capture probabilities for the following main reasons. We use a finite energy step approximation in the Monte Carlo method and the method considers only close approaches with the planets. The differences can be taken into account statistically and it is possible to calculate the correction factors for the capture probability and dynamical lifetime in the Monte Carlo method. Both corrections depend on the inclination and on the value of the minimum energy step. The capture probabilities of the short-period comets originating in the Oort Cloud are calculated by the corrected Monte Carlo method and compared with published results.     

Algorithms for Stellar Perturbation Computations on Oort Cloud Comets

We investigate different approximate methods of computing the perturbations on the orbits of Oort cloud comets caused by passing stars, by checking them against an accurate numerical integration using Everhart’s RA15 code. The scenario under study is the one relevant for long-term simulations of the cloud’s response to a predefined set of stellar passages. Our sample of stellar encounters simulates those experienced by the Solar System currently, but extrapolated over a time of 10¹⁰ years. We measure the errors of perihelion distance perturbations for high-eccentricity orbits introduced by several estimators – including the classical impulse approximation and Dybczyński’s (1994, Celest. Mech. Dynam. Astron. 58, 1330–1338) method – and we study how they depend on the encounter parameters (approach distance and relative velocity). We introduce a sequential variant of Dybczyński’s approach, cutting the encounter into several steps whereby the heliocentric motion of the comet is taken into account. For the scenario at hand this is found to offer an efficient means to obtain accurate results for practically any domain of the parameter space.     

Population of the Scattered Kuiper Belt

We study how the internal structure of dark halos is affected if cold dark matter particles are assumed to have a large cross section for elastic collisions. We identify a cluster halo in a large cosmological N-body simulation and resimulate its formation with progressively increasing resolution. We compare the structure found in the two cases in which dark matter is treated as collisionless or as a fluid. For the collisionless case, the overall ellipticity of the cluster, the central density cusp, and the amount of surviving substructure are all similar to those found in earlier high-resolution simulations. Collisional dark matter results in a cluster that is more nearly spherical at all radii, has a steeper central density cusp, and has less-but still substantial-surviving substructure. As in the collisionless case, these results for a "fluid" cluster halo are expected to carry over approximately to smaller mass systems. The observed rotation curves of dwarf galaxies then argue that self-interacting dark matter can only be viable if intermediate cross sections produce structure that does not lie between the extremes we have simulated.     

Long-term evolution of Oort Cloud comets: capture of comets

We test different possibilities for the origin of short-period comets captured from the Oort Cloud. We use an efficient Monte Carlo simulation method that takes into account non-gravitational forces, Galactic perturbations, observational selection effects, physical evolution and tidal splittings of comets. We confirm previous results and conclude that the Jupiter family comets cannot originate in the spherically distributed Oort Cloud, since there is no physically possible model of how these comets can be captured from the Oort Cloud flux and produce the observed inclination and Tisserand constant distributions. The extended model of the Oort Cloud predicted by the planetesimal theory consisting of a non-randomly distributed inner core and a classical Oort Cloud also cannot explain the observed distributions of Jupiter family comets. The number of comets captured from the outer region of the Solar system are too high compared with the observations if the inclination distribution of Jupiter family comets is matched with the observed distribution. It is very likely that the Halley-type comets are captured mainly from the classical Oort Cloud, since the distributions in inclination and Tisserand value can be fitted to the observed distributions with very high confidence. Also the expected number of comets is in agreement with the observations when physical evolution of the comets is included. However, the solution is not unique, and other more complicated models can also explain the observed properties of Halley-type comets. The existence of Jupiter family comets can be explained only if they are captured from the extended disc of comets with semimajor axes of the comets   a <5000 au  . The original flattened distribution of comets is conserved as the cometary orbits evolve from the outer Solar system era to the observed region.     

Coupling dynamical and collisional evolution of small bodies: II. Forming the Kuiper belt, the Scattered Disk and the Oort Cloud

The Oort Cloud, the Kuiper belt and the Scattered Disk are dynamically distinct populations of small bodies evolving in the outer regions of the Solar System. Whereas their collisional activity is now quiet, gravitational interactions with giant planets may have shaped these populations both dynamically and collisionally during their formation. Using a hybrid approach [Charnoz, S., Morbidelli, A., 2003. Icarus 166, 141-166], the present paper tries to couple the primordial collisional and dynamical evolution of these three populations in a self-consistent way. A critical parameter is the primordial size-distribution. We show that the initial planetesimal size distribution that allows an effective mass depletion of the Kuiper belt by collisional grinding, would decimate also the population of comet-size bodies that end in the Oort Cloud and, in particular, in the Scattered Disk. As a consequence, the Oort Cloud and the Scattered Disk would be too anemic, by a factor 20 to 100, relative to the estimates achieved from the observation of the fluxes of long period and Jupiter family comets, respectively. For these two reservoirs to have a sufficient number of comets, the initial size distribution in the planetesimal disk had to be such that the mass depletion by collisional erosion of the Kuiper belt was negligible. Consequently the current mass deficit of the Kuiper belt needs to be explained by dynamical mechanisms.     

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.     

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.     

Few Comments on the Relation Between the Initial Proto-planetary Disc Model and the Oort Cloud Formation

A fraction of small bodies from the once existing proto-planetary disc was ejected, by the giant planets, to large heliocentric distances and start to build the comet Oort cloud. Considering four models of initial proto-planetary disc, we attempt to roughly map a dependence between the initial disc’s structure and some properties of the Oort cloud. We find that it is difficult to construct the proto-planetary disc if (i) the amount of heavy chemical elements in Jupiter and Saturn is as high as currently accepted and (ii) the total mass of the minimum-mass solar nebula is assumed to be lower than $\approx0.05\,\hbox{M}_{\odot}.$ The behaviour of the Oort cloud formation does not crucially depend on the initial disc model. Some differences in its structure are obvious: since the cloud is known to be filled mainly by Uranus and Neptune, the efficiency of its formation is higher when the initial amount of particles in the Uranus-Neptune region is relatively higher. A significantly large number of Jupiter Trojans in our simulation appears, however, only in the case of the initially non-gapped disc, with the particles situated also close to the Jupiter’s orbit.