ESD模块计算ISC速率手册算例计算出错求指教

是aweia · 发表于 Post on 2022-2-10 00:33:13

我在重复ORCA5.0.2手册关于ISC速率计算的内容（8.15.4.2节）计算过程：首先需要计算获得anthracene分子的S1态和T1态的.hess文件
我的输入文件如下：
S1:
!B3LYP DEF2-TZVP(-F)  Opt freq
%tddft
nroots 5
IROOT 1
end
* xyzfile 0 1 anthracene.xyz

T1：
! B3LYP DEF2-TZVP(-F)  Opt freq
%TDDFT NROOTS 5
IROOTMULT TRIPLET
END
* xyzfile 0 1 anthracene.xyz

计算完成后将两个.hess文件及T1计算生成的.xyz用在下一步计算中
然后按照手册的example计算ISC速率，输入文件如下：
!B3LYP DEF2-TZVP(-F) TIGHTSCF ESD(ISC)
%TDDFT
NROOTS 5
SROOT 1
TROOT 1
TROOTSSL 0
DOSOC TRUE
END
%ESD
ISCISHESS "ANT_S1.hess"
ISCFSHESS "ANT_T1.hess"
USEJ TRUE
DOHT TRUE
TEMP 77
DELE 11548
END
* XYZFILE 0 1 ANT_T1.xyz
$NEW_JOB
!B3LYP DEF2-TZVP(-F) TIGHTSCF ESD(ISC)
%TDDFT
NROOTS 5
SROOT 1
TROOT 1
TROOTSSL -1
DOSOC TRUE
END
%ESD
ISCISHESS "ANT_S1.hess"
ISCFSHESS "ANT_T1.hess"
USEJ TRUE
DOHT TRUE
TEMP 77
DELE 11548
END
* XYZFILE 0 1 ANT_T1.xyz
...（这三个点是手册中就有的，我计算时没有写进去，是否是我漏掉了什么输入参数）
计算时报错如下：
------------------------------------------------------------------------------
                        ORCA EXCITED STATE DYNAMICS
------------------------------------------------------------------------------

Requested calculation:                         ...ISC
Coordinate system:                            ...Deloc. Internal (BAKER)
Lineshape function:                            ...Gaussian
Excited state PES:                            ...adiabatic hessian
Apply Duschinky rotation:                      ...yes
Use operator derivatives:                      ...yes
Convert derivatives between states:          ...no
Use same frequencies:                         ...no
Multiroot calculation:                         ...no
Single mode calculation:                      ...no
Selected mode calculation:                   ...no

Initial state hessian:                         ANT_S1.hess
Reading initial state hessian:                ...done
Final state hessian:                         ANT_T1.hess
Reading final state hessian:                   ...done
[file orca_tools/qcmat1.cpp, line 103]: Stored DataType 808463920 does not match expected 8 in OrcaObject_ReadHeader!

ORCA finished by error termination in ORCA_ESD
Calling Command: /public/home/bin/ORCA/orca_5_0_2_linux_x86-64_shared_openmpi411/orca_esd way.esdinp.tmp  way /public/home/anthracene/way2
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run

xu123123 · 发表于 Post on 2025-10-31 16:17:18

wzkchem5 发表于 2025-10-31 10:23
可以，5.0.3及以后的版本的ISC计算应该都没有bug

好的谢谢老师

wzkchem5 · 发表于 Post on 2025-10-31 10:23:13

xu123123 发表于 2025-10-30 22:28
老师您好，请问orca6.0.1的版本能计算 ISC速率吗

可以，5.0.3及以后的版本的ISC计算应该都没有bug

xu123123 · 发表于 Post on 2025-10-30 22:28:07

wzkchem5 发表于 2022-2-10 01:37
ORCA 5.0.2计算ISC有bug，参见https://orcaforum.kofo.mpg.de/viewtopic.php?f=54&t=8493
解决方法是退回 ...

老师您好，请问orca6.0.1的版本能计算 ISC速率吗

qian. · 发表于 Post on 2024-11-14 10:23:59

wzkchem5 发表于 2024-11-14 10:17
ORCA手册里给了换算系数，https://www.faccts.de/docs/orca/6.0/manual/contents/general.html#units-and ...

好的，谢谢

wzkchem5 · 发表于 Post on 2024-11-14 10:17:19

qian. 发表于 2024-11-14 02:49
你好，我想请教一下原子单位与波数的换算，我算了半天都不对

ORCA手册里给了换算系数，https://www.faccts.de/docs/orca/ ... -conversion-factors

qian. · 发表于 Post on 2024-11-14 09:49:59

是aweia 发表于 2022-2-17 01:44
您好我最近仍然在计算ISC速率感谢您的指导现在可以进行计算了但仍有疑惑想向您请教多次打扰请您见谅 ...

你好，我想请教一下原子单位与波数的换算，我算了半天都不对

wzkchem5 · 发表于 Post on 2022-8-10 17:35:23

waker 发表于 2022-8-10 09:49
那老师用TDDFT计算T2的频率的关键词还是吗？
%TDDFT NROOTS 5
IROOTMULT TRIPLET

原则上有两种做法，一种是以S0为参考态算T2，一种是以T1为参考态算T2。
如果是前者的话，在你输入的基础上把IROOT设为2，同时注意自旋多重度是1而不是3。
如果是后者的话，自旋多重度要设为3，IROOTMULT TRIPLET要删掉（因为当基态是UKS时，TDDFT不能指定自旋多重度），IROOT设为多少需要根据具体体系确定，因为当以T1为参考态时，第一激发态未必是T2，也可能是五重态。但是这个计算方法有自旋污染，所以除非特殊情况（比如S0波函数不稳定）否则不要用

waker · 发表于 Post on 2022-8-10 16:49:29

wzkchem5 发表于 2022-8-10 14:51
S1和T2的ISC只可能涉及S1的hess文件和T2的hess文件（严格的处理两者都需要，VG等近似方法只需要两者之一 ...

那老师用TDDFT计算T2的频率的关键词还是吗？
%TDDFT NROOTS 5
IROOTMULT TRIPLET

wzkchem5 · 发表于 Post on 2022-8-10 14:51:05

waker 发表于 2022-8-10 04:03
不好意思，我想问一下关于S1与T2的ISC计算的hess 文件和 T1的hess 文件是同一个吗？如果不是如何求得？谢 ...

S1和T2的ISC只可能涉及S1的hess文件和T2的hess文件（严格的处理两者都需要，VG等近似方法只需要两者之一），不涉及T1的hess文件。需要单独用TDDFT计算S1和/或T2的频率

waker · 发表于 Post on 2022-8-10 11:03:03

是aweia 发表于 2022-2-18 21:50
非常感谢您的指导，我补充了S1与T2的ISC计算后得出的结果和手册接近了，非常感谢您。

不好意思，我想问一下关于S1与T2的ISC计算的hess 文件和 T1的hess 文件是同一个吗？如果不是如何求得？谢谢

是aweia · 发表于 Post on 2022-2-18 21:50:39

wzkchem5 发表于 2022-2-17 18:35
SOC计算不要和opt freq写在一起，不然你优化的既不是S1态结构也不是T1态结构，而是旋量态的结构。应该先 ...

非常感谢您的指导，我补充了S1与T2的ISC计算后得出的结果和手册接近了，非常感谢您。

wzkchem5 · 发表于 Post on 2022-2-17 18:35:54

是aweia 发表于 2022-2-17 04:31
谢谢我将所有的输入输出文件一起打包了上传给您非常感谢

SOC计算不要和opt freq写在一起，不然你优化的既不是S1态结构也不是T1态结构，而是旋量态的结构。应该先做opt freq，再用一个单独的输入文件在优化好的结构下做SOC单点计算。
剩下的问题我晚些时候再看一下

是aweia · 发表于 Post on 2022-2-17 17:09:18

wzkchem5 发表于 2022-2-17 01:58
这几个计算的完整输入文件（包括分子坐标）可否给一下？我看看能不能重复出这些问题。

您好关于我的ISC速率结果和手册对不上的问题，我又看了一下手册，应该是我当时看的时候脑子糊涂理解错了，手册的ISC速率也包含了S1与T2态间的系间穿越，而我没算这部分，以及S1与T1态的第三个子态（TROOTSSL 1)之间的部分。所以这部分应该是我的问题，就不麻烦您了，我这就去补充这部分的计算，就请教您一下这个SOCME的计算的问题就行，给您带来的困扰真的非常抱歉。

是aweia · 发表于 Post on 2022-2-17 11:31:04

wzkchem5 发表于 2022-2-17 01:58
这几个计算的完整输入文件（包括分子坐标）可否给一下？我看看能不能重复出这些问题。

谢谢我将所有的输入输出文件一起打包了上传给您非常感谢

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] ESD模块计算ISC速率手册算例计算出错求指教

评分 Rate

浏览过的版块

[ORCA] ESD模块计算ISC速率 手册算例计算出错求指教

评分 Rate

浏览过的版块

[ORCA] ESD模块计算ISC速率手册算例计算出错求指教