结合案例图分析，我把这个项目的完整过程、参数设置、结果解读和踩坑经验整理如下。
客户要预测两个中性小分子之间的非共价相互作用，用于指导超分子自组装设计。他关心的是：最优间距是多少、最优相对取向如何、相互作用能有多强。
我采用的是势能面扫描方案，软件用 Gaussian 16 做 DFT-D3 计算（B97-D/def2-TZVP），并用 SAPT(DFT) 做能量分解作为对照。之所以选 DFT-D3，是因为中性分子间主要是色散作用，普通 GGA 会严重低估；SAPT 则能把总能拆成静电、诱导、色散和交换排斥。
参数设置：固定一个分子，另一个分子沿质心距 R 和取向角 θ 扫描。R 范围 0.25–0.90 nm，步长 0.01–0.02 nm；θ 范围 0–180°。每个点做单点能计算，最后用二维插值得到势能面。
第一张案例图是相互作用能随分子间距 R 的变化曲线。横轴是质心距 R（nm），纵轴是相互作用能（kJ·mol⁻¹）。当 R<0.34 nm 时，电子云重叠产生强排斥（>4 kJ·mol⁻¹）；在 R≈0.34 nm 处达到能量最低点 E_min=-4.0 kJ·mol⁻¹；R>0.34 nm 后能量逐渐趋近于 0，体现长程吸引尾巴。这个势阱深度属于中等偏弱相互作用，和典型氢键（10–30 kJ·mol⁻¹）比要弱，但足以驱动低温下的自组装。
第二张案例图是二维势能等高线图，横轴是间距 R，纵轴是取向角 θ。能量最低区域（深蓝色）出现在 R≈0.34–0.40 nm、θ≈20–60° 的范围内，说明两个分子之间有一个较宽的最佳取向窗口，而不是严格的单点最优。这个”能量谷较宽”的特征对自组装是有利的，意味着分子在形成聚集体时对取向偏差有一定容忍度。
第一，色散修正必须加。我用 B3LYP/6-31G(d) 先扫了一遍，得到的最优距离是 0.42 nm、E_min=-1.2 kJ·mol⁻¹，明显偏浅；换成 B97-D3 后变成 0.34 nm、-4.0 kJ·mol⁻¹，和 SAPT 分解结果一致。
第二，BSSE 校正不能忘。在没有 counterpoise 校正时，E_min 被高估约 0.8 kJ·mol⁻¹；校正后势阱深度才可靠。
第三，扫描步长要足够细。R 方向如果步长大于 0.05 nm，会漏掉势阱最优点；我用 0.01 nm 在势阱附近加密。
以后做分子间相互作用预测，我会采用”DFT-D3 快速扫描 + SAPT 能量分解”的组合：前者给出势能面，后者解释相互作用本质。交付时同时给 1D 势能曲线和 2D 角度-距离图，让客户既能看到最优构型，又能设计柔性自组装体系。
计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验 需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。
