Radial distribution of production rates, loss rates and densities corresponding to ion masses ≤40 amu in the inner coma of Comet Halley: composition and chemistry

Abstract

In this paper we have studied the chemistry of C, H, N, O, and S compounds corresponding to ions of masses ≤40 amu in the inner coma of the Comet 1P/Halley. The production rates, loss rates, and ion mass densities are calculated using the Analytical Yield Spectrum approach and solving coupled continuity equation controlled by the steady state photochemical equilibrium condition. The primary ionization sources in the model are solar EUV photons, photoelectrons, and auroral electrons of the solar wind origin. The chemical model couples ion-neutral, electron-neutral, photon-neutral and electron-ion reactions among ions, neutrals, electrons, and photons through over 600 chemical reactions. Of the 46 ions considered in the model the chemistry of 24 important ions (viz., CH₃OH⁺₂, H₃CO⁺, NH⁺₄, H₃S⁺, H₂CN⁺, H₂O⁺, NH⁺₃, CO⁺, C₃H⁺₃, OH⁺, H₃O⁺, CH₃OH⁺, C₃H⁺₄, C₂H⁺₂, C₂H⁺, HCO⁺, S⁺, CH⁺₃, H₂S⁺, O⁺, C⁺, CH⁺₄, C⁺₂, and O⁺₂) are discussed in this paper. At radial distances <1000 km, the electron density is mainly controlled by 6 ions, viz., NH⁺₄, H₃O⁺, CH₃OH⁺₂, H₃S⁺, H₂CN⁺, and H₂O⁺, in the decreasing order of their relative contribution. However, at distances >1000 km, the 6 major ions are H₃O⁺, CH₃OH⁺₂, H₂O⁺, H₃CO⁺, C₂H⁺₂, and NH⁺₄; along with ions CO⁺, OH⁺, and HCO⁺, whose importance increases with further increase in the radial distance. It is found that at radial distances greater than ~1000 km (±500 km) the major chemical processes that govern the production and loss of several of the important ions in the inner coma are different from those that dominate at distances below this value. The importance of photoelectron impact ionization, and the relative contributions of solar EUV, and auroral and photoelectron ionization sources in the inner coma are clearly revealed by the present study. The calculated ion mass densities are compared with the Giotto Ion Mass Spectrometer (IMS) and Neutral Mass Spectrometer (NMS) data at radial distances 1500, 3500, and 6000 km. There is a reasonable agreement between the model calculation and the Giotto measurements. The nine major peaks in the IMS spectra between masses 10 and 40 amu are reproduced fairly well by the model within a factor of two inside the ionopause. We have presented simple formulae for calculating densities of the nine major ions, which contribute to the nine major peaks in the IMS spectra, throughout the inner coma that will be useful in estimating their densities without running the complex chemical models.