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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 结构优化 l502错误
51/1返回列表
查看: 3043  |  回复: 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的回帖

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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 aliya0719 的主题更新
51/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 286求调剂 +3 lemonzzn 2026-03-16 5/250 2026-03-16 20:43 by lemonzzn
[文学芳草园] 伙伴们，祝我生日快乐吧 +17 myrtle 2026-03-10 26/1300 2026-03-16 18:32 by 青橙Ln
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7 Ncdx123456 2026-03-13 8/400 2026-03-16 12:15 by karry wen
[考研] 344求调剂 +3 knight344 2026-03-16 3/150 2026-03-16 09:42 by 无际的草原
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 274求调剂 +4 时间点 2026-03-13 4/200 2026-03-15 15:29 by Rambo13
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 255求调剂 +3 李嘉慧， 2026-03-12 4/200 2026-03-14 16:58 by 有只狸奴
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 0805，333求调剂 +3 112253525 2026-03-10 3/150 2026-03-13 23:42 by JourneyLucky
[考研] 材料与化工（0856）304求B区调剂 +6 邱gl 2026-03-12 7/350 2026-03-13 23:24 by 邱gl
[考研] 279求调剂 +3 Dizzy123@ 2026-03-10 3/150 2026-03-13 23:02 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 求b区学校调剂 +3 周56 2026-03-11 3/150 2026-03-13 16:20 by JourneyLucky
[考研] 【0856】化学工程（085602）313 分，本科学科评估A类院校化学工程与工艺，诚求调剂 +7 小刘快快上岸 2026-03-11 7/350 2026-03-13 16:06 by ruiyingmiao
[考研] 工科材料085601 279求调剂 +8 困于星晨 2026-03-12 10/500 2026-03-13 15:42 by ms629
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生，2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
关闭
请填处理意见
关闭