在使用ReaxFF进行正式的分子动力学模拟之前,一般会对体系进行弛豫(relax)或者能量最小化(或者叫做结构优化)。
弛豫或者能量最小化的目的,是为了得到一个平衡、均匀的反应物,防止初始结构不合理,而导致分子动力学模拟结果不可靠。那么什么时候使用能量最小化,什么时候使用弛豫呢?实际上所有的情况都可以使用弛豫代替能量最小化。能量最小化对MD而言意义不大,因为最终因为MD过程,都会偏离平衡位置。
对于体系中全部是强的化学键,例如硅单晶,里面原子和原子的相互作用都很强,那么能量最小化是可以很容易收敛的。但是如果体系中存在范德华力、氢键,那么这样的体系进行能量最小化很难收敛(而且收敛到能量最小值,本身对MD也没有多大意义),分子个数越多越难收敛,这种情况下就只能用弛豫的方式处理初始结构。
因此在MD模拟之前,建议都使用弛豫即可。