Molecular dynamics simulation of temperature-induced structural changes in cristobalite, coesite and amorphous silica

• Published in: Journal of non-crystalline solids Vol. 354; no. 2-9; pp. 181 - 187
• Main Authors: Takada, Akira, Richet, P., Catlow, C.R.A., Price, G.D.
• Format: Journal Article Conference Proceeding
• Language: English
• Published: Amsterdam Elsevier B.V 15.01.2008; Elsevier
• Subjects: Condensed matter: structure, mechanical and thermal properties; Disordered solids; Exact sciences and technology; Molecular dynamics; Other amorphous solids; Physics; Short-range order; Silica; Structure of solids and liquids; crystallography; 62.50.+p; 02.70.Ns; 64.70.−p; Molecular dynamics; Short-range order; 61.43.Fs; Silica; Melting; Phase transformations; Molecular dynamics method; Amorphous state; Bond lengths; Molecular rotation; Glass transition; Molecular dynamics; Silica; Short-range order; Temperature effects; 61.43.Fs; 62.50.+p; 64.70.-p; 02.70.Ns; Amorphous state structure
• ISSN: 0022-3093; 1873-4812
• DOI: 10.1016/j.jnoncrysol.2007.07.062

Structural changes in cristobalite, coesite, and amorphous silica have been investigated from 100 to 3000K by molecular dynamics simulations. These have been made with ‘soft potentials’ that yield melting and glass transition temperatures in fair agreement with experimental data. The simulated density of cristobalite increases abruptly at the α–β transition, that of coesite increases gradually, whereas that of amorphous silica goes through a density maximum around 2000K. In the three phases, the complex temperature-induced structural changes have been analyzed in terms of variations of the abundances of four ‘structons’ which represent different kinds of Si–O bond lengths and torsional angles between two SiO4 units. These changes in cristobalite, coesite and amorphous silica can then be rationalized in terms of loosening of network and of dynamical rotation of SiO4 units. The relation between structural changes and ring statistics is also discussed.