SAP‑X2C: 간단하면서도 정확한 2‑성분 상대론적 해밀토니안

초록

SAP‑X2C는 Lehtola의 원자 전위 중첩(SAP) 모델을 이용해 2전자 그림변환(picture‑change) 효과를 포함한 X2C 해밀토니안을 제안한다. 1‑전자 X2C(1eX2C)의 저비용·구현 용이성을 유지하면서, 전자‑전자 상호작용에 대한 보정과 크기 집약적(picture‑change) 특성을 확보한다. 테스트 결과는 SAP‑X2C가 4‑성분 Dirac‑Hartree‑Fock(DCHF) 및 기존 AMF‑X2C 방법과 비교해 에너지·스핀‑오빗 분할·구조 파라미터에서 높은 정확도를 보이며, 특히 대규모·주기적 시스템에 적용하기에 적합함을 보여준다.

상세 분석

본 논문은 상대론적 전자구조 계산에서 2‑성분(2C) 방법의 핵심 과제인 2전자 그림변환(picture‑change, 2ePC) 문제를 해결하고자 한다. 기존 1eX2C는 1‑전자 연산자를 블록 대각화하여 스칼라 및 스핀‑오빗 효과를 저비용으로 얻지만, 2ePC를 무시하고 핵전위의 발산으로 인해 크기 집약성이 결여된다. 이를 보완하기 위해 저자들은 Lehtola가 제안한 SAP(Superposition of Atomic Potentials)를 도입한다. SAP는 원자별 유효 핵전하 함수를 가우시안 형태로 표현해 짧은 거리에서 급격히 감소하는 특성을 가지며, 이는 전자‑전자 상호작용을 원자 수준에서 근사적으로 포함한다.

SAP‑X2C는 1eX2C와 동일한 변환 행렬 U를 사용하되, 핵전위 V와 스몰 스핀오빗 행렬 W를 SAP 전위(V_SAP, W_SAP)로 교체한다. 여기서 V_SAP = V + V_e, W_SAP = W + W_e이며 V_e·W_e는 세 중심 2전자 가우시안 적분을 통해 계산된다. 변환 후에는 중복 계산을 방지하기 위해 V_e를 최종 해밀토니안에서 빼준다. 이 과정은 기존 1eX2C 코드에 3‑중심 적분 모듈만 추가하면 구현이 가능하므로, 구현 난이도와 계산 비용이 크게 증가하지 않는다.

이론적 측면에서 SAP‑X2C는 다음 두 가지 중요한 장점을 제공한다. 첫째, 전자 전위가 원자별로 제한적이면서도 초다항식 감소(super‑polynomial decay)를 만족하므로, 무한 격자에 대한 절대 수렴성을 보장한다. 이는 주기적 시스템이나 대규모 클러스터에 적용할 때 필수적인 속성이다. 둘째, SAP 전위에 포함된 전자‑전자 스크리닝 효과가 2ePC를 근사적으로 반영하므로, 스핀‑오빗 상호작용과 스칼라 상대론적 효과가 동시에 개선된다.

계산 실험에서는 다양한 경량·중량 원소를 포함한 분자군과 Og₂와 같은 초중량 이원자를 대상으로 HF 에너지, 스핀오빗 에너지, 핵심 껍질의 스핀‑오빗 분할, 평형 결합 거리 및 조화 진동수를 비교하였다. 결과는 1eX2C가 Z가 커질수록 오차가 급증하는 반면, SAP‑X2C는 전반적으로 수 mEh 수준의 오차로 억제된다. 특히 Og₂의 바인딩 스핀오빗 레벨에서는 SAP‑X2C가 1eX2C 대비 1~2 orders of magnitude 작은 오차를 보였다. AMF‑X2C와 eamfX2C는 전반적으로 더 높은 정확도를 보이지만, 이들 방법은 원자별 4C Dirac‑Hartree‑Fock 계산을 사전 수행해야 하는 비용이 크다. 반면 SAP‑X2C는 사전 원자 계산이 필요 없으며, 구현이 간단하고 블랙박스 형태로 대규모 계산에 바로 적용할 수 있다.

크기 집약성 테스트에서는 Xe 결정 조각을 점진적으로 확대하면서 총 에너지 오차를 조사했으며, SAP‑X2C는 시스템 크기가 증가함에 따라 오차가 선형적으로 감소하는 이상적인 스케일링을 보였다. 이는 전자 전위가 원자별로 국소화되어 전체 시스템에 대한 그림변환이 독립적으로 적용될 수 있음을 의미한다.

요약하면, SAP‑X2C는 1eX2C의 저비용·단순성을 유지하면서 2ePC와 크기 집약성을 동시에 해결한 최적의 2‑성분 상대론적 해밀토니안이다. 기존 AMF‑X2C 계열과 비교해 약간의 정확도 차이는 존재하지만, 구현 용이성, 사전 계산 필요성 감소, 그리고 대규모·주기적 시스템에 대한 적용 가능성 측면에서 큰 장점을 제공한다.