体系带有一个正电荷，读取中性时的波函数SCF不收敛的问题

zhoumeimei · 发表于 Post on 2021-9-1 09:28:19

体系带有一个正电荷，算过中性的同结构，相同的计算方法泛函基组，读取的中性时的波函数。SCF一直不收敛，请问这种情况有没有什么好的解决办法？

zhoumeimei · 发表于 Post on 2021-9-1 16:41:08

wzkchem5 发表于 2021-9-1 16:32
加KDIIS不会改变，改格点、thresh/tcut会改变，所以可以先尝试KDIIS

好的，非常感谢

wzkchem5 · 发表于 Post on 2021-9-1 16:32:15

zhoumeimei 发表于 2021-9-1 09:30
谢谢老师的回复
我先试一下这些方法，看看能不能收敛。但是由于我算的是一个体系，考究每个体系的能量， ...

加KDIIS不会改变，改格点、thresh/tcut会改变，所以可以先尝试KDIIS

zhoumeimei · 发表于 Post on 2021-9-1 16:30:44

wzkchem5 发表于 2021-9-1 15:44
如果你算的是一个很大的闭壳层团簇的自由基阳离子的话（从你xyz文件的名字判断的），估计没啥好办法。可以 ...

谢谢老师的回复
我先试一下这些方法，看看能不能收敛。但是由于我算的是一个体系，考究每个体系的能量，改变这些参数会影响能量的值吧

wzkchem5 · 发表于 Post on 2021-9-1 15:52:29

zhoumeimei 发表于 2021-9-1 04:27
没有直接算过，读中性的波函数是以为以它为初猜会计算量会小一点，本来原子数就挺多的，带电更不好收敛

第一次算阳离子的时候，读中性波函数不是问题，因为默认初猜PModel相当于以中性原子密度的叠加为初猜密度，其实也是个中性体系，而且比你收敛的中性体系还要差。所以第一次算阳离子，读中性体系的波函数是有可能节省计算量的。但是第一次不收敛以后，后续计算不应该再从中性开始，而应该以前一次计算的最后一步的波函数开始，除非有证据表明前一次的波函数定性错误。
另外注意没有“带电更不好收敛”这个结论。开壳层比闭壳层难收敛，带电荷越负越难收敛，所以你的阳离子比中性体系难收敛不是因为带电，而是因为是开壳层。很多人长久以来只研究那些中性情况下是闭壳层的分子，然后观察到它的自由基阳离子、自由基阴离子都收敛困难，于是以讹传讹说带电一定会导致收敛变困难，其实只要试一个中性情况是开壳层的体系，就会发现根本不是这回事

wzkchem5 · 发表于 Post on 2021-9-1 15:44:28

如果你算的是一个很大的闭壳层团簇的自由基阳离子的话（从你xyz文件的名字判断的），估计没啥好办法。可以先试试加KDIIS，如果不行再加大格点，再不行就把thresh和tcut设小一些

zhoumeimei · 发表于 Post on 2021-9-1 15:19:14

sobereva 发表于 2021-9-1 12:16
根本就不会小。二者电子结构相差那么多，还指望中性时是多好的初猜。
分清楚阳离子和阴离子，明明电子数 ...

嗯嗯好的，以后这样的错误不会载犯了，谢谢老师

sobereva · 发表于 Post on 2021-9-1 12:16:24

zhoumeimei 发表于 2021-9-1 11:27
没有直接算过，读中性的波函数是以为以它为初猜会计算量会小一点，本来原子数就挺多的，带电更不好收敛

根本就不会小。二者电子结构相差那么多，还指望中性时是多好的初猜。
分清楚阳离子和阴离子，明明电子数越少往往越容易收敛，你当前的行为完全是自找麻烦。

zhoumeimei · 发表于 Post on 2021-9-1 11:27:11

zjxitcc 发表于 2021-9-1 09:56
那你为啥要读取中性时的波函数。。。直接默认初猜试了没。。

没有直接算过，读中性的波函数是以为以它为初猜会计算量会小一点，本来原子数就挺多的，带电更不好收敛

zjxitcc · 发表于 Post on 2021-9-1 09:56:05

那你为啥要读取中性时的波函数。。。直接默认初猜试了没。。

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

[ORCA] 体系带有一个正电荷，读取中性时的波函数SCF不收敛的问题

浏览过的版块