1
|
Graf D, Thom AJW. Corrected density functional theory and the random phase approximation: Improved accuracy at little extra cost. J Chem Phys 2023; 159:174106. [PMID: 37921249 DOI: 10.1063/5.0168569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2023] [Accepted: 10/16/2023] [Indexed: 11/04/2023] Open
Abstract
We recently introduced an efficient methodology to perform density-corrected Hartree-Fock density functional theory [DC(HF)-DFT] calculations and an extension to it we called "corrected" HF DFT [C(HF)-DFT] [Graf and Thom, J. Chem. Theory Comput. 19 5427-5438 (2023)]. In this work, we take a further step and combine C(HF)-DFT, augmented with a straightforward orbital energy correction, with the random phase approximation (RPA). We refer to the resulting methodology as corrected HF RPA [C(HF)-RPA]. We evaluate the proposed methodology across various RPA methods: direct RPA (dRPA), RPA with an approximate exchange kernel, and RPA with second-order screened exchange. C(HF)-dRPA demonstrates very promising performance; for RPA with exchange methods, on the other hand, we often find over-corrections.
Collapse
Affiliation(s)
- Daniel Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| | - Alex J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, England
| |
Collapse
|
2
|
Hermann J, Stöhr M, Góger S, Chaudhuri S, Aradi B, Maurer RJ, Tkatchenko A. libMBD: A general-purpose package for scalable quantum many-body dispersion calculations. J Chem Phys 2023; 159:174802. [PMID: 37933783 DOI: 10.1063/5.0170972] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2023] [Accepted: 10/17/2023] [Indexed: 11/08/2023] Open
Abstract
Many-body dispersion (MBD) is a powerful framework to treat van der Waals (vdW) dispersion interactions in density-functional theory and related atomistic modeling methods. Several independent implementations of MBD with varying degree of functionality exist across a number of electronic structure codes, which both limits the current users of those codes and complicates dissemination of new variants of MBD. Here, we develop and document libMBD, a library implementation of MBD that is functionally complete, efficient, easy to integrate with any electronic structure code, and already integrated in FHI-aims, DFTB+, VASP, Q-Chem, CASTEP, and Quantum ESPRESSO. libMBD is written in modern Fortran with bindings to C and Python, uses MPI/ScaLAPACK for parallelization, and implements MBD for both finite and periodic systems, with analytical gradients with respect to all input parameters. The computational cost has asymptotic cubic scaling with system size, and evaluation of gradients only changes the prefactor of the scaling law, with libMBD exhibiting strong scaling up to 256 processor cores. Other MBD properties beyond energy and gradients can be calculated with libMBD, such as the charge-density polarization, first-order Coulomb correction, the dielectric function, or the order-by-order expansion of the energy in the dipole interaction. Calculations on supramolecular complexes with MBD-corrected electronic structure methods and a meta-review of previous applications of MBD demonstrate the broad applicability of the libMBD package to treat vdW interactions.
Collapse
Affiliation(s)
- Jan Hermann
- Department of Mathematics and Computer Science, FU Berlin, 14195 Berlin, Germany
| | - Martin Stöhr
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Szabolcs Góger
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Shayantan Chaudhuri
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen, 28359 Bremen, Germany
| | - Reinhard J Maurer
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
3
|
A Cost Effective Scheme for the Highly Accurate Description of Intermolecular Binding in Large Complexes. Int J Mol Sci 2022; 23:ijms232415773. [PMID: 36555413 PMCID: PMC9780852 DOI: 10.3390/ijms232415773] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2022] [Revised: 11/23/2022] [Accepted: 12/07/2022] [Indexed: 12/15/2022] Open
Abstract
There has been a growing interest in quantitative predictions of the intermolecular binding energy of large complexes. One of the most important quantum chemical techniques capable of such predictions is the domain-based local pair natural orbital (DLPNO) scheme for the coupled cluster theory with singles, doubles, and iterative triples [CCSD(T)], whose results are extrapolated to the complete basis set (CBS) limit. Here, the DLPNO-based focal-point method is devised with the aim of obtaining CBS-extrapolated values that are very close to their canonical CCSD(T)/CBS counterparts, and thus may serve for routinely checking a performance of less expensive computational methods, for example, those based on the density-functional theory (DFT). The efficacy of this method is demonstrated for several sets of noncovalent complexes with varying amounts of the electrostatics, induction, and dispersion contributions to binding (as revealed by accurate DFT-based symmetry-adapted perturbation theory (SAPT) calculations). It is shown that when applied to dimeric models of poly(3-hydroxybutyrate) chains in its two polymorphic forms, the DLPNO-CCSD(T) and DFT-SAPT computational schemes agree to within about 2 kJ/mol of an absolute value of the interaction energy. These computational schemes thus should be useful for a reliable description of factors leading to the enthalpic stabilization of extended systems.
Collapse
|
4
|
Nagy PR, Gyevi-Nagy L, Lőrincz BD, Kállay M. Pursuing the basis set limit of CCSD(T) non-covalent interaction energies for medium-sized complexes: case study on the S66 compilation. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2109526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Péter R. Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - László Gyevi-Nagy
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Balázs D. Lőrincz
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| | - Mihály Kállay
- Faculty of Chemical Technology and Biotechnology, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, Budapest, Hungary
- ELKH-BME Quantum Chemistry Research Group, Budapest, Hungary
| |
Collapse
|
5
|
Villot C, Ballesteros F, Wang D, Lao KU. Coupled Cluster Benchmarking of Large Noncovalent Complexes in L7 and S12L as Well as the C 60 Dimer, DNA-Ellipticine, and HIV-Indinavir. J Phys Chem A 2022; 126:4326-4341. [PMID: 35766331 DOI: 10.1021/acs.jpca.2c01421] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
In this work, we report the benchmark binding energies of the seven complexes within the L7 data set, six host-guest complexes from the S12L data set, a C60 dimer, the DNA-ellipticine intercalation complex, and the largest system of the study, the HIV-indinavir system, which contained 343 atoms or 139 heavy atoms. The high-quality values reported were obtained via a focal point method that relies on the canonical form of second-order Møller-Plesset theory and the domain-based local pair natural orbital scheme for the coupled cluster with single double and perturbative triple excitations [DLPNO-CCSD(T)] extrapolated to the complete basis set (CBS) limit. The results in this work not only corroborate but also improve upon some previous benchmark values for large noncovalent complexes albeit at a relatively steep cost. Although local CCSD(T) and the largely successful fixed-node diffusion Monte Carlo (FN-DMC) have been shown to generally agree for small- to medium-size systems, a discrepancy in their reported binding energy values arises for large complexes, where the magnitude of the disagreement is a definite cause for concern. For example, the largest deviation in the L7 data set was 2.8 kcal/mol (∼10%) on the low end in C3GC. Such a deviation only grows worse in the S12L set, which showed a difference of up to 10.4 kcal/mol (∼25%) by a conservative estimation in buckycatcher-C60. The DNA-ellipticine complex also generated a disagreement of 4.4 kcal/mol (∼10%) between both state-of-the-art methods. The disagreement between local CCSD(T) and FN-DMC in large noncovalent complexes shows that it is urgently needed to have the canonical CCSD(T), the Monte Carlo CCSD(T), or the full configuration interaction quantum Monte Carlo approaches available to large systems on the hundred-atom scale to solve this dilemma. In addition, the performances of cheaper popular computational methods were assessed for the studied complexes with respect to DLPNO-CCSD(T)/CBS. r2SCAN-3c, B97M-V, and PBE0+D4 work well in large noncovalent complexes in this work, and GFN2-xTB performs well in π-π stacking complexes. B97M-V is the most reliable computationally efficient approach to predicting noncovalent interactions for large complexes, being the only one to have binding errors within the so-called 1 kcal/mol "chemical accuracy". The benchmark interaction energies of these host-guest complexes, molecular materials, and biological systems with electronic and medicinal implications provide crucial reference data for the improvement of current and future lower-cost methods.
Collapse
Affiliation(s)
- Corentin Villot
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Danyang Wang
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284 United States
| |
Collapse
|
6
|
Czernek J, Brus J, Czerneková V. A computational inspection of the dissociation energy of mid-sized organic dimers. J Chem Phys 2022; 156:204303. [DOI: 10.1063/5.0093557] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The gas-phase value of the dissociation energy ( D0) is a key parameter employed in both experimental and theoretical descriptions of noncovalent complexes. The D0 data were obtained for a set of mid-sized organic dimers in their global minima which was located using geometry optimizations that applied ample basis sets together with either the conventional second-order Møller–Plesset (MP2) method or several dispersion-corrected density-functional theory (DFT-D) schemes. The harmonic vibrational zero-point (VZP) and deformation energies from the MP2 calculations were combined with electronic energies from the coupled cluster theory with singles, doubles, and iterative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit to estimate D0 with the aim of inspecting values that were most recently measured, and an analogous comparison was performed using the DFT-D data. In at least one case (namely, for the aniline⋯methane cluster), the D0 estimate that employed the CCSD(T)/CBS energies differed from experiment in the way that could not be explained by a possible deficiency in the VZP contribution. Curiously, one of the DFT-D schemes (namely, the B3LYP-D3/def2-QZVPPD) was able to reproduce all measured D0 values to within 1.0 kJ/mol from experimental error bars. These findings show the need for further measurements and computations of some of the complexes. In order to facilitate such studies, the physical nature of intermolecular interactions in the investigated dimers was analyzed by means of the DFT-based symmetry-adapted perturbation theory.
Collapse
Affiliation(s)
- Jiří Czernek
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovsky Square 2, 162 06 Praha 6, The Czech Republic
| | - Jiří Brus
- Institute of Macromolecular Chemistry, Czech Academy of Sciences, Heyrovsky Square 2, 162 06 Praha 6, The Czech Republic
| | - Vladimíra Czerneková
- Institute of Physics, Czech Academy of Sciences, Na Slovance 2, 182 21 Praha 8, The Czech Republic
| |
Collapse
|
7
|
Gray M, Herbert JM. Comprehensive Basis-Set Testing of Extended Symmetry-Adapted Perturbation Theory and Assessment of Mixed-Basis Combinations to Reduce Cost. J Chem Theory Comput 2022; 18:2308-2330. [PMID: 35289608 DOI: 10.1021/acs.jctc.1c01302] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Hybrid or "extended" symmetry-adapted perturbation theory (XSAPT) replaces traditional SAPT's treatment of dispersion with better performing alternatives while at the same time extending two-body (dimer) SAPT to a many-body treatment of polarization using a self-consistent charge embedding procedure. The present work presents a systematic study of how XSAPT interaction energies and energy components converge with respect to the choice of Gaussian basis set. Errors can be reduced in a systematic way using correlation-consistent basis sets, with aug-cc-pVTZ results converged within <0.1 kcal/mol. Similar (if slightly less systematic) behavior is obtained using Karlsruhe basis sets at much lower cost, and we introduce new versions with limited augmentation that are even more efficient. Pople-style basis sets, which are more efficient still, often afford good results if a large number of polarization functions are included. The dispersion models used in XSAPT afford much faster basis-set convergence as compared to the perturbative description of dispersion in conventional SAPT, meaning that "compromise" basis sets (such as jun-cc-pVDZ) are no longer required and benchmark-quality results can be obtained using triple-ζ basis sets. The use of diffuse functions proves to be essential, especially for the description of hydrogen bonds. The "δ(Hartree-Fock)" correction for high-order induction can be performed in double-ζ basis sets without significant loss of accuracy, leading to a mixed-basis approach that offers 4× speedup over the existing (cubic scaling) XSAPT approach.
Collapse
Affiliation(s)
- Montgomery Gray
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
8
|
Shee J, Loipersberger M, Rettig A, Lee J, Head-Gordon M. Regularized Second-Order Møller-Plesset Theory: A More Accurate Alternative to Conventional MP2 for Noncovalent Interactions and Transition Metal Thermochemistry for the Same Computational Cost. J Phys Chem Lett 2021; 12:12084-12097. [PMID: 34910484 PMCID: PMC10037552 DOI: 10.1021/acs.jpclett.1c03468] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/14/2023]
Abstract
Second-order Møller-Plesset theory (MP2) notoriously breaks down for π-driven dispersion interactions and dative bonds in transition metal complexes. Herein, we investigate three physically justified forms of single-parameter, energy-gap dependent regularization which can yield high and transferable accuracy for a variety of noncovalent interactions (including S22, S66, and L7 test sets) and (mostly closed shell) transition metal thermochemistry. Regularization serves to damp overestimated pairwise additive contributions, renormalizing first-order amplitudes such that the effects of higher-order correlations are incorporated. The optimal parameter values for the noncovalent and transition metal sets are 1.1, 0.7, and 0.4 for κ, σ, and σ2 regularizers, respectively. However, such regularization slightly degrades the accuracy of conventional MP2 for some small-molecule test sets, most of which have relatively large average frontier energy gaps. Our results suggest that appropriately regularized MP2 models may improve double hybrid density functionals, at no additional cost over conventional MP2.
Collapse
Affiliation(s)
- James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthias Loipersberger
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
9
|
Carter-Fenk K, Lao KU, Herbert JM. Predicting and Understanding Non-Covalent Interactions Using Novel Forms of Symmetry-Adapted Perturbation Theory. Acc Chem Res 2021; 54:3679-3690. [PMID: 34550669 DOI: 10.1021/acs.accounts.1c00387] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Although sometimes derided as "weak" interactions, non-covalent forces play a critical role in ligand binding and crystal packing and in determining the conformational landscape of flexible molecules. Symmetry-adapted perturbation theory (SAPT) provides a framework for accurate ab initio calculation of intermolecular interactions and furnishes a natural decomposition of the interaction energy into physically meaningful components: semiclassical electrostatics (rigorously obtained from monomer charge densities), Pauli or steric repulsion, induction (including both polarization and charge transfer), and dispersion. This decomposition helps to foster deeper understanding of non-covalent interactions and can be used to construct transferable, physics-based force fields. Separability of the SAPT interaction energy also provides the flexibility to construct composite methods, a feature that we exploit to improve the description of dispersion interactions. These are challenging to describe accurately because they arise from nonlocal electron correlation effects that appear for the first time at second order in perturbation theory but are not quantitatively described at that level.As with all quantum-chemical methods, a major limitation of SAPT is nonlinear scaling of the computational cost with respect to system size. This cost can be significantly mitigated using "SAPT0(KS)", which incorporates monomer electron correlation by means of Kohn-Sham (KS) molecular orbitals from density functional theory (DFT), as well as by an "extended" theory called XSAPT, developed by the authors. XSAPT generalizes traditional dimer SAPT to many-body systems, so that a ligand-protein interaction (for example) can be separated into contributions from individual amino acids, reducing the cost of the calculation below that of even supramolecular DFT while retaining the accuracy of high-level ab initio quantum chemistry.This Account provides an overview of the SAPT0(KS) approach and the XSAPT family of methods. Several low-cost variants are described that provide accuracy approaching that of the best ab initio benchmarks yet are affordable enough to tackle ligand-protein binding and sizable host-guest complexes. These variants include SAPT+aiD, which uses ab initio atom-atom dispersion potentials ("+aiD") in place of second-order SAPT dispersion, and also SAPT+MBD, which incorporates many-body dispersion (MBD) effects that are important in the description of nanoscale materials. Applications to drug binding highlight the size-extensive nature of dispersion, which is not a weak interaction in large systems. Other applications highlight how a physics-based analysis can sometimes upend conventional wisdom regarding intermolecular forces. In particular, careful reconsideration of π-π interactions makes clear that the quadrupolar electrostatics (or "Hunter-Sanders") model of π-π stacking should be replaced by a "van der Waals model" in which conformational preferences arise from a competition between dispersion and Pauli repulsion. Our analysis also suggests that molecular shape, rather than aromaticity per se, is the key factor driving strong stacking interactions. Looking forward, we anticipate that XSAPT-based methods can play a role in screening of drug candidates and in materials design.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, United States
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
10
|
Doran AE, Qiu DL, Hirata S. Monte Carlo MP2-F12 for Noncovalent Interactions: The C 60 Dimer. J Phys Chem A 2021; 125:7344-7351. [PMID: 34433271 DOI: 10.1021/acs.jpca.1c05021] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A scalable stochastic algorithm is presented that can evaluate explicitly correlated (F12) second-order many-body perturbation (MP2) energies of weak, noncovalent, intermolecular interactions. It first transforms the formulas of the MP2 and F12 energy differences into a short sum of high-dimensional integrals of Green's functions in real space and imaginary time. These integrals are then evaluated by the Monte Carlo method augmented by parallel execution, redundant-walker convergence acceleration, direct-sampling autocorrelation elimination, and control-variate error reduction. By sharing electron-pair walkers across the supermolecule and its subsystems spanned by the joint basis set, the statistical uncertainty is reduced by one to 2 orders of magnitude in the MP2 binding energy corrected for the basis-set incompleteness and superposition errors. The method predicts the MP2-F12/aug-cc-pVDZ binding energy of 19.1 ± 4.0 kcal mol-1 for the C60 dimer at the center distance of 9.748 Å.
Collapse
Affiliation(s)
- Alexander E Doran
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - David L Qiu
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
11
|
Ballesteros F, Dunivan S, Lao KU. Coupled cluster benchmarks of large noncovalent complexes: The L7 dataset as well as DNA-ellipticine and buckycatcher-fullerene. J Chem Phys 2021; 154:154104. [PMID: 33887937 DOI: 10.1063/5.0042906] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
In this work, benchmark binding energies for dispersion-bound complexes in the L7 dataset, the DNA-ellipticine intercalation complex, and the buckycatcher-C60 complex with 120 heavy atoms using a focal-point method based on the canonical form of second-order Møller-Plesset theory (MP2) and the domain based local pair natural orbital scheme for the coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] extrapolated to the complete basis set (CBS) limit are reported. This work allows for increased confidence given the agreement with respect to values recently obtained using the local natural orbital CCSD(T) for L7 and the canonical CCSD(T)/CBS result for the coronene dimer (C2C2PD). Therefore, these results can be considered pushing the CCSD(T)/CBS binding benchmark to the hundred-atom scale. The disagreements between the two state-of-the-art methods, CCSD(T) and fixed-node diffusion Monte Carlo, are substantial with at least 2.0 (∼10%), 1.9 (∼5%), and 10.3 kcal/mol (∼25%) differences for C2C2PD in L7, DNA-ellipticine, and buckycatcher-C60, respectively. Such sizable discrepancy above "chemical accuracy" for large noncovalent complexes indicates how challenging it is to obtain benchmark binding interactions for systems beyond small molecules, although the three up-to-date density functionals, PBE0+D4, ωB97M-V, and B97M-V, agree better with CCSD(T) for these large systems. In addition to reporting these values, different basis sets and various CBS extrapolation parameters for Hartree-Fock and MP2 correlation energies were tested for the first time in large noncovalent complexes with the goal of providing some indications toward optimal cost effective routes to approach the CBS limit without substantial loss in quality.
Collapse
Affiliation(s)
- Francisco Ballesteros
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| | - Shelbie Dunivan
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| | - Ka Un Lao
- Department of Chemistry, Virginia Commonwealth University, Richmond, Virginia 23284, USA
| |
Collapse
|
12
|
Interactions between large molecules pose a puzzle for reference quantum mechanical methods. Nat Commun 2021; 12:3927. [PMID: 34168142 PMCID: PMC8225865 DOI: 10.1038/s41467-021-24119-3] [Citation(s) in RCA: 45] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2020] [Accepted: 06/02/2021] [Indexed: 02/05/2023] Open
Abstract
Quantum-mechanical methods are used for understanding molecular interactions throughout the natural sciences. Quantum diffusion Monte Carlo (DMC) and coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] are state-of-the-art trusted wavefunction methods that have been shown to yield accurate interaction energies for small organic molecules. These methods provide valuable reference information for widely-used semi-empirical and machine learning potentials, especially where experimental information is scarce. However, agreement for systems beyond small molecules is a crucial remaining milestone for cementing the benchmark accuracy of these methods. We show that CCSD(T) and DMC interaction energies are not consistent for a set of polarizable supramolecules. Whilst there is agreement for some of the complexes, in a few key systems disagreements of up to 8 kcal mol-1 remain. These findings thus indicate that more caution is required when aiming at reproducible non-covalent interactions between extended molecules.
Collapse
|
13
|
Morales-Silva MA, Jordan KD, Shulenburger L, Wagner LK. Frontiers of stochastic electronic structure calculations. J Chem Phys 2021; 154:170401. [PMID: 34241059 DOI: 10.1063/5.0053674] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
In recent years there has been a rapid growth in the development and application of new stochastic methods in electronic structure. These methods are quite diverse, from many-body wave function techniques in real space or determinant space to being used to sum perturbative expansions. This growth has been spurred by the more favorable scaling with the number of electrons and often better parallelization over large numbers of central processing unit (CPU) cores or graphical processing units (GPUs) than for high-end non-stochastic wave function based methods. This special issue of the Journal of Chemical Physics includes 33 papers that describe recent developments and applications in this area. As seen from the articles in the issue, stochastic electronic structure methods are applicable to both molecules and solids and can accurately describe systems with strong electron correlation. This issue was motivated, in part, by the 2019 Telluride Science Research Center workshop on Stochastic Electronic Structure Methods that we organized. Below we briefly describe each of the papers in the special issue, dividing the papers into six subtopics.
Collapse
Affiliation(s)
- Miguel A Morales-Silva
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94551, USA
| | - Kenneth D Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Luke Shulenburger
- HEDP Theory Department, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Lucas K Wagner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|