PNNL

Abstract

The sensitivity of resulting structures to starting configurations and quench algorithms were characterized using molecular dynamics (MD) simulations. The classical potential model introduced by Damodaran, Rao, and Rao (DRR) Phys. Chem. Glasses 31, 212 (1990) for lead silicate glass was used. Glasses were prepared using five distinct initial configurations and four glass forming algorithms. In previous MD work of bulk lead silicate glasses the ability of this potential model to provide good structural results were established by comparing to experimental results. Here the sensitivity of the results to the simulation methodology and the persistence of clustering with attention to details of molecular structure are determined.