计算化学公社

标题: 关于开壳层单重态过渡态的计算 [打印本页]

作者
Author: 一骑星尘 时间: 2022-4-3 19:05
标题: 关于开壳层单重态过渡态的计算
本帖最后由一骑星尘于 2022-4-3 20:27 编辑

最近在算两个自由基之间的反应，反应前有两个单电子，反应后没有单电子，想问下大家对于这种情况下，找过渡态是不是应该设置整体的自旋多重度为1，然后在基组前面加u，以及guess=mix关键词进行计算？还有没有需要注意的点呢？谢谢大家

作者
Author: 一骑星尘 时间: 2022-4-3 19:06
还是要自旋多重度为1和3的情况下都进行寻找？

作者
Author: sobereva 时间: 2022-4-3 20:02

一骑星尘发表于 2022-4-3 19:06
还是要自旋多重度为1和3的情况下都进行寻找？

有别人回复之前若需要对帖子进行修改、补充，应直接编辑原帖，不要通过回帖进行补充，这点在置顶的新社员必读贴里明确说了。

作者
Author: sobereva 时间: 2022-4-3 20:05
“基组前面加u”开玩笑
先确认这个反应到底有没有过渡态，很多情况俩自由基反应根本没过渡态
把具体反应物贴出来，免得别人猜着回答

作者
Author: 一骑星尘 时间: 2022-4-3 20:28

sobereva 发表于 2022-4-3 20:05
“基组前面加u”开玩笑
先确认这个反应到底有没有过渡态，很多情况俩自由基反应根本没过渡态
把具体反应 ...

好的社长，我已经把反应式图片放上了，怎么样才能看这个反应有没有过渡态呢？如果想找的话应该怎么找呢？谢谢

作者
Author: wzkchem5 时间: 2022-4-3 22:28

一骑星尘发表于 2022-4-3 13:28
好的社长，我已经把反应式图片放上了，怎么样才能看这个反应有没有过渡态呢？如果想找的话应该怎么找呢？ ...

先查文献，了解碳自由基和Ag(II)-OR键反应一般是什么机理，是形成C-Ag键再还原消除，还是按SH2机理碳自由基从Ag原子的背面进攻氧原子。然后按文献结论来摆放反应物，作为找过渡态的初猜。
此外注意，实验文章画的F3CO-Ag(II)-F这样的物种，其实未必是monomeric的，很可能是多聚物。应该查文献了解一下二价银的氟化物在溶液里一般是什么形式（必须找计算文章，或者测了聚集态的实验文章。那些只证明了Ag和F元素比例是1:2就说是AgF2的文章不要信），然后按文献结论构建反应物结构。产物Ag(I)F同理，几乎不可能是单分子的，而是(AgF)n多聚体团簇，这个n应该是多少、Ag和F的构型如何，乃至团簇表面有没有配溶剂、配了多少个溶剂分子，都需要查文献确定。查不到文献就用molclus（http://bbs.keinsci.com/thread-577-1-1.html）自己做构型搜索计算。

作者
Author: 一骑星尘 时间: 2022-4-4 09:50

wzkchem5 发表于 2022-4-3 22:28
先查文献，了解碳自由基和Ag(II)-OR键反应一般是什么机理，是形成C-Ag键再还原消除，还是按SH2机理碳自由 ...

好的，如果是先形成C-Ag键再进行还原消除，似乎也存在着一个整体上两个单电子到无单电子的一个过程，对于这种过程的过渡态应该如何计算呢？自旋多重度应该如何设置？用什么关键词呢？

作者
Author: fatpig 时间: 2022-4-4 12:29
本帖最后由 fatpig 于 2022-4-4 12:30 编辑

如果是我的话：
a. 先做triplet的单点，stable=opt
b. 用+u读triplet的guess，做singlet (open-shelled singlet, OSS)的单点, stable=opt
c. 确认OSS在stable之后存在且能量确实比triplet低(当然高一点也无妨，可能优化后变低)
d. OSS单点读上一步的guess，算freq，确认振动虚频是否是需要的
e. ts, readfc做opt

作者
Author: wzkchem5 时间: 2022-4-4 15:14

一骑星尘发表于 2022-4-4 02:50
好的，如果是先形成C-Ag键再进行还原消除，似乎也存在着一个整体上两个单电子到无单电子的一个过程，对于 ...

自旋多重度设为1，但是初猜要用broken symmetry初猜

作者
Author: fwz3888 时间: 2022-4-8 16:26
一般来说，先一个stable=opt任务拿到BS的波函数，后续的任务读这个波函数在做
不过建议check这个反应是Sh2还是自由基进攻Ag中心的过程
看上去像是nku某老师的体系hhh

作者
Author: 火星哥 时间: 2022-5-31 20:42
可以用QST2，虽然很可耻。但我一试就出来了

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