| 查看: 1083 | 回复: 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 |
» 猜你喜欢
推荐给教师的一本单词书《金鱼单词讲义：从26个拉丁字母到106万个英语单词》
已经有32人回复
-
北京工业大学材料科学与工程学院《院士团队》招收工程专业博士研究生（计算优先）
已经有1人回复
-
物理化学论文润色/翻译怎么收费?
已经有158人回复
-
核磁分析软件MestReNova打开文件时报错
已经有0人回复
-
在职博后不能申请博后基金了，那么在职博后意义何在？
已经有2人回复
-
青岛大学化学化工学院分子测量学研究院2026年招收博士研究生
已经有0人回复
-
香港科技大学（广州）诚招电催化方向博士生（2026秋入学）
已经有0人回复
-
求助Cu2+1O的CIF文件(PDF: 05-0667)
已经有0人回复
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
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
5楼2009-12-23 09:41:22