计算化学公社

标题: CP2K的thermochemistry和gaussian的对应关系求助 [打印本页]

作者
Author: wgbigfatcat 时间: 2020-5-16 15:29
标题: CP2K的thermochemistry和gaussian的对应关系求助
用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  ？？？？？？？

作者
Author: sobereva 时间: 2020-5-16 18:03
对

作者
Author: wgbigfatcat 时间: 2020-5-16 21:58
thx，sob。

作者
Author: CCCC 时间: 2020-5-18 17:12
你好，请问可以看看你的计算输入文件吗？

作者
Author: wgbigfatcat 时间: 2020-5-18 18:02

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

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

作者
Author: CCCC 时间: 2020-5-19 10:55

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

好的，谢谢分享

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

作者
Author: sobereva 时间: 2021-2-9 09:01

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

明显不是。这是自由能

作者
Author: wuandsun 时间: 2022-1-14 19:38

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

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

作者
Author: sobereva 时间: 2022-1-15 07:58

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）

作者
Author: elpa 时间: 2022-5-10 23:04
请问楼主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

作者
Author: ChemG 时间: 2022-5-11 09:28

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

Electronic energy (U) [kJ ...

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

作者
Author: elpa 时间: 2022-5-11 17:41

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

谢谢，真的可以

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