A theory for the reliable prediction of the EPR g tensor for paramagnetic defects in solids is presented. It is based on density functional theory and on the gauge including projector augmented wave approach to the calculation of all-electron magnetic response. The method is validated by comparison with existing quantum chemical and experimental data for a selection of diatomic radicals. We then perform the first prediction of EPR g tensors in the solid state and find the results to be in excellent agreement with experiment for the E 0 1 and substitutional phosphorus defect centers in quartz. DOI: 10.1103/PhysRevLett.88.086403 PACS numbers: 71.15. -m, 61.72.Bb, 76.30. -v Electron paramagnetic resonance (EPR), also known as electron spin resonance, is the most powerful spectroscopic technique for the study of paramagnetic defects in solids. Indeed, defect centers are often named directly after their EPR spectra. Applications of EPR extend to any situation where there are unpaired electrons, including the understanding of reactions involving free radicals in both biological and chemical contexts or the study of the structure and spin state of transition metal complexes.EPR spectra of spin 1͞2 centers are made up of two contributions: (i) the hyperfine parameters, which can be computed from the ground state spin density, and have been used to connect theoretical studies of defects to available experimental data [1 -6], and (ii) the g tensor. Only recently have there been attempts to calculate the g tensor in molecules from first principles using density functional theory (DFT) [7,8]. However, these approaches are valid only for finite systems and, thus, are not useful for the calculation of the g tensor for paramagnetic defects in solids, except possibly within a cluster approximation. In the absence of a predictive scheme, experimentally determined g tensors are, of necessity, interpreted in terms of their symmetry alone, leaving any remaining information unexploited. A reliable, first-principles approach to the prediction of g tensors in solids, in combination with structural and energetic calculations, would access this information, and could be used for an unequivocal discrimination between competing microscopic models proposed for defect centers. In this Letter we describe an approach for the calculation of the g tensor in extended systems, using periodic boundary conditions and supercells.In a previous paper [9] we have shown how to compute the all-electron magnetic linear response, in finite and extended systems, using DFT and pseudopotentials. To achieve this we introduced the gauge including projector augmented wave (GIPAW) method, which is an extension of Blöchl's projector augmented wave (PAW) method [10]. In Ref.[9], we used GIPAW to compute the NMR chemical shifts in molecules and solids. Here we apply the GIPAW approach to the first-principles prediction of EPR g tensors for paramagnetic defects in solids. We validate our theory and implementation for diatomic radicals, for which both all-electron...