제5주기 원소를 위한 고정밀 상관 및 전체 에너지 벤치마크: ccECP와 CBS 한계

초록

본 연구는 Rb‑Xe 원소군에 대해 상관 일관성 유효핵전위(ccECP)를 사용하여 CI, CC, QMC 등 다양한 고급 전자구조 방법으로 총 에너지와 상관 에너지를 CBS 한계까지 정확히 추정한다. 또한 고정노드 DMC 오차와 핵심 전자 제거가 전자밀도에 미치는 영향을 kinetic energy 분석을 통해 정량화하였다.

상세 분석

본 논문은 제5주기 원소(Rb‑Xe)를 대상으로 ccECP가 제공하는 valence‑only Hamiltonian의 정확성을 다각도로 검증한다. 먼저, 각 원소에 대해 aug‑cc‑p(C)VnZ (n = D, T, Q, 5, 6) 기반의 상관 일관성 기반 Gaussian 기저함수를 구축하고, HF 에너지와 상관 에너지를 각각 식(4)와 식(5)의 지수·다항식 형태로 CBS 한계로 외삽한다. HF 외삽은 a·exp(‑bn) 형태로, 상관 에너지 외삽은 (n+3/8)⁻³·α + (n+3/8)⁻⁵·β 형태를 사용해 과대추정 위험을 최소화하였다.

전자상관 계산은 CISD, RCCSD(T)/UCCSD(T), CCSDT(Q), 그리고 가능한 경우 FCI까지 수행했으며, valence 전자수가 6 이하인 경우는 전부 FCI를 적용하였다. valence가 많아 CCSDT(Q) 계산이 불가능한 경우에는 UCCSD(T)와 CCSDT(Q) 사이의 상관 에너지 비율 η를 이용해 고차 기저에서의 CCSDT(Q) 에너지를 추정하였다. 이러한 절차는 전체 원소에 대해 1 % 이내의 오차로 valence 상관 에너지를 제공한다.

QMC 측면에서는 단일 결정자 Slater‑Jastrow 파동함수를 사용해 DMC 계산을 수행했으며, 시간 단계 τ = 0.02, 0.01, 0.005, 0.0025 Ha⁻¹에 대해 선형 외삽으로 τ→0 한계를 얻었다. T‑moves 알고리즘을 적용해 비국소 ccECP AREP의 로컬라이제이션 오류를 변분적으로 상한으로 처리하였다. 결과적으로 고정노드 오차는 원자 전자배치가 새로운 각운동량 채널(p, d 등)을 처음 채울 때 급격히 증가함을 확인했으며, 4d 전이금속(Y‑Cd)에서는 3d 대비 노드가 더 부드러워 FN 오차가 감소하는 경향을 보였다.

또한, CISD/FCI 파동함수에서 얻은 kinetic energy ⟨−½∇²⟩를 CBS 한계로 외삽함으로써 전자밀도 국소화 정도를 정량화했다. ECP는 virial 정리를 만족하지 않으므로 kinetic‑to‑potential 비율이 원소마다 크게 변동한다. 핵심 전자를 제거한 후에도 kinetic energy는 원자별 고유값을 갖으며, 이는 ccECP가 전자밀도 재배치를 어떻게 유도하는지를 보여준다. QMC에서 Jastrow 최적화가 진행될수록 kinetic energy가 감소하는 경향이 관찰돼, 변분적 최적화가 전자밀도 꼬리 부분을 과소평가할 수 있음을 시사한다.

전반적으로 본 연구는 ccECP와 상관 일관성 기저를 결합한 체계적 CBS 외삽이 높은 정확도의 원자 총 에너지와 상관 에너지를 제공함을 입증하고, DMC 고정노드 편향과 kinetic energy 변화를 통해 ECP가 전자밀도에 미치는 미세한 효과까지 정량화하였다. 이러한 데이터베이스는 향후 다원자 시스템에서 ccECP 기반 고정밀 계산을 수행할 때 신뢰할 수 있는 기준점으로 활용될 수 있다.