在热解、积碳等反应模拟中,经常会出现产物分布与实验不符合的情况,实际上这对反应力场来说,得到与实验产物分布完全一致才是非常困难的。根本原因可以从2方面考虑(以下内容仅供参考,尤其第二点缺乏足够证据支持):
ReaxFF的力场并不能完全与实际一致,实际上ReaxFF产生的活化能往往与实际情况相差相当大,甚至大小顺序也不一致,这直接导致不同方向的反应速率常数不同,从而导致产物分布不同
ReaxFF的动力学过程与实际实验相比,扩散随机性较差,因此经常出现大量自由基无法中和(这一点可以考虑充分延长模拟时间),这与实验直观相比非常不同,实验中很难有大量残留自由基。由于自由基扩撒不均匀,因此自由基的结合,可能就会产生与实验不同的结合方式,从而产生与实验不符合的产物。
第一个方面,或许可以通过优化力场得到解决,但这个方面工作量较大,对研究者的理论素养要求非常高。第二个方面,用户可以尝试加入fbMC,增大随机性;或者使用纳米反应器,充分搅拌,但这两种方法的有效性有待验证,仅供参考。