란초스 알고리즘을 활용한 디랙 방정식 해법

초록

본 논문은 외부 쿨롱 퍼텐셜 하에서 상대론적 전자를 기술하는 디랙 방정식의 고유값을 란초스 알고리즘으로 구한다. 반복적으로 삼중대각 행렬을 구축하고 각 단계마다 대각화를 수행해 수렴하는 고유값을 얻으며, 연속 스펙트럼에 의해 발생하는 가짜 해를 효과적으로 식별한다.

상세 분석

디랙 방정식은 스핀‑½ 입자의 상대론적 동역학을 기술하는 4×4 행렬 미분 방정식으로, 연속 스펙트럼과 바운드 스테이트가 동시에 존재한다는 특성 때문에 수치적 고유값 해법이 까다롭다. 저자들은 이러한 문제에 란초스(Lanczos) 알고리즘을 적용함으로써 효율적인 삼중대각(tridiagonal) 행렬 표현을 얻는다. 란초스 과정은 초기 벡터를 선택한 뒤 해밀토니안 작용을 반복 적용해 정규 직교 기저를 생성하고, 이때 생성되는 삼중대각 행렬은 원래 무한 차원의 디랙 해밀토니안을 저차원 부분공간에 투사한 형태가 된다. 각 반복 단계마다 이 삼중대각 행렬을 전형적인 QR 알고리즘 등으로 대각화하면 현재까지 구축된 부분공간에 대한 근사 고유값과 고유벡터를 얻을 수 있다.

특히 저자들은 “반복마다 대각화” 전략을 채택해, 수렴 판정 기준을 고유값 변화량과 고유벡터 정규화 오차로 설정하였다. 바운드 상태(예: 수소 원자와 같은 쿠롱 포텐셜)에서는 고유값이 급격히 수렴하는 반면, 연속 스펙트럼에 해당하는 영역에서는 고유값이 불안정하게 움직이며 가짜 해(spurious solutions)를 만든다. 이러한 가짜 해는 란초스 과정에서 발생하는 수치적 오버플로우나 정규화 손실, 그리고 연속 상태를 부분공간에 강제로 포함시키면서 생기는 비물리적 모드이다. 저자들은 고유벡터의 정규화 정도와 기대값 ⟨ψ|H|ψ⟩의 실수성, 그리고 파동함수의 급격한 진동 패턴을 분석함으로써 가짜 해를 자동으로 구분하는 알고리즘을 제시한다.

또한, 삼중대각 행렬의 크기를 점진적으로 늘려가며 메모리 사용량을 최소화하고, 각 단계에서 필요한 연산량을 O(N) 수준으로 유지한다는 점이 큰 장점이다. 이는 전통적인 직접 대각화 방법이 O(N³) 복잡도를 갖는 것과 대비된다. 저자들은 수치 실험을 통해 1s, 2p 등 저에너지 바운드 상태는 10~12자리 정확도로 재현되며, 고에너지 연속 영역에서도 스펙트럼 밀도와 페르미 골든 레일을 정확히 추정한다는 결과를 보고한다.

이와 같은 접근법은 디랙 방정식뿐 아니라, 양자 전기역학(QED)에서의 외부 전자기장 문제, 그래핀과 같은 2차원 Dirac 물질, 그리고 핵물리학에서의 중성자-핵 상호작용 등 다양한 분야에 확장 가능성을 시사한다.