24小时热门版块排行榜    

CyRhmU.jpeg
查看: 825  |  回复: 6
当前主题已经存档。

liujl

铜虫 (初入文坛)

[交流] 【求助】irc的结构与ts的结构相近

对ts进行优化和频率计算正是所需进攻的mode而且只有一个虚频,但跑irc五个点后就收敛了。其结构与ts的结构相近
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 10-16 22:12
相近但不相等,IRC的结果并不是反应物或产物的最稳定性结构,应属亚稳定结构,也许还存在其它中间产物或其它过渡态。看文章可知,一个反应过程可能存在多个过渡态。
2楼2009-10-16 16:26:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mengfc

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖交流
跑irc五个点后就收敛了?这应该是你设置它跑几个点吧?
3楼2009-10-16 17:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liujl

铜虫 (初入文坛)

我没有设置maxpoints
irc(rcfc,reverse,stepsize=10,maxcyc=100)
跑完irc后的势能曲线基本上时平的,而且和过渡态的能量一样
4楼2009-10-17 09:26:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

majun04

木虫 (著名写手)

majun

★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+2,VIP+0):谢谢 10-17 10:17
把第五个点的结构取出来,优化,并算频率,看看是不是能量最小点,如果是能量最小点,才说明你的IRC是走到底了,如果不是的话,说明有问题。
乐观、执着、勤奋
5楼2009-10-17 09:57:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liujl

铜虫 (初入文坛)

在irc的输出文件中显示了minimu found。
因为第五个点的结构和能量都与过渡态差不多,所以再优化第五个点还有必要吗?
6楼2009-10-17 10:20:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yytsnake

至尊木虫 (职业作家)

★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
lei0736(金币+4,VIP+0):谢谢 10-17 10:22
引用回帖:
Originally posted by liujl at 2009-10-17 09:26:
我没有设置maxpoints
irc(rcfc,reverse,stepsize=10,maxcyc=100)
跑完irc后的势能曲线基本上时平的,而且和过渡态的能量一样

没做过IRC计算,但从LZ的论述可知只能算6个点,资料是这样描述的:
例如:
# MP2/6-311G** IRC(RCFC,Forward,Stepsize=N, maxpoints=M)geom=check guess=check
若无指定,缺省的MaxPoints是正、逆方向各计算6个点。
Stepsize=N,指定计算步长,默认为N=10,指定合适的步长和点数,是计算成与否的因素之一。
IRC某点处不收敛
      一般可以有三种解决办法:
      a. 用更为严格的收敛判据来进行IRC计算
IRC(RCFC, Tight,Stepsize=5,maxpoints=N+M)
      b. 改变IRC计算的步长, 小步长可以得到更为精确的反应路径。
IRC(RCFC, Stepsize=2,maxpoints=N+M)
       c. 先做IRC到出错点之前的某一点结束,再用CalcFc方式Restart. 输入如下:
IRC(CalcFC,Restart,Stepsize=5,maxpoints=N+M)
        如果上面方法都不行, 可以试试IRC(CalcAll)
7楼2009-10-17 10:20:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 liujl 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见