[논문 리뷰] Multi-State Formulation of the Frozen-Density Embedding Quasi-Diabatization Approach

Patrick Eschenbach, Denis G. Artiukhin외 1명|arXiv (Cornell University)|2021. 09. 04.

Spectroscopy and Quantum Chemical Studies참고 문헌 98인용 수 10

한 줄 요약

이 논문은 냉각 밀도 임bedded 준다이아벳화(이하 FDE-diab) 방법의 다중 상태 확장을 제안하며, 대규모 분자 체계에서 스핀 밀도, 전자 결합 및 여기 상태 에너지를 정확하고 효율적으로 계산할 수 있도록 한다. 하위계 FDE 계산에서 유도된 준다이아벳 상태를 구성하고 일반화된 고유값 문제를 해결함으로써, GGA 함수형(PW91 등)을 사용하여 CASSCF 수준의 정확도를 확보하면서도 계산 비용을 크게 감소시킨다. 이로 인해 스핀 기대값 ⟨S²⟩ ≈ 0.8 a.u.을 얻고, 스핀 분포 오차는 1.5% 이하로 유지된다.

ABSTRACT

We present a multi-state implementation of the recently developed FDE-diab methodology [J. Chem. Phys., 148 (2018), 214104] in the Serenity program. The new framework extends the original approach such that any number of charge-localized quasi-diabatic states can be coupled, giving an access to calculations of ground and excited state spin-density distributions as well as to excitation energies. We show that it is possible to obtain results similar to those from correlated wave function approaches such as the complete active space self-consistent field method at much lower computational effort. Additionally, we present a series of approximate computational schemes, which further decrease the overall computational cost and systematically converge to the full FDE-diab solution. The proposed methodology enables computational studies on spin-density distributions and related properties for large molecular systems of biochemical interest.

연구 동기 및 목표

• 두 상태 FDE-diab 방법을 더 큰 분자 체계 및 더 복잡한 시스템에서 다중 준다이아벳 상태를 다룰 수 있도록 확장하기 위해.
• 분산된 스핀 밀도를 가진 체계에서 기초 상태 및 여기 상태의 스핀 밀도 분포, 전자 결합 및 여기 상태 에너지를 정확히 계산할 수 있도록 하기 위해.
• 완전한 FDE-diab 해에 체계적으로 수렴하는 근사 계산 기법을 개발하여 계산 비용을 감소시키기 위해.
• DNA 기저 삼중체 및 벤젠 올리거머 루비스탈 카이온에서 고수준 CASSCF 기준 데이터 및 KS-DFT와의 비교를 통해 방법의 타당성을 검증하기 위해.

제안 방법

• 방법은 이전 FDE 계산에서 구한 고정된 전자 밀도를 가진 하위계 슬래터 결정로 구성된 n개의 준다이아벳 상태를 직접 곱으로 구성한다.
• 각 아디아벳 상태는 이 n개의 준다이아벳 상태의 선형 조합으로 표현되며, 계수는 일반화된 고유값 문제 H B = S B E를 풀어 구한다.
• FDE 체계는 비가환 운동 에너지 포텐셜을 사용하여 전하 및 스핀 밀도를 국소화함으로써 표준 KS-DFT의 자기 상호작용 오차를 피한다.
• 각 하위계에 대해 분자 궤도를 계산하기 위해 고정 전자 밀도를 가진 콜린-쇼우 방정식(KSCED)을 사용한다.
• 하나의 하위계당 분자 궤도 수를 제한함으로써 근사 기법을 도입한다 (예: FDE-diab(i,8,8)), 이는 완전한 해에 대한 체계적 수렴을 가능하게 한다.
• 다양한 근사 프로토콜을 통합할 수 있는 통일된 표기법을 제안하여, 융통성 있고 비용 효율적인 계산을 가능하게 한다.

실험 결과

연구 질문

• RQ1 FDE-diab 방법은 초분자 체계에서 다수의 프래그먼트를 통해 스핀 분산을 기술할 수 있도록 두 개 이상의 준다이아벳 상태로 일반화될 수 있는가?
• RQ2 다중 상태 FDE-diab 접근법은 CASSCF 기준 데이터와 비교해 스핀 밀도, 스핀 분포 및 여기 상태 에너지를 얼마나 정확하게 재현하는가?
• RQ3 단순한 GGA 함수형(PW91)을 사용할 경우 방법의 성능은 어떠한가? 스핀 오염 및 정확도 측면에서 하이브리드 함수형과 비교해 볼 때 어떻게 되는가?
• RQ4 근사 계산 기법을 통해 비용을 얼마나 줄일 수 있으며, 동시에 정확도를 유지하면서 완전한 FDE-diab 해에 체계적으로 수렴할 수 있는가?
• RQ5 이 방법은 π-스택된 벤젠 올리고머 및 DNA 기저 삼중체와 같은 큰 생화학적으로 관련된 체계에 적용될 수 있는가?

주요 결과

• 다중 상태 FDE-diab 방법은 [G1TG2]•+에서 CASSCF 기준 데이터와 우수한 일치를 보이며, 스핀 분포 오차는 1.5% 이하, 여기 상태 에너지 오차는 0.01 eV 미만이다.
• GGA 함수형 PW91를 사용할 경우, 스핀 기대값 ⟨S²⟩ ≈ 0.8 a.u.을 도출하여 스핀 오염이 매우 미약함을 나타내며, 이는 하이브리드 함수형에서 관찰되는 높은 값(≈1.1 a.u.)과 대비된다.
• 벤젠 올리거머 루비스탈 카이온의 경우, FDE-diab(i,i,i)에서 i=6을 사용할 경우, 전체 FDE-diab(8,8,8) 기준값과 1% 이내의 차이로 스핀 밀도 분포를 재현하며 체계적 수렴을 입증한다.
• FDE-diab(4,8,8) 및 FDE-diab(8,4,8) 기법은 기준값과 1% 이내의 스핀 분포를 생성하지만, FDE-diab(4,4,4)는 최대 7%의 오차를 보이지만 i 증가에 따라 빠르게 수렴한다.
• FDE-diab(i,8,8), FDE-diab(8,i,8), FDE-diab(i,i,i)와 같은 근사 기법은 계산 비용을 크게 감소시키며, 완전한 해에 체계적으로 수렴할 수 있다.
• 이 방법은 표준 KS-DFT보다 뛰어나며, 특히 하이브리드 함수형을 사용할 경우 스핀 과도한 분산 및 높은 스핀 오염 문제를 악화시키는 것을 방지한다.