计算化学公社

标题: 求助：CoPc和石墨烯的dos计算 [打印本页]

作者
Author: lyyylyyyl 时间: 2026-1-7 21:18
标题: 求助：CoPc和石墨烯的dos计算
各位老师同学大家好！
最近在计算CoPc和富勒烯（180C）放一起，间隔3.5A的结构的dos，用的是两步法，先scf计算，k点设置的111；再把scf计算生成WAVECAR和CHGCAR文件放入第二步计算的文件夹，k点设置为331。算出来后用vaspkit生成d带中心，发现d带中心为正，不知道是方法和参数哪里有问题，望各位帮忙解惑，感谢！
这是scf的INCAR：# control
SYSTEM=Cu2O_Au
SYMPREC=1E-10
ISTART=1
ICHARG=1
#IWAVE=1
LCHARG=.T.
LWAVE=.T.
LREAL=Auto
NPAR=4

#SCF
IALGO=48
EDIFF=1E-5
ISMEAR=0
SIGMA=0.1
PREC=Normal
NELM=150
NELMIN=10
ENCUT=520
LMAXMIX=4

#OPT
EDIFFG=-0.02
NSW=0
IBRION=2 # 0: MD, 1: Newton, 2: CG
ISIF=2 # 3 for bulk, 2 for surface and others
# IDIPOL=3 # 1: x, 2:y, 3, z, 4: all
POTIM=0.1

这是dos第二步的INCAR：# control
SYSTEM=Cu2O_Au
SYMPREC=1E-10
ISTART=1
ICHARG=11
#IWAVE=1
LCHARG=.T.
LWAVE=.T.
LREAL=Auto
NPAR=4

#SCF
IALGO=48
EDIFF=1E-5
ISMEAR=-5
SIGMA=0.1
PREC=Normal
NELM=150
NELMIN=10
ENCUT=520
LMAXMIX=4

#OPT
EDIFFG=-0.02
NSW=0
IBRION=2 # 0: MD, 1: Newton, 2: CG
ISIF=2 # 3 for bulk, 2 for surface and others
# IDIPOL=3 # 1: x, 2:y, 3, z, 4: all
POTIM=0.1

#DOS
LORBIT=11
NBANDS=704
EMIN=-30
EMAX=10
NEDOS=4001
#LDAU=T
#LDAUTYPE=2
#LDAUL=-1 2 -1
#LDAUU= 0 2 0
#LDAUJ= 0 0 0
#LMAXMIX=4

#spin- polar
ISPIN=2
#MAGMOM=8*0 6*0 1*0 2*2

作者
Author: shinnashiori 时间: 2026-1-10 15:09
1.  你是直接搭好的结构直接去跑了自洽没有做结构优化吗？
2.  NSW = 0 时候， IBRION = -1 应该要这样设置吧
3.  你自洽得到 WAVCAR 和 CHGCAR 文件了，做 dos 的时候应该把 LWAV 和 LCHA 给关上了
4.  有 Co 的话，最好还是要加上磁矩和 U 值吧

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