Learning from the density to correct total energy and forces in first principle simulations.
TL;DR: A new molecular simulation framework is proposed that combines the transferability, robustness, and chemical flexibility of an ab initio method with the accuracy and efficiency of a machine learning model.
read more
Abstract: We propose a new molecular simulation framework that combines the transferability, robustness, and chemical flexibility of an ab initio method with the accuracy and efficiency of a machine learning model. The key to achieve this mix is to use a standard density functional theory (DFT) simulation as a preprocessor for the atomic and molecular information, obtaining a good quality electronic density. General, symmetry preserving, atom-centered electronic descriptors are then built from this density to train a neural network to correct the baseline DFT energies and forces. These electronic descriptors encode much more information than local atomic environments, allowing a simple neural network to reach the accuracy required for the problem of study at a negligible additional cost. The balance between accuracy and efficiency is determined by the baseline simulation. This is shown in results where high level quantum chemical accuracy is obtained for simulations of liquid water at standard DFT cost or where high level DFT-accuracy is achieved in simulations with a low-level baseline DFT calculation at a significantly reduced cost.
read more
Chat with Paper
AI Agents for this Paper
Find similar papers on Google Scholar, PubMed and Arxiv
Write a critical review of this paper
Analyze citations of this paper to find unaddressed research gaps
Figures
![FIG. 10: a) and b) Radial distribution functions and c) vibrational density of states. MLCF is used to accelerate DFT calculations by correcting for both functional and basis-set related convergence errors, combining MLCF with the time step mixing method proposed by Anglada et al. [58]. Interpolation between a quick and dirty PBE(SZ) simulation and a slow and more accurate vdW-cx(DZP) simulations. n indicates the number of fast steps needed before correcting with a well converged step.](/figures/figure10-1-3y1oss5pch7k.png)
FIG. 10: a) and b) Radial distribution functions and c) vibrational density of states. MLCF is used to accelerate DFT calculations by correcting for both functional and basis-set related convergence errors, combining MLCF with the time step mixing method proposed by Anglada et al. [58]. Interpolation between a quick and dirty PBE(SZ) simulation and a slow and more accurate vdW-cx(DZP) simulations. n indicates the number of fast steps needed before correcting with a well converged step. 
FIG. 2: Energies of water hexamers with respect to the prism isomer, which is correctly determined to be the most stable structure by the MLCF. 
FIG. 11: Speed-up obtained by mixing PBE(SZ) with vdW-cx(DZP) for finite size systems compared to the speed-up that would be achieved for an algorithm that scales strictly cubic in the number of oribtals. The color gradient indicates the method’s reliability to reproduce the reference method for a given value of n. 
FIG. 4: Root mean squared errors (RMSE) of ∆-SchNet (top) and ∆-WACSF (bottom) models trained for different cutoff radii (see legend) on water monomers, dimers and trimers and tested on random clusters of increasing size. The errors are normalized by dividing them through the RMSE of the underlying baseline method (vdW-cx) and compared to those of MLCF (purple line) which can also be found in Tab. II. 
TABLE III: Optimal hyperparameters determined through cross-validation for SchNet. The hyperparameters that were optimized were the the dimension of the embedding space denoted as ”features” and the number of interaction layers. An ”x” in the column named ∆ denotes that the model was trained to correct vdW-cx. 
TABLE IV: Optimal hyperparameters determined through cross-validation for WACSF. The hyperparameters that were optimized were the the number of angular and radial basis functions and the number of nodes in the three layer fully connected neural network. An ”x” in the column named ∆ denotes that the model was trained to correct vdW-cx.
Citations
Mean-Field Density Matrix Decompositions
TL;DR: In this paper, the authors introduce new and robust decompositions of mean-field Hartree-Fock (HF) and Kohn-Sham density functional theory (KS-DFT) relying on the use of localized molecular orbitals and physically sound charge population protocols.
Generalizing deep learning electronic structure calculation to the plane-wave basis
Xiaoxun Gong,Steven G. Louie,Wenhui Duan,Yong Xu +3 more
TL;DR: Researchers propose a real-space reconstruction method to convert plane-wave DFT results to atomic-orbital basis, enabling deep learning electronic structure calculations in the plane-wave basis with high accuracy and efficiency, bridging a longstanding gap between methods.
References
Generalized Gradient Approximation Made Simple
TL;DR: A simple derivation of a simple GGA is presented, in which all parameters (other than those in LSD) are fundamental constants, and only general features of the detailed construction underlying the Perdew-Wang 1991 (PW91) GGA are invoked.
Density‐functional thermochemistry. III. The role of exact exchange
TL;DR: In this article, a semi-empirical exchange correlation functional with local spin density, gradient, and exact exchange terms was proposed. But this functional performed significantly better than previous functionals with gradient corrections only, and fits experimental atomization energies with an impressively small average absolute deviation of 2.4 kcal/mol.
The SIESTA method for ab initio order-N materials simulation
José M. Soler,Emilio Artacho,Julian D. Gale,Alberto García,Javier Junquera,Javier Junquera,Pablo Ordejón,Daniel Sánchez-Portal +7 more
TL;DR: In this paper, a selfconsistent density functional method using standard norm-conserving pseudopotentials and a flexible, numerical linear combination of atomic orbitals basis set, which includes multiple-zeta and polarization orbitals, was developed and implemented.
Van der Waals Density Functional for General Geometries
TL;DR: The proposed functional, which includes van der Waals forces in a seamless fashion, is applied to rare gas and benzene dimers, where it is shown to give a realistic description.
5.6K
Nosé-Hoover chains : the canonical ensemble via continuous dynamics
TL;DR: In this paper, a modification of the Nose-Hoover dynamics is proposed which includes not a single thermostat variable but a chain of variables, Nose chains, which gives the canonical distribution where the simple formalism fails.