밀도범함수이론

양자화학 계산물리 마지막 수정: 2026-07-03 22:03:51

밀도범함수이론 Density Functional Theory
약칭	DFT
분야	양자화학 × 계산물리 × 재료과학
기본 변수	전자 밀도 $n(\mathbf{r})$
이론적 기반	호엔베르크-콘 정리
대표 소프트웨어	VASP, Quantum ESPRESSO

1. 개요[편집]

파동함수는 3N차원이라 못 푼다. 그런데 밀도는 3차원이잖아? — 이 한 문장으로 노벨상을 받았다.

밀도범함수이론(Density Functional Theory, DFT)은 다전자계의 성질을 파동함수 대신 전자 밀도 $n(\mathbf{r})$ 를 기본 변수로 삼아 기술하는 양자역학 이론이다. 슈뢰딩거 방정식을 직접 풀려면 전자 $N$ 개짜리 계에서 $3N$ 차원 파동함수를 다뤄야 하는 차원의 저주에 부딪히지만, 전자 밀도는 전자가 몇 개든 언제나 3차원 공간의 함수다.¹ 이 발상의 전환이 계산 비용을 천문학적으로 낮췄다.

오늘날 재료과학, 촉매화학, 신약설계, 반도체 물성 예측 등에서 가장 널리 쓰이는 제일원리(first-principles) 계산 도구이며, 논문 인용 통계에서 물리·화학 분야를 통틀어 압도적 상위권을 차지한다. “이론적으로는 엄밀하고, 실제로는 근사에 의존하는” 묘한 이중성이 이 이론의 매력이자 함정이다.

2. 호엔베르크-콘 정리[편집]

DFT의 수학적 정당성은 1964년 피에르 호엔베르크와 월터 콘이 증명한 두 개의 정리에서 나온다.

제1정리 — 바닥상태의 전자 밀도 $n(\mathbf{r})$ 는 외부 퍼텐셜 $V_{\text{ext}}(\mathbf{r})$ 를 (상수 차이를 빼면) 유일하게 결정한다. 즉 밀도만 알면 원리적으로 계의 모든 것을 알 수 있다. 파동함수라는 어마어마한 정보량이 사실은 3차원 밀도 하나로 압축 가능하다는 충격적인 주장이다.²
제2정리 — 에너지를 밀도의 범함수 $E[n]$ 로 쓸 수 있으며, 진짜 바닥상태 밀도에서 이 범함수가 최솟값을 갖는다. 즉 변분 원리가 성립한다.

문제는 이 정리들이 “그런 범함수가 존재한다”는 것만 보장할 뿐, 그게 어떻게 생겼는지는 알려주지 않는다는 것. 존재는 하는데 정체는 모르는, 수학자들이 좋아하는 그런 종류의 결과다.

3. 콘-샴 방정식[편집]

1965년 월터 콘과 루제이 샴은 실용적인 돌파구를 냈다. 상호작용하는 진짜 전자계를, 같은 밀도를 내놓는 상호작용 없는 가상의 보조계로 치환하는 것이다. 이 보조계의 각 전자는 유효 퍼텐셜 $V_{\text{eff}}$ 속에서 움직이는 단일입자 방정식을 따른다:

\left[ -\frac{\hbar^2}{2m}\nabla^2 + V_{\text{eff}}(\mathbf{r}) \right] \psi_i(\mathbf{r}) = \varepsilon_i \psi_i(\mathbf{r})

여기서 밀도는 채워진 궤도들의 합으로 재구성된다:

n(\mathbf{r}) = \sum_{i}^{\text{occ}} |\psi_i(\mathbf{r})|^2

그런데 $V_{\text{eff}}$ 가 밀도 $n$ 에 의존하고, $n$ 은 다시 방정식을 풀어야 나온다. 그래서 콘-샴 방정식은 자기무모순장(self-consistent field, SCF) 반복으로 푼다. 초기 밀도를 찍고, 방정식을 풀어 새 밀도를 얻고, 수렴할 때까지 뺑뺑이를 돈다. 전산유체역학에서 잔차 그래프가 내려가길 기도하는 것과 정확히 같은 정신적 고통이 여기에도 있다.

4. 교환상관 범함수[편집]

콘-샴의 트릭에는 대가가 있다. 진짜 계와 보조계의 에너지 차이 — 전자들이 서로를 피하며 생기는 교환(exchange)과 상관(correlation) 효과 — 이 전부 하나의 항 $E_{xc}[n]$ 에 밀어넣어진다. 그리고 이 교환상관 범함수의 정확한 형태는 아무도 모른다.³ DFT 계산의 정확도는 사실상 이 범함수를 얼마나 잘 근사하느냐에 전부 달려 있다.

근사 등급	이름	사용하는 정보	특징
1단계	LDA	국소 밀도만	단순, 결합에너지 과대평가
2단계	GGA	밀도 + 밀도 구배	PBE 등, 가성비 표준
상위	하이브리드	GGA + 정확한 교환 일부	B3LYP 등, 정확하지만 비쌈

이른바 “야코프의 사다리(Jacob’s ladder)“라 불리는 이 위계는 위로 올라갈수록 정확해지지만 비용도 오른다. 그리고 어떤 사다리를 올라도 밴드갭 과소평가, 반데르발스 힘 부재, 강상관계 전자 실패 같은 고질병에서 완전히 자유롭지는 못하다. DFT를 “제일원리”라 부르지만, 교환상관 범함수라는 반쯤 경험적인 요소가 섞여 있어 순수주의자들은 이 명칭에 눈살을 찌푸리기도 한다.

5. 계산 비용과 소프트웨어[편집]

콘-샴 방정식을 푸는 표준 구현의 계산 비용은 대략 전자 수 $N$ 의 세제곱, 즉

\text{cost} \sim O(N^3)

에 비례한다. 궤도들을 서로 직교화하는 과정이 이 삼차 스케일링의 주범이다. 수백~수천 원자까지는 현실적이지만, 수만 원자 이상의 계는 여전히 벅차서 선형 스케일링 DFT나 준경험적 기법으로 우회한다.

주요 소프트웨어는 다음과 같다.

VASP — 평면파 기반의 상용 표준. 재료과학계에서 사실상 공용어. 라이선스가 비싸기로도 유명하다.
Quantum ESPRESSO — 오픈소스 평면파 코드. 공짜지만 입력 파일 앞에서 사용자의 정신력이 시험당한다.
Gaussian, ORCA — 분자 화학용 국소 기저함수 코드.
이들 대부분이 원자핵 근처의 급격한 파동함수를 부드럽게 대체하는 유사퍼텐셜(pseudopotential) 기법을 쓴다.

6. 재료·촉매에서의 위상[편집]

DFT는 실험 없이 물질의 성질을 예측하는 “계산 현미경” 역할을 한다. 새로운 배터리 전극 후보 물질의 전압을 미리 계산하고, 촉매 표면에서 반응물이 얼마나 잘 붙는지(흡착 에너지)를 따지고, 미지의 결정 구조가 안정한지를 판정한다. 머티리얼스 프로젝트(Materials Project) 같은 대규모 물성 데이터베이스는 수십만 개 물질에 대한 DFT 계산 결과를 쌓아두고 있으며, 최근에는 이 데이터로 기계학습 퍼텐셜을 훈련시켜 분자동역학의 정확도와 속도를 동시에 잡으려는 시도가 활발하다.

요컨대 DFT는 슈뢰딩거 방정식의 이상과 현실적 계산 가능성 사이에서 인류가 찾아낸 가장 성공적인 타협점이다. 완벽하진 않지만, 없으면 재료 연구가 멈춘다.

7. 관련 문서[편집]

이 발상 자체는 1927년 토머스-페르미 모형까지 거슬러 올라간다. 다만 당시엔 너무 조악해서 화학 결합조차 제대로 설명하지 못했고, 호엔베르크-콘 정리가 나오고서야 엄밀한 이론적 토대를 얻었다. ↩
직관적으로는 말이 안 되는 것 같다. 파동함수는 $3N$ 차원 함수인데 그 모든 정보가 3차원 밀도에 담겨 있다니? 하지만 밀도의 뾰족한 봉우리(cusp) 위치가 원자핵 위치를, 봉우리 높이가 원자번호를 알려주는 식으로, 밀도는 생각보다 많은 것을 품고 있다. ↩
만약 정확한 교환상관 범함수를 아는 사람이 있다면, 그는 다전자 슈뢰딩거 방정식을 사실상 정확히 푼 것이 된다. 즉 이 범함수의 정확한 형태를 찾는 것은 원래 문제만큼 어렵다. 우리는 문제를 푼 게 아니라 어려움을 한 항에 몰아넣었을 뿐이다. ↩