计算化学公社

标题: CP2K怎样生成焓熵自由能一类的文件 [打印本页]

作者
Author: Tiwmdh 时间: 2023-4-18 19:08
标题: CP2K怎样生成焓熵自由能一类的文件
请问一下，我在用CP2K时候，在out文件里面始终没有找到这种数值相关的内容，这个是需要.inp文件里面有什么设置吗，还是需要什么特殊的软件

作者
Author: sobereva 时间: 2023-4-19 00:46
CP2K做振动分析，然后Shermo载入其输出文件，按下文即可得到各种热力学量，非常简单
使用Shermo结合量子化学程序方便地计算分子的各种热力学数据
http://sobereva.com/552（http://bbs.keinsci.com/thread-17494-1-1.html）

对于周期性体系，记得settings.ini里imode设1去掉平动和转动贡献

如果是小晶胞的晶体，诸如CaO，则更适合用phonopy算热力学量，远远不需要对超胞做6*Natom次非常昂贵的受力计算，不过需要先安装，而且操作步骤远比用Shermo多。

作者
Author: Tiwmdh 时间: 2023-4-19 08:52

sobereva 发表于 2023-4-19 00:46
CP2K做振动分析，然后Shermo载入其输出文件，按下文即可得到各种热力学量，非常简单
使用Shermo结合量子化 ...

感谢sob老师

作者
Author: Tiwmdh 时间: 2023-4-19 15:58

sobereva 发表于 2023-4-19 00:46
CP2K做振动分析，然后Shermo载入其输出文件，按下文即可得到各种热力学量，非常简单
使用Shermo结合量子化 ...

sob老师，我还有个问题想问下，我用Multiwfn将输出的inp文件改为vibrational——analytic模式，然后模拟运算之后得到的振动能的贡献为0，这个是什么情况呢

作者
Author: sobereva 时间: 2023-4-20 04:24

Tiwmdh 发表于 2023-4-19 15:58
sob老师，我还有个问题想问下，我用Multiwfn将输出的inp文件改为vibrational——analytic模式，然后模拟 ...

看不懂你的意思。从out文件明显可见当前任务都还没算完

作者
Author: Tiwmdh 时间: 2023-4-20 15:00

sobereva 发表于 2023-4-20 04:24
看不懂你的意思。从out文件明显可见当前任务都还没算完

抱歉sob老师，表述不清，就是我按照步骤求得的vibratio项所有的数值都是0，输入文件是Multiwfn生成的，如果没跑完，那我再找找别的问题

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