The electronegativity equalization method and the split charge equilibration applied to organic systems: Parametrization, validation, and comparison

被引:83
作者
Verstraelen, Toon
Van Speybroeck, Veronique [1 ]
Waroquier, Michel
机构
[1] Univ Ghent, Ctr Mol Modeling CMM, B-9000 Ghent, Belgium
关键词
ab initio calculations; electronegativity; molecular force constants; organic compounds; perturbation theory; polarisability; POLARIZABLE FORCE-FIELD; ATOMIC CHARGES; AB-INITIO; MOLECULAR-MECHANICS; DIPOLE-MOMENT; PARAMETERIZATION; PROTEINS; SIMULATIONS; ZEOLITES; VACUUM;
D O I
10.1063/1.3187034
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
An extensive benchmark of the electronegativity equalization method (EEM) and the split charge equilibration (SQE) model on a very diverse set of organic molecules is presented. These models efficiently compute atomic partial charges and are used in the development of polarizable force fields. The predicted partial charges that depend on empirical parameters are calibrated to reproduce results from quantum mechanical calculations. Recently, SQE is presented as an extension of the EEM to obtain the correct size dependence of the molecular polarizability. In this work, 12 parametrization protocols are applied to each model and the optimal parameters are benchmarked systematically. The training data for the empirical parameters comprise of MP2/Aug-CC-pVDZ calculations on 500 organic molecules containing the elements H, C, N, O, F, S, Cl, and Br. These molecules have been selected by an ingenious and autonomous protocol from an initial set of almost 500 000 small organic molecules. It is clear that the SQE model outperforms the EEM in all benchmark assessments. When using Hirshfeld-I charges for the calibration, the SQE model optimally reproduces the molecular electrostatic potential from the ab initio calculations. Applications on chain molecules, i.e., alkanes, alkenes, and alpha alanine helices, confirm that the EEM gives rise to a divergent behavior for the polarizability, while the SQE model shows the correct trends. We conclude that the SQE model is an essential component of a polarizable force field, showing several advantages over the original EEM.
引用
收藏
页数:19
相关论文
共 73 条
[31]   Special issue on polarization [J].
Jorgensen, William L. .
JOURNAL OF CHEMICAL THEORY AND COMPUTATION, 2007, 3 (06) :1877-1877
[32]   Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids [J].
Jorgensen, WL ;
Maxwell, DS ;
TiradoRives, J .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1996, 118 (45) :11225-11236
[33]   AN ALGORITHM FOR DRAWING GENERAL UNDIRECTED GRAPHS [J].
KAMADA, T ;
KAWAI, S .
INFORMATION PROCESSING LETTERS, 1989, 31 (01) :7-15
[34]   COMPUTER AIDED DESIGN OF EXPERIMENTS [J].
KENNARD, RW ;
STONE, LA .
TECHNOMETRICS, 1969, 11 (01) :137-&
[35]  
Lebedev V.I., 1999, Doklady Math, V59, P477
[36]   MFI fingerprint: How pentasil-induced IR bands shift during zeolite nanogrowth [J].
Lesthaeghe, David ;
Vansteenkiste, Peter ;
Verstraelen, Toon ;
Ghysels, An ;
Kirschhock, Christine E. A. ;
Martens, Johan A. ;
Van Speybroeck, Veronique ;
Waroquier, Michel .
JOURNAL OF PHYSICAL CHEMISTRY C, 2008, 112 (25) :9186-9191
[37]   Semi-empirical atomic charges for use in computational chemistry of molecular sieves [J].
Louwen, JN ;
Vogt, ETC .
JOURNAL OF MOLECULAR CATALYSIS A-CHEMICAL, 1998, 134 (1-3) :63-77
[38]   A polarizable force field and continuum solvation methodology for modeling of protein-ligand interactions [J].
Maple, JR ;
Cao, YX ;
Damm, WG ;
Halgren, TA ;
Kaminski, GA ;
Zhang, LY ;
Friesner, RA .
JOURNAL OF CHEMICAL THEORY AND COMPUTATION, 2005, 1 (04) :694-715
[39]   Parameterization of the electronegativity equalization method based on the charge model 1 [J].
Menegon, G ;
Shimizu, K ;
Farah, JPS ;
Dias, LG ;
Chaimovich, H .
PHYSICAL CHEMISTRY CHEMICAL PHYSICS, 2002, 4 (24) :5933-5936
[40]   Note on an approximation treatment for many-electron systems [J].
Moller, C ;
Plesset, MS .
PHYSICAL REVIEW, 1934, 46 (07) :0618-0622