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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 量化中出错的log文件
61/1返回列表
查看: 1812  |  回复: 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的回帖

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

beefly

专家顾问 (职业作家)

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

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

fwenling0506

木虫 (小有名气)
做个记号，嘿嘿
回复此楼
6楼2011-09-04 10:34:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 我家小屁孩 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭