17O NMR parameters, both the chemical shifts and the quadrupolar parameters, are calculated for SiO2 polymorphs using density functional theory with the generalized gradient-corrected PBE functional. The gauge including projector augmented wave (GIPAW) method (Pickard, C. J.; Mauri, F. Phys. Rev. B2001, 63, 245101) ensures the reproduction of all electron results while using computationally efficient pseudopotentials. The use of plane-waves permits fully converged calculations to be performed on structures containing 144 atoms in the unit cell, without the need to resort to the cluster approximation. The calculated NMR parameters of cristobalite, quartz, coesite, and faujasite are in excellent agreement with experimental data. This demonstrates that density functional theory is able to reproduce with high accuracy the 17O NMR parameters in SiO2 systems. This precision is used to assign the spectrum of the zeolite ferrierite. The data calculated for SiO2 are used to confirm that no simple correlation between the chemical shift and Cq NMR parameters and Si-O-Si angle exists, emphasizing the importance of predictive theories in this field.

13C, 14N, 15N, 17O, and 35Cl NMR parameters, including chemical shift tensors and quadrupolar tensors for 14N, 17O, and 35Cl, are calculated for the crystalline forms of various amino acids under periodic boundary conditions and complemented by experiment where necessary. The 13C shift tensors and 14N electric field gradient (EFG) tensors are in excellent agreement with experiment. Similarly, static 17O NMR spectra could be precisely simulated using the calculation of the full chemical shift (CS) tensors and their relative orientation with the EFG tensors. This study allows correlations to be found between hydrogen bonding in the crystal structures and the 17O NMR shielding parameters and the 35Cl quadrupolar parameters, respectively. Calculations using the two experimental structures for L-alanine have shown that, while the calculated isotropic chemical shift values of 13C and 15N are relatively insensitive to small differences in the experimental structure, the 17O shift is markedly affected.

This paper presents results of first-principles calculations of nuclear magnetic resonance (NMR) parameters: the chemical shielding tensor and the electric field gradient tensor, of some crystalline and amorphous sodium silicate systems. The calculations have been performed using the recently introduced gauge including projector augmented wave (GIPAW) method, which was especially devised for periodic systems. It provides an attractive alternative to the cluster approximation, used in the previous NMR theoretical studies of silicates systems. Moreover, within the GIPAW formalism, amorphous systems can be efficiently described via a supercell approach as demonstrated in this work. Five reference crystalline compounds of known structure (α-quartz, α-cristobalite SiO2, and the sodium silicates Na2SiO3, α-Na2Si2O5, and β-Na2Si2O5) and two molecular dynamics models of the sodium tetrasilicate glass Na2Si4O9 (NS4) have been studied. The NS4 glass models were generated by a combination of classical and Car-Parrinello molecular dynamics simulations. The good agreement of the simulated 29Si MAS NMR and 17O,23Na 3Q-MAS NMR spectra with the corresponding experimental data demonstrates the accuracy of the GIPAW method. Using these numerically generated data, we have also been able to gain insight into the correlation between NMR parameters and local structural features.

We apply density functional theory (DFT) to the calculation of the (17)O NMR parameters in Ca and Mg oxides and aluminosilicates. We study the accuracy of the Perdew, Burke, and Ernzerhof (PBE) generalized-gradient approximation to DFT in the description of these systems and the origin of the experimentally observed large dependence of the (17)O chemical shift on the alkaline earth ion. We find that (i) the partially covalent nature of the Ca-O bond has a huge impact on the O chemical shifts. The Ca-O covalence alone explains why in Ca oxides and aluminosilicates the (17)O chemical shifts are much more deshielded than those of the corresponding Mg compounds. (ii) The Ca-O covalence is overestimated by the PBE functional. Thus PBE-DFT is not able to reproduce the measured (17)O NMR parameters in Ca oxide and Ca aluminosilicates. (iii) It is possible to correct for the PBE-DFT deficiency in a simple and transferable way and to predict very accurate (17)O NMR parameters. Such accuracy allows us to assign the (17)O NMR spectra of two important model systems: the grossite aluminate (CaAl(4)O(7)) and the wollastonite (CaSiO(3)) silicate.