如何看待《关于高斯计算精度参数的常见误解》这篇文章的观点

北大-陶豫 · 发表于 Post on 2022-1-15 10:55:02

关于高斯计算精度参数的常见误解https://mp.weixin.qq.com/s/EC3ZO8oz5A6UnQE9_Q99yA
文中的观点可取吗？实际计算中还有哪些要注意的呢？

Dolantin · 发表于 Post on 2024-6-1 11:09:31

sobereva 发表于 2024-6-1 10:20
找个有代表性的结构，用ultrafine积分格点精度重新计算，做个对比怼回去、证实他的揣测是完全多余的就完 ...

游客，本帖隐藏的内容需要积分高于 25 才可浏览，您当前积分为 0

sobereva · 发表于 Post on 2024-6-1 10:20:39

Dolantin 发表于 2024-6-1 01:41
**** 本内容被作者隐藏 ****

找个有代表性的结构，用ultrafine积分格点精度重新计算，做个对比怼回去、证实他的揣测是完全多余的就完了

Dolantin · 发表于 Post on 2024-6-1 01:41:13

sobereva 发表于 2022-1-18 09:14
Gaussian的几何优化默认收敛限对于一般目的已经足够了，在我来看是得当的。如果有人鼓吹Gaussian默认收敛 ...

游客，本帖隐藏的内容需要积分高于 25 才可浏览，您当前积分为 0

北大-陶豫 · 发表于 Post on 2022-2-16 16:39:27

喵星大佬发表于 2022-2-16 12:37
M06-2X这个破玩意，cphf=grid=fine相当有必要

居然这么久了还给我的帖子挖坟，ピエン

chrinide · 发表于 Post on 2022-2-16 13:38:54

喵星大佬发表于 2022-2-16 12:37
M06-2X这个破玩意，cphf=grid=fine相当有必要

参数太多约束太多，这个泛函可以部分的视为病态模型，原始的SCAN也是，各种数值问题，我记得用PSI4曾经测试SCAN的时候，好些个含重元素体系（X2C下）用过（500-974）都
不能完全收敛，居然要到（750-1202）才能完全收敛。

喵星大佬 · 发表于 Post on 2022-2-16 12:37:47

M06-2X这个破玩意，cphf=grid=fine相当有必要

wzkchem5 · 发表于 Post on 2022-1-18 16:10:34

sobereva 发表于 2022-1-18 01:53
我经常发现ORCA基于DIIS commutator判断收敛的时候，哪怕用strongSCF甚至tightSCF，达到收敛时能量收敛精 ...

从ORCA 5.0.3开始，UHF/UKS的DIIS收敛精度会变高一些，因为以前是用alpha和beta的commutator的和来判断收敛，所以在两个commutator异号的时候收敛会比较粗糙，现在改成了两者都必须收敛

冰释之川 · 发表于 Post on 2022-1-18 13:28:28

liyuanhe211 发表于 2022-1-18 12:58
相比于“Reply to comment to”，这叫“Correction of”更合理

233333

liyuanhe211 · 发表于 Post on 2022-1-18 12:58:08

mfdsrax2 发表于 2022-1-18 09:38
作者更新了《Reply to comment to “关于高斯计算精度参数的常见误解”》
https://mp.weixin.qq.com/s/kPd ...

相比于“Reply to comment to”，这叫“Correction of”更合理

mfdsrax2 · 发表于 Post on 2022-1-18 09:38:22

作者更新了《Reply to comment to “关于高斯计算精度参数的常见误解”》
https://mp.weixin.qq.com/s/kPdcqZ4qFbn1YTeLRUiS4w

sobereva · 发表于 Post on 2022-1-18 09:14:49

量化小菜鸡发表于 2022-1-17 22:48
我在想要Gaussian的默认收敛限要是都不行，那ORCA的结果不是按这个逻辑更差。。。

Gaussian的几何优化默认收敛限对于一般目的已经足够了，在我来看是得当的。如果有人鼓吹Gaussian默认收敛限太宽，甚至觉得tight才应该是默认的，则是非常偏激、很不成熟的，对各种类型体系都做过很多研究的经验丰富的研究者一般是不会认同的。

题外话：作为程序开发者不能无脑地把收敛限设得过严，比如非要试图保证所有情况下虚频严格清零什么的，而必须权衡考虑绝大多数应用场景下结果的精度、耗时、达到收敛的难度。如果把比如int=superfine、opt=tight当成默认关键词，虽然会给极个别用户的极个别研究带来好处（比如对柔性非常强的体系出现虚频的概率可能降低为原先的<1/10），却会给绝大多数用户的体验带来显著的坏处（平均耗时和遇到几何优化不收敛的情况增加N倍）。所以对于收敛限这个问题，在开发程序，以及写经验文章给出建议的时候，切不可光从自己的角度看待问题，更得从绝大多数程序使用者的角度全面、理性看待。

sobereva · 发表于 Post on 2022-1-18 08:53:47

wzkchem5 发表于 2022-1-18 00:10
这个不好比，因为ORCA默认的检查SCF收敛的方法是看总能量和单电子能量，当开启DIIS的时候是看DIIS commut ...

我经常发现ORCA基于DIIS commutator判断收敛的时候，哪怕用strongSCF甚至tightSCF，达到收敛时能量收敛精度都还不是很高，容易心里没底。我觉得要是始终也能判断能量收敛精度就好了，毕竟这是算单点时候特别关心的，比如strongSCF时至少<1E-6 Hartree。

量化小菜鸡 · 发表于 Post on 2022-1-18 00:12:45

wzkchem5 发表于 2022-1-18 00:10
这个不好比，因为ORCA默认的检查SCF收敛的方法是看总能量和单电子能量，当开启DIIS的时候是看DIIS commut ...

好的，谢谢老师。

wzkchem5 · 发表于 Post on 2022-1-18 00:10:31

量化小菜鸡发表于 2022-1-17 16:20
好的老师，我印象之前ORCA的TightOpt才对应高斯默认opt的MaxD要求。SCF来说ORCA就更低了，VeryTightSCF对 ...

这个不好比，因为ORCA默认的检查SCF收敛的方法是看总能量和单电子能量，当开启DIIS的时候是看DIIS commutator，当开启SOSCF的时候是看轨道梯度，而高斯是看密度矩阵。
你有兴趣的话可以用高斯和ORCA跑同一个体系，看最后一步的能量变化、密度矩阵的RMS和MAX变化，看看到底差多少。我们因为众所周知的原因没法用高斯，没办法自己来做这个事情

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Gaussian/gview] 如何看待《关于高斯计算精度参数的常见误解》这篇文章的观点