ORCA 第二步计算 M06-2X 单点计算错误

liyuanhe211

在ORCA中，希望优化某结构，随后使用M06-2X计算一个结构的单点能。

计算路径：

1. %pal nprocs 8 end
   
2. %maxcore 1500
   
3. 
   
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. *
   
9. 
   
10. %pal nprocs 3 end
    
11. %maxcore 1500
    
12. %scf SOSCFStart 0.00033 end
    
13. 
    
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
   
6. 
   
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
   
3. 
   
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
   
9. 
   
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
   
6. 
   
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
   
6. 
   
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
   
6. 
   
7. 
   
8. INFO   : the flag for use of LIBINT has been found!
   
9. 
   
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. ================================================================================
    
24. 
    
25.                        ****************************
    
26.                        * Single Point Calculation *
    
27.                        ****************************
    
28. 
    
29. ---------------------------------
    
30. CARTESIAN COORDINATES (ANGSTROEM)
    
31. ---------------------------------
    
32.   C      1.395014    2.320874   -0.205175
    
33. ......

复制代码

有没有办法解决，使其连续计算（除了写一个.sh分成两个输入文件）？

sobereva

算是bug，不辨是非，一看输入文件里又有freq又有meta-GGA，于是立马报错说没有解析Hessian
还是分成两个输入文件吧。一个任务一个输入文件，在管理计算数据上也比较清楚。

liyuanhe211

sobereva 发表于 2015-10-27 15:52
算是bug，不辨是非，一看输入文件里又有freq又有meta-GGA，于是立马报错说没有解析Hessian
还是分成两个输 ...

把Bug报给ORCAforum了
谢谢Sob