计算化学公社

标题: GAMESS算能量分解对锕系元素换哪些基组计算？ [打印本页]

---

作者
Author: xyz0101    时间: 2023-9-21 22:02
标题: GAMESS算能量分解对锕系元素换哪些基组计算？
在用GAMESS算能量分解，体系中含有锕系元素U，选用SBKJC基组，不支持86号之后的元素，该换那些基组计算呢？
报错如下： ERROR. SBKJC BASES DO NOT EXTEND PAST RN.
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT

---

作者
Author: sobereva    时间: 2023-9-21 22:08
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 “EDA-锕系元素” 改了，以后务必注意

---

作者
Author: sobereva    时间: 2023-9-21 22:10
用sobEDA省事得多，速度还更快。镧系通常用SDD即可，更好的def2-TZVP也可以
使用sobEDA和sobEDAw方法做非常准确、快速、方便、普适的能量分解分析
http://sobereva.com/685（http://bbs.keinsci.com/thread-39446-1-1.html）

---

作者
Author: zjxitcc    时间: 2023-9-21 22:18


可以对U元素用SDD。使用MOKIT的frag_guess_wfn小程序可自动调用Gaussian做计算，并生成GKS-EDA或LMO-EDA输入文件，内含写好的EDA关键词、坐标、基组和赝势数据，不用自己去下载SDD数据然后手动编辑GAMESS自定义基组；还内含高斯收敛的轨道，极大加速GAMESS的EDA计算收敛。例如我们写一个文件U.gjf

1. %mem=40GB
   
2. %nprocshared=20
   
3. #p UPBE1PBE/genecp guess(fragment=2) em=GD3
   
4. 
   
5. {gks}
   
6. 
   
7. 0 2 0 2 0 1
   
8. U(fragment=1)   0.0   0.0   0.0
   
9. H(fragment=2)   0.0   0.0   2.0
   
10. ...
    
11. 
    
12. H C O 0
    
13. def2TZVP
    
14. ****
    
15. U 0
    
16. SDD
    
17. ****
    
18. 
    
19. U 0
    
20. SDD
    
21. 
    

复制代码

提交给frag_guess_wfn小程序计算

1. frag_guess_wfn U.gjf >U.out 2>&1

复制代码

算完即生成U.inp，可直接提交给GAMESS或XEDA进行EDA计算。更多示例见https://jeanwsr.gitlab.io/mokit-doc-mdbook/chap5-3.html。

---

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

Powered by Discuz! X3.3