计算化学公社

标题: ORCA优化过程中能量变化 [打印本页]

作者
Author: yaochuang 时间: 2016-4-20 17:52
标题: ORCA优化过程中能量变化
最近用ORCA优化一个结构，最后把单点能提出来看了看，发现很奇怪，在能量较低的时候没有收敛，反而在能量升高一点的时候收敛了，这是不是就可以说明这个最终结构不是最优结构呢？
(, 下载次数 Times of downloads: 61)

作者
Author: sobereva 时间: 2016-4-20 18:13
可以这么说。
ORCA的优化不像高斯那么稳健

作者
Author: yaochuang 时间: 2016-4-20 20:10

sobereva 发表于 2016-4-20 18:13
可以这么说。
ORCA的优化不像高斯那么稳健

sob老师，那针对这种情况有没有好的办法解决呢？

作者
Author: sobereva 时间: 2016-4-20 20:13

yaochuang 发表于 2016-4-20 20:10
sob老师，那针对这种情况有没有好的办法解决呢？

取最低能量结构，减小步长上限再优化看看。

作者
Author: yaochuang 时间: 2016-4-21 10:59

sobereva 发表于 2016-4-20 20:13
取最低能量结构，减小步长上限再优化看看。

谢谢sob老师，我尝试了，将maxstep 设置的更小，并且用了更加严格的收敛限，但是最终的结果还是一样，能量会升高然后收敛，而且我比较了两种方法得到的结构，基本完全一样（求了一下RMSD为零）。而且有时候orca在判断收敛时会忽略掉能量的变化项，所以会不会是orca软件本身在定义收敛时对能量的要求就没有那么严格呢？

作者
Author: sobereva 时间: 2016-4-21 17:38

yaochuang 发表于 2016-4-21 10:59
谢谢sob老师，我尝试了，将maxstep 设置的更小，并且用了更加严格的收敛限，但是最终的结果还是一样，能 ...

关键词怎么写的？

作者
Author: yaochuang 时间: 2016-4-21 18:49

sobereva 发表于 2016-4-21 17:38
关键词怎么写的？

关键词是这样的：
!b3lyp/g def2-svp rijcosx def2-svp/C def2-svp/J tightopt
!pal6
%maxcore 1000
%geom
maxstep 0.1
end

最初用的是：
!b3lyp/g def2-svp rijcosx def2-svp/C def2-svp/J opt
!pal6
%maxcore 1000

作者
Author: sobereva 时间: 2016-4-21 19:28
不要写def2-svp/C def2-svp/J，应当写def2-SVP/JK
你也看看轨迹，考察具体结构是怎么变化的

作者
Author: yaochuang 时间: 2016-4-22 14:54
本帖最后由 yaochuang 于 2016-4-22 15:48 编辑

sobereva 发表于 2016-4-21 19:28
不要写def2-svp/C def2-svp/J，应当写def2-SVP/JK
你也看看轨迹，考察具体结构是怎么变化的

sob老师，我用def2-SVP/JK代替了原来的辅助基组，结果还是一样。并测试了两个结构的RMSD结果是0.007 埃（这应该是由辅助基组改变造成的）。
从轨迹上看就是基团间的小幅的调整，没有看出可疑的地方。
对了，我还想问一下对于def2-svp/C def2-svp/J和def2-svp/jk这两种辅助基组有没有优劣比较呢？因为记得在哪儿好像看到过说分别指定Coulomb和correlation辅助基组会更好一些，不知道是不是这样？
。 (, 下载次数 Times of downloads: 69)

作者
Author: sobereva 时间: 2016-4-22 18:19

yaochuang 发表于 2016-4-22 14:54
sob老师，我用def2-SVP/JK代替了原来的辅助基组，结果还是一样。并测试了两个结构的RMSD结果是0.007 埃（ ...

那就在优化中用精确计算的Hessian
/C是给后HF计算用的，DFT计算无用。
/J是给纯泛函用的，杂化泛函应当用/JK。

作者
Author: yaochuang 时间: 2016-4-24 08:30

sobereva 发表于 2016-4-22 18:19
那就在优化中用精确计算的Hessian
/C是给后HF计算用的，DFT计算无用。
/J是给纯泛函用的，杂化泛函应当 ...

sob老师，我添加了每一步都计算Hess矩阵，但是结果还是这样：
(, 下载次数 Times of downloads: 45)

这是我的输入：
!b3lyp/g def2-svp rijcosx def2-svp/JK tightopt tightscf
!pal6
%maxcore 1000
%geom
maxstep 0.1
Recalc_Hess 1
end

作者
Author: sobereva 时间: 2016-4-24 09:28

yaochuang 发表于 2016-4-24 08:30
sob老师，我添加了每一步都计算Hess矩阵，但是结果还是这样：

那你索性还是先用高斯搞搞吧，要不然虽然ORCA更快，翻来覆去折腾时间也都浪费了，还不如直接用高斯

作者
Author: yaochuang 时间: 2016-4-24 10:53

sobereva 发表于 2016-4-24 09:28
那你索性还是先用高斯搞搞吧，要不然虽然ORCA更快，翻来覆去折腾时间也都浪费了，还不如直接用高斯

嗯，好的！谢谢sob老师耐心的解答。

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