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.     

Composition of cometary dust from polarization spectra

The wavelength dependence of the polarization (“polarization spectra”) of cometary dust is discussed. It is shown that, in the case of large phase angles, the wavelength dependence of the polarization is mainly controlled by the complex refractive index of the particle material, whereas the spectral dependence of the intensity is also sensitive to the size of the particles. This suggests that observations of “polarization spectra” may determine the composition of cometary dust. An attempt is made to find the composition of the cometary dust material by comparing the observed polarimetric data with laboratory measurements of complex refractive indices of possible cometary constituents. Silicates, graphite, metals, organics, water ice and their mixtures are considered. It is shown that astronomical silicate must be the most abundant constituent of cometary dust in the range of heliocentric distances from 0.8 to 1.8 AU, whereas the volume fraction of pure graphite or pure metals is less then 1%. A substance similar to that of F-type asteroids may be present in comets. There is evidence for an organic material that is being destroyed between heliocentric distances of 0.8–1.8 AU.     

CCD Polarimetric Imaging Of Comet Hale-Bopp (C/1995 O1)

Polarimetric observations of the light scattered by dust have been carried out at Pic-du-Midi Observatory with the 2 m telescope in June and September–October 1996, and at Haute-Provence Observatory with the 0.80 m telescope in April 1997. They cover a total number of 11 nights and a large (6.9°–47.7°) phase angle range. The spatial resolution allows to underline structures in the coma, as well in the brightness images as in the polarization maps, with a correlation between the regions of bright structures and the regions of higher polarization. A clear difference appears between the sunward and antisunward side, with higher polarization on the antisunward side. The phase angle coverage allows us to obtain a polarimetric phase curve for the whole coma and to compare it with other cometary phase curves. The degree of polarization is higher for Hale-Bopp than for the comets previously observed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Cometary dust properties retrieved from polarization observations: Application to C/1995 O1 Hale-Bopp and 1P/Halley

The analysis of the polarized light scattered by cometary dust particles provides information on the physical properties of the solid component of cometary comae for C/1995 O1 Hale-Bopp and 1P/Halley. A model of light scattering by a size distribution of aggregates of up to 256 submicron-sized grains (spherical or spheroidal) mixed with single spheroidal particles has been developed, with its parameters adjusted to fit the phase angle and wavelength dependence of the polarization observations. The particles are built of two materials: a non-absorbing silicates-type material and a more absorbing organic-type material. The model reproduces accurately the inversion angle and the positive branch of the polarization phase curves from the visible to the near-infrared spectral domains. A negative branch of the polarization phase curves appears in our model, although the negative branch is not deep enough to reproduce accurately the observations. Significant differences are shown between the two comets, with dominance of small grains in the coma of Comet C/1995 O1 Hale-Bopp, well fitted by a distribution of the volume-equivalent diameter, a, following a^−3.0 with a lower cutoff around 0.20 μm and an upper cutoff of at least 40 μm. For 1P/Halley, the size distribution follows a^−2.8 with a lower cutoff around 0.26 μm and an upper cutoff of about 38 μm. The relative amount of organic-type particles is larger for 1P/Halley while the amount of aggregates, significant for both comets, is larger for C/1995 O1 Hale-Bopp.     

A Negative Branch of Polarization for Comets and Atmosphereless Celestial Bodies and the Light Scattering by Aggregate Particles

Optical observations of comets and atmosphereless celestial bodies show that a change of sign of the linear polarization of scattered light from negative to positive at phase angles less than 20° is typical of the cometary coma, as well as of the regolith of Mercury, the Moon, planetary satellites, and asteroids. To explain a negative branch of polarization, this research suggests a unified approach to the treatment of cometary-dust particles and regolith grains as aggregate forms. A composite structure of aggregate particles resulting in the interaction of composing structural elements (monomers) in the light-scattering process is responsible for the negative polarization at small phase angles, if the monomer sizes are comparable to the wavelength. The characteristics of single scattering of light calculated for aggregates of this kind turned out to be close to the properties observed for cometary dust. Unlike the cometary coma, the regolith is an optically semi-infinite medium, where the interaction between particles is significant. To find the reflectance characteristics of regolith, the radiative-transfer equation should be solved for a regolith layer. In this case, the interaction between scatterers can be modeled to a certain extent by representing the regolith grains as aggregate structures consisting of several or many elements. Although real regolith grains are much larger than the particles considered here, laboratory measurements have shown that it is precisely the surface irregularities comparable to the wavelength that cause a negative branch of polarization. The main observed features of the phase and spectral dependence of the linear polarization of light scattered from comets and atmosphereless celestial bodies, which are due to the difference of the elementary scatterers in composition, size, and structure, can be successfully explained using the aggregate model of particles.     

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

Spectropolarimetry of comets Austin and Churyumov-Gerasimenko

Spectropolarimetric observations from 5000 to 8000 Å have been obtained for comets P/Austin (1982g) and P/Churyumov-Gerasimenko (1982f). The observations were spaced over phase angles of 50–125° for comet Austin and 10–40° for comet Churyumov-Gerasimenko. The use of spectropolarimetry allowed an evaluation of continuum polarization without molecular line contamination. Especially for comet Churyumov-Gerasimenko, the curve of polarization versus phase angle resembles curves for asteroids, where the polarization is negative (electric vector maximum parallel to the scattering plane) for phase angles less than 20° and the most negative polarization is from ?1 to ?2%. The negative polarization at backscattering angles may be due to multiple scattering in agglomerated grains, as assumed for asteroids, or to Mie scattering by small dielectric particles. If multiple scattering is important in comet dust, polarization measurements may imply a low albedo, less than 0.08. The polarization of comet Austin remained steady during a large change in the dust production rate. Both comets increased continuum flux by a factor of 2 near perihelion. The continuum of comet Churyumov-Gerasimenko had the shape of the solar spectrum with derivations less than 5%. The equivalent width of spectral features of C₂, NH₂, and O varied as r^?2.     

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.     

Polarization of Light Scattered by Solar System Bodies and the Aggregate Model of Dust Particles

The present study considers the dependence of characteristics of light scattering by aggregate particles on the refractive index, size, and number of spherical particles composing the aggregate, as well as on the structure and porosity of the cluster. The parameters were varied in sufficiently wide ranges to let a coherent picture of the polarimetric properties of relatively small aggregate particles emerge (the size parameter of the aggregate is less than 10). It was shown that, in the framework of the aggregate model, the behavior of polarization phase curves observed for both comets and regolith surfaces can be explained. The modeling carried out confirms that the sizes of the cometary dust particles are larger than the wavelength. However, the grains forming the cometary dust particles or the regolith (or details of the particle surface) have a size less than 0.3–0.5 m. This agrees with estimates obtained by other methods. The determining role in the formation of the polarization phase curve is played by the structure of the external layer of the clusters. The appearance of the negative branch of polarization and its shape substantially depend on the effectiveness of the interference of multiply scattered waves and on the interaction in the near field at these phase angles. Interference and interaction in the near field in turn are determined by the sizes of elementary scatterers and the structure of the ensemble. If the number of constituent particles in the aggregate is larger than several tens, its role in the formation of the negative branch of polarization is minor, while the influence on the polarization maximum position is rather substantial. The polarimetric data alone cannot provide a unique estimate of the refractive index: the brightness measurements must be invoked as well. For a more complete quantitative interpretation of the observations, the scattering matrix of aggregates comparable in size to or larger than the wavelength must be calculated in the short- and long-wavelength ranges, which still encounters serious theoretical and technical difficulties. Moreover, in order to obtain unique results, it is obvious that the spectral range of observations must be extended and that other types of measurements, such as spectroscopic ones, must also be used.     

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

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

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.     

Polarimetric study of levitating dust aggregates with the PROGRA experiment

In order to interpret polarimetric remote observations of solar system dust clouds (e.g. cometary coma dust), laboratory measurements are needed. Three samples composed of aggregates are studied: crystallized enstatite, pyrogenic alumina and titanium oxide. The new version of the PROGRA² instrument allows to obtain polarimetric images of the samples under levitation. The dependence of polarization with phase angle and particle size is studied, as well as the effect of the porosity of the particles. Values of polarization at small phase angles are also discussed. The polarization near 90° decreases when the agglomerate size increases and when the porosity increases.     

Aperture Polarimetry And Photometry Of Comet Hale-Bopp

We present results of polarimetric and photometric observations of bright comet C/1995 O1 (Hale-Bopp) obtained at the 0.7 m telescope of Kharkov University Observatory from June 18, 1996 to April 24, 1997. The IHW and HB comet filters were used. The C₂ and C₃ production rates for Hale-Bopp are more than one order of magnitude larger and the dust production rates are more than two orders of magnitude larger than the Halley ones at comparable distances. Hence, Hale-Bopp was one of the most dusty comets. The average UC-BC and BC-RC colours of the dust were −0.02 and 0.13 mag, respectively. The polarization of comet Hale-Bopp at small phase angles of 4.8–13.0° was in good agreement with the date for comet P1/Halley at the same phase angles in spite of the fact that the heliocentric distances of comments differed nearly twice. However, at intermediate phase angles of 34–49° the polarization of comet Hale-Bopp was significantly larger than the polarization of the other dusty comets. It is the first case of such a large difference found in the continuum polarization of comets. The wavelength dependence of polarization for Hale-Bopp was steeper than for other dusty comets. The observed degree of polarization for the anti-sunward side of the coma was permanently higher than that for the sunward shell side. The polarization phase dependence of Hale-Bopp is discussed and compared with the polarization curves for other dusty comets. The peculiar polarimetric properties of comet Hale-Bopp are most likely caused by an over-abundance of small or/and absorbing dust particles in the coma. This revised version was published online in July 2006 with corrections to the Cover Date.     

A Monte Carlo Code to Compute Energy Fluxes in Cometary Nuclei

A Monte Carlo model designed to compute both the input and output radiation fields from spherical-shell cometary atmospheres has been developed. The code is an improved version of that by H. Salo (1988, Icarus76, 253-269); it includes the computation of the full Stokes vector and can compute both the input fluxes impinging on the nucleus surface and the output radiation. This will have specific applications for the near-nucleus photometry, polarimetry, and imaging data collection planned in the near future from space probes. After carrying out some validation tests of the code, we consider here the effects of including the full 4×4 scattering matrix in the calculations of the radiative flux impinging on cometary nuclei. As input to the code we used realistic phase matrices derived by fitting the observed behavior of the linear polarization as a function of phase angle. The observed single scattering linear polarization phase curves of comets are fairly well represented by a mixture of magnesium-rich olivine particles and small carbonaceous particles. The input matrix of the code is thus given by the phase matrix for olivine as obtained in the laboratory plus a variable scattering fraction phase matrix for absorbing carbonaceous particles. These fractions are 3.5% for Comet Halley and 6% for Comet Hale-Bopp, the comet with the highest percentage of all those observed.The errors in the total input flux impinging on the nucleus surface caused by neglecting polarization are found to be within 10% for the full range of solar zenith angles. Additional tests on the resulting linear polarization of the light emerging from cometary nuclei in near-nucleus observation conditions at a variety of coma optical thicknesses show that the polarization phase curves do not experience any significant changes for optical thicknesses τ?0.25 and Halley-like surface albedo, except near 90° phase angle.     

Physical properties of cometary and interplanetary dust

In the absence of numerous in situ studies, physical properties of cosmic dust may be derived from observations of their light scattering and thermal properties, through numerical simulations making use of realistic assumptions. Estimations about cometary and interplanetary dust composition, structure, size, as well as about their light scattering and thermal properties, are first summarized. We then present and discuss the numerical simulations we have performed with different types of particles: core-mantle submicron-sized elongated grains (having contributed to the formation of cometary dust), fractal aggregates of such grains (found in cometary comae and in the interplanetary dust cloud), and fractal aggregates of large dust grains (found in cometary dust trails).A very satisfactory fit to the numerous polarimetric observations of comet Hale-Bopp is obtained for a mixture with about 33–60% of organics in mass, with a power law size distribution with an index of (−3) and a radius of 20 μm for the upper cut-off. For the less-constrained polarimetric observations of interplanetary dust near 1 AU, a fit is obtained for a mixture with about 40% of organics in mass, with a similar size distribution and a radius of about 50 μm for the upper cut-off. The ensemble of results obtained for the interplanetary dust strongly suggest that its light scattering and thermal properties stem from the presence of compact and fluffy particles, with compositions ranging from silicates to more absorbing materials, whose contribution decreases with decreasing distance to the Sun.     

Spectropolarimetry of the Deep Impact target Comet 9P/Tempel 1 with HiVIS

High resolution spectropolarimetry of the Deep Impact target, Comet 9P/Tempel 1, was performed during the impact event on July 4th, 2005 with the HiVIS spectropolarimeter and the AEOS 3.67-m telescope on Haleakala, Maui. We observed atypical polarization spectra that changed significantly in the few hours after the impact. The polarization of scattered light as a function of wavelength is very sensitive to the size and composition (complex refractive index) of the scattering particles as well as the scattering geometry. As opposed to most observations of cometary dust, which show an increase in the linear polarization with the wavelength (at least in the visible domain and for phase angles greater than about 30, a red polarization spectrum) observations of 9P/Tempel 1 at a phase angle of 41° beginning 8 min after impact and centered at 6:30 UT showed a polarization of 4% at 650 nm falling to 3% at 950 nm. The next observation, centered an hour later showed a polarization of 7% at 650 nm falling to 2% at 950 nm. This corresponds to a spectropolarimetric gradient, or slope, of −0.9% per 1000 Å 40 min after impact, decreasing to a slope of −2.3% per 1000 Å an hour and a half after impact. This is an atypical blue polarization slope, which became more blue 1 h after impact. The polarization values of 4 and 7% at 650 nm are typical for comets at this scattering angle, whereas the low polarization of 2 and 3% at 950 nm is not. We compare observations of Comet 9P/Tempel 1 to that of a typical comet, C/2004 Machholz, at a phase angle of 30° which showed a typical red slope, rising from 2% at 650 nm to 3% at 950 nm in two different observations (+1.0 and +0.9% per 1000 Å).     

The Dust in Comet C/1999 S4 (LINEAR) during Its Disintegration: Narrow-Band Images, Color Maps, and Dynamical Models

Comet C/1999 S4 was observed with the 2m-telescopes of the Bulgarian National Observatory and Pik Terskol Observatory, Northern Caucasus, Russia, at the time of its disintegration. Maps of the dust brightness and color were constructed from images obtained in red and blue continuum windows, free from cometary molecular emissions. We analyze the dust environment of Comet C/1999 S4 (LINEAR) taking into account the observed changes apparent in the brightness images and in plots of Afρ profiles as function of the projected distance ρ from the nucleus. We also make use of the syndyne-synchrone formalism and of a Monte Carlo model based on the Finson-Probstein theory of dusty comets. The brightness and color of individual dust particles, which is needed to derive theoretical brightness and color maps of the cometary dust coma from the Monte Carlo model, is determined from calculations of the light scattering properties of randomly oriented oblate spheroids. In general, the dust of Comet C/1999 S4 (LINEAR) is strongly reddened, with reddening values up to 30%/1000 Å in some locations. Often the reddening is higher in envelopes further away from the nucleus. We observed two outbursts of the comet with brightness peaks on July 14 and just before July 24, 2000, when the final disintegration of the comet started. During both outbursts an excess of small particles was released. Shortly after both outbursts the dust coma “turns blue.” After the first outburst, the whole coma was affected; after the second one only a narrow band of reduced color close to the tail axis was formed. This difference is explained by different terminal ejection speeds, which were much lower than normal in case of the second outburst. In particular in the second, final outburst the excess small particles could originate from fragmentation of “fresh” larger particles.     

Radio observations of Jupiter-family comets

Radio observations from decimetric to submillimetric wavelengths are now a basic tool for the investigation of comets. Spectroscopic observations allow us: (i) to monitor the gas production rate of the comets, by directly observing the water molecule, or by observing secondary products (e.g., the OH radical) or minor species (e.g., HCN); (ii) to investigate the chemical composition of comets; (iii) to probe the physical conditions of cometary atmospheres: kinetic temperature and expansion velocity. Continuum observations probe large-size dust particles and (for the largest objects) cometary nuclei.Comets are classified from their orbital characteristics into two separate classes: (i) nearly isotropic, mainly long-period comets and (ii) ecliptic, short-period comets, the so-called Jupiter-family comets (JFCs). These two classes apparently come from two different reservoirs, respectively, the Oort cloud and the trans-Neptunian scattered disc. Due to their different history and—possibly—their different origin, they may have different chemical and physical properties that are worth being investigated.The present article reviews the contribution of radio observations to our knowledge of the JFCs. The difficulty of such a study is the commonly low gas and dust productions of these comets. Long-period, nearly isotropic comets from the Oort cloud are better known from Earth-based observations. On the other hand, JFCs are more easily accessed by space missions. However, unique opportunities to observe JFCs are offered when these objects come by chance close to the Earth (like 73P/Schwassmann-Wachmann 3 in 2006), or when they exhibit unexpected outbursts (as did 17P/Holmes in 2007).About a dozen JFCs were successfully observed by radio techniques up to now. Four to ten molecules were detected in five of them. No obvious evidence for different properties between JFCs and other families of comets is found, as far as radio observations are concerned.     

Modeling of meteoroid streams: The velocity of ejection of meteoroids from comets (a review)

An analytical review of the models of ejection of meteoroids from cometary nuclei is presented. Different formulas for the ejection velocity of meteoroids and the corresponding parameters are discussed and compared with the use of comet Halley and the Geminids meteoroid stream as examples. The ejection velocities obtained from observations of the dust trails of comets are discussed, and the values for comets 2P/Encke, 4P/Faye, 17P/Holmes, 22P/Kopff, and 67P/Churyumov-Gerasimenko are compared to the velocities yielded by Whipple’s model. The uncertainty intervals of the results are estimated.