skip to main content


Title: Quantum chemical accuracy from density functional approximations via machine learning
Abstract

Kohn-Sham density functional theory (DFT) is a standard tool in most branches of chemistry, but accuracies for many molecules are limited to 2-3 kcal ⋅ mol−1with presently-available functionals. Ab initio methods, such as coupled-cluster, routinely produce much higher accuracy, but computational costs limit their application to small molecules. In this paper, we leverage machine learning to calculate coupled-cluster energies from DFT densities, reaching quantum chemical accuracy (errors below 1 kcal ⋅ mol−1) on test data. Moreover, density-basedΔ-learning (learning only the correction to a standard DFT calculation, termedΔ-DFT ) significantly reduces the amount of training data required, particularly when molecular symmetries are included. The robustness ofΔ-DFT  is highlighted by correcting “on the fly” DFT-based molecular dynamics (MD) simulations of resorcinol (C6H4(OH)2) to obtain MD trajectories with coupled-cluster accuracy. We conclude, therefore, thatΔ-DFT  facilitates running gas-phase MD simulations with quantum chemical accuracy, even for strained geometries and conformer changes where standard DFT fails.

 
more » « less
Award ID(s):
1856165
PAR ID:
10305325
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Volume:
11
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. With dual goals of efficient and accurate modeling of solvation thermodynamics in molten salt liquids, we employ ab initio molecular dynamics (AIMD) simulations, deep neural network interatomic potentials (NNIP), and quasichemical theory (QCT) to calculate the excess chemical potentials for the solute ions Na + and Cl − in the molten NaCl liquid. NNIP-based molecular dynamics simulations accelerate the calculations by 3 orders of magnitude and reduce the uncertainty to 1 kcal mol −1 . Using the Density Functional Theory (DFT) level of theory, the predicted excess chemical potential for the solute ion pair is −178.5 ± 1.1 kcal mol −1 . A quantum correction of 13.7 ± 1.9 kcal mol −1 is estimated via higher-level quantum chemistry calculations, leading to a final predicted ion pair excess chemical potential of −164.8 ± 2.2 kcal mol −1 . The result is in good agreement with a value of −163.5 kcal mol −1 obtained from thermo-chemical tables. This study validates the application of QCT and NNIP simulations to the molten salt liquids, allowing for significant insights into the solvation thermodynamics crucial for numerous molten salt applications. 
    more » « less
  2. Rational design of molecular chelating agents requires a detailed understanding of physicochemical ligand–metal interactions in solvent phase. Computational quantum chemistry methods should be able to provide this, but computational reports have shown poor accuracy when determining absolute binding constants for many chelating molecules. To understand why, we compare and benchmark static- and dynamics-based computational procedures for a range of monovalent and divalent cations binding to a conventional cryptand molecule: 2.2.2-cryptand ([2.2.2]). The benchmarking comparison shows that dynamics simulations using standard OPLS-AA classical potentials can reasonably predict binding constants for monovalent cations, but these procedures fail for divalent cations. We also consider computationally efficient static procedure using Kohn–Sham density functional theory (DFT) and cluster-continuum modeling that accounts for local microsolvation and pH effects. This approach accurately predicts binding energies for monovalent and divalent cations with an average error of 3.2 kcal mol −1 compared to experiment. This static procedure thus should be useful for future molecular screening efforts, and high absolute errors in the literature may be due to inadequate modeling of local solvent and pH effects. 
    more » « less
  3. Abstract

    Reaction of {LiC6H2−2,4,6‐Cyp3⋅Et2O}2(Cyp=cyclopentyl) (1) of the new dispersion energy donor (DED) ligand, 2,4,6‐triscyclopentylphenyl with SnCl2afforded a mixture of the distannene {Sn(C6H2−2,4,6‐Cyp3)2}2(2), and the cyclotristannane {Sn(C6H2−2,4,6‐Cyp3)2}3(3).2is favored in solution at higher temperature (345 K or above) whereas3is preferred near 298 K. Van't Hoff analysis revealed the3to2conversion has a ΔH=33.36 kcal mol−1and ΔS=0.102 kcal mol−1 K−1, which gives a ΔG300 K=+2.86 kcal mol−1, showing that the conversion of3to2is an endergonic process. Computational studies show that DED stabilization in3is −28.5 kcal mol−1per {Sn(C6H2−2,4,6‐Cyp3)2unit, which exceeds the DED energy in2of −16.3 kcal mol−1per unit. The data clearly show that dispersion interactions are the main arbiter of the3to2equilibrium. Both2and3possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results).

     
    more » « less
  4. Abstract

    Reaction of {LiC6H2−2,4,6‐Cyp3⋅Et2O}2(Cyp=cyclopentyl) (1) of the new dispersion energy donor (DED) ligand, 2,4,6‐triscyclopentylphenyl with SnCl2afforded a mixture of the distannene {Sn(C6H2−2,4,6‐Cyp3)2}2(2), and the cyclotristannane {Sn(C6H2−2,4,6‐Cyp3)2}3(3).2is favored in solution at higher temperature (345 K or above) whereas3is preferred near 298 K. Van't Hoff analysis revealed the3to2conversion has a ΔH=33.36 kcal mol−1and ΔS=0.102 kcal mol−1 K−1, which gives a ΔG300 K=+2.86 kcal mol−1, showing that the conversion of3to2is an endergonic process. Computational studies show that DED stabilization in3is −28.5 kcal mol−1per {Sn(C6H2−2,4,6‐Cyp3)2unit, which exceeds the DED energy in2of −16.3 kcal mol−1per unit. The data clearly show that dispersion interactions are the main arbiter of the3to2equilibrium. Both2and3possess large dispersion stabilization energies which suppress monomer dissociation (supported by EDA results).

     
    more » « less
  5. 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