| 查看: 1111 | 回复: 4 | ||||
| 当前主题已经存档。 | ||||
[交流]
【求助】gaussian IRC 问题
|
||||
|
输入文件如下： %chk=abc #IRC=(forward,phase=(1,2),calcfc,maxcycle=100,stepsize=7,maxpoints=300) ub3lyp+G(d,p) scf=xqc iop(2/16=3)geom=allcheck guess=read 在我的abc.chk文件中，是优化好的做了频率分析的TS结构 文件跑了7个小时，可以用gaussian view看出最后的IRC跑到的结构就是我要的minimal point，就差最后的几步。结果因为 Number of steps exceeded,NStep=100,而跑死了 我试过其他的步长，shepsize=7，9，10，11，12，13，14，15都不行 哪位大侠给指点一下，谢谢啦 |
回复此楼» 收录本帖的淘帖专辑推荐
 Gaussian |
» 猜你喜欢
304求调剂
已经有5人回复
-
1
已经有21人回复
-
物理化学论文润色/翻译怎么收费?
已经有161人回复
-
0703 调剂
已经有5人回复
-
上海华东理工大学机械与动力工程学院先进材料与制造研究组招聘公告
已经有0人回复
-
304求调剂
已经有6人回复
-
西南某211招27届博士
已经有30人回复
-
大连工业杰青/长江团队-生物质材料-储能电池方向招收2026级博士生
已经有0人回复
-
304求调剂
已经有1人回复
-
中药学306，求调剂
已经有0人回复
hzfish
金虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1515.8
- 散金: 10
- 红花: 3
- 帖子: 604
- 在线: 358.1小时
- 虫号: 753222
- 注册: 2009-04-20
- 性别: GG
- 专业: 理论和计算化学
2楼2009-12-22 20:49:11
yytsnake
至尊木虫 (职业作家)
- 应助: 357 (硕士)
- 贵宾: 0.01
- 金币: 18204.3
- 散金: 5
- 红花: 36
- 帖子: 3189
- 在线: 444.6小时
- 虫号: 92066
- 注册: 2005-09-08
- 性别: GG
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+2,VIP+0):感谢转载 12-23 09:40
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+2,VIP+0):感谢转载 12-23 09:40
|
资料上如是说： IRC计算的选项 IRC计算的坐标系选项 MassWeighted 质量加权坐标 Internal 内坐标 Cartersian 直角坐标 IRC计算的力常数选项 RCFC 从Checkpoint读取直角坐 标力常数矩阵 CalcFC 第一点计算力常数矩阵 CalcAll 所有点都计算力常数矩阵(少用，耗时) IRC计算的路径选项： Forward 正向IRC计算 Reverse 逆向IRC计算 Stepsize=N IRC计算的步长 单位是0.01 amu1/2Bohr，缺省是10 Maxpoints=M IRC计算的点数， 缺省时6 你弄Maxpoints=300，是不是要从最高点到终点计算300 个点？即从过渡态到近产物区产生300个结构变化？是不是嫌计算机有点懒？ [ Last edited by yytsnake on 2009-12-22 at 21:27 ] |
3楼2009-12-22 21:26:04
llhxaustin
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 4494.5
- 红花: 7
- 帖子: 2749
- 在线: 174小时
- 虫号: 143679
- 注册: 2005-12-22
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 12-23 09:40
小木虫(金币+0.5):给个红包，谢谢回帖交流
yjcmwgk(金币+3,VIP+0): 12-23 09:40
| 哦，Number of steps exceeded,NStep=100，说明在走IRC路径上的一个步长时，已经超过100圈优化了，即，上一点你已经找到，而找下一点时已优化了100圈。此时程序跳出是因为超出了默认步数。个人觉得，没必要在做下去了，从上一点优化得到后，如果20圈内得不到下一点，一般优化的几何就已经远远偏离IRC路径了，所以，即使你得到了一个收敛的几何，一般也没有太多意义。我以前做时，一般都设定optcyc=20,对于原子数较多的，也就设定optcyc=40。如果有可能，建议你用内坐标优化一下TS，再用内坐标方式做IRC，收敛性要好些。 |
4楼2009-12-23 08:37:07
5楼2009-12-23 09:41:22
