500个原子的吸附模型CP2K感觉算起来好慢，不知道这个速度正常不正常

ZeeWong

我是在超算上面装的CP2K 2023
之前一直用的第一次CP2K培训的案例文件，能量部分的计算时间和卢老师给的文件大差不差。
这次自己弄了一个体系，500个Fe原子吸附一个有机分子，500多个原子，用的GFN1-xtb，算的速度好慢，7个小时结构优化才走了一步（超算给的16核）。
不知道是我计算设置的问题还是还是软件安装的问题。

sobereva

绝对不要用PSOLVER ANALYTIC。一定不要用太老版本的Multiwfn。较新版本的Multiwfn对二维周期性体系一律用MT，比ANALYTIC强多了。

slab完全没必要建得这么大，极度浪费时间。用小点slab并结合k点效率高得多

Fe的slab是具有磁性的，参考北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）里“磁性体系的计算”那部分的讲解

GFN1-xTB跑这种体系是否理想不一定，如果没有类似体系的研究文章做支撑，可能会被审稿人质疑。但凡DFT算得动就用DFT

Fe slab是导体，smearing显然必须开。你当前几何优化每一步做了巨量次数的outer SCF才收敛，显然效率太低了。一定要注意监控

下文仔细看看
解决CP2K中SCF不收敛的方法
http://sobereva.com/665（http://bbs.keinsci.com/thread-37196-1-1.html）

如果用的psmp应当用popt，并且当前体系很大，多给点核并行

ZeeWong

sobereva 发表于 2023-7-9 03:42
绝对不要用PSOLVER ANALYTIC。一定不要用太老版本的Multiwfn。较新版本的Multiwfn对二维周期性体系一律用MT ...

谢谢卢老师！！！我再试试！