이중 여기 상태를 위한 아우프바우 억제 결합군 이론

초록

아우프바우 억제 결합군(ASCC) 방식을 이중 여기(excited) 상태에 적용하기 위해 파동함수 초기화 전략을 고안하고 구현하였다. 상태 평균 CASSCF를 시작점으로 하여, 가장 큰 구성들을 0차 파동함수에 정확히 재현하도록 설계했으며, 계산 비용은 기존 CCSD와 동등한 O(N⁵) 수준을 유지한다. 단일 이중 여기 구성에 의해 지배되는 상태와, 글리옥살 등에서 두 개의 주요 이중 여기 구성이 혼재하는 경우 모두 평균 0.15 eV, 최대 0.3 eV 정도의 오차를 보였다. 이는 EOM‑CCSD(4–6 eV)와 EOM‑CCSDT(0.4–0.8 eV)와 비교해 현저히 개선된 결과이다.

상세 분석

본 논문은 기존의 상태특이적 결합군 이론인 Aufbau Suppressed Coupled Cluster(ASCC)를 이중 여기(excited) 상태에 일반화하는 데 초점을 맞춘다. ASCC는 기본(ground‑state) Aufbau 결정자를 억제하고 목표 상태의 주요 구성만을 0차 파동함수에 남기기 위해, 탈여(De‑excitation) 연산자 S†를 이용해 두 개의 지수(exp) 연산자를 도입한다. 1‑CSF(단일 구성) 이중 여기의 경우, T⁽⁰⁾=c₁S + c₂S² 형태의 초기 클러스터 연산자를 선택하고, η와 c₁, c₂를 해석적으로 결정해 S²|Φ₀⟩가 목표 이중 여기 결정자를 정확히 재현하도록 만든다. 이때 η=√2, c₁=√2, c₂=0이 최적값이며, 결과적으로 0차 파동함수는 목표 상태와 동일해진다.

다중 구성(2‑CSF) 경우에는 두 개의 탈여 연산자 S₁₁†, S₁₂†와 가중치 γ, ω를 도입하고, T⁽⁰⁾에 S₁₁, S₁₁², S₁₂, S₁₂², S₁₁S₁₂의 선형 결합을 포함한다. 6개의 방정식을 풀어 γ, ω 및 c₁…c₅를 결정함으로써, 복합적인 이중 여기 파동함수를 정확히 초기화한다.

그 후, 전체 클러스터 연산자를 T = T_N + T_M + T_P 로 분해하고, T_P에 0차 초기값을, T_N·T_M은 1차 이하로 제한한다. 이는 H⁽⁰⁾와 H⁽¹⁾(Fock 연산자와 전자 상호작용의 0차·1차 부분)로 구성된 유사 변환 해밀리안 Ḣ를 이용해 에너지 및 진폭 방정식을 도출한다. 특히, 4개 이상의 비주요(non‑primary) 인덱스를 가진 혼합 진폭은 H⁽⁰⁾와 T⁽¹⁾의 교환항만이 기여함을 보이며, 이를 통해 O(N⁵) 비용으로 충분히 수렴 가능한 절단 스킴을 정당화한다.

계산 테스트에서는 상태 평균 CASSCF(활성 공간 최소화) 파동함수를 시작점으로 삼아, (i) 단일 이중 여기 구성으로 지배되는 베이징, 포름알데히드, 아세톤 등 소분자와, (ii) 글리옥살 및 유사 분자에서 두 개의 주요 이중 여기 구성이 혼재하는 경우를 조사했다. 결과는 평균 0.15 eV, 최악 0.30 eV의 오차를 보였으며, 이는 EOM‑CCSD(4–6 eV)와 EOM‑CCSDT(0.4–0.8 eV) 대비 크게 개선된 수치이다. 또한 Δ‑CC와 달리 다중 구성 상태에서도 일정 수준의 정확도를 유지한다는 점이 강조된다.

한계점으로는 현재 2‑CSF 이상 다중 구성에 대한 일반화가 미완성이며, 활성 공간 선택에 대한 민감도가 여전히 존재한다는 점이다. 향후 연구에서는 더 복잡한 다중 구성 이중 여기(예: 3‑CSF 이상)와 큰 시스템에 대한 효율적인 활성 공간 자동 선택, 그리고 삼중 및 사중 여기 효과를 포함한 고차 상관 효과를 결합하는 방안을 모색할 필요가 있다.