An Extended Charge Equilibration Method.

We present a method for estimating partial atomic charges that uses all of the measured ionization energies (first, second, third, etc.) for every atom in the periodic table. We build on the charge equilibration (Qeq) method of Rappé and Goddard (which used only the first ionization energies) but reduce the number of ad hoc parameters from at least one for every type of atom to just two global parameters: a dielectric strength and a modified parameter for hydrogen atoms. Periodic electrostatic interactions are calculated via Ewald sums, and the partial charges are determined by simultaneously solving a system of linear equations; no iteration is required. We compare the predicted partial atomic charges of this extended charge equilibration (EQeq) scheme against plane-wave density-functional theory derived charges determined via the REPEAT method for 12 diverse metal-organic frameworks (MOFs). We also compare EQeq charges against ChelpG charges calculated using nonperiodic MOF fragments, as well as against Qeq charges as implemented in Accelrys Materials Studio. We demonstrate that for the purpose of ranking MOFs from best to worst for carbon capture applications, EQeq charges perform as well as charges derived from electrostatic potentials, but EQeq requires only a tiny fraction of the computational cost (seconds vs days for the MOFs studied). The source code for the EQeq algorithm is provided.

[1]  Rachel B. Getman,et al.  Review and analysis of molecular simulations of methane, hydrogen, and acetylene storage in metal-organic frameworks. , 2012, Chemical reviews.

[2]  U. Suter,et al.  Atomic Charges for Classical Simulations of Polar Systems , 2004 .

[3]  C. Ania,et al.  H2, N2, CO, and CO2 sorption properties of a series of robust sodalite-type microporous coordination polymers. , 2006, Inorganic chemistry.

[4]  Michael O’Keeffe,et al.  Exceptional chemical and thermal stability of zeolitic imidazolate frameworks , 2006, Proceedings of the National Academy of Sciences.

[5]  Darrin M. York,et al.  A chemical potential equalization method for molecular simulations , 1996 .

[6]  Johann Gasteiger,et al.  Electronegativity equalization: application and parametrization , 1985 .

[7]  R. Bader.  Atoms in molecules , 1990 .

[8]  David S Sholl,et al.  Chemically Meaningful Atomic Charges That Reproduce the Electrostatic Potential in Periodic and Nonperiodic Materials. , 2010, Journal of chemical theory and computation.

[9]  Toon Verstraelen,et al.  The electronegativity equalization method and the split charge equilibration applied to organic systems: parametrization, validation, and comparison. , 2009, The Journal of chemical physics.

[10]  Christopher E. Wilmer,et al.  Nanoscale forces and their uses in self-assembly. , 2009, Small.

[11]  David S. Sholl,et al.  Progress, Opportunities, and Challenges for Applying Atomically Detailed Modeling to Molecular Adsorption and Transport in Metal−Organic Framework Materials , 2009 .

[12]  C. Wilmer,et al.  Towards rapid computational screening of metal-organic frameworks for carbon dioxide capture: Calculation of framework charges via charge equilibration , 2011 .

[13]  Chongli Zhong,et al.  Estimation of Framework Charges in Covalent Organic Frameworks Using Connectivity-Based Atom Contribution Method , 2010 .

[14]  Sankar Nair,et al.  Finding MOFs for highly selective CO2/N2 adsorption using materials screening based on efficient assignment of atomic point charges. , 2012, Journal of the American Chemical Society.

[15]  Jackson,et al.  Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. , 1992, Physical review. B, Condensed matter.

[16]  Parameter determination for the charge equilibration method including third- and fourth-order terms applied to non-metallic compounds , 2010 .

[17]  H. Fjellvåg,et al.  An in situ high-temperature single-crystal investigation of a dehydrated metal-organic framework compound and field-induced magnetization of one-dimensional metal-oxygen chains. , 2005, Angewandte Chemie.

[18]  J. Marrot,et al.  A breathing hybrid organic-inorganic solid with very large pores and high magnetic characteristics. , 2002, Angewandte Chemie.

[19]  M. O'keeffe,et al.  Design and synthesis of an exceptionally stable and highly porous metal-organic framework , 1999, Nature.

[20]  Omar M Yaghi,et al.  Effects of functionalization, catenation, and variation of the metal oxide and organic linking units on the low-pressure hydrogen adsorption properties of metal-organic frameworks. , 2006, Journal of the American Chemical Society.

[21]  Kresse,et al.  Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. , 1996, Physical review. B, Condensed matter.

[22]  A. Rappé,et al.  Toward an Understanding of Zeolite Y as a Cracking Catalyst with the Use of Periodic Charge Equilibration , 1996 .

[23]  W. C. Lineberger,et al.  Binding energies in atomic negative ions , 1975 .

[24]  T. Darden,et al.  Molecular dynamics simulations of biomolecules: long-range electrostatic effects. , 1999, Annual review of biophysics and biomolecular structure.

[25]  Kevin L. Shaw,et al.  Charge–charge interactions influence the denatured state ensemble and contribute to protein stability , 2000, Protein science : a publication of the Protein Society.

[26]  Joachim Sauer,et al.  Molecular models in ab initio studies of solids and surfaces: from ionic crystals and semiconductors to catalysts , 1989 .

[27]  A. Matzger,et al.  Dramatic tuning of carbon dioxide uptake via metal substitution in a coordination polymer with cylindrical pores. , 2008, Journal of the American Chemical Society.

[28]  Lo,et al.  A chemically functionalizable nanoporous material , 1999, Science.

[29]  Toon Verstraelen,et al.  Assessment of Atomic Charge Models for Gas-Phase Computations on Polypeptides. , 2012, Journal of chemical theory and computation.

[30]  Donald E. Williams,et al.  Representation of the molecular electrostatic potential by a net atomic charge model , 1981 .

[31]  R. Snurr,et al.  Using molecular simulation to characterise metal-organic frameworks for adsorption applications. , 2009, Chemical Society reviews.

[32]  T. Verstraelen,et al.  The Significance of Parameters in Charge Equilibration Models. , 2011, Journal of chemical theory and computation.

[33]  Michael O'Keeffe,et al.  Systematic Design of Pore Size and Functionality in Isoreticular MOFs and Their Application in Methane Storage , 2002, Science.

[34]  Randall Q Snurr,et al.  Screening of metal-organic frameworks for carbon dioxide capture from flue gas using a combined experimental and modeling approach. , 2009, Journal of the American Chemical Society.

[35]  A. Matzger,et al.  Gas and liquid phase adsorption in isostructural Cu3[biaryltricarboxylate]2 microporous coordination polymers. , 2011, Chemical communications.

[36]  R. Levy,et al.  Computer simulations with explicit solvent: recent progress in the thermodynamic decomposition of free energies and in modeling electrostatic effects. , 1998, Annual review of physical chemistry.

[37]  P. Kollman,et al.  An approach to computing electrostatic charges for molecules , 1984 .

[38]  B. Grzybowski,et al.  Precision assembly of oppositely and like-charged nanoobjects mediated by charge-induced dipole interactions. , 2010, Nano letters.

[39]  J. Gasteiger,et al.  ITERATIVE PARTIAL EQUALIZATION OF ORBITAL ELECTRONEGATIVITY – A RAPID ACCESS TO ATOMIC CHARGES , 1980 .

[40]  Wilfried J. Mortier,et al.  Electronegativity-equalization method for the calculation of atomic charges in molecules , 1986 .

[41]  W. Goddard,et al.  Charge equilibration for molecular dynamics simulations , 1991 .

[42]  A. Matzger,et al.  Porous crystal derived from a tricarboxylate linker with two distinct binding motifs. , 2007, Journal of the American Chemical Society.

[43]  T. Darden,et al.  Toward the Accurate Modeling of DNA: The Importance of Long-Range Electrostatics , 1995 .

[44]  M. Hirscher,et al.  Hydrogen adsorption in a nickel based coordination polymer with open metal sites in the cylindrical cavities of the desolvated framework. , 2006, Chemical communications.

[45]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[46]  Tom K Woo,et al.  Electrostatic Potential Derived Atomic Charges for Periodic Systems Using a Modified Error Functional. , 2009, Journal of chemical theory and computation.

[47]  Susan B. Sinnott,et al.  Charge optimized many-body potential for the Si/SiO2 system , 2007 .

[48]  R. A. Nistor,et al.  A generalization of the charge equilibration method for nonmetallic materials. , 2006, The Journal of chemical physics.

[49]  C. E. Moore.  Ionization Potentials and Ionization Limits Derived from the Analyses of Optical Spectra , 1970 .