Читать книгу Encyclopedia of Glass Science, Technology, History, and Culture - Группа авторов - Страница 308

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 SiO₂ 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 SiO₄ 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 SiO₂ 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 SiO₂ units, to yield a value of 1.99 g/cm³ that proved to be much lower than the actual 2.20 g/cm³.

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 SiO₂, 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 SiO₂, silicates, or phosphates will be mentioned but special interest will be paid to the structure of B₂O₃ 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.