Water dissolution in albite melts:: Constraints from ab initio NMR calculations

Hartree-Fock and B3LYP NMR calculations were performed at the 6-311+G(2df,p) level on cluster models representing albite glasses using B3LYP/6 to 31G* optimized geometries. Calculation results on several well-known crystalline materials, such as low albite and KHSi₂O_5, were used to check the accuracy of the calculation methods.Calculated ²⁹Si-NMR results on clusters that model protonation of Al-O-Si linkages and the replacement of Na⁺ by H⁺ indicate a major increase in Si-O(H) bond length and a 5 ppm difference in δ_iso for ²⁹Si compared to that for anhydrous albite glass. The calculated δ_iso of ²⁷Al in such linkages agrees with the experimental data, but shows an increase in C_q that cannot be fully diminished by H-bonding to additional water molecules. This protonation model is consistent with both experimental ¹⁷O NMR data and the major peak of ¹H-NMR spectra. It cannot readily explain the existence of the small peak in the experimental ¹H spectra around 1.5 ppm. Production of the depolymerized units Al Q³]-O-H upon the dissolution of water is not consistent with ²⁷Al, ¹H, or ¹⁷O NMR experimental results. Production of Si Q³]-O-H is consistent with all of the experimental ¹⁷O and ¹H-NMR data; such units can produce both the major peak at 3.5 ppm and the small peak at 1.5 ppm in ¹H spectra, either with or without hydrogen bonding. This species, however, cannot produce the main features of ²⁹Si spectra.It is concluded that although neither protonation nor the production of Si Q³]-O-H alone is consistent with the available experimental data, the combination of these two processes is consistent with available experimental NMR data.