Analyzing the errors of DFT approximations for compressed water systems

Journal Article · Mon Jul 07 04:00:00 EDT 2014 · Journal of Chemical Physics

DOI:https://doi.org/10.1063/1.4885440· OSTI ID:22308740

Alfè, D. ^[1]; Bartók, A. P.; Csányi, G. ^[2]

---

1. Department of Earth Sciences, UCL, London WC1E 6BT (United Kingdom)
2. Department of Engineering, University of Cambridge, Cambridge CB2 1PZ (United Kingdom)

We report an extensive study of the errors of density functional theory (DFT) approximations for compressed water systems. The approximations studied are based on the widely used PBE and BLYP exchange-correlation functionals, and we characterize their errors before and after correction for 1- and 2-body errors, the corrections being performed using the methods of Gaussian approximation potentials. The errors of the uncorrected and corrected approximations are investigated for two related types of water system: first, the compressed liquid at temperature 420 K and density 1.245 g/cm{sup 3} where the experimental pressure is 15 kilobars; second, thermal samples of compressed water clusters from the trimer to the 27-mer. For the liquid, we report four first-principles molecular dynamics simulations, two generated with the uncorrected PBE and BLYP approximations and a further two with their 1- and 2-body corrected counterparts. The errors of the simulations are characterized by comparing with experimental data for the pressure, with neutron-diffraction data for the three radial distribution functions, and with quantum Monte Carlo (QMC) benchmarks for the energies of sets of configurations of the liquid in periodic boundary conditions. The DFT errors of the configuration samples of compressed water clusters are computed using QMC benchmarks. We find that the 2-body and beyond-2-body errors in the liquid are closely related to similar errors exhibited by the clusters. For both the liquid and the clusters, beyond-2-body errors of DFT make a substantial contribution to the overall errors, so that correction for 1- and 2-body errors does not suffice to give a satisfactory description. For BLYP, a recent representation of 3-body energies due to Medders, Babin, and Paesani [J. Chem. Theory Comput. 9, 1103 (2013)] gives a reasonably good way of correcting for beyond-2-body errors, after which the remaining errors are typically 0.5 mE{sub h} ≃ 15 meV/monomer for the liquid and the clusters.

• OSTI ID: 22308740
• Journal Information: Journal of Chemical Physics, Journal Name: Journal of Chemical Physics Journal Issue: 1 Vol. 141; ISSN JCPSA6; ISSN 0021-9606
• Country of Publication: United States
• Language: English

Related Subjects

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY
APPROXIMATIONS
CORRECTIONS
DENSITY
DENSITY FUNCTIONAL METHOD
ERRORS
LIQUIDS
MOLECULAR DYNAMICS METHOD
MONTE CARLO METHOD
NEUTRON DIFFRACTION
SIMULATION
WATER