|
|
本帖最后由 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
|
|