计算化学公社

标题: 求助：CP2K计算H2分子的频率与VASP中计算结果不一致 [打印本页]

作者
Author: guyingbing 时间: 2022-6-11 14:38
标题: 求助：CP2K计算H2分子的频率与VASP中计算结果不一致
大家好！我是在学习CP2K的初学者。现在尝试使用cp2k计算H2分子的能量。采用的基组是：DZVP-MOLOPT-SR-GTH-q1 赝势： GTH-PBE-q1。CUTOFF为1000。收敛后，我发现H-H之间的键长为：0.721 A。低于实际的H-H的键长（0.75A），也与vasp计算得到的数值不同。此外，计算频率后，发现采用CP2K优化得到有3个振动，其波数分别为：106.856689 cm-1 118.894653 cm-1 5025.259903 cm-1。然而，采用vasp计算得到的波数分别为：49.622957 cm-1 56.075428 cm-1 4343.01995 cm-1。此外，分别采用Shermo 与vaspkit计算了相应的热力学校正量。我发现热力学校正量存在较大的差距（如下图）。能请问一下是不是由于cp2k计算过程中设定错误造成的？能请各位老师帮忙看一下吗？vasp的计算参数为：ENCUT = 800，ISPIN = 1。

作者
Author: sobereva 时间: 2022-6-11 17:09
没具体输入文件没法说
CP2K算H2没理由用这么小的基组，应当用尽可能好的TZV2P-MOLOPT-GTH，或者算分子体系常用的诸如def2-TZVPP

作者
Author: guyingbing 时间: 2022-6-11 17:38
感谢老师的回复。我上传了我的输入和输出文件。请您审阅。我打算利用CP2K计算水分子在相关氧化物上的吸附能。所以采用了BASIS_MOLOPT中的DZVP-MOLOPT-SR-GTH-q1。

作者
Author: sobereva 时间: 2022-6-11 18:05

guyingbing 发表于 2022-6-11 17:38
感谢老师的回复。我上传了我的输入和输出文件。请您审阅。我打算利用CP2K计算水分子在相关氧化物上的吸附能 ...

换成更好的基组再说

你算的是水分子又不是氢分子，对氢分子描述差的基组不代表对水分子描述差

你的输入文件里还净是无意义的设置，诸如UKS T，体系又不是开壳层的写它干嘛。按下文几键就能产生又干净又合理的输入文件
使用Multiwfn非常便利地创建CP2K程序的输入文件
http://sobereva.com/587（http://bbs.keinsci.com/thread-21668-1-1.html）

作者
Author: guyingbing 时间: 2022-6-11 19:57
本帖最后由 guyingbing 于 2022-6-11 20:09 编辑

感谢sob老师的指点。通过换用基组TZV2PX-MOLOPT-GTH-q1。现在计算出的热力学校正量已经和vasp计算的基本一致，仅仅小数点后第二位有差距。具体结果如图：请老师审阅。

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