CP2K的thermochemistry和gaussian的对应关系求助

wgbigfatcat · 发表于 Post on 2020-5-16 15:29:04

用CP2K计算了一个小分子。得到它的化学热力学量如下，请假它们分别跟gaussian的对应是否正确？

VIB|                      NORMAL MODES - THERMOCHEMICAL DATA
VIB|
VIB|             Symmetry number:                                              1
VIB|             Temperature [K]:                                           298.15
VIB|             Pressure [Pa]:                                        101325.00


VIB|             Electronic energy (U) [kJ/mol]:                -412530.21514465             <===E,tot
VIB|             Zero-point correction [kJ/mol]:                      286.47778013                <===ZPE
VIB|             Entropy [kJ/(mol K)]:                                  0.24795733                      <===S,tot
VIB|             Enthalpy correction (H-U) [kJ/mol]:                299.11877477             <===H,corr
VIB|             Gibbs energy correction [kJ/mol]:                   225.19029687             <===G,corr
VIB|             Heat capacity [kJ/(mol*K)]:                            0.06693652

那么Sum of electronic and thermal free energies = E,tot + G,corr  ？？？？？？？

sobereva · 发表于 Post on 2020-5-16 18:03:52

wgbigfatcat · 发表于 Post on 2020-5-16 21:58:40

thx，sob。

CCCC · 发表于 Post on 2020-5-18 17:12:32

你好，请问可以看看你的计算输入文件吗？

wgbigfatcat · 发表于 Post on 2020-5-18 18:02:18

CCCC 发表于 2020-5-18 17:12
你好，请问可以看看你的计算输入文件吗？

我的计算结果有点问题。解决后发。

CCCC · 发表于 Post on 2020-5-19 10:55:45

wgbigfatcat 发表于 2020-5-18 18:02
我的计算结果有点问题。解决后发。

好的，谢谢分享

tuzhidingdong · 发表于 Post on 2020-11-1 20:10:35

请教一下，算单点能的时候就是算这个吗？Sum of electronic and thermal free energies = E,tot + G,corr

sobereva · 发表于 Post on 2021-2-9 09:01:54

tuzhidingdong 发表于 2020-11-1 20:10
请教一下，算单点能的时候就是算这个吗？Sum of electronic and thermal free energies = E,tot + G,corr

明显不是。这是自由能

wuandsun · 发表于 Post on 2022-1-14 19:38:46

wgbigfatcat 发表于 2020-5-18 18:02
我的计算结果有点问题。解决后发。

请问解决了吗，可以发一下输入文件吗，谢谢

sobereva · 发表于 Post on 2022-1-15 07:58:01

wuandsun 发表于 2022-1-14 19:38
请问解决了吗，可以发一下输入文件吗，谢谢

振动分析的输入文件用Multiwfn直接就能产生
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

Shermo可以直接基于CP2K和Gaussian振动分析任务的输出文件计算热力学量，输出是统一的，所以都不需要考虑对应关系问题了
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

elpa · 发表于 Post on 2022-5-10 23:04:25

请问楼主input文件是怎么设置打印热力学数据的？为什么我的输出有两个NaN。

Electronic energy (U) [kJ/mol]:                -10342166.27894069
Zero-point correction [kJ/mol]:                   6763.50262481
Entropy [kJ/(mol K)]:                                        NaN
Enthalpy correction (H-U) [kJ/mol]:                7455.51038792
Gibbs energy correction [kJ/mol]:                            NaN
Heat capacity [kJ/(mol*K)]:                            3.91701734

ChemG · 发表于 Post on 2022-5-11 09:28:11

elpa 发表于 2022-5-10 23:04
请问楼主input文件是怎么设置打印热力学数据的？为什么我的输出有两个NaN。

Electronic energy (U) [kJ ...

可以用Shermo处理振动分析的out文件，最终格式与Gaussian一致。之前也遇到过这种情况，Shermo就可以显示出来

elpa · 发表于 Post on 2022-5-11 17:41:33

ChemG 发表于 2022-5-11 09:28
可以用Shermo处理振动分析的out文件，最终格式与Gaussian一致。之前也遇到过这种情况，Shermo就可以显示 ...

谢谢，真的可以

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

[CP2K] CP2K的thermochemistry和gaussian的对应关系求助

浏览过的版块