求助如何对shermo计算得出的热力学量进行压力矫正的问题

shibanfish · 发表于 Post on 2026-4-24 17:14:48

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

sobereva · 发表于 Post on 2026-4-24 17:19:41

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

shibanfish · 发表于 Post on 2026-4-24 19:58:40

本帖最后由 shibanfish 于 2026-4-25 20:49 编辑

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

shibanfish · 发表于 Post on 2026-4-25 22:00:54

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

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

sobereva · 发表于 Post on 2026-4-26 05:04:04

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

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

shibanfish · 发表于 Post on 2026-4-26 11:24:54

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

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

sobereva · 发表于 Post on 2026-4-26 15:14:45

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

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

shibanfish · 发表于 Post on 2026-4-26 15:36:55

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[辅助/分析程序] 求助如何对shermo计算得出的热力学量进行压力矫正的问题