24小时热门版块排行榜    

查看: 1508  |  回复: 7

yoyo103108

铜虫 (小有名气)

[求助] 优化收敛失败 已有4人参与

使用g09计算,总是出现error in  internal coordinate system,怎么调整结构都优化不出来,请各位给我支支招吧~
回复此楼

» 猜你喜欢

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

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

wbf3ng

金虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +1
说清楚点,用的什么方法基组,什么体系?
一般先用低精度的初步优化然后再精细算

[ 发自小木虫客户端 ]
2楼2014-07-16 12:54:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
试一试opt=cartesian
3楼2014-07-16 13:01:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sutingting

铜虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
zhou2009: 金币+2 2014-07-17 21:12:47
通常不收敛的情况有好几种,你应该查看一下你的输出文件末行,看一下是哪一种错误,诸如L103,L502,L9999等,小木虫上有很多关于解决这类问题的。http://muchong.com/bbs/viewthread.php?tid=1079428&fpage=1&target=blank
加油
4楼2014-07-17 08:24:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

卡开发发

专家顾问 (著名写手)

Ab Initio Amateur

小基组下先opt=gdiis做个预先优化,也许能够有改善。
不一定挂在论坛,计算问题问题欢迎留言。
5楼2014-07-17 09:24:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoyo103108

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by wbf3ng at 2014-07-16 12:54:45
说清楚点,用的什么方法基组,什么体系?
一般先用低精度的初步优化然后再精细算

我要优化一个过渡态,
%chk=2-3_03_routeA1_ub3lyp.chk
%mem=2800MB
%nprocshared=12
# opt=(calcfc,ts,noeigen,gdiis,maxstep=10) ub3lyp/6-31g(d) geom=connectivity

结果优化后:GradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGradGrad
Berny optimization.
Using GDIIS optimizer.
NTrRot=    -1 NTRed=   497 NAtoms=    78 NSkip=   269 IsLin=F
Error in internal coordinate system.
Error termination via Lnk1e in /public/home/ciac_wu_4/software//g09/l103.exe at Fri Jul 18 04:15:46 2014.
Job cpu time:  9 days  3 hours 25 minutes 19.2 seconds.
File lengths (MBytes):  RWF=  10114 Int=      0 D2E=      0 Chk=     52 Scr=      1
6楼2014-07-19 10:19:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)

【答案】应助回帖

引用回帖:
6楼: Originally posted by yoyo103108 at 2014-07-19 10:19:22
我要优化一个过渡态,
%chk=2-3_03_routeA1_ub3lyp.chk
%mem=2800MB
%nprocshared=12
# opt=(calcfc,ts,noeigen,gdiis,maxstep=10) ub3lyp/6-31g(d) geom=connectivity

结果优化后:GradGradGradGradGradGr ...

因为你用了内坐标,内坐标中的角度必须小于180度。

而你的结构里有些角度接近了180度,这些角度在优化过程中算某一些单点时可能会超过180度,所以报错。

在一开始的结构里可以看到某些很接近180度的角度。

你可以手动调整这些角度令它们离180度远一点。

或者手写内坐标,添加虚原子令180度变成90度。

或者直接在opt的命令加上cart使用cartesian坐标,但计算花的时间会比用内坐标多。
7楼2014-07-21 16:35:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yoyo103108

铜虫 (小有名气)

引用回帖:
7楼: Originally posted by fatpig8832 at 2014-07-21 16:35:00
因为你用了内坐标,内坐标中的角度必须小于180度。

而你的结构里有些角度接近了180度,这些角度在优化过程中算某一些单点时可能会超过180度,所以报错。

在一开始的结构里可以看到某些很接近180度的角度。
...

你好,能不能给我个您的QQ号,请教您一下~你能帮我看看我这个结构怎么改啊~刚接触guassian 不太会~谢谢~
8楼2014-07-30 11:01:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yoyo103108 的主题更新
信息提示
请填处理意见