[논문 리뷰] Quantum Monte Carlo activation barrier for hydrogen dissociation on copper to unprecedented accuracy

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

Advanced Chemical Physics Studies참고 문헌 11인용 수 4

한 줄 요약

이 연구는 유체역학적 시뮬레이션을 사용하여 Cu(111) 상에서 H₂ 해리의 활성화 장벽을 매우 정확하게 계산한다. DFT로 생성한 기하구조를 입력으로 사용하고, 유한한 크기 보정 및 고정노드 오차 보정을 적용한다. 최종 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 프레임워크를 수립하기.

제안 방법

• H₂가 Cu(111)에 흡착된 상태의 초기 기하구조를 평면파 기저에서 PBE 기능을 사용한 DFT 계산을 통해 생성하며, 연장 및 해리된 구조를 포함한다.
• QMC 시뮬레이션에 사용하기 위해 DFT PBE 오비탈에서 유도된 시험 파동함수를 사용하고, 정확도 향상을 위해 선형 Jastrow 최적화를 수행한다.
• 통계적 오차를 대규모 샘플링을 통해 최소화하기 위해 확산 몽테카를로(DMC)를 사용하여 고정밀도 기초 상태 에너지를 계산한다.
• 유한한 크기 효과는 3×3×1 k-점 메esh에서의 꼬임 평균화를 통해 보정하며, 전이 상태와 점 渐진 상태 사이에서 오차를 약 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 이내의 정밀도를 달성할 수 있음을 입증하며, 새로운 기준을 설정한다.