Jorgensen William L, Jensen Kasper P, Alexandrova Anastassia N
Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520-8107.
J Chem Theory Comput. 2007 Nov 1;3(6):1987-1992. doi: 10.1021/ct7001754.
Variations in hydrogen-bond strengths are investigated for complexes of nine para-substituted phenols (XPhOH) with a water molecule and chloride ion. Results from ab initio HF/6-311+G(d, p) and MP2/6-311+G(d, p)//HF/6-311+G(d, p) calculations are compared with those from the OPLS-AA and OPLS/CM1A force fields. In the OPLS-AA model, the partial charges on the hydroxyl group of phenol are not affected by the choice of para substituent, while the use of CM1A charges in the OPLS/CM1A approach does provide charge redistribution. The ab initio calculations reveal a 2.0-kcal/mol range in hydrogen-bond strengths for the XPhOH⋯OH(2) complexes in the order X = NO(2) > CN > CF(3) > Cl > F > H >OH >CH(3) > NH(2). The pattern is not well-reproduced with OPLS-AA, which also compresses the variation to 0.7 kcal/mol. However, the OPLS/CM1A results are in good accord with the ab initio findings for both the ordering and range, 2.3 kcal/mol. The hydrogen bonding is, of course, weaker with XPhOH as acceptor, the order for X is largely inverted, and the range is reduced to ca. 1.0 kcal/mol. The substituent effects are found to be much greater for the chloride ion complexes with a range of 11 kcal/mol. For quantitative treatment of such strong ion-molecule interactions the need for fully polarizable force fields is demonstrated.
研究了9种对取代苯酚(XPhOH)与水分子和氯离子形成的配合物中氢键强度的变化。将从头算HF/6 - 311 + G(d, p)和MP2/6 - 311 + G(d, p)//HF/6 - 311 + G(d, p)计算结果与OPLS - AA和OPLS/CM1A力场的结果进行了比较。在OPLS - AA模型中,苯酚羟基上的部分电荷不受对位取代基选择的影响,而在OPLS/CM1A方法中使用CM1A电荷确实会导致电荷重新分布。从头算计算表明,XPhOH⋯OH(2)配合物的氢键强度范围为2.0千卡/摩尔,顺序为X = NO(2) > CN > CF(3) > Cl > F > H > OH > CH(3) > NH(2)。OPLS - AA不能很好地重现这种模式,它还将变化压缩到0.7千卡/摩尔。然而,OPLS/CM1A的结果在顺序和范围上与从头算结果都很吻合,为2.3千卡/摩尔。当然,当XPhOH作为受体时氢键较弱,X的顺序基本颠倒,范围减小到约1.0千卡/摩尔。对于氯离子配合物,取代基效应要大得多,范围为11千卡/摩尔。对于这种强离子 - 分子相互作用的定量处理,证明了需要完全可极化的力场。