Quantum dynamics on massively parallel computers: Efficient numerical implementation for preconditioned linear solvers and eigensolvers

Wenwu Chen, Bill Poirier

Research output: Contribution to journalArticlepeer-review

16 Scopus citations

Abstract

The eigenvalue/eigenvector and linear solve problems arising in computational quantum dynamics applications (e.g. rovibrational spectroscopy, reaction cross-sections, etc.) often involve large sparse matrices that exhibit a certain block structure. In such cases, specialized iterative methods that employ optimal separable basis (OSB) preconditioners (derived from a block Jacobi diagonalization procedure) have been found to be very efficient, vis-à-vis reducing the required CPU effort on serial computing platforms. Recently,1,2 a parallel implementation was introduced, based on a nonstandard domain decomposition scheme. Near-perfect parallel scalability was observed for the OSB preconditioner construction routines up to hundreds of nodes; however, the fundamental matrixvector product operation itself was found not to scale well, in general. In addition, the number of nodes was selectively chosen, so as to ensure perfect load balancing. In this paper, two essential improvements are discussed: (1) new algorithm for the matrixvector product operation with greatly improved parallel scalability and (2) generalization for arbitrary number of nodes and basis sizes. These improvements render the resultant parallel quantum dynamics codes suitable for robust application to a wide range of real molecular problems, running on massively parallel computing architectures.

Original languageEnglish
Pages (from-to)825-846
Number of pages22
JournalJournal of Theoretical and Computational Chemistry
Volume9
Issue number5
DOIs
StatePublished - Oct 2010

Keywords

  • Quantum dynamics
  • block Jacobi
  • chemical physics
  • eigensolver
  • linear solver
  • matrixvector product
  • parallel computing
  • preconditioning
  • sparse matrix

Fingerprint Dive into the research topics of 'Quantum dynamics on massively parallel computers: Efficient numerical implementation for preconditioned linear solvers and eigensolvers'. Together they form a unique fingerprint.

Cite this