求助Gaussian单点能计算命令关键词怎么输才正确

www...

高斯小白 现在开始做计算 以下是计算油酸分子单点能的计算文件内容 出现的问题是一提交任务就显示计算失败，是不是我的关键词输入错误啊。请各位大神解惑，不胜感激。

#p m062x/def2tzvp scrf=(smd,solvent=water) geom=allcheck Guess=read
em=gd3

Title Card Required

0 1
O                  9.66962600    1.53350900   -0.06664200
O                 10.76972100   -0.40736000    0.23412000
H                 11.52308800    0.21288900    0.19755700
C                  9.64388200    0.32760400    0.08198400
C                  8.40607500   -0.53609300    0.12277800
C                  7.10997100    0.25341400   -0.04868800
H                  8.50843300   -1.30225500   -0.65690900
H                  8.41277100   -1.08297100    1.07479800
H                  7.04405000    1.01537500    0.73783300
H                  7.14111100    0.79702900   -1.00107100
C                  5.87077500   -0.64418700   -0.00526300
C                  4.56295100    0.13430300   -0.17402600
H                  5.84851400   -1.19062200    0.94868600
H                  5.94452500   -1.40733200   -0.79348700
C                  3.31861000   -0.75714300   -0.13162600
H                  4.49066300    0.89743600    0.61438600
H                  4.58640600    0.68095000   -1.12790900
H                  3.39016900   -1.51985500   -0.92073000
H                  3.29601500   -1.30471100    0.82186500
C                  2.01110300    0.02166800   -0.29882200
C                  0.76358500   -0.87750500   -0.25173300
H                  2.02968800    0.56612400   -1.25381800
H                  1.93199600    0.78420200    0.48826500
C                 -0.51309200   -0.11272400   -0.46163900
H                  0.86419800   -1.65068000   -1.02895000
H                  0.72730800   -1.40378600    0.71137300
H                 -0.61225000    0.39556000   -1.42384400
C                 -1.49667400    0.00891300    0.43411700
C                 -2.77318900    0.77403600    0.22461900
H                 -1.39740200   -0.49937900    1.39631800
C                 -4.02144500   -0.12402700    0.27396800
H                 -2.87246400    1.54806900    1.00118000
H                 -2.73749700    1.29944000   -0.73899900
H                 -3.94399800   -0.88725100   -0.51266700
H                 -4.03901100   -0.66790000    1.22936200
C                 -5.32872300    0.65554000    0.10824000
C                 -6.57472600   -0.23364900    0.15365500
H                 -5.30700700    1.20241600   -0.84576300
H                 -5.39784500    1.41933300    0.89663500
C                 -7.88315600    0.54482800   -0.01239000
H                 -6.50436500   -0.99777200   -0.63426200
H                 -6.59595500   -0.78027800    1.10784200
H                 -7.95303200    1.30938800    0.77527600
H                 -7.86165100    1.09124700   -0.96675700
C                 -9.12966300   -0.34355100    0.03350100
C                -10.43875800    0.43432000   -0.13318700
H                 -9.15199600   -0.88983400    0.98808100
H                 -9.06021200   -1.10864300   -0.75386800
C                -11.67836200   -0.46246700   -0.08535400
H                -10.50704300    1.19866400    0.65368800
H                -10.41560100    0.97941700   -1.08732900
H                -12.59995700    0.11808300   -0.20697000
H                -11.64960300   -1.21609900   -0.88198900
H                -11.74179200   -0.99513400    0.87154900

1 4 2.0
2 3 1.0 4 1.0
3
4 5 1.0
5 6 1.0 7 1.0 8 1.0
6 9 1.0 10 1.0 11 1.0
7
8
9
10
11 12 1.0 13 1.0 14 1.0
12 15 1.0 16 1.0 17 1.0
13
14
15 18 1.0 19 1.0 20 1.0
16
17
18
19
20 21 1.0 22 1.0 23 1.0
21 24 1.0 25 1.0 26 1.0
22
23
24 27 1.0 28 2.0
25
26
27
28 29 1.0 30 1.0
29 31 1.0 32 1.0 33 1.0
30
31 34 1.0 35 1.0 36 1.0
32
33
34
35
36 37 1.0 38 1.0 39 1.0
37 40 1.0 41 1.0 42 1.0
38
39
40 43 1.0 44 1.0 45 1.0
41
42
43
44
45 46 1.0 47 1.0 48 1.0
46 49 1.0 50 1.0 51 1.0
47
48
49 52 1.0 53 1.0 54 1.0
50
51
52
53
54

wal

#p 6-31g(d) scrf=(smd,solvent=water) geom=connectivity m062x/def2tzvp
这个突兀的6-31g(d)是什么玩意 删了再跑看看
另外问报错最好把输出文件中的报错信息贴上来

wal

m062x好像也不支持BJ阻尼形式的色散校正，如果删了6-31g(d)还不行就把色散校正改成em=gd3

双黄蛋

怎么写了两个基组？

双黄蛋

wal 发表于 2024-10-22 14:58
m062x好像也不支持BJ阻尼形式的色散校正，如果删了6-31g(d)还不行就把色散校正改成em=gd3

之前试过明尼苏达确实不支持BJ

www...

wal 发表于 2024-10-22 14:58
m062x好像也不支持BJ阻尼形式的色散校正，如果删了6-31g(d)还不行就把色散校正改成em=gd3

谢谢解惑，我把6-31g（d）删了，还是跑不出来，em=gd3也改了，也还是跑不出来，我上传了计算失败的文件，麻烦大佬看看。

wal

www... 发表于 2024-10-22 20:05
谢谢解惑，我把6-31g（d）删了，还是跑不出来，em=gd3也改了，也还是跑不出来，我上传了计算失败的文件， ...

geom=allcheck Guess=read
你新加的这两个命令的意思是从%chk指向的检查点文件里读取坐标和初猜信息，你如果写了就不用再提供几何坐标了，但是%chk指向的文件必须要对。可以是你结构优化之后产生的chk文件。
如果你是想正常用gjf文件里的坐标去计算，就把这两个关键词删了

另外，你把geom=connectivity删了，几何坐标后面的键连关系就没用了，也给删了再跑试试

sobereva

%chk指定的路径不存在
提示写明了
Structure from the checkpoint file:  "/public1/home/t0s000841/software-t0s000841/gaussian16-test/HOLinputE1.chk"
FileIO operation on non-existent file.

别人回复后，不要编辑主贴，否则别人之前回复的内容全都作废，还会令其它读者迷惑，好像之前的人回答得不对似的

如果死活搞不懂，索性甭写geom=allcheck Guess=read，或者chk文件路径换个短点的免得写错。


根据前面别人的回复，输入文件里的硬伤太多，建议参加北京科音初级量子化学培训班（http://www.keinsci.com/workshop/KEQC_content.html）一次性系统地学个明明白白，免得踩各种坑白浪费时间，或者任务跑起来了实际上结果都是错的还浑然不知。

www...

sobereva 发表于 2024-10-23 12:55
%chk指定的路径不存在
提示写明了
Structure from the checkpoint file:  "/public1/home/t0s000841/sof ...

好的，谢谢sobereva老师