Evaluating charge equilibration methods to generate electrostatic fields in nanoporous materials

Authors: Daniele Ongari1*, Peter G. Boyd1, Amber K. Mace1, Berend Smit1*, Ozge Kadioglu2, Seda Keskin2

  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
  • Corresponding authors emails: daniele.ongari@epfl.ch, berend.smit@epfl.ch

DOI10.24435/materialscloud:2018.0017/v1 (version v1, submitted on 08 November 2018)

[There are newer versions. Click here to view the latest version v2]

How to cite this entry

Daniele Ongari, Peter G. Boyd, Amber K. Mace, Berend Smit, Ozge Kadioglu, Seda Keskin, Evaluating charge equilibration methods to generate electrostatic fields in nanoporous materials, Materials Cloud Archive (2018), doi: 10.24435/materialscloud:2018.0017/v1.


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.

Materials Cloud sections using this data

No Explore or Discover sections associated with this archive entry.


File name Size Description
MD5MD5: 8f035ccb155ed7073a91f6dfe026ad19
29.3 MiB Inputs for Qeq methods, adsorption results and charge comparison by atom type.


Files and data are licensed under the terms of the following license: Creative Commons Attribution 4.0 International.

External references

Journal reference
Manuscript under revision


partial charges Qeq charge equilibration methods MARVEL

Version history

29 November 2018

08 November 2018 [This version]