Zhou Zhimin, Qin Bin, Li Shenggang, Sun Yuhan
CAS Key Laboratory of Low-Carbon Conversion Science and Engineering, Shanghai Advanced Research Institute, Chinese Academy of Sciences, 100 Haike Road, Shanghai 201210, China.
Phys Chem Chem Phys. 2021 Jan 28;23(3):1888-1895. doi: 10.1039/d0cp05947a.
In this work, we performed density functional theory (DFT)-based microkinetic simulations to elucidate the reaction mechanism of methanol synthesis on two of the most stable facets of the cubic In2O3 (c-In2O3) catalyst, namely the (111) and (110) surfaces. Our DFT calculations show that for both surfaces, it is difficult for the H atom adsorbed at the remaining surface O atom around the O vacancy (Ov) active site to migrate to an O adsorbed at the Ov due to the very high energy barrier involved. In addition, we also find that the C-O bond in the bt-CO2* chemisorption structure can directly break to form CO with a lower energy barrier than that in its hydrogenation to the COOH* intermediate in the COOH route. However, our microkinetic simulations suggest that for both surfaces, CO2 deoxygenation to form CO in both pathways, namely the COOH and CO-O routes, are kinetically slower than methanol formation under typical steady state conditions assuming a CO2 conversion of 10% and a CO selectivity of 1%. Although these results agree with previous experimental observations at relatively low reaction temperature, where methanol formation dominates, they cannot explain the predominant formation of CO at relatively high reaction temperature. We tentatively attribute this to the simplicity of our microkinetic model as well as possible structural changes of the catalyst at relatively high reaction temperature. Furthermore, although the rate-determining step (RDS) from the degree of rate control (DRC) analysis is usually consistent with that judged from the DFT calculated energy barriers, for CO2 hydrogenation to methanol over the (111) surface, our DRC analysis suggests homolytic H2 dissociation to be the rate-controlling step, which is not apparent from the DFT-calculated energy barriers. This indicates that CO2 conversion and methanol selectivity over the (111) surface can be further enhanced if homolytic H2 dissociation can be accelerated for instance by introducing transition metal dopants as already shown by some experimental observations.
在本工作中,我们基于密度泛函理论(DFT)进行了微观动力学模拟,以阐明甲醇在立方相In2O3(c-In2O3)催化剂两个最稳定晶面,即(111)和(110)表面上合成的反应机理。我们的DFT计算表明,对于这两个表面,由于涉及的能垒非常高,吸附在氧空位(Ov)活性位点周围剩余表面氧原子上的H原子很难迁移到吸附在Ov处的氧原子上。此外,我们还发现,在bt-CO2化学吸附结构中的C-O键能够直接断裂形成CO,其能垒低于在COOH路径中加氢生成COOH中间体的能垒。然而,我们的微观动力学模拟表明,对于这两个表面,在假设CO2转化率为10%和CO选择性为1%的典型稳态条件下,通过COOH和CO-O两条路径将CO2脱氧形成CO的动力学都比甲醇生成的动力学要慢。尽管这些结果与之前在相对较低反应温度下的实验观察结果一致,即在该温度下甲醇生成占主导,但它们无法解释在相对较高反应温度下CO的主要生成情况。我们初步将此归因于我们微观动力学模型的简单性以及在相对较高反应温度下催化剂可能的结构变化。此外,尽管从速率控制程度(DRC)分析得出的速率决定步骤(RDS)通常与从DFT计算能垒判断出的一致,但对于(111)表面上CO2加氢生成甲醇的反应,我们的DRC分析表明均裂H2解离是速率控制步骤,这在DFT计算的能垒中并不明显。这表明,如果例如通过引入过渡金属掺杂剂来加速均裂H2解离,如一些实验观察所示,那么(111)表面上的CO2转化率和甲醇选择性可以进一步提高。