교환-상관 범함수

양자화학 계산물리 마지막 수정: 2026-07-08 05:17:22

교환-상관 범함수 Exchange-correlation functional
분야	밀도범함수이론 × 양자화학
역할	DFT의 모든 미지 물리를 떠안는 항
주요 등급	LDA, GGA, meta-GGA, 하이브리드
비유	야곱의 사다리(Jacob's Ladder)

1. 개요[편집]

DFT는 정확한 이론이다. 단 하나, 아무도 그 정확한 범함수를 모른다는 점만 빼면.

교환-상관 범함수(exchange-correlation functional, XC functional)는 밀도범함수이론(DFT)에서 전자들 사이의 복잡한 양자역학적 상호작용 중 고전적으로 다룰 수 없는 부분 — 즉 교환(exchange) 에너지와 상관(correlation) 에너지 — 을 전자밀도 $n(\mathbf{r})$ 의 함수로 표현한 항이다. 코-샴(Kohn-Sham) DFT에서 전체 에너지는 다음처럼 쓰인다.

E[n] = T_s[n] + E_{\text{ext}}[n] + E_H[n] + E_{xc}[n]

앞의 세 항(무상호작용 운동에너지, 외부 퍼텐셜, 고전 하트리 반발)은 정확히 계산할 수 있다. 문제는 마지막 $E_{xc}[n]$ . 호엔베르크-코-샴 정리는 이 범함수가 존재한다고 보장하지만, 그 구체적 형태는 알려주지 않는다.¹ 그래서 DFT의 실전은 “얼마나 좋은 $E_{xc}$ 근사를 고르느냐”의 예술이 된다. 하트리-폭 방법이 교환을 정확히 다루되 상관을 통째로 버리는 것과 정반대로, DFT는 교환과 상관을 근사적으로나마 한꺼번에 담으려 한다.

2. 왜 근사가 필요한가[편집]

정확한 다체 파동함수를 풀면 전자 상관을 다 담을 수 있지만, 그건 지수적으로 비싼 양자 몬테카를로나 CCSD(T) 급 방법의 몫이다. DFT의 매력은 3차원 밀도 하나만 다루면 되기 때문에 수백~수천 원자를 감당할 수 있다는 것. 이 저렴함의 대가가 바로 “미지의 $E_{xc}$ 를 근사로 때운다”는 파우스트적 거래다.

교환 에너지는 파울리 배타 원리에서 오는, 같은 스핀 전자끼리 서로 피하는 효과다. 상관 에너지는 그 너머, 모든 전자가 서로의 순간 위치를 눈치보며 움직이는 동역학적 회피 효과를 말한다. 정확한 XC를 알면 DFT는 원리적으로 정확하다. 그래서 지난 반세기 계산과학의 상당 부분이 “더 좋은 범함수 만들기” 경쟁이었다.

3. 야곱의 사다리[편집]

존 퍼듀(John Perdew)는 XC 범함수의 발전을 야곱의 사다리(Jacob’s Ladder)에 비유했다. 땅(하트리 근사, 화학적 정확도 0)에서 천국(정확한 XC, 화학적 정확도)으로 올라가는 사다리이며, 각 단은 더 많은 정보를 재료로 넣는다.²

3.1. 1단: LDA (국소밀도근사)[편집]

국소밀도근사(Local Density Approximation, LDA)는 각 점의 XC 에너지 밀도를 그 점의 밀도값과 똑같은 밀도를 가진 균일 전자가스(uniform electron gas)의 값으로 근사한다.

E_{xc}^{\text{LDA}}[n] = \int n(\mathbf{r})\, \varepsilon_{xc}\big(n(\mathbf{r})\big)\, d\mathbf{r}

놀랍게도 이 조잡한 근사가 고체의 구조·격자상수는 꽤 잘 맞춘다. 대신 결합에너지를 과대평가(overbinding)하는 고질병이 있다.

3.2. 2단: GGA (일반화 기울기 근사)[편집]

일반화 기울기 근사(Generalized Gradient Approximation, GGA)는 밀도값뿐 아니라 그 기울기 $\nabla n$ 까지 재료로 쓴다. 밀도가 얼마나 빠르게 변하는지를 반영해 분자 결합에너지가 크게 개선됐다. PBE(Perdew-Burke-Ernzerhof)와 BLYP가 대표 주자로, VASP를 비롯한 고체·표면 계산의 사실상 기본값이다.

3.3. 3단: meta-GGA[편집]

여기에 운동에너지 밀도 $\tau$ 나 밀도의 라플라시안까지 넣으면 meta-GGA가 된다. SCAN 범함수가 유명하며, 다양한 결합 종류를 한 범함수로 골고루 맞추려는 야심작이다.

3.4. 4단: 하이브리드[편집]

하이브리드 범함수(hybrid functional)는 DFT 교환에 하트리-폭 방법의 정확한 교환(exact exchange) 일부를 섞는다. 예컨대 그 유명한 B3LYP는 다음처럼 조합한다.

E_{xc}^{\text{B3LYP}} = (1-a)E_x^{\text{LDA}} + a E_x^{\text{HF}} + b\,\Delta E_x^{\text{GGA}} + \cdots

정확한 교환을 섞으면 자기상호작용 오차(self-interaction error)가 줄어 밴드갭·반응 장벽 예측이 좋아진다. B3LYP는 한때 계산화학 논문의 절대다수를 점령했을 정도로 국룰이었다.³ 고체용으로는 HSE06(단거리 스크리닝 하이브리드)이 밴드갭 계산의 표준이다. 대신 HF 교환이 들어가면 평면파 기저 계산에서 비용이 훅 뛴다.

3.5. 5단: 이중 하이브리드 이상[편집]

가상 궤도함수(virtual orbital)의 정보까지 끌어와 MP2 상관을 섞는 이중 하이브리드(double-hybrid, 예: B2PLYP)가 사다리 꼭대기 근처다. 정확하지만 비싸다. 사다리를 오를수록 정확도와 비용이 함께 오르는 것은 이 바닥의 만고불변 법칙.

4. 고르는 법과 함정[편집]

“어떤 범함수를 쓰죠?”는 계산화학 초심자의 영원한 질문이자, 사수가 가장 답하기 싫어하는 질문이다.⁴ 정답은 “계와 물성에 따라 다르다”.

반데르발스/분산력: 순수 GGA는 분산력을 아예 못 본다. DFT-D3 같은 분산 보정을 반드시 얹거나 vdW-DF 계열을 써야 한다.
밴드갭: LDA/GGA는 반도체 밴드갭을 30~100% 과소평가한다. 이 “밴드갭 문제”는 HSE06 같은 하이브리드나 GW 보정으로 완화한다. 밴드 구조 계산 참고.
강상관계(전이금속 산화물 d/f 전자): 표준 XC로는 금속으로 잘못 예측하기 일쑤. DFT+U 보정이 흔히 동원된다.
반응 장벽: 순수 DFT는 장벽을 낮게 본다. 하이브리드가 유리.

명심할 점 — 사다리를 높이 올라간다고 항상 더 정확한 건 아니다. 특정 물성에 맞춰 파라미터를 피팅한 범함수는 다른 물성에서 배신할 수 있다. 그래서 검증 및 확인의 정신으로, 실험이나 고급 이론(CCSD(T))에 대해 벤치마크한 뒤 쓰는 게 정도다.

5. 관련 문서[편집]

이게 DFT의 가장 아름답고도 잔인한 지점이다. 정확한 답이 존재한다는 건 증명됐는데, 그게 뭔지는 아무도 모른다. 마치 로또 1등 번호가 존재한다는 걸 아는 것과 그걸 맞히는 것의 차이랄까. ↩
창세기의 야곱이 꿈에서 본 하늘까지 닿는 사다리에서 따왔다. 퍼듀 본인이 신학적 비유를 좋아한 건지, 아니면 정확한 범함수가 정말 “천국”만큼 도달 불가능해 보였던 건지는 알 수 없다. ↩
한때 “B3LYP/6-31G(d)로 돌렸습니다”는 계산화학 논문의 정체성 그 자체였다. 지금은 분산 보정 없는 B3LYP를 그냥 쓰면 리뷰어에게 지적당하는 시대가 됐다. 국룰도 세월 앞에 장사 없다. ↩
왜냐하면 정직한 답은 “일단 몇 개 돌려보고 실험이랑 제일 잘 맞는 걸 골라라”인데, 이걸 대놓고 말하면 커브 피팅(curve fitting)처럼 들리기 때문이다. 실제로 종종 그렇기도 하다. ↩