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

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

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

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

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

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

不要写def2-svp/C def2-svp/J，应当写def2-SVP/JK
你也看看轨迹，考察具体结构是怎么变化的

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   

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

关键词怎么写的？

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

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

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

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

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

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

可以这么说。
ORCA的优化不像高斯那么稳健