计算化学公社

标题: ORCA finished by error termination in CIS [打印本页]

作者
Author: BinWang 时间: 2025-9-26 18:13
标题: ORCA finished by error termination in CIS
本帖最后由 BinWang 于 2025-9-26 18:38 编辑

各位老师：
   我用orca6.1计算一个体系的时候有如下报错，我的mpi版本是mpirun (Open MPI) 4.1.6,我换了好几种%maxcore的值进行计算，都报了同样的错误，我的服务器内存有512G，足够开销。同时，我换了说明书中推荐的泛函还是报同样的错误。我的输入文件和输出文件见附件，哪位老师能不能帮我看一眼有没有可能是别的错误？
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 input_eg.cisinp.tmp input_eg
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run

作者
Author: wzkchem5 时间: 2025-9-26 21:21
同一台服务器是否曾经成功计算过其他分子？
能成功计算的分子，和这个分子最大的区别是什么？

作者
Author: BinWang 时间: 2025-9-27 09:23
本帖最后由 BinWang 于 2025-9-27 09:26 编辑

wzkchem5 发表于 2025-9-26 21:21
同一台服务器是否曾经成功计算过其他分子？
能成功计算的分子，和这个分子最大的区别是什么？

非常感谢王老师的回复，这个同一台机器上成功算过不同的分子。区别是这个分子双自由基特征明显，用普通泛函算的话居然出现了S1的能量低于S0的情况，还有用普通泛函计算的话，这个分子的T1略低于S0。之前直接用def2-svp基组，貌似单点能可以成功收敛。

作者
Author: BinWang 时间: 2025-9-27 09:58
本帖最后由 BinWang 于 2025-9-27 10:07 编辑

wzkchem5 发表于 2025-9-26 21:21
同一台服务器是否曾经成功计算过其他分子？
能成功计算的分子，和这个分子最大的区别是什么？

王老师您好！我试了一下，改成这样能正常结束，但我还是希望能用原来设定的基组完成计算：
! wb97x-d3bj def2-svp  NoTRAH

%maxcore 12000

%pal
  nproc 8
end

%tddft
  sf true
  nroots 5
  end

%cpcm smd true
      SMDsolvent "water"
end

  *xyzfile 1 3 ring_o.xyz

输出文件如下：

the weight of the individual excitations are printed if larger than 1.0e-02

UHF/UKS reference: multiplicity estimated based on rounded <S**2> value, RELEVANCE IS LIMITED!

(SPIN-FLIP GROUND STATE)
STATE  1:  E= 0.022064 au    0.600 eV    4842.6 cm**-1 <S**2> = 0.777954 Mult 3
150a -> 194b  :    0.022102 (c= -0.14866893)
156a -> 194b  :    0.012658 (c= -0.11250666)
162a -> 194b  :    0.019322 (c= -0.13900330)
163a -> 194b  :    0.012680 (c=  0.11260362)
164a -> 194b  :    0.021937 (c= -0.14811263)
165a -> 194b  :    0.012085 (c= -0.10992978)
166a -> 194b  :    0.010013 (c=  0.10006253)
167a -> 194b  :    0.025020 (c= -0.15817858)
170a -> 194b  :    0.025632 (c=  0.16009862)
175a -> 194b  :    0.205788 (c=  0.45363853)
175a -> 195b  :    0.018706 (c=  0.13676902)
176a -> 194b  :    0.116147 (c=  0.34080300)
180a -> 194b  :    0.055148 (c=  0.23483603)
181a -> 194b  :    0.010426 (c= -0.10210565)
183a -> 194b  :    0.048321 (c= -0.21982008)
185a -> 194b  :    0.019573 (c=  0.13990470)
192a -> 194b  :    0.048279 (c=  0.21972581)
193a -> 194b  :    0.152779 (c= -0.39086978)

STATE  2:  E= 0.037641 au    1.024 eV    8261.3 cm**-1 <S**2> = 1.142444 Mult 3
175a -> 194b  :    0.011120 (c= -0.10545377)
176a -> 195b  :    0.024811 (c= -0.15751632)
189a -> 194b  :    0.015821 (c= -0.12578025)
189a -> 195b  :    0.056765 (c=  0.23825498)
190a -> 195b  :    0.014376 (c=  0.11990131)
191a -> 195b  :    0.019036 (c=  0.13796981)
192a -> 194b  :    0.024818 (c=  0.15753656)
192a -> 195b  :    0.126412 (c= -0.35554478)
193a -> 194b  :    0.101845 (c= -0.31913217)
193a -> 195b  :    0.468311 (c=  0.68433269)

STATE  3:  E= 0.068212 au    1.856 eV 14970.9 cm**-1 <S**2> = 1.087365 Mult 3
131a -> 194b  :    0.020665 (c=  0.14375397)
143a -> 194b  :    0.018306 (c=  0.13529806)
146a -> 194b  :    0.010705 (c= -0.10346356)
152a -> 194b  :    0.072983 (c=  0.27015339)
158a -> 194b  :    0.022192 (c=  0.14896819)
159a -> 194b  :    0.024676 (c=  0.15708491)
162a -> 194b  :    0.018836 (c= -0.13724268)
163a -> 194b  :    0.120151 (c= -0.34662764)
163a -> 195b  :    0.014716 (c= -0.12131119)
165a -> 194b  :    0.119353 (c=  0.34547552)
165a -> 195b  :    0.011904 (c=  0.10910483)
167a -> 194b  :    0.067707 (c=  0.26020647)
168a -> 194b  :    0.011673 (c= -0.10804194)
175a -> 194b  :    0.033293 (c= -0.18246249)
176a -> 194b  :    0.048002 (c=  0.21909339)
180a -> 194b  :    0.019549 (c= -0.13981813)
192a -> 194b  :    0.022072 (c=  0.14856495)
193a -> 194b  :    0.070550 (c= -0.26561182)
193a -> 195b  :    0.023010 (c= -0.15168958)
195a -> 194b  :    0.021403 (c=  0.14629627)

STATE  4:  E= 0.078839 au    2.145 eV 17303.2 cm**-1 <S**2> = 0.445288 Mult 1
150a -> 194b  :    0.017905 (c= -0.13381089)
156a -> 194b  :    0.011534 (c= -0.10739486)
162a -> 194b  :    0.052323 (c= -0.22874248)
175a -> 194b  :    0.016509 (c=  0.12848590)
183a -> 194b  :    0.062113 (c= -0.24922388)
184a -> 194b  :    0.011510 (c= -0.10728658)
186a -> 194b  :    0.010591 (c= -0.10291401)
189a -> 194b  :    0.038896 (c=  0.19722036)
190a -> 194b  :    0.026989 (c=  0.16428247)
192a -> 194b  :    0.035906 (c= -0.18948800)
192a -> 195b  :    0.027471 (c= -0.16574238)
193a -> 194b  :    0.304959 (c=  0.55223136)
193a -> 195b  :    0.107042 (c=  0.32717329)
195a -> 194b  :    0.103903 (c=  0.32233970)

STATE  5:  E= 0.085906 au    2.338 eV 18854.2 cm**-1 <S**2> = 1.042405 Mult 3
162a -> 194b  :    0.012102 (c=  0.11001073)
181a -> 194b  :    0.010164 (c=  0.10081680)
183a -> 194b  :    0.010997 (c=  0.10486760)
185a -> 195b  :    0.012834 (c= -0.11328731)
187a -> 194b  :    0.010301 (c=  0.10149191)
191a -> 194b  :    0.018884 (c= -0.13742009)
192a -> 194b  :    0.034887 (c=  0.18678081)
193a -> 194b  :    0.021250 (c= -0.14577476)
194a -> 194b  :    0.082825 (c= -0.28779382)
195a -> 194b  :    0.535849 (c=  0.73201717)
195a -> 195b  :    0.076653 (c= -0.27686218)

Storing amplitudes in GBW file ...
... GBW file created

-----------------------------
TD-DFT/TDA-EXCITATION SPECTRA
-----------------------------

作者
Author: BinWang 时间: 2025-9-27 11:04

wzkchem5 发表于 2025-9-26 21:21
同一台服务器是否曾经成功计算过其他分子？
能成功计算的分子，和这个分子最大的区别是什么？

王老师您好:
   我用def2svp基组优化的时候也报同样的错，输入文件如下：
  1 ! wb97x-d3bj def2-svp opt
2 !NoTRAH
3
4 %maxcore 30000
5
6 %pal
7  nproc 8
8 end
9
10 %tddft
11  sf true
12  nroots 5
13  end
14
15 %cpcm smd true
16       SMDsolvent "water"
17 end
18
19  *xyzfile 1 3 ring_o.xyz
20
21
输出文件如下：
CIS/TD-DFT TOTAL ENERGY
-----------------------

E(SCF)  =  -2964.856177861 Eh
DE(CIS) =    0.022909315 Eh (Root  1)
----------------------------- ---------
E(tot)  =  -2964.833268546 Eh

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

Reference wavefunction          ... UHF
State of interest                ... 1
Basis dimension                   ... 932
DFT-XC contribution needed       ... YES
RI-approximation                ... YES
Do Coulomb                      ... YES
Do Exchange                      ... YES
Need Hartree-Fock Exchange       ... YES
facj, facx                      ...  1.000, -0.167
Conventional                      ... NO
Range separation                ... YES
Triplet state                   ... NO
Spin Flip                         ... YES
CIS-vectors being read from       ... input_eg.cis
=> Number of Vectors found: 5
=> Orbital Window: 61...195 -> 196...931
                  61...193 -> 194...931

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 input_eg.cisinp.tmp input_eg
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run

t_eg.log" 5903L, 252901C

作者
Author: wzkchem5 时间: 2025-9-27 12:52

BinWang 发表于 2025-9-27 09:23
非常感谢王老师的回复，这个同一台机器上成功算过不同的分子。区别是这个分子双自由基特征明显，用普通泛 ...

成功计算的分子和这个分子在大小上有区别吗

作者
Author: BinWang 时间: 2025-9-27 14:30
本帖最后由 BinWang 于 2025-9-27 14:33 编辑

wzkchem5 发表于 2025-9-27 12:52
成功计算的分子和这个分子在大小上有区别吗

王老师您好！我现在的体系84个原子。我唯一跑成功的例子是说明书里面苯环双自由基的例子。我将说明书里面的例子改成wb97x-d3bj计算，报出同样错误，以下是输入文件：
!DEF2-TZVPD OPT wb97x-d3bj
%maxcore  2500
%pal nprocs 8 end
%tddft
SF true
NRoots 5
end
* xyz 0 3
C -1.39113 0.00000 0.00000
C 0.69557 1.20476 0.00000
C -0.69557 1.20476 0.00000
C -0.69557 -1.20476 0.00000
C 0.69557 -1.20476 0.00000
C 1.39113 0.00000 0.00000
H -1.24291 2.15278 0.00000
H -1.24291 -2.15278 0.00000
H 1.24291 -2.15278 0.00000
H 1.24291 2.15278 0.00000
*
以下是输出文件：

Reading the CIS file                   ... done
Norm of the CI vector                   ... 1.000000000
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 eg.cisinp.tmp eg
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run

我目前打算用说明书里面一模一样的方法添加溶剂效应计算一下我当前的体系。为了保持泛函一致性，我非常希望用wb97x-d3.

作者
Author: wzkchem5 时间: 2025-9-27 16:31

BinWang 发表于 2025-9-27 14:30
王老师您好！我现在的体系84个原子。我唯一跑成功的例子是说明书里面苯环双自由基的例子。我将说明书里面 ...

试试用libxc里的wb97x-d3泛函，输入文件书写方法在orca手册的libxc相关章节里有

作者
Author: BinWang 时间: 2025-9-27 17:43

wzkchem5 发表于 2025-9-27 16:31
试试用libxc里的wb97x-d3泛函，输入文件书写方法在orca手册的libxc相关章节里有

非常感谢王老师

，我用说明书里推荐的方法搭配def2-svp，设定Maxcores=12000，以及nproc=16的情况下，目前已经优化了两轮了。到后面我会用优化好的结构进行wb97x-d3的计算。

作者
Author: BinWang 时间: 2025-9-28 10:43

wzkchem5 发表于 2025-9-27 16:31
试试用libxc里的wb97x-d3泛函，输入文件书写方法在orca手册的libxc相关章节里有

王老师您好：我的体系直接用说明书推荐的泛函 hyb_gga_xc_bhandhlyp +def2-svp基组能正常收敛，输入文件如下：
! DEF2-SVP OPT
%maxcore 12000
%pal
  nproc 16
end
%method
  Method DFT
  Functional hyb_gga_xc_bhandhlyp
end
%tddft
  sf true
  nroots 5
  end

%cpcm smd true
      SMDsolvent "water"
end

  *xyzfile 1 3 ring_o.xyz
当我换成自定义基组的时候总是报错,输入文件如下：
!  OPT NoRI

%maxcore 30000

%pal
nproc 16
end
%basis
  NewGTO H "6-31G**" end
  NewGTO C "6-31G**" end
  NewGTO N "6-31G**" end
  NewGTO O "6-31G**" end
  NewGTO S "6-31G**" end
  NewGTO Ir
  S 3
  1       15.2937090          -1.6750339
  2       13.5736820             2.3934648
  3       5.8162740          -1.6504979
  S 1
  1       1.1955210             1.0000000
  S 1
  1       0.5657760             1.0000000
  S 1
  1       0.1405270             1.0000000
  S 1
  1       0.0508200             1.0000000
  S 1
  1       0.0150000             1.0000000
  P 2
  1       8.6697960             0.8677142
  2       6.2456140          -1.8144826
  P 2
  1       1.9663840             0.2709454
  2       1.0202060             0.7540443
  P 1
  1       0.4065800             1.0000000
  P 1
  1       0.0901160             1.0000000
  P 1
  1       0.0273640             1.0000000
  D 4
  1       3.6994850          -0.5973104
  2       3.3617430             0.6112252
  3       1.1745250             0.4732939
  4       0.4701760             0.5755302
  D 1
  1       0.1576780             1.0000000
  D 1
  1       0.0500000             1.0000000
  end
  NewECP Ir "SDD" end
  NewGTO 7 "6-31++G**" end
  NewGTO 83 "6-31++G**" end
end
%method
Method DFT
Functional hyb_gga_xc_bhandhlyp
end

%tddft
sf true
nroots 5
end

%cpcm smd true
   SMDsolvent "water"
end

*xyzfile 1 3 ring_o.xyz
我尝试关闭RI加速，设置过不同的maxcores数值，都失败了。输出文件如下：
Total run time:    898.045 sec

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

E(SCF)  =  -2964.444196441 Eh
DE(CIS) =    0.015789632 Eh (Root  1)
----------------------------- ---------
E(tot)  =  -2964.428406809 Eh

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

Reference wavefunction          ... UHF
State of interest                ... 1
Basis dimension                   ... 952
DFT-XC contribution needed       ... YES
RI-approximation                ... NO
Do Coulomb                      ... YES
Do Exchange                      ... YES
Need Hartree-Fock Exchange       ... YES
facj, facx                      ...  1.000, -0.500
Conventional                      ... NO
Range separation                ... NO
Triplet state                   ... NO
Spin Flip                         ... YES
CIS-vectors being read from       ... input_eg.cis
=> Number of Vectors found: 5
=> Orbital Window: 61...195 -> 196...951
                  61...193 -> 194...951

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 16  /public/home/wk/soft/orca610/bin/orca_cis_mpi input_eg.cisinp.tmp input_eg
[file orca_tools/qcmsg.cpp, line 394]:
  .... aborting the run
还有一个玄学的地方是同一个体系，原子完全相同，结构不同的中间体用SF-TDDFT优化，有极个别能收敛。有些就报CIS错误。针对这些问题您有什么建议吗？

作者
Author: zjxitcc 时间: 2025-9-28 11:17
你能直接上传xyz文件吗，你的inp文件里没给坐标，用的是引用xyz文件，读者看到也没法测试，仅凭ORCA那点报错信息很难看出什么。

比如说就针对10L你说的问题，其对应的xyz, inp, out文件放在一起打包压缩上传，前面的问题也暂时也可以不用看了。

作者
Author: BinWang 时间: 2025-9-28 13:49
本帖最后由 BinWang 于 2025-9-28 13:51 编辑

zjxitcc 发表于 2025-9-28 11:17
你能直接上传xyz文件吗，你的inp文件里没给坐标，用的是引用xyz文件，读者看到也没法测试，仅凭ORCA那点报 ...

非常感谢老师的回复，这是和实验合作的项目。实验的人怕发出来以后被抢发，所以目前他们不愿意上传分子结构。荧光探针这块竞争蛮激烈的。真的是一言难尽，同一个课题组都出现了抢发的事情。。。。以后遇到我自己的课题，我会尽量的上传分子结构。

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