thermo_pw如何计算吸附分子的吉布斯自由能？

Exodus

这是我目前的做法：准备优化后的*CO文件，在thermo_control中采用what='scf_disp',在ph_control中写入nat_todo=2，并在末尾处标记C,O的原子编号。在out文件中可以找到*CO在gamma点的频率 ，同时在热力学数据中还输出了如下信息 ，第二列似乎是U的校正，我验证了U=ZPE+CvdT，第三列应该是Helmholtz自由能的校正，我验证了也是相等的A=U-TS。
我的问题是：（1）由G=H-TS可知G=E+ZPE+CpdT-TS,我在一篇文章中看到 ，请问我可以按该公式将Cv换算为Cp吗？原文链接Lattice dynamics and thermophysical properties of h.c.p. Os and Ru from the quasi-harmonic approximation - IOPscience
（2）在输出的热力学数据中热容的贡献似乎是考虑的slab的，在298K时，CvdT=0.091Ry，这个数据是否太大了，我是否还应该计算slab的CvdT,然后在用*CO的减去slab的，以此得到*CO真实的CvdT。
谢谢各位~

卡开发发

1、可以，一些体胀系数不大的固体的Cp和Cv比较接近。
2、不需要。实际上我们可能会更关注吸附过程的自由能变，然后dG=G(CO@Slab)-G(CO)-G(Slab)之后这个值才有比较大小的价值，注意CO可能要手工求其平动转动对自由能贡献并进行加和。另外你声子频率的计算还有虚频，会对最终数据有影响。

Exodus

卡开发发 发表于 2021-6-24 15:10
1、可以，一些体胀系数不大的固体的Cp和Cv比较接近。
2、不需要。实际上我们可能会更关注吸附过程的自由能 ...

好的，谢谢您！