English
 
Help Privacy Policy Disclaimer
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange

MPS-Authors
/persons/resource/persons216825

Neese,  Frank
Lehrstuhl für Theoretische Chemie, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Wegelerstrasse 12, D-53115 Bonn, Germany;
Research Department Neese, Max Planck Institute for Bioinorganic Chemistry, Max Planck Society;

/persons/resource/persons216802

Becker,  Ute
Research Department Wieghardt, Max Planck Institute for Bioinorganic Chemistry, Max Planck Society;

External Resource
No external resources are shared
Fulltext (restricted access)
There are currently no full texts shared for your IP range.
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available
Citation

Neese, F., Wennmohs, F., Hansen, A., & Becker, U. (2009). Efficient, approximate and parallel Hartree–Fock and hybrid DFT calculations. A ‘chain-of-spheres’ algorithm for the Hartree–Fock exchange. Chemical Physics, 356(1-3), 98-109. doi:10.1016/j.chemphys.2008.10.036.


Cite as: https://hdl.handle.net/21.11116/0000-0008-3344-2
Abstract
In this paper, the possibility is explored to speed up Hartree–Fock and hybrid density functional calculations by forming the Coulomb and exchange parts of the Fock matrix by different approximations. For the Coulomb part the previously introduced Split-RI-J variant (F. Neese, J. Comput. Chem. 24 (2003) 1740) of the well-known ‘density fitting’ approximation is used. The exchange part is formed by semi-numerical integration techniques that are closely related to Friesner’s pioneering pseudo-spectral approach. Our potentially linear scaling realization of this algorithm is called the ‘chain-of-spheres exchange’ (COSX). A combination of semi-numerical integration and density fitting is also proposed. Both Split-RI-J and COSX scale very well with the highest angular momentum in the basis sets. It is shown that for extended basis sets speed-ups of up to two orders of magnitude compared to traditional implementations can be obtained in this way. Total energies are reproduced with an average error of <0.3 kcal/mol as determined from extended test calculations with various basis sets on a set of 26 molecules with 20–200 atoms and up to 2000 basis functions. Reaction energies agree to within 0.2 kcal/mol (Hartree–Fock) or 0.05 kcal/mol (hybrid DFT) with the canonical values. The COSX algorithm parallelizes with a speedup of 8.6 observed for 10 processes. Minimum energy geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angels from their canonical values. These developments enable highly efficient and accurate self-consistent field calculations including nonlocal Hartree–Fock exchange for large molecules. In combination with the RI-MP2 method and large basis sets, second-order many body perturbation energies can be obtained for medium sized molecules with unprecedented efficiency. The algorithms are implemented into the ORCA electronic structure system.