Uus/pMeC6H4-/キ 发表于 2024-11-29 09:49
&FORCE_EVAL/&DFT/&CHARGE是设系统总电荷的关键词。但这只是上面列的输入文件的一大堆问题之一，只改这个 ...

好嘞，明白，谢谢老师，我也没敢多算，还在学习参数

462855 发表于 2024-11-29 08:57
请老师，是哪个参数设置的电荷减1呢？现在的结果来看是应该对电荷进行设置的，只是一直没找到

&FORCE_EVAL/&DFT/&CHARGE是设系统总电荷的关键词。但这只是上面列的输入文件的一大堆问题之一，只改这个远不够。

一楼列的这篇文献我也简单看过。CP2K跑分子动力学的常见计算级别没有那么高，也就纯泛函PBE-D3结合2-zeta档次的DZVP-MOLOPT-SR-GTH以及设定恰当的截断能而已，虽然粗糙但计算速度够快；就算不直接用这个级别的结果，也至少可以确认流程和设置是否合理、让自己清楚知道各关键词的意义。而文献涉及的杂化泛函有更多讲究，参考社长博文http://sobereva.com/690或者科音培训班的讲义。

总之，如果没有低计算级别的AIMD计算经验，不要贸然模仿文献操作投入计算资源挑战高计算级别，否则只会像一楼那样浪费大把时间。

Uus/pMeC6H4-/キ 发表于 2024-11-28 20:54
这个水溶液体系是怎么构建的，有没有先几何优化和预平衡？Nose-Hoover热浴及其大到500的时间常数是怎么选的 ...

请老师，是哪个参数设置的电荷减1呢？现在的结果来看是应该对电荷进行设置的，只是一直没找到

这个水溶液体系是怎么构建的，有没有先几何优化和预平衡？Nose-Hoover热浴及其大到500的时间常数是怎么选的？为什么对只有一个氢氧根离子没有阳离子的电荷不平衡体系，不设体系总电荷-1却要设UKS T？MAX_MEMORY 37000如此巨大，符合实际运行CP2K的计算机的内存容量吗？

虽然体系不算复杂，但是刚开始用CP2K跑AIMD就要复现文献，还是用的revPBE0杂化泛函加TZV2P这么大的基组，实在是迈步大了点。文献的计算设置合不合理是要仔细考量的，就算合理也未见得可以完全照搬到自己的计算条件。推荐参加一个北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）了解背后的各种门道。