计算化学公社

标题: 有关stable关键词的疑问 [打印本页]

作者
Author: 2238663 时间: 2020-9-24 09:58
标题: 有关stable关键词的疑问
本帖最后由 2238663 于 2020-9-24 11:02 编辑

各位老师好！我在对一个配合物进行电荷分布计算时，采用pop关键词，但由于SCF难以收敛，参考了Sob老师的《解决SCF不收敛问题的方法》一文。
我在单点分析任务中采用了SCF=conver=6，文中提到这样放宽收敛对单点能的计算基本不产生影响，可以大胆使用，但是我用了stable测出我的波函数不稳定……
对于电荷分布的计算，使用该关键词结果是否还可信呢？

作者
Author: 喵星大佬 时间: 2020-9-24 10:30
所以和stable关键词有什么关系

作者
Author: sobereva 时间: 2020-9-24 10:34
当前问题和stable关键词没有丝毫联系
conver=6的波函数一般来说可用

作者
Author: 2238663 时间: 2020-9-24 11:02

sobereva 发表于 2020-9-24 10:34
当前问题和stable关键词没有丝毫联系
conver=6的波函数一般来说可用

我加了stable关键词说我的波函数不稳定我不知道电荷分布结果是否可信

作者
Author: 喵星大佬 时间: 2020-9-24 11:03

2238663 发表于 2020-9-24 11:02
我加了stable关键词说我的波函数不稳定我不知道电荷分布结果是否可信

不稳定就stable=opt变稳定咯

作者
Author: 2238663 时间: 2020-9-24 11:47

喵星大佬发表于 2020-9-24 11:03
不稳定就stable=opt变稳定咯

恩…体系比较大是直接用坐标计算单点任务不涉及结构优化

作者
Author: sobereva 时间: 2020-9-24 12:05

2238663 发表于 2020-9-24 11:47
恩…体系比较大是直接用坐标计算单点任务不涉及结构优化

本来stable=opt和结构优化也没有任何联系

作者
Author: wzkchem5 时间: 2020-9-24 12:28

2238663 发表于 2020-9-24 11:47
恩…体系比较大是直接用坐标计算单点任务不涉及结构优化

stable=opt的opt指的是优化波函数使其变得稳定，和单独写的opt关键字（结构优化）没有任何关系
而且并不是说不设scf=conver=6就一定能稳定，不要把不稳定怪罪到把conver放松到6上，不稳定性主要是因为你的初猜和最后的稳定波函数有定性差别，比如空间对称性不同，或该破缺的对称性（空间或自旋对称性）没有破缺

作者
Author: zjxitcc 时间: 2020-9-24 12:47

2238663 发表于 2020-9-24 11:02
我加了stable关键词说我的波函数不稳定我不知道电荷分布结果是否可信

波函数不稳定，要看是哪种类型的不稳定，如果是内部不稳定（internal instability），布局分析不可信。

如果你用了RDFT，且显示的是RHF -> UHF不稳定，那就不好讲了。此时若你稳定的UDFT波函数能量只比RDFT能量低一点点，且<S**2>非常接近于理想值，那么两个波函数的布局分析结果是很接近的，是可信的。若你稳定的UDFT波函数能量比RDFT能量低很多，那RDFT结果不怎么可信。

作者
Author: 2238663 时间: 2020-9-24 13:17
先谢谢各位老师的指导了。对波函数稳定性的测试是internal instability，请教一下：我下一步使用stable=opt后，波函数稳定后，布局分析产生的电荷值可信吗？

作者
Author: wzkchem5 时间: 2020-9-24 14:25

2238663 发表于 2020-9-24 13:17
先谢谢各位老师的指导了。对波函数稳定性的测试是internal instability，请教一下：我下一步使用stable=opt ...

只能说是不会因为波函数稳定性的问题导致不可信，因为决定电荷是否可信的还有其他因素，比如你的体系是否有强相关，是否容易受到自相互作用误差的影响，等等。
强相关体系有instability的概率比较大，所以反过来讲，出现instabililty的时候建议检查一下体系是否有强相关特征，如果有明显的强相关，就需要用多参考态方法，DFT算出来的无论如何是不可信的。

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