AUTHOR=Gong Kai , Yang Kengran , White Claire E. TITLE=Density functional modeling of the binding energies between aluminosilicate oligomers and different metal cations JOURNAL=Frontiers in Materials VOLUME=10 YEAR=2023 URL=https://www.frontiersin.org/journals/materials/articles/10.3389/fmats.2023.1089216 DOI=10.3389/fmats.2023.1089216 ISSN=2296-8016 ABSTRACT=

Interactions between negatively charged aluminosilicate species and positively charged metal cations are critical to many important engineering processes and applications, including sustainable cements and aluminosilicate glasses. In an effort to probe these interactions, here we have calculated the pair-wise interaction energies (i.e., binding energies) between aluminosilicate dimer/trimer and 17 different metal cations Mⁿ⁺ (Mⁿ⁺ = Li⁺, Na⁺, K⁺, Cu⁺, Cu²⁺, Co²⁺, Zn²⁺, Ni²⁺, Mg²⁺, Ca²⁺, Ti²⁺, Fe²⁺, Fe³⁺, Co³⁺, Cr³⁺, Ti⁴⁺ and Cr⁶⁺) using a density functional theory (DFT) approach. Analysis of the DFT-optimized structural representations for the clusters (dimer/trimer + Mⁿ⁺) shows that their structural attributes (e.g., interatomic distances) are generally consistent with literature observations on aluminosilicate glasses. The DFT-derived binding energies are seen to vary considerably depending on the type of cations (i.e., charge and ionic radii) and aluminosilicate species (i.e., dimer or trimer). A survey of the literature reveals that the difference in the calculated binding energies between different Mⁿ⁺ can be used to explain many literature observations associated with the impact of metal cations on materials properties (e.g., glass corrosion, mineral dissolution, and ionic transport). Analysis of all the DFT-derived binding energies reveals that the correlation between these energy values and the ionic potential and field strength of the metal cations are well captured by 2nd order polynomial functions (R² values of 0.99–1.00 are achieved for regressions). Given that the ionic potential and field strength of a given metal cation can be readily estimated using well-tabulated ionic radii available in the literature, these simple polynomial functions would enable rapid estimation of the binding energies of a much wider range of cations with the aluminosilicate dimer/trimer, providing guidance on the design and optimization of sustainable cements and aluminosilicate glasses and their associated applications. Finally, the limitations associated with using these simple model systems to model complex interactions are also discussed.