计算化学公社

标题: 改变k点后还能继续使用之前的波函数和电荷密度么？ [打印本页]

作者
Author: ra2ghgzh 时间: 2017-9-14 09:45
标题: 改变k点后还能继续使用之前的波函数和电荷密度么？
我先用一个较低的k点优化结构并得到了WAVECAR和CHGCAR，然后我将k点调高并设置ISTART=1, ICHARG=1。
我自己试了一下，发现k点调高后ISTART=1, ICHARG=1得到的结果和ISTART=0, ICHARG=2得到的结果结构上完全一样，能量之差1E-5个eV。
能不能确定这样做是没问题的？

作者
Author: 卡开发发 时间: 2017-9-14 14:38
如果提升k点后ISTART=0和ISTART=1做构型优化完全一样，能量差别也很小，那么没有关系。但一般情况下不推荐这样做，有时候有可能导致构型稍有差异（肉眼未必能看出来）。一般而言，ISTART=1和ICHARG=11会通过波函数构造密度或直接读取密度来进行非自洽迭代（一般通过Harris functional的方法之类的），如果收敛标准也不足够高，有可能会收敛不到真正的基态上。

不太建议一开始随便去取一个低的k点，最好做一下收敛性测试，k点的收敛一般很快。如果体系需要的k很大，可以加一些展宽来加速k点收敛以减小k点使用，但不能太大，自洽算能量的时候减小展宽增大k点即可。

作者
Author: hakuna 时间: 2017-9-15 14:12
其实吧，这样做也没什么不可以，俺在旧社会时曾经这么做过，测试好K点后，接下来的计算，先用较少的k点，快速获得波函数，然后在此基础上正常K点下计算。当时条件差，优化个结构都是以天为单位计算，先少量k点先算得波函数，再正常k点计算，目的就是省点时间。
现在大家都富豪了，没人愿意这么“多此一举“了

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