含混合基组的多步输入文件的写法求助

likjvau

求助，初学Gaussion，在计算掺杂分子结构时，想要分布进行激发态计算，先使用B3LYP，再使用PBE1PBE，无奈手中资源太少，想请教一下大家，多步输入的输入文件怎么写：
下面是我写的，运行报错的输入文件

%chk=\home\test\g09\Bi.chk
#p td=(singlets,nstates=25) b3lyp/genecp geom=connectivity

Title Card Required

0 1
Si                -0.09180200    3.33846800   -0.03012100
Si                 1.41093900    0.73113200    0.14736800
Si                -1.60303700    0.75316700   -0.04365600
O                 -1.42808000    2.38540500   -0.17749700
O                  1.22810500    2.36775100    0.12450200
O                 -0.11082800    0.08155800    0.12297600
O                 -2.55396600    0.49860900    1.27796900
H                 -2.70453400   -0.41349200    1.54939800
O                  2.22637000    0.22637400   -1.19542700
O                  2.20222400    0.32164500    1.52986700
H                  2.81405400   -0.42557000    1.45552400
O                 -2.25589800    0.08946100   -1.40405300
H                 -3.15976400    0.33839800   -1.62660500
O                 -0.33410100    4.27026100    1.30708200
H                  0.40730500    4.80068300    1.61938400
O                  0.16089000    4.25025600   -1.37821500
H                 -0.47316600    4.95185200   -1.56262500
Bi                 3.67435600   -1.19468700   -1.43486100

1 4 1.5 5 2.0 14 1.5 16 2.0
2 5 2.0 6 1.5 9 1.5 10 2.0
3 4 1.5 6 2.0 7 1.5 12 1.5
4
5
6
7 8 1.0
8
9 18 1.0
10 11 1.0
11
12 13 1.0
13
14 15 1.0
15
16 17 1.0
17
18

H O Si 0
6-31+G**
****
Bi     0
LANL2DZ
****

Bi     0
LANL2DZ



--Link1--
%chk=\home\test\g09\Bi+.chk
#p td=(singlets,nstates=25) pbepbe/genecp geom=connectivity guess=read geom=allcheck

H O Si 0
6-31+G**
****
Bi     0
LANL2DZ
****

Bi     0
LANL2DZ





geom=allcheck已经表示从chk文件中读取电荷、多重度、分子结构了，所以都没写，第二个任务中也按要求只把基组写上去了
但运行报错，输出如下：

%chk=\home\test\g09\Bi+-td25-08d-R.chk
----------------------------------------------------------------------
#p td=(singlets,nstates=25) pbepbe/genecp geom=connectivity guess=read
  geom=allcheck
----------------------------------------------------------------------
1/29=7,38=1,57=2/1;
2/12=2,40=1/2;
3/5=7,11=9,14=-4,16=1,17=8,25=1,30=1,74=1009,116=-2/1,2,8,3;
4/5=1/1;
5/5=2,38=6/2;
8/6=1,10=1,108=25/1;
9/41=25,42=1,48=1,70=2/14;
6/7=2,8=2,9=2,10=2/1;
99/5=1,9=1/99;
Leave Link    1 at Fri Oct 16 14:04:37 2020, MaxMem=           0 cpu:         0.2
(Enter /home/test/g09/l101.exe)
Structure from the checkpoint file:  ""
-------------------
Title Card Required
-------------------
Charge =  0 Multiplicity = 1
Z-Matrix found in file.
Si,0,-0.091802,3.338468,-0.030121
Si,0,1.410939,0.731132,0.147368
Si,0,-1.603037,0.753167,-0.043656
O,0,-1.42808,2.385405,-0.177497
O,0,1.228105,2.367751,0.124502
O,0,-0.110828,0.081558,0.122976
O,0,-2.553966,0.498609,1.277969
H,0,-2.704534,-0.413492,1.549398
O,0,2.22637,0.226374,-1.195427
O,0,2.202224,0.321645,1.529867
H,0,2.814054,-0.42557,1.455524
O,0,-2.255898,0.089461,-1.404053
H,0,-3.159764,0.338398,-1.626605
O,0,-0.334101,4.270261,1.307082
H,0,0.407305,4.800683,1.619384
O,0,0.16089,4.250256,-1.378215
H,0,-0.473166,4.951852,-1.562625
Bi,0,3.674356,-1.194687,-1.434861
Wanted an integer as input.
Found a string as input.
H O Si 0                                                                        
?
Error termination via Lnk1e in /home/test/g09/l101.exe at Fri Oct 16 14:04:37 2020.




希望好心人能帮我指点一下，很简单的问题，自学实在是缺乏资料。。

zjxitcc

你的报错是写了geom=connectivity但没写连接关系引起的，可以全不写。把所有的geom=connectivity及坐标底下的连接关系删去，这些都没用。

likjvau

zjxitcc 发表于 2020-10-16 15:11
你的报错是写了geom=connectivity但没写连接关系引起的，可以全不写。把所有的geom=connectivity及坐标底下 ...

已经根据您的建议修改，计算成功了，谢谢指点

喵星大佬

我看了一下啊，你这个Bi旁边总共配了一个O加一个变了形的二氧化硅，末端还加了H。。。。。

应该不能描述你想算得体系吧。。。。

如果是周期性的Bi参杂SiO2，建议趁早QE/VASP/Abinit/CASTEP

sobereva

论坛里的学习资源已经算是海量了，好好利用首页google框

LANL2DZ给Bi用完全莫名其妙，比其它元素差太远了
直接都用def2-SVP就完事了

如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，避免有任何歧义，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“含混合基组的多步输入文件的写法”改了，以后务必注意

天气很好

chk文件可以没有吗？ 我一般不输入chk文件这一行

sobereva

天气很好 发表于 2020-12-29 18:04
chk文件可以没有吗？ 我一般不输入chk文件这一行

当然可以
但不留着chk文件，以后做波函数分析、续算都没法算