CP2K给192个原子的Na2CO3超胞做振动分析崩溃退出

北大-陶豫

求助，我用CP2K给192个原子的Na2CO3超胞做振动分析，有限差分算到第75个原子的时候崩溃退出了。我现在脚本是这样的：

#SBATCH -n 96
source /public1/soft/modules/module.sh
module load cp2k/2023
mpirun -np 96 cp2k.popt -o xxx-freq.out xxx-freq.inp

inp 里设置了 NPROC_REP 24（结构优化时用的24个核，#SBATCH -n 24 和 mpirun -np 24，能正常跑完）。输入文件见附件，理论级别 PBE-D3 / pob-DZVP-rev2。
试过将 mpirun -np 96 改为其它核数（#SBATCH -n 96 这一行的核数没有修改），72核、48核也会爆掉但爆掉之前有限差分能多算几个原子，24核能一直不报错，但会卡在第90个左右原子的有限差分，之后好几天都没有算到下一个原子，我就把任务杀了。

我也尝试了用 psmp 方法，但是进入第一轮SCF的时候就会崩溃退出。我用 psmp 方法的脚本如下：

#SBATCH -n 96
source /public1/soft/modules/module.sh
module load cp2k/2023
export OMP_NUM_THREADS=24
mpirun -np 4 cp2k.psmp -o xxx-freq.out xxx-freq.inp

应该怎么办比较好呢？