28 28 Tichy, L. and Ticha, H. (2000). Remark on the glass forming ability in GexSe(1‐x) and AsxSe(1‐x) systems. J. Non Cryst. Solids 261: 277–281.
29 29 Gjersing, E.L., Sen, S., and Youngman, R.E. (2010). Mechanistic understanding of the effect of rigidity percolation on structural relaxation in supercooled germanium selenide liquids. Phys. Rev. B 82: 014203–014205.
30 30 Tatsumisago, M., Halfpap, B.L., Green, J.L. et al. (1990). Fragility of Ge‐As – Se glass‐forming liquids in relation to rigidity percolation, and the Kauzmann paradox. Phys. Rev. Lett. 64: 1549–1552.
31 31 Senapati, U. and Varshneya, A.K. (1996). Viscosity of chalcogenide glass forming liquids: an anomaly in the strong and fragile classification. J. Non Cryst. Solids 197: 210–218.
32 32 Naumis, G.G. (2006). Glass transition phenomenology and flexibility: an approach using the energy landscape formalism. J. Non Cryst. Solids 352: 4865–4870.
33 33 Mauro, J.C., Yue, Y., Ellison, A.J. et al. (2009). Viscosity of glass‐forming liquids. Proc. Natl. Acad. Sci. U. S. A. 106: 19780–19784.
34 34 Smedskjaer, M.M., Mauro, J.C., Sen, S., and Yue, Y. (2012). Quantitative design of glassy materials using the temperature dependent constraint theory. Chem. Mater. 22: 5358–5365.
35 35 Zhang, C., Hu, L., Yue, Y., and Mauro, J.C. (2010). Fragile to strong transition in metallic glass forming liquids. J. Chem. Phys. 133: 014508/7.
Note
1 Reviewers:P. Lucas, Materials Science and Engineering, University of Arizona, Tucson, AZ, USAJ. C. Mauro, Materials Science and Engineering, The Pennsylvania State University, PA 16802, USA
2.8 Atomistic Simulations of Glass Structure and Properties
Akira Takada
Research Center, Asahi Glass Co. Ltd., Hazawa‐cho, Yokohama, Japan
1 Introduction
Following Zachariasen's epoch‐making continuous random network (CRN) model [1], the first generation of atomistic models of glass relied on craft construction of ball‐and‐stick representations of atomic structures. Perhaps the most successful model was constructed for SiO2 glass by Bell and Dean [2] who patiently assembled manually rods and balls to compose 188 tetracoordinated units with a total of 614 atoms. The disordered linkage of the SiO4 groups, which were assumed to be rigid, did satisfy Zachariasen's structural rules.
From the angles and distances they measured, Bell and Dean also determined the pair distribution functions (PDF) for Si─Si, Si─O, and O─O, from which they derived a radial distribution function (RDF) that was in good agreement with the experimental data over their full range of definition, i.e. out to a distance of about five times the Si─O bond length. Yielding in particular an average O─Si─O bridging angle of 153°, this model enabled the three‐dimensional atomic configuration of SiO2 glass to be visualized, but it suffered from an obviously large arbitrariness in its handmade combination of building blocks. In addition, the model did not lend itself to estimations of physical properties. An exception was the density, which was calculated for an internal portion of the model system consisting of 72 SiO2 units, to yield a value of 1.99 g/cm3 that proved to be much lower than the actual 2.20 g/cm3.
Atomistic simulations have rendered such mechanical models obsolete as they readily provide not only atomic coordinates but also predict physical properties for glass and melts of any composition under a variety of temperature, pressure, or energy conditions. Besides, the accuracy and versatility of these calculations have been improving steadily, thanks to faster computing processors, more efficient algorithms, and bigger systems investigated. Originally, these simulations were mainly developed to give exact solutions to problems in statistical mechanics which would otherwise have been intractable for complex states of matter such as liquids, glasses, or aggregates.
The first simulations were made with the Monte‐Carlo (MC) method (e.g. [3]), which had been devised in the 1930s to solve general mathematical and statistical problems. It took advantage of the first electronic computers to sample the configurations of the system according to Boltzmann statistics, and weight them evenly when calculating the associated properties of interest. Because of this reliance on Boltzmann statistics, however, only equilibrium states can be investigated in MC simulations. At the cost of much computational complexity, the advantage of molecular dynamics (MD) simulations (e.g. [4]) is to characterize at every moment the state of the system by the positions and momenta of its constituting atoms and, thus, to account for the dynamics of the system whether in equilibrium or not. For SiO2, the first MD simulation performed by Woodcock et al. [5], for instance, dealt with the anomalous properties of the melt and assigned increases in diffusion coefficients to pressure‐induced structural changes. Following this pioneering study on a real material, a great many MD simulations have been carried out on glass/melt systems since that time.
Interatomic potentials are critical ingredients in both MC and MD simulations (e.g. [4]). As a preamble, it is thus appropriate to begin this chapter with a brief review of the manner in which they are determined within the general framework of numerical simulations. The principles of MC and MD simulations will then be presented. Finally, simulations made on amorphous oxide glasses will illustrate the interest and diversity of results that can be obtained with MD simulations and their complementarity with those of experimental studies. Amorphous SiO2, silicates, or phosphates will be mentioned but special interest will be paid to the structure of B2O3 glass that actually represents stringent tests of simulations as it requires accurate descriptions of both short‐ and medium‐range order to be understood. Other important applications of atomistic simulation to transport properties are discussed in Chapter 4.6.
2 Basics of Numerical Simulations
2.1 General Features
Atomistic simulations rely on statistical mechanical models. As such they may be performed within three main kinds of statistical ensembles. The canonical NPT ensemble (constant number of atoms, pressure, and temperature) is chiefly used in heating or quenching cycles. In contrast, the micro‐canonical NVE ensemble (constant number of atoms, volume, and energy) is primarily used when properties are calculated within the precise framework of statistical mechanics. As for the grand canonical μVT ensemble (constant chemical potential, volume, and temperature), it is typically used to investigate chemical equilibrium in systems that can exchange energy and matter with a reservoir.
For an isolated macroscopic system made up of a very large number N of atoms, each having three degrees of freedom, the microscopic state at a given instant is completely specified by the values of 3N coordinates r(i), collectively denoted by r N , and 3N momenta p(i), denoted similarly by p N . The values of the variables r N and p N define a point in a 6N‐dimensional space, called the phase space, symbolized by ΓN. If H(r N , p N ) is the Hamiltonian of the system, the path followed by this point in the phase space is determined by Hamilton's equations:
(1)
(2)
where i = 1,…,N. In principle, 6N coupled equations subjected to 6N initial conditions should be solved to specify the values of all r(i) and p(i) at a given time.
The