Assessing the domain-based local pair natural orbital (DLPNO) approximation for non-covalent interactions in sizable supramolecular complexes

13 May 2024, Version 2
This content is a preprint and has not undergone peer review at the time of posting.

Abstract

The titular DLPNO approximation {\red is the most widely used method} for extending correlated wave function models to large molecular systems, yet its fidelity for intermolecular interaction energies has not been thoroughly vetted. Non-covalent interactions are sensitive to tails of the electron density, involving parts of the wave function that are far from the nuclei and may be discarded in some local correlation treatments. Meanwhile, the accuracy of the DLPNO approximation is known to deteriorate as molecular size increases, and questions have been raised regarding the accuracy of benchmark calculations for large van der Waals complexes. We test the DLPNO approximation at the level of second-order Moller-Plesset perturbation theory (MP2) in systems containing as many as 240 atoms, for which canonical MP2 calculations can be performed for comparison. For small dimers with <~ 10 heavy atoms, DLPNO-MP2 interaction energies are within 3% of canonical values. However, the approximation is quite poor for larger systems unless the results are extrapolated to the limit where the threshold for discarding PNOs is taken to zero. For a sequence of nanoscale graphene dimers up to (C96H24)2, extrapolated DLPNO-MP2 interaction energies agree with canonical values to within 1%, independent of system size, provided that the basis set does not contain diffuse functions. The presence of diffuse basis functions causes oscillatory behavior as a function of the PNO threshold, making it impossible to extrapolate the results in a meaningful way. The significance of these results for higher-level coupled-cluster benchmarks is also addressed, via calculations in small dimers.

Keywords

local correlation
dispersion
van der Waals
benchmarks
quantum chemistry

Supplementary materials

Title
Description
Actions
Title
Additional data
Description
A significant amount of additional testing and analysis.
Actions
Title
Geometries
Description
Cartesian coordinates for test systems created as part of this work.
Actions
Title
Raw data
Description
Single-point energies for all systems studied
Actions

Comments

Comments are not moderated before they are posted, but they can be removed by the site moderators if they are found to be in contravention of our Commenting Policy [opens in a new tab] - please read this policy before you post. Comments should be used for scholarly discussion of the content in question. You can find more information about how to use the commenting feature here [opens in a new tab] .
This site is protected by reCAPTCHA and the Google Privacy Policy [opens in a new tab] and Terms of Service [opens in a new tab] apply.