밴드 구조 계산

양자화학 계산물리 마지막 수정: 2026-07-06 04:37:14

밴드 구조 계산 Band Structure Calculation
분야	고체물리 × 양자화학 × 계산물리
핵심 원리	블로흐 정리, 역격자, 브릴루앙 영역
주요 방법	밀도범함수이론(DFT), 평면파 + 유사퍼텐셜
대표 소프트웨어	VASP, Quantum ESPRESSO, WIEN2k

1. 개요[편집]

결정 하나가 금속인지 반도체인지 절연체인지는, 이 그림 한 장의 밴드갭이 있느냐 없느냐로 갈린다.

밴드 구조 계산(Band Structure Calculation)은 주기적 결정 안에서 전자가 가질 수 있는 에너지 준위를 파동벡터 $\mathbf{k}$ 의 함수 $E_n(\mathbf{k})$ 로 구하는 계산 기법이다. 이 에너지-운동량 관계, 즉 밴드 구조는 물질이 전기를 얼마나 잘 통하는지, 빛을 어떤 색으로 흡수·방출하는지, 열전 특성이 어떤지 같은 거의 모든 전자적 성질의 뿌리다.

원자 하나에서는 전자가 이산적인 에너지 준위를 갖지만, 원자가 아보가드로 수만큼 규칙적으로 모여 결정을 이루면 이 준위들이 폭을 가진 “띠(band)“로 퍼진다. 채워진 띠(valence band)와 빈 띠(conduction band) 사이의 금지된 에너지 간격이 바로 밴드갭이다. 반도체 산업 전체가 이 밴드갭 하나를 손보는 일에 매달려 있다고 해도 과언이 아니다.¹

2. 블로흐 정리[편집]

밴드 구조 계산의 출발점은 블로흐 정리(Bloch’s theorem)다. 주기 포텐셜 $V(\mathbf{r}) = V(\mathbf{r} + \mathbf{R})$ 안에서 슈뢰딩거 방정식의 해는 다음 형태를 가진다.

\psi_{n\mathbf{k}}(\mathbf{r}) = e^{i\mathbf{k}\cdot\mathbf{r}}\, u_{n\mathbf{k}}(\mathbf{r}), \qquad u_{n\mathbf{k}}(\mathbf{r}) = u_{n\mathbf{k}}(\mathbf{r} + \mathbf{R})

즉 파동함수는 평면파에 격자와 같은 주기를 갖는 함수 $u$ 를 곱한 꼴이다. 이 정리 덕분에 무한히 큰 결정 전체를 풀 필요 없이, 단위 세포 하나만 풀고 파동벡터 $\mathbf{k}$ 로 색인하면 된다. 무한을 유한으로 바꿔 주는 이 마법이 고체물리 계산을 가능하게 하는 핵심이다.

3. 역격자와 브릴루앙 영역[편집]

실공간 격자 벡터 $\mathbf{a}_i$ 에 대응해, 역공간에는 $\mathbf{b}_i \cdot \mathbf{a}_j = 2\pi\,\delta_{ij}$ 를 만족하는 역격자 벡터 $\mathbf{b}_i$ 가 정의된다. 파동벡터 $\mathbf{k}$ 는 이 역공간에 산다. 역격자의 위그너-자이츠 세포가 곧 제1 브릴루앙 영역(first Brillouin zone)이며, 물리적으로 서로 다른 모든 $\mathbf{k}$ 는 이 영역 안에 다 담긴다. 밴드 구조를 그릴 때 흔히 보는 가로축(Γ–X–M–Γ 같은 경로)은 이 브릴루앙 영역의 고대칭점들을 잇는 여정이다.

3.1. k-점 샘플링[편집]

물리량(전자 밀도, 전체 에너지 등)은 브릴루앙 영역 전체에 대한 적분으로 얻어지는데, 이를 유한 개의 $\mathbf{k}$ -점에서의 합으로 근사한다. Monkhorst-Pack 격자가 표준적인 샘플링 방식이다. $\mathbf{k}$ -점을 촘촘히 할수록 정확하지만 계산 비용이 비례해 늘어난다. 특히 금속은 페르미 면 근처에서 점유가 급변해 매우 조밀한 샘플링을 요구하는 반면, 밴드갭이 큰 절연체는 성긴 격자로도 수렴한다.²

4. 기저와 유사퍼텐셜[편집]

$u_{n\mathbf{k}}$ 를 수치적으로 표현하려면 기저 함수가 필요하다. 가장 널리 쓰이는 것은 평면파(plane-wave) 기저로, 주기성과 궁합이 좋고 정확도를 컷오프 에너지 하나로 체계적으로 조절할 수 있다는 장점이 있다.

문제는 원자핵 근처에서 파동함수가 격렬하게 진동해, 이를 평면파로 표현하려면 어마어마한 수의 평면파가 필요하다는 것이다. 그래서 화학 결합에 거의 관여하지 않는 내각(core) 전자를 매끄러운 유효 포텐셜로 대체하는 유사퍼텐셜(pseudopotential) 기법을 쓴다. 이렇게 하면 필요한 평면파 수가 극적으로 줄어든다. 전각(all-electron)을 그대로 다루는 LAPW 계열(WIEN2k)은 더 정확하지만 그만큼 비싸다.

5. 계산 방법과 DFT의 한계[편집]

밴드 구조 계산의 실질적 표준은 밀도범함수이론(DFT)이다. DFT는 다전자 문제를 유효 단일입자 문제(Kohn-Sham 방정식)로 환원해, 자기무모순장(self-consistent field) 반복으로 밴드를 얻는다. 계산이 감당 가능하다는 실용성 덕분에 물질 스크리닝의 주력이 되었다.

다만 유의할 함정이 하나 있다. 표준 DFT(LDA, GGA)는 밴드갭을 체계적으로 과소평가한다. 실리콘의 실측 밴드갭은 1.1 eV인데 LDA는 0.5 eV 안팎을 내놓는 식이다.³ 이는 Kohn-Sham 고유값을 참된 준입자 에너지로 오해한 데서 오는 구조적 한계로, 계산을 촘촘히 한다고 해결되지 않는다. 더 정확한 밴드갭이 필요하면 하이브리드 범함수(HSE)나 슈뢰딩거 방정식의 다체 이론인 GW 근사, 혹은 고전적 하트리-폭 방법 기반 접근을 동원해야 한다.

6. 주요 소프트웨어[편집]

VASP — 평면파 + 유사퍼텐셜(PAW). 상용이지만 물질 계산 분야의 사실상 표준.
Quantum ESPRESSO — 오픈소스 평면파 코드. 라이선스 부담 없이 널리 쓰인다.
WIEN2k — LAPW 전각 방식. 자화·초정밀 계산에서 신뢰받는다.

7. 관련 문서[편집]

밴드갭이 없으면 금속, 있으면(대략 4 eV 이하) 반도체, 크면(4 eV 이상) 절연체로 대략 나뉜다. 물론 경계는 흐릿하고, 반도체와 절연체의 구분은 사실 정도의 문제다. ↩
금속의 페르미 면 적분을 안정화하려고 점유수를 인위적으로 부드럽게 퍼뜨리는 스미어링(Fermi/Gaussian/Methfessel-Paxton) 기법을 쓴다. 온도가 없는 계에 가짜 온도를 넣는 셈인데, 결과 에너지는 이 스미어링 폭을 0으로 외삽해 보정한다. ↩
이 “밴드갭 과소평가 문제”는 DFT를 처음 배우는 이들이 반드시 데는 통과의례다. LDA 밴드갭을 실측인 줄 알고 논문을 쓰다가 심사자에게 혼나는 것이 계산 재료과학의 오랜 전통이다. ↩