量子化学模拟如果从”用什么方法”开始讨论，通常会在前十分钟陷入术语泥潭。但换一个视角——从”我要解决什么问题”出发——路线反而清晰了。
封闭壳层有机分子、化学键不涉及明显多参考效应、不需要精确的热化学精度：HF加MP2是起步方案。HF把交换能算准了而相关能完蛋，MP2把相关能补了一阶。这个组合对几何结构预测的误差通常在0.02埃以内，频率的偏差一般不超过5%。超过两百个原子的体系就别想MP2了，老老实实回到DFT。
开壳层体系、过渡金属、涉及到自由基反应的场景：DFT是现阶段性价比最高的选择。B3LYP在有机化学里依然是通用工作马，但如果体系有较多电荷分离态，M06-2X和ωB97X-D这些含长程校正的泛函会明显优于传统杂化泛函。本质上来讲，HF交换成分的比例是决定性的——B3LYP是20%，M06-2X是54%。一个直观的经验判断：如果你的体系沿某个键长方向有电荷转移特征（比如D-A型推拉体系），交换成分偏高的泛函大概率给出更合理的势能面。
但这里有一个需要仔细权衡的细节：高HF比例的泛函在热化学上的表现并不总是更好。ωB97X-D的原子化能误差在某些过渡金属配合物上反而比B3LYP大——这是因为HF交换项在描述静态相关方面天然短板，HF比例越高，对含有显著多组态特征的体系的描述越偏差。做计算之前先对自己的体系做一个基本的诊断：算一下HOMO-LUMO gap，如果gap小于2eV，说明体系可能有一定的多参考特征，此时DFT泛函选择的优先级不再是HF比例高低，而是看这个泛函在处理近简并态时的经验表现。
高精度热化学要求——反应能垒或结合能的误差要求在1kcal/mol以内：这就不是DFT能稳定达成的领域了。耦合簇方法（CCSD(T)）是gold standard，但计算量随体系大小的O(N⁷)标度让人头疼。做过一个含18个重原子的反应路径，DFT跑了3小时，CCSD(T)单点跑了将近一整天。这个时间差就是量子化学里最残酷的现实：精度和计算量之间不存在平滑过渡，而是一条断崖式的分界线。现在的折中方案是用DLPNO-CCSD(T)——局部自然轨道近似把标度压到接近线性，对中等体系来说是在可接受的时间内拿到接近CCSD(T)级别的结果。
最后值得提醒的是，任何方法都有失效区。近简并、强关联、涉及过渡金属的多组态问题——这时候UHF可能的自旋污染、DFT的离域误差、耦合簇的单参考假设都面临不同程度的失效。不要在一个失效的方法上反复调参，换一个更适合的方法才是正解。曾经见过一个博士生在一个有明显双自由基特征的体系上用B3LYP反复调基组，调了三个月，最后换成CASPT2一周就出了合理结果。选对方法的收益远大于在错方法上调参。
还有一点关于基组叠加误差（BSSE）的说明：很多人以为BSSE只在弱相互作用体系（如范德华复合物）里重要，但实际上在计算反应能垒时，过渡态和反应物之间基组完备性的差异也会引入系统性能量偏差。用counterpoise校正对反应路径上的每个驻点分别做BSSE评估，在论文的方法部分交代清楚校正前后的对比数据——这个透明度的提升比校正本身带来的数值变化更有价值。
https://www.keyanxueshu.com/category/dft/
