Combined quantum chemical and modeling study of CO hydrogenation on water ice

¹ Departament de Química, Universitat Autònoma de Barcelona, 08193 Bellaterra ( Cerdanyola del Vallès), Spain
e-mail: Albert.Rimola@uab.cat
² Astrochemistry Laboratory and The Goddard Center for Astrobiology, Mailstop 691, NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD 20770, USA
³ Dipartimento di Chimica, and NIS – Centre for Nanostructured Interfaces and Surfaces, Università di Torino, via P. Giuria 7, 10125 Torino, Italy
⁴ Dipartimento di Chimica, Biologia e Biotecnologie, Università di Perugia, via Elce di Sotto 8, 06123 Perugia, Italy
⁵ Univ. Grenoble Alpes, IPAG, 38000 Grenoble, France
⁶ CNRS, IPAG, 38000 Grenoble, France

Received: 22 April 2014
Accepted: 6 October 2014

Abstract

Context. Successive hydrogenation reactions of CO on interstellar icy grain surfaces are considered one of the most efficient mechanisms in interstellar environments for the formation of H₂CO and CH₃OH, two of the simplest organic molecules detected in space. In the past years, several experimental and theoretical works have been focused on these reactions, providing relevant information both at the macroscopic and atomic scale. However, several questions still remain open, such as the exact role played by water in these processes, a crucial aspect because water is the dominant constituent of the ice mantles around dust grain cores.

Aims. We here present a quantum chemical description of the successive H additions to CO both in the gas phase and on the surfaces of several water clusters.

Methods. The hydrogenation steps were calculated by means of accurate quantum chemical methods and structural cluster models consisting of 3, 18, and 32 water molecules.

Results. Our main result is that the interaction of CO and H₂CO with the water cluster surfaces through H-bonds with the O atoms increases the C−O polarization, thus weakening the C−O bond. Consequently, the C atoms are more prone to receiving H atoms, which in turn lowers the energy barriers for the H additions compared to the gas-phase processes. The calculated energy barriers and transition frequencies associated with the reaction coordinate were adopted as input parameters in our numerical model of the surface chemistry (GRAINOBLE) to simulate the distribution of the H₂CO and CH₃OH ice abundances (with respect to water). Our GRAINOBLE results based on the energy barriers and transition frequencies for the reactions on the 32 water molecule cluster compare well with the observed abundances in low-mass protostars and dark cores.