Finite-size supercell correction schemes for charged defect calculations

Various schemes for correcting the finite-size supercell errors in the case of charged defect calculations are analyzed and their performance for a series of defect systems is compared. We focus on the schemes proposed by Makov and Payne (MP), Freysoldt, Neugebauer, and Van de Walle (FNV), and Lany and Zunger (LZ). The role of the potential alignment is also assessed. We demonstrate a connection between the defect charge distribution and the potential alignment, which establishes a relation between the MP and FNV schemes. Calculations are performed using supercells of various sizes and the corrected formation energies are compared to the values obtained by extrapolation to infinitely large supercells. For defects with localized charge distributions, we generally find that the FNV scheme improves upon the LZ one, while the MP scheme tends to overcorrect except for point-charge-like defects. We also encountered a class of defects, for which all the correction schemes fail to produce results consistent with the extrapolated values. This behavior is found to be caused by partial delocalization of the defect charge. We associate this effect to hybridization between the defect state and the band-edge states of the host. The occurrence of defect charge delocalization also reflects in the evolution of the defect Kohn-Sham levels with increasing supercell size. We discuss the physical relevance of the latter class of defects.

[1]  Alex Zunger,et al.  Accurate prediction of defect properties in density functional supercell calculations , 2009 .

[2]  A. Zunger,et al.  Many-body GW calculation of the oxygen vacancy in ZnO , 2009, 0910.2962.

[3]  A. Pasquarello,et al.  Defect energy levels in density functional calculations: alignment and band gap problem. , 2008, Physical review letters.

[4]  R. L. Johnston,et al.  Temperature and pressure variation of the refractive index of diamond. , 1977, Applied optics.

[5]  A. Pasquarello,et al.  Defect levels through hybrid density functionals: Insights and applications , 2011 .

[6]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[7]  T. Bredow,et al.  Color centers in NaCl by hybrid functionals , 2010, 1006.3935.

[8]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[9]  Payne,et al.  Periodic boundary conditions in ab initio calculations. , 1995, Physical review. B, Condensed matter.

[10]  Jianmin Tao,et al.  Tests of a ladder of density functionals for bulk solids and surfaces , 2004 .

[11]  P. P. Ewald.  Die Berechnung optischer und elektrostatischer Gitterpotentiale , 1921 .

[12]  G. Fischer.  Numerical data and functional relationships in science and technology , 1987 .

[13]  Fabien Bruneval,et al.  Understanding and correcting the spurious interactions in charged supercells , 2011 .

[14]  A. L. McWhorter,et al.  Polarization and Intensity of Raman Scattering from Plasmons and Phonons in Gallium Arsenide , 1967 .

[15]  R. Nieminen,et al.  Convergence of supercell calculations for point defects in semiconductors: Vacancy in silicon , 1998 .

[16]  F. Bechstedt,et al.  Linear optical properties in the projector-augmented wave methodology , 2006 .

[17]  Kth,et al.  Finite-size scaling as a cure for supercell approximation errors in calculations of neutral native defects in InP , 2004, cond-mat/0512306.

[18]  Schultz,et al.  Charged local defects in extended systems , 2000, Physical review letters.

[19]  G. Kresse,et al.  Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set , 1996 .

[20]  Zhang,et al.  Chemical potential dependence of defect formation energies in GaAs: Application to Ga self-diffusion. , 1991, Physical review letters.

[21]  Alex Zunger,et al.  Assessment of correction methods for the band-gap problem and for finite-size effects in supercell defect calculations: Case studies for ZnO and GaAs , 2008 .

[22]  Fabiano Corsetti,et al.  System-size convergence of point defect properties: The case of the silicon vacancy , 2010, 1010.3921.

[23]  R. Nieminen,et al.  Density-functional calculations of defect formation energies using supercell methods: Defects in diamond , 2005 .

[24]  Boris Kozinsky,et al.  Electrostatics in periodic boundary conditions and real-space corrections , 2007, 0709.4647.

[25]  Nicholas D. M. Hine,et al.  Supercell size scaling of density functional theory formation energies of charged defects , 2009 .

[26]  Alfredo Pasquarello,et al.  Assessing the accuracy of hybrid functionals in the determination of defect levels: Application to the As antisite in GaAs , 2011 .

[27]  Alfredo Pasquarello,et al.  Ab initio molecular dynamics in a finite homogeneous electric field. , 2002, Physical review letters.

[28]  Sadao Adachi,et al.  Optical Constants of Crystalline and Amorphous Semiconductors , 1999 .

[29]  G. Kresse,et al.  Ab initio molecular dynamics for liquid metals. , 1993 .

[30]  Jörg Neugebauer,et al.  Electrostatic interactions between charged defects in supercells , 2011 .

[31]  N. Marzari,et al.  Electrostatics in periodic boundary conditions and real-space corrections (vol 77, 115139, 2008) , 2011 .

[32]  G. Pacchioni,et al.  Theoretical description of hole localization in a quartz Al center: The importance of exact electron exchange , 2000 .

[33]  C. Walle,et al.  First-principles calculations for defects and impurities: Applications to III-nitrides , 2004 .

[34]  S. S. Mitra,et al.  Temperature Dependence of Infrared Dispersion in Ionic Crystals LiF and MgO , 1966 .

[35]  K. Burke,et al.  Rationale for mixing exact exchange with density functional approximations , 1996 .

[36]  A. Pasquarello,et al.  Effect of improved band-gap description in density functional theory on defect energy levels in α-quartz , 2007 .

[37]  Baroni,et al.  Ab initio calculation of the macroscopic dielectric constant in silicon. , 1986, Physical review. B, Condensed matter.

[38]  F. Giustino,et al.  Theory of atomic-scale dielectric permittivity at insulator interfaces , 2005 .

[39]  安達 定雄,et al.  Optical constants of crystalline and amorphous semiconductors : numerical data and graphical information , 1999 .

[40]  A. Pasquarello,et al.  Band-edge problem in the theoretical determination of defect energy levels: The O vacancy in ZnO as a benchmark case , 2011, 1204.4127.

[41]  C. Castleton,et al.  Managing the supercell approximation for charged defects in semiconductors: Finite-size scaling, charge correction factors, the band-gap problem, and the ab initio dielectric constant , 2006 .

[42]  C. Freysoldt,et al.  Fully ab initio finite-size corrections for charged-defect supercell calculations. , 2009, Physical review letters.

[43]  R. Fox,et al.  Classical Electrodynamics, 3rd ed. , 1999 .

[44]  M Leslie,et al.  The energy and elastic dipole tensor of defects in ionic crystals calculated by the supercell method , 1985 .