计算化学公社

标题: 基态复合物优化不收敛 [打印本页]

作者
Author: dingliping1108 时间: 2016-7-14 20:02
标题: 基态复合物优化不收敛
请教各位老师优化一个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ｓ吗？？
有什么好一点的方法让他收敛呢？
谢谢各位老师

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

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

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

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

作者
Author: dingliping1108 时间: 2016-7-15 11:03

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？

作者
Author: dingliping1108 时间: 2016-7-15 11:19

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

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

作者
Author: plus 时间: 2016-7-15 11:24

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

加上d3bj试试

作者
Author: sobereva 时间: 2016-7-15 12:34

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

Grid4且GRIDX4

作者
Author: liyuanhe211 时间: 2016-7-15 13:47

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

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

作者
Author: zhanfei 时间: 2016-7-15 14:09
本帖最后由 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这样同步提高可能更好

作者
Author: dingliping1108 时间: 2016-7-15 15:32
本帖最后由 dingliping1108 于 2016-7-15 15:35 编辑

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

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

作者
Author: dingliping1108 时间: 2016-7-15 15:33

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

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

作者
Author: dingliping1108 时间: 2016-7-15 15:37

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

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

作者
Author: liyuanhe211 时间: 2016-7-15 16:48

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

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

作者
Author: sobereva 时间: 2016-7-15 19:28

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

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

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

作者
Author: sobereva 时间: 2016-7-15 19:30

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

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

作者
Author: dingliping1108 时间: 2016-7-15 20:06

liyuanhe211 发表于 2016-7-15 16:48
一般说震荡就是说能量曲线震荡。

明白了谢谢

作者
Author: dingliping1108 时间: 2016-7-15 20:12

sobereva 发表于 2016-7-15 19:30
一起写
/JK可以得到更好数值精度，收敛也可能更容易。你先增加格点精度，不行再尝试/JK

好的谢谢sob老师

作者
Author: dingliping1108 时间: 2016-7-19 10:38

sobereva 发表于 2016-7-15 19:30
一起写
/JK可以得到更好数值精度，收敛也可能更容易。你先增加格点精度，不行再尝试/JK

再请教sob老师，如果过渡态优化，优化刚开始精确计算了一遍hess，后面优化了几十圈之后，人为的对优化后的结构稍作修改，这样修改后的结构继续进行过渡态优化可不可以调用之前已经算好的hess？？还是需要重新计算？

作者
Author: sobereva 时间: 2016-7-20 12:17

dingliping1108 发表于 2016-7-19 10:38
再请教sob老师，如果过渡态优化，优化刚开始精确计算了一遍hess，后面优化了几十圈之后，人为的对优化后 ...

需要重新算。
不同结构下Hessian明显不同

作者
Author: dingliping1108 时间: 2016-7-20 20:31
本帖最后由 dingliping1108 于 2016-7-20 20:36 编辑

sobereva 发表于 2016-7-20 12:17
需要重新算。
不同结构下Hessian明显不同

谢谢sob老师但是如果我没有手动改它的结构，优化达到最大步数，需要继续优化是不是就可以调用以前的hess？？

作者
Author: sobereva 时间: 2016-7-21 00:11

dingliping1108 发表于 2016-7-20 20:31
谢谢sob老师但是如果我没有手动改它的结构，优化达到最大步数，需要继续优化是不是就可以调用以前的hess ...

是

作者
Author: dingliping1108 时间: 2016-7-21 09:01

sobereva 发表于 2016-7-21 00:11
是

谢谢sob老师指导

作者
Author: Theor_Comp 时间: 2019-9-12 09:58

sobereva 发表于 2016-7-15 19:30
一起写
/JK可以得到更好数值精度，收敛也可能更容易。你先增加格点精度，不行再尝试/JK

这样的话，用杂化泛函做RIJCOSX时，辅助基组不是必须用/JK吗？

作者
Author: sobereva 时间: 2019-9-12 11:05

Theor_Comp 发表于 2019-9-12 09:58
这样的话，用杂化泛函做RIJCOSX时，辅助基组不是必须用/JK吗？

完全错误
/JK是给RIJK用的
RIJCOSX只需要/J辅助基组。老版本ORCA也可以用/JK代替/J用，但显然不是“必须”

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