计算化学公社

标题: ORCA 的 %base 不能过长 [打印本页]
作者
Author: liyuanhe211    时间: 2016-11-21 11:31
标题: ORCA 的 %base 不能过长
本帖最后由 liyuanhe211 于 2016-11-21 16:23 编辑

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

用ORCA做双杂化计算，路径如下，计算正常结束但结果不合理。
  1. !RI-PWPB95 D3 RIJCOSX def2-TZVP def2-TZVP/J def2-TZVP/C grid4 tightSCF CDouble GCP(DFT/TZ)
复制代码

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

  2. -------------------------   ----------------
  3. Dispersion correction           -0.028678878
  4. -------------------------   ----------------

  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. ------------------   -----------------

  16. -------------------------   --------------------
  17. FINAL SINGLE POINT ENERGY      -930.367301940789
复制代码
  1. Timings for individual modules:

  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 子程序找不到中间文件了。

(, 下载次数 Times of downloads: 22)



几个教训：








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

经验指什么？容易自动化、需要重复数遍以上的工作全用程序实现→_→想实现什么就写什么呗。
另外做一个小lib，用来读各种输入、输出文件，留好接口，做好可扩展，想要新功能往上装。
GUI用的PyQt5
作者
Author: 我本是个娃娃    时间: 2016-11-21 14:39
李大博士又出实用的技巧了，撒花！
作者
Author: amani73986    时间: 2016-11-22 08:21
liyuanhe211 发表于 2016-11-21 13:46
经验指什么？容易自动化、需要重复数遍以上的工作全用程序实现→_→想实现什么就写什么呗。
另外做一个 ...

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



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3