Context:Ab initio molecular dynamics simulations were carried out to study the formation pathways to complex organic molecules when a OH + projectile hit an interstellar dust grain covered only by methanol molecules. The selected target material is a methanol cluster formed by ten units (CH 3 OH) 10 . Aims: The focus of this paper is the process where methanediol CH 2 (OH) 2 and hydroxymethyl CH 2 OH + , both key organic intermediate molecules, were involved in the formation mechanisms of stable complex organic molecules (COMs). Methods: We performed Born-Oppenheimer (ab initio) molecular dynamics (BOMD) simulations under the hybrid functional of Head-Gordon ωB97X-D. We used the initial kinetic impact energy of 10, 12, 15, 18, 20, and 22 eV. Results: We corroborate that CH 2 (OH) 2 and CH 2 OH + are the main precursors to form molecules such as methoxymethanol CH 3 OCH 2 OH, the formyl radical HCO, the Criegee biradical CH 2 OO, and formaldehyde H 2 CO and its elusive HCOH isomer. We discuss the mechanism formation of these complex organic molecules. We compare the formation pathways with previous theoretical results where both key intermediates are present. The pathways in some cases go through CH 2 (OH) 2 or undergo by CH 2 OH + . Conclusions: We confirm that CH 2 (OH) 2 and CH 2 OH + play a key role on the path to the formation of abundant H 2 CO. These mechanisms can give insight into alternative pathways relevant to understanding experimental processes with key steps within those precursors.