想请教类似的问题。我也正在学习cp2k计算吉布斯自由能，我准备练习的是纯Pt体系和Pt吸附CO后的体系，目前CP2K讲义有两种方法，一种使用Shermo，看着挺方便，但是好像适合大一些的晶体，Pt体系和Pt吸附CO后的体系扩胞成四层，晶胞在12埃左右，适合用Shermo计算吉布斯自由能吗？如果用第二种方法：CP2K结合phonopy算吉布斯自由能，但是看了也不明白所得结果怎么算吉布斯自由能。phonopy算的结果中可以找到亥姆霍兹自由能校正，熵，等容热容，内能校正，ZPE 。可不可以也像楼主一样将亥姆霍兹自由能校正当做吉布斯自由能校正量，再加上ZPE+电子的能量，这三项就是体系的吉布斯自由能？如果Pt吸附CO后的体系中亥姆霍兹自由能校正不能当做吉布斯自由能校正量，又该怎么算呢？

sobereva 发表于 2025-6-17 04:44
1 是

2 什么叫“造成这个原因”，哪个原因，把话说明白

明白了，谢谢sob老师

sobereva 发表于 2025-6-17 04:39
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

收到，谢谢sob老师。

1 是

2 什么叫“造成这个原因”，哪个原因，把话说明白
记住热容是依赖于温度的，对T的积分没法简化为常数乘以温度。热力学量的计算公式看下文介绍的Shermo程序手册附录部分
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “CP2K结合phonopy做自由能校正” 改了，以后务必注意！