TY - JOUR

T1 - Self-consistent full-potential total-energy Korringa-Kohn-Rostoker band-structure method

T2 - Application to silicon

AU - Bei Der Kellen, S.

AU - Oh, Yoonsik

AU - Badralexe, E.

AU - Freeman, A. J.

N1 - Copyright:
Copyright 2015 Elsevier B.V., All rights reserved.

PY - 1995

Y1 - 1995

N2 - The self-consistent full-potential total-energy Korringa-Kohn-Rostoker (KKR) coupled channel equation method for full-potential electronic band-structure determinations is applied to the case of silicon. For efficiency, the coupled channel equations are partially decoupled by using the point group symmetry of each basis atom, and solved on a linear radial mesh by using the Coulomb functions close to the singularity of the potential at the origin. In the KKR energy search, Chebyshev approximations for the scattering matrices as functions of energy save solving the coupled channel equations at many different energies. The results include the self-consistent anisotropic potential, the charge density, the energy band structure, the equilibrium lattice constant, and the bulk modulus. As expected, the energy bands agree well with results from other band-structure methods. We find that the cutoff of the spherical multipole expansion of potential and charge density must be chosen independently of the cutoff for the scattering matrices that determines the dimension of the secular matrix. In the case of silicon, we find that an lmax of 8 is sufficient for the potential and charge density, while the scattering matrices need at least lmax=6. The equilibrium lattice constant converges to the experimental value with increasing lmax and the bulk modulus converges to 96 GPa at lmax=6.

AB - The self-consistent full-potential total-energy Korringa-Kohn-Rostoker (KKR) coupled channel equation method for full-potential electronic band-structure determinations is applied to the case of silicon. For efficiency, the coupled channel equations are partially decoupled by using the point group symmetry of each basis atom, and solved on a linear radial mesh by using the Coulomb functions close to the singularity of the potential at the origin. In the KKR energy search, Chebyshev approximations for the scattering matrices as functions of energy save solving the coupled channel equations at many different energies. The results include the self-consistent anisotropic potential, the charge density, the energy band structure, the equilibrium lattice constant, and the bulk modulus. As expected, the energy bands agree well with results from other band-structure methods. We find that the cutoff of the spherical multipole expansion of potential and charge density must be chosen independently of the cutoff for the scattering matrices that determines the dimension of the secular matrix. In the case of silicon, we find that an lmax of 8 is sufficient for the potential and charge density, while the scattering matrices need at least lmax=6. The equilibrium lattice constant converges to the experimental value with increasing lmax and the bulk modulus converges to 96 GPa at lmax=6.

UR - http://www.scopus.com/inward/record.url?scp=0342568540&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0342568540&partnerID=8YFLogxK

U2 - 10.1103/PhysRevB.51.9560

DO - 10.1103/PhysRevB.51.9560

M3 - Article

AN - SCOPUS:0342568540

VL - 51

SP - 9560

EP - 9568

JO - Physical Review B

JF - Physical Review B

SN - 0163-1829

IS - 15

ER -