liyuanhe211 发表于 2017-3-10 11:31 老师你好,我也有个类似的三维立体的扫描图,是通过对一个分子同时扫描的两个二面角得来的,但是我不理解这个图有什么意义,然后第二个问题是什么时候采用两个柔性扫描同时进行的方法。 |
珊珊来迟 发表于 2017-3-10 11:19 缩短扫描就把步长改为负值就行。 不过有个情况。就是irc达不到的态 opt就能够达到,因为irc和opt用的优化算法不一样。opt好像能够跨过一些小的能垒。而IRC就不行。那么对于这种因为irc算法产生的‘准中间体’,能不能直接忽略?@liyuanhe |
本帖最后由 liyuanhe211 于 2017-3-10 11:32 编辑 珊珊来迟 发表于 2017-3-10 10:42 (为了方便理解,下面的描述很多地方不准确)请想象一个二维势能面(我没有好的图),一维是O-O键键长 x,另一维是“H的位置” y,将反应式(1)的底物(是一个极小点)设为原点(0,0),产物设为第一象限上的某个极小点(x,y),你需研究的是这两个极小点之间的路径。 如果存在协同机理,且是最优的,那么直接的(0,0)-->(x,y)的过程是最优的,此时如果保证优化总是能得到最小点,那么在柔性扫描 x 这个变量时,就应该观察到 y 这个变量也同步变化,并基本同时到达终点 (x,y),在这个路径中有一个极大点(一阶鞍点),是协同过渡态; 如果分布机理是较优的,那么反应路径有 (0,0)-->(x,0)-->(x,y) 或 (0,0)-->(0,y)-->(x,y)两种,此时柔性扫描 x 变量(并保证得到最小点),就应该观察到 y 变量基本不变化,此时就可以推测其分布机理是较优的,并且在柔性扫描图中能找到 O-O 键形成/断裂步骤的过渡态,然后观察扫描过程中 y 是接近0还是终点可以判断分布时哪一步在前,再找H转移这步的过渡态。(当然该反应中我猜如果分步、肯定是H转移在前、O-O断裂在后) 做两边扫描只是为了满足得到最小点这个要求,相当于优化时用两种初猜优化比较能量,取最优的。 |
| 参与人数Participants 1 | eV +1 | 收起 理由Reason |
|---|---|---|
|
| + 1 |
ggdh 发表于 2017-3-10 11:10 另外这个中间体是o-o键断过之后的结构吧,那就是个双自由基结构了? |
ggdh 发表于 2017-3-10 11:10 恩恩,我会试一下这个方法。扫描的时候我只会由初始结构向键长拉伸的方向扫描,请问怎么从该结构向缩短的方向扫描啊? |
| 把你向反应物irc终点取出来,做一个optimization。看看 oo键连起来没有。如果没有,那么你得到了一个中间体。然后从这个中间体出发,缩短oo键扫描。看看能不能找到oo键先断裂的过渡态。 |
liyuanhe211 发表于 2017-3-10 09:47 谢谢老师的解答!这种方法学生还是第一次学习到!又能get到一个新技能了~不过老师对于您说的方法我还有几点不太明白的地方 1.将两次扫描的结果合并是要借助什么技术手段把两次扫描的结果做成一个曲线来观察还是直接分开看两个scan结果就可以?最高点优化是指取出最高点结构为初猜做opt=ts吧?最高点是指哪个结果曲线的最高点呢? 2.“取能量较低的点画成曲线” 老师可否再解释详细一下,是指两次scan合并后的结果中取出能量较低的点即可么? 3.根据上述结果怎么就能判断出是协同还是分步?还是没太懂。。 4.老师最后说的是指O-O与N-C-H同时扫描,看有无能量最高点么 |
|
本帖最后由 liyuanhe211 于 2017-3-10 09:49 编辑 IRC和振动动画矛盾时以IRC为准(只要IRC做的正确)。 可以使初始结构中,H连于C上,然后对O-O键键长做柔性扫描;然后使初始结构,H连于N上,对O-O扫描,两次扫描所取的O-O键键长值相同;将两次扫描结果合并,取能量较低的点画成曲线,辅以最高点的优化,很可能能看出是协同的还是分布的;分两次扫描去最小是为了防止柔性扫描带到局域极小去。 我没有查过、但之前确实没见过这种协同的自由基过程,要是想完全确定是否确实存在协同过渡态(除了尝试各种初猜做过渡态优化),可以做O-O键和N-C-H键角的二维柔性扫描,但耗时较高。 |
| 参与人数Participants 1 | eV +1 | 收起 理由Reason |
|---|---|---|
|
| + 1 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-22 00:09 , Processed in 0.175115 second(s), 26 queries , Gzip On.