Methods for parallel computation of SCF NMR chemical shifts by GIAO method: Efficient integral calculation, multi-Fock algorithm, and pseudodiagonalization

by Wolinski, K.; Haacke, R.; Hinton, J. F.; Pulay, P.

We implemented our gauge-including atomic orbital (GIAO) NMR chemical shielding program on a workstation cluster, using the parallel virtual machine (PVM) message-passing system. On a modest number of nodes, we achieved close to Linear speedup. This program is characterized by several novel features. It uses the new integral program of Wolinski that calculates integrals in vectorized batches, increases efficiency, and simplifies parallelization. The self-consistent field (SCF) step includes a multi-Fock algorithm, i.e., the simultaneous calculation of several Fock matrices with the same integral set, increasing the efficiency of the direct SCF procedure. The SCF diagonalization step, which is difficult to parallelize, has been replaced by pseudo-diagonalization. The latter, widely used in semiempirical programs, becomes important in ab initio type calculations above a certain size, because the ultimate scaling of the diagonalization step is steeper than that of integral computation. Examples of the calculation of the NMR shieldings in large systems at the SCF level are shown. Parallelization of the density functional code is underway.

Journal
Journal of Computational Chemistry
Volume
18
Issue
6
Year
1997
Start Page
816-825
URL
https://dx.doi.org/10.1002/(sici)1096-987x(19970430)18:6<816::aid-jcc7>3.0.co;2-v
ISBN/ISSN
1096-987X; 0192-8651
DOI
10.1002/(sici)1096-987x(19970430)18:6<816::aid-jcc7>3.0.co;2-v