[논문 리뷰] Understanding Discrepancies of Wavefunction Theories for Large Molecules

Tobias Schäfer, Andreas Irmler외 2명|arXiv (Cornell University)|2024. 07. 01.

Quantum optics and atomic interactions인용 수 10

한 줄 요약

이 논문은 CCSD(T) 이론이 다체 섭동 급수의 절단으로 인해 큰 비공유 분자 복합체를 과도하게 결합시킨다는 점을 규명하고, 이를 보정한 CCSD(cT) 방법을 제안한다. 이는 132원자까지의 시스템에서 확산 양자 몽테카를로(DMC)와의 일치를 이끌어내는 화학적 정확도(1 kcal/mol 이내)를 달성한다. 평면파 기저와 체계적인 상관관계 처리를 통해 저자들은 CCSD(T)와 DMC 사이의 오랫동안 지속된 모순을 해결하였으며, 이는 대규모 재료 및 약물 설계 응용 분야에서 신뢰할 수 있는 벤치마크 예측을 가능하게 한다.

ABSTRACT

Quantum mechanical many-electron calculations can predict properties of atoms, molecules and even complex materials. The employed computational methods play a quintessential role in many scientifically and technologically relevant research fields. However, a question of paramount importance is whether approximations aimed at reducing the computational complexity for solving the many-electron Schrödinger equation, are accurate enough. Here, we investigate recently reported discrepancies of noncovalent interaction energies for large molecules predicted by two of the most widely-trusted many-electron theories: diffusion quantum Monte Carlo and coupled-cluster theory. We are able to unequivocally pin down the source of the puzzling discrepancies and present modifications to widely-used coupled-cluster methods needed for more accurate noncovalent interaction energies of large molecules on the hundred-atom scale. This enhances the reliability of predictions from quantum mechanical many-electron theories across a wide range of critical applications, including drug design, catalysis, and the innovation of new functional materials, such as those for renewable energy technologies.

연구 동기 및 목표

• 큰 분자 복합체에서 CCSD(T)와 DMC 사이의 비공유 상호작용 에너지에 대한 오랫동안 지속된 모순을 해결하기 위해.
• 특히 수백 개의 원자를 포함하는 큰 시스템에서 이러한 모순의 근본 원인을 규명하기 위해.
• CCSD(T)의 과도한 결합 문제를 수정하는 수정된 결합 클러스터 방법인 CCSD(cT)를 개발하고 검증하기 위해.
• 벤치마크 수준의 결과를 얻기 위해 평면파 기저 집합과 완전한 상관관계 처리를 사용하는 신뢰할 수 있고 체계적으로 개선 가능한 계산 프레임워크를 구축하기 위해.
• 약물 설계, 촉매 작용 및 재생 가능 에너지 재료 분야의 대규모 응용 분야에서 정확하고 고신뢰도의 예측을 가능하게 하기 위해.

제안 방법

• 원자 중심의 가우시안 함수에서 기인하는 기저 집합의 불완전성과 편향을 제거하기 위해 평면파 기저 집합을 사용하였다.
• 지난 연구에서 사용된 局부 상관관계 근사치를 피하기 위해, 완전하고 비국소적인 상관관계 처리를 수행한 캐논리컬 CCSD(T) 구현을 사용하였다.
• 단순하고 체계적인 보정을 통해 섭동 3중 자극을 조정함으로써 과도한 결합을 수정하는 수정된 결합 클러스터 방법인 CCSD(cT)를 도입하였다.
• LNO-CCSD(T) 데이터에서 유도된 선형 회귀 모델(CCSD(cT)-fit)을 사용하여 상호작용 에너지를 추정하였으며, 전체 CCSD(T) 및 DMC 결과와의 검증을 통해 타당성을 입증하였다.
• CCSD(T) 및 CCSD(cT)와의 비교를 위해 오차 구간을 포함한 대규모 DMC 계산을 수행하여 여러 구현 간의 일관성을 확보하였다.
• 특히 C60@[6]CPPA와 같이 높은 복잡성을 지닌 시스템에 대해 신뢰도를 평가하기 위해 체계적 오차 추정 및 상호 검증을 수행하였다.

실험 결과

연구 질문

• RQ1 CCSD(T)와 DMC가 모두 매우 정확한 것으로 간주되지만, 큰 분자 복합체에서 비공유 상호작용 에너지가 상당히 다를 수 있는 이유는 무엇인가?
• RQ2 이 모순은 CCSD(T)의 근사치, 예를 들어 섭동 급수의 절단 때문인가, 아니면 DMC의 유한한 크기 효과나 기저 집합 제한 때문인가?
• RQ3 CCSD(T)에서 관찰된 큰 비공유 시스템에서의 과도한 결합 문제를 해결하기 위해 수정된 CCSD 방법인 CCSD(cT)가 가능할 수 있는가?
• RQ4 CCSD(cT) 방법이 다양한 큰 분자 시스템에서 DMC 결과와 비교해 화학적 정확도(1 kcal/mol 이내)를 얼마나 잘 달성할 수 있는가?
• RQ5 C60@[6]CPPA와 같이 100개 이상의 원자를 포함하는 복잡한 시스템에 적용했을 때 CCSD(cT) 방법의 강건성은 어느 정도인가?

주요 결과

• 큰 분자에서 비공유 상호작용 에너지에 대한 CCSD(T)와 DMC 사이의 모순은 주로 CCSD(T)에서 다체 섭동 급수의 절단으로 인해 발생하며, 이는 과도한 결합을 초래한다.
• CCSD(T)의 보정된 형태인 CCSD(cT)는 과도한 결합 오차를 감소시키며, 테스트된 모든 시스템, 특히 132원자인 C60@[6]CPPA 복합체에서 DMC와의 일치를 1 kcal/mol 이내로 달성한다.
• 평행 이격 벤젠 이량체의 경우, 평면파 기반 CCSD(T) 결과가 가우시안 기저 집합 기반 기준값과 잘 일치하여 계산 프레임워크의 타당성을 검증하였다.
• LNO-CCSD(T) 데이터에서 유도된 CCSD(cT)-fit 모델은 뛰어난 예측 능력을 보이며, GGG, GCGC 및 C2C2PD 시스템에서 오차가 0.5 kcal/mol 이내이다.
• 여러 사례에서 DMC 결과의 오차 구간을 초월하는 DMC 추정치의 모순이 관찰되었으며, 이는 큰 시스템에서 DMC 결과의 불확실성을 시사한다.
• 남아 있는 작은 모순(예: C60@[6]CPPA의 경우 1.1 kcal/mol)은 DMC, CCSD(cT)-fit 또는 LNO-CCSD(T) 오차 추정의 불확실성 때문일 것이며, CCSD(cT)의 기본적인 결함은 아니다.