计算化学公社

标题: opt freq [打印本页]

作者
Author: 小菜鸟求关照 时间: 2014-11-11 15:53
标题: opt freq
请问：我用的相同的基组和方法，初始构型一样，一种是算opt 一种是freq 为什么两种结果优化得到的构型不一样呢只算opt的出现4个yes
算opt和freq的只前两个出现yes 这是为什么呢求指教？

作者
Author: 378024572 时间: 2014-11-11 19:31
为什么opt和freq不一起用呢？

作者
Author: sobereva 时间: 2014-11-11 21:52
如果你是指第一种情况只写opt，第二种情况写opt freq，那么第二种情况最后只有两个YES是因为在opt那一步用的Hessian不是精确的，而freq的时候则用的Hessian总是精确的。写opt(calcall) freq使得opt时也用精确Hessian才能绝对确保freq最后显示的也是四个yes，当然这样在优化过程中成本很高，故也可以考虑不用calcall而把opt的收敛限设成tight，这样使得优化出的结构更精准，也往往能使freq之后NO出现的可能性更少。

作者
Author: lyfchem 时间: 2014-11-12 09:18
两个YES是因为在opt那一步用的Hessian不是精确的？那为什么往往opt+freq计算的倒数第二圈是4个yes呢？

作者
Author: 小菜鸟求关照 时间: 2014-11-12 09:21

378024572 发表于 2014-11-11 19:31
为什么opt和freq不一起用呢？

因为一起用的时候freq后面出现两个no我想看一下是不是结构的问题

作者
Author: 小菜鸟求关照 时间: 2014-11-12 09:22

sobereva 发表于 2014-11-11 21:52
如果你是指第一种情况只写opt，第二种情况写opt freq，那么第二种情况最后只有两个YES是因为在opt那一步用 ...

opt的收敛限设成tight 是写成这样的关键词吗：opt（tight）

作者
Author: 小菜鸟求关照 时间: 2014-11-12 09:22

lyfchem 发表于 2014-11-12 09:18
两个YES是因为在opt那一步用的Hessian不是精确的？那为什么往往opt+freq计算的倒数第二圈是4个yes呢？

我那种也是这样情况很想知道为什么

作者
Author: sobereva 时间: 2014-11-12 09:23

lyfchem 发表于 2014-11-12 09:18
两个YES是因为在opt那一步用的Hessian不是精确的？那为什么往往opt+freq计算的倒数第二圈是4个yes呢？

你所说的倒数第二圈是opt最后一步。正是因为这里全yes了所以优化才收敛了。这步判断是否yes是根据优化时用的近似Hessian判断的，而最后一次是freq之后用精确Hessian判断的，结论不同正因为Hessian用的不同。

作者
Author: lyfchem 时间: 2014-11-12 09:27
哦如果用opt(calcall)计算相当于opt也是用精确Hessian 但是成本太高没必要这样是吗老师

作者
Author: sobereva 时间: 2014-11-12 09:35

lyfchem 发表于 2014-11-12 09:27
哦如果用opt(calcall)计算相当于opt也是用精确Hessian 但是成本太高没必要这样是吗老师

是的。
实际上，就算最后freq出来不是都yes，但NO的项里离收敛限偏差不太远，也没关系，结构是可用的。

作者
Author: lyfchem 时间: 2014-11-12 09:39
是的这样可以减少很多计算时间

作者
Author: 小菜鸟求关照 时间: 2014-11-17 09:35

sobereva 发表于 2014-11-12 09:23
你所说的倒数第二圈是opt最后一步。正是因为这里全yes了所以优化才收敛了。这步判断是否yes是根据优化时 ...

我发现加上opt（tight）以后就没有freq后面的converged了，搜索的converged直接就是complete上面那一个，这是为什么呢？这样计算的结果准确吗？

作者
Author: sobereva 时间: 2014-11-17 09:53

小菜鸟求关照发表于 2014-11-17 09:35
我发现加上opt（tight）以后就没有freq后面的converged了，搜索的converged直接就是complete上面那一个， ...

光从你的叙述上我弄不清楚。
freq过程之后肯定会判断一次YES/NO，这和写不写tight没关系。

作者
Author: 小菜鸟求关照 时间: 2014-11-17 10:23

sobereva 发表于 2014-11-17 09:53
光从你的叙述上我弄不清楚。
freq过程之后肯定会判断一次YES/NO，这和写不写tight没关系。

..skipping...
      Item             Value    Threshold  Converged?
Maximum Force          0.000000    0.000015    YES
RMS    Force          0.000000    0.000010    YES
Maximum Displacement    0.000007    0.000060    YES
RMS    Displacement    0.000001    0.000040    YES
Predicted change in Energy=-9.705411D-14
Optimization completed.

这是最后一个converged 之前这个不是在倒数第二个吗？为什么现在变成最后一个了呢

作者
Author: sobereva 时间: 2014-11-17 10:57
直接把输出文件压缩一下发上来（如果压缩后>5MB，用网盘）

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