skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Speeding up equation of motion coupled cluster theory with the chain of spheres approximation

Journal Article · · Journal of Chemical Physics
DOI:https://doi.org/10.1063/1.4939844· OSTI ID:22493654
; ;  [1]
  1. Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr (Germany)

In the present paper, the chain of spheres exchange (COSX) approximation is applied to the highest scaling terms in the equation of motion (EOM) coupled cluster equations with single and double excitations, in particular, the terms involving integrals with four virtual labels. It is found that even the acceleration of this single term yields significant computational gains without compromising the desired accuracy of the method. For an excitation energy calculation on a cluster of five water molecules using 585 basis functions, the four virtual term is 9.4 times faster using COSX with a loose grid than using the canonical implementation, which yields a 2.6 fold acceleration for the whole of the EOM calculation. For electron attachment calculations, the four virtual term is 15 times and the total EOM calculation is 10 times faster than the canonical calculation for the same system. The accuracy of the new method was tested using Thiel’s test set for excited states using the same settings and the maximum absolute deviation over the whole test set was found to be 12.945 cm{sup −1} (59 μHartree) for excitation energies and 6.799 cm{sup −1} (31 μHartree) for electron attachments. Using MP2 amplitudes for the ground state in combination with the parallel evaluation of the full EOM equations in the manner discussed in this paper enabled us to perform calculations for large systems. Electron affinity values for the two lowest states of a Zn protoporphyrine model compound (224 correlated electrons and 1120 basis functions) were obtained in 3 days 19 h using 4 cores of a Xeon E5-2670 processor allocating 10 GB memory per core. Calculating the lowest two excitation energies for trans-retinal (114 correlated electrons and 539 basis functions) took 1 day 21 h using eight cores of the same processor and identical memory allocation per core.

OSTI ID:
22493654
Journal Information:
Journal of Chemical Physics, Vol. 144, Issue 3; Other Information: (c) 2016 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA); ISSN 0021-9606
Country of Publication:
United States
Language:
English