Coma morphology of Jupiter-family comets

We present a general review of cometary coma morphology, with specific reference to how it used in studies of Jupiter-family comets. We introduce the most common features that are seen in gas and dust observations, and summarize of how they are used to infer the properties of the nucleus and coma. We also expand the discussion to cover other topics relating to morphology, including the general shape of the coma (characterized by radial gradient profiles) and spatial maps of the color, albedo and polarization of the dust. We address the pros and cons of the different approaches used in the interpretation and analysis of the features. Finally, we review the results obtained for specific comets and compare the Jupiter-family comets to those from other classes.     

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

Present and Future Cometary Science with the IRAM Plateau de Bure Interferometer

Interferometric observations are essential to probe the molecular emission in the inner cometary atmospheres and study the outgassing from the nucleus. Mapping the continuum emission can provide information about the dust and/or nucleus properties. We present here a summary of the observations of the dust and gas coma of comet 17P/Holmes and nuclear observations of 8P/Tuttle, both carried out with the IRAM interferometer at Plateau de Bure (PdBI) in 2007–2008. The observations of these two comets demonstrate the ability of the PdBI in terms of cometary science. In the near future, several improvements will be made (new receivers at 0.8 mm, a new wide-band correlator) allowing more frequent and more detailed studies of comets. On the long term, NOEMA, an expansion project, may add up to six antennas to the Plateau de Bure Interferometer, and increase the baseline lengths. Such an instrument would offer a complement to ALMA to track comets of the northern hemisphere with about half the sensitivity of ALMA for continuum studies.     

A comparison of the continuum spectra of four comets

The continuum spectra of comets carry information concerning the physical and chemical properties of solid coma grains. Although it is not feasible to use the continuum spectra to uniquely characterize the solid grains, variations among the continua of different comets may reveal subtle differences in their respective grain populations. We have taken and reduced optical spectra of four comets in the wavelength range 3700–7300 Å using a single observing system and reduction procedure. The continua all appear reddened with respect to the solar spectrum. The amount of reddening is consistent with a prevalence of ～2-μm-sized grains in all four comets, if the refractive indices of the grains are approximately equal to those of terrestrial rocks. Significant color differences were measured among the comets. Different intrinsic grain properties are suggested since the scattering geometries were very similar. The amount of reddening does not appear to be correlated with the amount of dust in the coma.     

Composition of Comets: Observations and Models

We analyze the chemical composition and abundances of comets based on in situ measurements of Comet 1P/Halley and remote sensing observations of several recent bright comets including Hale-Bopp (C/1995 O1) and Hyakutake (C/1996 B2), in light of the elemental abundances of the solar system. Nitrogen is underabundant in comets relative to the solar system because nitrogen tends to be in N₂, which is chemically relatively inert. While many details remain uncertain, some gross features are emerging. The abundance of water : silicates: carbonaceous molecules (CO, CO₂, and hydrocarbons) by mass is approximately 1 : 1 : 1. Furthermore, the mass abundance of ice : dust (silicates and hydrocarbon polycondensates) is about1 : 1. We compare a list of identified comet molecules with molecules detected in the interstellar medium, although a comparison with their relative abundances, particularly in the ice phase, would be more meaningful. However, ice-phase abundances are not yet available. One can expect a variation of the abundances of carbon-bearing molecules in comets to be associated with their place of origin in the solar nebula. However, we also note that comets are heterogeneous. Thus, observed differences may be related to the place of origin, heterogeneity of the nucleus, or acquired through evolution. The molecular and elemental compositions of the coma are most likely not the same as those in the nucleus. This is particularly true for volatile ices and their gases and for the dust-to-ice and dust-to-gas ratios. Analyses must carefully consider the three sources of gas: Water from the surface of the nucleus, gases more volatile than water from the interior of the nucleus, and gases from the sublimation of the dust distributed in the coma. Topography on the surface of the nucleus may cause important evolutionary differences in the dust-to-gas mass ratio. Relatively inactive areas on the surface of the nucleus are probably associated with convex topography. Gas sublimated from convex areas (hills and mountains) diverges more strongly relative to gas sublimated from concave areas, which can entrain dust more efficiently. Thus, the entrainment of dust from convex areas is poor and dust may fall back to the surface of the nucleus creating a dust mantle, which further inhibits outgassing.     

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.     

Dust coma of Comet C/1999 S4 (LINEAR): imaging polarimetry during nucleus disruption

Imaging and polarimetry of Comet C/1999 S4 (LINEAR) during its disruption provide information about the physical properties of the scattering dust particles, and some insight into the structure of the nucleus. A significant decrease in the brightness was noticed, together with a drastic change in the shape of the dust coma. The whole-coma polarization increased, which was typical of a comet with a near 27 percent maximum in polarization, the increase being comparable to previous observations for comets suffering a limited fragmentation. An important gradient in the intensity on the solar side corresponds to the ejected material. The degree of polarization in this region is higher than generally observed in jets and it increases with time as the nucleus (or its fragments) breaks up and ejects relatively large and compact particles. In the surrounding coma, these large particles are fragmented on short time-scales, indicated by the decrease of polarization. These results suggest that the fragile nucleus was not, as far as the physical properties of the dust are concerned, differentiated, and that it was possibly built of primordial cometesimals originating from the same formation region.     

Hale-Bopp: What Makes a Big Comet Different? Coma Dynamics: Observations and Theory

Comet Hale-Bopp was the largest comet by almost any definition, observed at least since the advent of modern observing techniques. In a more typical comet both the chemical and dynamical influences of collisional processes are limited by the short time a parcel of gas sublimated from the nucleus remains in the dense part of the coma. The resulting large size of the collisional coma in comet Hale-Bopp had important consequences on the dynamics of the coma, which in turn has important consequences on how observations are interpreted with standard models. Measured velocities of typical gas species (mostly the observed radicals) as well as dust were larger than normal comets. Conversely, velocities of super thermal atomic hydrogen were smaller than normal because of the samecollisional processes. Furthermore, as a consequence, dust particles, which are dragged by the outflowing gas, were also accelerated to larger velocities. Such larger velocities are not simply an interesting curiosity in their own right, because nearly all observations of dust and gas are interpreted with models of the coma that depend directly on some measurement or assumption with regard to velocity. In this presentation both observations and theory regarding the dynamical conditions in the coma of comet Hale-Bopp are summarized.     

Coma Morphology of Three Non-periodic Comets

The coma morphology and short-term evolution was investigated of three non-periodic comets in retrograde orbits, C/2001 Q4 (NEAT), C/2002 T7 (LINEAR), and C/2003 K4 (LINEAR). All three comets display distinct coma features, which were very different from one comet to the next and remained rather constant in shape during the observational period. A single, broad feature perpendicular to the sun-tail direction dominated the coma of C/2003 K4 in all used filters (B,V,R,I), whereas the coma of Comet C/2002 T7 exhibited different features in blue and red filters. C/2001 Q4 showed rather complex coma morphology with clear short-term variability in coma brightness. Therefore, these non-periodic comets neither show a featureless coma nor any similarities of the features detected. The overall distribution of coma material was investigated from the shape of radial coma profiles averaged around the comet nucleus. For C/2001 Q4 and C/2002 T7, the slopes fitted to the linear part of these profiles are flatter in the blue than in the red, which can be explained by the presence of coma gas. For C/2003 K4 no such difference is indicated in the May observations (r = 2.3 AU), while in July (r = 1.7 AU) the profiles in the B-filter are flatter than in V, R, and I, hence gas contamination was relevant at least in the B filter. The R and I filter images were used to determine approximate Afρ values of each comet as a function of time.     

Understanding measured water rotational temperatures and column densities in the very innermost coma of Comet 73P/Schwassmann–Wachmann 3 B

Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J., 685, 659?677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J., 732) to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann–Wachmann 3 obtained with a very high spatial resolution of ～30 km using IRCS/Subaru in May 2006 (Bonev, B.P., Mumma, M.J., Kawakita, H., Kobayashi, H., Villanueva, G.L. [2008]. Icarus, 196, 241?248). For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10–20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.     

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.     

Nitrogen Sulfide In Comets Hyakutake (C/1996 B2) And Hale–Bopp (C/1995 O1)

The chemistry of both nitrogen and sulfur presents interesting problems in comets.In this paper, we use a model of cometary comae with gas-phase chemical kineticsand gas dynamics to predict molecular abundances in the inner coma region for twoof the brightest comets in the past 20 years, Hyakutake (C/1996 B2) and Hale–Bopp(C/1995 O1). In this progress report we concentrate on the gas-phase chemistry of thenitrogen sulfide (NS) radical at a heliocentric distance of 1 AU to study the abundanceof NS using a detailed photo and chemical reaction network with over 100 species andabout 1000 reactions. The results are compared with recent observations of CometHale–Bopp and reveal that conventional gas-phase reactions schemes do not produceNS in sufficient quantities to explain the observations. We plan to continue therefinement of the model to improve agreement with observational constraints.     

Thermal modeling of cometary nuclei

A new model of the sublimation of volatile ices from a cometary nucleus has been developed which includes the effects of diurnal heating and cooling, rotation period and pole orientation, and thermal properties of the ice and subsurface layers. The model also includes the contribution from coma opacity, scattering, and thermal emission, where the properties of the coma are derived from the integrated rate of volatile production by the nucleus. The model is applied to the specific case of the 1986 apparition of Halley's comet. It is found that the generation of a cometary dust coma actually increases the total energy reaching the Halley nucleus. This results because of the significantly greater geometrical cross section of the coma as compared with the bare nucleus, and because the coma provides an essentially isotropic source of multiply scattered sunlight and thermal emission over the entire nucleus surface. For Halley, the calculated coma opacity is approximately 0.2 at 1 AU from the Sun, and 1.2 at perihelion (0.587 AU). At 1 AU this has little effect on dayside temperatures (maximum ≈200°K) but raises nightside temperatures (minimum ≈150°K) by about 40°K. At perihelion the higher opacity results in a nearly isothermal nucleus with only small diurnal and latitudinal temperature variations. The general surface temperature is 205°K with a maximum of 209°K at local noon on the equator. Some possible consequences of the results with respect to the generation of nongravitational forces, observed volatile production rates for comets, and cometary lifetimes against sublimation are discussed.     

Chemical and physical properties of gas jets in comets: I. Monte Carlo model of an inner cometary coma

We describe a 3-dimensional, time-dependent Monte Carlo model developed to analyze the chemical and physical nature of a cometary gas coma. Our model includes the necessary physics and chemistry to recreate the conditions applicable to Comet Hale-Bopp when the comet was near 1 AU from the Sun. Two base models were designed and are described here. The first is an isotropic model that emits particles (parents of the observed gases) from the entire nucleus; the second is a jet model that ejects parent particles solely from discrete active areas on the surface of the comet nucleus, resulting in coma jets. The two models are combined to produce the final model, which is compared with observations. The physical processes incorporated in both base models include: (1) isotropic ejection of daughter molecules (the observed gases) in the parent's frame of reference, (2) solar radiation pressure, (3) solar insolation effects, (4) collisions of daughter products with other molecules in the coma, and (5) acceleration of the gas in the coma. The observed daughter molecules are produced when a parent decays, which is represented by either an exponential decay distribution (photodissociation of the parent gas) or a triangular distribution (production from a grain extended source). Application of this model to the analysis the OH, C₂ and CN gas jets observed in the coma of Comet Hale-Bopp is the focus of the accompanying paper [Lederer, S.M., Campins, H., Osip, D.J., 2008. Icarus, in press (this issue)].     

Is Hale–Bopp A ‘Juvenile’ Comet – Study Of The Energy Balance To Explain The Vapor Flux Of Volatiles From The Surface

Sublimation of minor gases from ices inside of a porous comet nucleus strongly depends on the effective energy input. Our model meant to describe the gas flux inside and out of the porous nucleus has been used to study the influence of physical and structural parameters on the effective energy input. We solve the conservation equations for H₂O and CO as the most abundant minor component of higher volatility under appropriate boundary conditions. From the calculations we obtain the gas flux from volatile, icy components inside the porous nucleus, temperature profiles, changes in relative chemical abundances, and the gas flux into the coma for each of the volatiles. We will show results from our calculations for a model comet in the orbit of Hale-Bopp (C/1995 O1). In this paper we focus on the energy balance at the surface. We will also relate measurements of molecule fluxes to available energies and try to provide hints about the evolutionary status of the comet. This revised version was published online in July 2006 with corrections to the Cover Date.     

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.     

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.     

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.     

Rydberg Matter as the dust atmosphere in comets: Spectroscopic and polarization signatures

The long-lived excited state of matter called Rydberg Matter (RM) may explain several spectroscopic features in space, like the diffuse interstellar bands (DIBs) and the unidentified infrared bands (UIR, UIB). RM is here used to interpret some previously unexplained or inconsistent features in comets: (1) line absorption in the emission from the nucleus, (2) IR band emission from the coma, (3) selective and variable molecular line emission from the coma, (4) extended sources of molecules in the coma, (5) degree of linear polarization of light scattered from the coma. (1) The unexplained IR absorptions observed in the flyby of the nucleus of the Comet Borrelly agree well with RM emission bands observed by stimulated emission in the laboratory. It is proposed that RM is the so-called ultrared matter or at least formed from it. (2) The IR bands previously attributed to silicate particles are shown to be better described by RM theory. Transitions in atoms in RM are shown to dominate. (3) The inverted RM medium will optically amplify light from molecular transitions in the comet comae, in agreement with observations that many of the molecular IR emission lines lie within the emission bands from RM, or much too close to Rydberg transitions to be accidental. (4) The unexplained extended sources observed, e.g., for CO are proposed to be due to release of molecules previously incorporated in the RM clusters at low temperature. Such clusters may be the very small particles observed by space probes. (5) Finally, the largely unexplained measurements of the degree of linear polarization of scattered sunlight from comets are explained as due to scattering by the planar RM clusters. Quantitative agreement is demonstrated.