计算化学公社

标题: 势能面圆锥交叉点如何计算 [打印本页]

作者
Author: Kikyou 时间: 2020-9-23 10:50
标题: 势能面圆锥交叉点如何计算
请教一下各位老师，我在看一篇文章的时候看到了这样的机理路径，作者给了一个圆锥交叉点，并在文章中这样描述：Our results support the notion that hydrogen transfer to form the oquinodimethane must occur in the triplet state; this is followed by intersystem crossing at the conical intersection point of E− Z isomerization to form the singlet (E)- or (Z)-dienol.
我想知道这个交叉点是怎么算出来的啊，是跟我们平时讲的MECP有什么区别嘛？
谢谢各位老师！

作者
Author: 喵星大佬 时间: 2020-9-23 11:01
opt=conical

作者
Author: zjxitcc 时间: 2020-9-23 11:02

喵星大佬发表于 2020-9-23 11:01
opt=conical

思想精髓确实是这个，只不过直接写这个并不能算

。。。

作者
Author: 喵星大佬 时间: 2020-9-23 11:05

zjxitcc 发表于 2020-9-23 11:02
思想精髓确实是这个，只不过直接写这个并不能算。。。

你要这么说也有道理。。。。不过找圆锥交叉点是个很复杂的事情。。。。而且不同自选多重度的应该叫MECP。。。。用sobMECP找吧，论坛搜一下帖子就有了，或者把嘴哥那篇教学发来？

作者
Author: sobereva 时间: 2020-9-23 11:15
单重态-单重态之间的极小交叉点习俗上叫圆锥交叉点，不同自旋态之间的习俗上叫MECP。当前这个情况叫做MECP更符合习俗。
sobMECP可以很容易计算
使用sobMECP程序结合Gaussian程序搜索极小能量交叉点
http://sobereva.com/286（http://bbs.keinsci.com/thread-865-1-1.html）

作者
Author: biogon 时间: 2020-9-23 11:23
https://mp.weixin.qq.com/s?__biz ... 4f1db962a4a23b29#rd
看这个就知道了

作者
Author: wzkchem5 时间: 2020-9-23 11:33
应该是写错了，是MECP。
不同自旋态的交叉只有在考虑旋轨耦合的情况下才是锥形交叉（交叉空间N-2维），否则就是普通的势能面交叉（交叉空间N-1维），只能说MECP而不能说锥形交叉点。
另外即使是MECP，这个图也有问题，S0态要么和T1相交2次，要么不相交，几乎不可能刚好相切；即使刚好相切，切点也几乎不可能刚好是过渡态。所以作者应该是少算了一个MECP，或者把一个本来只涉及S0态的反应途径连到MECP上去了。

作者
Author: Kikyou 时间: 2020-9-23 12:13

喵星大佬发表于 2020-9-23 11:05
你要这么说也有道理。。。。不过找圆锥交叉点是个很复杂的事情。。。。而且不同自选多重度的应该叫MECP。 ...

大佬求一下你们说的这个教学教程

作者
Author: Kikyou 时间: 2020-9-23 12:14

wzkchem5 发表于 2020-9-23 11:33
应该是写错了，是MECP。
不同自旋态的交叉只有在考虑旋轨耦合的情况下才是锥形交叉（交叉空间N-2维），否 ...

嗯嗯，谢谢老师回复。我跟文章作者联系下问问她吧~

作者
Author: 喵星大佬 时间: 2020-9-23 12:36

Kikyou 发表于 2020-9-23 12:13
大佬求一下你们说的这个教学教程

点zjx下面的那个推广链接，里面有不少MECP相关的

作者
Author: Kikyou 时间: 2020-9-23 14:05

喵星大佬发表于 2020-9-23 12:36
点zjx下面的那个推广链接，里面有不少MECP相关的

嗯嗯看到了，谢谢！！

作者
Author: sobereva 时间: 2020-9-24 10:57
再特别强调一下，你若要算文中的MECP（文中说圆锥交叉那是缪称，前面已经有人提了），应当用的是sobMECP算S-T之间的MECP，在DFT下即可实现，极为容易。而若用opt=conical来找MECP，你必须用CASSCF，这要复杂麻烦得多得多。

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