[論文レビュー] 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

ひとこと要約

本研究では、DFTで生成された構造を入力として用い、有限サイズ補正および固定ノード誤差低減を施した、非常に高精度な量子モンテカルロ（QMC）シミュレーションを用いて、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フレームワークを確立すること。

提案手法

• 平面波基底関数を用いたPBE汎関数によるDFT計算により、H₂がCu(111)に吸着した状態、および伸長・解離した状態の初期構造を生成した。
• QMCシミュレーションに使用する試行波動関数は、DFT PBE軌道から得られ、線形ジャストロウ最適化により精度を向上させた。
• 統計誤差を大規模なサンプリングにより最小限に抑えるために、拡散モンテカルロ（DMC）法を用いて高精度な基底状態エネルギーを計算した。
• 有限サイズ効果は、3×3×1 kポイントグリッドを用いたねじれ平均化により補正され、遷移状態と漸近的状態との間で誤差が約90％低減された。
• 固定ノード近似に起因する系統的誤差は±1.5 kcal/molと推定され、統計誤差バーに加算された。
• 周期的系に適したCASINOコードを用い、系のスケーリングに効率的な立方スプライン（ブリップ）展開を採用した。

実験結果

リサーチクエスチョン

• RQ1 高精度な電子相関を考慮したQMC法を用いた場合、Cu(111)表面上のH₂解離の活性化エネルギー障壁はどの程度か？
• RQ2 有限サイズ効果および固定ノード誤差は、表面反応のQMCシミュレーションの精度にどのように影響するか？
• RQ3 QMCは、DFTよりも実験的分子ビームデータを高い精度で再現できるか？
• 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未満の精度で予測可能であり、新たなベンチマークを確立した。