Abstract
This paper investigates the applicability of the bond-order, variable-charge (BOVC) force fields of the Charge-Optimized Many-Body (COMB10), ReaxFFSiOH2O, and ReaxFFSiOGSI, for molecular dynamics (MD) simulations of crystalline SiO2. The calculated lattice constants and densities of the four SiO2 polymorphs, quartz, cristobalite, coesite, and stishovite, are compared to experimental values. Additionally, the calculated pair distribution and bond-angle distribution functions and the α–β transition for quartz, the most stable low-energy polymorph, are compared to experimental results. The simulations with the COMB10 force field accurately predict the properties of the SiO2 polymorphs, except the α-cristobalite, and the quartz α–β transition. The results with ReaxFFSiOH2O and ReaxFFSiOGSI accurately predict the properties of the SiO2 polymorphs, except the stishovite, but those with ReaxFFSiOH2O inaccurately predict the quartz α–β transition.
Published Version
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a call