A full implementation of the response iteration scheme for density functional calculations

This is the authors' abstract. We don't add key points for this paper.

We describe the implementation of the response algorithm for solving the Kohn-Sham equations and equations of similar structure that appear frequently in calculations of the structure of inhomogeneous many-body systems. The algorithm solves directly for the (spin) density. We have implemented the method for arbitrary geometries in a way that does not need the explicit calculation of unoccupied states. We apply the method to the two sample molecules C6H6 and C60, as well as to the metal clusters Mg30 and Na40. Starting from very rough guesses for the initial electron density, convergence is reached significantly faster than with previous update algorithms.