Developing atomistic glass models using potential-free Monte Carlo method: From simple to complex structures

We propose here a method to generate random networked amorphous structure using only readily available short-range properties like bond lengths, bond angles and connectivity of the constituents. This method is a variant of Monte-Carlo (MC) method wherein the basic constituents of an amorphous network i.e. rigid polyhedral units are connected randomly obeying certain steric constraints. The algorithm is designed to reproduce the medium-range order universally observed in glasses. The method somewhat resembles the reverse MC (RMC) method where a random move of an atom inside a box is accepted or rejected depending upon whether it decreases or increases the deviation from the experimentally observed features. However unlike RMC, this method does not demand large experimental sets of scattering data which in most cases is a priori not available for glasses. It rather relies on the stochasticity of MC method to produce glassy structures. The algorithm is first validated against SiO2 glass structure by comparing with the available structures from other methods and experimental data. The method is then extended for developing more complex Iron Phosphate Glass (IPG) structures and a comparison with existing models of IPG developed using quench-from-melt scheme implemented in classical Molecular Dynamics (MD) reveals that the method is extensible to complex glasses also. This study addresses the often-neglected issue of non-availability of correct starting structures in simulating glasses using MD or Density Functional Theory (DFT).Comment: 32 pages, 11 figures, journal articl