Susceptibilities for first principles band structures

We present a parallel implementation of a new method for calculating the unenhanced susceptibility proposed by us recently. Our implementation uses the first principles LMTO band structure within the tight binding approach as input to calculate the joint density of states. The susceptibility is then obtained by integrating over the product of the joint density of states and a Lindhard function. Our program, which has a simple friendly user interface, runs on the PC with a quadputer board, a Meiko Surface running CSTools powered either by T800 or i860 compute boards and the Intel iPSC/860 hypercube in Daresbury. Our method incorporates the troublesome matrix elements naturally and our results on Pd and Ni show that the decrease in χ( q ) as we go away from the Brillouin zone centre is due mainly to the matrix elements rather than to the band energies.