关于高斯计算报错L602

NoleTovey

大家好，我是一个新手，最近想用Gaussian09来优化Pd团簇，第一次计算就出现了L602报错，之后我看了高斯常见报错的帖子，按照里面的建议增加了pop=readradii,但仍出现了L602错误，不太明白为什么，希望可以得到指点。
下面是我的输入文件：
%nprocshared=1
%mem=500MW
%chk=Pd2.chk
# opt b3lyp/lanl2dz guess=save
pop=(mk,hirshfeld,nbo,savenbos,savemulliken,ReadRadii) geom=connectivity polar
Pd2
0 1
Pd                -3.12582054    1.08462005    0.00000000
Pd                -0.87417951    1.04441218    0.00000000
1 2 1.0
2
Pd 1.6


下面是输出文件：
Merz-Kollman atomic radii used.
Read replacement radii for each type of atom:
Atom Element Radius
    1    46    0.00
GetVDW:  no radius for atom   1 atomic number  46.
Error termination via Lnk1e in /home/stud06/g09/l602.exe at Fri Feb 14 14:20:09 2020.
Job cpu time:       0 days  0 hours  0 minutes 15.4 seconds.
File lengths (MBytes):  RWF=     72 Int=      0 D2E=      0 Chk=      2 Scr=      2

NoleTovey

输入文件复制的有些问题，上个图片

snljty

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

sobereva

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

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

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

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

ABetaCarw

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

lijiayisjtu

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

这关键词很多都是多余

NoleTovey

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

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

NoleTovey

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

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

多谢老师指教