Ab-initio calculation of intrinsic spin Hall conductivity (SHC) generally requires a strict convergence criterion and a dense k-point mesh to sample the Brillouin zone, making its convergence challenging and time-consuming. Here we present a scheme for efficiently and accurately calculating SHC based on maximally localized Wannier function (MLWF). The quantities needed by the Kubo formula of SHC are derived in the space of MLWF and it is shown that only the Hamiltonian, the overlap and the spin operator matrices are required from the initial ab-initio calculation. The computation of these matrices and the interpolation of Kubo formula on a dense k-point mesh can be easily achieved. We validate our results by prototypical calculations on fcc Pt and GaAs, which demonstrate that the Wannier interpolation approach is of high accuracy and efficiency. Calculations of α-Ta and β-Ta show that SHC of β-Ta is 2.7 times of α-Ta, while both have the opposite sign relative to fcc Pt and are an order of magnitude smaller than Pt. The calculated spin Hall angle of −0.156, is quite consistent with previous experiment on β-Ta, further suggesting intrinsic contribution may dominate in β-Ta. Our approach could facilitate large-scale SHC calculations, and may benefit the discovery of materials with high intrinsic SHC. arXiv:1810.07637v1 [cond-mat.mtrl-sci] 17 Oct 2018 (0) nk to