GAMESS算能量分解对锕系元素换哪些基组计算？

xyz0101 · 发表于 Post on 2023-9-21 22:02:16

在用GAMESS算能量分解，体系中含有锕系元素U，选用SBKJC基组，不支持86号之后的元素，该换那些基组计算呢？
报错如下： ERROR. SBKJC BASES DO NOT EXTEND PAST RN.
EXECUTION OF GAMESS TERMINATED -ABNORMALLY- AT

zjxitcc · 发表于 Post on 2023-9-21 22:18:19

本帖最后由 zjxitcc 于 2023-9-22 10:56 编辑

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

%mem=40GB
%nprocshared=20
#p UPBE1PBE/genecp guess(fragment=2) em=GD3
{gks}
0 2 0 2 0 1
U(fragment=1) 0.0 0.0 0.0
H(fragment=2) 0.0 0.0 2.0
...
H C O 0
def2TZVP
****
U 0
SDD
****
U 0
SDD

复制代码

提交给frag_guess_wfn小程序计算

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

复制代码

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

sobereva · 发表于 Post on 2023-9-21 22:10:52

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

sobereva · 发表于 Post on 2023-9-21 22:08:05

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[GAMESS-US] GAMESS算能量分解对锕系元素换哪些基组计算？