ORCA计算闭壳层阳离子的单点能，N-1态计算正常，N+1态报错

faylovesnow

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!                        FATAL ERROR ENCOUNTERED               !!!
!!!                        -----------------------               !!!
!!!                           BLAS-ERROR                         !!!
!!!              INCOMPATIBLE MATRICES OR VECTORS FOUND          !!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
--------------------------------------------------------------------------
Primary job  terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun detected that one or more processes exited with non-zero status, thus causing
the job to be terminated. The first process to do so was:


  Process name: [[26859,1],19]
  Exit code:    62
--------------------------------------------------------------------------
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run




ORCA finished by error termination in SCF
Calling Command: mpirun -np 36  /home/yzh/software/orca504/orca_scf_mpi 805228-31-3_sp.gbw b 805228-31-3_sp
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run




如题所述，ORCA计算+1电荷闭壳层阳离子的单点能，N-1态计算正常，N+1态报错。请各位老师看看是哪里出了问题。
以下是INP文件内容：
! DLPNO-CCSD(T) tightPNO RIJK cc-pVTZ cc-pVTZ/JK cc-pVTZ/C tightSCF noautostart miniprint nopop CPCM
%maxcore  10000
%pal nprocs  36 end
%CPCM
smd true
SMDsolvent "water"
End
* xyz   0   2
  C   -1.00406433677227     -0.44464228437832     -0.70086563040498
  C   0.30393529492131     -0.27625616748073     -1.01892473784876
  C   -0.00149610852417      0.46628334852234      1.04015404531691
  N   -1.17096208467125      0.02659929271298      0.58715890848283
  H   -1.82230455577689     -0.85690487248229     -1.26369935101363
  H   0.85266761347034     -0.50455981368858     -1.91521105151881
  H   0.17891922000365      0.89271571247415      2.01231806130703
  N   0.91413780441911      0.28954399727176      0.08369292924342
  C   -2.42956465227257      0.03087581446326      1.33377026973398
  H   -2.25365409148362      0.48485958031023      2.30438162527478
  H   -3.16671096901293      0.60850283761543      0.78086912352061
  H   -2.77376728443910     -0.99331334525001      1.45871930883628
  C   2.30496028886682      0.70904835442837      0.16202860679543
  H   2.61325619054600      0.71152582982866      1.20732817342895
  H   2.92102852051098     -0.01300137250481     -0.36998054434060
  C   2.58373013393250      2.09041899600993     -0.41174916539625
  O   3.71208998247670      2.50745827857752     -0.51946256686014
  O   1.48195598232918      2.75027732967683     -0.74354001733250
  C   1.66810437759881      4.08637883293682     -1.27712234799640
  H   2.17310318086040      4.70868586621397     -0.54146117583784
  H   2.25196570149843      4.03806219135023     -2.19386200171954
  H   0.66701279151859      4.45506459339226     -1.47477846167073
*

wzkchem5

估计是硬盘满了，或者偶然的硬件问题。重新算一遍可能就好了

faylovesnow

wzkchem5 发表于 2024-6-21 00:55
估计是硬盘满了，或者偶然的硬件问题。重新算一遍可能就好了

不是硬盘的问题呢，硬盘4T的。在该计算级别下，N+1态算了很多遍了，一直报错。换 RI-PWPB95-D3(BJ)/def2-QZVPP 计算没问题。

zjxitcc

faylovesnow 发表于 2024-6-21 08:14
不是硬盘的问题呢，硬盘4T的。在该计算级别下，N+1态算了很多遍了，一直报错。换 RI-PWPB95-D3(BJ)/def2- ...

（1）硬盘4T，不代表4T都能用，ORCA能用的空间是当前分区下剩余空间大小。
（2）既然坐标都给出来了，说明体系没什么保密性，干脆上传压缩后的输出文件。

faylovesnow

zjxitcc 发表于 2024-6-21 08:51
（1）硬盘4T，不代表4T都能用，ORCA能用的空间是当前分区下剩余空间大小。
（2）既然坐标都给出来了，说 ...

感谢大牛回帖，我通过论坛学习，在关键词里加入notrah 解决了问题，好像是默认20次迭代不收敛后自动开启了trah导致的报错。我加入notrah后30轮收敛了。

wzkchem5

faylovesnow 发表于 2024-6-21 02:28
感谢大牛回帖，我通过论坛学习，在关键词里加入notrah 解决了问题，好像是默认20次迭代不收敛后自动开启 ...

所以说明问问题报错信息要给全，如果你多给几十行输出，我们发现是在进入TRAH之后报错的，就更可能猜到真实报错原因

faylovesnow

wzkchem5 发表于 2024-6-21 15:10
所以说明问问题报错信息要给全，如果你多给几十行输出，我们发现是在进入TRAH之后报错的，就更可能猜到真 ...

是的，下次我会注意，把所有信息提供一下。