基态复合物优化不收敛

dingliping1108 · 发表于 Post on 2016-7-14 20:02:44

请教各位老师优化一个70多个原子的复合物，用SVP优化收敛了，命令如下! B3LYP RIJCOSX def2-SVP def2-SVP/J TightSCF Opt freq
%geom Calc_Hess true
end

但是将SVP优化后的结构换大基组TZVP继续优化一直不收敛

用TZVP优化时，刚开始给%geom Calc_Hess true这个命令优化，到最大步数后读之前的hess用下面这个命令
! B3LYP RIJCOSX def2-TZVP def2-TZVP/J TightSCF Opt freq
%geom
InHess Read
InHessName "20160701.hess"
end一直不收敛，想人为动一下内部的原子结构，动完之后可以读之前的hess吗？还是需要重新计算一下hess？因为计算hess要好久
在基态结构优化的时候，需要给这个命令%geom Calc_Hess true精确计算hesｓ吗？？
有什么好一点的方法让他收敛呢？
谢谢各位老师

sobereva · 发表于 Post on 2016-7-14 21:36:41

def2-SVP优化精度已经足够了（特别是对于多达70个原子的情况），通常没必要def2-TZVP优化

另外，你用杂化泛函，辅助基组应当用/JK。

碰见不好收敛的情况，适当增加积分格点质量再试。至于精确计算Hessian之类的之后再说。

plus · 发表于 Post on 2016-7-15 11:02:49

优化结构可以用def2-tzvp(-f)，效率会高点。
你看下是为什么不收敛，是震荡吗？

dingliping1108 · 发表于 Post on 2016-7-15 11:03:44

sobereva 发表于 2016-7-14 21:36
def2-SVP优化精度已经足够了（特别是对于多达70个原子的情况），通常没必要def2-TZVP优化

另外，你用杂 ...

请教sob老师，您的意思是不是先增加积分格点质量，实在不行再精确计算hess？
%method
Grid 0 # product grid
   1 # Lebedev 50 (not recommended)
   2 # Lebedev 110 points (default for SCF iterations)
   3 # Lebedev 194 points (more accurate)
   4 # Lebedev 302 points (default for FinalGrid)
   5 # Lebedev 434 points (large)
   6 # Lebedev 590 points (larger)
   7 # Lebedev 770 points (very large)
end
是不是用上面那个命令增加积分格点？从1加到7？

dingliping1108 · 发表于 Post on 2016-7-15 11:19:39

plus 发表于 2016-7-15 11:02
优化结构可以用def2-tzvp(-f)，效率会高点。
你看下是为什么不收敛，是震荡吗？

怎么看是不是震荡啊？优化复合物的时候，它总是跑偏。两个反应物跑远了。

plus · 发表于 Post on 2016-7-15 11:24:32

dingliping1108 发表于 2016-7-15 11:19
怎么看是不是震荡啊？优化复合物的时候，它总是跑偏。两个反应物跑远了。

加上d3bj试试

sobereva · 发表于 Post on 2016-7-15 12:34:53

dingliping1108 发表于 2016-7-15 11:03
请教sob老师，您的意思是不是先增加积分格点质量，实在不行再精确计算hess？
%method
Grid 0 # product ...

Grid4且GRIDX4

liyuanhe211 · 发表于 Post on 2016-7-15 13:47:31

dingliping1108 发表于 2016-7-15 11:19
怎么看是不是震荡啊？优化复合物的时候，它总是跑偏。两个反应物跑远了。

提取每步优化后的能量，作图观察。

zhanfei · 发表于 Post on 2016-7-15 14:09:18

本帖最后由 zhanfei 于 2016-7-15 14:13 编辑

dingliping1108 发表于 2016-7-15 11:03
请教sob老师，您的意思是不是先增加积分格点质量，实在不行再精确计算hess？
%method
Grid 0 # product ...

这么多原子精确hessen算好慢，我在24核机器上试过，ORCA有算多少轮重新计算hessen的参数。另外感觉COSX似乎是默认GRIDX有问题,至少我有时候关闭密度拟合NORI或者大幅提高密度拟合格点GRIDX9时候几何收敛容易很多，ORCA自洽和最终计算格点精度不同GRID4+FinalGRID6这样同步提高可能更好

dingliping1108 · 发表于 Post on 2016-7-15 15:32:31

本帖最后由 dingliping1108 于 2016-7-15 15:35 编辑

sobereva 发表于 2016-7-15 12:34
Grid4且GRIDX4

是要 Grid4和GRIDX4一起写吗？再请教sob老师，为什么杂化泛函，辅助基组应当用/JK？因为之前我一直用的是/J辅助基组，如果要重新计算的话，每一个都需要重新计算

dingliping1108 · 发表于 Post on 2016-7-15 15:33:02

plus 发表于 2016-7-15 11:24
加上d3bj试试

您好请问b3bj是什么意思啊？

dingliping1108 · 发表于 Post on 2016-7-15 15:37:13

liyuanhe211 发表于 2016-7-15 13:47
提取每步优化后的能量，作图观察。

每步计算后的单点能吗？是不是如果能量是震荡的，复合物就是震荡的？

liyuanhe211 · 发表于 Post on 2016-7-15 16:48:03

dingliping1108 发表于 2016-7-15 15:37
每步计算后的单点能吗？是不是如果能量是震荡的，复合物就是震荡的？

一般说震荡就是说能量曲线震荡。

sobereva · 发表于 Post on 2016-7-15 19:28:56

dingliping1108 发表于 2016-7-15 15:33
您好请问b3bj是什么意思啊？

DFT-D色散校正的使用
http://sobereva.com/210

如果体系跟色散作用没太大联系的话，没必要加d3bj，否则反倒还有可能略微有碍优化收敛

sobereva · 发表于 Post on 2016-7-15 19:30:01

dingliping1108 发表于 2016-7-15 15:32
是要 Grid4和GRIDX4一起写吗？再请教sob老师，为什 ...

一起写
/JK可以得到更好数值精度，收敛也可能更容易。你先增加格点精度，不行再尝试/JK

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

[ORCA] 基态复合物优化不收敛

浏览过的版块