ORCA中NumGrad不输出gradient文件的困惑

winnerwill

ORCA 5.0.2 版本，
发现仅仅在于NumGrad和EnGrad关键词的不同，前者没有输出gradient文件，后者有(.engrad文件)。
怎么控制NumGrad的文件输出呢？还是它的文件名是什么tmp文件被最后删除了？
找了手册没有找到线索。

关键词的唯一不同：

1. head -n 1 *inp
   
2. ==> molengrad.inp <==
   
3. ! B3LYP def2-SVP EnGrad
   
4. 
   
5. ==> mol.inp <==
   
6. ! B3LYP def2-SVP
   
7. 
   
8. ==> molnumgrad.inp <==
   
9. ! B3LYP def2-SVP NumGrad
   
10. 
    

复制代码

最后的输出文件：

1. ls -tr1
   
2. molengrad.inp
   
3. molengrad.densities
   
4. molengrad.gbw
   
5. molengrad.engrad
   
6. molengrad_property.txt
   
7. molengrad.log
   
8. molengrad.molden
   
9. mol.inp
   
10. mol.densities
    
11. mol.gbw
    
12. mol_property.txt
    
13. mol.log
    
14. mol.molden
    
15. molnumgrad.inp
    
16. molnumgrad.densities
    
17. molnumgrad.gbw
    
18. molnumgrad_property.txt
    
19. molnumgrad.log
    
20. molnumgrad.molden
    

复制代码

wzkchem5

要算数值梯度，需要既写NumGrad又写EnGrad。
NumGrad的意思是“如果算梯度的话，那么算数值梯度”，而不是“算数值梯度”。这一点经常被人误解（尤其是和NumFreq的意思不统一），确实造成了一些不便，但是出于历史原因，程序后来也没有再修改。

winnerwill

wzkchem5 发表于 2022-10-28 22:36
要算数值梯度，需要既写NumGrad又写EnGrad。
NumGrad的意思是“如果算梯度的话，那么算数值梯度”，而不是 ...

还真是，谢谢大佬！
本以为NumGrad和EnGrad是互斥的，然后试过Energy + NumGrad 就不行。

winnerwill

wzkchem5 发表于 2022-10-28 22:36
要算数值梯度，需要既写NumGrad又写EnGrad。
NumGrad的意思是“如果算梯度的话，那么算数值梯度”，而不是 ...

这个传统跟NumFreq不一样吧？

wzkchem5

winnerwill 发表于 2022-10-28 16:10
这个传统跟NumFreq不一样吧？

没错，我说的就是这个意思，我也不知道为什么。不过现在很多用户用习惯了，也不好再改成和NumFreq一样了。