youxx189 发表于 2023-10-24 11:12 是有可能的，所以结构优化收敛后，应当再次检验波函数稳定性。 |
sobereva 发表于 2018-1-9 10:41 老师您好 请问假如按第一种方法的话 用stable=opt得到了稳定的波函数 ，再读取这一步的chk后优化结构，会不会优化后的结构的波函数又不稳定了啊 |
wzkchem5 发表于 2023-7-5 19:55 我理解了 谢谢老师的解答 |
本帖最后由 wzkchem5 于 2023-7-5 12:56 编辑 youxx189 发表于 2023-7-5 09:34 但凡是你不能仅靠化学常识就确定波函数稳定的情况，都需要加stable=opt。和高自旋态、自旋极化单重态没有必然关系。 换句话说，如果你不能从化学常识判断肯定不需要加，那就是需要加；但是如果完全不知道什么情况下波函数肯定稳定，那就应该先去学习什么情况下波函数一定稳定 |
| 您好 请教一下 是不是高自旋态和自旋极化单重态都需要加stable=opt呢 |
wzkchem5 发表于 2023-6-17 04:07 我明白了，感谢老师的耐心解答！ |
zh7729 发表于 2023-6-16 14:49 对，不一样。只要最后得到的结构既收敛了、波函数又稳定即可，达到这一目的所经历的过程不重要。换句话说，优化收敛以后必须至少做一个stable；优化前是否做stable=opt只影响计算时间而不影响计算正确性。具体而言，如果所有结构的波函数都是稳定的，而你一开始做了stable=opt，会浪费稳定性分析时间；如果所有结构的波函数都不是稳定的，而你一开始没做stable=opt，会既浪费结构优化时间又浪费稳定性分析时间，因为优化收敛后做stable=opt发现波函数不稳定后，必须继续做结构优化到收敛。至于初始结构和收敛结构一个波函数稳定、另一个不稳定的情况，什么时候一开始做stable=opt省时，什么时候不做省时，这个你自己照理类推，我不一一分类讨论了 |
wzkchem5 发表于 2022-9-11 18:35 老师，如果按照您提到的在结构优化前是否进行波函数稳定性检测，是不是可以分为两种方法，第一种方法是在结构优化前用stable进行波函数稳定性测试，发现波函数不稳定，然后做stable=opt，得到稳定的波函数，再读取含有稳定波函数的check文件进行结构优化；第二种方法是先做结构优化，再做stable发现波函数不稳定，然后做stable=opt得到稳定波函数，再读取含有稳定波函数的check文件再次进行结构优化。因为做stable=opt时的结构不同，所以这两种方法通过stable=opt得到的稳定波函数是不是不一样？那读取不同的稳定波函数，再进行结构优化后得到的结果是不是也不一样呢？谢谢老师 |
wzy 发表于 2022-9-11 17:46 对于电子结构复杂、收敛到不稳定波函数概率较大的体系，优化前应当对优化任务所基于的波函数做稳定性测试，优化后对收敛结构下收敛的波函数再做一次稳定性测试，最为稳妥。前者避免初始波函数是错的，导致做无意义的优化，后者确保最终结果是有意义的。 |
wzy 发表于 2022-9-11 10:46 结构优化结束时必须加，因为后续计算一般用到的是结构优化之后的波函数，而非结构优化之前的波函数，所以显然必须保证结构优化后的波函数是稳定的，那就必须在结构优化后加。 但是这不意味着结构优化之前一定不用加。如果做distortion-interaction analysis等计算，需要用到未优化的结构下的能量或者波函数，那么优化前的波函数也必须确保稳定，此时就是至少要做两次稳定性测试才能保证结果完全可靠。 此外如果你预期这个体系的波函数不稳定的概率比较大（比如大于50%），那么建议在优化前也测试稳定性（但不是必需的，如果不做，可能会浪费计算时间，但不会导致结果不可用），原因是如果优化前的波函数是稳定的，那么优化后的波函数大概率也是稳定的（虽然不必然是稳定的），所以如果波函数不稳定的概率比较大，那么只优化后做稳定性测试就会导致很可能需要优化2次，而预先做稳定性测试就很大概率只需要优化1次，又因为稳定性测试比优化快，所以总的来说很可能可以省时间。相比之下，如果波函数很可能是稳定的，那么不管是否预先做稳定性测试，都很可能只需要优化1次，这样就白白浪费了一次稳定性测试的计算时间。不过再次强调一遍，这种情况下，不管优化前是否做稳定性测试，都不影响结果的可用性，至多只是浪费计算时间，所以搞错了也没关系。前面说的distortion-interaction analysis等情况下是不能搞错的。 |
sobereva 发表于 2018-1-9 10:08 老师，我有一个疑问，stable是刚开始的时候加上还是结构优化结束后用stable来测试波函数稳定性？ |
zjxitcc 发表于 2021-10-22 12:34 非常感谢各位老师的耐心解答 |
sakura1109 发表于 2021-10-22 03:16 找到全局最稳定的SCF解是NP-hard问题，如果有一种算法能在多项式时间内找到任意分子的全局最稳定解，那么你就证明了P=NP，可以去领Clay Institute的百万美元大奖了。 需要指数时间找到全局最稳定解的算法目前已经有了，但是极其耗时，目前做的极限可能也就10个电子。所以绝大部分情况下，不能严格确保你找到的是全局最稳定解，只能靠以下办法来增加全局最稳定的概率：（1）做stable=opt，确保得到局部最稳定解；（2）找到多个局部最稳定解，选其中能量最低者（但由于以上原因，注定需要指数时间才能保证把所有局部最稳定解都找全）；（3）分析SCF收敛的波函数，比如看布居数、自旋密度等，确保结果符合化学常识。 |
sakura1109 发表于 2021-10-22 10:16 对于scf多解的体系稳定点不一定是波函数空间的全局极小点，而且局部极小点是有意义的 scf=qc收敛到波函数空间的鞍点上都是正常事，这玩意也就是起到让波函数收敛的作用罢了，不能保证收敛到的是极小点 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-26 18:28 , Processed in 0.286917 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites