物理化学学报 >> 2010, Vol. 26 >> Issue (06): 1651-1656.doi: 10.3866/PKU.WHXB20100615

量子化学及计算化学 上一篇    下一篇

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

丁元法, 张跃, 张大海, 李仲平   

  1. 北京航空航天大学材料科学与工程学院, 空天材料与服役教育部重点实验室, 北京 100191; 航天材料及工艺研究所, 先进功能复合材料技术国防科技重点实验室, 北京 100076
  • 收稿日期:2010-01-11 修回日期:2010-03-10 发布日期:2010-05-28
  • 通讯作者: 张跃 E-mail:zhangy@buaa.edu.cn

Atomic Charge Transfer and Ensemble Selection inMolecular Dynamics Simulation of Vitreous Silica

DING Yuan-Fa, ZHANG Yue, ZHANG Da-Hai, LI Zhong-Ping   

  1. Key Laboratory of Aerospace Materials and Performance (Ministry of Education), School of Materials Science and Engineering, Beihang University, Beijing 100191, P. R. China; National Key Laboratory of Advanced Functional Composite Materials Technology, Aerospace Research Institute of Materials and Processing Technology, Beijing 100076, P. R. China
  • Received:2010-01-11 Revised:2010-03-10 Published:2010-05-28
  • Contact: ZHANG Yue E-mail:zhangy@buaa.edu.cn

摘要:

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

关键词: 分子动力学, 石英玻璃, 电荷转移, 系综, 空隙

Abstract:

In this paper, we introduced charge transfer of Si and O atoms for molecular dynamics simulation of SiO2. The effect of atomic charge transfer on vitreous silica simulation was investigated by Morse potential. Results show that atomic charge transfer not only affects the density of vitreous silica models, but also affects the interatomic nearest distance directly. The result comparison between the NPT and NVT ensembles indicates that the interatomic nearest-distances in the two ensembles are similar with the same atomic charge transfer but larger voids are found in the NVT model, which can account for the void structure in real glass. We clarify why simulations overestimate densities compared to experimental data. We also describe a better method of molecular dynamics modeling for vitreous silica. This method resolves the conflict between maintaining the interatomic nearest distance and keeping the model density as the atomic charge transfer changes. It also describes the long-range disordered structure of density non-uniformity and larger voids well. Furthermore, the result of atomic self-diffusion coefficient shows the significance of a larger void structure in the study of diffusion properties in vitreous silica.

Key words: Molecular dynamics, Vitreous silica, Charge transfer, Ensemble, Void

MSC2000: 

  • O641