计算化学公社

标题: orca查找锥形交叉点失败 [打印本页]

作者
Author: BinWang 时间: 2025-9-16 19:44
标题: orca查找锥形交叉点失败
本帖最后由 BinWang 于 2025-9-17 16:43 编辑

各位老师：
   我用orca计算一个体系的锥形交叉点老是失败，哪位老师能不能帮我看一下什么原因？这是输入文件： ! ci-opt  pbe0 def2-SVP NORI  miniprint
%maxcore 30000
#%pal
#nproc 8
#end
!PAL8
%tddft
  sf true
  nroots 5
  iroot 3
  jroot 1
  followiroot true
  # IROOTMULT TRIPLET
  # TDA false
end
%geom
maxstep 0.05
trust -0.05
end
%cpcm smd true
   SMDsolvent "water"
end
*xyzfile 0 3 refer.xyz
这个是输出的报错：
------------------------------------------------------------------------------------------
  1-3A  ->  2-3A 3.495698 28194.7 354.7 0.00001  -0.00138  -0.00319  -0.00000
  1-3A  ->  3-1A 4.014915 32382.5 308.8 0.00331  -0.38861  -0.89494  -0.00000
  1-3A  ->  4-3A 4.293602 34630.2 288.8 -0.00072 0.00000 0.00000  -0.09231
  1-3A  ->  5-3A 4.564809 36817.7 271.6 -0.00216  -0.00000  -0.00001 0.12617

Total run time:    50.636 sec

         *** ORCA-CIS/TD-DFT FINISHED WITHOUT ERROR ***

Saving reference transition density.

-----------------------
CIS/TD-DFT TOTAL ENERGY
-----------------------

E(SCF)  = -907.359555578 Eh
DE(CIS) =    0.047919672 Eh (Root  3)
----------------------------- ---------
E(tot)  = -907.311635906 Eh

                        ------------------
                           SINGLET ROOT 3
                        ------------------

--------------------------
TD-DFT/TDA RELAXED DENSITY
--------------------------

Reference wavefunction          ... UHF
State of interest                ... 3
Basis dimension                   ... 322
DFT-XC contribution needed       ... YES
RI-approximation                ... NO
Do Coulomb                      ... YES
Do Exchange                      ... YES
Need Hartree-Fock Exchange       ... YES
facj, facx                      ...  1.000, -0.250
Conventional                      ... NO
Range separation                ... NO
Triplet state                   ... NO
Spin Flip                         ... YES
CIS-vectors being read from       ... CI_orca.cis
=> Number of Vectors found: 5
=> Orbital Window: 21...63 -> 64...321
                  21...61 -> 62...321

Reading the CIS file                   ... done
Norm of the CI vector                   ... 1.000000000
Reading the CPCM file                   ... done
Making the CIS (pseudo)densities       ... done
Making operators G(D) and Q(T)          ... (SHARK-ok)
ORCA finished by error termination in CIS
Calling Command: mpirun -np 8  /public/home/wk/soft/orca610/bin/orca_cis_mpi CI_orca.cisinp.tmp CI_orca
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run

我把输入文件跟输出文件都上传了，大家有什么建议？
这个体系用GAMMES是能收敛的，但是GAMMES有两个缺点：1，锥形交叉点搜索不包含溶剂效应。2，有些泛函GAMMES不支持，比如wb97xd泛函（orca有wb97xd3代替），麻烦哪位老师帮我想想办法。我的orca版本是6.1，openmpi是(Open MPI) 4.1.6版本

作者
Author: wzkchem5 时间: 2025-9-16 20:54
S态和T态只有MECP，没有锥形交叉点

作者
Author: BinWang 时间: 2025-9-17 08:36
本帖最后由 BinWang 于 2025-9-17 09:03 编辑

wzkchem5 发表于 2025-9-16 20:54
S态和T态只有MECP，没有锥形交叉点

非常感谢王老师！

作者
Author: BinWang 时间: 2025-9-17 12:12

wzkchem5 发表于 2025-9-16 20:54
S态和T态只有MECP，没有锥形交叉点

非常感谢王老师的回复！我改成这样行不行啊？改成这样没有出现报错，得到的结果是S0/S1的CI吗？
! ci-opt  pbe0 def2-SVP NORI  miniprint
%maxcore 30000
#%pal
#nproc 8
#end
!PAL8
%tddft
  sf true
  nroots 5
  iroot 2
  jroot 1
  followiroot true
  # IROOTMULT TRIPLET
  # TDA false
end
%geom
maxstep 0.05
trust -0.05
end
%cpcm smd true
   SMDsolvent "water"
end
*xyzfile 0 3 refer.xyz

作者
Author: wzkchem5 时间: 2025-9-17 12:42

BinWang 发表于 2025-9-17 12:12
非常感谢王老师的回复！我改成这样行不行啊？改成这样没有出现报错，得到的结果是S0/S1的CI吗？
! ci-o ...

需要检查激发态组成才能确认。jroot 1一般对应的是S0没错，但iroot 2可能对应的是T1的Ms=0组分

作者
Author: BinWang 时间: 2025-9-17 14:22

wzkchem5 发表于 2025-9-17 12:42
需要检查激发态组成才能确认。jroot 1一般对应的是S0没错，但iroot 2可能对应的是T1的Ms=0组分

非常感谢王老师！

作者
Author: BinWang 时间: 2025-9-17 16:14

wzkchem5 发表于 2025-9-17 12:42
需要检查激发态组成才能确认。jroot 1一般对应的是S0没错，但iroot 2可能对应的是T1的Ms=0组分

王老师您好！这个帖子（https://bane-dysta.top/posts/CIopt/#mrsf-tddft）里面有用orca在SF-TDDFT下计算S0/S1的例子，我之前是按照这个例子做的，出现了报错。这个帖子里面确定激发态组成的方法正确吗？我最初的报错也是重复这个例子出来的。

作者
Author: wzkchem5 时间: 2025-9-17 16:29

BinWang 发表于 2025-9-17 16:14
王老师您好！这个帖子（https://bane-dysta.top/posts/CIopt/#mrsf-tddft）里面有用orca在SF-TDDFT下计算 ...

报错内容是什么？用的ORCA版本是什么？分子结构和这个例子严格一致吗？

作者
Author: BinWang 时间: 2025-9-17 16:45

wzkchem5 发表于 2025-9-17 16:29
报错内容是什么？用的ORCA版本是什么？分子结构和这个例子严格一致吗？

非常感谢王老师的回复，我已经把计算的结构，输入文件，输出文件传到帖子最上面了。麻烦您看一眼。

作者
Author: wzkchem5 时间: 2025-9-17 16:55
是你改过还是我一开始看错了？可能我看到你写了iroot和jroot，下意识以为你算的是S态和T态的锥形交叉点。
如果是一个spin flip计算，写一个iroot，一个jroot，那么只要态指认正确了就没问题。按<S**2>指认态就行了。
这样的话报错可能是内存问题，例如maxcore给太多了

作者
Author: BinWang 时间: 2025-9-17 16:58
本帖最后由 BinWang 于 2025-9-17 16:59 编辑

wzkchem5 发表于 2025-9-17 16:55
是你改过还是我一开始看错了？可能我看到你写了iroot和jroot，下意识以为你算的是S态和T态的锥形交叉点。
...

我这个输入文件没改过，非常感谢王老师

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