计算化学公社

标题: 用QST2寻找过渡态时反应物和产物自旋多重度不一致该如何处理？ [打印本页]

作者
Author: unforgivable 时间: 2017-9-25 11:32
标题: 用QST2寻找过渡态时反应物和产物自旋多重度不一致该如何处理？
各位前辈好，在用QST2寻找过渡态的时候，如果反应物和产物自旋多重度不一致该如何处理？比如，HNCNH+C——>2HCN，这种反应该如何设置自选多重度呢？谢谢！

作者
Author: xujian 时间: 2017-9-25 14:21
一步反应前后整体自旋多重度怎么会变化？

作者
Author: unforgivable 时间: 2017-9-25 16:31

xujian 发表于 2017-9-25 14:21
一步反应前后整体自旋多重度怎么会变化？

这个问题我也不太清楚，只是想计算看看反应势垒有多大。

作者
Author: sobereva 时间: 2017-9-25 21:41
自旋多重度必须一致，否则整个过程不在同一个势能面上，没法找过渡态。如果说反应过程发生了势能面交叉，你应该分别考察不同自旋多重度势能面的情况，并考察势能面极小交叉点（MECP）之类。可参考Acc. Chem. Res. 2000, 33, 139-145等文章。

作者
Author: unforgivable 时间: 2017-9-26 11:08

sobereva 发表于 2017-9-25 21:41
自旋多重度必须一致，否则整个过程不在同一个势能面上，没法找过渡态。如果说反应过程发生了势能面交叉，你 ...

好的，谢谢sob老师！

作者
Author: 刘小英雄steven 时间: 2018-5-22 17:36

sobereva 发表于 2017-9-25 21:41
自旋多重度必须一致，否则整个过程不在同一个势能面上，没法找过渡态。如果说反应过程发生了势能面交叉，你 ...

老师，我的衬底是负载O原子的石墨烯，如果先吸附气体NO，而后NO与衬底的O原子反应生成NO2。这个吸附的过程就会导致衬底和最后的吸附构型的基态自旋多重度发生变化。因为气体的吸附我们也必须加入反应路径图里，这是不是就产生了势能面的交叉？那么Eads的大小（计算时衬底，气体，吸附构型均为各自基态）与后面的反应路径（自旋多重度不变）的能量变化大小就没有可比性了吗？这种情况又该如何解决呢？

作者
Author: Abner 时间: 2020-3-30 09:13

刘小英雄steven 发表于 2018-5-22 17:36
老师，我的衬底是负载O原子的石墨烯，如果先吸附气体NO，而后NO与衬底的O原子反应生成NO2。这个吸附的过 ...

我也遇到了这个问题，请问问题解决了吗？可不可以交流一下

作者
Author: 刘小英雄steven 时间: 2020-3-30 11:14
本帖最后由刘小英雄steven 于 2020-3-30 11:17 编辑

Abner 发表于 2020-3-30 09:13
我也遇到了这个问题，请问问题解决了吗？可不可以交流一下

sub与sub+gas的自旋变化在反应台阶图的吸附能中就已经体现了，后续按照吸附后的基态自旋反应就可以了。

作者
Author: Abner 时间: 2020-3-30 17:32

刘小英雄steven 发表于 2020-3-30 11:14
sub与sub+gas的自旋变化在反应台阶图的吸附能中就已经体现了，后续按照吸附后的基态自旋反应就可以了。

不好意思，新手，您说的我还是不太理解，可以详细解释一下吗？十分感谢

作者
Author: 刘小英雄steven 时间: 2020-3-31 18:00

Abner 发表于 2020-3-30 17:32
不好意思，新手，您说的我还是不太理解，可以详细解释一下吗？十分感谢

举个例子，你图上是含N煤焦zig边缘吸附了O2。无论char(N)+O2这个体系的自旋多重度和原本的char(N)的基态自旋是否一样，你的反应台阶图都已经考虑了自旋变化了，原因是计算氧气Eads时已经考虑了自旋变化。后续一定要按照char(N)+O2这一体系（即吸附后）的基态自旋，继续计算生成NOX的路径（我猜你是做氮氧化物在煤焦表面生成机理）。

作者
Author: Abner 时间: 2020-3-31 18:22

刘小英雄steven 发表于 2020-3-31 18:00
举个例子，你图上是含N煤焦zig边缘吸附了O2。无论char(N)+O2这个体系的自旋多重度和原本的char(N)的基态 ...

我懂了，首选谢谢您，我的方向和您说的差不多。我现在还有几个问题想请教一下您：
1 char（N）+O2这个体系的自旋多重度如何确定？
2 我用傻瓜式的方法，把每个自旋多重度都计算了一遍，可是发现O2吸附不上去，这个是因为什么？
3还有就是计算O2吸附时的第一步是用过渡态搜索opt=（calcfc，ts，noeigen）还是只用opt结构优化？
如您能够回答，不胜感激,谢谢，急求，再次感谢

作者
Author: 刘小英雄steven 时间: 2020-3-31 18:51

Abner 发表于 2020-3-31 18:22
我懂了，首选谢谢您，我的方向和您说的差不多。我现在还有几个问题想请教一下您：
1 char（N）+O2这个体 ...

1.基态自旋多重度大部分都是1or2，但为确保计算正确，需要算3~4个不同自旋多重度opt能量，取最低的做基态。2，多算几个是必要的，吸附不上要看具体报错原因。尤其建议调整吸附构型，比如氧气摆远一些，离N原子3埃左右试试。3，氧气吸附可以直接opt，因为你的重头戏是氧气和char(N)反应，o2吸附过程是否有ts不重要。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)