计算化学公社

标题: 求助Gaussian中几何优化收敛限的问题 [打印本页]

作者
Author: qinchen 时间: 2023-12-28 11:20
标题: 求助Gaussian中几何优化收敛限的问题
在Author: sobereva 的文章“Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法(http://sobereva.com/278)”一文中3.1节提到的“值得一提的是优化过程中即便位移没收敛，但只要受力小于收敛限的100倍也自动被Gaussian算作收敛。”，请问：这里的收敛限指的是位移的收敛限还是力的收敛限呢？

不能直接给作者留言，没有权限，所以发在这里了，期待Author: sobereva 及其他前辈的回复哦！谢谢！！

作者
Author: mfdsrax2 时间: 2023-12-28 11:22
但只要受力小于收敛限的100倍也自动被Gaussian算作收敛

这说的还不够清楚吗？

作者
Author: qinchen 时间: 2023-12-28 15:28

mfdsrax2 发表于 2023-12-28 11:22
但只要受力小于收敛限的100倍也自动被Gaussian算作收敛

这说的还不够清楚吗？

Maximum Force          0.000401    0.000450    YES
RMS    Force          0.000197    0.000300    YES
Maximum Displacement    0.002231    0.001800    NO
RMS    Displacement    0.001325    0.001200    NO

文中前面提到：
笔者认为，如果数值比默认情况的收敛限高不太多，不超过2倍，且没有虚频出现，一般还是可以接受的。

后面说：
值得一提的是优化过程中即便位移没收敛，但只要受力小于收敛限的100倍也自动被Gaussian算作收敛。

我能力有限，那我觉得我对这里指代的理解是有点不清楚，因为不清楚所以来请教大家，谢谢哦！

另外：“位移or力”，1-2个字
      “这说的还不够清楚吗？”9个字

作者
Author: xuhj199508 时间: 2023-12-28 15:53

qinchen 发表于 2023-12-28 15:28
Maximum Force 0.000401 0.000450 YES
RMS Force 0.000197 0.0 ...

“如果数值比默认情况的收敛限高不太多，不超过2倍，且没有虚频出现，一般还是可以接受的。”这个指代的是freq结果中的最大位移和均方根位移不超过收敛限的2倍。
“但只要受力小于收敛限的100倍也自动被Gaussian算作收敛。”这个指代的是opt过程中即使最大位移和均方根位移没达到收敛限，但是最大受力小于收敛限的100倍也算优化正常结束。

作者
Author: qinchen 时间: 2023-12-28 17:02

xuhj199508 发表于 2023-12-28 15:53
“如果数值比默认情况的收敛限高不太多，不超过2倍，且没有虚频出现，一般还是可以接受的。”这个指代的 ...

谢谢您的回复哦！！

前面的内容理解，后面的部分有疑惑，，那也就是说“但是最大受力小于收敛限的100倍”这句指的应该是“最大受力小于最大受力的收敛限的100倍，对吗？再次感谢！！！

作者
Author: sobereva 时间: 2023-12-29 05:06
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并清楚、准确反映出帖子具体内容，避免有任何歧义和含糊性，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题 ““Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法(http://sobereva....” 改了，以后务必注意！下次将删帖+扣分处理。

作者
Author: sobereva 时间: 2023-12-29 05:08

qinchen 发表于 2023-12-28 17:02
谢谢您的回复哦！！

前面的内容理解，后面的部分有疑惑，，那也就是说“但是最大受力小于收敛限的100 ...

当前考察的是最大受力，此语境下“收敛限”不是指“最大受力的收敛限”还能指什么？

作者
Author: 哲人 时间: 2023-12-29 19:26
你可以去翻翻那个out文件看看

作者
Author: hebrewsnabla 时间: 2023-12-29 19:45
此处的“100倍”应该是 0.01倍的意思……或许这个让你产生了疑惑。

作者
Author: qinchen 时间: 2024-1-2 16:05

hebrewsnabla 发表于 2023-12-29 19:45
此处的“100倍”应该是 0.01倍的意思……或许这个让你产生了疑惑。

谢谢

作者
Author: qinchen 时间: 2024-1-2 16:05

hebrewsnabla 发表于 2023-12-29 19:45
此处的“100倍”应该是 0.01倍的意思……或许这个让你产生了疑惑。

谢谢

作者
Author: qinchen 时间: 2024-1-2 16:06

sobereva 发表于 2023-12-29 05:06
如置顶的新社员必读贴、论坛首页的公告栏、版头的红色大字非常明确所示，求助帖必须在帖子标题明确体现出此 ...

谢谢

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