计算化学公社
标题: 计算气相热力学参数遇到的问题 [打印本页]
作者Author: 雷鑫欣 时间: 2025-8-14 15:01
标题: 计算气相热力学参数遇到的问题
我在利用shermo计算热力学参数时，shermo通过读取高斯计算得到的结构优化与频率计算的输出文件，即可生成热力学数据。但是像计算SF6，SF5的热力学数据时与数据库的差异很小，但是在计算SF,S2时与数据库的差异就大。而且在计算SF+,SF-,S2+,S+,S-这样的离子时，得到的热力学数据与数据库的差异也大（SF+的结构优化与频率计算就是在SF的基础上改了电子数与多重度在高斯中进行计算的，所以的离子都是这么算的），到底是什么原因呢。Shermo是通过频率信息进行计算的，那像S,F这样的原子没有振动频率，如果要计算热力学参数该怎么计算呢。以及S+,F+,SF+,SF-,S2+,S-,F-这样的离子该怎样计算呢。恳请老师解答。
热力学数据库的数据参照的是这个数据库：NASA Thermo Build (beta)
作者Author: sobereva 时间: 2025-8-14 20:28
先交代清楚你说的热力学数据具体指什么
作者Author: 雷鑫欣 时间: 2025-8-20 15:13
CP(cal/mol/K) S(cal/mol/K) -(G° - H°298.15)/T H° - H°298.15
是算的这四个数据，恳请sob老师解答
作者Author: sobereva 时间: 2025-8-21 03:15
注意检查净电荷、自旋多重度确保正确，确保得到的是稳定的波函数，以及确保计算用的级别合理
旋轨耦合会造成其中一些体系简并的基态能级发生分裂，让Shermo恰当考虑这点结果会更好（参考Shermo手册Consideration of degeneracy of ground state and spin-orbit coupling splitting部分）
作者Author: 雷鑫欣 时间: 2025-8-21 13:53
多谢老师
作者Author: 雷鑫欣 时间: 2025-9-9 15:16
老师我想请教一下，我在shermo文件夹中的settings中进行了修改，加入了SOC部分也就是（您之前说的旋轨耦合会造成其中一些体系简并的基态能级发生分裂这一方面），但是运行shermo后shermo为什么无法读取到这部分的信息呢。settings的编写如下：
E= -795.455115
prtvib= 1
T=300,30000,100
P= 1.0
sclZPE= 0.97
sclheat= 1.0
sclS= 1.0
sclCV= 1.0
ilowfreq= 2
ravib= 100
imode= 0
conc= 0
outshm= 0
PGlabel= ?
defmass= 3
modmass
ElecLev 0.0 1
ElecLev 35.0 1
ElecLev 72.0 1
（注意： ElecLev这块的数据并不准，只是为了验证shermo能否读取这部分信息）恳请老师解答
作者Author: sobereva 时间: 2025-9-9 19:51
不要自己乱改settings.ini的内容。需要改什么就改什么，其它保留默认不动
Shermo根本没有ElecLev这种写法
仔细看启动Shermo时屏幕上的提示有没有warning说没找到settings.ini
作者Author: 雷鑫欣 时间: 2025-9-10 09:03
好的，谢谢老师
| 欢迎光临 计算化学公社 (http://bbs.keinsci.com/) |
Powered by Discuz! X3.3 |