24小时热门版块排行榜    

查看: 3184  |  回复: 12

yurivivi

新虫 (小有名气)

[求助] ts找到过渡态,可以把做irc两端的点提出来优化作为反应物和产物的稳定结构么? 已有2人参与

已知反应物和产物的结构,用的是ts的办法找到的过渡态,然后做irc各做了40个点,发现确实往反应物和产物的方向跑了。所以停掉,用各自的最后一个点拿出来再优化,得到了反应物和产物的结构。但是这个结构和我之前单独优化的结构构型是有差别的,角度啊,二面角什么的。请问各位前辈可以用irc的结果这么找反应物和产物么?
ps:算的双分子反应,所以是把A+B放在作为一个输入文件优化的,上面用irc的点优化走出来的构型和之前单独优化的构型不一样。(这个构型主要指两个分子的相对空间位置不一样,A分子和B分子的构型都是没差别的。)还有一个问题:用不同的方法和不同的基组,优化出来的A+B的构型也不一样,这个是正常的么?
希望各位前辈给予指点,完全的新手,也没人指导,都是自己瞎琢磨的。
回复此楼

» 收录本帖的淘帖专辑推荐

计算化学

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

清新俊逸

木虫 (正式写手)

A+B放在一起优化出来的不是反应物,而是势垒前面的一个范德华势阱。

[ 发自手机版 http://muchong.com/3g ]

» 本帖已获得的红花(最新10朵)

8楼2014-04-21 18:09:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yurivivi

新虫 (小有名气)

送红花一朵
引用回帖:
8楼: Originally posted by 清新俊逸 at 2014-04-21 18:09:41
A+B放在一起优化出来的不是反应物,而是势垒前面的一个范德华势阱。

那么请问,这样优化是合理的么?还是各自分开来优化呢?我主要想去计算反应的活化能。
9楼2014-04-23 10:40:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

wangsihang

木虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
yurivivi: 金币+3, 有帮助 2014-04-23 10:41:10
不同的基组和方法优化出来的构型会有差异的。如果仅仅是空间位置有差异或者角度有所不同应该是可以忽略的。
2楼2014-04-19 22:09:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yypcnbe

木虫 (小有名气)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
yurivivi: 金币+2, 有帮助 2014-04-23 10:41:19
那个两个分子的相对空间位置不一样,你是不是指的是互为镜像结构,可以在gview里面转换镜像
3楼2014-04-20 17:39:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yurivivi

新虫 (小有名气)

引用回帖:
3楼: Originally posted by yypcnbe at 2014-04-20 17:39:38
那个两个分子的相对空间位置不一样,你是不是指的是互为镜像结构,可以在gview里面转换镜像

不是镜像的结构。比如A+B-ts-C+D:A与B搁在一个输入文件里优化之后,用gview打开来,有一个构型。C和D也同样放在一个输入文件里优化,然后也优化出一个最小能量的构型。然后,我是用TS的方法找的过渡态,接下来做irc,但是如果把irc走出来的两端提出来再次分别进行反应物A+B和产物C+D的优化,发现这时候的优化出来的构型和前面优化出来的构型有差别,能量也有差别。这个差别里面仅仅就是A和B分子之间位置的区别,以及C和D之间分子位置的区别,单个的分子构型都是一样的。不知道我说滴明白不,最近这语言组织能力真不咋地啊,望给予指正。
4楼2014-04-20 20:01:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yurivivi

新虫 (小有名气)

引用回帖:
2楼: Originally posted by wangsihang at 2014-04-19 22:09:57
不同的基组和方法优化出来的构型会有差异的。如果仅仅是空间位置有差异或者角度有所不同应该是可以忽略的。

因为之前查看了一些帖子,比如说A和B,有人说这两个分子应该放在一起优化,因为有范德华力的存在(这个我不确定,不太懂)。所以两个分子的构型位置跟进攻位点的是有关系的,那这样的话,是不是应该A和B的相对位置在一个输入文件里应该是固定的呢?前提是同一个方法和同一个基组。后来,用irc走出来的两端试试优化构型,本来以为这样找出来反应物和产物的构型应该是要更准确些。结果两个分子之间的这个位置的差距还不小(这个都是用的同一个方法同一个基组),所以我就疑惑了,到底哪个优化出来的构型才是正确的呢?望前辈给予指正。
5楼2014-04-20 20:07:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

TINGLI_广

木虫 (小有名气)

请问,您用的是09还是03呢~谢谢
6楼2014-04-21 09:25:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yurivivi

新虫 (小有名气)

引用回帖:
6楼: Originally posted by TINGLI_广 at 2014-04-21 09:25:36
请问,您用的是09还是03呢~谢谢

09的
7楼2014-04-21 15:59:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

清新俊逸

木虫 (正式写手)

引用回帖:
9楼: Originally posted by yurivivi at 2014-04-23 10:40:45
那么请问,这样优化是合理的么?还是各自分开来优化呢?我主要想去计算反应的活化能。...

计算反应的活化能,肯定不能放在一起优化,活化能是反应物到势垒的能量差,一般还要考虑零点振动能。
至于怎么优化反应物,看你使用的是什么量子化学方法。
如果是size consistent的方法,就可以分开来优化,分别计算A和B的能量和零点能,反应物的能量就是两者之和。
如果不是size consistent的方法,比如MRCI,就要放在一起优化,并且固定A和B两者之间的距离,一般放到20个Boh日左右就可以了。因为不满足size consistent的方法,E(A+B)不等于E(A)+E(B)。
对于size consistent,翻翻教科书吧。

[ 发自手机版 http://muchong.com/3g ]
10楼2014-04-23 11:42:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yurivivi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-15 3/150 2026-09-18 02:40 by cITaGg3p5edr
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-18 02:30 by cITaGg3p5edr
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-18 01:54 by cITaGg3p5edr
[考博] 申博发邮件 +6 Lmengk 2026-09-15 7/350 2026-09-17 23:36 by Lmengk
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-17 23:16 by cuPZTDXS3VOd
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-17 23:05 by cuPZTDXS3VOd
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 9/450 2026-09-17 23:05 by cuPZTDXS3VOd
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 14:05 by ObVzqQQh1rrE
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-17 13:53 by ObVzqQQh1rrE
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:27 by ObVzqQQh1rrE
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-17 12:39 by ObVzqQQh1rrE
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-17 04:04 by k06BKNblrGRH
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-17 03:16 by Tql5LQhh5rLK
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +6 s3fFTmArrBt6 2026-09-14 6/300 2026-09-17 02:38 by Tql5LQhh5rLK
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 LwdutQ8HoqWP 2026-09-13 7/350 2026-09-17 02:02 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-16 22:26 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-16 19:15 by rRvxEj3zorCy
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:00 by i3KBRD19ERkG
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-16 03:14 by DgNGHc3h5tPl
信息提示
请填处理意见