orca_loc模块报错

任嘉人amber

各位老师好，想问一下这个错误是什么原因呢，重复的是CASSCF tutorial的例子
input是内容是
PPD_cas66.gbw # input orbitals
LOC.gbw # output orbitals
0 # operator:alpha here (CASSCF), 1 for beta
0 # orbital window: first orbital to be localized e.g. first active
25 # orbital window: last orbital to be localized e.g. last active
50 # maximum number of iterations
1 # localization method: 1=PIPEK-MEZEY

任嘉人amber

shenzp 发表于 2023-10-31 13:29
按照ORCA手册orca_loc部分的说明，input文件应该是这样的顺序：

好的好的，谢谢老师

shenzp

按照ORCA手册orca_loc部分的说明，input文件应该是这样的顺序：

1. input.gbw           // input orbitals
   
2. output.gbw        // output orbitals
   
3. 10                       // orbital window: first orbital to be localized e.g. first active
   
4. 15                       // orbital window: last orbital to be localized e.g. last active
   
5. 0                         // localization method: 1=PM,2=FB,3=IAOIBO,4=IAOBOYS,5=NEWBOYS,6=AHFB
   
6. 0                         // operator: 0 for alpha, 1 for beta。对于CASSCF/ROHF/RHF，应该用0。
   
7. 128                     // maximum number of iterations
   
8. 1e-6                   // convergence tolerance of the localization functional value
   
9. 0.0                     // relative convergence tolerance of the localization functional value
   
10. 0.95                   // printing thresh to call an orbital strongly localized
    
11. 0.85                   // printing thresh to call an orbital bond-like
    
12. 2                       // printlevel
    
13. 1                       // use Cholesky Decomposition (0=false, 1=true)
    
14. 1                       // randomize seed for localization (0=false, 1=true)
    

复制代码

后边几个选项可以不写。报错是因为第五行应该是局域化方法，但是没有25这个选项。CASSCF tutorial好像主要是根据ORCA 4写的，有可能是ORCA版本更新后输入文件的顺序变了

另外，个人感觉Multiwfn的轨道局域化功能更方便，可以看一下Multiwfn的轨道定域化功能的使用以及与NBO、AdNDP分析的对比。