A modified version of the method for calculating integrals over the Brillouin zone, which relies on the use of the orthogonal function series, is discussed. It is shown that adoption of the system of Chebyshev’s polynomials makes it possible, while retaining all the advantages of the original scheme, to overcome a number of its disadvantages and to improve the numerical stability in the case of high-degree polynomial approximation, to reduce the number of parameters in the method, and to rule out the development of regions with negative values of the density of states.
Read full abstract