와니에 함수

1. 개요

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

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

2. 정의

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

$$w_{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}$$

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

3. 위상 자유도와 최대 국소화

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

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

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

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

4. 어디에 쓰는가

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

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

5. 밀도범함수이론과의 관계

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

6. 관련 문서

• 밴드 구조 계산 · 평면파 기저
• 상태밀도
• 밀도범함수이론 · 콘-샴 방정식
• 고속 푸리에 변환
• 긴밀결합 모델

7. Footnotes

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

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

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

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