CP2K的k点测试不收敛 - 计算化学公社

CP2K默认的产生初猜波函数的方法是孤立原子密度的叠加，在整个任务一开始时会对每个KIND对应的原子在孤立状态下做计算，用PRINT_LEVEL MEDIUM时可以看到具体细节。计算孤立原子所用的电子组态可以通过&KIND里的MAGNETIZATION（定义alpha电子数减beta电子数）或者&BS（可以具体控制各个主量子数的各个角量子数上的alpha和beta电子数）来设置。初猜波函数中的孤立原子的组态越接近整个体系收敛的波函数中原子的实际组态则SCF越容易收敛、SCF不收敛的可能性越低。

如果体系有磁性，遇到SCF不收敛时一定别忘了尝试用MAGNETIZATION或&BS尽可能令初猜波函数中原子组态与实际情况接近，这对SCF是否能够顺利收敛有关键性影响。如果当前体系在文献（可用google搜索来找）或高通量数据库（如Material Project）中已经被计算过，可以参考那里面报道的原子自旋磁矩或原子组态来设。有时同一种元素在体系里处于不止一种化学环境，差异很大时（比如明显处于不同价态）可分成多个&KIND分别设置。