Gaussian对Fe(III)长链有机酸配合物结构优化相关问题请教

wzm12天空

本人在做分子1：硬脂酸铁（Fe(III)(OOC18H35)3）和分子2：HFe(II)(OOC18H35)2分子的相关计算。最开始做结构优化时计算级别为PBE0/6-311G* em=gd3bj，两个分子都能收敛。但接下来想计算同分子配体C自由基时对第一个分子相同级别下C1，C2都能正常收敛，C3以后就一直有-0.8左右的虚频，试了sob帖子：http://sobereva.com/278(Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法)中的除calcall（计算非常耗时）的其余方法都无法消除虚频。所以尝试换了6-311G**基组做优化，消除了虚频。由于我下一步要计算单点能，所以就把分子1和2也在6-311G**级别下重新优化。但是出现了问题：分子1可以正常收敛，但是分子2始终无法正常收敛（计算过程中查看log文件发现位移始终里比收敛限高两个数量级），期间试过GDIIS和maxstep=3/5,notrust关键词帮助收敛，由于离收敛太远没有尝试calcall相关关键词。
请教：1. 为什么6-311G*能够收敛，但给H加极化（6-311G**）却不能收敛？有什么比较好的办法帮助收敛吗？
         2. 结构优化用不同的基组（6-311G*和6-311G**），单点用统一基组做计算可行吗？对结果影响大吗？
附件fehst3_opt.gjf和fehst3_opt.log为正在计算的分子2的相关文件。

sobereva

1 纯属巧合。只要用恰当的关键词帮助收敛，二者都必定能收敛。能解决问题的关键词绝对不会超过下文所列范畴，反复认真一个字一个字阅读
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

“离收敛太远没有尝试calcall相关关键词” 这句话完全没有逻辑。怎么会以为离收敛限较远时calcall、recalc就不起作用？

2 原则上不行

wzm12天空

sobereva 发表于 2023-7-13 12:39
1 纯属巧合。只要用恰当的关键词帮助收敛，二者都必定能收敛。能解决问题的关键词绝对不会超过下文所列范畴 ...

谢谢sob老师
那像图中这样的分子，您给推荐几个有助于收敛的关键词吧，谢谢！

sobereva

wzm12天空 发表于 2023-7-13 13:05
谢谢sob老师
那像图中这样的分子，您给推荐几个有助于收敛的关键词吧，谢谢！

把结构尽可能摆得满足对称性，让程序判断出点群
之后考虑maxstep=3,notrust,recalc=3

wzm12天空

sobereva 发表于 2023-7-13 13:13
把结构尽可能摆得满足对称性，让程序判断出点群
之后考虑maxstep=3,notrust,recalc=3

谢谢sob老师，我试试

wzm12天空

sobereva 发表于 2023-7-13 12:39
1 纯属巧合。只要用恰当的关键词帮助收敛，二者都必定能收敛。能解决问题的关键词绝对不会超过下文所列范畴 ...

老师，我这个结构优化用PBE0/6-311G**; 单点能计算用PBE0/def2-TZVP级别可以吗？

sobereva

wzm12天空 发表于 2023-7-13 14:42
老师，我这个结构优化用PBE0/6-311G**; 单点能计算用PBE0/def2-TZVP级别可以吗？

基组可以。
应当带上色散校正。PBE0对色散作用描述很差。考虑色散作用后，很长的烃链会倾向于聚拢

wzm12天空

sobereva 发表于 2023-7-14 10:50
基组可以。
应当带上色散校正。PBE0对色散作用描述很差。考虑色散作用后，很长的烃链会倾向于聚拢

好的老师，谢谢。另外对于泛函，我查了下ωB97M-V计算能量是比较精确的，但是Gaussian没有内置。请问老师目前来说这个泛函可以通过手动设置参数实现Gaussian计算吗？我用的事G16.C01版本

sobereva

wzm12天空 发表于 2023-7-14 10:54
好的老师，谢谢。另外对于泛函，我查了下ωB97M-V计算能量是比较精确的，但是Gaussian没有内置。请问老师 ...

没法通过自定义方式使用

wzm12天空

sobereva 发表于 2023-7-14 11:18
没法通过自定义方式使用

明白了，谢谢老师的解答