计算化学公社

标题: 含混合基组的多步输入文件的写法求助 [打印本页]

作者
Author: likjvau 时间: 2020-10-16 15:05
标题: 含混合基组的多步输入文件的写法求助
求助，初学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.

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

作者
Author: zjxitcc 时间: 2020-10-16 15:11
你的报错是写了geom=connectivity但没写连接关系引起的，可以全不写。把所有的geom=connectivity及坐标底下的连接关系删去，这些都没用。

作者
Author: likjvau 时间: 2020-10-16 15:33

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

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

作者
Author: 喵星大佬 时间: 2020-10-16 15:36
我看了一下啊，你这个Bi旁边总共配了一个O加一个变了形的二氧化硅，末端还加了H。。。。。

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

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

作者
Author: sobereva 时间: 2020-10-16 21:02
论坛里的学习资源已经算是海量了，好好利用首页google框

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

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

作者
Author: 天气很好 时间: 2020-12-29 18:04
chk文件可以没有吗？我一般不输入chk文件这一行

作者
Author: sobereva 时间: 2021-1-2 20:14

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)