평면파 기저와 cutoff energy | 조우현의 연구 노트

무한 Fourier 급수를 유한 기저로 자르고 정확도와 계산비를 수렴시키기

DFT · PART 3 · 연속적인 양자 방정식을 유한한 계산으로 바꾸는 동안03 / 07

Kohn-Sham 방정식은 연속 공간의 미분방정식입니다. 컴퓨터는 유한한 숫자만 다룰 수 있으므로 오비탈을 유한한 기저의 계수로 바꾸어야 합니다.

평면파를 몇 개까지 포함해야 하며 cutoff를 높인다는 것은 공간적으로 무엇을 뜻할까요?

이번 장의 답주기 함수는 역격자 평면파로 전개하고 $|\mathbf k+\mathbf G|^2/2\le E_{cut}$인 성분만 남깁니다. cutoff가 높을수록 짧은 길이 변화를 표현하지만 기저 수와 FFT 비용이 증가하므로 관심 물성을 직접 수렴시켜야 합니다.

CONTINUOUS ORBITAL → FINITE COEFFICIENTS

01Bloch 오비탈$\psi_{n\mathbf k}$

02평면파 전개$\sum_G c_Ge^{i(k+G)r}$

03cutoff 구$E_G\le E_{cut}$

$\frac{|\mathbf k+\mathbf G|^2}{2}\le E_{cut}$

cutoff는 역공간에서 포함할 평면파 구의 반지름이자 실공간에서 표현 가능한 최소 길이의 척도입니다.

01왜 고체에서 평면파가 편리한가

평면파는 주기 경계조건을 자동으로 만족하고 서로 직교하며, 운동에너지 연산자가 대각 형태가 됩니다. 원자 위치가 바뀌어도 기저 자체가 원자에 붙어 움직이지 않아 힘 계산이 비교적 체계적이고 기저를 cutoff 하나로 수렴시킬 수 있습니다.

$$ \psi_{n\mathbf k}(\mathbf r)=\frac{1}{\sqrt\Omega}\sum_{\mathbf G}c_{n\mathbf k}(\mathbf G)e^{i(\mathbf k+\mathbf G)\cdot\mathbf r} $$

$\Omega$는 셀 부피이고 계수 $c$가 수치적으로 구할 미지수입니다.

단점은 핵 근처에서 전전자 파동함수가 빠르게 진동해 매우 많은 평면파가 필요하다는 점입니다. pseudopotential과 PAW가 이 급격한 변화를 효율적으로 처리합니다.

02kinetic-energy cutoff와 기저 크기

평면파의 운동에너지는 $|\mathbf k+\mathbf G|^2/2$이므로 일정 에너지 아래 성분만 포함하는 자연스러운 절단을 사용합니다. 3차원에서 cutoff가 커질수록 포함되는 평면파 수는 대략 $E_{cut}^{3/2}$에 비례해 증가합니다.

전자 오비탈 cutoff와 전하 밀도 FFT grid cutoff가 별도로 있을 수 있습니다. 밀도는 오비탈 곱이므로 더 높은 Fourier 성분을 포함하며 ultrasoft·PAW augmentation charge는 더 촘촘한 grid가 필요할 수 있습니다.

cutoff 증가	좋아지는 것	비용·주의
오비탈 cutoff	파동함수·힘·응력 표현	행렬 크기와 FFT 증가
밀도 cutoff/grid	전하·퍼텐셜의 짧은 파장	메모리와 통신 증가
초격자 부피	같은 cutoff에서 평면파 수 증가	원자 수보다 더 빠른 비용 증가 가능

03에너지보다 관심량을 수렴시키기

cutoff를 단계적으로 높이며 총에너지/atom만 보는 것은 시작일 뿐입니다. 상간 에너지 차이, force, stress, phonon처럼 미분에 민감한 양은 더 높은 cutoff가 필요할 수 있습니다. 서로 다른 원소를 비교할 때는 가장 단단한 원소의 데이터셋 요구치를 기준으로 봅니다.

권장 cutoff 검사

구조를 고정하고 같은 k-grid에서 cutoff만 변화
관심 에너지 차이와 최대 힘·응력 변화를 기록
production 값보다 한 단계 높은 cutoff에서 결론이 유지되는지 확인
pseudopotential 종류를 바꾸면 검사를 처음부터 다시 수행

04이 글의 용어 사전

본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.

기저함수basis function: 미지의 파동함수를 선형결합으로 표현할 때 사용하는 알려진 함수입니다.
평면파plane wave: 일정한 파동벡터를 갖는 복소 지수 함수로 주기계 Fourier 기저를 이룹니다.
kinetic-energy cutoff: 평면파 기저에 포함할 최대 운동에너지를 정하는 절단값입니다.
FFTfast Fourier transform: 실공간과 역공간 표현을 빠르게 변환하는 알고리즘입니다.
augmentation charge: ultrasoft·PAW에서 핵 근처의 정확한 전하 정보를 복원하기 위한 국소 보정 전하입니다.
Pulay stress: 셀 크기가 변할 때 유한 평면파 기저가 불연속적으로 바뀌며 생기는 인공 응력 오차입니다.

05참고자료와 더 읽을거리

Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
Blöchl, Projector augmented-wave method (1994) — 평면파 효율성과 핵 주변의 전전자 정보를 잇는 PAW 방법의 원 논문입니다.
Perdew, Burke & Ernzerhof, Generalized Gradient Approximation Made Simple (1996) — 재료 계산에서 널리 쓰이는 PBE-GGA 교환-상관 범함수의 원 논문입니다.

NEXT DFT NOTE핵심전자, pseudopotential과 PAW →