计算化学公社

标题: cp2k做aimd这个速度正常吗？ [打印本页]

作者
Author: Alset 时间: 2024-10-18 16:32
标题: cp2k做aimd这个速度正常吗？
老师们好，

目前我在采用cp2k进行aimd模拟，体系是包含8个水分子的冰，使用泛函为vdw-df2，基组为TZV2P-MOLOPT-GTH，ENCUT为800 Ry，k点为 6 6 4（由vaspkit产生的) ，使用32个核，并使用cp2k.potp -np 32 并行。模拟速度大约是1分钟/fs，如果想要模拟10 ps，大概需要一周的时间。这速度似乎远远比论坛上所说到的 “cp2k可以跑上百原子的aimd到几十上百皮秒” 要慢得多，所以我想向大家请教下，这里较慢的速度可能是我使用了较精确的基组泛函，还是其他因素所导致的呢？

附件的压缩包中包括了输入、输出、以及提交任务时所使用的脚本。

作者
Author: sobereva 时间: 2024-10-18 16:51
直接用CP2K基于标准的MP方法产生k点就完了，不需要借助其它工具
CP2K里没有叫ENCUT的东西

你的模拟体系完全不合理，盒子尺寸这么小，动力学行为会受到严重的自相互作用影响，结果根本没意义。用足够大的盒子尺寸时（如>=12埃）对此体系根本不需要考虑k点，此时小几百个原子，AIMD跑100步后每一步平均在普通双路服务器上花<6秒钟属于正常耗时。
基组完全没必要这么大。TZVP-MOLOPT-SR-GTH足矣
对于不控压的情况，400 Ry的CUTOFF就够，至多用到600也绰绰有余了

建议模拟前先参考前人的水的AIMD模拟文章，比如Chem. Sci., 2017, 8, 3554–3565，就是用的CP2K跑的

作者
Author: Alset 时间: 2024-10-18 17:03

sobereva 发表于 2024-10-18 16:51
直接用CP2K基于标准的MP方法产生k点就完了，不需要借助其它工具
CP2K里没有叫ENCUT的东西

好的老师，谢谢！

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