Universal Solvation Model Based on the Generalized Born Approximation with Asymmetric Descreening

被引:118
作者
Marenich, Aleksandr V.
Cramer, Christopher J. [1 ]
Truhlar, Donald G.
机构
[1] Univ Minnesota, Dept Chem, Minneapolis, MN 55455 USA
基金
美国国家科学基金会;
关键词
QUANTUM-MECHANICAL CALCULATIONS; REACTION FIELD MODEL; FREE-ENERGIES; DENSITY FUNCTIONALS; CONTINUUM MODEL; BASIS-SET; THERMOCHEMICAL KINETICS; MOLECULAR-MECHANICS; ATOMIC CHARGES; ACCURATE;
D O I
10.1021/ct900312z
中图分类号
O64 [物理化学(理论化学)、化学物理学];
学科分类号
070304 ; 081704 ;
摘要
We present a new self-consistent reaction field continuum solvation model based on the generalized Born (GB) approximation for the bulk electrostatic contribution to the free energy of solvation. The new model improves on the earlier SM8 model by using the asymmetric descreening algorithm of Grycuk to treat dielectric descreening effects rather than the Coulomb field approximation; it will be called Solvation Model 8 with asymmetric descreening (SM8AD). The SM8AD model is applicable to any charged or uncharged solute in any solvent or liquid medium for which a few key descriptors are known, in particular dielectric constant, refractive index, bulk surface tension, and acidity and basicity parameters. It does not require the user to assign molecular mechanics types to an atom or a group; all parameters are unique and continuous functions of geometry. This model employs a single set of parameters (solvent acidity-dependent intrinsic Coulomb radii for the treatment of bulk electrostatics and solvent description-dependent atomic surface tensions coefficients for the treatment of nonelectrostatic and shortrange electrostatic effects). The SM8AD model was optimized over 26 combinations of theoretical levels including various basis sets (MIDI!, 6-31 G(star), 6-31 + G(star), 6-31 + G(star star), 6-31 G(star star), cc-pVDZ, DZVP, 6-31B(star)) and electronic structure methods (M05-2X, M05, M06-2X, M06, M06-HF, M06-L, mPW1 PW, mPWPW, B3LYP, HF). It may be used with confidence with any level of electronic structure theory as long as self-consistently polarized Charge Model 4 or other self-consistently polarized charges compatible with CM4 charges are used, for example, CM4M charges can be used. With M05-2X/6-31G(star), the SM8AD model achieves a mean unsigned error of 0.6 kcal/mol on average over 2 560 solvation free energies of tested aqueous and nonaqueous neutral solutes and a mean unsigned error of 3.9 kcal/mol on average over 332 solvation free energies of aqueous and nonaqueous ions.
引用
收藏
页码:2447 / 2464
页数:18
相关论文
共 113 条
[1]  
Abraham M.H., 1994, QUANTITATIVE TREATME, P83
[3]   HYDROGEN-BONDING .31. CONSTRUCTION OF A SCALE OF SOLUTE EFFECTIVE OR SUMMATION HYDROGEN-BOND BASICITY [J].
ABRAHAM, MH .
JOURNAL OF PHYSICAL ORGANIC CHEMISTRY, 1993, 6 (12) :660-684
[4]   HYDROGEN-BONDING .7. A SCALE OF SOLUTE HYDROGEN-BOND ACIDITY BASED ON LOG K-VALUES FOR COMPLEXATION IN TETRACHLOROMETHANE [J].
ABRAHAM, MH ;
GRELLIER, PL ;
PRIOR, DV ;
DUCE, PP ;
MORRIS, JJ ;
TAYLOR, PJ .
JOURNAL OF THE CHEMICAL SOCIETY-PERKIN TRANSACTIONS 2, 1989, (06) :699-711
[5]   Exchange functionals with improved long-range behavior and adiabatic connection methods without adjustable parameters:: The mPW and mPW1PW models [J].
Adamo, C ;
Barone, V .
JOURNAL OF CHEMICAL PHYSICS, 1998, 108 (02) :664-675
[6]  
[Anonymous], 2017, J MOL STRUCT, DOI DOI 10.1016/J.MOLSTRUC.2017.03.014
[7]  
Babu CS, 1999, CHEM PHYS LETT, V310, P225, DOI 10.1016/S0009-2614(99)00714-9
[8]   First principles implementation of solvent effects without outlying charge error [J].
Baldridge, K ;
Klamt, A .
JOURNAL OF CHEMICAL PHYSICS, 1997, 106 (16) :6622-6633
[9]   Interpreting the Coulomb-field approximation for generalized-Born electrostatics using boundary-integral equation theory [J].
Bardhan, Jaydeep P. .
JOURNAL OF CHEMICAL PHYSICS, 2008, 129 (14)
[10]   Computation of protein pK's values by an integrated density functional theory/Polarizable Continuum Model approach [J].
Barone, V ;
Improta, R ;
Rega, N .
THEORETICAL CHEMISTRY ACCOUNTS, 2004, 111 (2-6) :237-245