프로레이트 구면좌표와 원자·동역학 균형 B 스플라인을 이용한 비분할 3차원 디랙 방정식 갈레르킨 해법

초록

본 논문은 원자·동역학 균형을 적용한 B‑스플라인 기저를 사용해 프로레이트 구면좌표계에서 두 중심 전자‑분자 시스템의 정적·동적 디랙 방정식을 비분할 갈레르킨 방식으로 풀고, 스펙트럼 오염을 억제하면서 높은 정확도의 에너지 스펙트럼과 외부 전자기장에 대한 파동함수 진화를 계산한다.

상세 분석

이 연구는 디랙 방정식의 수치해석에서 가장 난제 중 하나인 스펙트럼 오염(스퓨리어스 상태) 문제를 원자 균형(Atomic Balance)과 동역학 균형(Kinetic Balance)이라는 두 가지 균형 원리를 동시에 적용함으로써 효과적으로 해결한다. 기존의 연산자 분할 방법은 시간 의존 해석에 강점이 있지만, 초기 상태(바운드 및 연속 상태)를 정확히 구하기엔 수렴 속도가 매우 느리다. 저자들은 변분 원리를 기반으로 한 레일리‑라츠(Rayleigh‑Ritz) 접근을 채택하고, 큰 스핀오르와 작은 스핀오르 사이의 관계를 Eq.(20)과 같이 정확히 반영하는 균형 기저를 설계한다. 여기서 동역학 균형은 L_KB = (1/2mc²)α·p 로 정의되고, 원자 균형은 L_AB = (1/2mc²−V_c)α·p 로 정의된다. 이러한 균형 연산자를 사용하면, 제한된 차원 공간에서의 유한 행렬 근사에서도 스퓨리어스 고유값이 발생하지 않음이 정리 3.1·3.2에 의해 보장된다.

수치적 구현에서는 B‑스플라인 기저를 선택한다. B‑스플라인은 높은 차수와 국소성을 동시에 제공해 복잡한 구형 좌표계에서도 효율적인 행렬 조립을 가능하게 한다. 특히 프로레이트 구면좌표(ξ∈