skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Extending density functional theory with near chemical accuracy beyond pure water
Abstract Density functional simulations of condensed phase water are typically inaccurate, due to the inaccuracies of approximate functionals. A recent breakthrough showed that the SCAN approximation can yield chemical accuracy for pure water in all its phases, but only when its density is corrected. This is a crucial step toward first-principles biosimulations. However, weak dispersion forces are ubiquitous and play a key role in noncovalent interactions among biomolecules, but are not included in the new approach. Moreover, naïve inclusion of dispersion in HF-SCAN ruins its high accuracy for pure water. Here we show that systematic application of the principles of density-corrected DFT yields a functional (HF-r 2 SCAN-DC4) which recovers and not only improves over HF-SCAN for pure water, but also captures vital noncovalent interactions in biomolecules, making it suitable for simulations of solutions.  more » « less
Award ID(s):
2154371
PAR ID:
10396687
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Nature Communications
Volume:
14
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; (, The journal of physical chemistry letters)
    Empirical fitting of parameters in approximate density functionals is common. Such fits conflate errors in the self-consistent density with errors in the energy functional, but density-corrected DFT (DC-DFT) separates these two. We illustrate with catastrophic failures of a toy functional applied to H2+ at varying bond lengths, where the standard fitting procedure misses the exact functional; Grimme’s D3 fit to noncovalent interactions, which can be contaminated by large density errors such as in the WATER27 and B30 data sets; and double-hybrids trained on self-consistent densities, which can perform poorly on systems with density-driven errors. In these cases, more accurate results are found at no additional cost by using Hartree–Fock (HF) densities instead of self-consistent densities. For binding energies of small water clusters, errors are greatly reduced. Range-separated hybrids with 100% HF at large distances suffer much less from this effect. 
    more » « less
  2. ; ; (, Physical Chemistry Chemical Physics)
    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
  3. ; ; ; (, Chemical Physics Reviews)
    Density functional theory (DFT) has been applied to modeling molecular interactions in water for over three decades. The ubiquity of water in chemical and biological processes demands a unified understanding of its physics, from the single molecule to the thermodynamic limit and everything in between. Recent advances in the development of data-driven and machine-learning potentials have accelerated simulation of water and aqueous systems with DFT accuracy. However, anomalous properties of water in the condensed phase, where a rigorous treatment of both local and non-local many-body (MB) interactions is in order, are often unsatisfactory or partially missing in DFT models of water. In this review, we discuss the modeling of water and aqueous systems based on DFT and provide a comprehensive description of a general theoretical/computational framework for the development of data-driven many-body potentials from DFT reference data. This framework, coined MB-DFT, readily enables efficient many-body molecular dynamics (MD) simulations of small molecules, in both gas and condensed phases, while preserving the accuracy of the underlying DFT model. Theoretical considerations are emphasized, including the role that the delocalization error plays in MB-DFT potentials of water and the possibility to elevate DFT and MB-DFT to near-chemical-accuracy through a density-corrected formalism. The development of the MB-DFT framework is described in detail, along with its application in MB-MD simulations and recent extension to the modeling of reactive processes in solution within a quantum mechanics/MB molecular mechanics (QM/MB-MM) scheme, using water as a prototypical solvent. Finally, we identify open challenges and discuss future directions for MB-DFT and QM/MB-MM simulations in condensed phases. 
    more » « less
  4. ; ; (, International Journal of Quantum Chemistry)
    Abstract We present a graph‐theoretic approach to adaptively compute many‐body approximations in an efficient manner to perform (a) accurate post‐Hartree–Fock (HF) ab initio molecular dynamics (AIMD) at density functional theory (DFT) cost for medium‐ to large‐sized molecular clusters, (b) hybrid DFT electronic structure calculations for condensed‐phase simulations at the cost of pure density functionals, (c) reduced‐cost on‐the‐fly basis extrapolation for gas‐phase AIMD and condensed phase studies, and (d) accurate post‐HF‐level potential energy surfaces at DFT cost for quantum nuclear effects. The salient features of our approach are ONIOM‐like in that (a) the full system (cluster or condensed phase) calculation is performed at a lower level of theory (pure DFT for condensed phase or hybrid DFT for molecular systems), and (b) this approximation is improved through a correction term that captures all many‐body interactions up to any given order within a higher level of theory (hybrid DFT for condensed phase; CCSD or MP2 for cluster), combined through graph‐theoretic methods. Specifically, a region of chemical interest is coarse‐grained into a set of nodes and these nodes are then connected to form edges based on a given definition of local envelope (or threshold) of interactions. The nodes and edges together define a graph, which forms the basis for developing the many‐body expansion. The methods are demonstrated through (a) ab initio dynamics studies on protonated water clusters and polypeptide fragments, (b) potential energy surface calculations on one‐dimensional water chains such as those found in ion channels, and (c) conformational stabilization and lattice energy studies on homogeneous and heterogeneous surfaces of water with organic adsorbates using two‐dimensional periodic boundary conditions. 
    more » « less
  5. ; (, Physical Chemistry Chemical Physics)
    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
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email