Atomistic Modeling of Methyl Formate and Glycolaldehyde Formation on Interstellar Dirty Ice Mantles via a “Radical + Ice” Mechanism

Methylformate (MF) and glycolaldehyde (GA) are two primogenital organic molecules detected in both cold and warm regions of the interstellar medium (ISM). Both gas-phase and grain-surface pathways have been proposed to explain their abundances, yet uncertainties remain, since prevailing grain-surface mechanisms favor the formation of GA over MF, which mismatch observations in different ISM regions. In this work, MF and GA synthetic reactions are atomistically modeled on surfaces containing variable H <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mrow></mml:mrow> <mml:mn>2</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$_2$</mml:annotation></mml:semantics> </mml:math> O and CO percentages (interstellar dirty ices), in which one of the reactants coming from the gas phase reacts with an icy CO, thus adopting the following two-step "radical + ice" mechanism: for MF, OCH <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mrow></mml:mrow> <mml:mn>3</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$_3$</mml:annotation></mml:semantics> </mml:math> + <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>CO</mml:mi> <mml:mrow><mml:mo>(</mml:mo> <mml:mi>ice</mml:mi> <mml:mo>)</mml:mo></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$${\mathrm{CO}}_{(\mathrm{ice})}$$</mml:annotation></mml:semantics> </mml:math> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>→</mml:mo></mml:mrow> <mml:annotation>$\to$</mml:annotation></mml:semantics> </mml:math> COOCH <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mrow></mml:mrow> <mml:mn>3</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$_3$</mml:annotation></mml:semantics> </mml:math> + H <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>→</mml:mo></mml:mrow> <mml:annotation>$\to$</mml:annotation></mml:semantics> </mml:math> HCOOCH <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mrow></mml:mrow> <mml:mn>3</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$_3$</mml:annotation></mml:semantics> </mml:math> ; for GA, CH <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mrow></mml:mrow> <mml:mn>2</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$_2$</mml:annotation></mml:semantics> </mml:math> OH + <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mi>CO</mml:mi> <mml:mrow><mml:mo>(</mml:mo> <mml:mi>ice</mml:mi> <mml:mo>)</mml:mo></mml:mrow> </mml:msub> </mml:mrow> <mml:annotation>$${\mathrm{CO}}_{(\mathrm{ice})}$$</mml:annotation></mml:semantics> </mml:math> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>→</mml:mo></mml:mrow> <mml:annotation>$\to$</mml:annotation></mml:semantics> </mml:math> COCH <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mrow></mml:mrow> <mml:mn>2</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$_2$</mml:annotation></mml:semantics> </mml:math> OH + H <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics><mml:mrow><mml:mo>→</mml:mo></mml:mrow> <mml:annotation>$\to$</mml:annotation></mml:semantics> </mml:math> HCOCH <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msub><mml:mrow></mml:mrow> <mml:mn>2</mml:mn></mml:msub> </mml:mrow> <mml:annotation>$_2$</mml:annotation></mml:semantics> </mml:math> OH. Calculations show that the first step presents an energy barrier (32-38 kJ mol <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msup><mml:mrow></mml:mrow> <mml:mrow><mml:mo>-</mml:mo> <mml:mn>1</mml:mn></mml:mrow> </mml:msup> </mml:mrow> <mml:annotation>$^{-1}$</mml:annotation></mml:semantics> </mml:math> for MF and 17-20 kJ mol <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msup><mml:mrow></mml:mrow> <mml:mrow><mml:mo>-</mml:mo> <mml:mn>1</mml:mn></mml:mrow> </mml:msup> </mml:mrow> <mml:annotation>$^{-1}$</mml:annotation></mml:semantics> </mml:math> for GA), while the second step is nearly barrierless. Although the energetics favor GA formation, the observed abundances are better explained by desorption phenomena rather than reaction barriers are argued. Specifically, the weaker binding energies of MF (16.8-46.1 kJ mol <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msup><mml:mrow></mml:mrow> <mml:mrow><mml:mo>-</mml:mo> <mml:mn>1</mml:mn></mml:mrow> </mml:msup> </mml:mrow> <mml:annotation>$^{-1}$</mml:annotation></mml:semantics> </mml:math> ) than GA (28.4-90.2 kJ mol <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:semantics> <mml:mrow><mml:msup><mml:mrow></mml:mrow> <mml:mrow><mml:mo>-</mml:mo> <mml:mn>1</mml:mn></mml:mrow> </mml:msup> </mml:mrow> <mml:annotation>$^{-1}$</mml:annotation></mml:semantics> </mml:math> ) support its higher abundance in the ISM.