Hydrogenation of CO to methanol over In-doped m-ZrO: a DFT investigation into the oxygen vacancy size-dependent reaction mechanism.

Selective methanol synthesis CO hydrogenation has been thoroughly investigated over defective In-doped m-ZrO using density functional theory (DFT). Three types of oxygen vacancies (Ovs) generated either at the top layer (O1_v and O4_v) or at the subsurface layer (O2_v) are chosen as surface models due to low Ov formation energy. Surface morphology reveals that O1_v has smaller oxygen vacancy size than O4_v. Compared with perfect In@m-ZrO, indium on both O1_v and O4_v is partially reduced, whereas the Bader charge of In on O2_v remains almost the same. Our calculations show that CO is moderate in adsorption energy (∼-0.8 eV) for all investigated surface models, which facilitates the formate pathway for both O1_v and O4_v. O2_v is not directly involved in CO methanolization but could readily transform into O1_v once CO/H feed gas is introduced. Based on the results, the synthesis of methanol from CO hydrogenation turns out to exhibit conspicuous vacancy size-dependency for both O1_v and O4_v. The reaction mechanism for small-sized O1_v is controlled by both the vacancy size effect and surface reducibility effect. Thus, HCOO* favors direct C-O bond cleavage (c-mechanism) before further hydrogenation to methanol, which is similar to the defective InO. The vacancy size effect is more competitive than the surface reducibility effect for large-sized O4_v. Therefore, HCOO* prefers protonation to HCOOH before C-O bond cleavage (p-mechanism) which is similar to the ZnO-ZrO solid solution. Furthermore, we also determined that stable-CHO*, which is too stable to be hydrogenated, originates from the O1_v surface. In contrast, CHO* with similar configuration is allowed to be further converted to methanol on O4_v. Overall, our findings offer a new perspective towards how reaction mechanisms are determined by the size of oxygen vacancies.