CP2K频率分析计算总是报错求助（NiFe-NC ORR反应中间体，共160个原子）

guoguoping199 · 发表于 Post on yesterday 21:13

cp2k版本：2025.1
计算节点配置：AMP EPYC 7452 x2，共64核心512线程，为租的超算节点，软件安装也是超算方面工程师安装配置。
输入文件使用Multiwfn_2026.1.12_bin_Win64生成
求助问题：使用该节点计算OPT和单点能都正常，但是频率分析总是报错。
任务描述：进行结构优化后，然后进行频率分析
计算结构：NiFe-NC ORR反应中间体，共160个原子
所用基组：GFN1-xTB
使用核心：16核心
最开始使用64核心计算报错，搜索发现好像是内存不够导致的，于是改为使用16核心计算，结果依然报错。请问各位老师，这是软件出现通信器泄露问题了还是我的输入文件有问题啊？
输入文件和输出文件如下：

F.inp (12.36 KB, 下载次数 Times of downloads: 1)

F.out (49.18 KB, 下载次数 Times of downloads: 1)

Uus/pMeC6H4-/キ · 发表于 Post on yesterday 21:52

振动分析任务不能只看主输出文件而忽略各个副本的输出文件。我用mpirun -np 4 -x OMP_NUM_THREADS=2 cp2k.psmp和较少内存来跑这个任务，最后异常终止并不是因为内存占用过多，而是有一个副本SCF不收敛。这就自然引出问题：为什么要把这个体系当成电荷0自旋多重度1的闭壳层的样子，用GFN1-xTB来跑？该理论方法在此类体系上未必经过测试，再者CP2K本来也尚未支持明确考虑了自旋极化的spGFN-xTB形式，不宜以换理论方式来逃避如常规DFT的SCF不收敛等问题（类似http://bbs.keinsci.com/thread-55778-1-1.html的问题）。

另外，介绍体系尺寸不能只说总共160个原子而不提1-158号原子完全固定这点。我觉得把所有基底原子冻结、只考虑被吸附原子是不合理的，对几何优化和振动分析皆如此。

guoguoping199 · 发表于 Post on yesterday 22:14

Uus/pMeC6H4-/キ发表于 2026-5-16 21:52
振动分析任务不能只看主输出文件而忽略各个副本的输出文件。我用mpirun -np 4 -x OMP_NUM_THREADS=2 cp2k.p ...

原来如此，谢谢老师

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

[CP2K] CP2K频率分析计算总是报错求助（NiFe-NC ORR反应中间体，共160个原子）