计算化学公社

 找回密码 Forget password
 注册 Register
Views: 247|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[综合讨论] 含H体系AIMD时间步长如何选择

[复制链接 Copy URL]

6

帖子

0

威望

280

eV
积分
286

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 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

494

帖子

11

威望

4231

eV
积分
4945

Level 6 (一方通行)

2#
发表于 Post on 2025-10-21 13:51:06 | 只看该作者 Only view this author
就用0.5fs,你的研究体系涉及大量的含H的反应,并且还是在高温下,小时间步长跑起来更稳妥。
动力学模拟本来就耗时,要想节约模拟时间可以尝试简化模型大小。
自由发挥,野蛮生长

6

帖子

0

威望

280

eV
积分
286

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2025-10-21 20:47:25 | 只看该作者 Only view this author
丁越 发表于 2025-10-21 13:51
就用0.5fs,你的研究体系涉及大量的含H的反应,并且还是在高温下,小时间步长跑起来更稳妥。
动力学模拟本 ...

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-24 05:40 , Processed in 0.265886 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list