求助 表面计算化学与CASTEP Thermaldynamic analysis 数据的疑问。

Toussaint

题主在计算表面吸附能和催化反应的计算，需要计算吸附分子反应前后吉布斯自由能变化，用到的公式和（求助电催化中吉布斯自由能的计算问题 - 量子化学 (Quantum Chemistry) - 计算化学公社 (keinsci.com)）中提到的公式相同，如图

落实到CASTEP的操作层面，有几个主要的问题希望有经验的老师可以帮忙解答一下：
1.在CASTEP中一般是用几何优化完的吸附模型（2*2的8层表面层+1个吸附分子，原子数>120）来计算声子谱吗？
2.计算phonons时，选energy和properties两个选项的主要区别是什么？一般推荐用哪个？
3.analysis -> phonons.castep 文件可以直接得到下图所示的 ZPE，free energy等信息，题主查资料大概明白ΔS部分，也就是振动对熵值的贡献是通过下面的公式（对phonon density积分）计算得到的，ZPE同理，是一个常数，那么图中的Enthalpy部分是如何得来的？是振动引起的焓值吗？不太明白其物理意义，而且为何在上图的公式ΔG=ΔE+ΔEZPE-TΔS中为什么没有出现？

4.计算ΔG=ΔE(total energy)+ΔEZPE（analysis得到的数据）+ΔFree energy(图中红线数据)正确吗？

IreneAlder

我也遇到了相同的问题，计算出来吸附吉布斯自由能随着温度升高反而越来越小