利用Gaussian计算含C、N、O、H的阴离子静电势时，如何合理设置电荷等参数

林森

化合物结构是含有C、N、O、H的-1价的阴离子，目标想计算该结构的静电势；方法为DFT-B3LYP，基组设置为6-311G+(d，p)；

当charge：-1，spin：singlet时，计算结束用Multiwfn得到各极值点的静电势均为相似的负数（范围约：-1~-10kcal/mol），VMD呈现的静电势图为全蓝色；

当charge：0，spin=doublet时，计算结束可以得到有正有负的极值点静电势（正值最大可达58.21、负值可达-38.02kcal/mol）；

因为希望可以通过计算该结构与另一结构的静电势来判断可能发生相互作用力的位点，所以极值点有正有负的静电势结果更符合预期的；

想请教各位老师，①设置charge：0，spin=doublet所得到的静电势结果是否合适？②对于charge：-1，spin：singlet时所得全为负值的结果，有无方法可以改进？谢谢各位老师指点解惑

sobereva

带净正、负电荷的，分子范德华表面静电势分别全为正、负是再正常不过的情况，别胡思乱想以为结果错误
体系实际是什么净电荷，就以什么净电荷做量子化学计算并用Multiwfn正常分析静电势
带净电荷的体系，只需要关注不同位置静电势的相对大小

林森

sobereva 发表于 2025-2-16 01:02
带净正、负电荷的，分子范德华表面静电势分别全为正、负是再正常不过的情况，别胡思乱想以为结果错误
体系 ...

谢谢社长，受益匪浅