A systematic theoretical study of the consecutive hydrogenation reactions of the CO molecule for the methanol synthesis catalyzed by different oxides of Zn, Ce and Ga is reported in this work. First, the CO hydrogenation with the formation of formyl species (HCO) was analyzed, followed by the successive hydrogenations that lead to formaldehyde (H2CO), methoxy (H3CO) and, finally, methanol (H3COH). The co-adsorption with H, in almost all the intermediate species, allows the corresponding hydrogenation reaction. Oxygen vacancies promote the reactivity in the generation of both formaldehyde and methoxy species. The formation of these species involves an important geometric difference between the initial and the final states, leading to high activation barriers. Comparing the surfaces studied in this work, we found that ZnO (0001)vacO has shown to be of a greater interest for methanol synthesis. However, the foregoing is not the most relevant of our results, but, instead, that the Brönsted Evans Polanyi (BEP) relationships between the initial or the final states and the transition states (TS) allowed to find a very good correlation between surface structure and reactivity.