求助高斯固定部分结构的自由能计算问题

guoguoCO2

高斯计算： 不固定2D-MOF结构（取MOF中的重复单元作为MOF）优化，吸附HMF后结构变形厉害，不再是接近平面的2D结构
                 所以固定2D-MOF结构优化（如图脚本，固定2D-MOF，不固定HMF，吸附HMF后的结构和MOF自身结构的自由能校准值变化大（蓝色框出来的），所以计算吸附前后结构的吸附能是-1000多KJ/mol（红色框出来的），不知道那里存在问题？求助老师解答，非常感谢！

sobereva

不要把输入文件叫做“脚本”

那叫自由能热校正量，不叫“自由能校准值”

不要把MOF所有原子都冻住，只把边缘的原子冻了，能维持住实际MOF的骨架形态就足够了

冻住的原子对Hessian没贡献，显然冻结必然会直接影响热力学校正量。求差的时候冻结设置必须完全相同，才能最大程度抵消冻结带来的的人为影响。