4.4 Article

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

Journal

MOLECULAR PHYSICS
Volume 113, Issue 13-14, Pages 1961-1977

Publisher

TAYLOR & FRANCIS LTD
DOI: 10.1080/00268976.2015.1025114

Keywords

COSX; Hartree-Fock; molecular Hessian; Density Functional Theory

Funding

  1. U.S. National Science Foundation [ACI-1047696]
  2. Max Planck Society
  3. Office of Advanced Cyberinfrastructure (OAC)
  4. Direct For Computer & Info Scie & Enginr [1047696] Funding Source: National Science Foundation

Ask authors/readers for more resources

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 similar to 3 days using 16 Intel (R) Xeon (R) 2.60GHz processors with the COSX algorithm having a net parallelisation scaling of 11.9 which is at least similar to 20 times faster than the calculation without the RIJCOSX approximation.

Authors

I am an author on this paper
Click your name to claim this paper and add it to your profile.

Reviews

Primary Rating

4.4
Not enough ratings

Secondary Ratings

Novelty
-
Significance
-
Scientific rigor
-
Rate this paper

Recommended

Article Chemistry, Multidisciplinary

Quantum computing in pharma: A multilayer embedding approach for near future applications

Robert Izsak, Christoph Riplinger, Nick S. Blunt, Bernardo de Souza, Nicole Holzmann, Ophelia Crawford, Joan Camps, Frank Neese, Patrick Schopf

Summary: Quantum computers are expected to be useful in simulating strongly correlated chemical systems, but careful selection of orbital spaces is crucial. This study presents a scheme for automatically selecting an active space and combines quantum phase estimation and variational quantum eigensolver algorithms to accurately describe the environment. The protocol outlined here can be applied to chemical systems of any size, including those beyond the capabilities of classical computation.

JOURNAL OF COMPUTATIONAL CHEMISTRY (2023)

Article Chemistry, Multidisciplinary

The SHARK integral generation and digestion system

Frank Neese

Summary: This paper describes the SHARK integral generation and digestion engine, which is based on the McMurchie/Davidson approach and utilizes an efficient BLAS algorithm. SHARK can handle various types of basis function integrals and features programming constructs that simplify workflows and avoid code duplication.

JOURNAL OF COMPUTATIONAL CHEMISTRY (2023)

Article Chemistry, Physical

Second quantisation for unrestricted references: formalism and quasi-spin-adaptation of excitation and spin-flip operators

Robert Izsak

Summary: A second quantized formulation for evaluating spin-dependent properties in the UHF basis is introduced in this article. The practical usage of this formulation is illustrated using the configuration interaction singles ansatz, and the results are compared to spin-restricted cases of the same ansatz. On a more tentative basis, the notion of quasi-spin-adaptation is also discussed, which ensures that the ansatz becomes spin-adapted as the unrestricted orbitals approach the restricted ones.

MOLECULAR PHYSICS (2023)

Editorial Material Chemistry, Physical

Comment on Canonical transcorrelated theory with projected Slater-type geminals [J. Chem. Phys. 136, 084107 (2012)]

Conner Masteran, Ashutosh Kumar, Nakul Teke, Bimal Gaudel, Takeshi Yanai, Edward F. Valeev

JOURNAL OF CHEMICAL PHYSICS (2023)

Article Chemistry, Physical

Efficient Construction of Canonical Polyadic Approximations of Tensor Networks

Karl Pierce, Edward F. Valeev

Summary: This paper addresses the problem of constructing a canonical polyadic (CP) decomposition for a tensor network rather than a single tensor. It demonstrates how leveraging the structure of the network during CP factor optimization can reduce the complexity of constructing an approximate CP representation. The utility of this technique is shown for approximating the order-4 Coulomb interaction tensor with two order-3 tensors using an approximate generalized square-root (SQ) factorization.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Chemistry, Physical

Massively parallel GPU enabled third-order cluster perturbation excitation energies for cost-effective large scale excitation energy calculations

Andreas Erbs Hillers-Bendtsen, Dmytro Bykov, Ashleigh Barnes, Dmitry Liakh, Hector H. Corzo, Jeppe Olsen, Poul Jorgensen, Kurt V. Mikkelsen

Summary: We have developed a massively parallel implementation of the CPS(D-3) excitation energy model based on cluster perturbation theory. The extended algorithm efficiently scales with increasing computational resources and allows for low time-to-solution calculations on large molecular systems. The CPS(D-3) excitation energies are shown to be a computationally efficient alternative to those obtained from the coupled-cluster singles and doubles model.

JOURNAL OF CHEMICAL PHYSICS (2023)

Article Chemistry, Physical

SparseMaps-A systematic infrastructure for reduced-scaling electronic structure methods. VI. Linear-scaling explicitly correlated N-electron valence state perturbation theory with pair natural orbital

Yang Guo, Fabijan Pavosevic, Kantharuban Sivalingam, Ute Becker, Edward F. F. Valeev, Frank Neese

Summary: In this work, a linear scaling explicitly correlated N-electron valence state perturbation theory (NEVPT2-F12) is developed to reduce the computational scaling of the conventional NEVPT2-F12 to near-linear scaling using the idea of a domain-based local pair natural orbital (DLPNO). The excitation energies predicted by DLPNO-NEVPT2-F12 are as accurate as the exact NEVPT2-F12 results for low-lying excited states of organic molecules. The new algorithm is applied to study cluster models of rhodopsin and can handle systems with large basis functions and active spaces.

JOURNAL OF CHEMICAL PHYSICS (2023)

Article Chemistry, Physical

Distributed memory, GPU accelerated Fock construction for hybrid, Gaussian basis density functional theory

David B. Williams-Young, Andrey Asadchev, Doru Thom Popovici, David Clark, Jonathan Waldrop, Theresa L. Windus, Edward F. Valeev, Wibe A. de Jong

Summary: With the increasing reliance on accelerator-based architecture in modern supercomputers, the development and optimization of electronic structure methods to utilize massively parallel resources has become a recent priority. This work presents distributed memory algorithms for evaluating the Coulomb and exact exchange matrices for hybrid Kohn-Sham DFT with Gaussian basis sets using direct density-fitted and seminumerical methods. The performance and scalability of these methods are demonstrated on systems ranging from a few hundred to over one thousand atoms using up to 128 NVIDIA A100 GPUs on the Perlmutter supercomputer.

JOURNAL OF CHEMICAL PHYSICS (2023)

Article Chemistry, Physical

Memory-Efficient Recursive Evaluation of 3-Center Gaussian Integrals

Andrey Asadchev, Edward F. Valeev

Summary: In order to improve the efficiency of Gaussian integral evaluation on modern accelerated architectures, FLOP-efficient Obara-Saika-based recursive evaluation schemes are optimized to reduce memory usage. The use of multiquantal recurrences is shown to save significant memory for evaluating 3-center 2-particle integrals. Other innovations include leveraging register memory and compile-time features of modern C++/CUDA. The performance of the proposed schemes is demonstrated for various types of integrals and the implementation is available in the open-source LibintX library.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Chemistry, Physical

Measuring Electron Correlation: The Impact of Symmetry and Orbital Transformations

Robert Izsak, Aleksei V. Ivanov, Nick S. Blunt, Nicole Holzmann, Frank Neese

Summary: In this article, the different measures of electron correlation in wave function theory, density functional theory, and quantum information theory are briefly reviewed. The focus is then placed on a more traditional metric based on dominant weights in the full configuration solution and its behavior regarding the choice of N-electron and one-electron basis. The impact of symmetry is discussed, and the usefulness of distinguishing between determinants, configuration state functions, and configurations as reference functions is emphasized.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Chemistry, Physical

Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules

Yuqi Wang, Yang Guo, Frank Neese, Edward F. Valeev, Wei Li, Shuhua Li

Summary: In this article, a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework are presented. These methods allow F12 calculations of large molecules on a single node. The authors used these methods to investigate the relative stability between extended and folded alkane, the stability of different structures of polyglycine, and the binding energies of host-guest complexes. The results demonstrate the promising potential of combining CIM with F12 methods for studying large molecules.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Chemistry, Physical

Direct Determination of Optimal Real-Space Orbitals for Correlated Electronic Structure of Molecules

Edward F. Valeev, Robert J. Harrison, Adam A. Holmes, Charles C. Peterson, Deborah A. Penchoff

Summary: Optimizing nearly exact orthonormal orbitals in real space using a multiresolution spectral element basis allows for evaluation of energy of arbitrary states of atoms and molecules. This method provides more accurate electronic energies compared to Gaussian AO bases, and is suitable for high-end models of correlated electronic states. The approach can optimize over 100 fully correlated numerical orbitals on a single computer node, with significant room for improvement.

JOURNAL OF CHEMICAL THEORY AND COMPUTATION (2023)

Article Chemistry, Multidisciplinary

Dimerization of confined Bronsted acids in enantioselective organocatalytic reactions

Ingolf Harden, Frank Neese, Giovanni Bistoni

Summary: This study investigates the aggregation process of chiral organocatalysts derived from phosphoric acids using high-level quantum mechanical calculations. The results show that the formation of catalyst dimers is possible and the aggregation process can significantly affect the stereo-controlling factors, reaction kinetics, and selectivity of the transformations.

CHEMICAL SCIENCE (2023)

Article Chemistry, Multidisciplinary

Dimerization of confined Brønsted acids in enantioselective organocatalytic reactions

Ingolf Harden, Frank Neese, Giovanni Bistoni

Summary: The formation of Bronsted acid aggregates in asymmetric organocatalytic reactions has a significant impact on the stereo-controlling factors of the transformations. This study uses high-level quantum mechanical calculations to investigate the influence of catalyst structure and reaction conditions on the spontaneity of the aggregation process, shedding light on its importance.

CHEMICAL SCIENCE (2023)

No Data Available