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

应该怎么办比较好呢？

北大-陶豫

喝酸奶的巧乐兹 发表于 2026-5-6 23:41
用psmp版本，这么写试试呢
mpirun -np 8 --cpus-per-proc 12 -x OMP_NUM_THREADS=12 cp2k.psmp -i xxx.inp ...

谢谢，我又试了一下，你那种方法不是很管用（捂脸），我是这样解决的：改为 module load cp2k/2023-openmpi，脚本如下：

#!/bin/bash
#SBATCH -p v6_384
#SBATCH -N 1
#SBATCH -n 96
#SBATCH -J Li2CO3-222-freq
#SBATCH -o jobid%j-%N.out
#SBATCH -e jobid%j-%N.error

source /public1/soft/modules/module.sh
module load cp2k/2023-openmpi

export OMP_NUM_THREADS=24
mpirun -np 4 cp2k.psmp -o xxx-freq.out xxx-freq.inp

喝酸奶的巧乐兹

用psmp版本，这么写试试呢
mpirun -np 8 --cpus-per-proc 12 -x OMP_NUM_THREADS=12 cp2k.psmp -i xxx.inp -o xxxx.out