关于离子液体对汞氧化过程电子转移计算求助

天东耳日

老师好，我目前在计算[Bmim]Cl对汞的氧化过程，我的实验结果是[Bmim]Cl在150℃下将Hg0氧化为Hg2+，我想用gaussian计算研究一下这个过程汞的电子转移。
我用PBE0-D3(BJ)方法，基组对阳离子用6-311g**，阴离子用6-311+g**，汞用SDD，体系总电荷为0，自旋多重度为1。
计算后发现Hg的MK电荷值为0.06，没有发生电子转移。
因为我之前只做过简单的吸附计算（如载硫活性炭吸附汞），如果有化学吸附，汞的电荷值肯定远大于0。
想请教一下老师，是我的输入文件的设置不适用于研究这一氧化过程？还是其他什么原因导致我算的不对？
正确研究这一氧化反应的计算过程应该是什么样的？

sobereva

原子电荷算出来的大小和你用的原子电荷计算方法直接相关，看《原子电荷计算方法的对比》（http://www.whxb.pku.edu.cn/CN/abstract/abstract27818.shtml）。一方面，别光用MK电荷说事，另一方面，别以为实际中Hg的电子转移量会很大，哪怕MK电荷算出来只有0.06也是有一定转移。不要以为原子电荷会像氧化态的数值那么夸张，看《使用Multiwfn通过LOBA方法计算氧化态》（http://sobereva.com/362）。

gjf里的关键词本身并没什么问题，但是没有考虑溶剂效应。