计算化学公社

标题: 体系带有一个正电荷，读取中性时的波函数SCF不收敛的问题 [打印本页]

作者
Author: zhoumeimei 时间: 2021-9-1 09:28
标题: 体系带有一个正电荷，读取中性时的波函数SCF不收敛的问题
体系带有一个正电荷，算过中性的同结构，相同的计算方法泛函基组，读取的中性时的波函数。SCF一直不收敛，请问这种情况有没有什么好的解决办法？

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

作者
Author: zhoumeimei 时间: 2021-9-1 11:27

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

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

作者
Author: sobereva 时间: 2021-9-1 12:16

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

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

作者
Author: zhoumeimei 时间: 2021-9-1 15:19

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

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

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

作者
Author: wzkchem5 时间: 2021-9-1 15:52

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

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

作者
Author: zhoumeimei 时间: 2021-9-1 16:30

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

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

作者
Author: wzkchem5 时间: 2021-9-1 16:32

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

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

作者
Author: zhoumeimei 时间: 2021-9-1 16:41

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

好的，非常感谢

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