함수의 위계 — 야곱의 사다리
DFT 함수는 “야곱의 사다리(Jacob's ladder)”라는 비유로 분류되곤 한다. 사다리를 오를수록 정확도와 비용이 함께 올라간다.
| 층 | 이름 | 의존성 | 대표 예 |
|---|---|---|---|
| 1 | LDA | 밀도 ρ만 | SVWN — 거의 사용 안 함 |
| 2 | GGA | ρ + ∇ρ | BP86, BLYP, PBE |
| 3 | 메타-GGA | + ∇²ρ 또는 τ | TPSS, M06-L, r²SCAN |
| 4 | 하이브리드 | + HF 교환 일부 | B3LYP, PBE0, M06-2X |
| 5 | 이중 하이브리드 | + MP2 상관 | B2PLYP, DSD-PBEP86 |
높은 층일수록 대체로 정확하지만 항상 그런 건 아니다. 잘 만든 메타-GGA(r²SCAN)가 평범한 하이브리드보다 더 정확할 때도 있다. 최근의 광범위한 벤치마크(GMTKN55, MGCDB84 등)에서는 아래 두 함수가 종합 1~2위를 차지한다.
- ωB97M-V: 메타-GGA + 장범위 분리 하이브리드 + VV10 비국소 상관. 종합 1위.
- ωB97X-V: 위와 유사하지만 GGA 기반. 약간 더 빠름.
두 함수 모두 ORCA에 내장돼 있어, ! wB97M-V def2-TZVP 한 줄로 끝난다.
분산력 보정은 거의 필수
순수 GGA와 표준 하이브리드는 분산력(London 분산력)을 잡지 못한다. 그래서 π-π 적층, 수소결합, vdW 복합체, 콘포머 안정성 같은 데서 큰 오차가 난다. Grimme의 D3·D4 보정을 같이 켜는 게 좋다.
| 키워드 | 설명 |
|---|---|
D4 | 최신. ORCA 6 권장. 더 많은 함수에 fit. |
D3BJ | 널리 검증됨. D4가 없을 때의 대안. |
D3ZERO | D3의 또 다른 damping. D3BJ보다 덜 정확. |
| (없음) | ωB97M-V, ωB97X-V는 VV10이 내장되어 있어 별도 보정 불필요. |
# 분산력 보정이 포함된 표준 사용
! B3LYP D4 def2-TZVP RIJCOSX def2/J
! PBE0 D4 def2-TZVP RIJCOSX def2/J
! r2SCAN D4 def2-TZVP def2/J # 메타-GGA는 RIJCOSX 불필요
# VV10 내장 — 별도 보정 없음
! wB97M-V def2-TZVP RIJCOSX def2/J
작업별 추천
| 작업 | 추천 함수 | 비고 |
|---|---|---|
| 구조 최적화 (유기, 일반) | r²SCAN-3c, B97-3c | 합성 방법 — 기저집합 자동. |
| 구조 최적화 (전이금속 착물) | TPSS D4, BP86 D4 | 순수 GGA가 발산을 잘 견딤. |
| 반응 장벽 | M06-2X D4, ωB97X-V | 높은 HF 비율이 장벽을 잘 묘사. |
| 비공유 결합 (수소결합·π-π) | ωB97M-V, DSD-PBEP86 D4 | 이중 하이브리드가 가장 정확. |
| NMR 화학 이동 | PBE0 D4, ωB97X D4 | 전용 기저집합 권장(pcSseg). |
| UV-Vis 흡수 (국소 들뜸) | B3LYP D4, PBE0 D4 | 표준 하이브리드로 충분. |
| UV-Vis 흡수 (전하 이동) | CAM-B3LYP, ωB97X-V | 장범위 분리가 필수. |
| EPR (g-tensor, A-tensor) | PBE0, B3LYP | 30% HF 정도가 잘 맞음. |
합성 방법 (3c 시리즈)
Grimme 그룹의 “3c” 시리즈는 “함수 + 기저집합 + D3·D4 + gCP 보정”을 하나로 묶은 합성 방법이다. 기저집합이 자동으로 정해져 쓰기가 아주 간단하고, 비용 대비 정확도가 뛰어나다.
| 키워드 | 실제 구성 | 적합한 작업 |
|---|---|---|
HF-3c | HF + minix(작은 DZ) + 3개 보정 | 거대 분자의 빠른 스크리닝 |
PBEh-3c | PBEh 하이브리드 + def2-mSVP | 중대형 분자 최적화·진동수 |
B97-3c | 수정 B97 GGA + def2-mTZVP | 중대형 분자, 빠른 구조 최적화 |
r²SCAN-3c | r²SCAN 메타-GGA + def2-mTZVPP | 가성비 최강 추천. 자유에너지·구조 모두. |
ωB97X-3c | 장범위 분리 하이브리드 + valence DZ | 유기·들뜸 상태에 강점 |
# r²SCAN-3c — 일상 작업의 첫 시도로 강력 추천
! r2SCAN-3c Opt Freq
* xyzfile 0 1 mol.xyz
매뉴얼과 여러 외부 벤치마크에서 r²SCAN-3c는 B3LYP D4/def2-TZVP 수준의 정확도를 절반도 안 되는 시간에 낸다고 보고됐다. 일단 r²SCAN-3c로 시작하고, 필요할 때만 더 비싼 방법으로 정밀화하는 워크플로우가 효율적이다.
이중 하이브리드
이중 하이브리드(double hybrid)는 HF 교환과 MP2 상관을 둘 다 포함해서 하이브리드보다 한 단계 더 정확하지만, MP2 비용이 붙어 N⁵ 스케일링이 된다. 단일점 정밀 에너지에 잘 맞는다.
| 키워드 | 특징 |
|---|---|
B2PLYP | 고전 이중 하이브리드. |
DSD-PBEP86 | SCS-MP2 계열. 종합 성능 우수. |
DSD-BLYP | 비공유 상호작용에 강점. |
wB97X-2 | 장범위 분리 + 이중 하이브리드. |
PWPB95 | 높은 정확도, 비용 효율적. |
# DLPNO를 활용한 이중 하이브리드 — 비용 대폭 절감
! DLPNO-B2PLYP D4 def2-TZVPP def2-TZVPP/C TightSCF
* xyzfile 0 1 mol.xyz
시스템별 추천
유기 화학
- 구조 최적화:
r²SCAN-3c또는B97-3c. - 정밀 단일점:
ωB97M-V/def2-TZVPP 또는DLPNO-CCSD(T)/cc-pVTZ. - 비공유 결합:
ωB97M-V가 종합 최선.
전이금속 화학
- 1차 후보:
TPSSh D4또는TPSS D4(3d 금속). - 스핀 상태 에너지: 메타-GGA(
TPSS,M06-L)가 HF 비율의 영향을 적게 받음. - 분광학(g-tensor):
PBE0 D4.
거대 분자 (단백질 활성자리, MOF)
r²SCAN-3c또는HF-3c로 시작.- QM/MM 결합 시
PBEh-3c가 안정적.
들뜬 상태 (TD-DFT)
- 국소 π → π*:
B3LYP,PBE0. - 전하 이동:
CAM-B3LYP,ωB97X-V가 필수. - 리이드버그 상태: 장범위 분리 + 확산 함수.
선택을 피해야 할 함수
아래 함수들은 역사적 이유로 자주 보이지만, 2020년대 양자화학에서는 권하지 않는다.
| 함수 | 대체 |
|---|---|
SVWN / LSDA | GGA 이상 사용 권장. |
B3LYP (분산 보정 없이) | 반드시 D4 또는 D3BJ와 함께. |
BLYP (분산 보정 없이) | 같은 이유로 분산 보정 필수. |
D2 분산 보정 | D4 또는 D3BJ로 대체. |
VWN3로 정의된 변종 | 요즘은 별 효과 없음. B3LYP/G는 Gaussian 호환용으로만. |
B3LYP D4가 지금 권장되는 최고의 함수는 아니지만, 문헌의 절대다수가 여전히 B3LYP로 쓰여 있어서 “비교 가능한 결과”를 얻는 데는 의미가 있다. 새 연구라면 ωB97M-V, r²SCAN-3c, DSD-PBEP86 같은 최신 함수와 B3LYP 결과를 같이 보고하면, 과거 문헌과의 일관성을 지키면서도 정확도를 높일 수 있다.