|
|
您好 我最近仍然在计算ISC速率 感谢您的指导现在可以进行计算了 但仍有疑惑想向您请教 多次打扰 请您见谅
按照手册计算ISC的说明 按照是否考虑Herzberg-Teller effect分为两种方法
1 计算初末态间的SOCME和两个态各自的.hess文件和末态的结构，一起作为计算ISC速率的输入信息。
2 分别计算初末态的.hess文件以及末态的结构，并将其一起作为计算ISC速率的输入信息。
我自己按照手册说明算了一下，发现几个问题想请教您（手册计算ISC所用的分子是anthracene，结构如图）
- 1 我计算的SOCME与手册的不一致：手册展示了方法一计算ISC那一步的输入文件如下：
ESD(ISC) NOITER
%ESD
ISCISHESSIAN "ANT_S1.hess"
ISCFSHESSIAN "ANT_T1.hess"
DELE 11548
SOCME 0.0, 2.33e-5 （原子单位，虚部转换为cm-1后应该为5.11）
END
* XYZFILE 0 1 ANT_T1.xyz
但我计算出的SOCME 无论实部还是虚部都是零，我的SOC计算输入主要部分如下
%tddft
nroots 5
dosoc true
tda false
printlevel 3
end
* xyzfile 0 1 opt.xyz
输出的矩阵为
--------------------------------------------------------------------------------
CALCULATED SOCME BETWEEN TRIPLETS AND SINGLETS
--------------------------------------------------------------------------------
Root <T|HSO|S> (Re, Im) cm-1
T S MS= 0 -1 +1
--------------------------------------------------------------------------------
1 0 ( 0.00 , -0.00) ( -0.00 , 0.00) ( -0.00 , -0.00)
1 1 ( 0.00 , -0.00) ( 0.00 , 0.00) ( 0.00 , -0.00)
1 2 ( 0.00 , 0.04) ( -0.00 , 0.00) ( -0.00 , -0.00)
1 3 ( 0.00 , -0.00) ( -0.00 , -0.00) ( -0.00 , 0.00)
1 4 ( 0.00 , 0.00) ( -0.00 , -0.00) ( -0.00 , 0.00)
1 5 ( 0.00 , 0.00) ( -0.00 , -0.00) ( -0.00 , 0.00)
想跟您请教这样做有没有什么不合适的地方。
- 2 按照方法二 我计算得到了S1和T1态的.hess文件以及T1态的结构，计算的方法如下：
| S1态： | T1态： |
!B3LYP DEF2-TZVP(-F) Opt freq %tddft
nroots 5
IROOT 1
end
...
| !B3LYP DEF2-TZVP(-F) Opt freq
%TDDFT
NROOTS 5
IROOTMULT TRIPLET
END
...
|
但是用计算出来的.hess文件和结构，并按照手册给出的输入文件进行计算得出的结果和手册大不一样（手册结果是十的八次方，我的结果还不到1),且我从S1和T1的计算结果中得到的DELE参数为11478cm-1，手册里为11548cm-1(我是按照两个态的输出文件中的电子能Electronic energy做差值得到的，不知道是否正确），所以只能是我的S1和T1的计算有问题，但我找不到问题在哪里，想请您指教一下。 - 3 按照方法二进行ISC的计算时，作业能最终结束并输出ORCA TERMINATED NORMALLY的信息并打印出isc速率（The calculated ISC rate constant is 2.162519e-01 s-1 with 0.00 from FC and 100.00 from HT），但在计算时，输出文件输出有报错信息且其始终保持在文件末尾的位置，如图，不知道这个会不会影响计算结果?
如上，问题有点多，多有打扰，真的非常抱歉。同时很期待您的指导，万分感谢。
|
|
发表于 Post on 2022-2-10 00:33:13
收藏 Add to favorites
楼主 Author
