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 条
[1]   Polarizable empirical force field for the primary and secondary alcohol series based on the classical drude model [J].
Anisimov, Victor M. ;
Vorobyov, Igor V. ;
Roux, Benoit ;
MacKerell, Alexander D., Jr. .
JOURNAL OF CHEMICAL THEORY AND COMPUTATION, 2007, 3 (06) :1927-1946
[2]  
[Anonymous], 2016, GAUSSIAN 16 REV B01
[3]   ATOM DIPOLE INTERACTION MODEL FOR MOLECULAR POLARIZABILITY - APPLICATION TO POLYATOMIC-MOLECULES AND DETERMINATION OF ATOM POLARIZABILITIES [J].
APPLEQUIST, J ;
CARL, JR ;
FUNG, KK .
JOURNAL OF THE AMERICAN CHEMICAL SOCIETY, 1972, 94 (09) :2952-+
[4]   Parametrizing a polarizable force field from ab initio data.: I.: The fluctuating point charge model [J].
Banks, JL ;
Kaminski, GA ;
Zhou, RH ;
Mainz, DT ;
Berne, BJ ;
Friesner, RA .
JOURNAL OF CHEMICAL PHYSICS, 1999, 110 (02) :741-754
[5]   A combined electronegativity equalization and electrostatic potential fit method for the determination of atomic point charges [J].
Berente, Imre ;
Czinki, Eszter ;
Naray-Szabo, Gabor .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 2007, 28 (12) :1936-1942
[6]  
BOYD S, 2006, CONVEX OPTIMIZATIO 3, pCH11
[7]   CHARMM - A PROGRAM FOR MACROMOLECULAR ENERGY, MINIMIZATION, AND DYNAMICS CALCULATIONS [J].
BROOKS, BR ;
BRUCCOLERI, RE ;
OLAFSON, BD ;
STATES, DJ ;
SWAMINATHAN, S ;
KARPLUS, M .
JOURNAL OF COMPUTATIONAL CHEMISTRY, 1983, 4 (02) :187-217
[8]  
Bultinck P, 2004, J PHYS CHEM A, V108, P10359, DOI [10.1021/jp046928l, 10.1021/jp0469281]
[9]   The electronegativity equalization method II: Applicability of different atomic charge schemes [J].
Bultinck, P ;
Langenaeker, W ;
Lahorte, P ;
De Proft, F ;
Geerlings, P ;
Van Alsenoy, C ;
Tollenaere, JP .
JOURNAL OF PHYSICAL CHEMISTRY A, 2002, 106 (34) :7895-7901
[10]   The electronegativity equalization method I: Parametrization and validation for atomic charge calculations [J].
Bultinck, P ;
Langenaeker, W ;
Lahorte, P ;
De Proft, F ;
Geerlings, P ;
Waroquier, M ;
Tollenaere, JP .
JOURNAL OF PHYSICAL CHEMISTRY A, 2002, 106 (34) :7887-7894