计算化学公社

标题: vasp计算态密度问题 [打印本页]

作者
Author: asl1994 时间: 2017-10-12 19:57
标题: vasp计算态密度问题
计算PuO2态密度时，和文献比较发现有一段完全消失了。图b和c来源于文献。图a是我自己算的，计算细节照搬图c的（除了泛函我用的PBE而别人用的PBE-sol）。现在问题是b和c尽管计算方法差别很大，包括泛函、+U值以及磁性设置都不同，但他们算的DOS曲线的形状大体差不多。但是我算出来的为什么在E>6eV的区域完全没有态密度呢？EMAX的设置应该没有问题的，也就是说那一段的DOS值确实都是零。请问这可能是什么原因导致的呢？我做DOS非自洽计算的INCAR如下： ISTART = 1
ICHARG = 11
ISPIN = 2
VOSKOWN = 1
LREAL = .F.
NCORE = 4
LSORBIT = .True.
LMAXMIX = 6
ISYM = 0
SAXIS = 0 0 1

LORBIT = 10

# Accuracy controls:
PREC = normal
ENCUT = 500

# Electronic loop:
ALGO = Normal # Normal (Davidson)
EDIFF = 1E-07
NELM = 100    # Maximum number of electronic SC steps (default:60)

# DOS related values:
ISMEAR = -5       # tetrahedron method with Blochl corrections (use a G-centered k-mesh)
SIGMA = 0.200000 # eV
NBANDS = 140
EMIN = -5
EMAX = 20
NEDOS = 2000

# Ironic relaxation:
IBRION = -1
ISIF = 2    # Cell fixed
NSW = 0 # Maximum number of ionic steps
POTIM = 0.500000 (scaling constant)
EDIFFG= -0.01

# DFT +U
LDAU = .T.
LDAUTYPE = 2
LDAUL = 3 -1
LDAUU = 6.35 0
LDAUJ = 0 0

# Write flags:
LCHARG = .True.
LWAVE = .F.
LELF = .F.
LVTOT = .F.
LVHAR = .F.

作者
Author: asl1994 时间: 2017-10-12 20:00
更正：图b和c的位置放反了

作者
Author: 卡开发发 时间: 2017-10-12 21:58
非自洽的NBANDS可能不够，这样更高的能带就出不来。

作者
Author: Jack 时间: 2017-10-13 09:50
请问，PDOS怎么归一化使total等于各个分波的加和，谢谢！

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