[논문 리뷰] Analytic Gradients and Geometry Optimization for Orbital-Optimized Pair Coupled Cluster Doubles

Saman Behjou, Iulia Emilia Brumboiu외 1명|arXiv (Cornell University)|2026. 03. 20.

Advanced Chemical Physics Studies인용 수 0

한 줄 요약

분석적 기울기 기반의 기하 최적화 엔진을 PyBEST 내에서 궤도 최적화된 쌍 CC 이중식(OOpCCD/AP1roG)으로 도입하고, geomeTRIC 및 TRIC 좌표와 인터페이스하며 참조 기하에 대해 검증되었습니다.

ABSTRACT

We introduce a reusable geometry-optimization engine in PyBEST for analytic, gradient-driven molecular structure optimization, with particular emphasis on orbital-optimized pair coupled-cluster doubles (OOpCCD/AP1roG). The engine interfaces PyBEST with the exttt{geomeTRIC} optimizer, combining analytic electronic-structure gradients from PyBEST with the translation-rotation-internal coordinate (TRIC) framework, step control, and convergence machinery provided by exttt{geomeTRIC}. Specifically, we present the first implementation of analytic OOpCCD nuclear gradients within a Lagrangian formalism. Our approach and implementation are generally applicable to any seniority-zero wavefunctions that feature orbital optimization and allow for the evaluation of response one- and two-particle reduced density matrices. Owing to the seniority-zero structure of pCCD and the orbital stationarity of the optimized reference, the resulting gradient equations are compact, minimizing the storage of the full two-particle reduced density matrix, and avoiding finite-difference differentiation of wavefunction parameters. Validation on representative closed-shell systems shows that the OOpCCD-based PyBEST- exttt{geomeTRIC} workflow converges robustly and reproduces reference equilibrium geometries and energies within tight tolerances. Most importantly, OOpCCD produces structural parameters that deviate by approximately 0.02 Å (0.01 Å) for bond lengths or less than 1$^\circ$ for bond angles from CCSD(F12c)(T*) (MP2) reference structures.

연구 동기 및 목표

• 강한 상관 관계 시스템에서 궤도 최적화가 필수인 경우 정확한 기하 최적화를 목표로 한다.
• Lagrangian 프레임워크 내에서 OOpCCD를 위한 해석적 핵 기울기를 개발하고 구현한다.
• PyBEST 기울기를 geomeTRIC/TRIC와 결합한 재사용 가능한 기하 최적화 워크플로우를 제공하여 강건한 수렴을 달성한다.
• 선임성 0 pCCD 구조와 궤도 최적화로 인한 저장 공간 감소 및 안정적인 성능 이점을 보여준다.

제안 방법

• 관계식(Lagrangian, Z-vector) 접근법을 사용하여 명시적으로 궤도 반응 항을 제거한 OOpCCD의 해석적 기울기를 형식화한다.
• 그래디언트들을 응답 1-RDM 및 2-RDM과의 수축으로 표현하되, 이는 희소하고 짝-구조적이다.
• 비용과 희소성의 균형을 맞추기 위해 AO-MO 변환_basis에서 도함수 적분을 평가한다.
• 해밀토니안을 직교 MO(OMO) 기저로 변환하고 대칭 연결 도함수를 이용해 적분을 처리한다.
• TRIC 내부 좌표계와 geomeTRIC를 활용하여 그래디언트 기반 최적화 루프에서 기하 업데이트, 스텝 제어 및 수렴을 수행한다.
• cc-pVDZ/cc-pVTZ를 사용한 RHF 및 OOpCCD 기울기와 함께 이원자 및 소분자에 대해 벤치마크를 수행한다.

Figure 1: Molecular test set used for benchmarking the performance of the analytic-gradient-driven geometry optimization engine within OOpCCD. The molecular structures are taken from Ref. 25 and re-optimized for OOpCCD. All molecular structures are visualized using the Jmol package. 1

실험 결과

연구 질문

• RQ1 해석적 OOpCCD 기울기가 명시적 궤도 반응 항 없이 Lagrangian 프레임워크를 통해 구현될 수 있는가?
• RQ2 OOpCCD로 최적화된 기하가 MP2 및 CCSD(F12c)(T*)와 비교하여 결합 길이와 각을 얼마나 정확하게 재현하는가?
• RQ3 PyBEST–geomeTRIC 워크플로우가 OOpCCD 기울기를 이용해 바닥상태 및 전이상태 구조의 최적화에 대해 견고한가?
• RQ4 혼합 AO-MO 기울기 평가를 사용할 때 OOpCCD의 저장 및 계산 비용 절감과 수치적 안정성 이점은 무엇인가?
• RQ5 내부 좌표 최적화(TRIC)가 OOpCCD 맥락에서 연성인 분자 시스템의 경우 어떻게 작용하는가?

주요 결과

분자	기 basis 세트	수치 PES 적합: E(r_e) [Eh]	수치 PES 적합: r_e [Å]	해석적: E_e [Eh]	해석적: r_e [Å]
BN	cc-pVDZ	-79.029999	1.2688	-79.029999	1.2688
BN	cc-pVTZ	-79.065094	1.2582	-79.065094	1.2582
C2	cc-pVDZ	-75.549523	1.2387	-75.549522	1.2387
C2	cc-pVTZ	-75.582371	1.2213	-75.582369	1.2213
CN+	cc-pVDZ	-91.803418	1.1657	-91.803774	1.1630
CN+	cc-pVTZ	-91.842591	1.1499	-91.842591	1.1499
CO	cc-pVDZ	-112.855529	1.1231	-112.855529	1.1231
CO	cc-pVTZ	-112.911671	1.1156	-112.911671	1.1156
F2	cc-pVDZ	-198.855477	1.5186	-198.857066	1.5190
F2	cc-pVTZ	-198.949747	1.4624	-198.949747	1.4624
N2	cc-pVDZ	-109.062713	1.1016	-109.062706	1.1016
N2	cc-pVTZ	-109.127740	1.0867	-109.127740	1.0867

• 해석적 OOpCCD 기울기가 Lagrangian 프레임워크 내에서 올바르게 구현되어 정지점에서의 편도함수와 동일한 기울기를 산출한다.
• OOpCCD 기울기 평가가 희소하고 짝-구조의 1-RDM 및 2-RDM에 의존하며 도함수 적분과의 MO 기반 수축을 가능하게 한다.
• 혼합 AO-MO 접근은 1-체항은 AO 기저에서, 2-체항은 MO 기저에서 계산하여 저장과 비용의 균형을 이룬다.
• 검증 결과 OOpCCD 기하 최적화가 참조 평형 기하 및 에너지를 엄격한 공차 내에서 재현하며, 결합 길이는 대략 CCSD(F12c)(T*)/MP2 참조와의 차이가 보통 ~0.02 Å(일부 경우 0.01 Å) 이내이다.
• 분자-원자 간의 이차적 결합 길이는 해석적 최적화에서 PES 적합값과 일치하는 경향이 있어 그래디언트의 올바름을 신뢰하게 한다.
• 전이상태 최적화와 일반 바닥상태 기하가 PyBEST–geomeTRIC 인터페이스를 사용해 성공적으로 타깃팅된다.

Figure 2: RMSEs of bond lengths (Å, blue bars) and bond angles (degrees, orange bars) for OOpCCD equilibrium geometries relative to CCSD(F12c)(T*) reference structures, with all OOpCCD calculations performed using the cc-pVDZ basis set. The right vertical axis refers to angular deviations (degrees),