k‑점 샘플링을 위한 SCDM‑k: 고체 전자구조의 국소 궤도 자동 생성법

초록

본 논문은 기존 SCDM 방법을 일반적인 k‑점 샘플링 체계에 확장한 SCDM‑k 알고리즘을 제안한다. 밀도 행렬의 선택된 열을 직접 이용해 게이지에 무관한 국소 궤도를 얻으며, 최적화 과정 없이 O(N log N) (N: k‑점 수) 복잡도로 정규 직교화된 Wannier‑like 함수들을 효율적으로 생성한다. 2·3 차원 모델 시스템을 통해 정확도와 스케일링을 검증하였다.

상세 분석

SCDM‑k는 전자 밀도 행렬이 실공간에서 지수적으로 감소한다는 ‘근접성 원리’를 이용한다. 기존 SCDM는 Γ‑점(단일 k‑점)에서만 적용 가능했으나, 고체에서는 Brillouin 구역을 다수의 k‑점으로 샘플링해야 한다. 여기서 핵심은 각 k‑점에 대해 얻어지는 Kohn‑Sham 궤도 ψ_{b,k}가 복소 위상(게이지)과, 경우에 따라는 유니터리 변환에 의해 자유롭게 변할 수 있다는 점이다. 기존의 Maximally Localized Wannier Functions(MLWF) 방식은 이러한 게이지를 사전에 ‘스무딩’해야 하며, 비선형 최적화 과정에서 지역 최소에 빠질 위험이 있다.

SCDM‑k는 먼저 전체 초셀(Ω_g)에서 모든 k‑점에 대한 밀도 행렬 P_k = Ψ_k Ψ_k† (Ψ_k는 n_b개의 밴드에 대한 파동함수 행렬)를 구성한다. 각 P_k는 주기성을 이용해 블록 토리 구조를 가지며, 전체 밀도 행렬 P = Σ_k P_k는 k‑점에 걸친 합으로 정의된다. 중요한 점은 P가 게이지 불변량이라는 사실이다; 즉, ψ_{b,k}에 임의의 위상·유니터리 변환을 적용해도 P는 변하지 않는다.

그 다음, 저차원 ‘로컬 슈퍼셀’ Ω_ℓ (N_ℓ1×N_ℓ2×N_ℓ3개의 유닛셀) 안에서 P의 열을 선택한다. 열 선택은 QR 분해에 컬럼 피벗을 적용하는 알고리즘으로 수행되며, 피벗된 열 인덱스 집합 C는 실질적인 국소 기저를 제공한다. 선택된 열은 원래 파동함수들의 선형 결합으로 표현되므로, 자동으로 Kohn‑Sham 불변 부분공간을 스팬한다.

선택된 열을 정규 직교화하면 최종적인 국소 궤도 φ_{μ}를 얻는다. 정규 직교화 단계는 FFT 기반의 컨볼루션을 이용해 O(N log N) 비용으로 수행된다. 여기서 N은 전체 k‑점 수이며, 로컬 슈퍼셀 크기 L_i는 정확도와 비용 사이의 트레이드오프 파라미터이지만, 실험적으로는 결과에 크게 민감하지 않다.

알고리즘의 복잡도는 크게 두 부분으로 나뉜다. (1) 비직교 국소 궤도 생성: 각 k‑점에 대해 P_k를 구성하고 QR 피벗을 수행하는 단계는 O(N)이다. (2) 직교화 및 전역 정규화: FFT를 이용한 전역 전파와 Gram‑Schmidt 과정이 O(N log N)으로 스케일한다. 따라서 수천에서 수만 개의 k‑점을 사용하는 대규모 고체 계산에서도 효율적으로 동작한다.

수치 실험에서는 2‑D와 3‑D 모델 포텐셜(예: 파라볼릭, 격자 포텐셜)에서 SCDM‑k가 얻은 국소 궤도가 전통적인 Wannier 함수와 거의 동일한 스프레드와 대칭성을 보였으며, 특히 게이지 스무딩이 필요 없는 점이 구현상의 편리함을 크게 향상시켰다. 또한, 복소 위상에 민감한 경우에도 안정적인 결과를 제공한다는 것이 강조된다.

요약하면, SCDM‑k는 (i) 게이지 독립성, (ii) 최적화 불필요, (iii) 선형·로그 선형 스케일링, (iv) 구현이 간단한 네 가지 장점을 갖는 새로운 국소 궤도 생성 방법이며, 대규모 주기적 시스템에 특히 적합한 도구이다.