Simulations of ice chemistry in cometary nuclei

Robin T. Garrod
DOI: https://doi.org/10.3847/1538-4357/ab418e
2019-09-26
Abstract:The first computational model of solid-phase chemistry in cometary nuclear ices is presented. An astrochemical kinetics model, MAGICKAL, is adapted to trace the chemical evolution in multiple layers of cometary ice, over a representative period of 5 Gyr. Physical conditions are chosen appropriate for "cold storage" of the cometary nucleus in the outer Solar System, prior to any active phase. The chemistry is simulated at a selection of static temperatures in the range 5 - 60 K, while the ice is exposed to the interstellar radiation field, inducing a photochemistry in the outer ice layers that produces significant formation of complex organic molecules. A treatment for the chemistry resulting from cosmic-ray bombardment of the ices is also introduced into the model, along with a new formulation for low-temperature photochemistry. Production of simple and complex molecules to depth on the order of 10~m or more is achieved, with local fractional abundances comparable to observed values in many cases. The production of substantial amounts of O$_2$ (and H$_2$O$_2$) is found, suggesting that long-term processing by high-energy cosmic rays of cometary ices in situ, over a period on the order of 1 Gyr, may be sufficient to explain the large observed abundances of O$_2$, if the overall loss of material from the comet is limited to a depth on the order of 10 m. Entry into the inner solar system could produce a further enhancement in the molecular content of the nuclear ices that may be quantifiable using this modeling approach.
Earth and Planetary Astrophysics,Astrophysics of Galaxies,High Energy Astrophysical Phenomena
What problem does this paper attempt to address?