Orca_5.0.3报错no orbitals were found in the gbw file解决过程记录，期望老师们指正

zhaoxiaoni

zjxitcc 发表于 2024-8-14 21:57
对于Am元素，Basis Set Exchange上的SARC-DKH2 就是ORCA程序内置的SARC-DKH-TZVP，可以这么验证：
写一 ...

感谢两位老师。那我在Gaussian中是直接用Basis Set Exchange上的SARC-DKH2还是用您上面提到的修改后的Am.gjf，如果是用Am.gjf，那应该怎么修改就可以用于Gaussian计算。

zhaoxiaoni

zjxitcc 发表于 2024-8-14 21:57
对于Am元素，Basis Set Exchange上的SARC-DKH2 就是ORCA程序内置的SARC-DKH-TZVP，可以这么验证：
写一 ...

老师，我在运行orca_2mkl Am -mkl这条命令的时候提示
Reading the input file Am.gbw      ...

WARNING!!!!!!!
Your GBWFile is either corrupt or from a different ORCA version!
Please be VERY careful with your calculation results!!!

Data_size, sizeof( TOrcaInfo ): 651728, 491312
[file orca_tools/qcmem.cpp, line 870]:  OUT OF MEMORY ERROR!

[file orca_tools/qcmem.cpp, line 870]:  OUT OF MEMORY ERROR!
我用的ORCA版本是ORCA-5.0.3, 看log文件是正常结束任务了的。

snljty2

zhaoxiaoni 发表于 2024-8-15 11:55
老师，我在运行orca_2mkl Am -mkl这条命令的时候提示
Reading the input file Am.gbw      ...

先确保你计算用的orca和转波函数用的orca_2mkl是来自同一个版本的看看。

zhaoxiaoni

snljty2 发表于 2024-8-15 14:46
先确保你计算用的orca和转波函数用的orca_2mkl是来自同一个版本的看看。

我是直接就在当前计算的文件夹下运行的orca_2mkl这个命令

snljty2

zhaoxiaoni 发表于 2024-8-15 16:55
我是直接就在当前计算的文件夹下运行的orca_2mkl这个命令

跟在哪个目录没关系，提交的脚本里面加入

1. dirname $(\which orca)

复制代码

命令，然后跑orca_2mkl之前看

1. dirname $(\which orca_2mkl)

复制代码

，看两者是否一致。

zjxitcc

zhaoxiaoni 发表于 2024-8-15 11:55
老师，我在运行orca_2mkl Am -mkl这条命令的时候提示
Reading the input file Am.gbw      ...

这是因为你混用了不同版本的ORCA所致，这与混用不同版本Gaussian类似，都是十分不好的习惯，解决办法在18L已经有人说了。要确保orca_2mkl和orca是同一个版本的，比如都是5.0.3。此贴2L和4L我也给出了详细解答http://bbs.keinsci.com/thread-47632-1-1.html

zjxitcc

zhaoxiaoni 发表于 2024-8-15 10:58
感谢两位老师。那我在Gaussian中是直接用Basis Set Exchange上的SARC-DKH2还是用您上面提到的修改后的Am. ...

前一个问题：都可以。
后一个问题：取决于你要做什么计算。产生的Am.gjf文件中写好了UHF的关键词，当然你还可以修改核数和内存等。这里举两个例子
（1）UHF计算，gjf文件前几行示例

1. %chk=Am.chk
   
2. %nprocshared=64
   
3. %mem=180GB
   
4. #p UHF/gen int(nobasistransform,DKH2) nosymm guess=cards stable=opt
   
5. ...(剩余内容不动)

复制代码

相当于用ORCA的UHF初猜，使用Gaussian做UHF计算，开启了DKH2哈密顿，并检验波函数稳定性。这个单原子的计算结果是E(UHF) = -30336.1908255 a.u.

（2）UB3LYP计算，gjf文件前几行示例

1. %chk=Am_GAU.chk
   
2. %nprocshared=64
   
3. %mem=180GB
   
4. #p UB3LYP/gen int(nobasistransform,DKH2) nosymm stable=opt
   
5. 
   
6. generated by utility mkl2gjf in MOKIT
   
7. 
   
8. 0 8
   
9. Am         0.00000000         0.00000000         0.00000000
   
10. 
    
11. 1 0
    
12. S   9  1.00
    
13. ... (基组数据)
    
14. (空行)
    
15. (空行)

复制代码

删掉底下轨道数据，这里我们做的是UB3LYP计算，用不到UHF初猜。