A Master-Worker Type Eigensolver for Molecular Orbital Computations

We consider a parallel method for solving generalized eigenvalue problems that arise from molecular orbital computations. We use a moment-based method that finds several eigenvalues and their corresponding eigenvectors in a given domain, which is suitable for masterworker type parallel programming models. The computation of eigenvalues using explicit moments is sometimes numerically unstable. We show that a Rayleigh-Ritz procedure can be used to avoid the use of explicit moments. As a test problem, we use the matrices that arise in the calculation of molecular orbitals. We report the performance of the application of the proposed method with several PC clusters connected through a hybrid MPI and GridRPC system.

[1]  P. Deuflhard,et al.  Large Scale Scientific Computing , 1987 .

[2]  Mardochée Magolu monga Made,et al.  Incomplete factorization-based preconditionings for solving the Helmholtz equation , 2001 .

[3]  Z. Bai.  Krylov subspace techniques for reduced-order modeling of large-scale dynamical systems , 2002 .

[4]  T. Sakurai,et al.  On Locating Clusters of Zeros of Analytic Functions , 1999 .

[5]  T. Sakurai,et al.  A projection method for generalized eigenvalue problems using numerical integration , 2003 .

[6]  K. Kitaura,et al.  Definition of molecular orbitals in fragment molecular orbital method , 2002 .

[7]  H. V. D. Vorst,et al.  A Petrov-Galerkin type method for solving Axk=b, where A is symmetric complex , 1990 .

[8]  Axel Ruhe.  Rational Krylov algorithms for nonsymmetric eigenvalue problems. II. matrix pairs , 1994 .

[9]  W. Arnoldi.  The principle of minimized iterations in the solution of the matrix eigenvalue problem , 1951 .

[10]  T. Sakurai,et al.  A perturbation result for generalized eigenvalue problems and its application to error estimation in a quadrature method for computing zeros of analytic functions , 2003 .

[11]  Tetsuya Sakurai,et al.  A Hybrid Parallel Method for Large Sparse Eigenvalue Problems on a Grid Computing Environment Using Ninf-G/MPI , 2005, LSSC.

[12]  Mitsuhisa Sato,et al.  A Parallel Method for Large Sparse Generalized Eigenvalue Problems by OmniRPC in a Grid Environment , 2004, PARA.

[13]  T. Sakurai,et al.  An error analysis of two related quadrature methods for computing zeros of analytic functions , 2003 .

[14]  T. Sakurai,et al.  A projection method for generalized eigenvalue problems , 2002 .