MS分析Si-O的原子间径向分布问题求助！

边锋小奶爸

大佬们，快救救我，谢谢大佬们。

我用MS分析Si-O的原子间径向分布函数的首峰总是在2.6埃左右，我看很多文献都是1.6埃，不知道咋回事。

具体是首先用MS中的amorphous cell tools建模导入成一定比例的Ca、Al、Si原子；然后用forcite中dynamics驰豫一定时间，再用analysis中的RDF进行分析。

设置基本上均参考文献，具体如下：
amorphous cell tolls设置：力场选用pcff，温度为800K，无定形结构数为10，密度为2.7 g/cm3，模型边长为23.67埃，其他默认。
forcite中dynamics设置：系综选用NVE，力场为pcff，温度为300K，总时间20ps，其他默认。
analysis设置如图1所示。

文献有：
[1]Molecular dynamics calculations for Mg3Al2Si3O12 (pyrope) and Ca3Al2Si3O12 (grossular) glass structures, Masayuki Okuno aY and Katsuyuki Kawamura
[2]CaO-Al2O3-SiO2-Fe2O3体系微晶玻璃分子动力学模拟的探索研究, 王明国
[3]CaO-SiO2-Al2O3-Na2O体系熔体微观结构与性质的分子动力学模拟, 肖成

AaronO_o

检查一下力场类型和电荷分配？
看您先前的模型描述，似乎原子间没有价键，那么我个人会考虑非价键的库伦作用和vdW作用计算是否存在问题。比如Si和O可能没有指派合适的电荷、中性的原子间没有库仑作用带来的引力，导致没有预期的相互作用？
另外这个问题应该放在“分子模拟”版块更好吧、也跟标签中的CASTEP和DMol3并没有关系、、、

边锋小奶爸

AaronO_o 发表于 2021-9-12 14:17
检查一下力场类型和电荷分配？
看您先前的模型描述，似乎原子间没有价键，那么我个人会考虑非价键的库伦作 ...

好的，多谢大佬。我再试试（主要在分子模拟那边发了没有人回复，这边人多点我试试）

sobereva

帖子标题里别随便带叹号，仔细看http://bbs.keinsci.com/thread-9348-1-1.html，给你去掉了，以后注意

若无绝对必要，不要在标题里用叹号，否则显得有攻击性或者显得很浮躁。

不是第一性原理计算相关的问题就不要发到第一性原理板块，给你移动了。
“在分子模拟那边发了没有人回复，这边人多点我试试”这绝对不是借口，你的做法纯粹是扰乱论坛秩序！别把论坛秩序不当回事。下次再这样直接删帖扣分处理

边锋小奶爸

sobereva 发表于 2021-9-12 20:34
帖子标题里别随便带叹号，仔细看http://bbs.keinsci.com/thread-9348-1-1.html，给你去掉了，以后注意

...

好的，不好意思