Use of site symmetry in supercell models of defective crystals: Polarons in CeO2

The authors thank R. Merkle and G. W. Watson for stimulating discussions. E. K. also acknowledges partial financial support from the Russian Science Foundation for the study of charged defects under the project 14-43-00052. A. C. also acknowledges financial support from the University of Latvia Foundation (Arnis Riekstins's "MikroTik" donation). E. K. and D. G. express their gratitude to the High Performance Computer Centre in Stuttgart (HLRS, project DEFTD 12939) for the provided computer facilities whereas R. A. E. thanks the St. Petersburg State University Computer Center for assistance in high-performance calculations.In supercell calculations of defective crystals, it is common to place a point defect or vacancy in the atomic position with the highest possible point symmetry. Then, the initial atomic structure is often arbitrary distorted before its optimization, which searches for the total energy minimum. In this paper, we suggest an alternative approach to the application of supercell models and show that it is necessary to preliminarily analyze the site symmetry of the split Wyckoff positions of the perfect crystal supercell atoms (which will be substituted or removed in defective crystals) and then perform supercell calculations with point defects for different possible site symmetries, to find the energetically most favorable defect configuration, which does not necessarily correspond to the highest site symmetry. Using CeO2 as an example, it is demonstrated that this use of the site symmetry of the removed oxygen atoms in the supercells with vacancies allows us to obtain all the possible atomic and magnetic polaron configurations, and predict which vacancy positions correspond to the lowest formation energies associated with small polarons. We give a simple symmetry based explanation for the existence of controversies in the literature on the nature of the oxygen vacancies in CeO2. In particular, the experimentally observed small polaron formation could arise for oxygen vacancies with the lowest Cs site symmetry, which exist in 3 x 3 x 3 and larger supercells. The results of first principles calculations using a linear combination of atomic orbitals and hybrid exchange-correlation functionals are compared with those from previous studies, obtained using a widely used DFT+U approach.Russian Science Foundation 14-43-00052; Saint Petersburg State University; University of Latvia Foundation (Arnis Riekstins's "MikroTik" donation); Institute of Solid State Physics, University of Latvia as the Center of Excellence has received funding from the European Union’s Horizon 2020 Framework Programme H2020-WIDESPREAD-01-2016-2017-TeamingPhase2 under grant agreement No. 739508, project CAMART