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

被引:2222
作者
Neese, Frank [1 ,2 ]
Wennmohs, Frank [1 ]
Hansen, Andreas [1 ]
Becker, Ute [2 ]
机构
[1] Univ Bonn, Inst Phys & Theoret Chem, Lehrstuhl Theoret Chem, D-53115 Bonn, Germany
[2] Max Planck Inst Bioanorgan Chem, D-45470 Mulheim An Der Ruhr, Germany
关键词
Hartree-Fock; Hybrid DFT; Density fitting; Numerical integration; Parallel electronic structure; ZETA-VALENCE QUALITY; AUXILIARY BASIS-SETS; PLESSET PERTURBATION-THEORY; GAUSSIAN-BASIS SETS; ELECTRONIC-STRUCTURE; DIRECT SCF; IDENTITY APPROXIMATION; LARGE MOLECULES; DENSITY; RESOLUTION;
D O I
10.1016/j.chemphys.2008.10.036
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
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 degrees 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. (C) 2008 Elsevier B.V. All rights reserved.
引用
收藏
页码:98 / 109
页数:12
相关论文
共 81 条
[41]   IMPROVEMENTS ON THE DIRECT SCF METHOD [J].
HASER, M ;
AHLRICHS, R .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1989, 10 (01) :104-111
[42]  
Helgaker T., 1995, MODERN ELECT STRUCTU, P725, DOI DOI 10.1016/J.CCR.2014.12.016
[43]   PRESCREENING OF 2-ELECTRON INTEGRAL DERIVATIVES IN SCF GRADIENT AND HESSIAN CALCULATIONS [J].
HORN, H ;
WEISS, H ;
HASER, M ;
EHRIG, M ;
AHLRICHS, R .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1991, 12 (09) :1058-1064
[44]   The impact of the resolution of the identity approximate integral method on modern ab initio algorithm development [J].
Kendall, RA ;
Fruchtl, HA .
THEORETICAL CHEMISTRY ACCOUNTS, 1997, 97 (1-4) :158-163
[45]   Pseudospectral time-dependent density functional theory [J].
Ko, Chaehyuk ;
Malick, David K. ;
Braden, Dale A. ;
Friesner, Richard A. ;
Martinez, Todd J. .
JOURNAL OF CHEMICAL PHYSICS, 2008, 128 (10)
[46]   An adaptive numerical integrator for molecular integrals [J].
Krack, M ;
Köster, AM .
JOURNAL OF CHEMICAL PHYSICS, 1998, 108 (08) :3226-3234
[47]   Range definitions for Gaussian-type charge distributions in fast multipole methods [J].
Kudin, KN ;
Scuseria, GE .
JOURNAL OF CHEMICAL PHYSICS, 1999, 111 (06) :2351-2356
[48]  
LAMBRECHT DS, 2005, J CHEM PHYS, V123
[49]   Fast evaluation of a linear number of local exchange matrices [J].
Liang, WZ ;
Shao, YH ;
Ochsenfeld, C ;
Bell, AT ;
Head-Gordon, M .
CHEMICAL PHYSICS LETTERS, 2002, 358 (1-2) :43-50
[50]   ONE-ELECTRON AND 2-ELECTRON INTEGRALS OVER CARTESIAN GAUSSIAN FUNCTIONS [J].
MCMURCHIE, LE ;
DAVIDSON, ER .
JOURNAL OF COMPUTATIONAL PHYSICS, 1978, 26 (02) :218-231