计算出的RESP2电荷有氢原子的电荷是负的，合理吗？

lufenghxu

H    15.848140   17.652703   56.269738   0.1276885005
H    14.867852   10.362404   55.852962   0.1448700466
H    13.745314   11.432494   55.014873   0.1448700466
H    15.188147   13.382005   55.474598   0.0850071003
H    15.400735   15.227866   56.568994   0.2284932668
H    14.044425   11.644233   59.137514   0.1459053741
H    15.044953   10.475765   58.261874   0.1459053741
H    11.037988   11.084650   55.354882   0.0460819627
H    11.226706   12.602799   56.222448   0.0460819627
H    21.135241   22.181999   56.868406   0.0746726079
H    20.548273   20.980308   58.020236   0.0746726079
H    21.786227   20.548532   56.828267   0.0746726079
H    18.886499   21.333500   56.173500  -0.0373506851
H    20.118517   20.968049   54.987160  -0.0373506851
上述是计算出的resp2电荷的chg文件中的一部分，但是有几个氢原子的resp2电荷是负值（倒数第一，第二行），出现这种情况合理吗？另外，溶剂环境下是选择smd还是cpcm？

sobereva

没必然问题，尤其是不给出体系结构的情况下

看下文了解溶剂模型之间的差别。老老实实用scrf默认的IEFPCM。Gaussian用户没事绝对甭用CPCM！！！下文明确说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）