石英玻璃分子动力学模拟中的原子电荷转移与系综选择

摘要

介绍了SiO2体系分子动力学模拟中的Si、O原子电荷转移问题;采用Morse势函数研究了原子电荷转移对石英玻璃模拟的影响,发现原子电荷转移在影响模型密度的同时,还直接影响着原子的最近邻距离.NPT和NVT系综下的模拟结果对比显示,系综对模型中原子最近邻情况影响不大,但在NVT系综下模拟结果表明实际玻璃中存在的较大的空隙结构,找到了以往模拟中密度结果偏高的原因,提出了一种较好的石英玻璃分子动力学建模的方法.该方法不但解决了在调整电荷时维持原子最近邻距离与保证模型密度之间的矛盾,而且可以很好地描述石英玻璃在远程结构上密度不均、存在较大空隙的无序结构.此外,原子自扩散系数的计算结果展示了空隙结构在石英玻璃扩散性质研究中的作用.