计算化学公社

 找回密码 Forget password
 注册 Register
Views: 540|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[VASP] vasp计算dos时电子步收敛很慢是什么原因?

[复制链接 Copy URL]

13

帖子

0

威望

95

eV
积分
108

Level 2 能力者

本帖最后由 zss2022 于 2023-4-1 14:20 编辑

16个原子的磁性体系GGA+U计算,第一次不设置能量范围时,计算DOS正常收敛;第二次设置能量范围(-1,1), 共1000个点,OSZICAR里的电子步收敛特别慢,64个电子步才收敛到0.01. 搞不懂是什么原因。下面贴出INCAR文件。KPOINTS取的(16,16,16)。
---------------------------INCAR---------------------------------
Global Parameters
ISTART = 1
ICHARG = 11         (Default: 2 if ISTART = 1;else: 0)
LREAL  = .FALSE.      (Projection operators: automatic)
ENCUT  =  700        (Cut-off energy for plane wave basis set, in eV)
PREC   =  Accurate       (Precision level)
ADDGRID = .TRUE.        (Increase grid; helps GGA convergence)
GGA_COMPAT = .FALSE.
ALGO = Very
# SYMPERC = 1E-08
LWAVE = .FALSE.
LCHARG = .FALSE.
NEDOS = 1000
EMIN = -1
EMAX = 1
ISYM = 0
NPAR = 8

Electronic Relaxation
ISMEAR =  -5            (Gaussian smearing; metals:1)
SIGMA = 0.2
NELM   =  90           (Max electronic SCF steps)
NELMIN =  6            (Min electronic SCF steps)
EDIFF  =  1E-08        (SCF energy convergence; in eV)

Magnetism+SOC
ISPIN  =  2            (Non-Spin polarised DFT)
MAGMOM = 4 -4 -4 4 12*0
LORBIT = 11

GGA+U
LDAU         = .TRUE.
LDAUTYPE     = 2
LDAUL        = 2 -1 -1
LDAUU        = 4.0 0.0 0.0
LDAUJ        = 0.0 0.0 0.0
LMAXMIX      = 4
LASPH        = .TRUE.

----------------------------------END------------------------------

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 20:24 , Processed in 5.767890 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list