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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 结构优化 l502错误
51/1返回列表
查看: 3290  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6GojgJvkDudM 2026-08-16 3/150 2026-08-17 09:49 by xLVPIRuSvCUe
[基金申请] 今天系统多次维护，明天很可能放榜！ +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +32 医学老男孩 2026-08-13 74/3700 2026-08-16 21:05 by 医学老男孩
[基金申请] 快农历七夕节了，轻松一下，男人悄悄话，女施主请不要进来。 +3 Tide man 2026-08-14 3/150 2026-08-16 17:47 by jurkat.1640
[基金申请] 我的国基提前知道中了，可是同事的操作让我实在接受不了，怎么会有这样的人 +11 家与远方 2026-08-10 16/800 2026-08-16 10:28 by ray43
[基金申请] 有时候，自然基金真的不能太认真 （我的申报经验） +10 majunge000 2026-08-11 12/600 2026-08-16 08:18 by xli1984
[精细化工] 招聘 金属平磨液，抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 各位道友，我要去昆明玩几天，回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 是这周出结果还是下周出结果？ +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[基金申请] 重要来源：本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了，巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 2019年青年基金涵评意见，大家看看几个A，几个B? +11 Tide man 2026-08-11 11/550 2026-08-13 07:35 by 撸猫猫
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 确定了，国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
信息提示
关闭
请填处理意见
关闭