밀도범함수 이론(Density Functional Theory, DFT)은 물질, 분자 내부에 전자가 들어있는 모양과 그 에너지를 양자역학으로 계산하기 위한 이론 및 계산 방법이다. 다전자(多電子) 계의 바닥상태 에너지와 전자 밀도를 전자 밀도 함수 n(r)을 기본 변수로 하여 기술한다. 전자 파동함수 대신 3차원 공간상의 전자 밀도만을 다루므로 계산 효율이 크게 향상되며, 컴퓨터를 사용하는 과학 계산 중 가장 널리 쓰이는 양자역학 계산 분야 중 하나이다.
역사
밀도범함수 이론의 기초는 1927년 L. H. 토마스(Thomas)와 엔리코 페르미(Enrico Fermi)가 독립적으로 제안한 통계적 모델(Thomas-Fermi model)에서 비롯되었다. 이 모델은 전자가 위상 공간에 균일하게 분포한다고 가정하여 운동 에너지를 전자 밀도의 함수로 표현하였으나, 분자 결합을 설명하지 못하는 등 정확도에 한계가 있었다. 1964년 피에르 호헨베르크(Pierre Hohenberg)와 월터 콘(Walter Kohn)은 호헨베르크-콘 정리(Hohenberg-Kohn theorem)를 발표하여 밀도범함수 이론에 엄밀한 이론적 기초를 제공하였다. 1965년 월터 콘과 루 전 샴(Lu Jeu Sham)은 콘-샴 방정식(Kohn-Sham equations)을 도입하여 실제 상호작용하는 전자계를 상호작용이 없는 가상의 전자계로 근사하는 실용적인 계산 체계를 확립하였다. 콘은 이 공로로 1998년 노벨 화학상을 수상하였다.
핵심 원리
호헨베르크-콘의 제1정리에 따르면, 바닥상태의 전자 밀도 n(r)이 주어지면 계의 외부 퍼텐셜(예: 핵 전하 분포)이 유일하게 결정되며, 따라서 모든 물리적 성질이 전자 밀도의 범함수(functional)로 표현될 수 있다. 제2정리는 총 에너지 범함수 E[n]이 바닥상태 밀도에서 최솟값을 가짐을 증명한다.
콘-샴 접근법에서는 상호작용하는 실제 전자계를 상호작용이 없는 가상의 전자계(콘-샴 계)로 대체한다. 이 가상 계의 단일-입자 궤도함수(콘-샴 오비탈)는 다음과 같은 콘-샴 방정식을 만족한다:
$$ \left[ -\frac{1}{2} abla^2 + V_{\text{eff}}(\mathbf{r}) \right] \phi_i(\mathbf{r}) = \varepsilon_i \phi_i(\mathbf{r}) $$
여기서 유효 퍼텐셜 $V_{\text{eff}}$는 외부 퍼텐셜, 하트리(Hartree) 퍼텐셜(전자 간 정전기적 상호작용), 그리고 교환-상관(exchange-correlation) 퍼텐셜의 합으로 구성된다. 교환-상관 퍼텐셜은 다체 효과를 집약적으로 표현하는 항으로, 그 정확한 형태는 알려져 있지 않으며 근사가 필요하다.
교환-상관 범함수의 근사
교환-상관 에너지 범함수 $E_{\text{XC}}[n]$의 근사 방법은 밀도범함수 이론의 정확도를 결정짓는 핵심 요소이다. 주요 근사법은 다음과 같다.
- 국소 밀도 근사(Local Density Approximation, LDA): 각 지점에서의 전자 밀도만을 사용하여 균일한 전자 기체의 교환-상관 에너지로 근사한다. 금속 및 조밀한 전자계에 적합하다.
- 일반화된 기울기 근사(Generalized Gradient Approximation, GGA): 전자 밀도와 그 기울기(∇n)를 함께 고려하여 비균질성을 보정한다. PBE(Perdew-Burke-Ernzerhof) 범함수가 대표적이다.
- 메타-GGA(Meta-GGA): 밀도의 2차 도함수(라플라시안) 또는 운동 에너지 밀도를 추가로 포함한다.
- 하이브리드 범함수(Hybrid functional): 하트리-폭(Hartree-Fock) 이론의 정확한 교환 에너지 일부를 혼합한다. B3LYP, PBE0 등이 대표적이다.
수학적 기초
호헨베르크-콘 정리는 다음과 같이 요약된다.
- 해밀토니언이 주어지면 파동함수를 구할 수 있고, 파동함수에서 밀도함수를 구할 수 있다. 역으로 (이론상) 밀도함수에서 파동함수를 구할 수도 있다. 즉, 바닥상태 밀도 함수와 외부 퍼텐셜은 일대일 대응 관계에 있다.
- 변분 원리를 이용하여 정확한 바닥상태 에너지를 계산할 수 있다.
- 외부 퍼텐셜을 제외한 나머지 항(운동 에너지 + 전자-전자 상호작용)은 보편적(universal)이다.
콘-샴 방정식은 자체-일관성 장(Self-Consistent Field, SCF) 방법으로 풀이한다. 초기 전자 밀도를 가정하고 유효 퍼텐셜을 계산한 후 콘-샴 방정식을 풀어 새로운 파동함수와 밀도를 얻는다. 이 과정을 수렴할 때까지 반복한다.
장점과 한계
밀도범함수 이론은 계산 효율성과 실용적 정확도의 균형으로 인해 분자 구조 최적화, 총에너지 계산, 결합 에너지, 전자 밀도, 상태 밀도(DOS), 밴드 구조, 진동수, 반응 경로 등 다양한 물리·화학적 성질을 계산하는 데 널리 사용된다. 응용 분야로는 재료과학(밴드 구조, 결함, 표면, 촉매), 화학(분자 구조·반응 경로·스펙트럼 예측), 나노과학, 배터리 전극 재료, 태양전지 재료 설계 등이 있다.
그러나 다음과 같은 한계도 존재한다. LDA와 GGA는 반도체 및 절연체의 밴드갭을 과소평가하는 경향이 있다. 강한 상관관계를 가진 계(Mott 절연체 등)에서는 표준 DFT가 실패할 수 있어 DFT+U, DMFT 등의 보완 기법이 필요하다. 또한 분자간 상호작용, 특히 반데르발스 힘(분산력)의 처리가 불완전하며, 전자 들뜬 상태(광학 스펙트럼 등) 계산은 기본 DFT로는 정확하지 않아 시간의존 밀도범함수 이론(TDDFT)이나 GW, BSE 등의 보완이 필요하다.
소프트웨어
밀도범함수 이론을 구현한 주요 소프트웨어로는 VASP, Quantum ESPRESSO, ABINIT, GPAW, Gaussian, NWChem, CP2K 등이 있다. 대부분의 소프트웨어는 유닉스 또는 리눅스 환경을 기반으로 하며, 포트란 소스 코드 형태로 제공되거나 상용 소프트웨어로 배포된다.
