Acknowledgements

It is a pleasure to thank Professor James Talman, University of Western Ontario, who kindly donated important routines for 3j, 6j, 9j symbols and a fast Hankel transform. The ongoing encouragement and advice of the SIESTA development team, particularly Alberto García, Institut de Ciència de Materials de Barcelona and Daniel Sanchez-Portal, Centro de Física de Materiales Donostia-San Sebastián, and of Luc Giraud, Inria-Bordeaux Sud-Ouest, is gratefully acknowledged. This work was supported by A.N.R. grant CIS 2007 NOSSI. --