原子间径向分布函数问题求助

边锋小奶爸

大佬们，我按照某论文（论文题目：CaO-Al2O3-SiO2微晶玻璃分子动力学模拟的探索研究）中所给条件计算出的原子间径向分布和该文献中差别很大。例如他计算出的Si-O的径向分布函数的首峰在1.6埃左右（其有文献支撑，结果应该没问题），而我是在2.6埃左右。文献中径向分布见图（1），本人计算径向分布见图（2）。弄了三四天了还是不知道原因，请教一下各位大佬，谢谢。具体操作过程如下（均与文献中所给相同）
1）利用amorphous cell tolls建立模型，模型中添加Ca、Al、Si、Fe原子的数量分别为60、120、180、16，再相应补充624个O原子。模型如图1-1所示。设置如下：力场选用pcff，温度为800K，无定形结构数为10，密度为2.7 g/cm3，模型边长为23.67埃，其他默认。
2）利用forcite中calculation的anneal板块模拟降温。模拟中的分段降温过程如图2-1所示，可见过程为先在6000K下驰豫20 ps，然后进行6000K-3000K-2000K-300K的退火，最后在300K下驰豫20 ps达到平衡。设置如下：起始驰豫和降温过程选用NVT系综，结尾驰豫选用NVE系综。积分步程为1fs，
3）利用forcite中analysis板块进行原子间径向分布分析，设置界面如图3-1所示。