CP2K做熔盐体系的AIMD偏慢，如何改进？（100个原子左右的金属熔盐体系）

duoconsole

各位老师好，我目前正基于文献重现熔盐体系的计算，过程中遇到内存不足与计算效率低下的问题，恳请各位帮忙分析是否为输入文件设置不当所致，具体情况如下：
一、模拟体系与整体流程

• 体系组成：含 36 个 Li⁺、25 个 K⁺、61 个 Cl⁻的熔盐体系（与文献一致）；
• 模拟流程：
  
  • 第一步：使用 CP2K 软件，采用 PIM 力场进行分子动力学（MD）模拟，时长 1ns；
  • 第二步：从第一步平衡后的 MD 轨迹中提取一帧结构，作为初始结构，进行从头算分子动力学（AIMD）模拟。
    
  
  

二、遇到的具体问题
模拟在超算平台上运行，主要出现两类问题：

• 频繁中断，需依赖 restart 续算：AIMD 每计算约 1 万步，就会触发报错（怀疑是内存不足），虽可通过 restart 文件继续计算，但频繁中断；
• 计算效率持续下降，耗时远超预期：初始阶段 AIMD 单步耗时10s左右，但计算至约 5 万步时，单步耗时已增至约 50 秒；按此效率，要完成文献中要求的 11 万步 AIMD，总耗时过长。
  

恳请各位老师帮忙排查输入文件是否存在参数设置不合理的情况，给出提升 AIMD 效率的建议，非常感谢！
附：（1）输入文件（自己借助multiwfn编写）；（2）参考文献（Nguyen H H, Gibson L D, Emerson M S, et al. Chlorine gas and anion radical reactivity in molten salts and the link to chlorobasicity[J]. Physical Chemistry Chemical Physics, 2025, 27(8): 4290-4297）.

ultramanm87

你AIMD效率逐渐下降可能编译有问题，内存有泄露啥的。
cp2k参数你可以尝试使用DZVP-MOLOPT-SR-GTH，测试下截断能能不能降一些。

sobereva

内存泄露问题必须解决

如果不是非得重复文献，基组可以降到DZVP-MOLOPT-GTH档次，REL_CUTOFF也可以适当减小点，并且建议用CSVR热浴