A novel parallel algorithm for large‐scale Fock matrix construction with small locally distributed memory architectures: RT parallel algorithm

We developed a novel parallel algorithm for large‐scale Fock matrix calculation with small locally distributed memory architectures, and named it the “RT parallel algorithm.” The RT parallel algorithm actively involves the concept of integral screening, which is indispensable for reduction of computing times with large‐scale biological molecules. The primary characteristic of this algorithm is parallel efficiency, which is achieved by well‐balanced reduction of both communicating and computing volume. Only the density matrix data necessary for Fock matrix calculations are communicated, and the data once communicated are reutilized for calculations as many times as possible. The RT parallel algorithm is a scalable method because required memory volume does not depend on the number of basis functions. This algorithm automatically includes a partial summing technique that is indispensable for maintaining computing accuracy, and can also include some conventional methods to reduce calculation times. In our analysis, the RT parallel algorithm had better performance than other methods for massively parallel processors. The RT parallel algorithm is most suitable for massively parallel and distributed Fock matrix calculations for large‐scale biological molecules with more than thousands of basis functions. © 2002 Wiley Periodicals, Inc. J Comput Chem 23: 1337–1346, 2002

[1]  Matt Challacombe,et al.  Linear scaling computation of the Fock matrix. V. Hierarchical Cubature for numerical integration of the exchange-correlation matrix , 2000 .

[2]  Thomas R. Furlani,et al.  Implementation of a parallel direct SCF algorithm on distributed memory computers , 1995, J. Comput. Chem..

[3]  Eric Schwegler,et al.  Linear scaling computation of the Fock matrix. II. Rigorous bounds on exchange integrals and incremental Fock build , 1997 .

[4]  Gustavo E. Scuseria,et al.  A quantitative study of the scaling properties of the Hartree–Fock method , 1995 .

[5]  Fumitoshi Sato,et al.  Distributed parallel processing by using the object-oriented technology in ProteinDF program for all-electron calculations on proteins , 2001 .

[6]  Ching-Hsing Yu,et al.  Ab Initio Geometry Determinations of Proteins. 1. Crambin , 1998 .

[7]  H. Fritsche,et al.  Reviews in Computational Chemistry IV , 1994 .

[8]  Toshikazu Ebisuzaki,et al.  A Highly Parallelized Special-Purpose Computer for Many-Body Simulations with an Arbitrary Central Force: MD-GRAPE , 1996 .

[9]  Michael E. Colvin,et al.  Parallel direct SCF for large-scale calculations , 1993 .

[10]  Kazutoshi Tanabe,et al.  Is large‐scale ab initio Hartree‐Fock calculation chemically accurate? Toward improved calculation of biological molecule properties , 1999 .

[11]  Mark S. Gordon,et al.  General atomic and molecular electronic structure system , 1993, J. Comput. Chem..

[12]  Takashi Amisaki,et al.  Development of MD Engine: High‐speed accelerator with parallel processor design for molecular dynamics simulations , 1999 .

[13]  Rick Stevens,et al.  Toward high‐performance computational chemistry: II. A scalable self‐consistent field program , 1996 .

[14]  James J. P. Stewart,et al.  Application of localized molecular orbitals to the solution of semiempirical self‐consistent field equations , 1996 .

[15]  Timothy G. Mattson.  Parallel Computing in Computational Chemistry , 1995 .

[16]  Robert J. Harrison,et al.  Toward high‐performance computational chemistry: I. Scalable Fock matrix construction algorithms , 1996 .

[17]  Jürgen Gauss,et al.  An unconventional scf method for calculations on large molecules , 1986 .

[18]  Marco Häser,et al.  Improvements on the direct SCF method , 1989 .

[19]  Haruki Nakamura,et al.  Ab initio MO study of the chlorophyll dimer in the photosynthetic reaction center. I. A theoretical treatment of the electrostatic field created by the surrounding proteins , 1997 .

[20]  Kenichiro Tsuda,et al.  Ab initio MO studies of interaction mechanisms of Protein Kinase C with cell membranes , 2001 .

[21]  Kazuaki Murakami,et al.  Development of MOE (molecular orbital calculation engine) , 1999 .

[22]  Itai Panas,et al.  ABINITIO METHODS FOR LARGE SYSTEMS , 1991 .

[23]  Fumitoshi Sato,et al.  Calculation of all-electron wavefunction of hemoprotein cytochrome c by density functional theory , 2001 .

[24]  J. Almlöf,et al.  Principles for a direct SCF approach to LICAO–MOab‐initio calculations , 1982 .

[25]  Stefan Brode,et al.  Parallel direct SCF and gradient program for workstation clusters , 1993, J. Comput. Chem..