[논문 리뷰] Ab initio instanton rate theory made efficient using Gaussian process regression

Gabriel Laude, Danilo Calderini외 2명|arXiv (Cornell University)|2018. 05. 07.

Advanced Chemical Physics Studies인용 수 4

한 줄 요약

이 논문은 정확한 양자 반응 속도 예측을 위해 필요한 고수준 전자구조 계산의 수를 극적으로 감소시키는 가우시안 프로세스 회귀(GPR) 기반의 아비니티 인스턴턴 이론을 제안한다. 인스턴턴 경로 주변의 국소적 포텐셜 에너지 표면을 GPR로 피팅함으로써, 표준 방법에 비해 최대 10배 적은 아비니티 평가를 통해 기준 결과에 1% 이내로 수렴하는 것을 달성하였으며, 이는 H + C2H6와 같은 반응에서 UCCSD(T)-F12b 수준에서 효율적이고 고정밀도의 반응 속도 예측을 가능하게 한다.

ABSTRACT

Ab initio instanton rate theory is a computational method for rigorously including tunnelling effects into calculations of chemical reaction rates based on a potential-energy surface computed on the fly from electronic-structure theory. This approach is necessary to extend conventional transition-state theory into the deep-tunnelling regime, but is also more computationally expensive as it requires many more ab initio calculations. We propose an approach which uses Gaussian process regression to fit the potential-energy surface locally around the dominant tunnelling pathway. The method can be converged to give the same result as from an on-the-fly ab initio instanton calculation but requires far fewer electronic-structure calculations. This makes it a practical approach for obtaining accurate rate constants based on high-level electronic-structure methods. We show fast convergence to reproduce benchmark H + CH4 results and evaluate new low-temperature rates of H + C2H6 in full dimensionality at a UCCSD(T)-F12b/cc-pVTZ-F12 level.

연구 동기 및 목표

• 아비니티 인스턴턴 이론의 계산 비용을 줄이기 위해, 전체 인스턴턴 경로를 따라 에너지, 기울기, 헤시안을 필요로 하는 고수준 전자구조 계산의 빈도를 감소시키는 것.
• 인스턴턴 경로를 따라 실시간으로 계산하는 아비니티 인스턴턴 계산의 비효율성을 해결하기 위해, 전이 상태 이론보다 더 비싼 계산 비용을 유발하는 많은 전자구조 평가가 필요로 하는 문제를 해결하는 것.
• 고정밀도 전자구조 방법을 사용하면서도 고가의 계산 횟수를 크게 줄일 수 있는 방법을 개발하여 높은 정확도를 유지하는 것.
• 복잡한 다원자 반응에 대해 아비니티 인스턴턴 이론을 실용적으로 적용할 수 있도록, 고수준 이론 수준에서 계산 가능성을 확보하는 것.

제안 방법

• 가우시안 프로세스 회귀(GPR)를 사용하여 인스턴턴 경로를 따라 샘플링된 아비니티 데이터 포인트에 기반해 서브스티튜션 포텐셜 에너지 표면(PES)을 구축한다.
• GPR 모델은 고수준 전자구조 계산(에너지, 기울기, 헤시안)에서 얻은 소량의 적응형 데이터 포인트를 기반으로 훈련되며, 이는 레인지-폴리머 인스턴턴 경로의 핵심 지점에서 계산된다.
• 예측 불확실성에 기반해 새로운 훈련 포인트를 적응적으로 선택함으로써 진정한 인스턴턴 경로와 속도 상수에 수렴하도록 보장한다.
• GPR 서브스티튜션은 표준 실시간 인스턴턴 방법보다 훨씬 적은 고수준 전자구조 평가를 통해 정확한 인스턴턴 속도 계산을 가능하게 한다.
• 수렴된 속도를 기준 전체 아비니티 인스턴턴 계산 결과와 비교하여 검증하였으며, 1% 이내의 일치를 보였다.
• 이 방법은 H + CH4 및 H + C2H6에 적용되었으며, UCCSD(T)-F12b/cc-pVTZ-F12 수준에서 약 6개의 헤시안만으로도 전차원적 인스턴턴 속도를 도출하였다.

실험 결과

연구 질문

• RQ1 가우시안 프로세스 회귀(GPR)를 사용하여 고수준 전자구조 계산의 수를 최소화하면서도 인스턴턴 경로를 따라 포텐셜 에너지 표면을 정확하게 근사화할 수 있는가?
• RQ2 표준 실시간 결과에 비해 1% 이내로 인스턴턴 속도가 수렴하기 위해 필요한 최소한의 아비니티 평가 횟수는 얼마인가?
• RQ3 기준 반응인 H + CH4에 대해 GPR 기반 가속 인스턴턴 방법은 전통적인 아비니티 인스턴턴 이론에 비해 정확도와 효율성 면에서 어떻게 비교되는가?
• RQ4 이 방법은 H + C2H6와 같은 복잡한 다원자 반응에 대해 고수준 이론 수준에서 실용적인 계산 비용으로 확장 가능한가?
• RQ5 GPR 프레임워크는 저수준 전자구조 데이터와 조합하여 고수준 계산 횟수를 더 줄일 수 있으며, 정확도를 손상시키지 않을 수 있는가?

주요 결과

• GPR 기반 가속 인스턴턴 방법은 표준 방법에 비해 필요한 전자구조 계산 수를 극히 적게 사용하면서도 H + CH4 반응에 대해 기준 전체 아비니티 인스턴턴 속도에 1% 이내로 수렴한다.
• H + C2H6 반응에 대해 이 방법은 UCCSD(T)-F12b/cc-pVTZ-F12 수준에서 약 6개의 헤시안만으로도 전차원적 인스턴턴 속도를 도출하였으며, 이는 고전적 TST 수준과 거의 동일한 효율성을 보였다.
• 인스턴턴 경로가 최소 에너지 경로에서 크게 이탈하더라도 높은 정확도를 유지함으로써, 복잡한 반응 메커니즘에서의 강건성을 입증하였다.
• GPR의 사용은 고리-폴리머 비드 수와 훈련 세트 크기 측면에서 체계적인 수렴을 가능하게 하여 신뢰할 수 있는 결과를 보장한다.
• 이 방법은 이전에는 많은 시스템에서 계산 비용이 지나치게 높아 실현 가능하지 않았던 고수준 전자구조 방법을 인스턴턴 이론에 적용할 수 있도록 하였다.
• 이 방법은 핵심 단계에 대해서만 효율적이고 고정밀도의 인스턴턴 계산을 가능하게 하여, 복잡한 반응 네트워크에서 터널링 효과를 체계적으로 포함할 수 있는 길을 열어주었다.