Structure and density of cometary nuclei

Abstract— Understanding the nature of the cometary nucleus remains one of the major problems in solar system science. Whipple's (1950) icy conglomerate model has been very successful at explaining a range of cometary phenomena, including the source of cometary activity and the nongravitational orbital motion of the nuclei. However, the internal structure of the nuclei is still largely unknown. We review herein the evidence for cometary nuclei as fluffy aggregates or primordial rubble piles, as first proposed by Donn et al. (1985) and Weissman (1986). These models assume that cometary nuclei are weakly bonded aggregations of smaller, icy‐conglomerate planetesimals, possibly held together only by self‐gravity. Evidence for this model comes from studies of the accretion and subsequent evolution of material in the solar nebula, from observations of disrupted comets, and in particular comet Shoemaker‐Levy 9, from measurements of the ensemble rotational properties of observed cometary nuclei, and from recent spacecraft missions to comets. Although the evidence for rubble pile nuclei is growing, the eventual answer to this question will likely not come until we can place a spacecraft in orbit around a cometary nucleus and study it in detail over many months to years. ESA's Rosetta mission, now en route to comet 67P/Churyumov‐Gerasimenko, will provide that opportunity.     

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)     

MASS FLUX OF ASTEROIDAL ORIGIN METEOROIDS ON PERIODIC COMET NUCLEI

We examine the potential contamination of cometary nuclei through impacts from asteroidal origin meteoroids. The paper uses a simple model and has the goal of determining whether asteroidal contamination is potentially significant. We assume a meteoroid power law mass distribution with index values in the range from s=1.83 to s=2.09. We used maximum and minimum models which we believe will bracket the true meteoroid mass distribution. We identify those comets which are expected to be most significantly contaminated, and find values of up to 3.6 kg of asteroidal meteoroid impact per square meter of the cometary surface per orbital revolution. This is less than the expected mass loss per perihelion passage for most comets. Therefore any remnant effects of the contamination will depend on the penetration depth of the meteoroids in the cometary nucleus, and possibly on the distribution of active and inactive areas on cometary nuclei. We present a simple model which suggests that even small meteoroids will embed relatively deeply into a cometary nucleus.     

Thermal destruction of cometary materials – application to activity of comets

Destruction mechanisms connected with thermodynamical behaviour of cometary material are reviewed with a special consideration of their effects on activity of comets. Consequences of thermal stresses which occur in the interior of a comet are discussed with reference to changes in the cometary brightness. Moreover, thermal destruction of grains placed in the head of the comet as well as on the surface of the nucleus is considered. It has been shown that the destruction of the cometary material can lead to an essential increase in the activity of the comet. Calculations have been carried out for a large assumed range of cometary parameters. The obtained simulated changes in the brightness of comets are consistent with the ones observed during the real variations and outbursts of brightness. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Dynamics of high‐temperature materials delivered by jets to the outer solar nebula

Abstract— The presence of high‐temperature materials in the Stardust collection that are isotopically similar to those seen in chondritic meteorites argues for the outward transport of materials from the hot, inner region of the solar nebula to the region where comets formed. A number of mechanisms have been proposed to be responsible for this transport, with a number of models being developed to show that such outward transport is possible. However, these models have not examined in detail how these grains are transported after they have been delivered to the comet formation region or how they may be distributed in the cometary nuclei that form. Here, the dynamical evolution of crystalline silicates injected onto the surface of the solar nebula as proposed by jet models for radial transport is considered. It is generally found that crystalline grains should be heterogeneously distributed within the population of comets and within individual cometary nuclei. In order to achieve a homogeneous distribution of such grains, turbulence must be effective at mixing the crystalline silicates with native, amorphous grains on fine scales. However, this turbulent mixing would serve to dilute the crystalline silicates as it would redistribute them over large radial distances. These results suggest that it is difficult to infer the bulk properties of Wild 2 from the Stardust samples, and that the abundance of crystalline grains in these samples cannot alone be used to rule out or in favor of any of the different radial transport models that have been proposed.     

Processing of cometary grains at the nucleus surface

Cometary material inevitably undergoes chemical changes before and on leaving the nucleus. In seeking to explain comets as the origin of many IDPs (interplanetary dust particles), an understanding of potential surface chemistry is vital. Grains are formed and transformed at the nucleus surface; much of the cometary volatiles may arise from the organic material. In cometary near-surface permafrost, one expects cryogenic chemistry with crystal growth and isotope. This could be the hydrous environment where IDPs form. Seasonal and geographic variations imply a range of environmental conditions and surface evolution. Interplanetary dust impacts and electrostatic forces also have roles in generating cometary dust. The absence of predicted cometary dust ‘envelopes’ is compatible with the wide range of particle structures and compositions. Study of IDPs would distinguish between this model and alternatives that see comets as aggregates of core-mantle grains built in interstellar clouds.     

Large Dust Grains Around Cometary Nuclei

Large amounts of particles ejected from the nucleus surface are present in the vicinity of the cometary nuclei when comets are near the Sun (at heliocentric distances ≤2 AU). The largest dust grains ejected may constitute a hazard for spatial vehicles. We tried to obtain the bounded orbits of those particles and to investigate their stability along several orbital periods. The model includes the solar and the cometary gravitational forces and the solar radiation pressure force. The nucleus is assumed to be spherical. The dust grains are also assumed to be spherical, and radially ejected. We include the effects of centrifugal forces owing to the comet rotation. An expression for the most heavy particles that can be lifted is proposed. Using the usual values adopted for the case of Halley’s comet, the largest grains that can be lifted have a diameter about 5 cm, and the term due to the rotation is negligible. However, that term increases the obtained value for the maximum diameter of the lifted grain in a significant amount when the rotation period is of the order of a few hours.     

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.     

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.     

On the formation of cometary nuclei in dense globules

Evolution of cometary orbits by planetary perturbations, weakly hyperbolic original orbits of comets calculated by Marsdenet al. (1978) are taken to indicate the interstellar origin of comets, and the possible formation of cometary nuclei in interstellar globules is discussed. The process is sedimentation of dust grains. It is shown that if a globule is at 40 K, its lifetime is sufficiently long to allow the sedimentation.     

Possibilities of using extreme observations to estimate cometary nucleus sizes

Not considering very rare in situ measurements of cometary nuclei, observations of comets at large heliocentric distances are the only direct source of our knowledge on their sizes. Observations of a cometary nucleus in pure reflected sunlight, at the time when coma is absent, are the way in which the nucleus size can be estimated. Probabilities that extreme observations represent non—active stages of cometary nuclei and also reliability of derived cometary nucleus sizes are investigated. Statistical analysis is based on a sample of 2842 photometric observations of 67 long-period comets observed at large heliocentric distances. For any long-period comet, there is a probability of 2:3 that the sizes derived on the basis of observations at extreme distances are in good agreement with the real nucleus sizes. For new comets in Oort's sense the probability is 3:4 independent of investigated arcs of orbits. For old comets a chance to estimate correct sizes is 1:2 but on the pre-perihelion arc only 1:3. It is also demonstrated that a premature start of activity prior to perihelion or a longer fading after perihelion is more frequent than a short-time isolated activity at large heliocentric distances.     

Protosolar nebula and composition of comets

Comets seem to be composed of matter, which is supposed to have the same molecular composition as protosolar nebula. Although there are no unbiased evidence that cometary nuclei retain the molecular composition inherited from the protosolar cloud, the observed properties of comets indicate that there is at least a resemblance between cometary composition and the material properties of dense interstellar clouds. Therefore the origin of comets could be searched in the cold stages of the protosolar nebula and molecular abundances of grain mantles in this nebula may be similar to those in the cometary dust. It is suggested that comets may contain pristine, virtually unaltered protosolar material and their study might be very relevant way to more information about processes in early stages of the solar nebula. Our knowledge about composition of the cometary nucleus is still relatively scarce, but we can partly deduce it from data obtained either by ground-based spectroscopy or by in situ mass spectrometry from space experiments. Most important were the discovery of fluffy CHON particles composed partly or even completely from compounds containing light elements. No consensus concerning the presence of interstellar pristine matter in comet has been reached from various approaches to determine the relationship between comets and interstellar grains. Most of these studies are based on infrared spectroscopy. Another method is the comparison on the chemical models of the protosolar nebula with the volatile compounds of the cometary nuclei. Both gas-phase and grain-surface chemistry are considered and initial gas-phase atomic abundances are assumed to be protosolar. The cometary matter is certainly not identical with the typical material of dense interstellar cool dense clouds, but it is closer to it than any other type of matter in solar system so far accessible to us. The data from comets combined with models of chemical evolution of matter in environment similar as prevailed the early stage of presolar nebula may at least impose constrains on the condition for comet formation. Here presented study is a preliminary contribution to such studies.     

Hamiltonian Formulation and Perturbations for Dust Motion Around Cometary Nuclei

In this paper we analyze the dynamical behavior of large dust grains in the vicinity of a cometary nucleus. To this end we consider the gravitational field of the irregularly shaped body, as well as its electric and magnetic fields. Without considering the effect of gas friction and solar radiation, we find that there exist grains which are static relative to the cometary nucleus; the positions of these grains are the stable equilibria. There also exist grains in the stable periodic orbits close to the cometary nucleus. The grains in the stable equilibria or the stable periodic orbits won’t escape or impact on the surface of the cometary nucleus. The results are applicable for large charge dusts with small area-mass ratio which are near the cometary nucleus and far from the Solar. It is found that the resonant periodic orbit can be stable, and there exist stable non-resonant periodic orbits, stable resonant periodic orbits and unstable resonant periodic orbits in the potential field of cometary nuclei. The comet gravity force, solar gravity force, electric force, magnetic force, solar radiation pressure, as well as the gas drag force are all considered to analyze the order of magnitude of these forces acting on the grains with different parameters. Let the distance of the dust grain relative to the mass centre of the cometary nucleus, the charge and the mass of the dust grain vary, respectively, fix other parameters, we calculated the strengths of different forces. The motion of the dust grain depends on the area-mass ratio, the charge, and the distance relative to the comet’s mass center. For a large dust grain (> 1 mm) close to the cometary nucleus which has a small value of area-mass ratio, the comet gravity is the largest force acting on the dust grain. For a small dust grain (< 1 mm) close to the cometary nucleus with large value of area-mass ratio, both the solar radiation pressure and the comet gravity are two major forces. If the a small dust grain which is close to the cometary nucleus have the large value of charge, the magnetic force, the solar radiation pressure, and the electric force are all major forces. When the large dust grain is far away from the cometary nucleus, the solar gravity and solar radiation pressure are both major forces.     

Are the cometary outbursts caused by impacts with interplanetary vagabonds probable?

The possibility of impacts and their results in relation to the cometary outbursts between comets and other small bodies in the solar system has been investigated. Taking into consideration certain physical features of cometary nuclei and impacting bodies, the probability of impacts of small bodies moving in the main asteroid belt with hypothetical comets which represent three types: Jupiter family comets, Halley family comets and long period comets has been computed. The probability of impacts between comets and meteoroids at large heliocentric distances has also been estimated. Potential consequences of these events in relation to outbursts of the cometary brightness have been discussed. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Heat content and evolution of cometary nuclei

In continuation of an earlier study of the influence of phase transitions on the thermal behavior of cometary nuclei, the heat flux into nuclei at various distances from the Sun before and after perihelion has been investigated for the isothermal case and for the fixed subsolar point. It turns out that this heat flux may be a large fraction of the incident solar heat input, so that the surface temperature and the associated rate of evaporation are lower than usually calculated. The effect is strongly dependent on the porosity of the nucleus. The surface temperature of the nucleus reaches a maximum after perihelion, as does the size of the coma, in agreement with several observations. The denser surface layers made either of ice or of dust may break away. An ideal, initially homogeneous and spherical nucleus cannot remain isothermal so that it must gradually develop considerable surface nonuniformities through localized phase changes, evaporation, and break-away. An explanation of the splitting of comets as far as 9 AU from the Sun is suggested in terms of heating of a CO₂-rich inclusion in a nucleus.     

Processing of cometary grains at the nucleus surface

Cometary material inevitably undergoes chemical changes before and on leaving the nucleus. In seeking to explain comets as the origin of many IDPs (interplanetary dust particles), an understanding of potential surface chemistry is vital. Grains are formed and transformed at the nucleus surface; much of the cometary volatiles may arise from the organic material. In cometary near-surface permafrost, one expects cryogenic chemistry with crystal growth and isotope. This could be the hydrous environment where IDPs form. Seasonal and geographic variations imply a range of environmental conditions and surface evolution. Interplanetary dust impacts and electrostatic forces also have roles in generating cometary dust. The absence of predicted cometary dust envelopes is compatible with the wide range of particle structures and compositions. Study of IDPs would distinguish between this model and alternatives that see comets as aggregates of core-mantle grains built in interstellar clouds.     

Probing the internal structure of the nuclei of comets

There is no direct evidence about the internal structure of cometary nuclei, which are mostly hidden by their gas and dust comae, and have not yet been orbited by any spacecraft. Their densities are low, typically of about 400 kg m⁻³ for 9P/Tempel 1 (that was impacted by the Deep Impact probe) and 67P/Churyumov-Gerasimenko (that is the target of the Rosetta mission). Such low densities are in favour of a high macro-porosity, or a high micro-porosity, or both. Observations of disruption or splitting of nuclei indeed suggest that some huge sub-nuclei or some meter-sized fragments could be the building blocks of comets. Analysis, from in-situ measurements and from remote light scattering observations, of the structure of the dust particles, which significantly consist of fluffy aggregates of submicron-sized grains, could be in favour of a fractal structure. However, the presence of huge icy grains in the innermost coma, and of flat layers on the surface of 9P/Tempel 1, are clues to the complexity of these objects, which have suffered drastic erosion phenomena on their elongated orbits. It is expected that the Rosetta mission will provide a fair understanding of the structure of the deep interior of the nucleus of 67P/Churyumov-Gerasimenko, thanks to the on-board CONSERT experiment.     

Direct statistical simulation of the near-surface layers of a cometary atmosphere. II: A nonspherical nucleus

The study presents the results of numerical simulations of mass-transfer processes in the near-surface layer of the cometary nucleus and in the inner part of the cometary atmosphere, which is formed under the action of solar radiation. The gas-kinetic model of the inner part of the cometary atmosphere surrounding a spherical nucleus (Skorov et al., 2004) is extended to the case of a nonspherical nucleus with axial symmetry. After high-resolution images of comets 19P/Borrelly and Wild 2 have been obtained by Deep Space 1 and Stardust spacecraft, such an extension seems to be vital and important. The nucleus and the inner part of the coma are closely related to each other because of the permanent exchange of energy and mass; therefore, they are modeled consistently. As in the first part of our study, the boundary conditions at the inner boundary of the simulation domain, which are necessary for gas-kinetic simulations, were determined from the self-consistent model of heat and mass transfer in a porous cometary nucleus that was developed earlier by the authors. The model took into account the volumetric character of the radiation absorption in a porous sublimating medium, the kinetic regime of the transport of sublimation products in the pores, and the backward gas fluxes from the coma due to intermolecular collisions. We considered different models of the nucleus structure that determined the effective gas production. Using the direct simulation Monte Carlo method, we computed the two-dimensional gas flow from a heterogeneous nonspherical cometary nucleus. The simulations were performed using the SMILE software. The parallel computer implementation of the software made it possible to calculate the spatial structure of the gas flow for the entire circumnucleus zone.     

On the rotational breakup of cometary nuclei and centaurs

Here we estimate the regions of stability, fragmentation, and destruction for cometary bodies versus rotational breakup in the radius-rotational period plane. By testing different plausible physical models of the cometary nucleus equation of state, we show that the plane is divided into 3 segments: the allowed, damaged, and forbidden regions. We then compare the location of well-observed comets with respect to the separation lines. The range of constituent material parameters from the literature for cometary nuclei are used to show that all the observed comets lie in the allowed region, except for Comet C/1995 O1 (Hale-Bopp), which resides in the damaged region (where the body is fractured and only held together gravitationally). We speculate that the extremely high activity demonstrated by Comet Hale-Bopp during the 1997 apparition may have been due to its highly fractured state. Comet Hyakutake, observed to emit fragments at perigee in 1996, may be near the boundary of the damaged region. Comet C/1999 S4 (LINEAR) was solidly in the rotationally allowed region, making its disintegration in July 2000 due to centrifugal forces unlikely. In contrast to the comets, the centaurs do not cluster in the allowed region, with the majority falling instead into the rotationally damaged and forbidden regions. The centaurs are only stable against breakup assuming much stronger solid water ice properties, strongly suggesting that on the whole, these bodies have different bulk physical properties than cometary nuclei.