Nuclear magnitudes and the size distribution of Jupiter family comets

We present a new catalog of absolute nuclear magnitudes of Jupiter family (JF) comets, which is an updated version of our previous catalog [Tancredi, G., Fernández, J.A., Rickman, H., Licandro, J., 2000. Astron. Astrophys. Suppl. Ser. 146, 73-90]. From the new catalog we find a linear cumulative luminosity function (CLF) of slope 0.54±0.05 for JF comets with q?2.5 AU. By considering this CLF combined with the few measured geometric albedos with their respective uncertainties, and assuming a canonical albedo of 0.035±0.012 for those comets with undetermined albedos, we derive a cumulative size distribution that follows a power-law of index −2.7±0.3. The slope is similar to that derived from some theoretical collisional models and from some populations of Solar System bodies like the trans-neptunian objects. We also discuss and compare our size distribution with those by other authors that have recently appeared in the literature. Some striking differences in the computed slopes are explained in terms of biases in the studied samples, the different weights given to the brightest members of the samples, and discrepancies in the values of a few absolute nuclear magnitudes. We also compute sizes and fractions of active surface area of JF comets from their estimated absolute nuclear magnitudes and their water production rates. With the outgassing model that we use, about 60% of the computed fractions f of active surface area are found to be smaller than 0.2, with one case (28P/Neujmin 1) of no more than 0.001, which suggests that JF comets may transit through stages of very low activity, or even dormancy. There is an indication that JF comets with radii RN?3 km have active fractions f?0.01, which might be due to the rapid formation of insulating dust mantles on larger nuclei.     

The size-frequency distribution of dormant Jupiter family comets

We estimate the total number and the slope of the size-frequency distribution (SFD) of dormant Jupiter family comets (JFCs) by fitting a one-parameter model to the known population. We first select 61 near-Earth objects (NEOs) that are likely to be dormant JFCs because their orbits are dynamically coupled to Jupiter [Bottke, W.F., Morbidelli, A., Jedicke, R., Petit, J., Levison, H.F., Michel, P., Metcalfe, T.S., 2002a. Icarus 156, 399-433]. Then, from the numerical simulations of Levison and Duncan [1997. Icarus 127, 13-32], we construct an orbit distribution model for JFCs in the NEO orbital element space. We assume an orbit-independent SFD for all JFCs, the slope of which is our unique free parameter. Finally, we compute observational biases for dormant JFCs using a calibrated NEO survey simulator [Jedicke, R., Morbidelli, A., Spahr, T., Petit, J., Bottke, W.F., 2003. Icarus 161, 17-33]. By fitting the biased model to the data, we estimate that there are ∼75 dormant JFCs with H<18 in the NEO region and that the slope of their cumulative SFD is −1.5±0.3. Our slope for the SFD of dormant JFCs is very close to that of active JFCs as determined by Weissman and Lowry [2003. Lunar Planet. Sci. 34. Abstract 2003]. Thus, we argue that when JFCs fade they are likely to become dormant rather than to disrupt and that the fate of faded comets is size-independent. Our results imply that the size distribution of the JFC progenitors—the scattered disk trans-neptunian population—either (i) has a similar and shallow SFD or (i^′) is slightly steeper and physical processes acting on the comets in a size-dependent manner creates the shallower active comet SFD. Our measured slope, typical of collisionally evolved populations with a size-dependent impact strength [Benz, W., Asphaug, E., 1999. Icarus 142, 5-20], suggests that scattered disk bodies reached collisional equilibrium inside the protoplanetary disk prior to their removal from the planetary region.     

The population of faint Jupiter family comets near the Earth

We study the population of faint Jupiter family comets (JFCs) that approach the Earth (perihelion distances q<1.3 AU) by applying a debiasing technique to the observed sample. We found for the debiased cumulative luminosity function (CLF) of absolute total magnitudes H₁₀ a bimodal distribution in which brighter comets (H₁₀?9) follow a linear relation with a steep slope α=0.65±0.14, while fainter comets follow a much shallower slope α=0.25±0.06 down to H₁₀∼18. The slope can be pushed up to α=0.35±0.09 if a second break in the H₁₀ distribution to a much shallower slope is introduced at H₁₀∼16. We estimate a population of about 10³ faint JFCs with q<1.3 AU and 10<H₁₀<15 (radii ∼0.1-0.5 km). The shallowness of the CLF for faint near-Earth JFCs may be explained either as: (i) the source population (the scattered disk) has an equally very shallow distribution in the considered size range, or (ii) the distribution is flattened by the disintegration of small objects before that they have a chance of being observed. The fact that the slope of the magnitude distribution of the faint active JFCs is very similar to that found for a sample of dormant JFCs candidates suggests that for a surviving (i.e., not disintegrated) object, the probability of becoming dormant versus keeping some activity is roughly size independent.     

Secular light curves of comets, II: 133P/Elst-Pizarro, an asteroidal belt comet

We present the secular light curve (SLC) of 133P/Elst-Pizarro, and show ample and sufficient evidence to conclude that it is evolving into a dormant phase. The SLC provides a great deal of information to characterize the object, the most important being that it exhibits outburst-like activity without a corresponding detectable coma. 133P will return to perihelion in July of 2007 when some of our findings may be corroborated. The most significant findings of this investigation are: (1) We have compiled from 127 literature references, extensive databases of visual colors (37 comets), rotational periods and peak-to-valley amplitudes (64 comets). 2-Dimensional plots are created from these databases, which show that comets do not lie on a linear trend but in well defined areas of these phase spaces. When 133P is plotted in the above diagrams, its location is entirely compatible with those of comets. (2) A positive correlation is found between cometary rotational periods and diameters. One possible interpretation suggest the existence of rotational evolution predicted by several theoretical models. (3) A plot of the historical evolution of cometary nuclei density estimates shows no trend with time, suggesting that perhaps a consensus is being reached. We also find a mean bulk density for comets of 〈ρ〉=0.52±0.06 g/cm³. This value includes the recently determined spacecraft density of Comet 9P/Tempel 1, derived by the Deep Impact team. (4) We have derived values for over 18 physical parameters, listed in the SLC plots, Figs. 6-9. (5) The secular light curve of 133P/Elst-Pizarro exhibits a single outburst starting at +42±4 d (after perihelion), peaking at LAG=+155±10 d, duration 191±11 d, and amplitude 2.3±0.2 mag. These properties are compatible with those of other low activity comets. (6) To explain the large time delay in maximum brightness, LAG, two hypothesis are advanced: (a) the existence of a deep ice layer that the thermal wave has to reach before sublimation is possible, or (b) the existence of a sharp polar active region pointing to the Sun at time = LAG, that may take the form of a polar ice cap, a polar fissure or even a polar crater. The diameter of this zone is calculated at ∼1.8 km. (7) A new time-age is defined and it its found that T-AGE = 80 cy for 133P, a moderately old comet. (8) We propose that the object has its origin in the main belt of asteroids, thus being an asteroid-comet hybrid transition object, an asteroidal belt comet (ABC), proven by its large density. (9) Concerning the final evolutionary state of this object, to be a truly extinct comet the radius must be less than the thermal wave depth, which at 1 AU is ∼250 m (at the perihelion distance of 133P the thermal wave penetrates only ∼130 m). Comets with radius larger than this value cannot become extinct but dormant. Thus we conclude that 133P cannot evolve into a truly extinct comet because it has too large a diameter. Instead it is shown to be entering a dormant phase. (10) We predict the existence of truly extinct comets in the main belt of asteroids (MBA) beginning at absolute magnitude ∼21.5 (diameter smaller than ∼190 m). (11) The object demonstrates that a comet may have an outburst of ∼2.3 mag, and not show any detectable coma. (12) Departure from a photometric R⁺² law is a more sensitive method (by a factor of 10) to detect activity than star profile fitting or spectroscopy. (13) Sufficient evidence is presented to conclude that 133P is the first member of a new class of objects, an old asteroidal belt comet, ABC, entering a dormant phase.     

Plasma environment of Jupiter family comets

As any comet nears the Sun, gas sublimes from the nucleus taking dust with it. Jupiter family comets are no exception. The neutral gas becomes ionized, and the interaction of a comet with the solar wind starts with ion pickup. This key process is also important in other solar system contexts wherever neutral particles become ionized and injected into a flowing plasma such as at Mars, Venus, Io, Titan and interstellar neutrals in the solar wind. At comets, ion pickup removes momentum and energy from the solar wind and puts it into cometary particles, which are then thermalised via plasma waves. Here we review what comets have shown us about how this process operates, and briefly look at how this can be applied in other contexts. We review the processes of pitch angle and energy scattering of the pickup ions, and the boundaries and regions in the comet-solar wind interaction. We use in-situ measurements from the four comets visited to date by spacecraft carrying plasma instrumentation: 21P/Giacobini-Zinner, 1P/Halley, 26P/Grigg-Skjellerup and 19P/Borrelly, to illustrate the process in action. While, of these, comet Halley is not a Jupiter class comet, it has told us the most about cometary plasma environments. The other comets, which are from the Jupiter family, give an interesting comparison as they have lower gas production rates and less-developed interactions. We examine the prospects for Rosetta at comet Churyumov-Gerasimenko, another Jupiter family comet where a wide range of gas production rates will be studied.     

Embedded star clusters and the formation of the Oort cloud: III. Evolution of the inner cloud during the Galactic phase

In a previous publication [Brasser, R., Duncan, M.J., Levison, H.F., 2006. Icarus 184, 59-82], models of the inner Oort cloud were built which included the effect of an embedded star cluster on cometary orbits about the Sun. The main conclusions of that paper were that the formation efficiency is about 10% and the median distance of the cloud to the Sun only depends on the mean density of gas and stars the Sun encountered. Here we report on the results of simulations which followed the ensuing dynamical evolution of these comet clouds in the current Galactic environment once the Sun left the embedded star cluster. The goal is to determine whether or not the dynamical influence of passing Galactic field stars and the Galactic tidal field is sufficient to replenish the current outer cloud (semi-major axis a>20,000 AU) with enough material from the inner cloud (a<20,000 AU). Since visible new comets come directly from the outer cloud, a mass estimate only exists for the latter, with a lower limit of 1 M_⊕ [Francis, P.J., 2005. Astrophys. J. 635, 1348-1361]. Knowing the amount of expansion of the inner cloud may therefore yield an estimate of the mass of said (unseen) inner cloud. Our results indicate that typically only 10% of the comets from the inner cloud land in the outer cloud and are bound after 4.5 Gyr. If one assumes that in the extreme case all or the majority of the current population of the outer cloud has come from the inner cloud, then a typical value of the mass of the inner cloud is about 10 M_⊕. The results of [Brasser, R., Duncan, M.J., Levison, H.F., 2006. Icarus 184, 59-82] showed that ∼10% of comets from the Jupiter-Saturn region were implanted in the inner Oort cloud, which implies an uncomfortably large value of about 100 M_⊕ for the mass of solids in the primordial Jupiter-Saturn region. This extreme case might be remedied in two says: either the effect of Giant Molecular Cloud complexes on the inner Oort cloud must be much more severe than originally thought, or there was a two-stage formation process for the Oort cloud, in which the outer cloud was largely populated by comets scattered once the Sun had left its primordial birth cluster.     

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.     

The origin and distribution of the Centaur population

We analyze the Centaur population as a group of objects with perihelion distances (q) of less than 30 AU and heliocentric distances outside the orbit of Jupiter, formed by objects entering this region from the Scattered Disk (SD). We perform a numerical integration of 95 real Scattered Disk Objects (SDOs) extracted from the Minor Planet Center database and of 905 synthetic SDOs compensating for observational biases. SDOs have in the Centaur zone a mean lifetime of 72 Myr, though this number falls with a decrease of q. After this incursion, 30% of them enter the zone interior to Jupiter's orbit. We find that the contribution to the Centaur population from the SD gives a total of ∼2.8×¹⁰⁸ Centaurs with a radius R>1 km. We also propose a model for the intrinsic distribution of orbital elements of Centaurs and their distance and apparent magnitude distribution.     

Optical properties of dust from Jupiter family comets

To try to define specific physical properties of the dust of Jupiter-family comets (JFCs), we compare the light scattered by them. Amongst the more than 1000 JFCs, less than 200 are numbered, 40 of them being rather bright. In the present work we use data from the latter. In situ observations of three nuclei show low albedo surfaces. The albedo of the dust particles in the coma is low, with generally a red colour. The A(α)fρ product is a measure of cometary activity and secular changes. Images of different regions (jets and fans) give indications on the nucleus rotation and position of the emitting areas, as compared to the position of the rotation axis. Differences in physical properties between the particles in different regions are pointed out by differences in the linear polarization of the scattered light and by spectral variations in brightness and polarization. Jupiter family comets are considered as dust-poor comets. Tails and trails’ studies give an estimation of the size distribution of the particles. However the dust production rates depend on the largest particles (up to centimetre size), which are mainly observed in the trails where large dark compact particles are found. These dark particles are also responsible for the high polarization in the inner most coma of some comets. The meaning, in terms of physical properties, of the linear polarization is discussed through different examples such as 2P/Encke, 9P/Tempel 1 or the fragments of 73P/Schwassmann-Wachmann 3. Cometary outbursts and splitting events show that the properties of the dust ejected from the interior of the nucleus are similar to the ones of more active comets (new or with larger semi-major axis).     

The nuclei of Jupiter family comets: A critical review of our present knowledge

Of the currently over 300 identified Jupiter family comets (JFCs), we have estimated nucleus sizes and shapes for fewer than 70 and have detailed nucleus observations arising from spacecraft fly-bys for just 3: 19P/Borrelly (Deep Space 1), 81P/Wild 2 (Stardust), and 9P/Tempel 1 (Deep Impact). These observations reveal similarities but also significant diversity. In this review, we make a critical assessment of our knowledge of JFC nuclei and suggest a priority list for observations of the nucleus of the JFC, 67P/Churyumov-Gerasimenko, the Rosetta target comet.     

A ∼25 K temperature of formation for the submicron ice grains which formed comets

Following the observations of ice grains in cometary comae and their size distributions, we reexamined experimentally our previous conclusion that the ice grains which agglomerated to form comet nuclei were formed at ∼25 K. The suggestion of a ∼25 K formation temperature was confirmed experimentally. Moreover, we suggest that these ice grains had to be of submicron size.     

Activity of comets at large heliocentric distances pre-perihelion

We present observational data for two long-period and three dynamically new comets observed at heliocentric distances between 5.8 to 14.0 AU. All of the comets exhibited activity beyond the distance at which water ice sublimation can be significant. We have conducted experiments on gas-laden amorphous ice samples and show that considerable gas emission occurs when the ice is heated below the temperature of the amorphous-crystalline ice phase transition (T∼137 K). We propose that annealing of amorphous water ice is the driver of activity in comets as they first enter the inner Solar System. Experimental data show that large grains can be ejected at low velocity during annealing and that the rate of brightening of the comet should decrease as the heliocentric distance decreases. These results are consistent with both historical observations of distant comet activity and with the data presented here. If observations of the onset of activity in a dynamically new comet are ever made, the distance at which this occurs would be a sensitive indicator of the temperature at which the comet had formed or represents the maximum temperature that it has experienced. New surveys such as Pan STARRS, may be able to detect these comets while they are still inactive.     

The internal structure of Jupiter family cometary nuclei from Deep Impact observations: The “talps” or “layered pile” model

We consider the hypothesis that the layering observed on the surface of Comet 9P/Tempel 1 from the Deep Impact spacecraft and identified on other comet nuclei imaged by spacecraft (i.e., 19P/Borrelly and 81P/Wild 2) is ubiquitous on Jupiter family cometary nuclei and is an essential element of their internal structure. The observational characteristics of the layers on 9P/Tempel 1 are detailed and considered in the context of current theories of the accumulation and dynamical evolution of cometary nuclei. The works of Donn [Donn, B.D., 1990. Astron. Astrophys. 235, 441-446], Sirono and Greenberg [Sirono, S.-I., Greenberg, J.M., 2000. Icarus 145, 230-238] and the experiments of Wurm et al. [Wurm, G., Paraskov, G., Krauss, O., 2005. Icarus 178, 253-263] on the collision physics of porous aggregate bodies are used as basis for a conceptual model of the formation of layers. Our hypothesis is found to have implications for the place of origin of the JFCs and their subsequent dynamical history. Models of fragmentation and rubble pile building in the Kuiper belt in a period of collisional activity (e.g., [Kenyon, S.J., Luu, J.X., 1998. Astron. J. 115, 2136-2160; 1999a. Astron. J. 118, 1101-1119; 1999b. Astrophys. J. 526, 465-470; Farinella, P., Davis, D.R., Stern, S.A., 2000. In: Mannings, V., Boss, A.P., Russell, S.S. (Eds.), Protostars and Planets IV. Univ. of Arizona Press, Tucson, pp. 1255-1282; Durda, D.D., Stern, S.J., 2000. Icarus 145, 220-229]) following the formation of Neptune appear to be in conflict with the observed properties of the layers and irreconcilable with the hypothesis. Long-term residence in the scattered disk [Duncan, M.J., Levison, H.F., 1997. Science 276, 1670-1672; Duncan, M., Levison, H., Dones, L., 2004. In: Festou, M., Keller, H.U., Weaver, H.A. (Eds.), Comets II. Univ. of Arizona Press, Tucson, pp. 193-204] and/or a change in fragmentation outcome modeling may explain the long-term persistence of primordial layers. In any event, the existence of layers places constraints on the environment seen by the population of objects from which the Jupiter family comets originated. If correct, our hypothesis implies that the nuclei of Jupiter family comets are primordial remnants of the early agglomeration phase and that the physical structure of their interiors, except for the possible effects of compositional phase changes, is largely as it was when they were formed. We propose a new model for the interiors of Jupiter family cometary nuclei, called the talps or “layered pile” model, in which the interior consists of a core overlain by a pile of randomly stacked layers. We discuss how several cometary characteristics—layers, surface texture, indications of flow, compositional inhomogeneity, low bulk density low strength, propensity to split, etc., might be explained in terms of this model. Finally, we make some observational predictions and suggest goals for future space observations of these objects.     

Stardust findings favor not only the planetary origin of comets but the underlying close-binary cosmogony of the Solar System as well

We analyze findings of the Stardust mission that brought to the Earth dust from the 81P/Wild 2 coma. Just as the data obtained in the Deep Impact mission to 9P/Tempel 1, they are at odds with the universally accepted condensation/sublimation comet paradigm. They fit rather well to the approach assuming ejection of nuclei of short-period comets from Moon-like bodies of the type of Galilean satellites in rare (six to seven events in 4.5 aeons) global explosions of their massive icy envelopes saturated by 2H₂ + O₂, products of the electrolysis of ice. This approach offers an explanation, in particular, for the jet activity of comets, which is sustained by combustion of the 2H₂ + O₂ + organics mixture ignited and complemented by the solar radiation. Combustion accounts also for other observations, in particular, the presence in the dust of products of high-temperature (800-900 K) metamorphism. The presence of minerals forming at still higher temperatures (∼1400-2000 K), just as the undoubtedly planetary origin of some long-period comets arriving from the joint planeto-cometary cloud beyond Neptune, forces one, however, to invoke the close-binary cosmogony of the Solar System, which three decades ago predicted the existence of such a cloud (in the recent decade, this prediction has been substantiated by the discovery there of many dwarf planets). This cosmogony is based on the modern understanding of the processes involved in the formation of multiple stellar systems and of their gas-dynamic evolution. It considers the Jupiter-Sun system as the limiting case of a binary star and uses it as a basis for explanation of all the known observations and for prediction of the new ones to come. It provides a plausible explanation, in particular, for both the origin and capture by the Earth of the Moon as a high-temperature condensate and the formation of the Galilean satellites, which also contain inclusions of refractory minerals in the ices of their envelopes.     

Meteoroid streams that trace to candidate dormant comets

In an effort to identify space mission targets of interest, the association of known meteoroid streams with Near-Earth Objects (NEOs) was investigated. In addition to updating previous searches to include NEOs discovered up to January 1, 2007, a new dissimilarity criterion based on dynamical arguments was applied to evaluate the likelihood of each candidate association. The new criterion is based on the fact that the few established cases, such as 2003 EH1 and the Quadrantid stream, involve parent bodies that fragmented in the most recent nutation cycle of their secular orbital evolution. In established cases, the statistics speak strongly of an association due to the lack of NEOs in the a, e, i phase space occupied by these showers. The newly proposed associations are much more uncertain, because the odds of chance associations greatly increase as orbital inclination of the showers decreases. Forty-two plausible candidate dormant comets were identified, that deserve further scrutiny. Both comet and stream typically lack sufficient data to prove the association. Most candidate parent bodies pertain to NEOs with an aphelion distance just short of Jupiter's orbit, a perihelion distance near Earth orbit, and an eccentricity in the range 0.5-0.8. Surprisingly many have , which means that most candidate parent bodies are dormant Jupiter family comets that have not yet fully decoupled from Jupiter. Establishing these associations can provide further evidence that (mostly) dormant comets break frequently, making this the dominant mechanism for replenishing the zodiacal cloud.     

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.     

Gas trapping in water ice at very low deposition rates and implications for comets

The effect of water ice formation temperature and rate of ice deposition on a cold plate on the amount of trapped argon (equivalent to CO), and the ratios of Ar/Kr/Xe trapped in the water ice were studied at 50, 27 and 22 K and at ice formation rates ranging over four orders of magnitude, from 10⁻¹ to 10⁻⁵ μm min⁻¹. Contrary to our previous conclusions that cometary ices were formed at 50-60 K, we now conclude that these ices were formed at about 25 K. At 25 K the enrichment ratios for Ar, Kr, and Xe remained the same as those at 50 K, reinforcing our suggestion of cometary contribution of these noble gases to the atmospheres of Earth and Mars.     

Organics-glued dust aggregates in Halley's coma and some implications of the snow line inward motion for the formation of comets

In view of the solar nebula models, organics-glued dust aggregates (whose disintegration resulted in the two phenomena found in Halley's coma, the dust boundary and small-scale dust structures) originated due to coagulation of iceless dust particles somewhere within the snow line, and then were incorporated into Halley's nucleus as a consequence of the snow line inward motion. This implies that two types of comets exist: outer comets, formed entirely beyond the snow line, and inner comets, similar to Halley, which are bodies intermediate between outer comets and primitive asteroids. The presence of large iceless dust aggregates in nuclei of inner comets constrains the inward drift velocity of meter-sized dust bodies, which in turn implies that the radial transport of water in the solar nebula was predominantly outward. It is shown that in nuclei of inner comets: both the upper mass limit of iceless dust aggregates and the ice mantle thickness increase with decreasing formation heliocentric distance, while the cumulative mass distribution index decreases; the lower limit of the mass index is ∼0.8, and the upper limit of the ice mantle thickness is ∼10⁻³ cm (∼200 times the interstellar value); the lower limit of the latent heat of organics in organic mantles of submicron particles increases toward small heliocentric distances; the recondensation of organics combined with the growth of dust bodies leads to a fractionation of organics within iceless dust aggregates; last accreted sub-units of an aggregate are always glued by organics with the lowest value of the latent heat, which somewhat exceeds 60 kJ/mol. Based on in situ observations at Halley, the parameters characterizing iceless dust aggregates in that comet are calculated. Finally, feasible observational tests of the conclusions drawn are discussed.     

Crystalline silicates in comets: How did they form?

Two processes have been proposed to explain observations of crystalline silicate minerals in comets and in protostellar sources, both of which rely on the thermal annealing of amorphous grains. First, high temperatures generated by nebular shock processes can rapidly produce crystalline magnesium silicate grains and will simultaneously produce a population of crystalline iron silicates whose average grain size is ∼10-15% that of the magnesium silicate minerals. Second, exposure of amorphous silicate grains to hot nebular environments can produce crystalline magnesium silicates that might then be transported outward to regions of comet formation. At the higher temperatures required for annealing amorphous iron silicates to crystallinity the evaporative lifetime of the grains is much shorter than a single orbital period where such temperatures are found in the nebula. Thermal annealing is therefore unable to produce crystalline iron silicate grains for inclusion into comets unless such grains are very quickly transported away from the hot inner nebula. It follows that observation of pure crystalline magnesium silicate minerals in comets or protostars is a direct measure of the importance of simple thermal annealing of grains in the innermost regions of protostellar nebulae followed by dust and gas transport to the outer nebula. The presence of crystalline iron silicates would signal the action of transient processes such as shock heating that can produce crystalline iron, magnesium and mixed iron-magnesium silicate minerals. These different scenarios result in very different predictions for the organic content of protostellar systems.