24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 3079  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

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.
回复此楼
已阅   回复此楼   关注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的回帖
查看全部 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的回帖

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 一志愿郑州大学材料与化工085600,求调剂 +34 吃的不少 2026-04-02 34/1700 2026-04-07 20:01 by lrll?l
[考研] 306求调剂 +3 15287505595 2026-04-03 3/150 2026-04-07 18:08 by 蓝云思雨
[考研] 081700学硕,323分,一志愿中国海洋大学求调剂学校 +19 披星河 2026-04-04 19/950 2026-04-07 15:00 by 上岸快快
[考研] 305分求调剂 +3 哈_哈_哈_哈_哈 2026-04-04 5/250 2026-04-07 14:49 by 哈_哈_哈_哈_哈
[考研] 328求调剂 +4 ghhh88888 2026-04-06 5/250 2026-04-07 14:45 by ghhh88888
[考研] 材料工程302分求调剂 +13 zyx上岸! 2026-04-04 13/650 2026-04-07 11:14 by 诗与自由
[考研] 319分085702安全工程求调剂 +6 rious 2026-04-05 6/300 2026-04-07 09:42 by jp9609
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 09:25 by 小黑不怕难
[考研] 08600生物与医药-327 +9 18755400796 2026-04-05 9/450 2026-04-06 22:35 by 52305043001
[考研] 专硕304找调剂,一线城市最好 +3 李lsl李 2026-04-05 3/150 2026-04-06 12:16 by ffffjjjj
[考研] 求调剂,一志愿郑州大学材料与化工专硕,英二数二342分,求老师收留 +19 v12abo 2026-04-02 21/1050 2026-04-06 09:29 by 蓝云思雨
[考研] 调剂 一志愿吉林大学357分 +5 .Starry. 2026-04-04 5/250 2026-04-06 09:28 by cql1109
[考研] 一志愿武汉理工大学-085601材料工程(专硕)-总分353求调剂 +3 2626262626li 2026-04-02 3/150 2026-04-06 09:08 by 无际的草原
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 313求调剂 +5 海日海日 2026-04-04 5/250 2026-04-05 15:52 by jndximd
[考研] 0835学硕299求调剂 08大类可接受 +5 useryy 2026-04-03 5/250 2026-04-04 20:07 by 蓝云思雨
[考研] 342求调剂 +3 Liang7111 2026-04-04 5/250 2026-04-04 19:47 by dongzh2009
[考研] 专硕085601求调剂 +7 suyifei 2026-04-03 8/400 2026-04-03 14:00 by 欣喜777
[考研] 379求调剂 +3 ?苦瓜不苦 2026-04-01 3/150 2026-04-01 20:09 by 暮云清寒
信息提示
请填处理意见