sobereva 发表于 2020-2-14 22:40
对这体系算电荷完全是多余的，双原子分子各个原子电荷精确为0

没事甭写一大堆乱七八糟的关键词，除非你 ...

多谢老师指教

snljty 发表于 2020-2-14 21:14
先说清楚你到底要算什么。那一大堆布局分析除了NBO没必要用Gaussian，也不嫌输出文件太长，用Multiwfn做非 ...

多谢指教，第一次拿高斯计算，有很多地方都不懂，你的回复对我很有帮助，谢谢

# opt b3lyp/lanl2dz guess=save
pop=(mk,hirshfeld,nbo,savenbos,savemulliken,ReadRadii) geom=connectivity polar

这关键词很多都是多余

说句风凉话，你的输入文件可以作为社长下次培训班讲的反面例子了

对这体系算电荷完全是多余的，双原子分子各个原子电荷精确为0

没事甭写一大堆乱七八糟的关键词，除非你真的非常清楚都是什么含义

对于算原子电荷的目的，用Multiwfn比用Gaussian强得多，看手册Multiwfn手册4.7节的计算例子。Multiwfn在计算拟合静电势电荷的时候都会直接给你建议半径，都不用自己补，省事得多还避免初学者弄错。

顺带一提，建议参加北京科音初级量子化学培训班（http://www.keinsci.com/workshop/KEQC_content.html）真正系统、正确地学一遍量化计算，免得胡算瞎算瞎耽误工夫

先说清楚你到底要算什么。那一大堆布局分析除了NBO没必要用Gaussian，也不嫌输出文件太长，用Multiwfn做非常方便，输出结果也好看。
lanl2dz这个基组现在看差的不行，给几何优化用都嫌够呛。就两个原子，何必呢，直接上def2-TZVP。
不知道你为什么要nproc=1，现在电脑再惨也有双核吧...双核并行效率还是挺好的。
polar如果你想算极化率，首先算极化率任务和几何优化任务不要一起用，其次这么小的基组尤其不带弥散函数算极化率是绝绝绝对不能接受的，保底也得def2-SVPD起。你的报错原因不是上面说的这个，但是先把上面的改了。
geom=connectivity和量化方法计算无关，对于非分子力学方法只是gjf和chk文件看图指定了连接方式的作用。后面额外的连接方式和别的额外输入信息的顺序还容易搞混。
两个原子而已也不影响作图，geom=connectivity和后面的几何连接关系都删了，省的麻烦。直接报错原因估计是半径那行数据应该在连接关系之前，不想查手册了，总之geom那个删了没影响还方便。但是别的问题一定要改正后再算不然结果没意义或者其他报错。