An algorithm for the multi-reference configuration interaction method on distributed memory architectures

We have reformulated the multi-reference configuration interaction algorithm for massively parallel architectures with distributed memory, such that the communication effort associated with the matrix element computation is independent of the number of units of the machine. Exploiting the special structure of the MRCI Hilbert space, this formulation allows the evaluation of all matrix elements in six distinct phases—independent of the number of electrons, the number of orbitals, the number of references and of the number of nodes of the computer used. Using an object-oriented approach, we have implemented and tested the a symmetry-adapted determinantal version of the algorithm in C++ on an IBM-SP2 computer. Our initial results for benchmark calculations on N2 and O3 for Hilbert spaces of up to 12×106 determinants indicate that the vector redistribution effort is negligible in actual calculations. The remaining difficulties to perform calculations for significantly larger Hilbert spaces are critically addr...

[1]  Robert J. Buenker,et al.  Individualized configuration selection in CI calculations with subsequent energy extrapolation , 1974 .

[2]  Robert J. Harrison,et al.  A massively parallel multireference configuration interaction program: The parallel COLUMBUS program , 1997 .

[3]  Michel Dupuis,et al.  Parallel computation of the Moller–Plesset second‐order contribution to the electronic correlation energy , 1988 .

[4]  B. Roos,et al.  A new method for large-scale Cl calculations , 1972 .

[5]  Rodney J. Bartlett,et al.  Multi-reference averaged quadratic coupled-cluster method: a size-extensive modification of multi-reference CI , 1993 .

[6]  Robert J. Gdanitz,et al.  The averaged coupled-pair functional (ACPF): A size-extensive modification of MR CI(SD) , 1988 .

[7]  Bernard Pullman,et al.  The World of Quantum Chemistry , 1974 .

[8]  Robert J. Buenker,et al.  Configuration interaction calculations for the N2 molecule and its three lowest dissociation limits , 1977 .

[9]  P. Löwdin.  Quantum Theory of Many-Particle Systems. III. Extension of the Hartree-Fock Scheme to Include Degenerate Systems and Correlation Effects , 1955 .

[10]  Péter G. Szalay,et al.  High-quality theoretical potential energy surface for Be2 by using the multireference averaged quadratic coupled-cluster (MR-AQCC) method and large basis sets , 1996 .

[11]  Leszek Meissner,et al.  Size-consistency corrections for configuration interaction calculations , 1988 .

[12]  Ernest R. Davidson,et al.  Configuration interaction calculations on the nitrogen molecule , 1974 .

[13]  G. Diercksen,et al.  Methods in Computational Molecular Physics , 1983 .

[14]  Martyn F. Guest,et al.  Parallelism in computational chemistry , 1993 .

[15]  Robert J. Buenker,et al.  Energy extrapolation in CI calculations , 1975 .

[16]  Robert J. Harrison,et al.  Approximating full configuration interaction with selected configuration interaction and perturbation theory , 1991 .

[17]  David E. Bernholdt,et al.  Orbital‐invariant second‐order many‐body perturbation theory on parallel computers: An approach for large molecules , 1995 .

[18]  R. Harrison,et al.  AB Initio Molecular Electronic Structure on Parallel Computers , 1994 .

[19]  Robert J. Harrison,et al.  Computational chemistry on the FPS-X64 scientific computers , 1987 .

[20]  Rodney J. Bartlett,et al.  Approximately extensive modifications of the multireference configuration interaction method: A theoretical and practical analysis , 1995 .

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

[22]  David E. Woon,et al.  Gaussian basis sets for use in correlated molecular calculations. IV. Calculation of static electrical response properties , 1994 .