ORCA 的 %base 不能过长

liyuanhe211

解决了ORCA中的一个奇怪的问题，提醒一下，也当做笔记。

用ORCA做双杂化计算，路径如下，计算正常结束但结果不合理。

1. !RI-PWPB95 D3 RIJCOSX def2-TZVP def2-TZVP/J def2-TZVP/C grid4 tightSCF CDouble GCP(DFT/TZ)

复制代码

检查后发现部分文件的gCP计算没有运行、gCP校正能为0，但不报错。

1. 
   
2. -------------------------   ----------------
   
3. Dispersion correction           -0.028678878
   
4. -------------------------   ----------------
   
5. 
   
6. input error
   
7. ERROR: Number of atoms in ORCA and otool_gcp are different:
   
8. ORCA: 52  gCP: 0
   
9. NOTIFICATION: Different basis set in ORCA and otool_gcp:
   
10. ORCA: 862  gCP: 133110635
    
11. If you are NOT using ECPs, check your basis set inputs!
    
12. ------------------   -----------------
    
13. gCP correction             0.000000000
    
14. ------------------   -----------------
    
15. 
    
16. -------------------------   --------------------
    
17. FINAL SINGLE POINT ENERGY      -930.367301940789

复制代码

1. Timings for individual modules:
   
2. 
   
3. Sum of individual times         ...     7964.166 sec (= 132.736 min)
   
4. GTO integral calculation        ...        8.221 sec (=   0.137 min)   0.1 %
   
5. SCF iterations                  ...     2360.074 sec (=  39.335 min)  29.6 %
   
6. MP2 module                      ...     5595.871 sec (=  93.265 min)  70.3 %
   
7.                              ****ORCA TERMINATED NORMALLY****
   
8. TOTAL RUN TIME: 0 days 2 hours 12 minutes 49 seconds 29 msec
   

复制代码

查询、调整关键词许久，才发现与其无关，如下图对比，只要将 %base 名称改短即可正常计算。应该是ORCA程序内部字符串长度的某些限制、导致 otool_gcp 子程序找不到中间文件了。

几个教训：

• ORCA的%base名称、输入文件名不要太长，出现奇怪结果先行简化。（目前试过的最长%base、且计算正常的文件为 64 字符）
• ORCA的中间程序出错可能也不会报错、而（看似）正常的输出结果，注意检查
• 自己写辅助程序的时候不要生成过长的文件名→_→
  

amani73986

又是Python写的辅助程序吗？能不能请您分享一下您在实际工作中，利用自己写程序提高效率的经验呢？非常感谢！
最近也经常做一些重复性的工作（批次构型的优化、单点以及后续波函数分析等），经常在VMD Gaussian ORCA Vasp 还有作图软件之间折腾，学习了python，想像您一样写一些脚本或者GUI帮助自己完成这些繁琐的工作（监视程序运行，处理结果并传递给下一个软件等）。谢谢。

liyuanhe211

amani73986 发表于 2016-11-21 12:47
又是Python写的辅助程序吗？能不能请您分享一下您在实际工作中，利用自己写程序提高效率的经验呢？非常感谢 ...

经验指什么？容易自动化、需要重复数遍以上的工作全用程序实现→_→想实现什么就写什么呗。
另外做一个小lib，用来读各种输入、输出文件，留好接口，做好可扩展，想要新功能往上装。
GUI用的PyQt5

我本是个娃娃

李大博士又出实用的技巧了，撒花！

amani73986

liyuanhe211 发表于 2016-11-21 13:46
经验指什么？容易自动化、需要重复数遍以上的工作全用程序实现→_→想实现什么就写什么呗。
另外做一个 ...

谢谢你的回复，我先学着自己写，遇到具体问题再请教！