24小时热门版块排行榜    

查看: 737  |  回复: 8
当前主题已经存档。

36du

银虫 (正式写手)

[交流] [求助]gaussian计算又出错,受伤!

哪为高手知道这是怎么回事啊?我算得是一个镁离子和三个甲醇的体系,电荷是2,结构优化断掉了,是我的初始结构输入的不多吗?还是其他那里错了?
%Chk=Untitled-1.chk
%mem=400MB
# opt rob3lyp/6-31+g guess=read geom=connectivity

CH3OH

2 1
C
O                  1              B1
H                  1              B2    2              A1
H                  1              B3    2              A2    3              D1
H                  1              B4    2              A3    3              D2
H                  2              B5    1              A4    3              D3
H                  2              B6    1              A5    3              D4
C                  2              B7    1              A6    3              D5
O                  8              B8    2              A7    1              D6
H                  8              B9    2              A8    1              D7
H                  8             B10    2              A9    1              D8
H                  9             B11    8             A10    2              D9
H                  2             B12    1             A11    9             D10
C                  8             B13    2             A12    1             D11
O                 14             B14    8             A13    2             D12
H                 14             B15    8             A14    2             D13
H                 14             B16    8             A15    2             D14
H                 15             B17   14             A16    8             D15
Mg                 9             B18    8             A17    2             D16

输出:
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(ROB+HF-LYP) =  -533.134008776     A.U. after  129 cycles
             Convg  =    0.1663D-01             -V/T =  1.9943
             S**2   =   0.0000
Annihilation of the first spin contaminant:
S**2 before annihilation     0.0000,   after     0.0000
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe at Thu May 08 21:40:37 2008.
Job cpu time:  0 days  1 hours 32 minutes 55.0 seconds.
File lengths (MBytes):  RWF=     22 Int=      0 D2E=      0 Chk=     10 Scr=      1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lixiaona158

木虫 (著名写手)

自由人

★ ★
lei0736(金币+2,VIP+0):谢谢
如果有guess=read geom=connectivity
那么下面就不应该输入坐标
而且 geom=connectivity 我不知道什么意思  我计算的时候加上会导致比收敛
2楼2008-05-19 15:43:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

erylingjet

铁杆木虫 (著名写手)

个人认为还是构型的问题,建议合理修改!
3楼2008-05-19 19:47:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huahua1981

铜虫 (小有名气)


lei0736(金币+1,VIP+0):谢谢
同意1楼的 命令不对
你的应该是 步数少了 加上optcyc=300 试试
4楼2008-05-20 08:46:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

suny2005

木虫 (小有名气)


lei0736(金币+1,VIP+0):谢谢
收敛失败,可以加scf=qc
不过楼主的输入文件也有问题应该去掉guess=read geom=connectivity
并且优化的时候应该算频率freq
5楼2008-05-20 08:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brover

捐助贵宾 (正式写手)

学术刺客


lei0736(金币+1,VIP+0):谢谢
结构坐标后面有连接表吗?  没有的话把geom=connectivity去掉试试吧
6楼2008-05-20 12:37:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jjf_sxnu

铜虫 (正式写手)


lei0736(金币+1,VIP+0):谢谢
看一看我的贴子吧:《【共享资料】How to handle a SCF convergence fail (Link 502)》
7楼2008-05-20 22:35:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xxylf

金虫 (正式写手)

请参考6楼的意见
8楼2008-05-21 07:42:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hairan

木虫 (著名写手)

瞌睡木虫

★ ★ ★
lei0736(金币+3,VIP+0):谢谢
2到5楼说得都不对
这个问题我想是rob3lyp关键词中的ro导致的
做限制性开壳层(restricted-open, ro)计算常常会遇到收敛的问题,lz先做一个rb3lyp试试看,感觉你的体系不需要要用ro的,自旋污染很严重吗?
guess=read geom=connectivity两个关键词和收敛不收敛没有必然关系。
据权威报纸调查,84%的青少年犯罪者有玩电子游戏的经历,故应该禁止经营电子游戏。据我们所知,100%的青少年犯罪者有穿鞋的经历,故制鞋厂应该停业。
9楼2008-05-21 16:06:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 36du 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见