在进行两个反应物的吸附能计算时，采用同样的命令，构建相近的模型但结果差异巨大

sommnus

      求助各位，我在进行计算环氧丙烷（PO）和环氧丁烷（BO）吸附在改性分子筛（Ce-SBA-15/KI）中时想进行两个反应物吸附能的计算。
      在计算过程中，催化剂采用PBE0方法和混合基组，其中I用sdd赝势，Ce用sdd下的大核版本MWB47，在统计能量的过程中统计EE + Thermal Free Energy Correction并将单位转换为kcal/mol。吸附能的计算采用E(吸附)=E（PO+catalyst)-E(PO)-E(Catalyst）。在计算完PO的吸附后，我们尝试用PO的结构近似构建BO吸附在分子筛中的结构，但计算完成后，二者的吸附能大相径庭。想求教一下各位产生这种情况的原因。
      附件为计算过程中产生的各种输入输出文件。

dnlx

1. 计算级别不一致
2. 单点能（EE）计算级别太低
3. 没有考虑色散校正
4. 如果改用3-zeta基组，需要考虑BSSE
5. 计算自由能的时候，不应该包括固体的平动和转动贡献。