本文对MgO-Al 2O 3-SiO 2 (MAS)微晶玻璃的高温体系的微观结构进行了分子动力学模拟,得出该体系微观体系结构特点以及SiO 2含量对其网络结构的影响。计算结果显示:在MAS高温体系中,Mg 2+的均方位移最大,其次是Al 3+,然后是O 2+,Si 4+的均方位移最小。MAS系微晶玻璃中Si和O之间的作用力最强,Al和O之间的作用力次之,Mg和O之间的结合力最弱。四面体数量由多到少依次是[SiO 4]、[AlO 4],随着SiO 2含量增多[SiO 4]数量逐渐增加,[AlO 4]数量随SiO 2含量的增加呈现先增加然后减少最终趋向平稳的现象。 In this paper, MgO-Al 2O 3-SiO 2 glass-ceramics were calculated by using a molecular dynamics method in this paper, and the microstructure characteristics of the system and the influence of SiO 2 content on its network structure were obtained. It is shown by the calculation results that: In the MAS system at high temperature, the mean square displacement (MSD) of Mg 2+ is the largest, followed by Al 3+, followed by O 2+ and the MSD of Si 4+ is the smallest. The binding force between Si and O in MAS glass-ceramics is the strongest, followed by the force between Al and O, and the binding force between Mg and O is the weakest. The number of tetrahedrons is from more to less, followed by [SiO 4] and [AlO 4]. With the increase of SiO 2 content, the number of [SiO 4] tetrahedrons gradually increased, and the number of [AlO 4] tetrahedrons increased at first and then decreased, and finally reached a state of equilibrium.
张永豪,张金龙,李省伟,崔凯旋,王健健
山东建筑大学,材料科学与工程学院,山东 济南
收稿日期:2019年1月31日;录用日期:2019年2月12日;发布日期:2019年2月19日
本文对MgO-Al2O3-SiO2 (MAS)微晶玻璃的高温体系的微观结构进行了分子动力学模拟,得出该体系微观体系结构特点以及SiO2含量对其网络结构的影响。计算结果显示:在MAS高温体系中,Mg2+的均方位移最大,其次是Al3+,然后是O2+,Si4+的均方位移最小。MAS系微晶玻璃中Si和O之间的作用力最强,Al和O之间的作用力次之,Mg和O之间的结合力最弱。四面体数量由多到少依次是[SiO4]、[AlO4],随着SiO2含量增多[SiO4]数量逐渐增加,[AlO4]数量随SiO2含量的增加呈现先增加然后减少最终趋向平稳的现象。
关键词 :MgO-Al2O3-SiO2系微晶玻璃,分子动力学模拟,微观结构,SiO2
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
微晶玻璃(玻璃陶瓷)多年来一直受到研究人员的广泛关注,其中MgO-Al2O3-SiO2 (MAS)系具有良好绝缘性和抗热震性、高机械强度、低介电常数等优异特性,可广泛应用于雷达天线罩、集成电路基片、硬盘基板和红外辐射材料等高新材料 [
在计算机的发展和多势体 [
在计算过程中建立模型的首要原则是保证体系的电中性,同时需要保持体系中较高的硅含量,这样才可以形成致密的网络结构。在综合考虑计算机的计算能力和能代表MAS系微晶玻璃的粒子数量情况下,本次模拟确定了系统中粒子量保持在1300左右 [
组号 | SiO2 | SiO2/% | 粒子总数 |
---|---|---|---|
1 | 195 | 48.1 | 1215 |
2 | 205 | 49.4 | 1245 |
3 | 215 | 50.6 | 1275 |
4 | 225 | 51.7 | 1305 |
5 | 235 | 52.8 | 1335 |
6 | 245 | 53.8 | 1365 |
7 | 255 | 54.8 | 1395 |
表1. 七组体系初始模拟的具体数据
本次计算的降温机制如图1所示,降温过程为先让计算模型在6000 K下运行20 ps,得到无序和充分熔融的体系,然后用40 ps的时间降温至2000 K,在2000 K弛豫20 ps;然后经过20 ps时间冷却到1800 K,最后在1800 K温度下弛豫20 ps,保证体系处于稳定的高温熔体状态。
图1. 计算中的温度变化图
图2是不同组各原子的均方位移(MSD)变化图。均方位移大小反映了体系中粒子的活跃程度,均方位移量越大,表明粒子在体系中越活跃,其他粒子对该粒子束缚力越小。反之,粒子在体系中位移量越小,表明粒子在体系中的活跃性差,其他粒子对该粒子的束缚能力强。从图中能够得出各组分中不同原子的均方位移量,它们的迁移速率关系均为Mg2+ > Al3+ > O2− > Si4+,这满足微晶玻璃内部粒子的运动特征。
图2. 不同组各粒子的均方位移(MSD)图
出现这种现象的原因是硅离子作为网络形成子参与了网络结构的构建,处于网络中间。镁离子处于游离态不参与网络结构组成,而铝离子作为网络中间体只有在特定条件下才会参与到结构中,有一部分不存在于网络结构之中,氧离子有部分作为桥氧处于网络内部,另外有部分作为非桥氧在网络的边缘活跃性高,所以氧离子的均方位移量比网络形成子硅离子大,但是由于体系中三氧的笼蔽效应导致了体系中氧离子的整体均方位移降低,所以铝离子的均方位移量比氧离子大。
由图3可以得出,高温体系中[SiO4]的数量比[AlO4]的数量多。当SiO2含量增多时,该体系中[SiO4]数量逐渐增多,这表明[SiO4]内部作用力最大,受其他元素的影响小,Si和O具有比Al和O更强的结合能力。而[AlO4]数量随SiO2含量的增加呈现先增加然后减少最终趋向平稳的现象,并且在SiO2含量达到50.6%时,在镁铝硅玻璃陶瓷网络里[AlO4]数量达到最大值161个,出现这种现象的原因是铝离子是以网络中间体形式处于玻璃陶瓷网络结构中,只有在特定的情况下铝离子才可以进入网络结构中成为网络中间体。随着SiO2含量继续不断增加,网络结构变的越来越完整,此时[AlO4]的数量几乎不再变化。
图3. 各组四面体数量示意图
第四组分中,Mg、Al、Si三种粒子数之比为2:4:5,这是堇青石Mg2Al4Si5O18中的原子比,因此是一种较稳定的配方,所以以此组进行分析MAS系微晶玻璃微观熔融结构特征。
由图4可以看出Si-O、Al-O、Mg-O对应的第一峰峰形尖锐程度明显递减,表明Si和O两者的相互作用力最大,组成的团簇结构最牢固,Al和O的相互作用力次之,Mg和O之间作用力最小。并且同种原子之间没有第一尖锐峰,证明同种原子之间没有团聚现象。
图4. 第四组各粒子间的径向函数分布图
本文利用Materials Studio软件对MgO-Al2O3-SiO2微晶玻璃的高温体系的微观结构进行了分子动力学模拟,计算分析了该体系微观结构特点和SiO2含量对其网络结构的影响,得出以下结论:
1) 不同离子在MAS微晶玻璃熔融网络结构中均方位移大小关系为Mg2+ > Al3+ > O2− > Si4+。
2) 在MAS微晶玻璃高温体系中,[SiO4]的数量最多,其次是[AlO4]的数量,随着SiO2含量增多[SiO4]数量逐渐增加,[AlO4]数量随SiO2含量的增加呈现先增加然后减少最终趋向平稳的现象。
3) MAS系微晶玻璃中Si和O之间的作用力最强,Al和O之间的作用力次之,Mg和O之间的结合力最弱。
张永豪,张金龙,李省伟,崔凯旋,王健健. MgO-Al2O3-SiO2系微晶玻璃微观结构的分子动力学模拟 Molecular Dynamics Simulation for the Micro-Structure of MgO-Al2O3-SiO2 Glass-Ceramics[J]. 材料科学, 2019, 09(02): 164-169. https://doi.org/10.12677/MS.2019.92021