计算化学公社

标题: 含H体系AIMD时间步长如何选择 [打印本页]

作者
Author: Miles 时间: 2025-10-21 12:56
标题: 含H体系AIMD时间步长如何选择
本帖最后由 Miles 于 2025-10-21 12:59 编辑

各位老师好，对于H2在金属/氧化物体系高温下的氧化反应AIMD模拟，涉及H2-->H2O以及H2和CO氧化生成烃类物质或自由基中间体类的反应，按照 @sobereva 老师的介绍 http://bbs.keinsci.com/forum.php?mod=redirect&goto=findpost&ptid=23337&pid=158130&fromuid=23518，含H的选0.5 fs更为合适，如果不涉及质子转移这样成键断键问题，1 fs也可以接受。

对于上述我提到的具体体系，涉及H的成键断键，如果能选0.5 fs更好，但是从计算时间角度考虑，想请教各位老师时间步长取1 fs的可行性：
1. 直接取1 fs时间步长，对于主要关注涉H化学反应的体系，这样设置后AIMD是否可以准确合理地模拟含H体系的动力学
2. 用D取代H，将H的质量设置为2，这样D-O振动频率变低，周期变长，时间步长取1 fs似乎更具可行性，但增加H质量带来的问题是，对于主要关注涉H化学反应的体系，人为增加H的质量，势必减缓H的动力学，造成模拟得到的反应路径可能与实际有偏差

综上，方案1和方案2哪个更合理呢，抑或是只能选择直接设置时间步长0.5 fs

作者
Author: 丁越 时间: 2025-10-21 13:51
就用0.5fs，你的研究体系涉及大量的含H的反应，并且还是在高温下，小时间步长跑起来更稳妥。
动力学模拟本来就耗时，要想节约模拟时间可以尝试简化模型大小。

作者
Author: Miles 时间: 2025-10-21 20:47

丁越发表于 2025-10-21 13:51
就用0.5fs，你的研究体系涉及大量的含H的反应，并且还是在高温下，小时间步长跑起来更稳妥。
动力学模拟本 ...

谢谢！决定用0.5 fs了
CP2K软件，体系180原子，应该不算大。目前算着有些SCF步难以收敛造成耗时较长（几百s甚至不收敛），一般的SCF在50 s~100 s左右就可以收敛。我继续优化下参数

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