A hierarchy of functionals — Jacob's ladder

DFT functionals are often organised along Perdew's "Jacob's ladder". Each rung adds more ingredients and (typically) more accuracy at higher cost.

Jacob's ladder of DFT functionals — from LDA to double hybrids
RungNameDepends onRepresentative
1LDAdensity ρ onlySVWN — almost never used
2GGAρ + ∇ρBP86, BLYP, PBE
3meta-GGA+ ∇²ρ or τTPSS, M06-L, r²SCAN
4Hybrid+ some HF exchangeB3LYP, PBE0, M06-2X
5Double hybrid+ MP2 correlationB2PLYP, DSD-PBEP86

Higher rungs are typically more accurate, but not always. A well-built meta-GGA (e.g. r²SCAN) can outperform a mediocre hybrid. The latest broad benchmarks (GMTKN55, MGCDB84, …) tend to rank these two functionals at the very top overall:

  • ωB97M-V: meta-GGA + range-separated hybrid + VV10 non-local correlation. Overall #1.
  • ωB97X-V: similar idea but GGA-based. Slightly faster.

Both are built into ORCA: ! wB97M-V def2-TZVP on one line is enough.

Dispersion corrections are essentially required

Pure GGAs and standard hybrids do not describe London dispersion. Without a correction, π–π stacking, hydrogen bonds, vdW complexes, and conformer energetics come out badly. Use one of Grimme's D3 / D4 corrections by default.

KeywordDescription
D4Newest. Recommended in ORCA 6. Fit for more functionals.
D3BJWidely validated. Fallback if D4 is unavailable for your functional.
D3ZEROAlternative D3 damping. Less accurate than D3BJ.
(none)ωB97M-V and ωB97X-V already include VV10, so no extra correction is needed.
# Standard usage with dispersion correction
! B3LYP D4 def2-TZVP RIJCOSX def2/J
! PBE0 D4 def2-TZVP RIJCOSX def2/J
! r2SCAN D4 def2-TZVP def2/J      # meta-GGA doesn't need RIJCOSX

# VV10 already included — no extra correction
! wB97M-V def2-TZVP RIJCOSX def2/J

Recommendations by task

TaskRecommended functionalNotes
Geometry optimisation (general organic)r²SCAN-3c, B97-3cComposite — basis included.
Geometry optimisation (TM complexes)TPSS D4, BP86 D4Pure GGAs are more forgiving of poor initial geometries.
Reaction barriersM06-2X D4, ωB97X-VHigher HF percentage gives better barriers.
Non-covalent interactionsωB97M-V, DSD-PBEP86 D4Double hybrids are the most accurate.
NMR chemical shiftsPBE0 D4, ωB97X D4Use dedicated NMR bases (pcSseg).
UV-Vis absorption (local excitations)B3LYP D4, PBE0 D4Standard hybrids are usually enough.
UV-Vis absorption (charge transfer)CAM-B3LYP, ωB97X-VRange-separation is mandatory.
EPR (g-tensor, A-tensor)PBE0, B3LYP~30 % HF tends to fit well.

Composite methods (the 3c series)

The Grimme group's "3c" methods bundle "functional + basis + D3/D4 + gCP correction" into a single keyword. The basis is fixed, the usage is therefore extremely simple, and the price/accuracy ratio is hard to beat.

KeywordWhat it really isBest for
HF-3cHF + minix (small DZ) + 3 correctionsQuick screening of very large molecules
PBEh-3cPBEh hybrid + def2-mSVPMedium-sized opt + frequencies
B97-3cModified B97 GGA + def2-mTZVPMedium-to-large molecules, fast optimisation
r²SCAN-3cr²SCAN meta-GGA + def2-mTZVPPTop recommendation. Geometries and free energies alike.
ωB97X-3cRange-separated hybrid + valence DZOrganic / excited states
# r²SCAN-3c — a strong first choice for everyday work
! r2SCAN-3c Opt Freq

* xyzfile 0 1 mol.xyz
r²SCAN-3c works well almost everywhere

The manual and many independent benchmarks report that r²SCAN-3c reaches B3LYP-D4/def2-TZVP accuracy in under half the time. Using r²SCAN-3c as the first attempt — and refining with more expensive methods only when needed — is an efficient workflow.

Double hybrids

Double hybrids include both HF exchange and MP2 correlation. They are one rung above the conventional hybrids in accuracy, but they inherit MP2's N⁵ scaling. Best suited to precision single-point energies.

KeywordNotes
B2PLYPThe original double hybrid.
DSD-PBEP86SCS-MP2-coupled. Excellent overall performance.
DSD-BLYPStrong for non-covalent interactions.
wB97X-2Range-separated + double hybrid.
PWPB95High accuracy, cost-efficient.
# DLPNO-accelerated double hybrid — much cheaper
! DLPNO-B2PLYP D4 def2-TZVPP def2-TZVPP/C TightSCF

* xyzfile 0 1 mol.xyz

Recommendations by system

Organic chemistry

  • Geometry optimisation: r²SCAN-3c or B97-3c.
  • Precision single point: ωB97M-V/def2-TZVPP or DLPNO-CCSD(T)/cc-pVTZ.
  • Non-covalent interactions: ωB97M-V is the overall best.

Transition-metal chemistry

  • First try: TPSSh D4 or TPSS D4 (3d metals).
  • Spin-state energetics: meta-GGAs (TPSS, M06-L) are less sensitive to the HF percentage.
  • Spectroscopy (g-tensor): PBE0 D4.

Very large systems (protein active sites, MOFs)

  • Start with r²SCAN-3c or HF-3c.
  • For QM/MM embeddings PBEh-3c tends to be robust.

Excited states (TD-DFT)

  • Local π → π*: B3LYP, PBE0.
  • Charge transfer: CAM-B3LYP, ωB97X-V are essential.
  • Rydberg states: range-separated + diffuse functions.

Functionals to avoid

For historical reasons these functionals are still common, but they are not recommended for new work in the 2020s:

FunctionalReplacement
SVWN / LSDAMove up to a GGA or higher.
B3LYP without dispersionAlways pair with D4 (or D3BJ).
BLYP without dispersionSame — dispersion is mandatory.
The D2 correctionUse D4 or D3BJ instead.
VWN3 variantsLittle effect today. B3LYP/G is for Gaussian compatibility only.
B3LYP is still the "reference"

B3LYP D4 is not the most accurate functional available today, but the overwhelming majority of the literature is written with B3LYP, so it still has value as a comparison baseline. For new work, reporting modern functionals (ωB97M-V, r²SCAN-3c, DSD-PBEP86) alongside B3LYP gives both literature consistency and higher accuracy at once.