计算化学公社

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

作者
Author: guoguoping199 时间: yesterday 21:13
标题: CP2K频率分析计算总是报错求助（NiFe-NC ORR反应中间体，共160个原子）
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核心计算，结果依然报错。请问各位老师，这是软件出现通信器泄露问题了还是我的输入文件有问题啊？
输入文件和输出文件如下：
(, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 1)

作者
Author: Uus/pMeC6H4-/キ 时间: 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号原子完全固定这点。我觉得把所有基底原子冻结、只考虑被吸附原子是不合理的，对几何优化和振动分析皆如此。

作者
Author: guoguoping199 时间: yesterday 22:14

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

原来如此，谢谢老师

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