计算化学公社

标题: 请问关于stable判断波函数稳定性的问题 [打印本页]

作者
Author: dingniu2 时间: 2021-5-3 09:23
标题: 请问关于stable判断波函数稳定性的问题
请问各位，关于波函数稳定性的判断一直有个疑问，stable关键词用于判断波函数稳定性，而stable=opt则既可以判断波函数稳定性，又可以优化至稳定的波函数，如果波函数稳则定两者耗时相等，不稳定时通常也需要得到稳定波函数，这种耗时也是必须要花费的。那么stable关键词岂不是没有什么用了？任何情况都用stable=opt岂不更好？

作者
Author: sobereva 时间: 2021-5-3 09:45
也不是
stable=opt寻找稳定波函数的代价比较高。stable判断不稳定时，用户可以通过别的代价更低的方式（如调整初猜）试图得到稳定的波函数。

作者
Author: zjxitcc 时间: 2021-5-3 10:07
本帖最后由 zjxitcc 于 2021-5-3 10:32 编辑

（几乎）任何情况下stable=opt都更比stable更好。虽然有sob老师的这句话“stable判断不稳定时，用户可以通过别的代价更低的方式（如调整初猜）试图得到稳定的波函数”，但是你想啊：你咋知道你构造的初猜能一次性得到稳定的波函数呢？即使你有很好的初猜，你还是得老老实实做一次stable=opt的计算；只做stable计算是不够的，因为运气好的话是稳定的，运气不好的话不稳定怎么办？拼运气么？再调整初猜么？人脑能想出多少种初猜呢？又能想出多少种稳定的初猜呢？

想初猜是简单的，但是一次性就能构造出稳定的初猜，这在d和f单电子居多的过渡金属里是几乎不可能实现的（尤其是双核、多核过渡金属），可能你构造的初猜很合理，但就是不稳定，还是会经过一次短暂的波函数（指迭代圈数少）优化，到达一个能量更低的波函数。还有其他情况：同一个自旋多重度下需要几种不同的初猜，那算前咋知道哪种初猜更合理、哪种能量最低呢？仅用stable万一不稳定再手动补交计算？所以说还是得靠stable=opt。

那在双自由基、有机分子中，确实经常你能通过构建初猜或简单的guess=mix一次性得到稳定的波函数，但这其实是后验性的，你算完才知道，算前无法预测，最多只能是一种经验性规律。所以在你知道答案的时候stable是优于stabe=opt；在你不知道答案的时候，stable=opt优于stable。

作者
Author: sobereva 时间: 2021-5-3 21:06
比直接用stable=opt有很多更廉价的得到最终能收敛到基态波函数的方式
比如可以在更小基组下做stable=opt，然后把收敛的波函数当大基组的初猜
对于发现是对称破缺单重态的时候用guess=mix，或者片段初猜波函数
根据stable输出的负本征值的组态信息，基于stable任务已经收敛的轨道对轨道占据方式进行调换
等等

stable=opt在较高计算级别是很昂贵的，如果可能有其它做法得到稳定波函数，建议先尝试其它的

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