求问stable，opt=stable，stable=opt的差异

xyyan1207

各位老师好，
图片中展示的是sob老师博文中提到的，我想问一下stable，opt=stable，stable=opt，这三个关键词之间有什么细节差异？
若优化，还需要再单独添加opt吗？ opt=stable是既优化又检测波函数稳定性吗？

zjxitcc

youxx189 发表于 2023-10-24 11:12
老师您好  请问假如按第一种方法的话  用stable=opt得到了稳定的波函数 ，再读取这一步的chk后优化结构， ...

是有可能的，所以结构优化收敛后，应当再次检验波函数稳定性。

youxx189

sobereva 发表于 2018-1-9 10:41
第一种第二种皆可。第三种则stable=opt完全没对优化的初猜波函数产生影响
搞清楚几个关键词的含义，该 ...

老师您好  请问假如按第一种方法的话  用stable=opt得到了稳定的波函数 ，再读取这一步的chk后优化结构，会不会优化后的结构的波函数又不稳定了啊

youxx189

wzkchem5 发表于 2023-7-5 19:55
但凡是你不能仅靠化学常识就确定波函数稳定的情况，都需要加stable=opt。和高自旋态、自旋极化单重态没有 ...

我理解了  谢谢老师的解答

wzkchem5

youxx189 发表于 2023-7-5 09:34
您好  请教一下  是不是高自旋态和自旋极化单重态都需要加stable=opt呢

但凡是你不能仅靠化学常识就确定波函数稳定的情况，都需要加stable=opt。和高自旋态、自旋极化单重态没有必然关系。
换句话说，如果你不能从化学常识判断肯定不需要加，那就是需要加；但是如果完全不知道什么情况下波函数肯定稳定，那就应该先去学习什么情况下波函数一定稳定

youxx189

您好  请教一下  是不是高自旋态和自旋极化单重态都需要加stable=opt呢   

zh7729

wzkchem5 发表于 2023-6-17 04:07
对，不一样。只要最后得到的结构既收敛了、波函数又稳定即可，达到这一目的所经历的过程不重要。换句话说 ...

我明白了，感谢老师的耐心解答！

wzkchem5

zh7729 发表于 2023-6-16 14:49
老师，如果按照您提到的在结构优化前是否进行波函数稳定性检测，是不是可以分为两种方法，第一种方法是在 ...

对，不一样。只要最后得到的结构既收敛了、波函数又稳定即可，达到这一目的所经历的过程不重要。换句话说，优化收敛以后必须至少做一个stable；优化前是否做stable=opt只影响计算时间而不影响计算正确性。具体而言，如果所有结构的波函数都是稳定的，而你一开始做了stable=opt，会浪费稳定性分析时间；如果所有结构的波函数都不是稳定的，而你一开始没做stable=opt，会既浪费结构优化时间又浪费稳定性分析时间，因为优化收敛后做stable=opt发现波函数不稳定后，必须继续做结构优化到收敛。至于初始结构和收敛结构一个波函数稳定、另一个不稳定的情况，什么时候一开始做stable=opt省时，什么时候不做省时，这个你自己照理类推，我不一一分类讨论了

zh7729

wzkchem5 发表于 2022-9-11 18:35
结构优化结束时必须加，因为后续计算一般用到的是结构优化之后的波函数，而非结构优化之前的波函数，所以 ...

老师，如果按照您提到的在结构优化前是否进行波函数稳定性检测，是不是可以分为两种方法，第一种方法是在结构优化前用stable进行波函数稳定性测试，发现波函数不稳定，然后做stable=opt，得到稳定的波函数，再读取含有稳定波函数的check文件进行结构优化；第二种方法是先做结构优化，再做stable发现波函数不稳定，然后做stable=opt得到稳定波函数，再读取含有稳定波函数的check文件再次进行结构优化。因为做stable=opt时的结构不同，所以这两种方法通过stable=opt得到的稳定波函数是不是不一样？那读取不同的稳定波函数，再进行结构优化后得到的结果是不是也不一样呢？谢谢老师

sobereva

wzy 发表于 2022-9-11 17:46
老师，我有一个疑问，stable是刚开始的时候加上还是结构优化结束后用stable来测试波函数稳定性？

对于电子结构复杂、收敛到不稳定波函数概率较大的体系，优化前应当对优化任务所基于的波函数做稳定性测试，优化后对收敛结构下收敛的波函数再做一次稳定性测试，最为稳妥。前者避免初始波函数是错的，导致做无意义的优化，后者确保最终结果是有意义的。

wzkchem5

wzy 发表于 2022-9-11 10:46
老师，我有一个疑问，stable是刚开始的时候加上还是结构优化结束后用stable来测试波函数稳定性？

结构优化结束时必须加，因为后续计算一般用到的是结构优化之后的波函数，而非结构优化之前的波函数，所以显然必须保证结构优化后的波函数是稳定的，那就必须在结构优化后加。
但是这不意味着结构优化之前一定不用加。如果做distortion-interaction analysis等计算，需要用到未优化的结构下的能量或者波函数，那么优化前的波函数也必须确保稳定，此时就是至少要做两次稳定性测试才能保证结果完全可靠。
此外如果你预期这个体系的波函数不稳定的概率比较大（比如大于50%），那么建议在优化前也测试稳定性（但不是必需的，如果不做，可能会浪费计算时间，但不会导致结果不可用），原因是如果优化前的波函数是稳定的，那么优化后的波函数大概率也是稳定的（虽然不必然是稳定的），所以如果波函数不稳定的概率比较大，那么只优化后做稳定性测试就会导致很可能需要优化2次，而预先做稳定性测试就很大概率只需要优化1次，又因为稳定性测试比优化快，所以总的来说很可能可以省时间。相比之下，如果波函数很可能是稳定的，那么不管是否预先做稳定性测试，都很可能只需要优化1次，这样就白白浪费了一次稳定性测试的计算时间。不过再次强调一遍，这种情况下，不管优化前是否做稳定性测试，都不影响结果的可用性，至多只是浪费计算时间，所以搞错了也没关系。前面说的distortion-interaction analysis等情况下是不能搞错的。

wzy

sobereva 发表于 2018-1-9 10:08
scf=qc或xqc和opt不可能有任何冲突，完全不相干
stable和opt别图省事搅合在一起用，不然程序糊涂自己也 ...

老师，我有一个疑问，stable是刚开始的时候加上还是结构优化结束后用stable来测试波函数稳定性？

sakura1109

zjxitcc 发表于 2021-10-22 12:34
不稳定的波函数 一般没有意义。稳定至少能确保在轨道系数的空间是个local minimum，总比不稳定的有意义多 ...

非常感谢各位老师的耐心解答

wzkchem5

sakura1109 发表于 2021-10-22 03:16
请问就算找到了最近的稳定点也不能保证是全局最小值吧？这样理解对吗？
如果无法确保找出来的是全局最稳 ...

找到全局最稳定的SCF解是NP-hard问题，如果有一种算法能在多项式时间内找到任意分子的全局最稳定解，那么你就证明了P=NP，可以去领Clay Institute的百万美元大奖了。
需要指数时间找到全局最稳定解的算法目前已经有了，但是极其耗时，目前做的极限可能也就10个电子。所以绝大部分情况下，不能严格确保你找到的是全局最稳定解，只能靠以下办法来增加全局最稳定的概率：（1）做stable=opt，确保得到局部最稳定解；（2）找到多个局部最稳定解，选其中能量最低者（但由于以上原因，注定需要指数时间才能保证把所有局部最稳定解都找全）；（3）分析SCF收敛的波函数，比如看布居数、自旋密度等，确保结果符合化学常识。

biogon

sakura1109 发表于 2021-10-22 10:16
请问就算找到了最近的稳定点也不能保证是全局最小值吧？这样理解对吗？
如果无法确保找出来的是全局最稳 ...

对于scf多解的体系稳定点不一定是波函数空间的全局极小点，而且局部极小点是有意义的
scf=qc收敛到波函数空间的鞍点上都是正常事，这玩意也就是起到让波函数收敛的作用罢了，不能保证收敛到的是极小点