stable检查波函数稳定性问题

Jerryluo · 发表于 Post on 2018-6-7 17:14:41

请问，用M06L/TZVP scf=xqc opt freq得到最终收敛的结构。想检查波函数的稳定性，再用关键词M06L/TZVP stable geom=check guess=read，还是用M06L/TZVP scf=xqc stable geom=check guess=read算？
另外，opt的每一步用的波函数都是一样的吗？
opt的每一步不一定都用了qc，如果最后一步没有用qc，还需要检查波函数稳定性吗？

weiyiran · 发表于 Post on 2022-9-20 09:19:25

sobereva 发表于 2022-9-19 23:42
如果对最终结构下最终的波函数用stable=opt时提示不稳定，就用这个任务产生的稳定波函数结合guess=read o ...

明白了！感谢社长耐心解答！

sobereva · 发表于 Post on 2022-9-19 23:42:06

weiyiran 发表于 2022-9-19 15:31
社长好！您提到“为了确信稳定，可以对最后结构再测试稳定性”。这时候如果不稳定，应该如何解决呢？因为 ...

如果对最终结构下最终的波函数用stable=opt时提示不稳定，就用这个任务产生的稳定波函数结合guess=read opt继续做优化。反复如此，直到最终优化的结构又无虚频波函数又稳定

weiyiran · 发表于 Post on 2022-9-19 15:31:45

sobereva 发表于 2018-6-8 02:56
没事别乱加xqc
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

社长好！您提到“为了确信稳定，可以对最后结构再测试稳定性”。这时候如果不稳定，应该如何解决呢？因为stable和opt是不能同时用的。我猜想，可以对最终结构进行stable=opt使得波函数稳定，然后通过guess=read opt再进行优化，然后再检验稳定性。不知道这样子是否正确呢？谢谢您！

Jerryluo · 发表于 Post on 2018-6-9 10:33:07

终于明白了，感谢大家的耐心解答！

sobereva · 发表于 Post on 2018-6-9 03:56:36

Jerryluo 发表于 2018-6-8 11:19
谢谢sob老师，我的体系是含有金属的体系，按照您的帖子尝试了其他方法都没能解决scf收敛问题才用qc的 ...

1 建议同时加上guess=read，这样测试的才相当于优化过程中一直维持的波函数的状态，如果不加这个，可能当前结构下用默认的初猜得到的波函数和优化过程中一直维持的并不相同

2 确切来说是当前几何结构下当前得到的波函数的稳定性。

3 每一步的结构不同，显然用上一步收敛的波函数当初猜，在当前几何结构下还要做SCF，波函数自然会变化
guess=read和opt没直接联系

4 不需要。测最终结构的就行了。

anhuixu · 发表于 Post on 2018-6-8 14:32:00

波函数的收敛和稳定性是不同的概念。

应该是M06L/TZVP stable geom=check guess=read吧，测试你最终得到的收敛波函数的稳定性。
M06L/TZVP stable 会根据你的最终结构猜测一个初始波函数，迭代使其收敛，然后测试该收敛波函数的稳定性。这时，如果初猜不合理，波函数可能不收敛。

按道理，每一步都应该测试波函数稳定性；但这样做太麻烦，一般认为从稳定波函数出发，之后就都稳定了。但为了确信稳定，可以再对最后结构再测试稳定性。

Jerryluo · 发表于 Post on 2018-6-8 11:19:00

sobereva 发表于 2018-6-8 02:56
没事别乱加xqc
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

谢谢sob老师，我的体系是含有金属的体系，按照您的帖子尝试了其他方法都没能解决scf收敛问题才用qc的，所以想测试波函数稳定性。
1、我也可以这样理解吗，M06L/TZVP scf=xqc opt freq 得到最终收敛的结构坐标文件，再用输入文件：

M06L/TZVP stable
****
0 1
坐标

来测试波函数稳定性

2、stable测的是当前几何结构波函数的稳定性？

3、另外，您说opt每一步的结构都不同，波函数自然不同。opt每一步初猜波函数自动用上一步收敛的波函数。
那么，如果用了guess=read，结构最终收敛。根据opt每一步初猜波函数自动用上一步收敛的波函数，整个opt每一步的波函数不是就一样了吗？
我是不是可以这样理解，用了guess=read，opt每一步波函数都一样。没用guess=read的话，opt每一步都不一样？

4、您说。opt每一步的结构都不同，波函数自然不同。如果用了scf=xqc，如何确定每一步的波函数都稳定，每一步都需要stable吗？

sobereva · 发表于 Post on 2018-6-8 02:56:28

没事别乱加xqc
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331

之后你guess=read本身就是读取当前级别收敛的波函数，显然再用scf=xqc这种有可能解决不收敛问题的关键词更是完全多余

opt每一步的结构都不同，波函数自然不同

opt每一步初猜波函数自动用上一步收敛的波函数。有可能某个态在初始结构下是稳定的，但到了最终结构就不稳定了。为了确信稳定，可以再对最后结构再测试稳定性。

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

[Gaussian/gview] stable检查波函数稳定性问题