Negative ion chemistry in the coma of comet 1P/Halley

Negative ions (anions) were identified in the coma of comet 1P/Halley during in situ Electron Electrostatic Analyzer measurements performed by the Giotto spacecraft in 1986. These anions were detected with masses in the range 7–110 amu, but with insufficient mass resolution to permit unambiguous identification. We present details of a new chemical‐hydrodynamic model for the coma of comet Halley that includes—for the first time—atomic and molecular anions, in addition to a comprehensive hydrocarbon chemistry. Anion number densities are calculated as a function of radius in the coma, and compared with the Giotto results. Important anion production mechanisms are found to include radiative electron attachment, polar photodissociation, dissociative electron attachment, and proton transfer. The polyyne anions C₄H^? and C₆H^? are found to be likely candidates to explain the Giotto anion mass spectrum in the range 49–73 amu. The CN^? anion probably makes a significant contribution to the mass spectrum at 26 amu. Larger carbon‐chain anions such as C₈H^? can explain the peak near 100 amu provided there is a source of large carbon‐chain‐bearing molecules from the cometary nucleus.     

The Near-Nuclear Coma of Comet Halley in March 1986

The cameras carried onboard the flyby missions to comet P/Halleyin 1986 imaged the near nuclear jet activity fromseveral spatial directions. The observed, very structured near nucleardust jets were considered at that timeas the result of dust emission from well localized active surface regions(without supporting 3-D model computations, however).Based on the first, recently developed 3-D gas dynamical model ofP/Halley's activity,we have been shown that jet features can be reproduced assuming ahomogeneous dusty icenucleus surface. The dust in the collisional near nuclear comais concentrated along the gas flow discontinuities resulting from thecomplicated surface orography, creating the visual impression ofdust jets. We present here the results of these calculations forthe near nucleus dust distributions,and we compare them with the direct observations made during thethree Halley flybys (Vega 1, Vega 2, and Giotto).     

Evidence for methane and ammonia in the coma of comet P/Halley

Methane and ammonia abundances in the coma of Halley are derived from Giotto IMS data using an Eulerian model of chemical and physical processes inside the contact surface to simulate Giotto HIS ion mass spectral data for mass-to-charge ratios (m/q) from 15 to 19. The ratio m/q = 19/18 as a function of distance from the nucleus is not reproduced by a model for a pure water coma. It is necessary to include the presence of NH3, and uniquely NH3, in coma gases in order to explain the data. A ratio of production rates Q(NH3)/Q(H2O) = 0.01-0.02 results in model values approximating the Giotto data. Methane is identified as the most probable source of the distinct peak at m/q = 15. The observations are fit best with Q(CH4)/Q(H2O) = 0.02. The chemical composition of the comet nucleus implied by these production rate ratios is unlike that of the outer planets. On the other hand, there are also significant differences from observations of gas phase interstellar material.     

The dust distribution within the inner coma of comet P/Halley 1982i: encounter by Giotto's impact detectors

Analysis of the data from Giotto's Dust Impact Detection System experiment (DIDSY) is presented. These data represent measurement of the size of dust grains incident on the Giotto dust shield along its trajectory through the coma of comet P/Halley on 1986 March 13/14. First detection occurred at some 287000 km distance from the nucleus on the inbound leg; the majority of the DIDSY subsystems remained operational after closest approach (604 km) yielding the last detection at about 202000 km from the nucleus. In order to improve the data coverage (and especially for the smallest grains, to approximately 10(-19) kg particle mass), data from the PIA instrument has been combined with DIDSY data. Flux profiles are presented for the various mass channels showing, to a first approximation, a 1/R2 flux dependence, where R is the distance of the detection point from the cometary nucleus, although significant differences are noted. Deviations from this dependence are observed, particularly close to the nucleus. From the flux profiles, mass and geometrical area distributions for the dust grains are derived for the trajectory through the coma. Groundbased CCD imaging of the dust continuum in the inner coma at the time of encounter is also used to derive the area of grains intercepted by Giotto. The results are consistent with the area functions derived by Giotto data and the low albedo of the grains deduced from infrared emission. For the close encounter period (-5 min to +5 min), the cumulative mass distribution function has been investigated, initially in 20 second periods; there is strong evidence from the data for a steepening of the index of the mass distribution for masses greater than 10(-13) kg during passage through dust jets which is not within the error limits of statistical uncertainty. The fluences for dust grains along the entire trajectory is calculated; it is found that extrapolation of the spectrum determined at intermediate masses (cumulative mass index alpha = 0.85) is not able to account for the spacecraft deceleration as observed by the Giotto Radio Science Experiment and by ESOC tracking operations. Data at large masses (>10(-8) kg) recently analysed from the DIDSY data set show clear evidence of a decrease in the mass distribution index at these masses within the coma, and it is shown that such a value of the mass index can provide sufficient mass for consistency with the observed deceleration. The total particulate mass output from the nucleus of comet P/Halley at the time of encounter would be dependent on the maximum mass emitted if this change in slope observed in the coma were also applicable to the emission from the nucleus; this matter is discussed in the text. The flux time profiles have been converted through a simple approach to modeling of the particle trajectories to yield an indication of nucleus surface activity. There is indication of an enhancement in flux at t approximately -29 s corresponding to crossing of the dawn terminator, but the flux detected prior to crossing of the dawn terminator is shown to be higher than predicted by simple modelling. Further enhancements corresponding to jet activity are detected around +190 s and +270 s.     

Solar wind interactions with Comet 19P/Borrelly

The Plasma Experiment for Planetary Exploration (PEPE) made detailed observations of the plasma environment of Comet 19P/Borrelly during the Deep Space 1 (DS1) flyby on September 22, 2001. Several distinct regions and boundaries have been identified on both inbound and outbound trajectories, including an upstream region of decelerated solar wind plasma and cometary ion pickup, the cometary bow shock, a sheath of heated and mixed solar wind and cometary ions, and a collisional inner coma dominated by cometary ions. All of these features were significantly offset to the north of the nucleus-Sun line, suggesting that the coma itself produces this offset, possibly because of well-collimated large dayside jets directed 8°-10° northward from the nucleus as observed by the DS1 MICAS camera. The maximum observed ion density was 1640 ion/cm³ at a distance of 2650 km from the nucleus while the flow speed dropped from 360 km/s in the solar wind to 8 km/s at closest approach. Preliminary analysis of PEPE mass spectra suggest that the ratio of CO⁺/H₂O⁺ is lower than that observed with Giotto at 1P/Halley.     

Challenging a Paradigm: Do We Need Active and Inactive Areas to Account for Near-Nuclear Jet Activity?

We briefly describe an advanced 3D gas dynamical model developed for the simulation of theenvironment of active cometary nuclei. The model canhandle realistic nucleus shapes and alternative physical models for the gas and dust production mechanism.The inner gas coma structure is computed by solving self-consistently(a) near to the surface the Boltzman Equation(b) outside of it, Euler or Navier-Stokes equations.The dust distribution is computed from multifluid ``zero-temperature' Euler equations,extrapolated with the help of a Keplerian fountain model.The evolution of the coma during the nucleus orbital and spin motion,is computed as a succession of quasi-steady solutions. Earlier versions of the model using simple,``paedagogic' nuclei have demonstrated that the surface orographyand the surface inhomogeneity contribute similarly to structuring the near-nucleusgas and dust coma,casting a shadow on the automatic attribution of such structures to ``active areas'.The model was recently applied to comet P/Halley, for whichthe nucleus shape is available. In the companion paper of this volume,we show that most near-nucleus dust structuresobserved during the 1986 Halley flybys are reproduced, assuming that the nucleus is strictly homogeneous. Here, we investigate the effect of shape perturbations and homogeneityperturbations. We show that the near nucleus gas coma structure is robust vis-a-vissuch effects. In particular, a random distribution of active and inactive areaswould not affect considerably this structure, suggesting that such areas,even if present, could not be easily identified on images of the coma.     

Infrared emission from P/Halley's dust coma during March 1986

2 to 20 micrometers photometry of the inner dust coma of comet Halley was obtained at the NASA IRTF on Mar 6.85, 12.8, 13.75, 17.7, and 24.8. Positions offset 10" were measured as well as the central brightness. The strength of the 10 micrometers emission feature was observed to vary with location in the coma. The infrared emission is in general agreement with the dust size distribution measured from the Vega and Giotto spacecraft. Mar 6.8, 17.7, and 24.8 corresponded to strong dust jet activity. The strength of the 10 micrometers silicate emission is shown to be a sensitive indicator of grain size and thus of jet activity. Dust production rate on March 13.75, 6 h before Giotto encounter, was approximately 10(7) gm s-1.     

Analysis of the morphological structures of comet 1P/Halley in the perihelion passages in 1910 and 1986

This work is based on a systematic analysis of images of comet 1P/Halley collected during its penultimate and ultimate approaches, i.e., in 1910 and 1986. This research has identified, characterized, classified, and compared tail structures of comet 1P/Halley, namely disconnection events (DEs), wavy structures, and solitons. The images of the comet during its 1910 passage, as illustrated in the Atlas of Comet Halley 1910 II (Donn et al. 1986), were compared with those of its approach in 1986 as illustrated in The International Halley Watch Atlas of Large‐Scale Phenomena (Brandt et al. 1992). Two onsets of DEs were discovered after the perihelion passage in 1910 with an average value of the corrected cometocentric velocity (Vc) of 57 ± 15 km s^?1. Ten onsets of DEs were discovered after the perihelion passage in 1986 with an average Vc equal to 130 ± 37 km s^?1. The mean value of the corrected wavelength λc of wavy structures in 1910 is equal to 1.7 ± 0.1 × 10⁶ km, as compared to 2.2 ± 0.2 × 10⁶ km in 1986. The mean value of the amplitude A of the wave in 1910 is equal to 1.4 ± 0.1 × 10⁵ km and 2.8 ± 0.5 × 10⁵ km in 1986. The goals of this research were to report the results obtained from the analysis of the P/Halley's images from 1910 and 1986, to provide empirical data for comparison, and to form the input for future physical/theoretical work.     

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.     

Photoelectric photometry of comets in the system of standard IHW filters and the special case of comet P/Halley

Magnitudes of comets P/Giacobini-Zinner (1984e), P/Halley (1982i), P/Hartley-Good (1985 1), and Thiele (1985 m) in the bandpasses of the standard IHW comet filters are presented. For comet P/Halley production rates for CN, C₃, C₂, and solids were derived. For the gaseous components these show a strong dependence on heliocentric distance. The dependence is less steep for the solids which may be due to relatively pronounced backscattering properties in case of comet P/Halley. During one night (1985 Dec. 22/23) intensity profiles along three sections through the coma of comet P/Halley were measured. Compared with theoretical profiles they show a global anisotropy of the coma and possibly local structure.     

On the visibility of nuclei of dusty comets

The problem of visibility of a cometary nucleus discussed in general terms for single scattering by dust grains. The ratio of radiatio scattered in the dust column above the surface and that reflected from the nucleus determines the visibility of features on the nuclear surface. A contrast parameter characterizing the ration of radiation foming from the nuclear surface and that of the nuclear vicinity describes the visibility of the full nucleus against the dust fore- and background. These quantities and the intensity distribution of scattered solar radiation across the nucleus and its vicinity are calculated for the case of comet Halley at a heliocentric distance of 0.9 AU after perihelion (Giotto encounter). The scattering calculations are based on an isotropic dust distribution derived from hydrodynamics gas-dust interactions resulting in a steep densiity increase right above the surface. For Newburn's nominal model of comet Halley, an optical depth of about 0.5 impairs the visibility of the nucleus somewhat.     

Influence of the magnetic field on the density distribution of solar wind protons and cometary ions in the shock layer ahead of cometary ionospheres

The “mass loading” of the solar wind by cometary ions produced by the photoionization of neutral molecules outflowing from the cometary nucleus plays a major role in the interaction of the solar wind with cometary atmospheres. In particular, this process leads to a decrease in the solar wind velocity with a transition from supersonic velocities to subsonic ones through the bow shock. The so-called single-fluid approximation, in which the interacting plasma flows are considered as a single fluid, is commonly used in modeling such an interaction. However, it is occasionally necessary to know the distribution of parameters for the components of the interacting plasma flows. For example, when the flow of the cometary dust component in the interplanetary magnetic field is considered, the dust particle charge, which depends significantly on the composition of the surrounding plasma, needs to be known. In this paper, within the framework of a three-dimensional magnetohydrodynamic model of the solar wind flow around cometary ionospheres, we have managed to separately obtain the density distributions of solar wind protons and cometary ions between the bow shock and the cometary ionopause (in the shock layer). The influence of the interplanetary magnetic field on the position of the point of intersection between the densities with the formation of a region near the ionopause where the proton density is essentially negligible compared to the density of cometary ions is investigated. Such a region was experimentally detected by the Vega-2 spacecraft when investigating Comet Halley in March 1986. The results of the model considered below are compared with some experimental data obtained by the Giotto spacecraft under the conditions of flow around Comets Halley and Grigg–Skjellerup in 1986 and 1992, respectively. Unfortunately, our results of calculations on Comet Churyumov–Gerasimenko are only predictive in character, because the trajectory of the Rosetta spacecraft, which manoeuvred near its surface for several months, is complex.     

H2O + ions in comets: models and observations

An improved magnetohydrodynamic (MHD) model with chemistry is presented. The analysis of the source and sink terms for H₂O ⁺ shows that for small comets up to 11% of water molecules are finally ionized. For large comets (such as Halley) this fraction decreases to less than 3%. From the MHD scaling laws a similarity law for the individual ion densities is deduced which takes into account that the mother molecules are depleted by dissociation. This is applied to H₂O ⁺ ions. Radial density profiles from model calculations, observations by Giotto near comet Halley, and ground based observations of three comets confirm this scaling law for H₂O ⁺ ions. From the similarity law for the density a scaling law for the column density is derived which is more convenient to apply for ground based observations. From these scaling laws methods are derived which allow the determination of the water production rate from the ground based images of the H₂O ⁺ ions. Finally, the two dimensional images of model column densities are compared with observations.     

BVR broadband photometry of comets 1P/Halley and 4P/Faye

We analyzed the visible broadband photometric data of comets 1P/Halley and 4P/Faye, obtained during their perihelion passages of 1986 and 1991, respectively at the Sanglokh Observatory (Tajikistan) and the European Southern Observatory (Chile). Applying the Lomb-Scargle periodogram in the V-band magnitudes and B-V color index, we find that the most probable periodicities are 79 ± 6 and 7.36 ± 0.04 days for 1P, and 6.1 ± 0.3 days for 4P. After comparing results of color and magnitude periodograms, we argue there is a systematic difference in the number of signals identified and the level of confidence of the same periodicity in the periodograms. Our results suggest the quest for periodicities in the color of the coma of active comets should be complementary to ones in magnitudes. We have verified that the distribution of the color B-V of Faye’s coma was invariable during and after the possible occurrence of a post-perihelion outburst. We verify a symmetry in the pre- and post-perihelion H₀ photoelectric absolute magnitude of the comet Halley. The same issues were not observed in the B-V color index. We verify that the absolute magnitude H₀ of the comet Halley differs from each other when calculated from the visual or photoelectric magnitudes, due to the section of the coma used to estimate these magnitudes. We also verified that this difference in the photometric aperture can compromise comparisons of B-V color distributions between active comets.     

Space missions to small bodies: asteroids and cometary nuclei

The knowledge of the physical and dynamical properties, distribution, formation, and evolution of small bodies is fundamental to understand how planet formation occurred and, even more importantly, if and how these objects have played a role in the apparition of life on Earth. In the last century, asteroids began to no longer appear as starlike points of light in our telescopes, but to be resolved worlds with distinctly measurable sizes, shapes, and surface morphologies. Only in the last 25 years, the exploration of small bodies by spacecraft has begun and revealed objects widely diverse in formation region, evolution and properties (e.g. shape, albedo density, gravity, regolith size distribution, and porosity). In this paper we will provide a chronological analysis of comet nuclei and asteroids as revealed by space missions. The real breakthrough began with the ESA Giotto mission in 1986 to the comet Halley, while the latest JAXA Hayabusa mission was devoted to hover above the small asteroid Itokawa with a touch-and-go for a sample return of asteroidal regolith. Comet and asteroid science stands at the threshold of a new exceptional era, with many new missions to be devoted to these widely diverse and still poorly known small bodies.     

Narrow-Band Photometry of Comet C/1995 O1 (Hale-Bopp)

Photometric observations of comet C/1995 O1 (Hale-Bopp) carried out at the Stará Lesná Observatory since February to April 1997 are analyzed and discussed. Emission band fluxes and continuum fluxes are presented, from which the total numbers of molecules in the columns of the coma encircled by diaphragms are calculated. The production rates are estimated from the conventional Haser model. We found that the photometric exponent of dust contribution two months prior perihelion was n = 5.2. The photometric exponent n of the cometary magnitude solely to the C₂ emission alone equals 3.3 and that of CN equals 2.5. These values can be explained by a fact that the maximums of production rates of the gases were reached between March 2and 12 and not at the perihelion as it is valid for dust. These results are compared with the values of 1P/Halley (1986 III) under the similar conditions, obtained with the same method and instrument. C/Hale-Bopp exhibited 4.1 times more molecules radiating the CN-emission than 1P/Halley in the same column of the coma. The continuum flux of C/Hale-Bopp was also very strong. The ratios (to 1P/Halley) are 94:1 (Cont. 484.5) and 74:1 (Cont. 365.0). The cometary colour was the same as that of the Sun. This revised version was published online in July 2006 with corrections to the Cover Date.     

On the evolution and activity of cometary nuclei

The thermal evolution of a spherical cometary nucleus (initial radius of 2.5 km), composed initially of very cold amorphous ice and moving in comet Halley's orbit, is simulated numerically for 280 revolutions. It is found that the phase transition from amorphous to crystalline ice constitutes a major internal heat source. The transition does not occur continuously, but in five distinct rounds, during the following revolutions: 1, 7, 40-41, 110-112, and 248-252. Due to the (slow) heating of the amorphous ice between crystallization rounds, the phase transition front advances into the nucleus to progressively greater depths: 36 m on the first round, and then 91 m, 193 m, 381 m, and 605 m respectively. Each round of crystallization starts when when the boundary between amorphous and crystalline ice is brought to approximately 15 m below the surface, as the nucleus radius decreases due to sublimation. At the time of crystallization, the temperature of the transformed ice rises to 180 K. According to experimental studies of gas-laden amorphous ice, a large fraction of the gas trapped in the ice at low temperatures is released. Whereas some of the released gas may find its way out through cracks in the crystalline ice layer, the rest is expected to accumulate in gas pockets that may eventually explode, forming "volcanic calderas." The gas-laden amorphous ice thus exposed may be a major source of gas and dust jets into the coma, such as those observed on comet Halley by the Giotto spacecraft. The activity of new comets and, possibly, cometary outbursts and splits may also be explained in terms of explosive gas release following the transition from amorphous to crystalline ice.     

Possible collisions of P/Halley and P/Machholz 1 with Jupiter and the terrestrial planets

The perturbed motion of comet Halley and comet Mackholz 1 1986 VIII was investigated within a time interval of about 20 millennia. The minimal distance of 0.043 AU between P/Halley and Venus may occur on April 4, 4868 AD. The distance of 0.036 AU between P/Halley and Jupiter will take place on April 1, 6616 AD. The orbit of P/Machholz 1 crosses the orbits of Mercury and Venus eight times, that of the Earth six or eight times, and the orbit of Mars four times per a period of advance of the argument of perihelion. A distance of about 0.06 AU between P/Machholz 1 and the Earth may take place in August 2576 AD and 5751 AD and in February 4770 AD. The minimal comet-Earth distance of 0.035 AU occurs on September 14, 5971 AD. The closest encounter between P/Machholz 1 and Jupiter at the distance of 0.098 AU may be in May 4499 AD. These results may be considered as a forecast of possible collisions.     

Astrometric positions of comets Halley and Giacobini-Zinner obtained at Fabra Observatory

We present 190 astrometric positions of comet Halley and 46 positions of comet Giacobini-Zinner obtained at Fabra Observatory during the 1985/1986 campaign by the astrometric network of the International Halley Watch.     

Observations of the GIOTTO target comet P/Grigg-Skjellerup at the Calar Alto Observatory

CCD images of comet P/Grigg-Skjellerup, obtained for astrometric purposes with the 3.5 m telescope at the Calar Alto Observatory/Spain, were used for an analysis of the activity status of the nucleus and for a search of faint coma structures. The nucleus was found essentially inactive beyond 2.7 AU solar distance both inbound and outbound (observations on 12–13 August, 1986, 21–23 October, 1986, 22 August, 1988, 18 October, 1988, 9 and 12 September, 1991 and 3 December, 1991). The coma of the comet was well developed in May and July 1987 with a diameter of at least 190 000 km on 24 May, 1987 and of at least 80 000 km on 24 July 1987. The coma showed a cone of diffuse brightness enhancement in the sunward hemisphere. The orientation of the cone axis changed from the Sun direction in May 1987 towards about North in July 1987, i.e., it was almost perpendicular to the projected Sun-nucleus line on the sky. The cone opening angle became smaller from about 100 in May to about 50 in July 1987. A weak and narrow plasma tail was found in the images of May 1987.