24小时热门版块排行榜    

查看: 741  |  回复: 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的回帖

erylingjet

铁杆木虫 (著名写手)

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

lixiaona158

木虫 (著名写手)

自由人

★ ★
lei0736(金币+2,VIP+0):谢谢
如果有guess=read geom=connectivity
那么下面就不应该输入坐标
而且 geom=connectivity 我不知道什么意思  我计算的时候加上会导致比收敛
2楼2008-05-19 15:43:31
已阅   回复此楼   关注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的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见