高斯计算电离能

项目过程：苯分子电离能的多方法收敛测试

客户需要苯分子的第一电离能，用于校准他自己开发的有机半导体电离能预测模型。他的要求很明确：结果要和实验值 9.243 eV 接近，方法要有系统性，能说明不同理论水平的偏差。

我用 Gaussian 做了从小到大多层次的计算：HF/Koopmans、MP2、B3LYP、CAM-B3LYP，一直到 CCSD(T)/CBS。基组统一用 cc-pVTZ，几何优化在每个方法下单独做。电离能分两种定义：垂直电离能（VIP，保持中性分子几何）和绝热电离能（AIP，阳离子也做几何优化）。

计算流程是：先优化中性苯分子 S₀ 态，再优化阳离子 D₀ 态，最后分别用两种几何做单点能计算。

案例分析：CCSD(T)/CBS 几乎完美复现实验

第一张案例图把中性 S₀ 和阳离子 D₀ 的势能面画在了一起。两条抛物线的最低点能量差就是绝热电离能，图上标注 AIP=9.24 eV；从中性最低点垂直向上到阳离子势能面的能量差是垂直电离能，VIP=9.26 eV。两者只差 0.02 eV，说明苯的几何弛豫很小，这是芳香环高度刚性的体现。

第二张案例图对比了不同方法的电离能。HF（Koopmans）给出约 9.42 eV，明显偏高；MP2 约 9.15 eV，偏低；B3LYP 约 9.20 eV；CAM-B3LYP 约 9.27 eV；CCSD(T)/CBS 约 9.25 eV，和实验值 9.243 eV 几乎重合。灰色柱是各种近似方法，红色柱是实验值，一眼就能看出 CCSD(T) 的可靠性。

总结教训：方法选择比基组增大更重要

第一，HF 用 Koopmans 定理估算电离能会系统偏高约 0.2 eV，因为 Koopmans 忽略了阳离子的弛豫和相关能。仅适合做量级估计，不能做定量。

第二，CCSD(T)/CBS 虽然贵，但对芳香烃电离能可以达到实验精度。如果客户预算有限，CAM-B3LYP/cc-pVTZ 是一个误差约 0.05 eV 的性价比方案。

第三，垂直和绝热电离能的差别不能忽略。对于柔性分子，AIP 和 VIP 可能差 0.3 eV 以上，本例因为苯环刚性所以差别很小。

项目经验：把多方法对比做成校准数据集

客户最后把这套计算结果作为训练集标签，用来评估他的机器学习模型。这个思路可以推广：对一类分子（芳香烃、噻吩衍生物等）批量算 CCSD(T)/CBS，再用 DFT 方法做快速预测，用误差分布来决定哪个泛函最适合这类体系。

这篇文章把我在这个项目里用到的计算流程、参数取舍和图上读数做了完整复盘。如果你也在做类似方向，可以参考我的思路：高斯计算电离能。