结合能/吸附能计算步骤及初猜选择问题

gyp

已经阅读谈谈分子间结合能的构成以及分解分析思想 - 思想家公社的门口：量子化学·分子模拟·二次元，了解结合能与相互作用能区别。
现在有些疑问：

计算步骤1：
手动摆A分子初猜，进行结构优化及单点计算得到优化后的结构A1及能量E(A1)
手动摆B分子初猜，进行结构优化及单点计算得到优化后的结构B1及能量E(B1)
手动把优化好的结构A1和B1摆在一起进行结构优化及单点计算，得到优化后的结构AB1及能量E(AB1)
计算Ebind1=E(AB1)-E(A1)-E(B1)

计算步骤2：
手动把未经优化的A、B分子摆在一起，进行结构优化及单点计算得到优化后的结构AB2及能量E(AB2)
把优化好的AB2结构中A分子坐标删除，保留B分子坐标，进行结构优化及单点计算得到优化后的结构B2及能量E(B2)
同样方法得到结构A2及能量E(A2)
计算Ebind2=E(AB2)-E(A2)-E(B2)

问题1：计算步骤2得到的结合能Ebind2是否有意义？我理解的结合能应该是按计算步骤1来的，是否正确？
问题2：相互作用能的计算，是在步骤2中，删除坐标后不做结构优化直接做单点计算还是应该在步骤1结束后再删除坐标进行单点计算（即A、B分子需不需要事先优化）
主要是觉得不同的初猜会优化到势能面不同的极小点导致能量有差异

Stardust0831

其实不需要这么复杂，各自独立建模和结构优化即可。你提到的“优化到势能面不同的极小点”是个很敏锐的观察，一般是需要做构象搜索来找出合理的结构。具体可以看社长的帖子：《使用molclus程序做团簇构型搜索和分子构象搜索》

一切从简的情况下，确实可以先优化团簇再删坐标优化单体，这种情况下涉及三次结构优化，这三次结构优化前都可以用半经验DFT（比如GFN2-xTB）先预优化一下再上正式的DFT,会省不少耗时，社长的molclus程序里也实现了现成的工作流。
记得加gCP或者选用足够完备的基组来处理计算时的BSSE问题。