[논문 리뷰] The S66 Noncovalent Interaction Benchmark Re-examined: Composite Localized Coupled Cluster Approaches

Emmanouil Semidalas, Golokesh Santra외 2명|arXiv (Cornell University)|2021. 11. 02.

Advanced Chemical Physics Studies인용 수 1

한 줄 요약

이 연구는 S66 기준값에서 비공유 상호작용(NCI) 에너지 계산을 위한 비용 효율적인 복합 국소화된 결합형 코풀러-클러스터(cLNO-CCSD(T))) 방법을 제안한다. 더 작은 기저함수에서의 외삽된 LNO-CCSD(T) 에너지와 더 넓은 임계값 계산에서의 보정항을 조합함으로써, 고정밀 기준값과의 RMS 편차가 0.05 kcal/mol 이하인 정확도를 달성한다. 이는 전체 vTight/vvTight 계산의 비용의 일부에 불과한 비용으로도 달성 가능한 수준이다.

ABSTRACT

The S66 non-covalent interactions are studied through localized coupled-cluster methods and general LNO-CCSD(T)-based composite schemes. Very small RMS deviations (\leq 0.05 kcal/mol) for the low-cost composite approaches from the SILVER reference interaction energies of S66 indicate that we can safely avoid carrying out the largest basis set calculations with veryVeryTight thresholds, and apply instead additivity corrections in smaller basis sets. Interestingly, the counterpoise corrections do not have an appreciable effect on the composite schemes. These findings may prove useful for intermolecular and intramolecular NCIs of larger systems.

연구 동기 및 목표

• 비공유 상호작용(NCI) 에너지 계산을 위한 저비용, 고정밀 복합 방법 개발
• 복합 LNO-CCSD(T) 방법이 고비용 vTight/vvTight 임계값 계산의 정확도에 도달하거나 초월할 수 있는지 평가
• 이러한 복합 방법에서 카운터포인트(CP) 보정의 필요성과 영향 평가
• 표준 고정밀 LNO-CCSD(T) 계산이 비용이 너무 많이 들기 때문에 수행이 어려운 대규모 시스템에 대한 계산 효율적인 대안 제공

제안 방법

• 복합 방법은 더 작은 기저함수(haVTZ, haVQZ)에서의 LNO-CCSD(T) 에너지와 더 좁은 임계값 계산에서의 보정항을 조합한다.
• 기저함수 카디널 수 T, Q, 5를 사용하여 관련 에너지와 SCF 에너지에 대해 이중점 외삽을 적용하며, α=3.22 및 α=3.0을 각각 사용한다.
• 복합 방법의 최적 선형 조합 계수는 SILVER 기준 상호작용 에너지와의 RMSD를 최소화함으로써 결정된다.
• 카운터포인트(CP) 보정을 적용하고, CP 및 비보정 결과의 평균인 '반-CP' 결과와 비교하여 정확도에 미치는 영향을 평가한다.
• 벽 시계(wall time) 측정을 통해 벤젠 이량체에서 다양한 임계값과 기저함수에 따른 계산 비용 스케일링을 평가한다.
• 복합 방법의 성능을 vTight 및 vvTight 임계값에서의 전체 LNO-CCSD(T) 계산과 비교한다.

실험 결과

연구 질문

• RQ1 복합 LNO-CCSD(T) 방법이 S66 NCIs에서 SILVER 기준값과의 RMS 편차를 0.05 kcal/mol 이하로 유지할 수 있는가?
• RQ2 더 작은 기저함수에서 더 좁은 임계값 계산으로부터의 보정항이 vTight 또는 vvTight 임계값을 사용한 대규모 기저함수 계산의 필요성을 효과적으로 대체할 수 있는가?
• RQ3 정확도에 미치는 영향이 미미한 상황에서, 이러한 복합 방법에 카운터포인트 보정을 포함하는 것이 필수적인가?
• RQ4 최고의 복합 방법이 전체 vTight 또는 vvTight LNO-CCSD(T) 계산에 비해 계산 비용 측면에서 어느 정도의 이점이 있는가?
• RQ5 카운터포인트 보정이 비현실적인 인트라분자 NCIs에 대해 이러한 복합 방법을 신뢰성 있게 적용할 수 있는가?

주요 결과

• Tight{T,Q} + 0.60[vvTight – Tight]/T 복합 방법에 반-CP 보정을 적용한 결과, RMSD가 0.042 kcal/mol으로 나타나, LNO-CCSD(T,vTight)/haV{T,Q}Z 및 LNO-CCSD(T,vTight)/haV{Q,5}Z와 동일한 정확도를 확보하였다.
• 가장 저비용 복합 방법인 Normal{T,Q} + 0.95[vTight – Normal]/T에 반-CP 보정을 적용한 결과, RMSD가 0.063 kcal/mol로 나타나, 더 비싼 LNO-CCSD(T,Tight)/haV{Q,5}Z를 초월하는 성능을 보였다.
• 카운터포인트 보정은 정확도 향상에 거의 기여하지 않으며, 반-CP 방법이 전반적인 CP 보정과 유사한 성능을 보였고, 양자 모두 복합 방법의 성능에 열등하였다.
• 가장 정밀한 복합 방법(RMSD = 0.035 kcal/mol)은 LNO-CCSD(T,vTight)/haVQZ의 1.5배, LNO-CCSD(T,vvTight)/haV5Z의 3.7배 더 저렴한 비용으로 수행되었다.
• 벽 시계 스케일링 결과, 각 임계값 증가(예: Normal → Tight → vTight) 시 계산 비용이 약 3배 증가함을 확인하여 고임계값 계산의 높은 비용을 확인하였다.
• 제안된 cLNO 방법은 전체 vTight 또는 vvTight 계산에 비해 계산 비용을 최대 3.7배까지 절감하면서도, 0.05 kcal/mol 이하의 정확도를 유지한다.