Valence electrons real-space energy

The relaxation of the structure in the KMC-DR method was done using an approach based on the density functional theory and linear combination of atomic orbitals implemented in the Siesta code [97]. The minimum basis set of localized numerical orbitals of Sankey type [98] was used for all atoms except silicon atoms near the interface, for which polarization functions were added to improve the description of the SiOx layer. The core electrons were replaced with norm-conserving Troullier-Martins pseudopotentials [99] (Zr atoms also include 4p electrons in the valence shell). Calculations were done in the local density approximation (LDA) of DFT. The grid in the real space for the calculation of matrix elements has an equivalent cutoff energy of 60 Ry. The standard diagonalization scheme with Pulay mixing was used to get a self-consistent solution. In the framework of the KMC-DR method, it is not necessary to perform an accurate optimization of the structure, since structure relaxation is performed many times. 

In our real-space core-valence separation, defines the boundary of the region around nucleus k that contains the core electrons associated with it. The pertinent energy components are ... 

---