计算化学公社

标题: 求助ORCA手册 蒽 分子 旋轨耦合计算问题 [打印本页]

作者
Author:
是aweia    时间: 2022-3-29 13:52
标题: 求助ORCA手册 蒽 分子 旋轨耦合计算问题
本帖最后由 是aweia 于 2022-3-29 17:08 编辑

根据ORCA手册,在计算蒽分子的ISC速率时,需要首先计算S1与对应三重态各子态间的旋轨耦合,以下为手册原文
As a simple example, one could compute the excited singlet and ground triplet geometries and Hessians foranthracene using TD-DFT, then compute the SOC matrix elements with a given triplet spin-sublevel by the same method (see the details below), maybe even using CASSCF, MRCI, STEOM-CCSD or some other theory level and finally obtain the ISC rates using an input such as:


!ESD(ISC) NOITER

%ESD ISCISHESSIAN "ANT_S1.hess"

ISCFSHESSIAN "ANT_T1.hess"

DELE 11548

SOCME 0.0, 2.33e-5

END

* XYZFILE 0 1 ANT_T1.xyz


另外手册中说明,ESD模块内SOCME应使用的是原子单位,:

The relevant matrix elements can be calculated from any method available in ORCA and inputed as SOCME Re,Im under %ESD

where Re and Im are its real and imaginary parts (in in atomic units!)

也就是说将上述手册ESD模块中SOCME的值2.33e-5为Eh单位 若转变为cm-1单位就是约为5.11 cm-1

但我实际算的蒽分子的S1到T1\T2d的旋轨耦合几乎为零,差别比较大

请问是我理解的有问题吗,还是我计算的有问题 以下是我的输入文件,请各位批评指正

! B3LYP/G TZVP miniprint tightSCF
%pal
   nprocs 64
end

%tddft
nroots 20
dosoc true
tda false
printlevel 10
end
* xyzfile 0 1 anthracene.xyz
同时我把输出文件上传








作者
Author:
sobereva    时间: 2022-3-29 14:37
如置顶的新社员必读贴和论坛的版头所示,上传超过500KB的文本型文件上传前必须先压缩再上传,以节约论坛空间、节约下载耗时。
作者
Author:
是aweia    时间: 2022-3-29 14:46
sobereva 发表于 2022-3-29 14:37
如置顶的新社员必读贴和论坛的版头所示,上传超过500KB的文本型文件上传前必须先压缩再上传,以节约论坛空 ...

好的  我已经将文件压缩重新上传 今后我会注意类似问题的
作者
Author:
是aweia    时间: 2022-3-29 15:49
请问大家有没有ORCA的软件交流群 因为课题组只有我一个人在用  时常遇到问题的时候感觉找不到组织   希望能找个群和大家一起交流
作者
Author:
冰释之川    时间: 2022-3-29 16:01
本帖最后由 冰释之川 于 2022-3-29 16:05 编辑

先重现帖子里的实例吧:

使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱
http://sobereva.com/462http://bbs.keinsci.com/thread-12114-1-1.html
作者
Author:
是aweia    时间: 2022-3-29 16:15
冰释之川 发表于 2022-3-29 16:01
先重现帖子里的实例吧:

使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱

感谢回复  我在计算之前已经读过且重现过sob这篇帖子的第一个计算实例:甲醛的计算,我的计算也是按照那个为参考的,结果出来却与ORCA手册相去甚远 ,不知道哪里出了问题
作者
Author:
zjxitcc    时间: 2022-3-29 16:42
是aweia 发表于 2022-3-29 15:49
请问大家有没有ORCA的软件交流群 因为课题组只有我一个人在用  时常遇到问题的时候感觉找不到组织   希望能 ...

ORCA有官方论坛https://orcaforum.kofo.mpg.de,每天都有不少人提问和回复,算活跃的(不过最近不知为啥似乎上不去了)
作者
Author:
是aweia    时间: 2022-3-29 17:06
zjxitcc 发表于 2022-3-29 16:42
ORCA有官方论坛https://orcaforum.kofo.mpg.de,每天都有不少人提问和回复,算活跃的(不过最近不知为啥 ...

好的 谢谢您
作者
Author:
冰释之川    时间: 2022-3-29 17:19
zjxitcc 发表于 2022-3-29 16:42
ORCA有官方论坛https://orcaforum.kofo.mpg.de,每天都有不少人提问和回复,算活跃的(不过最近不知为啥 ...

刚测试了一下,能上去啊
作者
Author:
Novice    时间: 2022-3-30 08:38
冰释之川 发表于 2022-3-29 17:19
刚测试了一下,能上去啊

那你的网和大多数都不太一样
作者
Author:
tiandikuoyuan    时间: 2022-3-30 08:57
手册里面说了这个方法计算蒽的ISC不合适,需要考虑HT effect;你计算SOC使用的是基态构型吗,可能需要使用T1的,你这结构S1和T2的能级差和手册上的结果对不上
作者
Author:
是aweia    时间: 2022-3-30 23:38
本帖最后由 是aweia 于 2022-3-30 23:45 编辑
tiandikuoyuan 发表于 2022-3-30 08:57
手册里面说了这个方法计算蒽的ISC不合适,需要考虑HT effect;你计算SOC使用的是基态构型吗,可能需要使用T ...

您好  非常感谢您的回复,非常乐意一起讨论这个问题
1、确实手册里提到蒽不适合这个从耦合矩阵元开始的方法,因为其中只计入了FC的部分,还需要考虑HT效应 ,对此我有个疑问,考虑HT效应后的方法和只考虑FC的方法得出的结果是需要相加还是只把包含HT的计算的结果当作最终结果,也就是说其实我疑问的地方在于考虑HT效应的方法里是否同时包含了FC效应2、我计算soc矩阵元时使用的是基态构型,按您所说我使用t1态的结构重算了一下,结果依然不对,希望能得到您的指导,谢谢。
3、我用我的结构分别在s1和t2态下做结构优化得到的能量符合手册描述,不知道您是怎么做的,我上传了我的s1和t2的out文件可以作为参考比较。

作者
Author:
wzkchem5    时间: 2022-3-31 00:02
是aweia 发表于 2022-3-30 16:38
您好  非常感谢您的回复,非常乐意一起讨论这个问题
1、确实手册里提到蒽不适合这个从耦合矩阵元开始的 ...

1. 是包括的,所以不用另做FC计算
我有模糊的印象这个体系手册报的结果好像和手册给的输入文件不完全一致,不过找不到原帖了,建议还是在orca论坛上发帖问一下。ESD模块是Bernardo de Souza负责的,他应该会在orca论坛上回你帖。
作者
Author:
tiandikuoyuan    时间: 2022-3-31 09:17
是aweia 发表于 2022-3-30 23:38
您好  非常感谢您的回复,非常乐意一起讨论这个问题
1、确实手册里提到蒽不适合这个从耦合矩阵元开始的 ...

我之前是用Gaussian优化激发态结构并计算频率,然后再转成orca识别的hess文件;orca计算激发态频率有点难度。后来就没有计算ISC速率了,我觉得最后只要速率数值和文献上一致就行了。
作者
Author:
是aweia    时间: 2022-3-31 10:19
tiandikuoyuan 发表于 2022-3-31 09:17
我之前是用Gaussian优化激发态结构并计算频率,然后再转成orca识别的hess文件;orca计算激发态频率有点难 ...

Gaussian 我有使用过 不过对于您所说的“Gaussian优化激发态结构并计算频率,然后再转成orca识别的hess文件”
具体是怎么转换呢 ,因为我现在做的计算体系较大 现在算频率非常慢,我觉得您的办法会很有帮助,非常期待得到您的解答,谢谢
作者
Author:
是aweia    时间: 2022-3-31 10:21
wzkchem5 发表于 2022-3-31 00:02
1. 是包括的,所以不用另做FC计算
我有模糊的印象这个体系手册报的结果好像和手册给的输入文件不完全一 ...

非常感谢您的指导
作者
Author:
tiandikuoyuan    时间: 2022-4-2 21:38
是aweia 发表于 2022-3-31 10:19
Gaussian 我有使用过 不过对于您所说的“Gaussian优化激发态结构并计算频率,然后再转成orca识别的hess文 ...

我编译出来的程序文件比较大,源代码见 (, 下载次数 Times of downloads: 20) ;我不保证结果,之前用的时候是没有问题的





欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3