24小时热门版块排行榜    

查看: 2326  |  回复: 5

RRAM

金虫 (小有名气)

[求助] 高斯优化l9999出现错误,请教各位大虾~ 已有5人参与

这是我的输入文件
%chk=l3.chk        
#T b3lyp/6-31G** OPT

l3                 

0 1
C                 -3.60813676   -1.01124205   -0.46947725
C                 -2.22741142   -0.81107030   -0.46900027
C                 -1.71037209    0.48426216   -0.47143006
C                 -2.57419774    1.58016586   -0.47553751
C                 -3.95458004    1.37996395   -0.47649321
C                 -4.47161709    0.08417430   -0.47306644
H                 -4.01557467   -2.03258018   -0.46711127
H                 -1.54692806   -1.67488229   -0.46539303
H                 -2.16629868    2.60139590   -0.47751197
H                 -4.63563800    2.24345903   -0.47973815
C                  3.38129393    0.02624117   -0.33983967
C                  4.76201927    0.22641292   -0.33936269
C                  5.27905860    1.52174538   -0.34179248
C                  4.41523295    2.61764908   -0.34589993
C                  3.03485065    2.41744717   -0.34685563
C                  2.51781360    1.12165752   -0.34342886
H                  2.97385602   -0.99509696   -0.33747369
H                  5.44250263   -0.63739908   -0.33575545
H                  4.82313201    3.63887912   -0.34787439
H                  2.35379269    3.28094225   -0.35010057
C                 -0.18632221    0.70532769   -0.47001595
C                  1.00243670    0.87775880   -0.46891294
C                 -5.99572044   -0.13652539   -0.47384560
S                  7.04066605    1.77696023   -0.33965423
H                  7.47975128    1.78572158    0.89453676
S                 -6.65234265   -1.79098412   -0.47015939
H                 -7.96139970   -1.74193138   -0.47812541
P                 -6.99200332    1.11747862   -0.47796131
H                 -7.83177723    2.17448749   -0.48143047   
出现了这样的错误Sum of Mulliken charges with hydrogens summed into heavy atoms =   0.00000
Electronic spatial extent (au):  <R**2>=          13696.9194
Charge=              0.0000 electrons
Dipole moment (field-independent basis, Debye):
    X=              1.5062    Y=             -0.7739    Z=              0.7853  Tot=              1.8666
Error termination request processed by link 9999.
Error termination via Lnk1e in d:\g09\l9999.exe at Sun Dec 28 00:05:51 2014.
Job cpu time:  0 days  8 hours 55 minutes 37.0 seconds.
File lengths (MBytes):  RWF=     77 Int=      0 D2E=      0 Chk=      7 Scr=      1
请教各位了,不胜感激
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ztl19780124

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
RRAM: 金币+2 2015-01-21 14:27:18
RRAM: 金币+2, 有帮助 2015-01-21 15:02:43
可以尝试
1.增大循环次数,比如maxcyc=500,
2.改变初始构型
3.减少步长,比如maxstep=10.
4.是否带电荷,即自旋多重度是否正确
我想第三个应该最有可能搞定,再不行就maxstep=1。
2楼2014-12-28 10:01:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhou2009

版主 (著名写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
RRAM: 金币+2, 有帮助 2015-01-21 15:03:12
看了输入文件,21C-22C为三键,3C-21C和22C-C16为单键,那三键的P电子必定要和苯上的P电子形成更大的共轭体系,怎么会是三键、单键分明呢?应形成大π键。
而且两个苯环是共平面还是垂直?也应分别算一下,比较一下二者的能量。当作为平面输入时,是不能自行优化为垂直的。
在方法上建议用RHF/6-31G*先算一下。
3楼2014-12-28 10:24:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chu5263821

铜虫 (小有名气)

【答案】应助回帖

★ ★
RRAM: 金币+2 2015-01-21 15:03:19
最后没有收敛
可以把计算错误的输出文件用guassview打开,再存为gjf文件,同样写好命令 再算一次
不悔梦归处 只恨太匆匆
4楼2014-12-31 10:48:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xieruobing

木虫 (著名写手)

【答案】应助回帖


RRAM: 金币+1, 有帮助 2015-01-21 15:03:36
Electronic spatial extent (au):  <R**2>=          13696.9194
Charge=              0.0000 electrons
可以提取优化结构,在基础上加大步长继续运算
5楼2015-01-10 11:51:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianchangke

金虫 (小有名气)

【答案】应助回帖


RRAM: 金币+1, 有帮助 2015-01-21 15:03:31
我计算中也经常遇见,就是还没算完就结束了的意思,把最后坐标结果截出来,用复制然后粘贴到新的gjf文件中继续提交计算就是了
生活就像剥洋葱,一次只能剥开一层有时还会流泪
6楼2015-01-13 08:43:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 RRAM 的主题更新
信息提示
请填处理意见