English
 
User Manual Privacy Policy Disclaimer Contact us
  Advanced SearchBrowse

Item

ITEM ACTIONSEXPORT

Released

Journal Article

Efficient implementation of the analytic second derivatives of Hartree–Fock and hybrid DFT energies: a detailed analysis of different approximations

MPS-Authors
/persons/resource/persons237549

Bykov,  Dmytro
Research Department Neese, Max Planck Institute for Chemical Energy Conversion, Max Planck Society;

/persons/resource/persons237672

Petrenko,  Taras
Research Department Neese, Max Planck Institute for Chemical Energy Conversion, Max Planck Society;

/persons/resource/persons216815

Izsák,  Róbert
Research Department Neese, Max Planck Institute for Chemical Energy Conversion, Max Planck Society;

/persons/resource/persons237622

Kossmann,  Simone
Research Department Neese, Max Planck Institute for Chemical Energy Conversion, Max Planck Society;

/persons/resource/persons216802

Becker,  Ute
Research Department Neese, Max Planck Institute for Chemical Energy Conversion, Max Planck Society;

/persons/resource/persons216825

Neese,  Frank
Research Department Neese, Max Planck Institute for Chemical Energy Conversion, Max Planck Society;

External Ressource
No external resources are shared
Fulltext (public)
There are no public fulltexts stored in PuRe
Supplementary Material (public)
There is no public supplementary material available
Citation

Bykov, D., Petrenko, T., Izsák, R., Kossmann, S., Becker, U., Valeev, E., et al. (2015). Efficient implementation of the analytic second derivatives of Hartree–Fock and hybrid DFT energies: a detailed analysis of different approximations. Molecular Physics, 113(13-14), 1961-1977. doi:10.1080/00268976.2015.1025114.


Cite as: http://hdl.handle.net/21.11116/0000-0007-8973-D
Abstract
In this paper, various implementations of the analytic Hartree–Fock and hybrid density functional energy second derivatives are studied. An approximation-free four-centre implementation is presented, and its accuracy is rigorously analysed in terms of self-consistent field (SCF), coupled-perturbed SCF (CP-SCF) convergence and prescreening criteria. The CP-SCF residual norm convergence threshold turns out to be the most important of these. Final choices of convergence thresholds are made such that an accuracy of the vibrational frequencies of better than 5 cm−1 compared to the numerical noise-free results is obtained, even for the highly sensitive low frequencies (<100–200 cm−1). The effects of the choice of numerical grid for density functional exchange–correlation integrations are studied and various weight derivative schemes are analysed in detail. In the second step of the work, approximations are introduced in order to speed up the computation without compromising its accuracy. To this end, the accuracy and efficiency of the resolution of identity approximation for the Coulomb terms and the semi-numerical chain of spheres approximation to the exchange terms are carefully analysed. It is shown that the largest performance improvements are realised if either Hartree–Fock exchange is absent (pure density functionals) and otherwise, if the exchange terms in the CP-SCF step of the calculation are approximated by the COSX method in conjunction with a small integration grid. Default values for all the involved truncation parameters are suggested. For vancomycine (176 atoms and 3593 basis functions), the RIJCOSX Hessian calculation with the B3LYP functional and the def2-TZVP basis set takes ∼3 days using 16 Intel® Xeon® 2.60GHz processors with the COSX algorithm having a net parallelisation scaling of 11.9 which is at least ∼20 times faster than the calculation without the RIJCOSX approximation.