Parallelization of four-component calculations. II. Symmetry-driven parallelization of the 4-Spinor CCSD algorithm

Given the importance of the Coupled-cluster (CC) method as an efficient and accurate way to take electron correlation into account, we extend the parallelization technique in the second part of this series also to the 4-Spinor CCSD algorithm implemented in the Dirac-Fock packages DIRAC and MOLFDIR. The present implementation is based on the availability of the transformed molecular two-electron integrals on an external storage medium. The linearity of the CC equations in these two-electron integrals is used in a parallelization strategy that is based on distribution of the two largest integral classes that carry three or four virtual spinor indices. The corresponding partial contributions to the T(1) and T(2) amplitudes are calculated on each node and added using Message Passing Interface (MPI) library calls. Although we did not employ a master/slave principle, one specific node was assigned to also perform the remaining serial parts of the algorithm. In the critical sections considerable savings in storage requirements and computer time could be achieved, and this allows for computations on larger systems in the framework of four-component theory.