计算化学公社

标题: s1-T1和s1-T2的soc和系间窜越速率能用什么软件去计算 [打印本页]

作者
Author: zhouyulin 时间: 2024-7-24 21:14
标题: s1-T1和s1-T2的soc和系间窜越速率能用什么软件去计算
各位老师好，我想请教一下momap能计算s1-T1和s1-T2的soc和系间窜越速率吗，我看到文献里面是用orca算的，想问一下momap是否可以

作者
Author: wzkchem5 时间: 2024-7-24 21:36
可以，搭配任何能计算SOCME的量化程序即可

作者
Author: zhouyulin 时间: 2024-7-24 22:17

wzkchem5 发表于 2024-7-24 21:36
可以，搭配任何能计算SOCME的量化程序即可

谢谢老师，按照老师您的说法，那么假设我在计算s1-T2的系间窜越速率的时候，首先通过gaussian分别得到s1激发态优化好的结构和t2激发态优化好的结构，再以这两个激发态优化好的结构文件为输入条件，分别用orca得到s1-t2的soc以及用momap计算得到s1-t2的evc.cart.dat文件，最后再通过momap计算得到s1-t2的输入文件。老师，您看这样是否合理

作者
Author: wzkchem5 时间: 2024-7-24 23:56

zhouyulin 发表于 2024-7-24 15:17
谢谢老师，按照老师您的说法，那么假设我在计算s1-T2的系间窜越速率的时候，首先通过gaussian分别得到s1 ...

除了优化以外还要算频率。其他看起来没问题

作者
Author: zhouyulin 时间: 2024-7-25 16:03
本帖最后由 zhouyulin 于 2024-7-25 16:05 编辑

wzkchem5 发表于 2024-7-24 23:56
除了优化以外还要算频率。其他看起来没问题

谢谢老师，不过我有个问题想请教老师您，我这边用orca计算s1-t1的soc的时候，发现文献给出的计算结果是0.0033cm-1,但是我这边用文献提供的orca计算结果来分析发现它给的结果是0.01cm-1，我想问一下是我计算有错误吗

作者
Author: wzkchem5 时间: 2024-7-25 17:09

zhouyulin 发表于 2024-7-25 09:03
谢谢老师，不过我有个问题想请教老师您，我这边用orca计算s1-t1的soc的时候，发现文献给出的计算结果是0. ...

如果SOCME这么小的话，结构优化收敛限/SCF收敛限稍微松一点就已经足够造成这么大的区别了。但另一方面，这也意味着Herzberg-Teller效应估计占主导，所以FC结构下SOCME大一点、小一点也对ISC/磷光速率没影响。

作者
Author: zhouyulin 时间: 2024-7-25 19:06

wzkchem5 发表于 2024-7-25 17:09
如果SOCME这么小的话，结构优化收敛限/SCF收敛限稍微松一点就已经足够造成这么大的区别了。但另一方面， ...

谢谢老师，只不过这个orca输出的out文件不是我算的，而是文献的support information给的。按照soc的计算方法，它的自旋轨道耦合应该是0.01cm-1，但是文献里给的确实0.0033cm-1。文献提供的orca输出文件的soc信息与文献里面在文章实际给出的soc值不一样，所以我才想请教一下老师您，这是我算的方法有问题吗还是文献的值有问题

作者
Author: wzkchem5 时间: 2024-7-25 19:50

zhouyulin 发表于 2024-7-25 12:06
谢谢老师，只不过这个orca输出的out文件不是我算的，而是文献的support information给的。按照soc的计算 ...

那估计是文献有问题。文献可能误以为算出SOC以后应该除以3，但是没有任何道理除以3。有的情况下需要把SOC除以根号3（或者说把SOC的平方除以3再开根号），但是文章如果算S1->T1 ISC速率常数的话，也不属于应当除以根号3的情形。当然还有一个更细节的问题就是orca只给了一位有效数字，文献不应该给两位。
参见https://www.faccts.de/docs/orca/ ... isc-general-typical

作者
Author: zhouyulin 时间: 2024-7-25 20:46
本帖最后由 zhouyulin 于 2024-7-25 20:49 编辑

wzkchem5 发表于 2024-7-25 19:50
那估计是文献有问题。文献可能误以为算出SOC以后应该除以3，但是没有任何道理除以3。有的情况下需要把SOC ...

而且我也用相同的方法来计算了其他几个化合物s1-t1的soc,也是用文献给的orca文件,结果除了NMe2-DCM-OH-NBD,Br-DCM-OH和benzo2-DCM-OH跟文献结果相同以外以外，其他的都存在差别，所以我才想是不是我算错了。

作者
Author: wzkchem5 时间: 2024-7-25 21:28

zhouyulin 发表于 2024-7-25 13:46
而且我也用相同的方法来计算了其他几个化合物s1-t1的soc,也是用文献给的orca文件,结果除了NMe2-DCM-OH-NB ...

那应该是文献的问题，否则不会有的数字能对上、有的数字对不上。如果是所有化合物你算出来都和文献差一样的倍数，那还可能是哪里有一个3或者根号3没有乘/除

作者
Author: 三月星 时间: 2025-5-12 17:49
可以问一下这篇文献是什么吗？

作者
Author: JohnCase 时间: 2025-9-24 19:05

wzkchem5 发表于 2024-7-24 23:56
除了优化以外还要算频率。其他看起来没问题

听起来MOMAP算系间窜越不是那么的黑盒化？

作者
Author: wal 时间: 2025-9-24 19:24

JohnCase 发表于 2025-9-24 19:05
听起来MOMAP算系间窜越不是那么的黑盒化？

为什么这么说好像主流ESD程序都可以吃hessian+soc+de算kisc

作者
Author: JohnCase 时间: 2025-9-24 22:24

wal 发表于 2025-9-24 19:24
为什么这么说好像主流ESD程序都可以吃hessian+soc+de算kisc

我是从user-friendly的角度，我看帖子描述似乎不是只需要写一份MOMAP输入文件就能达到“端到端一键计算ISC”的目的。

作者
Author: wal 时间: 2025-9-24 22:34

JohnCase 发表于 2025-9-24 22:24
我是从user-friendly的角度，我看帖子描述似乎不是只需要写一份MOMAP输入文件就能达到“端到端一键计算IS ...

是的但其实不仅MOMAP 就算是(我觉得)最user-friendly的ORCA也要好做好几个计算才能算这个

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