24小时热门版块排行榜    

查看: 1036  |  回复: 6

ww1987

铁杆木虫 (小有名气)

[交流] 【求助】计算错误 帮忙看看 已有5人参与

我计算C80(D5d)构型的分子 计算方法用B3LYP/6-31G(d)水平  算到一半出现这种现象 请帮忙!
  
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB+HF-LYP) =  -3048.38575494     A.U. after  129 cycles
             Convg  =    0.7361D-06             -V/T =  2.0094
             S**2   =   0.0000
KE= 3.020073824779D+03 PE=-3.325628686430D+04 EE= 1.408728660716D+04
Convergence failure -- run terminated.
Error termination via Lnk1e in C:\G03W\l502.exe at Sat Jul 03 19:20:23 2010.
Job cpu time:  0 days 10 hours 15 minutes 11.0 seconds.
File lengths (MBytes):  RWF=    897 Int=      0 D2E=      0 Chk=     98 Scr=      1
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sunhaitao

木虫 (小有名气)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+1):鼓励交流!! 2010-07-05 15:44:04
不收敛,搜索论坛帖子
一点一点进步!
2楼2010-07-05 15:15:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heyo_123

至尊木虫 (职业作家)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1):谢谢讨论 2010-07-14 09:38:42
nosymm
3楼2010-07-05 15:43:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ter20

铁杆木虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1):谢谢讨论 2010-07-14 09:38:47
加iop(5/13=1)试试
4楼2010-07-05 17:21:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ter20

铁杆木虫 (正式写手)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
heyo_123(金币+2):鼓励交流!! 2010-07-05 17:34:41
网上的建议:
1. 检查是否有初始文件错误
在命令行中加入 %kJob L301 or %kJob L302
如果通过则一般初始文件ok。
常见初级错误:
a. 自旋多重度错误
b. 变量赋值为整数
c. 变量没有赋值或多重赋值
d. 键角小于等于0度,大于等于180度
e. 分子描述后面没有空行
f. 二面角判断错误,造成两个原子距离过近
g. 分子描述一行内两次参考同一原子,或参考原子共线
2. SCF(自洽场)不收敛
则一般是L502错误
省却情况做64个cycle迭代(G03缺省128 cycles)
a. 修改坐标,使之合理
b. 改变初始猜 Guess=Huckel 或其他的,看Guess关键词。
c. 增加叠代次数 SCFCYC=N (对小分子作计算时最好不要增加,很可能结构不合理)
d. iop(5/13=1)这样忽略不收敛,继续往下做。
3. 分子对称性改变
a. 修改坐标,强制高对称性或放松对称性
b. 给出精确的、对称性确定的角度和二面角。 如CH4的角度给到109.47122
c. 放松对称性判据 Symm=loose
d. 不做对称性检查 iop(2/16=1) (最好加这个选项)
iop(2/16=2) 则保持新的对称性来计算

4. Opt时收敛的问题
a. 修改坐标,使之合理
b. 增加叠代次数optcyc=N

5. 优化过渡态,若势能面太平缓,则不好找到。
iop(1/8=10) 默认30(下一个结构和该结构的差别0.3Å),可改成10。如果每一步都要用到小的步长,应该加opt(notrustupdate)

6. 在CI(组态)方法中
如QCISD(T),CCSD(T),CID方法中,省却最大循环50,若出错(L913错误)
解决方法:#P QCISD(maxcyc=N) 注:N≤512

7. 优化过渡态
opt=TS (给出过渡态)
opt=qst2 (给出反应物和产物)
opt=qst3 (给出反应物和产物和过渡态)
a. 用G03时的出错 opt=ts 必须加FC (force constant)
写法:opt=(TS, calcFc)
or opt=(TS,calchffc)计算HF力常数,对QCISD,CCSD等方法用;
or opt=(TS,modRedundant) (最好写这个)
b. 如果计算采用QCISD计算(不好计算FC)
则写为QCISD opt=(TS, calcHFFC) (用HF计算FC)

8. 无法写大的Scratch文件RWF
a. 劈裂RWF文件
%rwf=loc1,size1,loc2,size2,……..,locN,-1
b. 改变计算方法 MP2=Direct可以少占硬盘空间
c. 限制最大硬盘 maxdisk=N GB,****MB,有些系统写2GB会出错,可以写2000MB

9. FOPT出错
原因是变量数与分子自由度数不相等。 可用POPT 或直接用OPT

10. 优化过渡态只能做一个STEP
原因是负本征数目不对 添加 iop(1/11)=1或者noeigentest
5楼2010-07-05 17:24:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhaokaimrt

新虫 (小有名气)


小木虫(金币+0.5):给个红包,谢谢回帖交流
不知道对不起。
6楼2010-07-05 18:35:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

beyondboyy

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
aylayl08(金币+1):谢谢讨论 2010-07-14 09:39:34
Error termination via Lnk1e in C:\G03W\l502.exe at Sat Jul 03 19:20:23 2010.
中的l502表示RHF。vhf,rohf自洽场迭代错误。
心静自然成!!
7楼2010-07-13 19:08:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ww1987 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见