A Wigner Matrix Based Convolution Algorithm For Matrix Elements in the LCAO Method

The linear combination of atomic orbitals (LCAO) method uses a small basis set in exchange for expensive matrix element calculations. The most efficient approximation for the matrix element calculations is the two-center approximation (2CA) in tight binding (TB). In the 2CA, a variety of matrix elements are neglected with only "two-center integrals" (2CI) remaining. The 2CI are calculated efficiently by rotating to symmetrical coordinates where the integral is parameterized. This makes TB fast in exchange for diminished transferability. An ideal electronic structure method has both the efficiency of TB and the transferability of ab-initio methods. In this work, I expand the full crystal potential into multipoles where the resulting matrix elements are transformed into the form of 2CI between high angular momentum functions. The usual Slater-Koster formulae for TB are limited to $l\leq3$; to enable efficient evaluation of the full crystal potential 2CI, I derive a Wigner matrix based convolution algorithm (WMCA) that works for arbitrary angular momentum. Given a suitable method for generating a local ab-initio Kohn-Sham potential, the algorithm for calculating matrix elements is applicable to fully ab-initio LCAO methods (this is the subject of forthcoming work). In this paper, I apply the WMCA to silicon using a model crystal potential.