CP2K 中 PBE/DZVP 优化后使用 PBE0/TZVP 进行计算报错

yyyyang

老师们好，

最近我使用 CP2K 在 PBE/DZVP-MOLOPT-SR-GTH 水平下优化了 Pt(111) 表面上氢分子的结构，并在能量收敛后得到了结构 B。随后我想在 PBE0/TZVP-MOLOPT-PBE-GTH 水平下对结构 B 进行单点能和频率计算。

计算单点能 output 报错： hfx_energy_potential.F:486。我尝试增大 cutoff energy，但仍然出现同样的报错。请问这种错误通常是什么原因导致的？

计算频率 没有输出的.mol文件计算就结束了。

另外，如果需要使用 PBE0 计算 铂表面的吸附能，有什么比较推荐的基组吗？
谢谢！

Uus/pMeC6H4-/キ

CP2K里杂化泛函带k点计算需要用&XC/&HF/&RI字段启用RI-HFx技术（为啥multiwfn创建的输入文件没考虑这点？），而且现实中即使只算单点也很难处理边长11埃、超过60个重金属原子的庞大高密度体系（除非有巨量的核心数与内存可调用），更别说要用有限差分法算6N个单点的振动分析任务了。

真想要PBE0算吸附能的话，可能还是像http://sobereva.com/540说的挖团簇模型在量子化学软件当孤立体系算比较可行（这样无需k点，但得自行按边界效应和吸附覆盖度要求等做团簇尺寸收敛性测试了）。

UW_0728.

CP2K的RI-HFXk，超过个位数原子根本甭考虑，不可能算得动

yyyyang

Uus/pMeC6H4-/キ 发表于 2026-3-9 18:04
CP2K里杂化泛函带k点计算需要用&XC/&HF/&RI字段启用RI-HFx技术（为啥multiwfn创建的输入文件没考虑这点？） ...

感谢回复。我已经使用 PBE0/def2-TZVP 计算过氢分子在 Pt 团簇上的吸附能，现在希望与 slab 模型的吸附能进行对比。但CP2K用 PBE/DZVP 计算得到的吸附能与之前结果差异较大，因此想尝试使用 PBE0 再做一次计算进行比较。请问如果使用 PBE/TZVP-MOLOPT-PBE-GTH 计算 slab 的吸附能，与团簇模型的结果进行对比，这样的差异是否可以接受？还是说为了保证结果的可比性，使用相同的泛函和基组进行计算会更合理？

UW_0728.

yyyyang 发表于 2026-3-9 18:58
感谢回复。我已经使用 PBE0/def2-TZVP 计算过氢分子在 Pt 团簇上的吸附能，现在希望与 slab 模型的吸附能 ...

就算吸附能而言，DZVP档次的基组不够用；PBE算吸附能还可以，但算不上好的。建议BEEF-vdW/TZV2P-MOLOPT(-SR)-GTH，或者便宜一些的RPBE-D3(BJ)或revPBE-D3(BJ)/TZV2P-MOLOPT(-SR)-GTH。

个人觉得这样的体系没必要用簇模型；如果一定要用并且要比较的话，确保团簇模型做了合理的收敛性测试并合理冻结了原子、避免了悬键

yyyyang

UW_0728. 发表于 2026-3-9 19:52
就算吸附能而言，DZVP档次的基组不够用；PBE算吸附能还可以，但算不上好的。建议BEEF-vdW/TZV2P-MOLOPT(-SR ...

收到！感谢回复！

sobereva

DZVP-MOLOPT-SR-GTH算吸附能必须考虑BSSE校正，否则误差巨大。如果不考虑BSSE校正，至少也要用不带SR的

参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）讲吸附能计算的一个例子的对比