求助：锕系离子水解结构收敛不合理

zyaks.liu

老师们好，我在计算锕系离子水解结构时候结构经常跑散至不合理，做气相优化得到结构，使用smd再优化结构跑散。寻求大佬建议。

图一是气相优化结构，图二是基于气相优化的初始结构，图三是利用guess=mix初猜的结构是第二步结果，拉开了一个水分子，图四是经过93steps优化后得到的变形结构
开头文件：气相下：
%nprocshared=28
%mem=56GB
%chk=Am3+8h2o-pbe0.chk
#p opt freq PBE1PBE/genecp nosymm int=noxctest scf=xqc guess=mix
溶剂化效应下：
%nprocshared=28
%mem=56GB
%chk=Am3+8h2o-pbe0.chk
#p opt freq PBE1PBE/genecp scrf=(smd,solvent=water) nosymm int=noxctest scf=xqc guess=mix
所有体系均能量收敛
是不是不应该用guess=mix？ 同时计算时间特别长，要好几百个核时，是不是scf=xqc也多余？
祝各位老师身体健康，工作顺利

wzkchem5

这么重的元素用赝势不可靠吧？应该用sf-X2C描述相对论效应

zyaks.liu

wzkchem5 发表于 2025-7-12 18:48
这么重的元素用赝势不可靠吧？应该用sf-X2C描述相对论效应

老师使用的是ECP60MWB赝势，外带参数
H-Komponente
1
2 1.00000000 0.00000000
S-H
1
2 16.77194770 541.53853851
P-H
1
2 9.78355277 191.08457137
D-H
1
2 9.03422660 146.40258923
F-H
1
2 7.12184827 60.39433314
G-H
1
2 13.41832253 -62.62618508

sobereva

没有完整输入文件，没法判断输入文件格式是否合理
guess=mix完全多余，别随便乱用。没事甭用xqc
参考下文的讨论
使用Gaussian做镧系金属配合物的量子化学计算
http://sobereva.com/581（http://bbs.keinsci.com/thread-21274-1-1.html）

问“计算时间特别长，要好几百个核时”这种问题前应当先自行查看收敛情况，弄清楚就是计算慢还是难收敛导致的

wzkchem5

zyaks.liu 发表于 2025-7-12 19:54
老师使用的是ECP60MWB赝势，外带参数
H-Komponente
1

这个赝势已知可以准确描述水合锕系离子吗？
你的体系用显式相对论方法真的算不动吗？如果算得动，那么即使赝势的精度可以，也最好用sf-X2C算。
如果我是审稿人，必问这两个问题

我本是个娃娃

Am离子的电荷是多少？如果是+3，那么要注意体系的multiplicity。