Mechanisms and energetics of calcium aluminosilicate glass dissolution through ab initio molecular dynamics-metadynamics simulations

The dissolution of silicate glasses has implications in diverse fields ranging from the immobilization of radioactive waste to the development of sustainable alternatives to Portland cement. Here, we used ab initio molecular dynamics simulations biased with well-tempered metadynamics to study Si-O-T bridge dissociation in calcium aluminosilicate glasses, crucial for understanding their dissolution. In a departure from the conventional Michalske-Freiman model, our findings reveal a nucleophilic substitution reaction mechanism characterized by a short-lived, 5-fold coordinated Si intermediate or transition state, depending on the Si bridge coordination, with a near-trigonal bipyramidal geometry. We find that the reorganization required for reaching this state causes the activation energy barriers to be dependent on the Si bridge coordination, with Si Q3 species serving as the rate-limiting step in the dissolution reaction. Our findings not only challenge long-standing theoretical models but also pave the way for more accurate and comprehensive frameworks for understanding the dissolution of silicate glasses in various applications.