Skip to Main content Skip to Navigation
Journal articles

A Parallel Iterative Method for Computing Molecular Absorption Spectra

Abstract : We describe a fast parallel iterative method for computing molecular absorption spectra within TDDFT linear response and using the LCAO method. We use a local basis of "dominant products" to parametrize the space of orbital products that occur in the LCAO approach. In this basis, the dynamical polarizability is computed iteratively within an appropriate Krylov subspace. The iterative procedure uses a a matrix-free GMRES method to determine the (interacting) density response. The resulting code is about one order of magnitude faster than our previous full-matrix method. This acceleration makes the speed of our TDDFT code comparable with codes based on Casida's equation. The implementation of our method uses hybrid MPI and OpenMP parallelization in which load balancing and memory access are optimized. To validate our approach and to establish benchmarks, we compute spectra of large molecules on various types of parallel machines. The methods developed here are fairly general and we believe they will find useful applications in molecular physics/chemistry, even for problems that are beyond TDDFT, such as organic semiconductors, particularly in photovoltaics.
Complete list of metadatas
Contributor : Olivier Coulaud <>
Submitted on : Tuesday, June 1, 2010 - 9:14:37 AM
Last modification on : Thursday, September 6, 2018 - 2:58:12 PM

Links full text




Peter Koval, Dietrich Foerster, Olivier Coulaud. A Parallel Iterative Method for Computing Molecular Absorption Spectra. Journal of Chemical Theory and Computation, American Chemical Society, 2010, 6 (9), pp.2654-2668. ⟨10.1021/ct100280x⟩. ⟨inria-00488048⟩



Record views