计算化学公社

标题: 关于 G09 中 Cartesian Force 和Internal Force收敛标准的疑惑 [打印本页]

作者
Author: LGC020305 时间: 2016-7-11 03:36
标题: 关于 G09 中 Cartesian Force 和Internal Force收敛标准的疑惑
本帖最后由 LGC020305 于 2016-7-11 09:04 编辑

RT,有几个概念不太清楚，希望大家帮忙解惑。
（1）G09中opt（关键字opt，不加其它选项）优化过程的Maximum Force， RMS Force 判断标准默认的是基于redundant internal coordinates吗？
（2）多原子体系的Cartesian Force 和Internal Force是一样数值吗？（好像大小不一样？）
（3）下面一个例子：opt+freq 最后一步力已经收敛，Cartesian Force 和Internal Force并不一样，Cartesian Force是0.004344856似乎这个job又没有收敛？
Maximum Force：  0.000145    0.000450    YES
RMS    Force： 0.000026    0.000300    YES
Cartesian Forces:  Max 0.004344856 RMS    0.000776246
遇到4个YES，而Cartesian Force又相对比较大该怎么处理？

作者
Author: 冰释之川 时间: 2016-7-11 08:19
默认情况下， OPT 收敛基于内坐标的，四个YES就完事了。其实我们更关注的是接下去Freq做出来有没有虚频（对于几何优化极小值点）

作者
Author: LGC020305 时间: 2016-7-11 09:06

冰释之川发表于 2016-7-11 08:19
默认情况下， OPT 收敛基于内坐标的，四个YES就完事了。其实我们更关注的是接下去Freq做出来有没有虚频（对 ...

谢谢你的回复。遇到opt+freq 4个YES，而Cartesian Force又相对比较大该怎么处理？

作者
Author: sobereva 时间: 2016-7-11 09:45
注意看读取cartesian force的地方和最后判断成四个YES的地方有无可比性。比如你取的是上一个结构乃至更早结构下的cartesian force，肯定比最后给出的冗余内坐标的力要大很多。
总的来说你不用管cartesian force的输出，就看最后对收敛判断的地方的数值就行了。

作者
Author: LGC020305 时间: 2016-7-11 10:48

sobereva 发表于 2016-7-11 09:45
注意看读取cartesian force的地方和最后判断成四个YES的地方有无可比性。比如你取的是上一个结构乃至更早结 ...

多谢。好像读取的cartesian force有问题。脚本出问题了

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