A Kinetic Monte Carlo Simulation Study of Synthesis Variables and Diffusion Coefficients in Early Stages of Silicate Oligomerization

Abstract

A kinetic Monte Carlo (kMC) approach combined with density functional theory (DFT) calculations is used to examine the effects of molecular diffusion and synthesis parameters (pH 7–12) as well as initial monomer concentration (0.01–1 mol L–1) for a silicate oligomerization model. To implement this approach, we have adapted the open source kMC SPPARKS software in order to simulate the early stages involved in zeolite formation and more generally kinetically driven reactional systems, using a variety of lattice models (fcc/octahedral/tetrahedral, i.e., fcc/oct/tet). First, these adaptations were validated with kinetically driven reactional systems of the “Lodka” model, providing an excellent match with the analytical solution of the reactive system. The calculations reveal that both the lattice complexity and the diffusion coefficients of species have an impact on the steady state concentration (ssc) of oligomers in solution. Second, the approach is further applied to the early stages of silicate oligomeriz...