An official website of the United States government
We report the implementation of a symmetry-adapted perturbation theory algorithm based on a density functional theory [SAPT(DFT)] description of monomers. The implementation adopts a density-fitting treatment of hybrid exchange–correlation kernels to enable the description of monomers with hybrid functionals, as in the algorithm by Bukowski, Podeszwa, and Szalewicz [Chem. Phys. Lett. 414, 111 (2005)]. We have improved the algorithm by increasing numerical stability with QR factorization and optimized the computation of the exchange–correlation kernel with its 2-index density-fitted representation. The algorithm scales as O( N 5 ) formally and is usable for systems with up to ∼3000 basis functions, as demonstrated for the C 60 –buckycatcher complex with the aug-cc-pVDZ basis set. The hybrid-kernel-based SAPT(DFT) algorithm is shown to be as accurate as SAPT(DFT) implementations based on local effective exact exchange potentials obtained from the local Hartree–Fock (LHF) method while avoiding the lower-scaling [ O( N 4 )] but iterative and sometimes hard-to-converge LHF process. The hybrid-kernel algorithm outperforms Hartree–Fock-based SAPT (SAPT0) for the S66 test set, and its accuracy is comparable to the many-body perturbation theory based SAPT2+ approach, which scales as O( N 7 ), although SAPT2+ exhibits a more narrow distribution of errors.
more »
« less
Optimization of damping function parameters for -D3 and -D4 dispersion models for Hartree–Fock based symmetry-adapted perturbation theory
Symmetry-adapted perturbation theory (SAPT) directly computes intermolecular interaction energy in terms of electrostatics, exchange-repulsion, induction/polarization, and London dispersion components. In SAPT based on Hartree–Fock (“SAPT0”) or based on density functional theory, the most time-consuming step is the computation of the dispersion terms. Previous work has explored the replacement of these expensive dispersion terms with simple damped asymptotic models. We recently examined [Schriber et al. J. Chem. Phys. 154, 234107 (2021)] the accuracy of SAPT0 when replacing its dispersion term with Grimme’s popular -D3 correction, reducing the computational cost scaling from O(N5) to O(N3). That work optimized damping function parameters for SAPT0-D3/jun-cc-pVDZ using estimates of the coupled-cluster complete basis set limit [CCSD(T)/CBS] on a 8299 dimer dataset. Here, we explore the accuracy of SAPT0-D3 with additional basis sets, along with an analogous model using -D4. Damping parameters are rather insensitive to basis sets, and the resulting SAPT0-D models are more accurate on average for total interaction energies than SAPT0. Our results are surprising in several respects: (1) improvement of -D4 over -D3 is negligible for these systems, even charged systems where -D4 should, in principle, be more accurate; (2) addition of Axilrod–Teller–Muto terms for three-body dispersion does not improve error statistics for this test set; and (3) SAPT0-D is even more accurate on average for total interaction energies than the much more computationally costly density functional theory based SAPT [SAPT(DFT)] in an aug-cc-pVDZ basis. However, SAPT0 and SAPT0-D3/D4 interaction energies benefit from significant error cancellation between exchange and dispersion terms.
more »
« less
- Award ID(s):
- 1955940
- PAR ID:
- 10566311
- Publisher / Repository:
- American Institute of Physics Publishing
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 161
- Issue:
- 11
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 114115
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
An impressive collection of accurate two-body interaction energies for small complexes has been assembled into benchmark databases and used to improve the performance of multiple density functional, semiempirical, and machine learning methods. Similar benchmark data on nonadditive three-body energies in molecular trimers are comparatively scarce, and the existing ones are practically limited to homotrimers. In this work, we present a benchmark dataset of 20 equilibrium noncovalent interaction energies for a small but diverse selection of 10 heteromolecular trimers. The new 3BHET dataset presents complexes that combine different interactions including π −π, anion−π, cation−π, and various motifs of hydrogen and halogen bonding in each trimer. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition of the two- and three-body interaction energies shows that 3BHET consists of electrostatics- and dispersion-dominated complexes. The nonadditive three-body contribution is dominated by induction, but its influence on the overall bonding type in the complex (as exemplified by its position on the ternary diagram) is quite small. We also tested the extended SAPT (XSAPT) approach which is capable of including some nonadditive interactions in clusters of any size. The resulting three-body dispersion term (obtained from the many-body dispersion formalism) is mostly in good agreement with the supermolecular CCSD(T)−MP2 values and the nonadditive induction term is similar to the three-body SAPT(DFT) data, but the overall three-body XSAPT energies are not very accurate as they are missing the first-order exchange terms.more » « less
-
Coupled cluster theory has had a momentous impact on the ab initio prediction of molecular properties, and remains a staple ingratiate in high-accuracy thermochemical model chemistries. However, these methods require inclusion of at least some connected quadruple excitations, which generally scale at best as 𝒪(𝑁9) with the number of basis functions. It is very difficult to predict, a priori, the effect correlation of past CCSD(T) on a given reaction energy. The purpose of this work is to examine cost-effective quadruple corrections based on the factorization theorem of the many-body perturbation theory that may address these challenges. We show that the 𝒪(𝑁7) factorized CCSD(TQf) method introduces minimal error to predicted correlation and reaction energies as compared to the 𝒪(𝑁9) CCSD(TQ). Further, we examine the performance of Goodson’s continued fraction method in the estimation of CCSDT(Q)Λ contributions to reaction energies as well as a “new” method related to %TAE[(T)] that we refer to as a scaled perturbation estimator. We find that the scaled perturbation estimator based upon CCSD(TQf)/cc-pVDZ is capable of predicting CCSDT(Q)Λ/cc-pVDZ contributions to reaction energies with an average error of 0.07 kcal mol–1 and an L2D of 0.52 kcal mol–1 when applied to a test-suite of nearly 3000 reactions. This offers a means by which to reliably “ballpark” how important post-CCSD(T) contributions are to reaction energies while incurring no more than CCSD(T) formal cost and a little mental math.more » « less
-
Second-order Møller–Plesset perturbation theory (MP2) provides a valuable alternative to density functional theory for modeling problems in organic and biological chemistry. However, MP2 suffers from known limitations in the description of van der Waals (London) dispersion interactions and reaction thermochemistry. Here, a spin-component-scaled, dispersion-corrected MP2 model (SCS-MP2D) is proposed that addresses these weaknesses. The dispersion correction, which is based on Grimme's D3 formalism, replaces the uncoupled Hartree–Fock dispersion inherent in MP2 with a more robust coupled Kohn–Sham treatment. The spin-component scaling of the residual MP2 correlation energy then reduces the remaining errors in the model. This two-part correction strategy solves the problem found in earlier spin-component-scaled MP2 models where completely different spin-scaling parameters were needed for describing reaction energies versus intermolecular interactions. Results on 18 benchmark data sets and two challenging potential energy curves demonstrate that SCS-MP2D considerably improves upon the accuracy of MP2 for intermolecular interactions, conformational energies, and reaction energies. Its accuracy and computational cost are competitive with state-of-the-art density functionals such as DSD-BLYP-D3(BJ), revDSD-PBEP86-D3(BJ), ωB97X-V, and ωB97M-V for systems with ∼100 atoms.more » « less
-
Symmetry-adapted perturbation theory (SAPT) is an ab initio approach that directly computes noncovalent interaction energies in terms of electrostatics, exchange repulsion, induction/polarization, and London dispersion components. Due to its high computational scaling, routine applications of even the lowest order of SAPT are typically limited to a few hundred atoms. To address this limitation, we report here the addition of electrostatic embedding to the SAPT (EE-SAPT) and ISAPT (EE-ISAPT) methods. We illustrate the embedding scheme using water trimer as a prototype example. Then, we show that EE-SAPT/EE-ISAPT can be applied for efficiently and accurately computing noncovalent interactions in large systems, including solvated dimers and protein–ligand systems. In the latter application, particular care must be taken to properly handle the quantum mechanics/molecular mechanics boundary when it cuts covalent bonds. We investigate various schemes for handling charges near this boundary and demonstrate which are most effective in the context of charge-embedded SAPT.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0219185
