赵 谦 1,祖 群 1,齐 亮 2,胡永杰 2,孙雪坤3,陈 阳 1
(1. 中材科技股份有限公司,中国南京 210012;2. 美国密歇根大学,美国密歇根州安娜堡 48109-236; 3. 美国大陆地科技有限公司,美国印第安纳波利斯 46033)
摘 要:采用 3 种原子间作用势函数,对不同含量的二元钠硅酸盐玻璃[(xNa2O(100–x)SiO2)]结构与性能进行分子动力学(MD) 模拟。结果表明,相比 Born-Mayer-Huggins (BMH)势和 Morse 势,Buckingham 势函数计算的密度和模量结果与实验数据更为接近。基于 Buckingham 势的 MD 计算结果表明:随 Na2O 含量增加,玻璃密度增大,剪切和杨氏模量降低;Na2O含量小于 20%时,体积弹性模量呈降低趋势,当 Na2O含量(摩尔分数)达到 20%后,随 Na2O 含量的增加而增大,变化趋势与实测数据相符。径向分布函数分析玻璃结构表明,随Na2O含量增加,玻璃短程有序性先增加后降低,与体积弹性模量变化呈负向对应关系,剪切模量和杨氏模量随玻璃结构中非桥氧的离子键增加而降低。