[논문 리뷰] Imaginary shift in CASPT2 nuclear gradient and derivative coupling theory

Jae Woo Park, Rachael Al-Saadon외 2명|arXiv (Cornell University)|2019. 04. 14.

Spectroscopy and Quantum Chemical Studies인용 수 4

한 줄 요약

이 논문은 가상 이동 정규화를 사용한 CASPT2에 대한 분석적 핵 구배 이론을 제안하며, 이는 이전의 실수 이동 구현을 확장한 것이다. 이 방법은 극히 소량의 계산 비용 증가로도 자극 상태 성질과 원뿔형 교차점 계산의 정확도를 향상시킨다. 주로 가상 이동 특화 밀도 행렬 항을 평가하는 데 기인한 계산 비용 증가가 발생하며, p-HBDI 및 FeP 모델을 통해 이를 입증하였다.

ABSTRACT

We report the analytical nuclear gradient theory for complete active space second-order perturbation theory (CASPT2) with imaginary shift, which is commonly used to avoid divergence of the perturbation expression. Our formulation is based on the Lagrangian approach and is an extension of the algorithm for CASPT2 nuclear gradients with real shift. The working equations are derived and implemented into an efficient parallel program. Numerical examples are presented for the ground- and excited-state geometries and conical intersections of a green fluorescent protein model chromophore, $p$-HBDI$^-$. We also report timing benchmarks with adenine, $p$-HBDI$^-$, and iron porphyrin. It is demonstrated that the energies and geometries obtained with the imaginary shift improve accuracy at a minor additional cost which is mainly associated with evaluating the effective density matrix elements for the imaginary shift term.

연구 동기 및 목표

• 내부 상태 발산을 방지하기 위해 가상 이동 정규화를 사용한 CASPT2에 대한 분석적 핵 구배 이론을 개발한다.
• 실수 이동의 한계를 해결하기 위해 더 견고하고 특이점이 없는 정규화 기법을 도입한다.
• 자극 상태 기하구조, 원뿔형 교차점, 수직/자극 상태 에너지를 정확하게 계산할 수 있도록 보장한다.
• 기존 CASPT2 구배 알고리즘에 가상 이동 항을 효율적으로 통합하여 계산 비용 증가를 최소화한다.
• p-HBDI 및 철 페오프린을 포함한 기준 시스템에서 메서드를 검증한다. 활성 공간 크기가 다양하게 설정된 경우를 대상으로 한다.

제안 방법

• 라그랑주 방법을 적응하여 가상 이동을 포함한 CASPT2에 대한 분석적 핵 구배를 유도하며, 이는 이전의 실수 이동 수식을 확장한 것이다.
• 에너지 분모의 정규화를 반영하기 위해 가상 이동 항에 대한 작동 수식을 유도하며, 이는 $\mathbf{d}^{(2)}_{\text{shift}}$ 및 $\tilde{y}_{I,M}^{\text{shift}}$를 포함한다.
• 병렬 계산을 활용하여 효율성을 높인 BAGEL 양자화학 패키지 내에 이 те오리 구현을 수행한다.
• 기존 코드에 최소한의 수정으로 $\lambda$-방정식 및 $Z$-벡터 방정식의 해법에 가상 이동 항을 통합한다.
• 특히 원뿔형 교차점 근처에서 중요한 특이점 방지를 위해 복소수 이동 $1/\Delta \to \Re[1/(\Delta + i\epsilon)]$을 사용한다.
• 새로운 항을 전통적 CASPT2 기여 항과 결합하여 계산 비용을 최적화하여, 큰 활성 공간에서는 비용 증가율을 약 15%로 줄였다.

실험 결과

연구 질문

• RQ1 CASPT2 핵 구배에서 가상 이동이 실수 이동에 비해 정확도를 어떻게 향상시키는가? 특히 원뿔형 교차점 근처에서의 성능을 중심으로 분석한다.
• RQ2 CASPT2 핵 구배 계산에 가상 이동 항을 포함할 경우 발생하는 계산 비용 증가율은 어느 정도인가?
• RQ3 활성 공간 크기가 증가함에 따라 가상 이동 방법의 성능 스케일링 특성은 어떻게 되는가?
• RQ4 실수 이동에 비해 가상 이동 방법이 이동 매개변수 $\epsilon$에 대해 얼마나 덜 민감한가?
• RQ5 가상 이동 수식은 자극 상태 기하구조 최적화 및 원뿔형 교차점 계산에서 수치적 안정성과 정확도를 유지할 수 있는가?

주요 결과

• 가상 이동 수식은 실수 이동에 비해 $\epsilon$의 변동에 덜 민감한 결과를 도출하여 자극 상태 계산의 견고성을 향상시킨다.
• 작은 활성 공간(예: 아데닌, p-HBDI)에서는 가상 이동 구배의 계산 비용이 실수 이동과 거의 동일하며, 추가 시간은 약 2–3%에 불과하다.
• 더 큰 활성 공간(예: CAS(10e,9o)를 가진 FeP)에서는 벽 시계 시간이 약 15% 증가(2947 s 대 3387 s)했으며, 주로 가상 이동 항 평가의 $O(N_{\text{act}}^9)$ 스케일링에 기인한다.
• FeP에서 가상 이동 항 $\mathbf{d}^{(2)}_{\text{shift}}$는 관련 밀도 행렬 계산 시간의 52%를 차지하며, 383초에 이르렀고, 반면 작은 시스템에서는 무시할 수 있을 정도로 낮은 시간(1초 미만)이 소요되었다.
• $\lambda$-방정식을 풀기 위한 추가 비용은 FeP에서 단지 11초 증가하여 반복 해법기의 영향을 최소화함을 시사한다.
• 이 방법은 성공적으로 공개된 BAGEL 양자화학 패키지에 통합되어 다구성 상태 양자화학 응용 분야에서 광범위한 사용을 가능하게 하였다.