| 查看: 1136 | 回复: 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 |
» 猜你喜欢
华南师范大学广州市光谱分析与功能探针重点实验室诚招青年英才和博士后岗位研究人员
已经有10人回复
-
求标准卡片Cu0.4In0.4Zn1.2S2 PDF#32-0340
已经有2人回复
-
物理化学论文润色/翻译怎么收费?
已经有115人回复
-
量子点电致发光器件（QLED）：蓝光量子点（QDs）和纳米氧化锌电子传输材料
已经有0人回复
-
七嗪类物质合成求助
已经有0人回复
-
需要一个圆偏振光诱导化学反应，请问使用圆偏振光作为光源，应该需要哪些设备搭建
已经有0人回复
-
半透性容器进行高湿或者低湿研究
已经有0人回复
-
求ENDNOTE破解软件安装包
已经有4人回复
-
河海大学 （211.双一流高校）钙钛矿光电实验室丁勇团队欢迎有志青年！
已经有3人回复
-
黑龙江省原子能研究院化学课题组招收2026级化学专业或化学相关专业硕士生
已经有0人回复
-
26博士申请
已经有0人回复
5楼2009-12-23 09:41:22
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
20