We present a massively parallel implementation of time-dependent density functional theory in real space, aimed at computing optical absorption spectra of realistic systems with hundreds of atoms from first principles. We provide details of the formalism and discuss its implementation, optimization, and efficient parallelization, as well as remaining limitations, in detail. The capabilities of the code are illustrated by calculations of optical properties of hydrogenated silicon quantum dots.
Bibliographical noteFunding Information:
We thank Zhongze Li for his help in parallelizing an early version of the coupling matrix construction code, and Yu Liang for his assistance in adapting the Incomplete Cholesky factorization preconditioning for use in this work. LK acknowledges the generous support of the Estelle Funk Foundation and the Delta Career Development Chair. We acknowledge support for this work by the National Science Foundation (under Grant ITR 0082094) and the Minnesota Supercomputing Institute.
Work supported by NSF grant NSF/ITR-0082094 and by the Minnesota Supercomputing Institute.