Korringa–Kohn–Rostoker method

The Korringa–Kohn–Rostoker method or KKR method is used to calculate the electronic band structure of periodic solids. In the derivation of the method using multiple scattering theory by Jan Korringa[1] and the derivation based on the Kohn and Rostoker variational method,[2] the muffin-tin approximation was used.[3] Later calculations are done with full potentials having no shape restrictions.[4][5]

Introduction

All solids in their ideal state are single crystals with the atoms arranged on a periodic lattice. In condensed matter physics, the properties of such solids are explained on the basis of their electronic structure. This requires the solution of a complicated many-electron problem, but the density functional theory of Walter Kohn makes it possible to reduce it to the solution of a Schroedinger equation with a one-electron periodic potential. The problem is further simplified with the use of group theory and in particular Bloch's theorem, which leads to the result that the energy eigenvalues depend on the crystal momentum ${\displaystyle {\bf {k}}}$ and are divided into bands. Band theory is used to calculate the eigenvalues and wave functions.

Many band theory methods have been proposed over the years. Some of the most widely used, such as the electronic structure programs VASP and WIEN2k, make use of approximations so that acceptable accuracy can be achieved with a minimum of computer resources. The KKR method is chosen when the primary goal is high accuracy.

The parameters obtained from reliable band-theory calculations are useful in theoretical studies of problems, such as superconductivity, for which the density functional theory does not apply.

Mathematical formulation

The KKR band theory equations for space-filling non-spherical potentials are derived in books[4][5] and in the article on multiple scattering theory.

The wave function near site ${\displaystyle j}$ is determined by the coefficients ${\displaystyle c_{l'm'}^{j}}$. According to Bloch's theorem, these coefficients differ only through a phase factor ${\displaystyle c_{l'm'}^{j}={e^{-i{\bf {k}}\cdot {{\bf {R}}_{j}}}}{c_{l'm'}}\left(E,{\bf {k}}\right)}$. The ${\displaystyle {c_{l'm'}}\left(E,{\bf {k}}\right)}$ satisfy the homogeneous equations

${\displaystyle \sum \limits _{l'm'}{{M_{lm,l'm'}}\left({E,{\bf {k}}}\right){c_{l'm'}}\left(E,{\bf {k}}\right)=0}}$,

where ${\displaystyle {M_{lm,l'm'}}\left({E,{\bf {k}}}\right)={m_{lm,l'm'}}\left(E\right)-{A_{lm,l'm'}}\left({E,{\bf {k}}}\right)}$ and ${\displaystyle {A_{lm,l'm'}}\left({E,{\bf {k}}}\right)=\sum \limits _{j}{e^{i{\bf {k}}\cdot {{\bf {R}}_{ij}}}}{g_{lm,l'm'}}\left({E,{{\bf {R}}_{ij}}}\right)}$.

The ${\displaystyle {m_{lm,l'm'}}\left(E\right)}$ is the inverse of the scattering matrix ${\displaystyle {t_{lm,l'm'}}\left(E\right)}$ calculated with the non-spherical potential for the site. As pointed out by Korringa,[1] Ewald derived a summation process that makes it possible to calculate the structure constants, ${\displaystyle {A_{lm,l'm'}}\left({E,{\bf {k}}}\right)}$. The energy eigenvalues of the periodic solid for a particular ${\displaystyle {\bf {k}}}$, ${\displaystyle {E_{b}}\left({\bf {k}}\right)}$, are the roots of the equation ${\displaystyle \det {\bf {M}}\left({E,{\bf {k}}}\right)=0}$. The eigenfunctions are found by solving for the ${\displaystyle {c_{l,m}}\left(E,{\bf {k}}\right)}$ with ${\displaystyle E={E_{b}}\left({\bf {k}}\right)}$. By ignoring all contributions that correspond to an angular momentum ${\displaystyle l}$ greater than ${\displaystyle {l_{\max }}}$, they have dimension ${\displaystyle {\left({{l_{\max }}+1}\right)^{2}}}$.

In the original derivations of the KKR method, spherically symmetric muffin-tin potentials were used. Such potentials have the advantage that the inverse of the scattering matrix is diagonal in ${\displaystyle l}$

${\displaystyle {m_{lm,l'm'}}=\left[{\alpha \cot {\delta _{l}}\left(E\right)-i\alpha }\right]{\delta _{l,l'}}{\delta _{m,m'}}}$,

where ${\displaystyle {\delta _{l}}\left(E\right)}$ is the scattering phase shift that appears in the partial wave analysis in scattering theory. The muffin-tin approximation is good for closely packed metals, but it does not work well for ionic solids like semiconductors. It also leads to errors in calculations of interatomic forces.

References

1. J. Korringa (1947). "On the calculation of the energy of a Bloch wave in a metal". Physica. XIII (6–7): 392–400. Bibcode:1947Phy....13..392K. doi:10.1016/0031-8914(47)90013-x.
2. W. Kohn, N. Rostoker (1954). "Solution of the Schrödinger Equation in Periodic Lattices with an Application to Metallic Lithium". Phys. Rev. 94 (5): 1111–1120. Bibcode:1954PhRv...94.1111K. doi:10.1103/physrev.94.1111.
3. W. Jones, N. H. March (1973). Theoretical Solid State Physics. Wiley and Sons - Dover Publications. ISBN 0-486-65015-4.
4. Jan Zabloudil; Robert Hammerling; Laszlo Szunyogh; Peter Weinberger (2010) [2005]. Electron Scattering in Solid Matter: A Theoretical and Computational Treatise (Softcover reprint of hardcover 1st ed. 2005 ed.). Springer. ISBN 978-3642061387.
5. Yang Wang; G. Malcolm Stocks; J. Sam Faulkner (2015). Multiple Scattering beta edition (Kindle Interactive ed.). Amazon. ASIN B015NFAN6M.