金属氧化物ELF图为啥只能看到氧看不到金属

qwe1832518773

用VASP算了一个Fe2O3表面的ELF，想用看一下Fe和O之间的成键情况，可是放到VESTA一看isosurfaces全在O原子上，Fe原子都没有，是我INCAR设置的问题吗？

也做过2D截面，但截出来的原子全是氧原子，用VESTA里有一个shift slice by mouse wheel反复查看，完全找不到Fe，就和上面这个3D是对应的，但是看别人文章里的ELF明明金属是可以出现的


是我INCAR设置的问题吗？INCAR如下：
ISTART = 0
ICHARG = 1
ENCUT = 400
EDIFF = 1E-4
ISYM = 0
NELMIN = 3
LREAL= Auto
ISMEAR = 0
SIGMA = 0.2
ALGO = Fast
VOSKOWN = 1
NPAR = 2

LWAVE = .FALSE.
LVDW = .TRUE.
LELF = .TRUE.
NGX = 108
NGY = 128
NGZ = 300

卡开发发

可能的原因是这些，尚不确定以下哪些因素是主导：
1、ISPIN是反铁磁性的，必须考虑自旋极化，且需要做磁矩初始猜测。
2、我不确定你PAW用的是何种，也许需要半芯态的PAW才能够一定程度显示出离子芯对ELF的贡献。
3、Fe2O3中的Fe可能需要进行DFT+U才有可能正确体现价电子的定域化特性，这方面PBE可能真的不太行。

cokoy

换个包含内层电子的赝势应该就好了。

qwe1832518773

前几天发不了言，后面经过尝试发现换成sv的赝势就好了，DFT+U也做了尝试，但+U对于描述ELF没啥用