[论文解读] Quantum Monte Carlo activation barrier for hydrogen dissociation on copper to unprecedented accuracy

Philip E. Hoggan|arXiv (Cornell University)|Nov 24, 2015

Advanced Chemical Physics Studies参考文献 11被引用 4

一句话总结

本研究采用高精度的量子蒙特卡罗（QMC）模拟，基于DFT生成的几何构型，计算Cu(111)表面H₂解离的活化能垒，应用有限尺寸校正与固定节点误差缓解方法。最终的QMC校正能垒为14.79 kcal/mol，统计误差为0.55 kcal/mol，与实验分子束测得的14.48 kcal/mol值高度一致，展示了在预测表面催化能垒方面前所未有的精确度。

ABSTRACT

Many chemical reactions involve bond-dissociation. This is also true for reactions at solid surfaces, in which the dissociation step is often limiting but facilitated in comparison to gas phase reaction channels. This work considers hydrogen dissociation. Reliable molecular beam results are available for this reaction at some copper surfaces. Heterogeneous catalysis by copper is simulated. It was investigated in our previous work since it is in many ways a prototype metal presenting a close-packed surface here. These hydrogen molecules are adsorbed at Cu(111) and fixed geometries on the dissociation reaction pathway for stretched and distant equilibrium H$_2$ are given by using Density Functional Theory (DFT) calculations in a plane wave basis. The PBE wave-functions at these bond-lengths serve as trial input for Quantum Monte Carlo (QMC) simulations of the ground states to obtain highly accurate correlated results for the associated activation barriers indicating the catalytic effect on this dissociation. This correlation varies as bonds dissociate, requiring its accurate evaluation. Finite size effects and fixed-node error are possible limitations to accuracy of this type of QMC study. We are able to limit fixed node error, using certain trial wave-functions. The finite size effect is considerable, although comparing two adsorbed geometries cancels about 90% with respect to clean surfaces. The pseudo-potential used to represent the atomic core of copper must also be determined carefully: we leave 11 active electrons but include the 3d shell in the pseudo-potential.

研究动机与目标

• 使用量子蒙特卡罗（QMC）方法，以前所未有的精度计算Cu(111)表面H₂解离的活化能垒。
• 克服DFT在描述化学键解离过程中电子相关性的局限性，尤其是在解离极限区域。
• 通过扭曲平均法与优化的试探波函数，减小QMC中的系统性误差，包括有限尺寸效应与固定节点误差。
• 将QMC结果与实验分子束数据进行验证，为表面催化模拟提供基准参考。
• 建立可靠的QMC框架，用于研究涉及强电子关联的多相催化中决速步的反应能垒。

提出的方法

• 采用PBE泛函在平面波基组下进行DFT计算，生成H₂吸附在Cu(111)表面的初始几何构型，包括拉伸与解离构型。
• 在QMC模拟中使用基于DFT PBE轨道构建的试探波函数，并通过线性Jastrow优化提升精度。
• 采用扩散蒙特卡罗（DMC）方法计算高精度基态能量，通过大规模采样最小化统计误差。
• 利用3×3×1 k点网格上的扭曲平均法校正有限尺寸效应，使过渡态与渐近态之间的误差降低约90%。
• 固定节点近似带来的系统性误差估计为±1.5 kcal/mol，并加入统计误差条中。
• 使用CASINO代码处理周期性体系，采用三次样条（blip）展开以实现与体系尺寸的高效缩放。

实验结果

研究问题

• RQ1 当使用高度相关联的QMC方法计算时，Cu(111)表面H₂解离的活化能垒是多少？
• RQ2 有限尺寸效应与固定节点误差如何影响QMC模拟在表面反应中的精度？
• RQ3 QMC能否比DFT更准确地再现Cu(111)表面H₂解离的实验分子束数据？
• RQ4 试探波函数与扭曲平均法在多大程度上可降低QMC在催化能垒模拟中的系统性误差？
• RQ5 在H₂解离过程中，电子相关性如何变化？QMC能否准确捕捉这一变化？

主要发现

• 原始QMC活化能垒为15.00 ± 0.55 kcal/mol，已施加有限尺寸与固定节点校正。
• 校正后最终的QMC能垒（c-QMC）为14.79 ± 0.55 kcal/mol，与实验值14.48 kcal/mol高度一致。
• QMC与实验结果的偏差仅0.31 kcal/mol，表明具有极高的精度与极低的系统性误差。
• 统计误差已降低至0.55 kcal/mol，总不确定性包含1.5 kcal/mol的系统性误差估计。
• 有限尺寸效应在过渡态与渐近几何构型之间降低了90%，3×3×1 k点网格有效实现了校正。
• 本研究证明，QMC可实现表面催化能垒的亚0.5 kcal/mol精度，树立了新的基准。