这里会显示出您选择的修订版和当前版本之间的差别。
后一修订版 | 前一修订版 | ||
adf:non-integer-occupations [2018/11/16 10:25] – 创建 liu.jun | adf:non-integer-occupations [2023/05/03 22:17] (当前版本) – liu.jun | ||
---|---|---|---|
行 1: | 行 1: | ||
======logfile报错:Error detected: The fragments need to have integer occupations====== | ======logfile报错:Error detected: The fragments need to have integer occupations====== | ||
- | 这个报错信息就是指现在碎片的电子占据出现了非整数,解决这个问题即可。可以看看电子占据方式,是不是想要的占据方式。 | + | 这个报错信息是指某个碎片的电子占据出现了非整数,而这是EDA分析禁止出现的。 |
+ | =====解决方法===== | ||
+ | 在该碎片的out文件中可以看到例如: | ||
+ | <code bash> | ||
- | 可以尝试人为指定占据:[[adf: | + | 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)是这条带自旋向上的能量最大值(本例中因为只计算了Gamma点,因此min up和max up是相等的) | ||
+ | * min(down)是这条带自旋向下的能量最大值(本例中因为只计算了Gamma点,因此min down和max down是相等的) | ||
+ | * occup(up)是自旋向上的占据电子数目 | ||
+ | |||
+ | 所以可以看到11、12、13、14、15能带上,分别占据了自旋向上的0.2电子,这样进行片段分析,计算NOCV是不正确的。但这个例子是极难解决的,因为这个体系对称性极高从11~15号能带,能量绝对地简并(可能是一个孤立的金属原子,有一个d电子),因此无论怎么指定(指定占据的格式参考:[[adf: | ||
+ | < | ||
+ | Occupations | ||
+ | 1 11 // 10 | ||
+ | END | ||
+ | </ | ||
+ | 或者 | ||
+ | < | ||
+ | Occupations | ||
+ | 1 10 0 1 // 10 | ||
+ | END | ||
+ | </ | ||
+ | 最后都不会收敛,例如到100次迭代时,可能这个电子随机地占据在某个轨道上。101次迭代,又随机在另一个轨道上。也许用户可以考虑这种方式反复计算,然后在SCM → View中查看对应的占据轨道,到底是哪个d轨道被占据了,然后将这个结果,用到后续的EDA-NOCV计算中。 | ||
+ | |||
+ | 碎片的电子占据成为整数之后,就可以进行整体计算。整体计算的时候,Mutilevel - fragment窗口,为片段选中之前计算得到的rkf文件(*.results/ |