Recursive hybrid tetrahedron method for Brillouin-zone integration

A recursive extension of the hybrid tetrahedron method for Brillouin-zone integration is proposed, allowing iterative tetrahedron refinement and significantly reducing the error from the linear tetrahedron method. The Brillouin-zone integral is expressed as a weighted sum on the initial grid, with integral weights collected recursively from the finest grid. Our method is capable of simultaneously handling multiple singularities in the integrand and thus may provide practical solutions to various Brillouin-zone integral tasks encountered in realistic calculations, including the computation of response and spectral function with superior sampling convergence. We demonstrate its effectiveness through numerical calculations of the density response functions of two model Hamiltonians and one real material system, the face-centered cubic cobalt.