计算化学公社

标题: MD收敛的帧，单独拿出来算SCF却不一定收敛 [打印本页]

作者
Author: a386514315 时间: 2024-6-16 11:30
标题: MD收敛的帧，单独拿出来算SCF却不一定收敛
本帖最后由 a386514315 于 2024-6-16 11:57 编辑

各位老师您们好！
本人在利用cp2k计算分子动力学后（20000步），抽出50帧单独算scf（每间隔400步抽一帧），看看是否收敛。但是实际上发现并不是所有的帧都收敛。因此想请教这是为什么？

1.正常收敛的帧，在收敛速度方面，在分子动力学中，往往经历一个OUTER_SCF就能收敛，而单独算scf则往往需要3个以上OUTER_SCF才能收敛。根据社长的帖子http://bbs.keinsci.com/thread-37196-1-1.html，MD过程中程序会自动利用之前一些结构的波函数外推出当前步的初猜波函数，因此这也合理。

2.但是，却出现了10帧左右数据，单独scf计算不能收敛，请问是否温度或者构型有关？

3. 参数设置方面，SCF和MD计算尽量保持一致：
理论方法为PBE，基组和赝势为DZVP-MOLOPT-SR-GTH
采用OT
Inner和Outer SCF的EPS_SCF为10E-5
CUTOFF、REL_CUTOFF分别为400,50
MD采用CSVR热浴，温度控制在300k。

输入和部分输出文件如附件，输入文件均采用Multiwfn生成。

作者
Author: sobereva 时间: 2024-6-16 19:13
AIMD过程每一步都会用ASPC方式基于之前多步的波函数信息外推产生本次用的初猜，质量比起直接CP2K重新产生的初猜好得多得多。也因此AIMD跑起来一定时间后平均每步的耗时比初期往往低得多

作者
Author: a386514315 时间: 2024-6-16 20:37

sobereva 发表于 2024-6-16 19:13
AIMD过程每一步都会用ASPC方式基于之前多步的波函数信息外推产生本次用的初猜，质量比起直接CP2K重新产生的 ...

感谢社长解答！请问这个问题如何解决呢？按照您的之前的帖子http://bbs.keinsci.com/thread-37196-1-1.html，请问我是提高CUTOFF，切换成对角化（设置对应的k点），以上两种方式我也试过了（而且担心和md跑出来的能量值会有差异），但也不收敛。还是说继续增加MAX_SCF的数量（现在inner是35，outer是40），谢谢您。

作者
Author: sobereva 时间: 2024-6-16 22:02

a386514315 发表于 2024-6-16 20:37
感谢社长解答！请问这个问题如何解决呢？按照您的之前的帖子http://bbs.keinsci.com/thread-37196-1-1.ht ...

如果要和AIMD过程的结果保持相同，CUTOFF之类系统性影响能量的设置不要改，按照博文说的尝试其它各种方法

作者
Author: a386514315 时间: 2024-6-16 22:20

sobereva 发表于 2024-6-16 22:02
如果要和AIMD过程的结果保持相同，CUTOFF之类系统性影响能量的设置不要改，按照博文说的尝试其它各种方法

好的，谢谢您，我再看看。

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