计算化学公社

标题: 求助MS中GULP模块的ReaxFF力场去模拟煤的低温氧化计算RDF没有峰怎么办 [打印本页]

作者
Author: wanglinchuan 时间: 2026-1-10 13:52
标题: 求助MS中GULP模块的ReaxFF力场去模拟煤的低温氧化计算RDF没有峰怎么办
用MS的GULP模块的ReaxFF力场进行分子动力学模拟煤的低温氧化，以下是我的具体操作，具体操作步骤如下，请帮我看一下这个模拟思路有问题没？还有一些具体的问题，求帮忙解答，谢谢。
1.导入实验数据构建的煤分子结构，进行COMPASS力场来进行一个简单的几何优化。
(, 下载次数 Times of downloads: 0)

2.使用AC模块的COMPASS力场进行盒子的构建，图中是一些参数设置。
(, 下载次数 Times of downloads: 1)

3.对构建的盒子取10帧中能量最低的一个，保存为.xsd文件在GULP模块中进行ReaxFF力场的几何优化。
(, 下载次数 Times of downloads: 1)

4.在GULP中进行Anneal退火操作，以下是设置的参数。
(, 下载次数 Times of downloads: 1)
5.使用GULP模块进行MD模拟，以下是设置的参数。
(, 下载次数 Times of downloads: 1)
6.取出.xtd文件标记氧气分子和煤分子中的羰基，进行下一步的RDF分析。
(, 下载次数 Times of downloads: 0)
7.使用脚本文件进行RDF扩散系数的计算，但是对于这个分析结果，总是没有所谓第一配位峰和第二配位峰，这个问题应该如何解决？
(, 下载次数 Times of downloads: 1)
以下是脚本文件的参数设置，还有跑出来的RDF分析。脚本参数中截断距离为RDF_cut = 13，Interval = 0.05。
(, 下载次数 Times of downloads: 1)
请各位老师帮忙看看我的这个操作有没有问题，这个峰的问题应该如何解决？

作者
Author: wanglinchuan 时间: 2026-1-13 18:08
是不是对于RDF分析的对象选择错误呀

欢迎光临计算化学公社 (http://bbs.keinsci.com/)