와니에 함수

편집 역사 토론
계산물리 고체역학 마지막 수정: 2026-07-23 04:38:26

1. 개요[편집]

와니에 함수(Wannier function)는 주기적 결정 안에서 넓게 퍼져 있는 밴드 구조 계산의 블로흐 상태(Bloch state)들을 한 밴드(또는 여러 밴드) 전체에 걸쳐 유니터리 푸리에 변환으로 결합해 얻은 실공간에 국소화된 궤도다. 블로흐 정리에 따르면 결정 속 전자 상태는 결정 전체에 퍼진 평면파 같은 함수인데, 화학자의 직관인 “원자에 붙은 국소 오비탈”과는 영 딴판이다. 와니에 함수는 이 둘 사이의 다리 — 퍼진 블로흐 상태를 국소화된 상으로 다시 묶어낸 것이다.1

핵심 발상은 위치-운동량이 푸리에 쌍인 것과 거의 같다. 블로흐 상태 ψnk(r)\psi_{n\mathbf{k}}(\mathbf{r})는 파수 k\mathbf{k}로 라벨된 퍼진 상태이고, 와니에 함수는 그 k\mathbf{k}들을 브릴루앙 영역 전체에 걸쳐 합성해 특정 격자점 R\mathbf{R}에 국소화시킨 실공간 상이다. 평면파 기저로 표현된 퍼진 상태와, 그 국소화된 짝 사이의 관계라고 보면 된다.

2. 정의[편집]

밴드 nn의 격자점 R\mathbf{R}에 놓인 와니에 함수는 블로흐 상태를 브릴루앙 영역(BZ)에서 적분해 정의한다.

wnR(r)=V(2π)3BZψnk(r)eikRdkw_{n\mathbf{R}}(\mathbf{r}) = \frac{V}{(2\pi)^3} \int_{\mathrm{BZ}} \psi_{n\mathbf{k}}(\mathbf{r})\, e^{-i\mathbf{k}\cdot\mathbf{R}}\, d\mathbf{k}

이렇게 만든 wnRw_{n\mathbf{R}}는 서로 다른 R\mathbf{R}nn에 대해 직교하며, 한 함수를 평행이동한 것이 다른 격자점의 함수가 된다. 즉 하나의 국소 궤도가 결정 전체에 격자 주기로 복제된 셈이다.

3. 위상 자유도와 최대 국소화[편집]

문제는 이 정의가 유일하지 않다는 것이다. 블로흐 상태에는 ψnkeiϕn(k)ψnk\psi_{n\mathbf{k}} \to e^{i\phi_n(\mathbf{k})}\psi_{n\mathbf{k}}라는 게이지 자유도(위상을 k\mathbf{k}마다 마음대로 곱해도 물리는 안 변함)가 있고, 여러 밴드가 얽히면 유니터리 회전 Umn(k)U_{mn}(\mathbf{k})까지 자유롭다. 이 게이지를 어떻게 고르느냐에 따라 와니에 함수의 모양 — 특히 얼마나 국소화되는지 — 가 완전히 달라진다. 게이지를 잘못 고르면 지수적으로 국소화되기는커녕 넓게 질질 퍼진 쓸모없는 함수가 나온다.

여기서 등장하는 게 최대 국소화 와니에 함수(maximally localized Wannier function, MLWF)다. 마르차리와 밴더빌트가 제안한 방법으로, 와니에 함수들의 총 퍼짐(spread)

Ω=n[wnr2wnwnrwn2]\Omega = \sum_n \left[ \langle w_n|r^2|w_n\rangle - |\langle w_n|\mathbf{r}|w_n\rangle|^2 \right]

을 최소화하도록 게이지 U(k)U(\mathbf{k})를 최적화한다.2 이렇게 얻은 MLWF는 화학 결합의 직관(결합 궤도, 고립쌍)과 놀랍도록 잘 대응한다. 실무에서는 Wannier90 같은 도구가 DFT 코드의 출력을 받아 이 최소화를 수행하는 것으로 널리 알려져 있다.

4. 어디에 쓰는가[편집]

와니에 함수가 인기 있는 이유는 결정 계산을 다루기 쉬운 국소 언어로 번역해주기 때문이다. 대표적 용도는 다음과 같다.

  • 긴밀결합 모델 구성: MLWF는 자연스러운 국소 기저라, 여기서 인접 격자점 간 겹침 적분을 뽑으면 그대로 긴밀결합(tight-binding) 모형의 호핑 파라미터가 된다. DFT의 방대한 밴드를 소수의 유효 궤도와 호핑으로 압축한 모형은 강상관계 등 후속 계산의 출발점이 된다.
  • 화학 결합 분석: DFT 밴드 구조는 퍼져 있어 “이 전자가 어느 결합에 있나”를 읽기 어렵지만, MLWF로 옮기면 결합 궤도·고립쌍 같은 익숙한 그림이 실공간에 그려진다. 상태밀도를 궤도별로 분해할 때도 유용하다.
  • 현대 분극 이론과 위상 밴드: 결정의 전기 분극은 사실 파동함수의 위상(베리 위상)에 담겨 있는데, 이것이 와니에 중심(Wannier center, 와니에 함수의 무게중심 위치)의 이동과 직접 연결된다. 위상 부도체의 지표를 와니에 중심의 흐름으로 추적하는 것도 이 덕분이다.3

5. 밀도범함수이론과의 관계[편집]

와니에 함수는 대개 밀도범함수이론이나 콘-샴 방정식 계산의 결과물인 블로흐 상태를 후처리해 얻는다. 즉 DFT가 밴드를 내놓으면, 그 밴드를 국소 궤도 언어로 재포장하는 사후 단계인 셈이다. 그래서 와니에 함수의 품질은 원 계산의 밴드 얽힘(entanglement) 처리와 게이지 최적화의 수렴에 크게 좌우된다.4 잘 수렴한 MLWF는 밴드 구조를 거의 완벽히 재현하는 초정확한 보간자로도 쓰여, 조밀한 k\mathbf{k}-격자 위 물성(예: 이상 홀 전도도)을 값싸게 계산하게 해준다.

6. 관련 문서[편집]

7. Footnotes[편집]

  1. 그레고리 와니에가 1937년에 도입했다. 당시엔 그저 형식적인 도구였는데, 21세기 들어 현대 분극 이론과 위상 물질 붐을 타고 계산물리의 주력 연장으로 화려하게 부활했다. 오래 묵힌 아이디어가 시대를 잘 만난 케이스.

  2. 퍼짐 Ω\Omega는 게이지에 안 변하는 부분(ΩI\Omega_I)과 변하는 부분(Ω~\tilde\Omega)으로 쪼개지는데, 최소화가 건드리는 건 후자뿐이다. 마르차리-밴더빌트(1997)의 이 분해가 MLWF 실용화의 핵심이었다.

  3. “전기 분극은 원래 잘 정의되지 않는다”는 충격적인 사실 — 무한 결정에서 분극은 절대값이 아니라 양자(quantum)만큼의 모호함을 가진 양이다. 이걸 와니에 중심(=베리 위상)으로 깔끔하게 정리한 게 킹-스미스·밴더빌트·레스타의 현대 분극 이론이다.

  4. 금속처럼 밴드가 다른 밴드와 얽혀 깔끔히 분리되지 않으면, 관심 에너지 창만 떼어내는 disentanglement 절차가 추가로 필요하다. 이 과정이 삐끗하면 국소화가 망가져 호핑 파라미터가 엉뚱하게 나온다.