Laboratory simulation of the physical processes occurring on and near the surfaces of comet nuclei

Abstract— Laboratory comet simulation experiments are discussed in the context of theoretical models and recent ground-based and spacecraft observations, especially the Giotto observations of P/Halley. The set-up of various comet simulation experiments is reviewed. A number of small-scale experiments have been performed in many laboratories since the early 1960s. However, the largest and most ambitious series of experiments were the comet simulation experiments known as KOSI (German = Kometen Simulation). These experiments were prompted by the appearance of Comet P/Halley in 1986 and in planning for the European Space Agency's Rossetta mission that was originally scheduled to return samples. They were performed between 1987 and 1993 using the German Space Agency's (DLR) space hardware testing facilities in Cologne. As with attempts to reproduce any natural phenomenon in the laboratory, there are deficiencies in such experiments while there are major new insights to be gained. Simulation experiments have enabled the development of methods for making comet analogues and for exploring the properties of such materials in detail. These experiments have provided new insights into the morphology and physical behavior of aggregates formed from silicate grains likely to exist in comets. Formation of a dust mantle on the surfaces and a system of ice layers below the mantle caused by chemical differentiation have been identified after the insolation of the artificial comet. The mechanisms for heat transfer between the comet's surface and its interior, the associated gas diffusion from the interior of the surface, and compositional, structural, and isotopic changes that occur near the surface have been described by modeling the experimental results. The mechanisms of the ejection of dust and ice grains from the surface and the importance of gas-drag in propelling grains have also been explored.     

All comets are born equal: infrared emission by dust as a key to comet nucleus composition

The infrared emission of various comets can be matched within the framework that all comets are made of aggregated interstellar dust. This is demonstrated by comparing results on Halley (a periodic comet), Borrelly (a Jupiter family short period comet), Hale-Bopp (a long period comet), and extra-solar comets in the β Pictoris disk. Attempts have been made to generalize the chemical composition of comet nuclei based on the observation of cometary dust and volatiles and the interstellar dust model. Finally, we deduce some of the expected dust and surface properties of comet Wirtanen from the interstellar dust model as applied to other comets.     

An idealized short-period comet model: Surface insolation,H2O flux,dust flux,and mantle evolution

A model of cometary activity is developed which integrates the feedback processes involving heat, gas, and dust transport, and dust mantle development. The model includes the effects of latitude, rotation, and spin axis orientation. Results are obtained for various grain size distributions, dust-to-ice ratios, and spin axis orientations. Attention is focused on the development, change of structure and distribution of dust mantles and their mutual interaction with ice surface temperature and gas and dust production. In this model the dust mantle controls the mechanism of gas transport not onlu by its effect on the temperature but, more importantly, by its own dynamic stability. Results suggest that an initially homogeneous short-period comet with a “cosmic” dust-to-water ice ratio, typical orbit, rotation rate, and grain size distribution would develop at most only a thin (<1 mm) cyclic mantle at all points on the nucleus. Such a fully developed temporary mantle would exist throughout the diurnal cycle only beyond ～4AU. Thus, cyclic behavior would be expected for such an idealized comet, at least for most of its lifetime. Long-term irreversible mantle development on comets with typical rotation rates was not found except regionally on Encke and also on objects with perihelia ?1.5 AU. Even in these cases, free silicate exists, after a few cycles, only as relatively rare large grains and agglomerates with radii ～1 cm scattered over exposed ice. Full mantle development would require hundreds to thousands of cycles. In the case of an initially homogeneous comet Encke, this slow incipient mantle development is shown to be the direct result of its peculiar axial orientation. High obliquity appears required for long-term mantle development for typical rotation rates and perihelia ?1.5 AU. Heat conduction into the nucleus for an incompletely mantled or bald comet has been found to be very important in maintaining relatively higher ice surface temperatures, and hence fluxes, during those portions of the diurnal and orbital cycles which would otherwise be cooler. It is also shown to be at least one cause of post perihelion brightness asymmetries, especially in lower obliquity comets. Maximum heliocentric distances at which 1-μm dust, sand, pebbles, cobbles, and boulders can be permanently ejected from the subsolar point by H₂O (CO₂) are (in AU): 6.9 (16.8), 5.2 (11.5), 1.8 (3.0), 0.21 (0.34) and 0.07 (0.11), respectively. A detailed anatomy of temperature, gas and dust fluxes vs latitude and longitude for a homogeneous rotating comet with fixed axis is given for comparison with future observations. Most H₂O flux histories deduced from brightness data are found to be in reasonable agreement with the model, allowing for uncertainty in radius and albedo. A clear exception is Encke. It is shown that the large discrepancy between Encke's observed and model predicted fluxes, based on radar cross section, can be used to evaluate the extent of exposed ice (<10%). The model is then used to place an active area so as to explain a reported sharp drop in flux on approach to the Sun at 0.78 AU. An active area or areas, <10% of the comet's surface, centered near 65°N latitude appears indicated. Although cyclic mantles are generally indicated for the set of parameters we used, our results show that a global mantle only 1 to 3 cm thick (depending on the orbit) consisting of a full range of grain sizes can cause irresversible evolution to a noncometary body. We investigated the long-term evolution of such a postulated initially thinly mantled cometary object. It was found that after the first few passes and until the end of its dynamic lifetime the object averaged <3 × 10^?12 g cm^?1 sec^?1 H₂O flux. Therefore, if cometary objects evolve into Apollo asteroids, ice should always be accessible within 10 m of the surface despite numerous close perihelion passages. The possible impact of factors not included in the model, such as initial inhomogeneities, coma scattering of radiation, and global redistribution of ejected silicate around the nucleus, are discussed.     

Modeling the Comet Nucleus Interior

Numerical simulation of the structure and evolution of a comet nucleus is reviewed both from the mathematical and from the physical point of view. Various mathematical procedures and approximations are discussed, and different attempts to model the physical characteristics of cometary material, such as thermal conductivity, or permeability to gas flow, are described. The evolution and activity of comets is shown to depend on different classes of parameters: Defining parameters, such as size and orbit, structural parameters, such as porosity and composition, and initial parameters, such as temperature and live radio isotope content. The latter are related to the formation of comets. Despite the large number of parameters, general conclusions, or common features, appear to emerge from the numerous model calculations — for different comets — performed to date. Thus, the stratified structure of comet nuclei, volatile depletion, and the role of crystallization of ice in cometary outbursts are discussed. Finally, an evolution model applied to comet C/1995 O1 Hale-Bopp — using different assumptions — is described and analysed in the light of observations.     

Models of P/Wirtanen nucleus: active regions versus non-active regions

From the current understanding we know that comet nuclei have heterogeneous compositions and complex structures. It is believed that cometary activity is the result of a combination of physical processes in the nucleus, like sublimation and recondensation of volatile ices, dust grains release, phase transition of water ice, depletion of the most volatile components in the outer layers and interior differentiation.The evolution of the comet depends on the sublimation of ices and the release of different gases and dust grains: the formation of a dust crust, the surface erosion and the development of the coma are related to the gas fluxes escaping from the nucleus. New observations, laboratory experiments and numerical simulations suggest that the gas and dust emissions are locally generated, in the so-called active regions. This localized activity is probably superimposed to the global nucleus activity. The differences between active and inactive regions can be attributed to differences in texture and refractory material content of the different areas.In this paper we present the results of numerical models of cometary nucleus evolution, developed in order to understand which are the processes leading to the formation of active and non-active regions on the cometary surface. The used numerical code solves the equations of heat transport and gas diffusion within a porous nucleus composed of different ices—such as water (the dominant constituent), CO₂, CO- and of dust grains embedded in the ice matrix.By varying the set of physical parameters describing the initial properties of comet P/Wirtanen, the different behaviour of the icy and dusty areas can be followed.Comet P/Wirtanen is the target of the international ROSETTA mission, the cornerstone ESA mission to a cometary nucleus. The successful design of ROSETTA requires some knowledge of comet status and activity: surface temperatures, amount of active and inactive surface areas, gas production rate and dust flux.     

Short-period Jupiter family comets after Stardust

The NASA Stardust mission has provided for laboratory study an extensive data set of cometary dust of known provenance (from comet 81P/Wild 2) yielding detailed insights into the composition of the comet. Combined with the results of data from other missions to short-period Jupiter family comets (JFC), this has greatly deepened the understanding of such objects. If depressions on the surface of comet 81P/Wild 2 are all taken as evidence of impact cratering, their number suggests a long occupancy in the outer region of the Solar System. The dust from comet 81P/Wild 2 has been shown to be heavily deficient in pre-Solar grains and rich in materials formed at high temperatures in the inner Solar System. Although it is too early to know if this is typical of JFC, it does argue for rapid and thorough mixing of materials in the disk on timescales related to comet formation, and may also suggest outward migration of small icy bodies after their formation. Thus, instead of providing mainly new knowledge of the pre-Solar materials expected to be rich in comets, Stardust and comet 81P/Wild 2 have instead focussed attention on large-scale transport processes during the critical period when cometary parent bodies were forming in the early Solar System.     

Large cometary grains – their destruction and changes in the luminosity of comets

The catastrophic thermodynamic destruction of large cometary heterogeneous grains lying on the surface of a comet nucleus is examined. The core–mantle grain-structure model is assumed. Grain fragmentation as an explanation of sudden changes in cometary brightness is proposed. The approach presented to the problem of cometary outbursts is a development of a previous author's paper. The proposed mechanism is based on the idea of thermodynamical destruction of heterogeneous cometary grains. Numerical simulations have been carried out for a wide range of values of physical characteristics of cometary material. The results obtained are consistent with observational data. The main conclusion of this paper is that thermodynamical fragmentation of large grains can explain variations in brightness and also outbursts of comets.     

On the importance of dust in cometary nuclei

The icy conglomerate model introduced by Whipple more than 40 years ago has been widely accepted in cometary science because it is able to describe numerous cometary phenomena. In this model comets are described as a conglomerate of ices and dust where the ices represent the major component. However, some recent observations seem to favour dust rich comets. The purpose of this paper is to summarize the observational facts supporting the dominance of refractories in comets and to discuss the consequences of a dust dominated nucleus for cometary physics.     

On the importance of dust in cometary nuclei

The icy conglomerate model introduced by Whipple more than 40 years ago has been widely accepted in cometary science because it is able to describe numerous cometary phenomena. In this model comets are described as a conglomerate of ices and dust where the ices represent the major component. However, some recent observations seem to favour dust rich comets. The purpose of this paper is to summarize the observational facts supporting the dominance of refractories in comets and to discuss the consequences of a dust dominated nucleus for cometary physics.     

A tale of two very different comets: ISO and MSX measurements of dust emission from 126P/IRAS (1996) and 2P/Encke (1997)

We present the characteristics of the dust comae of two comets, 126P/IRAS, a member of the Halley family (a near-isotropic comet), and 2P/Encke, an ecliptic comet. We have primarily used mid- and far-infrared data obtained by the ISOPHOT instrument aboard the Infrared Space Observatory (ISO) in 1996 and 1997, and mid-infrared data obtained by the SPIRIT III instrument aboard the Midcourse Space Experiment (MSX) in 1996. We find that the dust grains emitted by the two comets have markedly different thermal and physical properties. P/IRAS's dust grain size distribution appears to be similar to that of fellow family member 1P/Halley, with grains smaller than 5 microns dominating by surface area, whereas P/Encke emits a much higher fraction of big (20 μm and higher) grains, with the grain mass distribution being similar to that which is inferred for the interplanetary dust population. P/Encke's dearth of micron-scale grains accounts for its visible-wavelength classification as a “gassy” comet. These conclusions are based on analyses of both imaging and spectrophotometry of the two comets; this combination provides a powerful way to constrain cometary dust properties. Specifically, P/IRAS was observed preperihelion while 1.71 AU from the Sun, and seen to have a 15-arcmin long mid-infrared dust tail pointing in the antisolar direction. No sunward spike was seen despite the vantage point being nearly in the comet's orbital plane. The tail's total mass at the time was about 8×10⁹ kg. The spectral energy distribution (SED) is best fit by a modified greybody with temperature T=265±15 K and emissivity ε proportional to a steep power law in wavelength λ: ε∝λ^−α, where α=0.50±0.20(2σ). This temperature is elevated with respect to the expected equilibrium temperature for this heliocentric distance. The dust mass loss rate was between 150-600 kg/s (95% confidence), the dust-to-gas mass loss ratio was about 3.3, and the albedo of the dust was 0.15±0.03. Carbonaceous material is depleted in the comet's dust by a factor of 2-3, paralleling the C₂ depletion in P/IRAS's gas coma. P/Encke, on the other hand, observed while 1.17 AU from the Sun, had an SED that is best fit by a Planck function with T=270±15 K and no emissivity falloff. The dust mass loss rate was 70-280 kg/s (95% confidence), the dust-to-gas mass loss ratio was about 2.3, and the albedo of the dust was about 0.06±0.02. These conclusions are consistent with the strongly curved dust tail and bright dust trail seen by Reach et al. (2000; Icarus 148, 80) in their ISO 12-μm imaging of P/Encke. The observed differences in the P/IRAS and P/Encke dust are most likely due to the less evolved and insolated state of the P/IRAS nuclear surface. If the dust emission behavior of P/Encke is typical of other ecliptic comets, then comets are the major supplier of the interplanetary dust cloud.     

The destruction of the cometary nucleus – The case of 73P/Schwassmann-Wachmann

The paper presents an analysis of the actual brightness change of comet 73P/Schwassmann-Wachmann, which took place in 1995. The consequence of a cometary outburst is the destruction of a fragment of its surface. This causes the emission of comet material from both the surface and from exposed subsurface layers. Therefore, the calculations take into account the scattering cross-sections that come from ice and dust particles. It was assumed that the dust particles are silicates which are characterized by high irregularity of their structure. This assumption is a consequence of the analysis of the results provided by the Rosetta mission to the comet 67P/Churyumov-Gerasimenko. The main factor determining the amplitude of a cometary outburst is the mass ejected as well as the loss of ice that holds the individual nucleus structures together. Consequently, this phenomenon can significantly contribute to the destruction and even decay of the cometary nucleus.     

Dust in cometary comae: Present understanding of the structure and composition of dust particles

In situ probing of a very few cometary comae has shown that dust particles present a low albedo and a low density, and that they consist of both rocky material and refractory organics. Remote observations of solar light scattered by cometary dust provide information on the properties of dust particles in the coma of a larger set of comets. The observations of the linear polarization in the coma indicate that the dust particles are irregular, with a size greater (on the average) than about 1 μm. Besides, they suggest, through numerical and experimental simulations, that both compact grains and fluffy aggregates (with a power law of the size distribution in the −2.6 to −3 range), and both rather transparent silicates and absorbing organics are present in the coma. Recent analysis of the cometary dust samples collected by the Stardust mission provide a unique ground truth and confirm, for comet 81P/Wild 2, the results from remote sensing observations. Future space missions to comets should, in the next decade, lead to a more precise characterization of the structure and composition of cometary dust particles.     

Millimeter-wave experiments for cometary space missions

Predicted brightness temperatures for a variety of cometary nucleus models, consisting of homogeneous layers comprised of mixtures of water ice and refractory grains, are presented as functions of wavelength. These illustrative spectra are computed using simple radiative transfer techniques adapted from modeling of terrestrial ice and snow fields. The computed millimeter-wave spectra are sensitive to the values of physically significant nucleus parameters such as crust thickness, the subsurface temperature gradient, and the boundary temperature of the sublimating surface. It appears that millimeter-wave sensing from an interplanetary spacecraft is an effective means for distinguishing between alternate models of the nucleus and for inferring the rough physical state of substrata; modern theories on the nature of the nucleus indicate that sublimation from the substrata provides the gas phase cometary volatiles that are actually observed from ground-based and Earth-orbiting instruments. Antenna beam dilution is a major obstacle for ground-based molecular spectral line radio observations (e.g., water and ammonia) of comets but a modest millimeter-wave radiometer system in the near vicinity of the comet would not be subject to this problem. Such a system can make definitivebservations of several candidate parent molecules in the gas phase and should contribute to the understanding of the physics of the inner coma.     

Comet Splitting – Observations and Model Scenarios

Splitting events affect cometary nuclei to a different level of severity ranging from complete disruption of the nucleus (e.g., C/1999 S4 LINEAR) to separation of major fragments (e.g., 73P/Schwassmann-Wachmann 3) and spill-offs of smaller boulders (e.g., C/2001 A2 LINEAR).Fragmentation of comets produces secondary products over a wide range of sizes (from cometesimals to sub-micron dust). It is detectable through the presence of fragments (with own comae and tails) in the coma of the parent nucleus, through outbursts in its activity and through arc-lets (“coma wings”)associated with fragments. The secondaries have different life times and show different non-gravitational forces. Nucleus splitting is also considered to generate whole families of comets (Kreutz group) or — if gravitational bound — multiple nuclei (e.g., C/1995 O1 Hale-Bopp). It may explain the striae phenomena seen in dust tails of bright comets (C/1995 O1 Hale-Bopp) and the detection of chains of impact craters onother bodies in the solar system. As process of significant mass loss it is relevant for the scenario of nucleus extinction, at the same time it also plays a role for the number statistics of existing (observable) comets and for the size distribution of comet nuclei. Various model scenarios for nucleus splitting are proposed: tidal disruption, rotational splitting, break-up due to internal gas pressure, fragmentation due to collision with other bodies. Only in one case, Comet D/1993 F1Shoemaker-Levy 9, the physical process of fragmentation could be undoubtedly identified. In any case, comet splitting provides important insights inthe internal structure, surface layering and chemistry of comet nuclei.     

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.     

Is the D/H Ratio in the Comet Coma Equal to the D/H Ratio in the Comet Nucleus?

We present a simple, semianalytic model of the vaporization of H₂O and HDO ice from a comet nucleus. We use this model to show that the flux of HDO relative to H₂O can be much higher, at times, than would be expected from the D/H ratio in the nuclear ice itself. This effect varies with position in the comet's orbit. It is negligible sufficiently near the Sun but could lead to erroneous interpretations of the primordial D/H ratio in cometary ice if measurements are made in other parts of the cometary orbit.     

The destruction of cometary grains and changes in the luminosity of comets

One explanation of the sudden changes in the brightness of comets is proposed based on the author's earlier suggestions involving the fragmentation of cometary grains. Within the inner coma, a core‐mantle model of the structure of grains is assumed. The proposed mechanism is a combination of electrostatic stress and thermodynamical fragmentation of the cometary grains water‐ice mantle. It has been shown that the vapour pressure of volatile inclusions placed in the waterice mantle of grains can increase sufficiently to cause their fragmentation. It takes place before grains can completely sublime into the vacuum away. Numerical calculations have been carried out for a large range of values of probable physical characteristics of cometary material. The proposed approach yields increases in cometary brightness consistent with observations of typical cometary outbursts. It is concluded that this approach can provide an explanation of the sudden change in activity of comets for a wide range of heliocentric distances (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Do comets have satellites?

The recent evidence that many minor planets may have satellites, together with recently iscovered physical, chemical, and lightcurve similarities between minor planets and comets, lead naturally to the question, “Might comets have satellites also?” This paper explores several puzzling features of comets which do not fit easily into conventional cometary models, but which can be satisfactorily explained if it is assumed that comets have a full range of gravitationally bound masses, from dust size to the size of the nucleus, in orbit around the principal nucleus. This discussion also implies a higher probability of destruction of a spacecraft near a comet than is usually assumed.     

Dust tail of the active distant Comet C/2003 WT42 (LINEAR) studied with photometric and spectroscopic observations

We present the study of dust environment of dynamically new Comet C/2003 WT42 (LINEAR) based on spectroscopic and photometric observations. The comet was observed before and after the perihelion passage at heliocentric distances from 5.2 to 9.5 AU. Although the comet moved beyond the zone where water ice sublimation could be significant, its bright coma and extended dust tail evidenced the high level of physical activity. Afρ values exceeded 3000 cm likely reaching its maximum before the perihelion passage. At the same time, the spectrum of the comet did not reveal molecular emission features above the reflected continuum. Reddening of the continuum derived from the cometary spectrum is nonlinear along the dispersion with the steeper slop in the blue region. The pair of the blue and red continuum images was analyzed to estimate a color of the comet. The mean normalized reflectivity gradient derived from the innermost part of the cometary coma equals to 8% per 1000 Å that is typical for Oort cloud objects. However, the color map shows that the reddening of the cometary dust varies over the coma increasing to 15% per 1000 Å along the tail axis. The photometric images were fitted with a Monte Carlo model to construct the theoretical brightness distribution of the cometary coma and tail and to investigate the development of the cometary activity along the orbit. As the dust particles of distant comets are expected to be icy, we propose here the model, which describes the tail formation taking into account sublimation of grains along their orbits. The chemical composition and structure of these particles are assumed to correspond with Greenberg’s interstellar dust model of comet dust. All images were fitted with the close values of the model parameters. According to the results of the modeling, the physical activity of the comet is mainly determined by two active areas with outflows into the wide cones. The obliquity of the rotation axis of the nucleus equals to 20° relative to the comet’s orbital plane. The grains occupying the coma and tail are rather large amounting to 1 mm in size, with the exponential size distribution of a^−4.5. The outflow velocities of the dust particles vary from a few centimeters to tens of meters per second depending on their sizes. Our observations and the model findings evidence that the activity of the nucleus decreased sharply to a low-level phase at the end of April–beginning of May 2007. About 190 days later, in the first half of November 2007 the nucleus stopped any activity, however, the remnant tail did not disappear for more than 1.5 years at least.     

On the Role of Dust Mass Loss in the Evolution of Comets and Dusty Disk Systems

We present the results of an ongoing optical/thermal infraredphotometric imaging survey of thedust mass loss and nuclear size ofcomets. We find an evolution with time of the kind of dust emittedfrom a comet's surface. Our results indicate that the mass loss ratefrom the short period comets alone is enough to supply the interplanetary dust (IPD) complexagainst losses. We conclude that the nature of the IPD cloud haschanged over time, with small particles dominating early, and largerparticles dominating in the present era. The rate of destruction of short period (SP)comets due to sublimation mass loss is found to be low compared to therate of surface mantle formation, but fast compared to the rate of their dynamical removal from the inner solar system.