Published November 29, 2018 | Version v2
Dataset Open

Evaluating charge equilibration methods to generate electrostatic fields in nanoporous materials

  • 1. Laboratory of Molecular Simulation, Institut des Sciences et Ingénierie Chimiques, Ecole Polytechnique Fédérale de Lausanne (EPFL), Rue de l'Industrie 17, CH-1951 Sion, Valais, Switzerland
  • 2. Department of Chemical and Biological Engineering, Koc University, Rumelifeneri Yolu, Sariyer, 34450 Istanbul, Turkey

* Contact person

Description

Charge equilibration (Qeq) methods can estimate the electrostatic potential of molecules and periodic frameworks by assigning point charges to each atom, using only a small fraction of the resources needed to compute density functional (DFT)-derived charges. This makes possible, for example, the computational screening of thousands of microporous structures to assess their performance for the adsorption of polar molecules. Recently, different variants of the original Qeq scheme were proposed to improve the quality of the computed point charges. One focus of this research was to improve the gas adsorption predictions in Metal Organic Frameworks (MOFs), for which many different structures are available. In this work, we review the evolution of the method from the original Qeq scheme, understanding the role of the different modifications on the final output. We evaluated the result of combining different protocols and set of parameters, by comparing the Qeq charges with high quality DFT-derived DDEC charges for 2338 MOF structures. We focused on the systematic errors that are attributable to specific atom types to quantify the final precision that one can expect from Qeq methods in the context of gas adsorption where the electrostatic potential plays a significant role, namely CO2 and H2S adsorption. In conclusion, both the type of algorithm and the input parameters have a large impact on the resulting charges and we draw some guidelines to help the user to choose the proper combination of the two for obtaining a meaningful set of charges. We show that, considering this set of MOFs, the accuracy of the original Qeq scheme is often still comparable with the most recent variants, even if it clearly fails in the presence of certain atom types, such as alkali metals. Version 2 of the record adds the data underlying the figures in CSV format.

Files

File preview

files_description.md

All files

Files (39.7 MiB)

Name Size
md5:a98bca5b9ac978007cf743fc65a132f8
227 Bytes Preview Download
md5:b03c264df9f4958ea25ceda7ca93a66b
39.7 MiB Preview Download

References

Journal reference (Article where the results are published.)
Daniele Ongari, Peter G Boyd, Ozge Kadioglu, Amber K. Mace, Seda Keskin, and Berend Smit, J. Chem. Theory Comput., Just Accepted Manuscript, doi: 10.1021/acs.jctc.8b00669