Fast and Flexible Coupled Cluster Implementation.

A new coupled cluster singles and doubles with triples correction, CCSD(T), algorithm is presented. The new algorithm is implemented in object oriented C++, has a low memory footprint, fast execution time, low I/O overhead, and a flexible storage backend with the ability to use either distributed memory or a file system for storage. The algorithm is demonstrated to work well on single workstations, a small cluster, and a high-end Cray computer. With the new implementation, a CCSD(T) calculation with several hundred basis functions and a few dozen occupied orbitals can run in under a day on a single workstation. The algorithm has also been implemented for graphical processing unit (GPU) architecture, giving a modest improvement. Benchmarks are provided for both CPU and GPU hardware.

[1]  Jack J. Dongarra,et al.  Accelerating GPU Kernels for Dense Linear Algebra , 2010, VECPAR.

[2]  R. Bartlett,et al.  A full coupled‐cluster singles and doubles model: The inclusion of disconnected triples , 1982 .

[3]  J. Pople,et al.  Variational configuration interaction methods and comparison with perturbation theory , 2009 .

[4]  Karol Kowalski,et al.  Efficient computer implementation of the renormalized coupled-cluster methods: The R-CCSD[T], R-CCSD(T), CR-CCSD[T], and CR-CCSD(T) approaches , 2002 .

[5]  R J Bartlett,et al.  Parallel implementation of electronic structure energy, gradient, and Hessian calculations. , 2008, The Journal of chemical physics.

[6]  M. Head‐Gordon,et al.  A fifth-order perturbation comparison of electron correlation theories , 1989 .

[7]  Jarek Nieplocha,et al.  Advances, Applications and Performance of the Global Arrays Shared Memory Programming Toolkit , 2006, Int. J. High Perform. Comput. Appl..

[8]  Mark S. Gordon,et al.  New Multithreaded Hybrid CPU/GPU Approach to Hartree-Fock. , 2012, Journal of chemical theory and computation.

[9]  J. Cizek.  On the Correlation Problem in Atomic and Molecular Systems. Calculation of Wavefunction Components in Ursell-Type Expansion Using Quantum-Field Theoretical Methods , 1966 .

[10]  Henry F. Schaefer,et al.  Accelerating the convergence of the coupled-cluster approach: The use of the DIIS method , 1986 .

[11]  Mark S. Gordon,et al.  Chapter 41 – Advances in electronic structure theory: GAMESS a decade later , 2005 .

[12]  Martin Head-Gordon,et al.  Quadratic configuration interaction. A general technique for determining electron correlation energies , 1987 .

[13]  Peter Pulay,et al.  Parallel Calculation of Coupled Cluster Singles and Doubles Wave Functions Using Array Files. , 2007, Journal of chemical theory and computation.

[14]  Curtis L. Janssen,et al.  An efficient reformulation of the closed‐shell coupled cluster single and double excitation (CCSD) equations , 1988 .

[15]  Wilfried Meyer,et al.  Theory of self‐consistent electron pairs. An iterative method for correlated many‐electron wavefunctions , 1976 .

[16]  David E. Bernholdt,et al.  Parallel computational chemistry made easier: The development of NWChem , 1995 .

[17]  Debashis Mukherjee,et al.  Applications of a non-perturbative many-body formalism to general open-shell atomic and molecular problems: calculation of the ground and the lowest π-π* singlet and triplet energies and the first ionization potential of trans-butadiene , 1977 .

[18]  Brett M. Bode,et al.  Uncontracted Rys Quadrature Implementation of up to G Functions on Graphical Processing Units. , 2010, Journal of chemical theory and computation.

[19]  I. Lindgren.  A coupled-cluster approach to the many-body perturbation theory for open-shell systems , 2009 .

[20]  F. Coester,et al.  Short-range correlations in nuclear wave functions , 1960 .

[21]  Alistair P. Rendell,et al.  A direct coupled cluster algorithm for massively parallel computers , 1997 .

[22]  Josef Paldus,et al.  Correlation Problems in Atomic and Molecular Systems. IV. Extended Coupled-Pair Many-Electron Theory and Its Application to the B H 3 Molecule , 1972 .

[23]  Tomasz Janowski,et al.  Efficient Parallel Implementation of the CCSD External Exchange Operator and the Perturbative Triples (T) Energy Calculation. , 2008, Journal of chemical theory and computation.

[24]  A.P. Rendell,et al.  Enabling the Efficient Use of SMP Clusters: The GAMESS/DDI Model , 2003, ACM/IEEE SC 2003 Conference (SC'03).

[25]  Rodney J. Bartlett,et al.  Coupled-cluster methods with internal and semi-internal triply and quadruply excited clusters: CCSDt and CCSDtq approaches , 1999 .

[26]  Alistair P. Rendell,et al.  Quantum chemistry on parallel computer architectures: coupled-cluster theory applied to the bending potential of fulminic acid , 1992 .

[27]  Mark S. Gordon,et al.  A Novel Approach to Parallel Coupled Cluster Calculations: Combining Distributed and Shared Memory Techniques for Modern Cluster Based Systems , 2007 .

[28]  Jon Baker,et al.  Recent developments in the PQS program , 2012 .

[29]  Benny G. Johnson,et al.  Kohn—Sham density-functional theory within a finite basis set , 1992 .

[30]  Alistair P. Rendell,et al.  Evaluation of the contribution from triply excited intermediates to the fourth-order perturbation theory energy on Intel distributed memory supercomputers , 1993 .