Assessment of periodic and cluster-in-vacuo models for first principles calculation of EPR parameters of paramagnetic defects in crystals: Rh2+ defects in NaCl as case study.

In order to find a reliable and efficient calculation scheme for electron paramagnetic resonance (EPR) spectroscopic parameters for transition metal complexes in ionic solids from first principles, periodic and finite cluster-in-vacuo density functional theory (DFT) simulations are performed for g tensors, ligand hyperfine tensors (A), and quadrupole tensors (Q) for Rh(2+)-related centers in NaCl. EPR experiments on NaCl:Rh single crystals identified three Rh(2+) monomer centers, only differing in the number of charge compensating vacancies in their local environment, and one dimer center. Periodic and cluster calculations, both based on periodically optimized structures, are able to reproduce experimentally observed trends in the ligand A and Q tensors and render very satisfactory numerical agreement with experiment. Taking also computation time into account as a criterion, a full periodic approach emerges as most appropriate for these parameters.The g tensor calculations, on the other hand, prove to be insufficiently accurate for model assessment. The calculations also reveal parameters of the complexes which are not directly accessible through experiments, in particular related to their geometry.