ORCA5.0.3报错:Rebuilding of the RPA subspace failed

cokie

各位老师好，我在使用ORCA 5.0.3 win版本程序计算自旋-轨道耦合常数时，发生如下报错

1. [file orca_tools/Tool-Linalg/qcrpasolv.cpp, line 785, Process 0]: Rebuilding of the RPA subspace failed because there are no old eigenvectors, aborting
   
2. 
   
3. [file orca_tools/Tool-Linalg/qcrpasolv.cpp, line 785, Process 10]: Rebuilding of the RPA subspace failed because there are no old eigenvectors, aborting
   
4. 
   
5. [file orca_tools/Tool-Linalg/qcrpasolv.cpp, line 785, Process 15]: Rebuilding of the RPA subspace failed because there are no old eigenvectors, aborting
   
6. 
   
7. [file orca_tools/Tool-Linalg/qcrpasolv.cpp, line 785, Process 13]: Rebuilding of the RPA subspace failed because there are no old eigenvectors, aborting
   
8. 
   
9. [file orca_tools/Tool-Linalg/qcrpasolv.cpp, line 785, Process 14]: Rebuilding of the RPA subspace failed because there are no old eigenvectors, aborting
   
10. 
    
11. 
    
12. job aborted:
    
13. [ranks] message
    
14. 
    
15. [0-15] process exited without calling finalize
    
16. 
    
17. ---- error analysis -----
    
18. 
    
19. [0-15] on DESKTOP-NFN1R15
    
20. C:\ORCA\orca_cis_mpi.exe ended prematurely and may have crashed. exit code 78
    
21. 
    
22. ---- error analysis -----
    
23. 
    
24. ORCA finished by error termination in CIS
    
25. Calling Command: mpiexec -np 16  C:\ORCA\orca_cis_mpi.exe PXD_soc.cisinp.tmp PXD_soc
    
26. [file orca_tools/qcmsg.cpp, line 465]:
    
27.   .... aborting the run

复制代码

我自己做过几个对比：更换泛函/基组（辅助基组）、增加nroots（之前看sob老师在某个帖子给过的可能有用的建议）、不使用TDA FALSE关键词、将并行核数改为12(每个核3500 MB，我的机器是16核64GB)均有上述报错。但我更换分子发现并没有上述报错，计算可以正常进行，所以应该也不是mpi.exe出了问题。但我不是很明白我的分子有什么特别之处。我把输入文件、输出文件和分子结构（用gaussian做过S1的优化）放下面了，希望能得到各位老师的帮助和解惑，非常感谢！

wzkchem5

不要只看这行报错信息，关键在于下面这行
Negative eigenvalue found in transformed Hamiltonian
这个代表波函数不稳定。导致RPA subspace rebuild的主要有两个可能的原因，一个是波函数不稳定，一个是内存不足，增加nroots、改基组、改核数基本只对后者有效。
打开TDA以后按理说不应该报这个错了，你确定报错信息也是这个？

cokie

wzkchem5 发表于 2024-3-26 16:48
不要只看这行报错信息，关键在于下面这行
Negative eigenvalue found in transformed Hamiltonian
这个代 ...

是的老师，不过不使用TDA FALSE的输出文件我手欠删了……就留了这最后一次。那请问老师这种情况是不是只能用gaussian做stable=opt，然后把chk文件转成wfn格式作为orca初猜才有可能解决呢？orca有没有类似stable=opt这种方式来找到自己稳定的基函数呢？

wzkchem5

cokie 发表于 2024-3-26 09:52
是的老师，不过不使用TDA FALSE的输出文件我手欠删了……就留了这最后一次。那请问老师这种情况是不是只 ...

orca可以写
%scf
stabperform true
end
等价于高斯的stable=opt。
另外注意，视情况而定，找到稳定波函数再做计算未必是最理想的解决方法，参见http://bbs.keinsci.com/forum.php ... 6&fromuid=21811

cokie

wzkchem5 发表于 2024-3-26 16:55
orca可以写
%scf
stabperform true

明白了，感谢老师指点，我根据您的建议和推荐的博文尝试一下。

另外就是orca好像没有现成的wb97x带D2矫正的版本，我如果想让激发能尽可能和gaussian中使用wb97xd匹配，是否选用wb97x会比wb97x-d3更好一些呢？还是说计算激发能，带不带D矫正的泛函对所计算的结果并没有什么明显的影响呢？

wzkchem5

cokie 发表于 2024-3-26 09:58
明白了，感谢老师指点，我根据您的建议和推荐的博文尝试一下。

另外就是orca好像没有现成的wb97xd的D2 ...

一般来说不加色散校正对激发能影响精确为0。但是wB97X系列泛函是例外，因为wB97X某种意义上类似M06-2X，考虑了一部分中程色散，所以但凡要加色散校正，就会发现最好要重新拟合wB97X部分的参数，以避免double-count色散部分。所以wB97X、wB97XD、wB97X-D3、wB97X-D3BJ等等的wB97X部分的参数都不一样，导致激发能结果有一些区别，不过一般区别不大。

cokie

wzkchem5 发表于 2024-3-26 17:04
一般来说不加色散校正对激发能影响精确为0。但是wB97X系列泛函是例外，因为wB97X某种意义上类似M06-2X， ...

明白了，非常感谢老师的帮助！