准备用CP2K计算一个150个原子结构内水分子扩散的MSD，cutoff和realcutoff怎么选择

yang543

准备计算一个150个原子体系的分子筛结构内，多种小分子的扩散系数，之前用vasp做过DFT计算，那些截断能和k点应该不能直接拿来用吧。所以就标准流程来说这个的cutoff和realcutoff应该怎么确定呢？

sobereva

那叫rel_cutoff

k点可以延用。但由于VASP跑不小的体系的计算效率比CP2K低太多太多，所以VASP跑AIMD时往往被迫用过小的盒子将就，而效率高得多的CP2K完全可以用更大盒子得到更可靠结果，此时跑AIMD的盒子一般已经大到不需要考虑k点。

cutoff、rel_cutoff拿不准就做收敛性测试，这样也最严格。懒得测试图省事就根据用的GTH赝势基组的最大GTF指数估计，北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/KFP）里面“能量的计算及相关问题” 那一节的 “k点的选择和收敛性测试” 和 “截断能的选择和收敛性测试” 部分专门详细讲了各种细节，并且还给了测试脚本。还专有一节“第一性原理分子动力学” 给了巨量可以参考的AIMD例子。