시간 의존성 두 성분 상대론적 EA/IP 방정식운동 결합군법: TD‑EA/IP‑EOMCC

초록

본 논문은 정확한 두‑성분(X2C) 상대론적 프레임워크와 결합된 가상시간/실시간 접근법을 제시한다. 전자 부착(EA) 및 이온화 전위(IP) EOMCCSD 수준에서 개방‑쉘 원자(알칼리 금속·할로겐)의 선형 흡수 스펙트럼을 실시간 전기쌍극자 자기상관함수의 푸리에 변환으로 얻는다. 초기 상태는 코프만(det) 상태를 가상시간 전파해 얻으며, 저에너지 상태와의 겹침이 초기 전파 길이에 미치는 영향을 분석한다.

상세 분석

이 연구는 두 가지 핵심 기술을 결합한다. 첫째, EA‑EOMCC와 IP‑EOMCC를 비입자 보존 연산자로 정의함으로써, 제한된 Hartree‑Fock(RHF) 기준파에서 스핀 적합(open‑shell) 상태를 정확히 기술한다. 클러스터 연산자 (\hat T) 와 비보존 연산자 (\hat R^{\pm1}) 를 이용해 (N±1) 전자 시스템의 좌·우 고유벡터를 구하고, 이들 사이의 bi‑orthonormal 관계를 유지한다. 둘째, 이러한 고유벡터를 실시간 전파에 바로 사용할 수 있도록, 가상시간 전파(imaginary‑time propagation)를 도입한다. 코프만(det) 상태(예: 가장 높은 점유 궤도 (i) 또는 가장 낮은 비점유 궤도 (a))를 (\exp(-\bar H \tau)) 로 전파하면, (\tau)가 충분히 클 때 최저 에너지 EA/IP 상태에 수렴한다. 이 과정은 전통적인 좌·우 고유값 문제를 푸는 대신, 에너지와 정규화를 모니터링하며 자동으로 최적 초기 상태를 얻는다.

시간 의존성은 (\hat U(t)=e^{i\bar H_N t}) 로 정의되며, 여기서 (\bar H_N=\bar H-E^{(N\pm1)}0)는 (N±1) 전자 기준에 대해 에너지 이동된 유사 변환 해밀리안이다. 전기쌍극자 연산자 (\mu\alpha) 를 유사 변환한 (\bar\mu_\alpha) 를 이용해 순간 함수 (|M_{\alpha}^{(N\pm1)}\rangle = \bar\mu_\alpha \hat R^{\pm1}0|\Phi\rangle) 와 (\langle\tilde M{\alpha}^{(N\pm1)}| = \langle\Phi|\hat L^{\pm1}0 \bar\mu\alpha) 를 구성한다. 이 순간 함수는 EA/IP‑EOMCCSD의 ‘시그마 빌드’와 동일한 구조를 가지므로 기존 코드에 최소한의 수정만으로 초기 상태를 생성할 수 있다. 실시간 전파는 Runge‑Kutta와 같은 고정 스텝 통합법으로 수행되며, 전기쌍극자 자기상관함수 (\langle\tilde M(t) | M(0)\rangle) 를 얻어 푸리에 변환 후 Lorentzian 폭 (\gamma) 를 적용해 스펙트럼을 얻는다.

핵심 발견은 저에너지 상태와 코프만(det) 상태 사이에 비무시적 겹침이 존재할 경우, 가상시간 전파가 충분히 길어야 정확한 절대 에너지와 피크 강도를 재현한다는 점이다. Na, K, Rb (EA)와 F, Cl, Br (IP) 원자에 대해 X2C‑SVPall‑2c 기저를 사용했으며, 실시간 스펙트럼은 전통적인 주파수 영역 EA/IP‑EOMCCSD 결과와 거의 일치한다. 또한, 코어 전자와 고에너지 가상궤도를 적절히 동결(frozen‑core)함으로써 계산 비용을 크게 절감하면서도 스핀‑오비트 결합 효과를 정확히 포착했다.

이 논문은 (i) 비입자 보존 EOMCC를 실시간 TD‑EOMCC에 적용하는 방법론적 틀, (ii) 가상시간 전파를 통한 초기 상태 생성이라는 실용적 전략, (iii) 정확한 두‑성분 상대론적 X2C 프레임워크와 결합해 스핀‑오비트가 강한 원자 스펙트럼을 효율적으로 계산할 수 있음을 입증한다는 점에서 의미가 크다.