计算化学公社

标题: opt freq和opt然后freq有区别么？ [打印本页]

作者
Author: fankui1990129 时间: 2015-6-13 11:42
标题: opt freq和opt然后freq有区别么？
我前几天计算一个结构，先opt，在freq，总共用了10小时左右吧。今天我进行了同样的计算，只不过用的是opt freq，已经计算24小时了，还没有结束。请问这是为什么呢？

两次用的其他设置是一模一样的。

作者
Author: 小范范1989 时间: 2015-6-13 18:28
是不是 1+1 不等于2？

作者
Author: sobereva 时间: 2015-6-13 20:06
要仔细阅读输出文件，比较一下，观看收敛情况就知道了。光看时间什么也不知道。
如果方法、基组以及其它关键词，包括初始结构都完全一样，很可能由于某种数值巧合，opt freq的优化过程没顺利收敛，反复震荡。

作者
Author: fankui1990129 时间: 2015-6-14 00:28
本帖最后由 fankui1990129 于 2015-6-14 01:16 编辑

sobereva 发表于 2015-6-13 20:06
要仔细阅读输出文件，比较一下，观看收敛情况就知道了。光看时间什么也不知道。
如果方法、基组以及其它关 ...

原来的计算结果被别人不小心删除了，所以要重算。

sob老师，有没有这样的可能。当时优化的时候没有优化到最佳的结构，只是得到能量极小值，而不是最小值。这种情况下，在算Freq时有虚频。而合并计算的时候，一边查看是否有虚频，一边优化，直至找到最小值未知。

（刚刚自学这部分，不了解具体的情况，一些理论细节也不懂，趁着这个机会好好了解一些）

作者
Author: fankui1990129 时间: 2015-6-14 00:50

小范范1989 发表于 2015-6-13 18:28
是不是 1+1 不等于2？

不知道啊，我也觉得比较好奇

作者
Author: sobereva 时间: 2015-6-14 03:01

fankui1990129 发表于 2015-6-14 00:28
原来的计算结果被别人不小心删除了，所以要重算。

sob老师，有没有这样的可能。当时优化的时候没有优 ...

不是。
用opt优化，只能得到离初猜最近的极小点结构。如果这个极小点结构恰好就是最小点，那你得到的就是最小点结构。
无论是最小点，还是能量比它高的极小点结构，都是没有虚频的。

opt freq是由高斯自动拆成opt和freq两个任务依次单独去做（注freq的时候用了guess=read），和手动先做opt再做freq并无区别。

作者
Author: sobereva 时间: 2015-6-14 03:09

fankui1990129 发表于 2015-6-14 00:50
不知道啊，我也觉得比较好奇

计算机以整数来计算1+1总是等于2的。但对于浮点数的加法，累加顺序不同结果往往会不同，这也是并行运算结果无法精确重现的原因之一。感兴趣的话建议看看此文
数值误差对计算化学结果重现性的影响
http://sobereva.com/88

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