24小时热门版块排行榜    

查看: 571  |  回复: 8
当前主题已经存档。
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

皓儿165

铜虫 (正式写手)

[交流] 【求助】两种方法算为什么有一种不收敛

我开始用的DFT算的,结果很好,正常结束,但是刚刚用MP2算,却不收敛,结果是这样的:Annihilation of the first spin contaminant:
S**2 before annihilation     0.8753,   after     0.7529
Convergence failure -- run terminated.
Error termination via Lnk1e in /home/g03/g03/l502.exe at Wed Jun  4 10:39:13 2008.
Job cpu time:  0 days  0 hours 29 minutes 43.7 seconds.
File lengths (MBytes):  RWF=     48 Int=      0 D2E=      0 Chk=      9 Scr=      1
求高手帮忙改正一下
回复此楼
做个聪明优雅的女人
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓儿165

铜虫 (正式写手)

引用回帖:
Originally posted by hairan at 2008-6-5 23:35:


直接用scf=qc,计算量太大

lz可以尝试使用scf=(maxcycle=80,xqc)
那个maxcycle之后的那个数值可以自己定,建议在50--100之间,如果这么多循环之后scf还不收敛,则跳转到qc做计算,如果在循环数之内收敛了 ...

万分感谢啊
做个聪明优雅的女人
9楼2008-06-06 09:39:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 9 个回答

jghe

金虫 (著名写手)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,化学区正在举办十大“耐思”版主评选,http://emuch.net/bbs/viewthread.php?tid=845527,欢迎你去投票支持计模版版主!
这种现象很正常啊,L502不收敛,改变结构或加一些关键字如,nosymm, scf=maxcycle=500,等等
2楼2008-06-04 20:06:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

皓儿165

铜虫 (正式写手)

可是以前用DFT算的收敛啊,是不是结构应该没多大问题,我也试过scf=maxcycle=500了,还是出现同样的错误。另外还有个问题,在文献中看得G3(MP2)-RAD方法是否就是指的G3MP2?
做个聪明优雅的女人
3楼2008-06-04 21:47:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lsl6585

银虫 (初入文坛)

★ ★
zzgyb(金币+2,VIP+0):谢谢你的参与,化学区正在举办十大“耐思”版主评选,http://emuch.net/bbs/viewthread.php?tid=845527,欢迎你去投票支持计模版版主!
引用回帖:
Originally posted by 皓儿165 at 2008-6-4 10:25:
我开始用的DFT算的,结果很好,正常结束,但是刚刚用MP2算,却不收敛,结果是这样的:Annihilation of the first spin contaminant:
S**2 before annihilation     0.8753,   after     0.7529
Convergence f ...

是在Linux系统上面跑的吗?不收敛的原因都很多种,你可以把grep "SCF D" *.log和egrep 'out of | SCF Done | Converged | NO | YES | exceeded' *.log | grep -v '\\\\'这两条命令的结果贴上来看一下,我们可以再一起分析一下,注意把这两条命令里面的*.log换成你自己的输出文件的名字
4楼2008-06-05 00:32:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见