这是本文档旧的修订版!
这个报错信息就是指现在碎片的电子占据出现了非整数,在该碎片的out文件中可以看到例如:
B A N D E N E R G Y R A N G E S The energies below are in Hartree. band min(up) max(up) occup(up) min(dn) max(dn) occup(dn) 10 -0.15579 -0.15579 1.00000 -0.14616 -0.14616 1.00000 11 -0.13304 -0.13304 0.20000 -0.10815 -0.10815 0.00000 12 -0.13304 -0.13304 0.20000 -0.10815 -0.10815 0.00000 13 -0.13304 -0.13304 0.20000 -0.10815 -0.10815 0.00000 14 -0.13304 -0.13304 0.20000 -0.10815 -0.10815 0.00000 15 -0.13304 -0.13304 0.20000 -0.10815 -0.10815 0.00000 16 -0.04274 -0.04274 0.00000 -0.03537 -0.03537 0.00000 17 -0.04274 -0.04274 0.00000 -0.03537 -0.03537 0.00000 18 -0.04274 -0.04274 0.00000 -0.03537 -0.03537 0.00000 19 0.05265 0.05265 0.00000 0.05668 0.05668 0.00000 20 0.16818 0.16818 0.00000 0.17820 0.17820 0.00000 21 0.16818 0.16818 0.00000 0.17820 0.17820 0.00000 22 0.16818 0.16818 0.00000 0.17820 0.17820 0.00000 23 0.16818 0.16818 0.00000 0.17820 0.17820 0.00000 24 0.16818 0.16818 0.00000 0.17820 0.17820 0.00000 25 0.46662 0.46662 0.00000 0.47856 0.47856 0.00000 26 0.79965 0.79965 0.00000 0.82037 0.82037 0.00000 27 0.79965 0.79965 0.00000 0.82037 0.82037 0.00000 28 0.79965 0.79965 0.00000 0.82037 0.82037 0.00000 ========================================================================================
其中
所以可以看到11、12、13、14、15能带上,分别占据了自旋向上的0.2电子,这样进行片段分析,计算NOCV是不正确的。
需要想办法让电子占据成为整数。可以尝试人为指定占据:BAND里面如何指定占据方式
这种情况下,碎片计算完毕,如果结果没有问题,就可以进行整体计算。整体计算的时候,Mutilevel - fragment窗口,为片段选中之前计算得到的rkf文件(*.results/band.rkf),这样整体计算就直接计算,不再重复片段计算了。