同一结构几何优化与频率计算所得电荷不同

Happypant

采用B3LYP/6-31+g(d,p)方法，用PCM模拟水环境，研究镁离子对腺嘌呤-胸腺嘧啶碱基对的质子转移的影响。反应物、过渡态、产物的几何优化与频率计算所得输出文件中Mulliken电荷有着不小差距。有的是用不同电脑算的，但是应该不会差这么多。我想知道是不是计算出了问题。sob老师请指教。

sobereva

带弥散函数时候算的Mulliken没意义，这是常识。看《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）
两个任务给出的Mulliken电荷都是相同的，只不过是你读错了位置而已。

写=pcm完全是多余的，这本来就是默认的

优化这种体系用2-zeta加弥散是极差的选择，用6-311G**更便宜还更好，仔细看
谈谈弥散函数和“月份”基组
http://sobereva.com/119

Happypant

sobereva 发表于 2019-5-11 11:31
带弥散函数时候算的Mulliken没意义，这是常识。看《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn ...

多谢sob老师。