[논문 리뷰] Is valence CCSD(T) enough for the binding of water clusters? The isomers of (H$_2$O)$_6$ and (H$_2$O)$_{20}$ as a case study

Golokesh Santra, Margarita Shepelenko외 2명|arXiv (Cornell University)|2023. 08. 11.

Molecular spectroscopy and chiralityChemistry인용 수 3

한 줄 요약

이 연구는 (H₂O)₆ 및 (H₂O)₂₀ 이성질체에서 수소결합을 고려할 때, 벨런스 CCSD(T)의 적정성을 평가하기 위해 후기-CCSD(T) 전자상관효과를 분석한다. 핵심 결과로, 벨런스 CCSD(T)를 초월한 보정 중에서 핵-벨런스 상관이 가장 중요한 것으로 나타났으며, 고차항 전이(excitations)는 더 조밀한 이성질체를 선호하고, 상대론적 효과와 대각 Born-Oppenheimer 보정은 거의 상쇄된다.

ABSTRACT

Benchmark calculations on noncovalent interactions typically exclude correlation effects beyond valence CCSD(T) owing to their steep computational cost scaling. In this work, we consider their importance for water clusters, specifically, eight isomers of (H$_2$O)$_6$ and four Wales-Hodges isomers of (H$_2$O)$_{20}$. Higher order connected triples, $T_3$--(T), reduce dissociation energies of the latter by about 0.4 kcal/mol, but this is more than compensated by an increase of up to 0.85 kcal/mol due to connected quadruple excitations. In general, higher-order correlation effects favor more compact isomers over more `spread-out' ones. We also consider additional small effects for balance: scalar relativistics reduce binding in (H$_2$O)$_{20}$ by ca. --0.4 kcal/mol, which fortuitously is compensated by the ca. 0.55 kcal/mol diagonal Born-Oppenheimer correction. Core-valence correlation has the greatest impact, at ca. 1.3 kcal/mol for the icosamer.

연구 동기 및 목표

• 수소결합 클러스터에서 벨런스 CCSD(T)가 결합에너지 예측에 충분한지 평가하기 위해.
• 특히 연결된 삼중항(connected triples, T₃–(T)) 및 사중항 전이 효과가 (H₂O)₆ 및 (H₂O)₂₀의 이성질체 안정성에 미치는 영향을 정량화하기 위해.
• 핵-벨런스 상관, 스칼라 상대론적 효과, 대각 Born-Oppenheimer 보정이 비공명 결합에너지에 미치는 영향을 평가하기 위해.
• 더 큰 클러스터에서 후기-CCSD(T) 효과를 이중체 근사로 충분히 모델링할 수 있는지 판단하기 위해.

제안 방법

• 기준 CCSD(T) 계산은 MOLPRO 2022.3를 사용하여 수행되었으며, 완전 기저집합 보간법을 통해 cc-pVnZ 및 aug-cc-pVnZ 기저집합을 활용하였다.
• 후기-CCSD(T) 기여는 CFOUR 및 MRCC를 통해 CCSDT, CCSDT(Q), CCSDTQ(5) 방법을 사용하여 평가되었으며, CCSDT(Q) 계산에는 NCC 모듈을 활용하였다.
• 핵-벨런스 상관은 aug-cc-pwCVnZ 기저집합을 사용하여 평가하였고, 스칼라 상대론적 보정은 2차 상대론적 Douglas-Kroll-Hess 방법으로 계산하였다.
• 대각 Born-Oppenheimer 보정(DBOC)은 CFOUR 2.1을 사용하여 계산하였으며, 다체 전개(MBE) 기법을 활용해 기여도를 분석하였다.
• (H₂O)₆ 및 (H₂O)₂₀의 구조는 각각 BEGDB 및 WATER27 데이터셋에서 확보하였으며, 이성질체 간 에너지 차이를 분석하였다.
• 기저집합 수렴성 및 보간 파rameter는 Karton (2022) 및 참고문헌 43과 같은 기존 문헌에서 유도되었다.

실험 결과

연구 질문

• RQ1 연결된 삼중항 및 사중항 전이와 같은 고차항 전자상관 효과가 수소결합 클러스터 이성질체의 상대적 안정성에 어떤 영향을 미치는가?
• RQ2 벨런스 CCSD(T)를 초월한 핵-벨런스 상관 효과가 (H₂O)₆ 및 (H₂O)₂₀의 결합에너지에 얼마나 큰 영향을 미치는가?
• RQ3 스칼라 상대론적 효과와 대각 Born-Oppenheimer 보정은 이성질체 간 결합에너지 차이를 유의미하게 변화시키며, 상호 보완적인 영향을 미치는가?
• RQ4 후기-CCSD(T) 효과를 이중체 근사로 근사할 경우, 더 큰 클러스터에서의 대부분의 후기-CCSD(T) 기여를 잘 반영할 수 있는가?

주요 결과

• 고차항 연결 삼중항 전이(T₃–(T))는 (H₂O)₂₀의 해리 에너지를 약 0.4 kcal/mol 감소시키며, 연결 사중항 전이 효과는 최대 0.85 kcal/mol 증가시켜 더 조밀한 이성질체를 선호한다.
• 후기-CCSD(T) 벨런스 상관 효과의 순기여는 결합에너지 증가를 유도하며, 사중항 전이 효과는 프리즘 구조와 같은 조밀한 이성질체를 더 안정하게 만든다.
• 핵-벨런스 상관은 (H₂O)₂₀의 결합에너지에 1.25–1.36 kcal/mol 기여하며, 이는 벨런스 CCSD(T)를 초월한 보정 중 가장 큰 기여를 하며, 이는 이전의 벤젠 이량체와 비교해도 현저히 크다.
• 스칼라 상대론적 보정은 (H₂O)₂₀의 결합에너지를 약 0.4 kcal/mol 감소시키며, 이는 대략 0.56 kcal/mol 증가시키는 대각 Born-Oppenheimer 보정에 의해 거의 완전히 상쇄된다.
• 이중체 근사는 후기-CCSD(T) 기여의 90% 이상을 포괄하므로, 더 큰 클러스터에서 고차항 항을 이중체 전개로 현실적으로 모델링할 수 있음을 시사한다.
• 핵-벨런스 상관의 삼체 기여는 무시할 수 없을 정도로 작다(0.009–0.020 kcal/mol), 또한 스칼라 상대론적 효과 및 DBOC 보정의 삼체 및 사체 항 기여도는 유의미하지 않다.