计算化学公社

标题: 求助：ORCA两种格式输入相同基组计算结果不同 [打印本页]

---

作者
Author: chands    时间: 2022-1-20 21:20
标题: 求助：ORCA两种格式输入相同基组计算结果不同


各位老师：
       我用ORCA 5.0.1在CCSD(T)/aug-cc-pVQZ-PP级别下计算Br2的单点能，用两种输入格式：
第一种

1. !  CCSD(T) aug-cc-pVQZ-PP tightSCF noautostart miniprint nopop
   
2. %maxcore  1000
   
3. %mdci
   
4.     maxiter 100
   
5. end
   
6. 
   
7. * xyz   0   1
   
8. Br                 0.00000000    0.00000000    1.14214800
   
9. Br                 0.00000000    0.00000000   -1.13885200
   
10. *
    

复制代码

第二种

1. !  CCSD(T)  tightSCF noautostart miniprint nopop
   
2. %maxcore  1000
   
3. %basis
   
4.     newGTO Br "aug-cc-pVQZ-PP" end
   
5. end
   
6. %mdci
   
7.     maxiter 100
   
8. end
   
9. 
   
10. * xyz   0   1
    
11. Br                 0.00000000    0.00000000    1.14214800
    
12. Br                 0.00000000    0.00000000   -1.13885200
    
13. *
    

复制代码

   结果第一种的单点能是-831.4347 Hartree, 第二种的单点能是-2391.8318 Hartree。请教
     （1）为什么结果差这么大？是不是第一种ORCA默认开了什么选项，第二种ORCA没有开默认选项？
     （2） 能够加一些关键词使得两种输入方式结果相同吗？
    我用newGTO命令是因为我的体系含有Br，想对Br用赝势基组，对其他原子用普通的基组，计算结果很奇怪，纠结了半天发现可能是Br的问题，就用Br2做个测试。

---

作者
Author: hebrewsnabla    时间: 2022-1-20 23:00
第一种是正确写法，第二种只指定了赝势基组，没有指定赝势。

---

作者
Author: sobereva    时间: 2022-1-20 23:02
第二个得定义赝势

---

作者
Author: chands    时间: 2022-1-21 02:18


谢谢两位老师指教，我加上赝势

1. !  CCSD(T)  tightSCF noautostart miniprint nopop
   
2. %maxcore  1000
   
3. %basis
   
4.     newGTO Br "aug-cc-pVQZ-PP" end
   
5.     newECP Br "SK-MCDHF-RSC" end
   
6. end
   
7. %mdci
   
8.     maxiter 100
   
9. end
   
10. 
    
11. * xyz   0   1
    
12. Br                 0.00000000    0.00000000    1.14214800
    
13. Br                 0.00000000    0.00000000   -1.13885200
    
14. *

复制代码

现在两种输入法结果一样了。

---

欢迎光临 计算化学公社 (http://bbs.keinsci.com/)

Powered by Discuz! X3.3