An ab initio calculation of ¹⁷O and ²⁹Si NMR parameters for SiO₂ polymorphs

Ab initio molecular orbital calculations (Hartree-Fock, HF and density functional theories, DFTs) have been carried out for SiO₂ polymorphs coesite, low cristobalite, and α-quartz, in order to investigate the reliability of this method for predicting ²⁹Si and ¹⁷O nuclear magnetic resonance (NMR) properties of silicates. Oxygen- and silicon-centered clusters consisting of one (1T) to three tetrahedral (3T) shells (one to four atomic shells), taken from real crystal structure, have been investigated. It is found that for reasonable predication of both the ²⁹Si and ¹⁷O chemical shifts (δ_i^Si and δ_i^O), the minimum cluster is one that gives the correct second neighbors to the nucleus of interest. Both the δ_i^Si and δ_i^O have reached convergence with respect to cluster size at the OH-terminated two tetrahedral (2T) shell (three atomic shells around Si and four atomic shells around O) model. At convergence, the calculated δ_i^Si values agree well (within ±1 ppm) with experimental data. The calculated ¹⁷O electric field gradient (EFG)-related parameters also agree with experimental data within experimental uncertainties. The calculation also reproduces small differences in δ_i^O for O sites with similar tetrahedral connectivities, but shows deviations up to about 10 ppm in relative difference for O sites with different tetrahedral connectivities. The poor performance for the latter is mainly due to the approximations of the HF method. Our study thus suggests that the ab initio calculation method is a reliable mean for predicting ²⁹Si and ¹⁷O NMR parameters for silicates. Such an approach should find application not only to well-ordered crystalline phases, but also to disordered materials, by combining with other techniques, such as the molecular dynamics simulation method.