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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 量化中出错的log文件
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 1936  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

我家小屁孩

铜虫 (小有名气)

[求助] 量化中出错的log文件

虫虫们好：
有劳见多识广的虫虫帮俺分析一下此log文件出错的原因并怎样改正，单线态已顺利完成，在算五线态时出错了，多谢了！！！！
Gaussian 09:  EM64L-G09RevA.01  8-May-2009
                17-Aug-2011
******************************************
%chk=/home/lyy/dqs/perpen.chk
%mem=2000MB
%nprocshared=4
Will use up to    4 processors via shared memory.
---------------------------------------------------
# opt freq=noraman b3lyp/6-31g(d) geom=connectivity
---------------------------------------------------
1/14=-1,18=20,19=15,26=3,38=1,57=2/1,3;
2/9=110,12=2,17=6,18=5,40=1/2;
3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3;
4//1;
5/5=2,38=5/2;
6/7=2,8=2,9=2,10=2,28=1/1;
7//1,2,3,16;
1/14=-1,18=20,19=15/3(2);
2/9=110/2;
99//99;
2/9=110/2;
3/5=1,6=6,7=1,11=2,16=1,25=1,30=1,71=1,74=-5/1,2,3;
4/5=5,16=3/1;
5/5=2,38=5/2;
7//1,2,3,16;
1/14=-1,18=20,19=15/3(-5);
2/9=110/2;
6/7=2,8=2,9=2,10=2,19=2,28=1/1;
99/9=1/99;
------
perpen
------
Symbolic Z-matrix:
Charge =  0 Multiplicity = 5
     （中间部分略）
(A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A) (A)
                 (A) (A) (A) (A) (A) (A) (A) (A)
The electronic state of the initial guess is 5-A.
Initial guess = 0.0000 = 0.0000 = 2.0000 = 6.0000 S= 2.0000
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.
EnCoef did   100 forward-backward iterations
EnCoef did    11 forward-backward iterations
EnCoef did     3 forward-backward iterations
EnCoef did     6 forward-backward iterations
EnCoef did     5 forward-backward iterations
EnCoef did     3 forward-backward iterations
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
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -4473.53934912     A.U. after  129 cycles
             Convg  =    0.2980D-04             -V/T =  2.0026
= 0.0000 = 0.0000 = 2.0000 = 6.0847 S= 2.0169
= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation     6.0847,   after     6.0019
DiagDN has N=   676 LTot=  2060 but NE2=     3 cannot use DSYEVD.
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/lyy/g09/l502.exe at Thu Aug 18 02:36:02 2011.
Job cpu time:  0 days 18 hours  8 minutes 33.1 seconds.
File lengths (MBytes):  RWF=    427 Int=      0 D2E=      0 Chk=     11 Scr=      1
小女是实验为主，量化知识浅薄，希望能得到高手虫虫们的指点！
回复此楼

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

高斯日常学习总结

» 猜你喜欢

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

beefly

专家顾问 (职业作家)

地沟油冶炼专家
如果当前结构的能量很重要，必算不可，用scf=qc

否则，找个小点的基组或者换个泛函，只要能收敛就行，做个预优化。然后用这个结构作为初始，做正式的优化计算
赞一下 回复此楼
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
5楼2011-08-24 21:15:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
zhou2009(金币+2): 2011-08-24 19:40:30
我家小屁孩(金币+10): 不好意思，我是小木虫新手，不怎么上，发金币也晚了点。很感谢你的！ 2011-09-02 17:39:48
无他，构型不好而已...

1. 改构型
2. 可用 scfcyc=500 增大scf计算的cycle，默认的是129
（SCF Done:  E(UB3LYP) =  -4473.53934912     A.U. after  129 cycles）
3. scf=qc，用了这个，速度很慢；所以一般用scf=qc算第一点即可，然后再读取第一点的chk文件继续往下做...
赞一下(1人) 回复此楼
2楼2011-08-24 11:06:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bubuweiying

金虫 (著名写手)
★ ★
zhou2009(金币+2): 2011-08-24 19:40:47
引用回帖:
2楼: Originally posted by fatpig8832 at 2011-08-24 11:06:23:
无他，构型不好而已...

1. 改构型
2. 可用 scfcyc=500 增大scf计算的cycle，默认的是129
（SCF Done:  E(UB3LYP) =  -4473.53934912     A.U. after  129 cycles）
3. scf=qc，用了这个，速度很慢；所以一般 ...

如果还不行，只能拿出杀手锏了，就是使用qc，但不建议直接使用，而是使用xqc关键词，比如scf=(maxcycle=80,xqc)，意思是如果scf正常计算(dc)在80个循环之内不收敛才进行昂贵的qc计算，因为scf不收敛多数在几个优化的过程中出现，无法判断哪一步优化的时候会出现scf不收敛，所以用xqc比纯粹使用qc要省时的多。
使用scf=qc，这个计算会慢，而且需要用stable关键字来测试结果是否波函数稳定。如果这个还不收敛，会提示L508错误，请问如何测试波函数稳定啊，是加上关键词scf=xqc stable 就可以么，即# opt freq rb3lyp/6-31+g(d) pop=(nbo,savenbo,full) scf=xqc stable 谢谢啊
赞一下(1人) 回复此楼
春天花会开！
3楼2011-08-24 11:43:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

★ ★
zhou2009(金币+2): 2011-08-24 19:41:04
我家小屁孩(金币+5): 3ks 2011-09-12 12:50:20
引用回帖:
3楼: Originally posted by bubuweiying at 2011-08-24 11:43:49:
如果还不行，只能拿出杀手锏了，就是使用qc，但不建议直接使用，而是使用xqc关键词，比如scf=(maxcycle=80,xqc)，意思是如果scf正常计算(dc)在80个循环之内不收敛才进行昂贵的qc计算，因为scf不收敛多数在几个优 ...

测试波函数稳定否，用stable=opt，只需算一个单点即可，不用算freq之类

例如
# hf/6-31g* stable=opt

若用scf=qc，也可以写成
# hf/6-31g* stable=opt scf=qc

如要在波函数稳定的基础上继续优化或算freq，可用--link1--
# hf/6-31g* stable=opt scf=qc
...
...
...
空一行
--link1--
%chk=...跟上面的一致...
# hf/6-31g* opt freq=noraman
guess=read geom=allcheck
赞一下(1人) 回复此楼
4楼2011-08-24 12:10:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 334分 一志愿武理 材料求调剂 +7 李李不服输 2026-03-26 7/350 2026-03-29 16:25 by 学员8dgXkO
[考研] 南京大学化学调剂 +6 景随风 2026-03-29 9/450 2026-03-29 15:27 by 唐沐儿
[考研] 300求调剂，材料科学英一数二 +9 leaflight 2026-03-24 9/450 2026-03-29 14:18 by wxiongid
[考研] 329求调剂 +10 钮恩雪 2026-03-25 10/500 2026-03-29 13:32 by peike
[考研] 309求调剂 +7 谁不是少年 2026-03-29 7/350 2026-03-29 13:17 by mumin1990
[考研] 本科双非材料，跨考一志愿华电085801电气，283求调剂，任何专业都可以 +6 芝士雪baoo 2026-03-28 8/400 2026-03-29 08:16 by 松花缸1201
[考研] 11408软件工程求调剂 +3 Qiu学ing 2026-03-28 3/150 2026-03-28 21:50 by zhq0425
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考研] 0856调剂 +5 求求让我有书读� 2026-03-26 6/300 2026-03-27 15:12 by caszguilin
[考研] 一志愿陕师大生物学071000，298分，求调剂 +5 SYA！ 2026-03-23 5/250 2026-03-27 09:29 by 不吃魚的貓
[考研] 351求调剂 +4 麦克阿磊 2026-03-24 4/200 2026-03-27 00:32 by wxiongid
[考研] 材料调剂 5+4 想要一壶桃花水 2026-03-25 10/500 2026-03-26 19:56 by 不吃魚的貓
[考研] 一志愿哈工大，085400，320，求调剂 +4 gdlf9999 2026-03-24 4/200 2026-03-25 23:01 by boxking200
[考研] 0854AI CV方向招收调剂 +4 章小鱼567 2026-03-23 4/200 2026-03-25 17:04 by CoderLoser
[考研] 296求调剂 +4 汪！？！ 2026-03-25 7/350 2026-03-25 16:41 by 汪！？！
[考研] 调剂 +4 13853210211 2026-03-24 4/200 2026-03-24 19:44 by ms629
[考研] 化工专硕求调剂 +3 question挽风 2026-03-24 3/150 2026-03-24 18:48 by jhhcooi
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
信息提示
关闭
请填处理意见
关闭