FFSiOH: a New Force Field for Silica Polymorphs and Their Hydroxylated Surfaces Based on Periodic B3LYP Calculations

A partial charge shell-ion model potential for silica polymorphs and their hydroxylated surfaces(FFSiOH) was parametrized in a self-consistent way using periodic B3LYP results for bulk R-cristobaliteand the (100) and (001) hydroxylated surfaces. The reliability of the new potentials was checked bycomparing structures, vibrational frequencies and relative phase stabilities of dense bulk silica polymorphs,namely R-quartz, R-cristobalite, R-tridymite, and Stishovite with both experimental and B3LYP data.The FFSiOH was also checked for computing structural and vibrational features of representative all-silica microporous materials, namely edingtonite, chabazite, and faujasite. As a last step, FFSiOH wasadopted to predict OH stretching vibrational frequencies and relative thermodynamic stability of themost common fully hydroxylated surfaces of the dense silica polymorphs, the (100) and (001) facesof all-silica edingtonite, the features of the local Si-defect in chabazite and sodalite known as (SiOH)4hydrogarnet and the geometries of H-bonded silanol groups of an amorphous silica surface. In all casesexcellent agreement resulted between FFSiOH and B3LYP periodic data and experimental data, whenavailable. The new FFSiOH force field opens up the molecular simulation of materials in which thesurface hydroxyl groups play a key role, as is the case for amorphous silica surfaces, all-silica zeoliteexternal surfaces, and the internal walls of mesoporous materials.