Constraints on the Formation Regions of Comets from their D:H Ratios

Studies of the D:H ratio in H₂O within the Solar nebula provide a relationship between the degree of enrichment of deuterium and the distance from the young Sun. In the context of cometary formation, such models suggest that comets which formed in different regions of the Solar nebula should have measurably different D:H ratios. We aim to illustrate how the observed comets can give information about the formation regions of the reservoirs in which they originated. After a discussion of the current understanding of the regions in which comets formed, simple models of plausible formation regions for two different cometary reservoirs (the Edgeworth–Kuiper belt and the Oort Cloud) are convolved with a deuterium-enrichment profile for the pre-solar nebula. This allows us to illustrate how different formation regions for these objects can lead to great variations in the deuterium enrichment distributions that we would observe in comets today. We also provide an illustrative example of how variations in the population within a source region can modify the resulting observational profile. The convolution of a deuterium-enrichment profile with examples of proto-cometary populations gives a feel for how observations could be used to draw conclusions on the formation region of comets which are currently fed into the inner Solar system from at least two reservoirs. Such observations have, to date, been carried out on only three comets, but future work with instruments such as ALMA and Herschel should vastly improve the dataset, leading to a clearer consensus on the formation of the Oort cloud and Edgeworth–Kuiper belt.     

New constraints on the delivery of cometary water and nitrogen to Earth from the N/N isotopic ratio

New independent constraints on the amount of water delivered to Earth by comets are derived using the ¹⁵N/¹⁴N isotopic ratio, measured to be roughly twice as high in cometary CN and HCN as in the present Earth. Under reasonable assumptions, we find that no more than a few percent of Earth’s water can be attributed to comets, in agreement with the constraints derived from D/H. Our results also suggest that a significant part of Earth’s atmospheric nitrogen might come from comets. Since the ¹⁵N/¹⁴N isotopic ratio is not different in Oort-cloud and Kuiper-belt comets, our estimates apply to the contribution of both types of objects.     

A survey of debris trails from short-period comets

We observed 34 comets using the 24 μm camera on the Spitzer Space Telescope. Each image contains the nucleus and covers at least ¹⁰⁶ km of each comet's orbit. Debris trails due to mm-sized or larger particles were found along the orbits of 27 comets; 4 comets had small-particle dust tails and a viewing geometry that made debris trails impossible to distinguish; and only 3 had no debris trail despite favorable observing conditions. There are now 30 Jupiter-family comets with known debris trails, of which 22 are reported in this paper for the first time. The detection rate is >80%, indicating that debris trails are a generic feature of short-period comets. By comparison to orbital calculations for particles of a range of sizes ejected over 2 yr prior to observation, we find that particles comprising 4 debris trails are typically mm-sized while the remainder of the debris trails require particles larger than this. The lower-limit masses of the debris trails are typically ¹⁰¹¹ g, and the median mass loss rate is 2 kg/s. The mass-loss rate in trail particles is comparable to that inferred from OH production rates and larger than that inferred from visible-light scattering in comae.     

The activities of comets related to their aging and origin

Two indices have been developed for the purpose of comparing the natures of various classes of comets. The first is the Activity Index (AI), measuring the inherent magnitude increase in brightness from great solar distances to maximum near perihelion. The second, or Volatility Index (VI), measures the variation in magnitude near perihelion. Tentative determinations of these two indices are derived from observations by Max Beyer over more than 30 years for long-period (L-P) and short-period (S-P) comets near perihelion and from other homogeneous sources. AI determinations are made for 32 long-period (L-P) comets and for 14 short-period (S-P). The range of values of AI is of the order of 3 to 10 magnitudes with a median about 6. An expected strong correlation with perihelion distance q, is found to vary as q ^–2.3. Residuals from a least-square solution (AI) are used for comparing comets of different orbital classes, the standard deviation of a single value of AI is only ±1^m.1 for L-P comets and ±1^m.2 for S-P comets.Among the L-P comets, 19 of period P larger than 10⁴ years yield AI = 0^m.27 ± 0^m.25 compared to 0^m.39 ± 0^m.26 for 13 of period between 10² years and 10⁴ years. This denies any fading with aging among the L-P comets. Also no systematic change with period occurs for the VI index, leading to the same conclusions. Weak correlations are found with the Gas/Dust ratio of comets. No correlations are found between the two indices, nor of either index with near-perihelion magnitudes or orbital inclination.The various data are consistent with a uniform origin for all types of comets, the nuclei being homogeneous on the large scale but quite diverse on a small scale (the order of a fraction of kilometer in extent). Small comets thus may sublimate away entirely, leaving no solid core, while huge comets may develop a less volatile core by radioactive heating and possibly become inactive like asteroids after many S-P revolutions about the Sun. When relatively new, huge comets may be quite active at great solar distances because of volatiles from the core that have refrozen in the outer layers.     

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

On the hypothesis of the eruption origin of near-parabolic comets

The hypothesis on the genetic connection of near-parabolic comets with Jupiter, Saturn, and the transPlutonian region (5–3000 AU) proposed by E.M. Drobyshevskii is considered. It has been shown that, on average, 5.6 comets per an area of 10⁶ AU² passed through the transPlutonian region during the whole history of observations. Six-hundred nineteen comets crossed the ecliptic at heliocentric distances ranging from 0 to 2 AU. As has been shown, from the total number of 945 near-parabolic comets, eight comets closely approached Jupiter and five closely approached Saturn. The Kreutz comets, 1277 objects, did not approach Jupiter closer than 3 AU. Their minimal distance to Saturn was 5.5 AU. The minimal distance of the Kreutz comets from the edge of the transPlutonian region was 28.8 AU. The analysis led to the conclusion that the concept on the origin of the near-parabolic comets suggested by Drobyshevskii is groundless.     

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.     

A note on the total mass of comets in the solar system

The assumption that the very low albedo determined for Halley's comet is typical of all short period comets, taken together with the assumption that the average sizes of long and short period comets are approximately equal, leads to an increase in the total mass of comets in the solar system by almost two orders of magnitude. If gravitational ejection from the Uranus - Neptune zone during the later phases of planet formation is indeed responsible for the classical Oort cloud between 10⁴–10¹⁵ AU, then the mass of comets in this transplanetary region during cosmogonie times has to exceed the combined masses of Uranus and Neptune by over an order of magnitude. Furthermore, if the recent arguments for as many as 10¹⁴ comets in an inner Oort cloud between ~40– 10⁴AU are valid, then the total mass of comets in the solar system approaches 2% of a solar mass.     

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.     

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.     

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

The lifetime of the OH radical in comets at 1 AU

It has been shown that the photochemical lifetime of OH in comets is a function of the comet's radial velocity. The calculated lifetime at 1 AU can vary between 6.9 × 10⁴ to 2.1 × 10⁵ sec for radial velocities that vary from ?58 to +59 km/sec. A comparison between the scale lengths observed for three comets and those calculated based upon the theoretical lifetime has been made. This comparison shows that in two of the comets the lifetime derived from the scale lengths is a factor of 1.7 larger than the theoretical lifetime. Suggestions are made about the origin of this discrepancy.     

Comet Grains: Their IR Emission and Their Relation to ISm Grains

Comets and the chondritic porous interplanetary dust particles (CP IDPs) that they shed in their comae are reservoirs of primitive solar nebula materials. The high porosity and fragility of cometary grains and CP IDPs, and anomalously high deuterium contents of highly fragile, pyroxene-rich Cluster IDPs imply these aggregate particles contain significant abundances of grains from the interstellar medium (ISM). IR spectra of comets (3–40 μm) reveal the presence of a warm (near-IR) featureless emission modeled by amorphous carbon grains. Broad andnarrow resonances near 10 and 20 microns are modeled by warm chondritic (50% Feand 50% Mg) amorphous silicates and cooler Mg-rich crystalline silicate minerals, respectively. Cometary amorphous silicates resonances are well matched by IRspectra of CP IDPs dominated by GEMS (0.1 μm silicate spherules) that are thought to be the interstellar Fe-bearing amorphous silicates produced in AGB stars. Acid-etched ultramicrotomed CP IDP samples, however, show that both the carbon phase (amorphous and aliphatic) and the Mg-rich amorphous silicate phase in GEMS are not optically absorbing. Rather, it is Fe and FeS nanoparticles embedded in the GEMS that makes the CP IDPs dark. Therefore, CP IDPs suggest significant processing has occurred in the ISM. ISM processing probably includes in He⁺ ion bombardment in supernovae shocks. Laboratory experiments show He⁺ ion bombardment amorphizes crystalline silicates, increases porosity, and reduces Fe into nanoparticles. Cometary crystalline silicate resonances are well matched by IR spectra of laboratory submicron Mg-rich olivine crystals and pyroxene crystals. Discovery of a Mg-pure olivine crystal in a Cluster IDP with isotopically anomalous oxygen indicates that a small fraction of crystalline silicates may have survived their journey from AGB stars through the ISM to the early solar nebula. The ISM does not have enough crystalline silicates (<5%), however, to account for the deduced abundance of crystalline silicates in comet dust. An insufficient source of ISMMg-rich crystals leads to the inference that most Mg-rich crystals in comets are primitive grains processed in the early solar nebula prior to their incorporation into comets. Mg-rich crystals may condense in the hot (～1450 K), inner zones of the early solar nebula and then travel large radial distances out to the comet-forming zone. On the other hand, Mg-rich silicate crystals may be ISM amorphous silicates annealed at ～1000 K and radially distributed out to the comet-forming zone or annealed in nebular shocks at ～5-10 AU. Determining the relative abundance of amorphous and crystalline silicatesin comets probes the relative contributions of ISM grains and primitive grains to small, icy bodies in the solar system. The life cycle of dust from its stardust origins through the ISM to its incorporation into comets is discussed.     

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.     

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.     

The origin of the CN radical in comets: A review from observations and models

The origin of CN radicals in comets is not completely understood so far. We present a study of CN and HCN production rates and CN Haser scale lengths showing that: (1) at heliocentric distances larger than 3 AU, CN radicals could be entirely produced by HCN photolysis; (2) closer to the Sun, for a fraction of comets CN production rates are higher than HCN ones whereas (3) in the others, CN distribution cannot be explained by the HCN photolysis although CN and HCN production rates seem to be similar. Thus, when the comets are closer than 3 AU to the Sun, an additional process to the HCN photolysis seems to be required to explain the CN density in some comets.The photolysis of HC₃N or C₂N₂ could explain the CN origin. But the HC₃N production rate is probably too low to reproduce CN density profile, even if uncertainties on its photolysis leave the place for all possible conclusions. The presence of C₂N₂ in comets is a reliable hypothesis to explain the CN origin; thus, its detection is a challenging issue. Since C₂N₂ is very difficult to detect from ground-based observations, only in situ measurements or space observations could determine the contribution of this compound in the CN origin.Another hypothesis is a direct production of CN radicals by the photo- or thermal degradation of complex refractory organic compounds present on cometary grains. This process could explain the spatial profile of CN inside jets and the discrepancy noted in the isotopic ratio ¹⁴N/¹⁵N between CN and HCN. Laboratory studies of the thermal and UV-induced degradation of solid nitrogenated compounds are required to model and validate this hypothesis.     

Transneptunian object 2003 UB 313 as a source of comets

The possibility of interrelation between long-period comets and 2003 UB 313, a recently discovered large Kuiper Belt body, is investigated. For this purpose, 78 objects crossing the plane of motion of this body at distances from 37.8 to 97.6 AU have been selected from 860 long-period comets. The overpopulation of comets with this characteristic is also considered. The plane of motion of 2003 UB 313 is compared with the orbital planes of other objects in number of comet crossings in the specified distance interval or in some parts of it. A statistically significant overpopulation of elliptic and intermediate comets with the corresponding orbital nodes has been established. Recently discovered and absolutely faint comets show the best effect in this sense. The same is also true for comets with osculating eccentricities e < 1. A similar result is also obtained for comets with “original” a ^?1 > 0.010000. It is hypothesized that the 2003 UB 313 family is present among the 78 comets. Four of them have aphelion distances from 37.8 to 97.6 AU. An ellipticity is traceable in the distribution of some of the 78 distant nodes. This may be considered as a further argument for the suggested hypothesis. Generally, the body 2003 UB 313 may be assumed to play a prominent role in injecting observable comets from the transneptunian region     

Capture of comets during the evolution of a star cluster and the origin of the Oort Cloud

It is generally assumed that the Solar System is surrounded by a swarm of comets, the so-called Oort Cloud, which contains approximately 10¹¹ members. The observed comets belong to a small subsection of the Cloud, and they have very elongated orbits. The origin of the Cloud is presently unclear. Here we consider the possibility that the comets were born in a star cluster together with the Sun. We follow the evolution of the star cluster with its embedded swarm of comets and calculate the rate at which stars accumulate stable comet companions. We conclude that if the Oort Cloud of comets was born in this process, then the present day density of comets in interstellar space has to be high, and that comets make a significant contribution to the overall mass density of the Galaxy.     

Dynamical evolution of the Oort cloud: I. A Monte Carlo simulation

We have studied the long-time dynamical evolution of a population of comets surrounding the Solar System at a large distance. Orbital changes are caused by random passing stars. We first emphasize the need for a new simulation because of the lack of completeness of previous analytical and numerical studies. Then the solar neighborhood is modeled by a sphere of 1 pc in radius, which stars cross at random in direction and distance. The geometry of the encounters allows us to compute the impulse gained by the star and the Sun, in the context of an impact approximation. Then we determine the change of orbital elements for a population of comets and follow the evolution of the frequency distribution for the five Keplerian elements. Clouds are selected in such a way that we test the two main hypotheses for the origin of the Oort cloud, and also the regions of stability in an aphelion-eccentricity diagram. We show that stellar perturbations randomize the cloud and prevent one from inferring the initial cloud configuration from the current distribution. Clouds are depleted by the diffusion of comets into the planetary regions, where they become planet-influenced comets or are ejected from the Solar System. The diffusion of aphelion toward interstellar regions proves to be the major source of cometary loss. Direct ejection to hyperbolic orbits amounts to 9% of the originally population over the age of the Solar System. Finally the current and original cloud populations are estimated at 1.8 × 10¹² and 2 × 10¹³ comets and we discuss these results.     

Comparison of X-rays from Comets LINEAR (C/1999 S4) and McNaught-Hartley (C/1999 T1)

The Chandra X-ray Observatory (CXO) observations of Comets McNaught-Hartley (MH) and LINEAR S4 (S4) have been processed in the same way to compare X-rays from those comets. The X-ray isophotes are crescent-like in S4 and more circular in MH because of the different phase angles (98° and 44°, respectively). The peak X-ray brightness is greater in S4 than that in MH by a factor of 1.5 and smaller by a factor of 1.7 after the correction for heliocentric distance. The X-ray luminosities of MH and S4 are equal to 8.6 and 1.4×10¹⁵ erg s⁻¹ inside the apertures of ρ=1.5 and 0.5×10⁴ km, respectively. (Brightness is 20% of the peak value at these ρ.) Efficiencies of X-ray excitation corrected to the solar wind flow are similar and equal to 4.3×10⁻¹⁴ erg AU^3/2 in both comets. This confirms the solar wind excitation of X-rays in comets. Spectra of the comets were extracted with a special care of the background correction and using an energy-dependent spectral resolution code. The MH spectrum consists of ten emissions instead of nine emissions in the previously published spectrum. The new emission at 307 eV fills in a strong minimum in the previous spectrum and removes the major difference between that spectrum and the synthetic spectrum. This emission is assigned to the C⁺⁴ and Mg⁺⁹ lines. The positions of the other emissions and their identification are similar to those in the previous spectrum. The S4 spectrum consists of eight emissions, and four emissions are the same as in MH. The line identification is given. Ion ratios in the solar wind have been extracted from the spectra. O⁺⁸/O⁺⁷ is equal to 0.29±0.04 and 0.14±0.02 in MH and S4, and this difference correlates with the higher solar wind speed in S4. Ne⁺⁹/O⁺⁷ is (15±6)×10⁻³ and (19±7)×10⁻³, and these are the first data on Ne⁺⁹ in the solar wind. C⁺⁶/O⁺⁷ is 0.7±0.2 in both MH and S4. X-ray spectroscopy of comets may be used as a diagnostic tool to study the solar wind composition.