Abstract:Formation energies of charged point defects in semiconductors are calculated using periodic supercells, which entail a divergence arising from long-range Coulombic interactions. The divergence is typically removed by the so-called jellium approach. Recently, Wu, Zhang and Pantelides [WZP, Phys. Rev. Lett. 119, 105501 (2017)] traced the origin of the divergence to the assumption that charged defects are formed by physically removing electrons from or adding electrons to the crystal, violating charge neutrality, a key principle of statistical mechanics that determines the Fermi level. An alternative theory was constructed by recognizing that "charged" defects form by trading carriers with the energy bands, whereby supercells are always charge-neutral so that no divergence is present and no ad-hoc procedures need to be adopted for calculations. Here we give a more detailed exposition of the foundations of both methods and show that the jellium approach can be derived from the statistical-mechanics-backed WZP definition by steps whose validity cannot be assessed a priori. In particular, the divergence appears when the charge density of band carriers is dropped, leaving a supercharged crystal. In the case of charged defects in two-dimensional (2D) materials, unphysical fields appear in vacuum regions. None of these pathological features are present in the reformulated theory. Finally, we report new calculations in both bulk and 2D materials. The WZP approach yields formation energies that differ from jellium values by up to ~1 eV. By analyzing the spatial distribution of wave functions and defect potentials, we provide insights into the inner workings of both methods and demonstrate that the failure of the jellium approach to include the neutralizing electron density of band carriers, as is the case in the physical system, is responsible for the numerical differences between the two methods.

Importance of elastic finite-size effects: neutral defects in ionic compounds

Point defect modeling in materials: coupling ab initio and elasticity approaches

Finite-size corrections to the free energy of crystalline solids

Imperfections are not 0 K: free energy of point defects in crystals

The continuum elastic and atomistic viewpoints on the formation volume and strain energy of a point defect

Molecular Dynamics Simulation of Coherent Interfaces in Fluorite Heterostructures

First-principles calculations for point defects in solids

Physics of band-filling correction in defect calculations of solid-state materials

Finite-Temperature Atomistic and Continuum Stress Fields of Coherent Precipitates with a Small Lattice Misfit

Linear Scaling DFT for defects in metals

Defect creation and Diffusion under electric fields from first-principles: the prototypical case of silicon dioxide

A realistic dimension-independent approach for charged defect calculations in semiconductors

Elastic Fields In Quantum Dot Structures With Arbitrary Shapes And Interface Effects

Defects from phonons: Atomic transport by concerted motion in simple crystalline metals

Finite Field Methods for the Supercell Modelling of Charged Insulator-Electrolyte Interfaces

Investigating finite-size effects in molecular dynamics simulations of ion diffusion, heat transport, and thermal motion in superionic materials

Lattice Green function for extended defect calculations: Computation and error estimation with long-range forces

Image charge interaction correction in charged-defect calculations

Comparison of charged-defect finite-size supercell correction methods in a general framework

Aperiodic defects in periodic solids

Formation energies of charged defects in two-dimensional materials -- resolution of long-standing difficulties