计算化学公社

标题: 求助：锕系离子水解结构收敛不合理 [打印本页]

作者
Author: zyaks.liu 时间: 2025-7-12 15:11
标题: 求助：锕系离子水解结构收敛不合理
本帖最后由 zyaks.liu 于 2025-7-12 15:15 编辑

老师们好，我在计算锕系离子水解结构时候结构经常跑散至不合理，做气相优化得到结构，使用smd再优化结构跑散。寻求大佬建议。
(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 1)
图一是气相优化结构，图二是基于气相优化的初始结构，图三是利用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也多余？
祝各位老师身体健康，工作顺利

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

作者
Author: zyaks.liu 时间: 2025-7-12 19:54

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

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

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

作者
Author: wzkchem5 时间: 2025-7-13 10:58

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

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

作者
Author: 我本是个娃娃 时间: 2025-7-13 20:50
Am离子的电荷是多少？如果是+3，那么要注意体系的multiplicity。

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