应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 结构优化 l502错误
51/1返回列表
查看: 2971  |  回复: 4
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

aliya0719

铜虫 (初入文坛)

[求助] 结构优化 l502错误

本人是gaussian计算新手，请大家帮忙看看问题出在哪里！先谢谢了！
二价钴离子和柠檬酸盐（CoC6H5O7-）的络合物，采用HF/Lanl2dz 方法进行结构优化，计算出错。输入结构结构摘自有机晶体结构数据库CCDC，删掉了部分原子。
输出Log文件报错：
The electronic state of the initial guess is 2-A.
of initial guess= 0.7500
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Restarting incremental Fock formation.
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Rare condition: small coef for last iteration:  0.000D+00
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Density matrix is not changing but DIIS error= 1.82D-06 CofLast= 6.60D-02.
The SCF is confused.
Error termination via Lnk1e in /home/ww/g03/l502.exe at Tue Aug  9 11:45:36 2011.

输入gjf文件：
%Chk=CoL-1.chk
%mem=20MW
%nproc=4
# opt hf/lanl2dz pop=reg

CoL-1

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

   B1             1.54757447
   B2             1.55032174
   B3             1.45131974
   B4             1.56119553
   B5             1.54115604
   B6             1.51060848
   B7             1.29391149
   B8             1.24379651
   B9             1.24125108
   B10            1.29058423
   B11            1.21474654
   B12            1.37796941
   B13            1.07819405
   B14            1.08521286
   B15            1.01270630
   B16            1.07936249
   B17            1.08146036
   B18            1.99529094
   A1           114.45175190
   A2           108.89312205
   A3           108.17772783
   A4           110.06425575
   A5           115.00470971
   A6           115.06575010
   A7           127.24593992
   A8           116.01138338
   A9           117.20194942
   A10          126.23457020
   A11          115.57267368
   A12          109.42217991
   A13          106.94866997
   A14          106.08583530
   A15          110.31524379
   A16          108.76707262
   A17          115.32436236
   D1            47.55401305
   D2           167.84367388
   D3           -70.87120092
   D4           164.91367592
   D5           -29.39594907
   D6           178.20658815
   D7            70.88776878
   D8          -108.01302269
   D9          -144.74386740
   D10           37.54708974
   D11           28.13346663
   D12          -87.63407589
   D13          -36.34965302
   D14           45.64915492
   D15          -73.72963063
   D16          -16.09209920

结构进行稍微调整后，又出现了不收敛的问题。下面是输出的log文件报错：
Convergence criterion not met.
SCF Done:  E(UHF) =  -897.994913496     A.U. after  129 cycles
             Convg  =    0.2426D-07             -V/T =  2.0642
             S**2   =   0.7620
Annihilation of the first spin contaminant:
S**2 before annihilation     0.7620,   after     0.7501
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/ww/g03/l502.exe at Tue Aug  9 12:08:41 2011.
回复此楼

» 收录本帖的淘帖专辑推荐

气态分子化学动力学

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2011-08-29 15:00:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

julyfly

银虫 (小有名气)
引用回帖:
2楼: Originally posted by luoronggen at 2011-08-30 17:00:20:
单点不收敛，可以增加循环次数，或者强制收敛。可以再#栏加上
scf=(xqc,maxcon=128,maxcyc=512,conver=8)，也可以将conver=8改为conver=6降低一下收敛标准。

谢谢你的帮忙！
但是，加上关键词scf=(xqc,maxcon=128,maxcyc=512,conver=6) 计算了，还是出错。 由不收敛回到了之前的出错：
  The SCF is confused.
Error termination via Lnk1e in /home/ww/g03/l502.exe
不知道二价钴离子是不是电子结构有什么特别的地方，分子的自选多重度不是2呢？
因为我从CCDC查了多个结构，每个结构计算都是出现上面那个错误
赞一下 回复此楼
爱拼才会赢，新的一年，敢于拼搏的一年
3楼2011-08-31 08:58:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

luoronggen

木虫 (正式写手)

【答案】应助回帖

aliya0719(金币+5): 2011-08-31 09:05:17
引用回帖:
1楼: Originally posted by aliya0719 at 2011-08-29 15:00:50:
本人是gaussian计算新手，请大家帮忙看看问题出在哪里！先谢谢了！
二价钴离子和柠檬酸盐（CoC6H5O7-）的络合物，采用HF/Lanl2dz 方法进行结构优化，计算出错。输入结构结构摘自有机晶体结构数据库CCDC，删掉了部 ...

单点不收敛，可以增加循环次数，或者强制收敛。可以再#栏加上
scf=(xqc,maxcon=128,maxcyc=512,conver=8)，也可以将conver=8改为conver=6降低一下收敛标准。
赞一下 回复此楼
贵在虚心坚持。
2楼2011-08-30 17:00:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zenghua63

新虫 (初入文坛)

【答案】应助回帖

aliya0719(金币+5): B1到B6都改成了1.4。。。，仍旧是同一个错。 2011-08-31 14:45:04
你把B1到B6的键长改小再试试。
赞一下 回复此楼
4楼2011-08-31 09:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aliya0719

铜虫 (初入文坛)
引用回帖:
4楼: Originally posted by zenghua63 at 2011-08-31 09:06:23:
你把B1到B6的键长改小再试试。

先谢谢你！
把B1到B6都改成了1.4多，仍旧是同一个错。
Density matrix is not changing but DIIS error= 2.35D-04 CofLast=-4.37D-02.
The SCF is confused.
Error termination via Lnk1e in /home/ww/g03/l502.
赞一下 回复此楼
5楼2011-08-31 14:46:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭