计算化学公社

标题: 金属氧化物ELF图为啥只能看到氧看不到金属 [打印本页]

作者
Author: qwe1832518773 时间: 2021-1-16 23:00
标题: 金属氧化物ELF图为啥只能看到氧看不到金属
用VASP算了一个Fe2O3表面的ELF，想用看一下Fe和O之间的成键情况，可是放到VESTA一看isosurfaces全在O原子上，Fe原子都没有，是我INCAR设置的问题吗？
(, 下载次数 Times of downloads: 31)
也做过2D截面，但截出来的原子全是氧原子，用VESTA里有一个shift slice by mouse wheel反复查看，完全找不到Fe，就和上面这个3D是对应的，但是看别人文章里的ELF明明金属是可以出现的
(, 下载次数 Times of downloads: 21)

是我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

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

作者
Author: cokoy 时间: 2021-1-19 11:30
换个包含内层电子的赝势应该就好了。

作者
Author: qwe1832518773 时间: 2021-1-22 14:12
前几天发不了言，后面经过尝试发现换成sv的赝势就好了，DFT+U也做了尝试，但+U对于描述ELF没啥用

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