ORCA 5.0.3优化激发态时报错：ORCA finished by error termination in CIS

j5888xm

各位老师，我在使用ORCA5.0.3优化激发态S1min时，报错如下，苦苦得不到解决办法：


[file orca_tools/Tool-Linalg/qcrpasolv.cpp, line 785, Process 3]: Rebuilding of the RPA subspace failed because there are no old eigenvectors, aborting

ORCA finished by error termination in CIS
Calling Command: mpirun -np 8  /export/home/test/app/orca_5_0_3_linux_x86-64_shared_openmpi411/orca_cis_mpi Carbazole_imine_S1_OPT.cisinp.tmp Carbazole_imine_S1_OPT
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run

我在计算基态小分子以及优化过渡态时，任务均能正常结束，但是计算激发态极小值优化时，总会报错，有时算了几个小时优化三十多步后报错，有时直接报错

服务器的内存是512GB，输入和输出均贴在下方，望求大佬们不吝赐教，十分感谢

sobereva

不要自己在标题里手写[ORCA]这种标签，直接选相应分类就完了。这次给你改了，以后注意

wzkchem5

“Rebuilding of the RPA subspace failed because there are no old eigenvectors”表示基态波函数不稳定。虽然可以找到稳定的波函数再继续做TDDFT结构优化，但是结果会有严重的自旋污染。所以此时最好重新评估一下所用泛函从理论角度讲是否适合这个体系，如果不适合，换泛函；如果所有泛函都不适合，考虑用spin flip TDA优化结构（当然如果算得动的话CASPT2等方法更好）；如果所有泛函都不适合，且体系不适合做spin flip TDA，也可以做普通TDA计算，但结果可靠性较低

j5888xm

sobereva 发表于 2023-9-8 23:22
不要自己在标题里手写[ORCA]这种标签，直接选相应分类就完了。这次给你改了，以后注意

好的老师，下次一定注意

j5888xm

wzkchem5 发表于 2023-9-9 02:13
“Rebuilding of the RPA subspace failed because there are no old eigenvectors”表示基态波函数不稳定 ...

嗯嗯老师，您的意思是这可能是PBE0并不适用这个体系是吗，那有没有您所推荐的泛函更适合激发态计算呢？我是想用TD先定性看一下激发态然后用CASSCF优化的，我想先看一下激发态结构

wzkchem5

j5888xm 发表于 2023-9-9 08:45
嗯嗯老师，您的意思是这可能是PBE0并不适用这个体系是吗，那有没有您所推荐的泛函更适合激发态计算呢？我 ...

看http://sobereva.com/272，针对你的体系判断

15689900879@163

wzkchem5 发表于 2023-9-9 02:13
“Rebuilding of the RPA subspace failed because there are no old eigenvectors”表示基态波函数不稳定 ...

老师，您好！我也遇到这个问题了。我的T1文件是在Gaussian 16中优化和频率计算正常结束的，但ORCA同样报错，Rebuilding of the RPA subspace failed because there are no old eigenvectors, aborting。

wzkchem5

15689900879@163 发表于 2024-3-21 02:10
老师，您好！我也遇到这个问题了。我的T1文件是在Gaussian 16中优化和频率计算正常结束的，但ORCA同样报 ...

你的高斯计算是怎么做的？用的UKS还是TDDFT？开了TDA了吗？做了broken symmetry吗？

15689900879@163

wzkchem5 发表于 2024-3-21 15:53
你的高斯计算是怎么做的？用的UKS还是TDDFT？开了TDA了吗？做了broken symmetry吗？

老师，您好！高斯命令# opt freq td=(nstates=5,triplet) def2tzvp em=gd3bj pbe1pbe 0 1。broken symmetry是手动破坏原始构象，并添加nosymm命令实现吗？老师，我是否应该分别尝试UKS、TDA和broken symmetry，重新优化构象，再计算SOC呢？

Excited State   1:      Triplet-A      2.1490 eV  576.93 nm  f=0.0000  <S**2>=2.000
      31 -> 33         0.26438
      32 -> 33         0.64627
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -385.913491864   
Copying the excited state density for this state as the 1-particle RhoCI density.

wzkchem5

15689900879@163 发表于 2024-3-22 02:21
老师，您好！高斯命令# opt freq td=(nstates=5,triplet) def2tzvp em=gd3bj pbe1pbe 0 1。broken symmet ...

broken symmetry指的是破坏自旋对称性，不是破坏空间对称性。做法参见http://sobereva.com/82，也可以在orca里用%scf brokensym 1,1 end之类的关键词得到broken symmetry波函数，之后可以传给高斯。
你给的高斯输出结果和orca计算根本就不是一个结构吧，你的orca计算的结构有12个氟原子，一个氟原子的绝对能量就将近-100 au了，激发态绝对能量不可能是-385.9au

15689900879@163

wzkchem5 发表于 2024-3-22 17:10
broken symmetry指的是破坏自旋对称性，不是破坏空间对称性。做法参见http://sobereva.com/82，也可以在o ...

是的，老师，我复制错文件了。感谢您的回复，我再研究研究，谢谢老师。

Excited State   1:      Triplet-A      0.9399 eV 1319.12 nm  f=0.0000  <S**2>=2.000
      63 -> 77         0.12676
      73 -> 77        -0.54453
      75 -> 77         0.72333
      73 <- 77        -0.36962
      75 <- 77         0.45910
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -1501.06774366   
Copying the excited state density for this state as the 1-particle RhoCI density.

526827797

各位老师好，我这也遇到了类似的报错，但是没有
“Rebuilding of the RPA subspace failed because there are no old eigenvectors”这句提示，求教如何解决

ORCA finished by error termination in CIS
Calling Command: mpirun -np 256 --allow-run-as-root /work/home/acrsgru5zp/apprepo/orca/5.0.4-openmpi411/app/orca_cis_mpi TAPT-F-TD.cisinp.tmp TAPT-F-TD
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run