title: Use of site symmetry in supercell models of defective crystals: polarons in CeO2.
abstract: 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 CeO2as 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 Cssite symmetry, which exist in 3 × 3 × 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.
authors.affiliation: Department of Chemistry, St. Petersburg State University, St. Petersburg, Russia.
authors.forename: R A
authors.initials: RA
authors.lastname: Evarestov
authors.affiliation: Department of Theoretical Physics and Computer Modelling, Institute of Solid State Physics, University of Latvia, Riga, Latvia. gryaznov@cfi.lu.lv.
authors.forename: D
authors.identifier: http://orcid.org/0000-0002-2894-4471
authors.initials: D
authors.lastname: Gryaznov
authors.affiliation: Department of Physical Chemistry of Solids, Max Planck Institute for Solid State Research, Stuttgart, Germany.
authors.forename: M
authors.initials: M
authors.lastname: Arrigoni
authors.affiliation: Department of Theoretical Physics and Computer Modelling, Institute of Solid State Physics, University of Latvia, Riga, Latvia. gryaznov@cfi.lu.lv and Department of Physical Chemistry of Solids, Max Planck Institute for Solid State Research, Stuttgart, Germany.
authors.forename: E A
authors.identifier: http://orcid.org/0000-0001-6662-1705
authors.initials: EA
authors.lastname: Kotomin
authors.affiliation: Department of Theoretical Physics and Computer Modelling, Institute of Solid State Physics, University of Latvia, Riga, Latvia. gryaznov@cfi.lu.lv.
authors.forename: A
authors.initials: A
authors.lastname: Chesnokov
authors.affiliation: Department of Physical Chemistry of Solids, Max Planck Institute for Solid State Research, Stuttgart, Germany.
authors.forename: J
authors.initials: J
authors.lastname: Maier