gaussian的ADMP模拟中加或不加FullScf的区别及与BOMD的区别

diaolanxinyu

下面给出的结果均在同一结构下的绝热的动力学模拟。
gaussian的ADMP模拟中加或不加FullScf的区别，见图1和图2，其中图1和图2的左边一列是未加FullScf的ADMP模拟，右边一列是加了FullScf的ADMP模拟




gaussian的BOMD的模拟结果如图3
1. 从图1和图2中可看出，图1的Electronic Kinetic energy差别最大，原因是加了FullScf，在每一点用 SCF 收敛的结果进行动力学计算，我的问题是这个SCF收敛是对什么进行迭代计算的？
总能量在整个模拟过程中没有完全守恒，但偏差都较小，这又是什么原因造成的呢？
2.从图3中可以看出，BOMD模拟结果没有Electronic Kinetic energy这一项，而且总能量在整个模拟过程中是完全守恒的，这符合绝热的MD模拟，那BOMD模拟中每一步又是对什么进行迭代计算的呢？
谢谢大家！

xiewenli

你好，我也对这两个关键字的用法和结果不太懂，可以解释一下吗？谢谢。

清微

sobereva 发表于 2018-1-23 16:53
和算法以及程序实现有关系。BOMD完成之后末尾有汇总，会给出每一步的时间，虽然不均匀，但作图没问题

BOMD每一步的时间有什么含义呢？是按着预设的步幅跑完的时间吗？

sobereva

清微 发表于 2018-1-23 16:25
请问老师，ADMP和BOMD中的参数StepSize含义大相径庭，前者是步长，单位是fs,后者是步幅，单位是Bohr，究竟 ...

和算法以及程序实现有关系。BOMD完成之后末尾有汇总，会给出每一步的时间，虽然不均匀，但作图没问题

清微

请问老师，ADMP和BOMD中的参数StepSize含义大相径庭，前者是步长，单位是fs,后者是步幅，单位是Bohr，究竟是什么原因呢？而BOMD中不能设定步长，那模拟结果怎么对时间作图？

sobereva

BOMD每一步就跟单点+force关键词一样，即每走一步都做SCF直到波函数收敛，然后计算受力。这样也没有电子在动力学过程中的动能的概念。

用了FullSCF则ADMP本质上也是每步迭代到波函数收敛了，这也就和BOMD一样没有动力学中的电子动能概念了。有动能的话就必然是没相对于当前几何结构收敛的波函数。

ADMP总能量那点波动很正常，毕竟计算时候步长是有限的，没法做到能量的绝对守恒。看上去BOMD总能量没变，是条直线，但注意坐标轴范围，远比ADMP的图大得多，有波动从图上也看不出来。