求助使用dmol3计算分子筛团簇时吉布斯自由能的计算

Goodkkk

本人做的是一个在周期性结构模型中优化后截出来的19T大小的带质子酸位点的团簇，在将除中间5T部分的活动区域固定笛卡尔坐标后，直接加入烯烃进行优化和过渡态搜索，但是TS opt得到的过渡态结果常有冻结原子部分的虚频（分子筛骨架的振动），但是活性区域只有一个虚频而且是我想要的虚频。现在审稿人要求补充吉布斯自由能的计算，我想求助这时候我是应该算Partial hessian还是在这种情况下继续计算吉布斯自由能。