计算化学公社

标题: ORCA 第二步计算 M06-2X 单点计算错误 [打印本页]
作者
Author: liyuanhe211    时间: 2015-10-27 15:17
标题: ORCA 第二步计算 M06-2X 单点计算错误
在ORCA中，希望优化某结构，随后使用M06-2X计算一个结构的单点能。

计算路径：

  1. %pal nprocs 8 end
  2. %maxcore 1500

  4. %scf SOSCFStart 0.00033 end
  5. ! opt freq B3LYP/G 6-31+G(d,p) D3 grid4 xyzfile
  6. * xyz 0   1
  7. [the structure]
  8. *

  10. %pal nprocs 3 end
  11. %maxcore 1500
  12. %scf SOSCFStart 0.00033 end

  14. ! M062X ma-def2-TZVP D3 Grid6
  15. * xyzfile 0 1 filename.xyz
复制代码

返回错误

  1. WARNING: Analytical Hessians are not yet implemented for meta-GGA functionals
  2.          Please use the numerical Hessian instead:
  3.          for frequencies  : !NumFreq
  4.          for optimisations: %geom Calc_Hess true; NumHess true end
  5.   ===> : Skipping actual calculation

  7. Error (ORCA_MAIN): ... aborting the run
  8. Error (ORCA_MAIN): ... aborting the run
复制代码

按照要求加入NumFreq（实际上并不想算Freq）

  1. ! NumFreq M062X ma-def2-TZVP D3 Grid6
复制代码

返回错误

  1. WARNING: Gradients needed for Numerical Frequencies
  2.   ===> : Setting RunTyp to EnGrad

  4. WARNING: Analytical Hessians are not yet implemented for meta-GGA functionals
  5.          Please use the numerical Hessian instead:
  6.          for frequencies  : !NumFreq
  7.          for optimisations: %geom Calc_Hess true; NumHess true end
  8.   ===> : Skipping actual calculation

  10. Error (ORCA_MAIN): ... aborting the run
  11. Error (ORCA_MAIN): ... aborting the run
复制代码

再加入 EnGrad，返回与第一次相同

  1. ! EnGrad NumFreq M062X ma-def2-TZVP D3 Grid6
复制代码
  1. WARNING: Analytical Hessians are not yet implemented for meta-GGA functionals
  2.          Please use the numerical Hessian instead:
  3.          for frequencies  : !NumFreq
  4.          for optimisations: %geom Calc_Hess true; NumHess true end
  5.   ===> : Skipping actual calculation

  7. Error (ORCA_MAIN): ... aborting the run
  8. Error (ORCA_MAIN): ... aborting the run
复制代码
而如果直接把输出的.xyz文件做输入，是可以运行的。
  1. #没有上一步
  2. %pal nprocs 3 end
  3. %maxcore 1500
  4. %base "Me2_34_Prod_step1"
  5. %scf SOSCFStart 0.00033 end

  7. ! M062x ma-def2-tzvp D3 Grid6
  8. * xyzfile 0 1 last_output.xyz
复制代码
输出：
  1. ================================================================================
  2.                                         WARNINGS
  3.                        Please study these warnings very carefully!
  4. ================================================================================
  5. Now building the actual basis set


  8. INFO   : the flag for use of LIBINT has been found!

  10. ================================================================================
  11.                                        INPUT FILE
  12. ================================================================================
  13. NAME = SP.inp
  14. |  1> %pal nprocs 3 end
  15. |  2> %maxcore 1500
  16. |  3> %base "Me2_34_Prod_step1"
  17. |  4> %scf SOSCFStart 0.00033 end
  18. |  5>
  19. |  6> ! M062x ma-def2-tzvp D3 Grid6
  20. |  7> * xyzfile 0   1 SP.xyz
  21. |  8>
  22. |  9>                          ****END OF INPUT****
  23. ================================================================================

  25.                        ****************************
  26.                        * Single Point Calculation *
  27.                        ****************************

  29. ---------------------------------
  30. CARTESIAN COORDINATES (ANGSTROEM)
  31. ---------------------------------
  32.   C      1.395014    2.320874   -0.205175
  33. ......
复制代码
有没有办法解决，使其连续计算（除了写一个.sh分成两个输入文件）？

作者
Author: sobereva    时间: 2015-10-27 15:52
算是bug，不辨是非，一看输入文件里又有freq又有meta-GGA，于是立马报错说没有解析Hessian
还是分成两个输入文件吧。一个任务一个输入文件，在管理计算数据上也比较清楚。
作者
Author: liyuanhe211    时间: 2015-10-27 17:00
sobereva 发表于 2015-10-27 15:52
算是bug，不辨是非，一看输入文件里又有freq又有meta-GGA，于是立马报错说没有解析Hessian
还是分成两个输 ...

把Bug报给ORCAforum了
谢谢Sob





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