Brillouin zone 적분을 유한한 점과 점유 함수로 정확하게 근사하기

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

주기적 고체의 전자 밀도와 에너지는 Brillouin zone 전체 k 상태의 합입니다. 컴퓨터는 연속 적분 대신 유한한 점만 계산하며, 금속은 Fermi 준위에서 점유가 급격히 바뀌어 더 어렵습니다.

몇 개의 k점을 어디에 두고, 금속의 불연속 점유를 어떻게 다뤄야 할까요?

이번 장의 답k-grid는 reciprocal-space 밀도와 대칭으로 수렴시키고, 금속은 적절한 smearing으로 점유와 적분을 부드럽게 합니다. smearing 폭은 단순 수렴 도구이자 전자 자유에너지 항이므로 물성에 맞춰 외삽·보고해야 합니다.

01Monkhorst-Pack grid와 대칭

균일한 k-grid는 각 역격자 방향을 일정 간격으로 나누어 만듭니다. 결정 대칭과 시간반전 대칭으로 동등한 점을 묶으면 irreducible Brillouin zone의 대표점과 가중치만 계산해도 전체 적분을 재구성할 수 있습니다.

격자 크기 $N_1\times N_2\times N_3$는 셀 길이와 역비례하도록 잡아 각 방향의 k 간격을 비슷하게 만듭니다. slab의 진공 방향처럼 실제 분산이 없는 방향은 보통 1점을 사용합니다.

$$ Q=\sum_{\mathbf k}^{IBZ}w_{\mathbf k}Q(\mathbf k),\qquad \sum_{\mathbf k}w_{\mathbf k}=1 $$
가중치는 생략된 대칭 동등점의 몫을 포함합니다.

02금속이 더 촘촘한 표본을 요구하는 이유

0 K에서 상태 점유는 Fermi 준위 아래 1, 위 0으로 급격히 바뀝니다. k가 조금만 변해도 밴드가 Fermi surface를 가로질러 점유가 바뀌므로 적분 integrand가 매끄럽지 않고 에너지와 힘이 k-grid에 민감합니다.

smearing은 step 점유를 Gaussian, Fermi-Dirac, Methfessel-Paxton 같은 매끄러운 함수로 바꿉니다. SCF와 적분을 안정화하지만 폭이 너무 크면 전자구조와 에너지를 왜곡합니다.

방법주요 용도주의
Gaussian/Fermi-Dirac금속 SCF·유한 전자온도entropy 보정과 폭 수렴
Methfessel-Paxton총에너지·힘의 빠른 적분DOS 그림에 음의 점유 모양 가능
Tetrahedron수렴된 정적 DOS·에너지충분한 grid와 금속 보정 필요

03grid와 smearing을 함께 수렴시키기

한 가지 smearing 폭에서 k-grid만 늘리면 폭 오차를 놓칠 수 있고, 거친 grid에서 폭만 줄이면 Fermi surface 표본 오차가 커집니다. grid를 촘촘히 하면서 여러 폭에서 extrapolation 또는 free-energy correction의 안정성을 봅니다.

실무 판단 순서
  • 금속·반도체·절연체의 점유 성격을 먼저 확인
  • 셀 크기에 맞는 k 간격으로 grid 계열을 생성
  • 총에너지 차이·힘·응력·자기모멘트를 함께 비교
  • 최종 DOS는 더 촘촘한 grid와 적절한 적분법으로 별도 계산

04이 글의 용어 사전

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

k-point
Brillouin zone에서 Bloch 상태와 적분을 표본화하는 결정 운동량 점입니다.
Monkhorst-Pack grid
Brillouin zone을 대칭적인 균일 격자로 표본화하는 표준 k-point 생성법입니다.
irreducible Brillouin zone
결정 대칭으로 서로 동등한 k점을 제거하고 남은 최소 대표 영역입니다.
Fermi surface
0 K 금속에서 점유 상태와 비점유 상태를 가르는 k공간의 면입니다.
smearing
Fermi 경계의 불연속 점유를 유한 폭의 매끄러운 함수로 바꾸는 수치 기법입니다.
전자 엔트로피electronic entropy
부분 점유 상태가 만드는 엔트로피로, 유한온도 전자 자유에너지 또는 smearing 보정에 나타납니다.

05참고자료와 더 읽을거리

  1. Monkhorst & Pack, Special points for Brillouin-zone integrations (1976) — 주기적 고체의 k-point 격자를 체계적으로 만드는 표준 방법입니다.
  2. Methfessel & Paxton, High-precision sampling for Brillouin-zone integration (1989) — 금속의 부분 점유와 Brillouin zone 적분을 다루는 대표적인 smearing 방법입니다.
  3. Sholl & Steckel, Density Functional Theory: A Practical Introduction — 평면파 DFT의 이론·수치 설정·재료 응용·AIMD와 정확도를 잇는 이 과목의 중심 교재입니다.
NEXT DFT NOTESCF 반복, density mixing과 수렴 실패 →