计算化学公社

标题: 求助如何对shermo计算得出的热力学量进行压力矫正的问题 [打印本页]

作者
Author: shibanfish 时间: 2026-4-24 17:14
标题: 求助如何对shermo计算得出的热力学量进行压力矫正的问题
@sobereva 卢老师，我之前用您开发的shermo计算一下CO在Pt表面氧化过程过程中所涉及的反应物、生成物分子在Pt表面的吸附吉布斯自由能变，发现shermo好像算不了低压时的热力学量(我研究的压强是在1E-8到10bar这个区间的，我输进去压强1E-8，shermo只能给出压强为0的数据)，请问我该如何用1E-8到1E-4bar这个区间的压强来对shermo计算出的压强为0的数据进行修正呢？麻烦老师了
图中的是气态CO的热力学量

作者
Author: sobereva 时间: 2026-4-24 17:19
此处输出信息只保留三维小数，自然表面上看不出差异来。内部实际计算的时候你设多少就按多少计算

作者
Author: shibanfish 时间: 2026-4-24 19:58
本帖最后由 shibanfish 于 2026-4-25 20:49 编辑

好的，多谢卢老师解答。但我还想再问一下，如果我要算固体表面吸附气体分子构型的热力学量的话，电子能在CP2K的输出文件的哪里能找到呢？如果要计算吸附过程的吉布斯自由能变的话，直接用吸附态分子计算出的自由能(这个能量由固定固体表面原子，只放开气体分子并做振动分析后得到)减去气态分子计算得到的自由能就行吗？希望老师能解答一下@sobereva

作者
Author: shibanfish 时间: 2026-4-25 22:00

sobereva 发表于 2026-4-24 17:19
此处输出信息只保留三维小数，自然表面上看不出差异来。内部实际计算的时候你设多少就按多少计算

另外，我还想问一下，如果要计算吸附过程的吉布斯自由能变的话，直接用吸附态分子计算出的自由能(这个能量由固定固体表面原子，只放开气体分子并做振动分析后得到)减去气态分子计算得到的自由能就行吗？麻烦老师解答一下

作者
Author: sobereva 时间: 2026-4-26 05:04

shibanfish 发表于 2026-4-25 22:00
另外，我还想问一下，如果要计算吸附过程的吉布斯自由能变的话，直接用吸附态分子计算出的自由能(这个能 ...

否
和被吸附分子有显著相互作用的固体表面原子不能冻住
哪怕你不考虑固体表面原子在吸附过程中结构改变，至少你还得再减去固体表面的自由能才行

作者
Author: shibanfish 时间: 2026-4-26 11:24

sobereva 发表于 2026-4-26 05:04
否
和被吸附分子有显著相互作用的固体表面原子不能冻住
哪怕你不考虑固体表面原子在吸附过程中结构改变 ...

好的，感谢卢老师的解答。我想知道，如果需要放开固体表面原子的话，只需要放开与气体分子接触的最外层固体原子就行是吗？还是说要把所有固体原子都放开？如果是只需要放开最外层固体原子的话，后面计算吸附吉布斯自由能变的时候只需要G吸附=G(放开最外层固体原子+吸附气体分子)-G(放开固体最外层原子)-G(气相气体分子)就行是吗？

作者
Author: sobereva 时间: 2026-4-26 15:14

shibanfish 发表于 2026-4-26 11:24
好的，感谢卢老师的解答。我想知道，如果需要放开固体表面原子的话，只需要放开与气体分子接触的最外层固 ...

吸附可能造成结构不可忽视改变的固体表面原子都得放开
是

作者
Author: shibanfish 时间: 2026-4-26 15:36

sobereva 发表于 2026-4-26 15:14
吸附可能造成结构不可忽视改变的固体表面原子都得放开
是

好的，感谢卢老师在百忙之中解答我的问题

欢迎光临计算化学公社 (http://bbs.keinsci.com/)