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

好的，非常感谢

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

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

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

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

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

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

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

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

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

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

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

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

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

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