24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1710  |  回复: 9

XYLF

金虫 (小有名气)


[交流] 关于L502 ERROR

前阵子看到有人提出Gaussian03W计算过程中出现了L502错误提醒计算终止的时候,可以通过提取计算最后一个优化坐标重建input文件重新开始计算!
请问这样的处理可行么?
如果跳开这个SCF计算队最终结果会有什么影响?
常规的处理L502 error 的方法是什么?
请指导~!
回复此楼

» 猜你喜欢

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

» 抢金币啦!回帖就可以得到:

查看全部散金贴

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

独孤狼

金虫 (正式写手)


★ ★
XYLF(金币+1):谢谢参与
gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-06-08 16:58:56
加上scf=xqc
2楼2011-06-08 15:29:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

XYLF

金虫 (小有名气)


引用回帖:
Originally posted by 独孤狼 at 2011-06-08 15:29:36:
加上scf=xqc

XQC表示?
不是还可以用nosym关键字么?有什么算法上的区别呢?
3楼2011-06-08 17:57:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤狼

金虫 (正式写手)



gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-06-08 21:38:20
nosym是消除轨道对称性
xqc是当一级scf未收敛的情况,再加入额外的scf=qc
qc是使用二次收敛scf程序
4楼2011-06-08 20:32:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

XYLF

金虫 (小有名气)


引用回帖:
Originally posted by 独孤狼 at 2011-06-08 20:32:01:
nosym是消除轨道对称性
xqc是当一级scf未收敛的情况,再加入额外的scf=qc
qc是使用二次收敛scf程序

那如果直接截取最后优化坐标重新输入INPUT文件坐标这种方法是否可行?
5楼2011-06-08 21:23:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

伟明

铁杆木虫 (著名写手)


★ ★
XYLF(金币+1):谢谢参与
gmy1990(金币+1): 非常感谢你对论坛的贡献,奖励一下! 2011-06-08 22:34:39
重新输入我认为就相当于增加循环次数了,可以加上scf=tight有助于收敛,楼上的 方法据说可行,我加scf=tight就收敛了
6楼2011-06-08 22:20:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

XYLF

金虫 (小有名气)


引用回帖:
Originally posted by 伟明 at 2011-06-08 22:20:12:
重新输入我认为就相当于增加循环次数了,可以加上scf=tight有助于收敛,楼上的 方法据说可行,我加scf=tight就收敛了

我想知道一个问题,如果直接跳掉这个L502-ERROR直接取最后坐标重新计算的话确实是相当增大的优化次数,但是在这个体系的这个坐标点重新计算的话,会不会相当于忽视了一个潜在的不收敛函数关系,这种潜在的可忽视的不收敛关系在之后的计算中会不会不断积加导致最终结果的较大偏差?
7楼2011-06-08 23:15:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyuanhe211

铜虫 (小有名气)



XYLF(金币+1):谢谢参与
opt=gdiis
8楼2011-06-08 23:45:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

XYLF

金虫 (小有名气)


引用回帖:
Originally posted by liyuanhe211 at 2011-06-08 23:45:55:
opt=gdiis

加了~用的是opt=(gdiis,maxcycle=100)但是还是偶尔报错L502,这个更GDIIS这个关键词没有关系吧!gdiis不是用于大分子体系修正的一个手段么?
9楼2011-06-08 23:55:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yongleli

木虫 (正式写手)



XYLF(金币+1):谢谢参与
引用回帖:
Originally posted by XYLF at 2011-06-08 23:55:30:
加了~用的是opt=(gdiis,maxcycle=100)但是还是偶尔报错L502,这个更GDIIS这个关键词没有关系吧!gdiis不是用于大分子体系修正的一个手段么?

scf=verytight。
G09里边默认的是tight。

还有什么vshift的之类也可以用用。

还有integral=(grid=untrafine,Acc2e=11)之类的也可以试试。
10楼2011-06-09 12:07:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 XYLF 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见