함수의 위계 — 야곱의 사다리

DFT 함수는 “야곱의 사다리(Jacob's ladder)”라는 비유로 분류되곤 한다. 사다리를 오를수록 정확도와 비용이 함께 올라간다.

DFT 함수의 야곱의 사다리 — LDA부터 이중 하이브리드까지, 각 단이 재료 하나씩 추가
이름의존성대표 예
1LDA밀도 ρ만SVWN — 거의 사용 안 함
2GGAρ + ∇ρ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가 없을 때의 대안.
D3ZEROD3의 또 다른 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, B3LYP30% HF 정도가 잘 맞음.

합성 방법 (3c 시리즈)

Grimme 그룹의 “3c” 시리즈는 “함수 + 기저집합 + D3·D4 + gCP 보정”을 하나로 묶은 합성 방법이다. 기저집합이 자동으로 정해져 쓰기가 아주 간단하고, 비용 대비 정확도가 뛰어나다.

키워드실제 구성적합한 작업
HF-3cHF + minix(작은 DZ) + 3개 보정거대 분자의 빠른 스크리닝
PBEh-3cPBEh 하이브리드 + def2-mSVP중대형 분자 최적화·진동수
B97-3c수정 B97 GGA + def2-mTZVP중대형 분자, 빠른 구조 최적화
r²SCAN-3cr²SCAN 메타-GGA + def2-mTZVPP가성비 최강 추천. 자유에너지·구조 모두.
ωB97X-3c장범위 분리 하이브리드 + valence DZ유기·들뜸 상태에 강점
# r²SCAN-3c — 일상 작업의 첫 시도로 강력 추천
! r2SCAN-3c Opt Freq

* xyzfile 0 1 mol.xyz
r²SCAN-3c는 거의 어디서나 잘 통한다

매뉴얼과 여러 외부 벤치마크에서 r²SCAN-3c는 B3LYP D4/def2-TZVP 수준의 정확도를 절반도 안 되는 시간에 낸다고 보고됐다. 일단 r²SCAN-3c로 시작하고, 필요할 때만 더 비싼 방법으로 정밀화하는 워크플로우가 효율적이다.

이중 하이브리드

이중 하이브리드(double hybrid)는 HF 교환과 MP2 상관을 둘 다 포함해서 하이브리드보다 한 단계 더 정확하지만, MP2 비용이 붙어 N⁵ 스케일링이 된다. 단일점 정밀 에너지에 잘 맞는다.

키워드특징
B2PLYP고전 이중 하이브리드.
DSD-PBEP86SCS-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 / LSDAGGA 이상 사용 권장.
B3LYP (분산 보정 없이)반드시 D4 또는 D3BJ와 함께.
BLYP (분산 보정 없이)같은 이유로 분산 보정 필수.
D2 분산 보정D4 또는 D3BJ로 대체.
VWN3로 정의된 변종요즘은 별 효과 없음. B3LYP/G는 Gaussian 호환용으로만.
B3LYP는 여전히 “표준 참조”

B3LYP D4가 지금 권장되는 최고의 함수는 아니지만, 문헌의 절대다수가 여전히 B3LYP로 쓰여 있어서 “비교 가능한 결과”를 얻는 데는 의미가 있다. 새 연구라면 ωB97M-V, r²SCAN-3c, DSD-PBEP86 같은 최신 함수와 B3LYP 결과를 같이 보고하면, 과거 문헌과의 일관성을 지키면서도 정확도를 높일 수 있다.