University of North Texas 
Witold Brostow a, Mieczyslaw Chybicki b, Robert Laskowski b, and Jaroslaw Rybicki b
a Department of Materials Science, University of North Texas, Denton, TX 76203-5310,
USA
b Department of Solid State Physics, Technical University of Gdansk, Narutowicza 11/12,
80-952 Gdansk, Poland
ABSTRACT
Voronoi and Delaunay tessillations are applied to pattern recognition of atomic environment and to investigation of the nonlocal order in molecular-dynamics (MD)-simulated materials. The method is applicable also to materials generated using other computer techniques such as Monte Carlo. The pattern recognition is based on an analysis of the shapes of the Voronoi polyhedron (VP). A procedure for contraction of short edges and small faces of the polyhedron is presented. It involves contraction to vertices of all edges shorter than a certain fraction x of the average edge length, with concomitant contraction of the associated faces. Thus, effects of fluctuations are eliminated, providing "true" values of the geometric coordination numbers f, both local averaged over the material. Nonlocal order analysis involves geometric relations between Delaunay simplexes. The methods proposed are used to analyze the structure of MD-simulated solid lead [J. Rybicki, W. Alda, S. Feliziani, and W. Sandowski, in Proceedings of the Conference on Intermolecular Interactions in Matter, edited by K. Sangwal, E. Jartych, and J.M. Olchowik (Technical University of Lublin, Lublin, 1995), p. 57: J. Rybicki, R. Laskowski, and S. Feliziani, Comput. Phys. Commun. 97, 185 (1997)] and germianium dioxide [T. Nanba, T. Miyaji, T. Takada, A. Osaka, Y. Minura, and I. Yosui, J. Non-Cryst. Solids 177, 131 (1994)]. For Pb the contraction results are independent of x. For the open structure of GeO2 there is an x dependence of the contracted structure, so that using several values of x is preferable. In addition to removing effects of thermal perturbation, in open structures the procedure also cleans the resulting VP from faces contributed by the second neighbors. The analysis can be combined with that in terms of the radial distribution g ( R ), making possible comparison of geometric coordination numbers with structural ones [W. Brostow, Chem. Phys. Lett. 49, 285 (1977)].
* Physical Review 1998, 57, 13448