计算化学公社

 找回密码 Forget password
 注册 Register
Views: 1891|回复 Reply: 10
打印 Print 上一主题 Last thread 下一主题 Next thread

[辅助/分析程序] s1-T1和s1-T2的soc和系间窜越速率能用什么软件去计算

[复制链接 Copy URL]

68

帖子

0

威望

458

eV
积分
526

Level 4 (黑子)

各位老师好,我想请教一下momap能计算s1-T1和s1-T2的soc和系间窜越速率吗,我看到文献里面是用orca算的,想问一下momap是否可以

1万

帖子

0

威望

9001

eV
积分
20757

Level 6 (一方通行)

2#
发表于 Post on 2024-7-24 21:36:56 | 只看该作者 Only view this author
可以,搭配任何能计算SOCME的量化程序即可
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

458

eV
积分
526

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2024-7-24 22:17:16 | 只看该作者 Only view this author
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的输入文件。老师,您看这样是否合理

1721830206919.jpg (44.97 KB, 下载次数 Times of downloads: 44)

momap计算kisc的输入脚本

momap计算kisc的输入脚本

1万

帖子

0

威望

9001

eV
积分
20757

Level 6 (一方通行)

4#
发表于 Post on 2024-7-24 23:56:56 | 只看该作者 Only view this author
zhouyulin 发表于 2024-7-24 15:17
谢谢老师,按照老师您的说法,那么假设我在计算s1-T2的系间窜越速率的时候,首先通过gaussian分别得到s1 ...

除了优化以外还要算频率。其他看起来没问题
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

458

eV
积分
526

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2024-7-25 16:03:56 | 只看该作者 Only view this author
本帖最后由 zhouyulin 于 2024-7-25 16:05 编辑
wzkchem5 发表于 2024-7-24 23:56
除了优化以外还要算频率。其他看起来没问题

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

1721894595673.jpg (15.56 KB, 下载次数 Times of downloads: 39)

文献给出的计算结果

文献给出的计算结果

1721894624487.jpg (21.75 KB, 下载次数 Times of downloads: 37)

文献提供的orca输出文件的计算结果

文献提供的orca输出文件的计算结果

DCM-OH.out

259.77 KB, 下载次数 Times of downloads: 3

文献提供的orca计算文件

1万

帖子

0

威望

9001

eV
积分
20757

Level 6 (一方通行)

6#
发表于 Post on 2024-7-25 17:09:38 | 只看该作者 Only view this author
zhouyulin 发表于 2024-7-25 09:03
谢谢老师,不过我有个问题想请教老师您,我这边用orca计算s1-t1的soc的时候,发现文献给出的计算结果是0. ...

如果SOCME这么小的话,结构优化收敛限/SCF收敛限稍微松一点就已经足够造成这么大的区别了。但另一方面,这也意味着Herzberg-Teller效应估计占主导,所以FC结构下SOCME大一点、小一点也对ISC/磷光速率没影响。
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

458

eV
积分
526

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2024-7-25 19:06:29 | 只看该作者 Only view this author
wzkchem5 发表于 2024-7-25 17:09
如果SOCME这么小的话,结构优化收敛限/SCF收敛限稍微松一点就已经足够造成这么大的区别了。但另一方面, ...

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

1万

帖子

0

威望

9001

eV
积分
20757

Level 6 (一方通行)

8#
发表于 Post on 2024-7-25 19:50:23 | 只看该作者 Only view this author
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
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

68

帖子

0

威望

458

eV
积分
526

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2024-7-25 20:46:58 | 只看该作者 Only view this author
本帖最后由 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跟文献结果相同以外以外,其他的都存在差别,所以我才想是不是我算错了。

1721911339196.jpg (11.66 KB, 下载次数 Times of downloads: 32)

数据对比

数据对比

1721911371337.jpg (22.36 KB, 下载次数 Times of downloads: 36)

1721911371337.jpg

1721911426280.jpg (21.72 KB, 下载次数 Times of downloads: 38)

1721911426280.jpg

Br-DCM-OH.out

255.31 KB, 下载次数 Times of downloads: 1

benzo2-DCM-OH.out

265.22 KB, 下载次数 Times of downloads: 0

NMe2-DCM-O-NBD.out

285.94 KB, 下载次数 Times of downloads: 0

1万

帖子

0

威望

9001

eV
积分
20757

Level 6 (一方通行)

10#
发表于 Post on 2024-7-25 21:28:47 | 只看该作者 Only view this author
zhouyulin 发表于 2024-7-25 13:46
而且我也用相同的方法来计算了其他几个化合物s1-t1的soc,也是用文献给的orca文件,结果除了NMe2-DCM-OH-NB ...

那应该是文献的问题,否则不会有的数字能对上、有的数字对不上。如果是所有化合物你算出来都和文献差一样的倍数,那还可能是哪里有一个3或者根号3没有乘/除
Zikuan Wang
山东大学光学高等研究中心 研究员
BDF(https://bdf-manual.readthedocs.io/zh_CN/latest/Introduction.html)、ORCA(https://orcaforum.kofo.mpg.de/index.php)开发团队成员
Google Scholar: https://scholar.google.com/citations?hl=zh-CN&user=XW6C6eQAAAAJ&view_op=list_works&sortby=pubdate
ORCID: https://orcid.org/0000-0002-4540-8734
主页:http://www.qitcs.qd.sdu.edu.cn/info/1034/1702.htm
本团队长期招收研究生,有意者可私信联系

2

帖子

0

威望

227

eV
积分
229

Level 3 能力者

11#
发表于 Post on 2025-5-12 17:49:58 | 只看该作者 Only view this author
可以问一下这篇文献是什么吗?

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2025-8-16 15:15 , Processed in 0.159429 second(s), 23 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list