Numerically improved computational scheme for the optical conductivity   tensor in layered systems
  The contour integration technique applied to calculate the optical conductivity tensor at finite temperatures in the case of layered systems within the framework of the spin-polarized relativistic screened Korringa-Kohn-Rostoker band structure method is improved from the computational point of view by applying the Gauss-Konrod quadrature for the integrals along the different parts of the contour and by designing a cumulative special points scheme for two-dimensional Brillouin zone integrals corresponding to cubic systems. 
A. Vernes, L. Szunyogh and P. Weinberger