求助MS中GULP模块的ReaxFF力场去模拟煤的低温氧化计算RDF没有峰怎么办

wanglinchuan

用MS的GULP模块的ReaxFF力场进行分子动力学模拟煤的低温氧化，以下是我的具体操作，具体操作步骤如下，请帮我看一下这个模拟思路有问题没？还有一些具体的问题，求帮忙解答，谢谢。
1.导入实验数据构建的煤分子结构，进行COMPASS力场来进行一个简单的几何优化。



2.使用AC模块的COMPASS力场进行盒子的构建，图中是一些参数设置。


3.对构建的盒子取10帧中能量最低的一个，保存为.xsd文件在GULP模块中进行ReaxFF力场的几何优化。



4.在GULP中进行Anneal退火操作，以下是设置的参数。

5.使用GULP模块进行MD模拟，以下是设置的参数。

6.取出.xtd文件标记氧气分子和煤分子中的羰基，进行下一步的RDF分析。

7.使用脚本文件进行RDF扩散系数的计算，但是对于这个分析结果，总是没有所谓第一配位峰和第二配位峰，这个问题应该如何解决？

以下是脚本文件的参数设置，还有跑出来的RDF分析。脚本参数中截断距离为RDF_cut  =  13，Interval = 0.05。

请各位老师帮忙看看我的这个操作有没有问题，这个峰的问题应该如何解决？

wanglinchuan

是不是对于RDF分析的对象选择错误呀