VASP基本参数设置
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
SYSTERM=hafnium oxide
ISTART=0
ICHARG=2
NWRITE=2
LWAVE=.False.
PREC=high???? 默认值Medium在4版本,在5版本normal
specify high --mannual 46
#LVTOT=.True.
#IDIPOL=4
ALGO=Fast
LREAL=Auto
#LELF=.TRUE.
LCHARG=.FALSE.
################electronic relaxation############
NELM=80
ENCUT=500?????
EDIFF=1E-05
#ISPIN=2
###################ion relaxtion################
NSW=
IBRION=2
ISIF=2
EDIFFG=-0.005
####################################################
ISMEAR=0
SIGMA=0.05
#LORBIT=11
#NEDOS=
#NBANDS =
ISYM=0
#NELECT
###################bader charge analysis#################### #LAECHG=.Ture.
#NGXF
#NGYF
#NGZF
###################NEB######################################## #ICHAIN=
#IMAGES=
#SPRING=
#LCLIMB=
#POTIM=
#IOPT=
vasp incar 主要参数设置
SYSTEM=B N ------作业说明语句,即任务的名字
IALGO=38 ------指定电子自洽计算的算法 38 Davidson 48 RMM-DIIS 算法
NELM=200 ------自洽迭代次数,一般默认40次
NELMIN=4 ------最小迭代次数,表面或动力学计算是应增大默认为2次
#NELMDL=-12 ------开始几步不自洽处理,默认IALGO=8时为-5,IALGO=48时为-12,其他0 #ISTART=0 ------是否使用已有波函数。有WAVECAR是默认为1使用,否则为0 新的作业
#ALGO=very_fast ----同IALGO
EDIFF=1E-05 ----总能计算允许误差,默认1E-04
ISMEAR=0 ------确定电子的部分占据数,半导体一般取小于0,金属大于1
LREAL=Auto ------投影操作是在实空间还是在倒空间,.FALSE.倒,.TRUE.实。AUTO实优化ISYM=0 ------是否具有对称性,默认1,具有,0 不具有。
#ICHARG=11 ------如何产生初始电荷密度,默认ISTART=0时取2,LCAO。否则0,从初始波函
#LDIAG = .FALSE. -----是否执行子空间对角化。IALGO=8时为.TRUE.
NELECT ------总电子数,如果系统不是电中性的就必须设置,所带电荷作为均一的背景电子气考虑。
# ! plot charge density
#LPARD = .TRUE. --------是否计算部分电荷。
#IBAND = 176 177 ----指定要计算的带
#KPUSE = 1 7 -----要考虑的K点
#LSEPB = .TRUE.
#LSEPK = .TRUE.
# ! change configuration
#FERWE = 159*1.0 1.0 0.0 80*0.0--电子占据数从INCAR中读入并在计算中不变
#FERDO = 159*1.0 0.0 1.0 80*0.0--电子自旋极化计算(按每个K点每个能带排列)
#SIGMA=0.2 -----------展开宽度,默认ISMEAR=1 SIGMA=0.2
#ISPIN=2 ------------指定电子体系的总自旋。
# ! optimize
IBRION=2 ------离子是否运动,-1不运动但做NSW外循环。0动力学模拟,1准牛顿法离子弛豫------2 CG法离子弛豫,3 采用衰减二阶运动方程离子弛豫,
ISIF=2 -----是否离子弛豫及是否改变原胞。
NSW=200 -----离子运动的步数,默认0
POTIM=0.2 ----IBRION=0时时间步长,=123作用在力上的比例系数,默认0.5
#NSW=0
#IBRION=-1
EDIFFG=-1E-02 ---离子弛豫运动的结束条件,默认为EDIFF*10
#PREC=HIGH ------控制ENCUT,FFT格子和实空间投影时RPOT的默认值。
# ! parallel calculation
#LPLANE = .TRUE. ---对带和平面波系数并行
#NPAR = 1 ---- 1 只对平面波系数并行。并行计算band的节点数,每一个节点计算一个
band当然可以提高并行效率,减少通讯量,不过貌似现在硬件的主要限制
还是内存,而这个选项的使用可能会大幅增加内存的需求。
#LSCALU = .FALSE. ---
#NSIM = 4 ----有N条带一起优化。
PREC=normal ------控制ENCUT,FFT格子和实空间投影时RPOT的默认值
VOSKOWN = 1 ---对交换关联泛函是否采用VWM的形式,默认0不用而用标准泛函,1采用对磁性有用
AMIX = 0.2 ---线性混合参数
BMIX = 0.00001 ---Kerker混合方法的截断波失
AMIX_MAG = 0.8 ---磁的线性混合参数
BMIX_MAG = 0.00001--磁的Kerker混合方法的截断波失
#ADDGRID = .TRUE.
LORBIT=11 ---是否输出投影波函数到PROCAR和PROOUT中
#LVTOT=.TRUE. ---是否输出总局域势。
LWAVE=.FALSE. ---是否输出波函数(WAVECAR)
LCHARG=.FALSE. ---是否输出电荷密度(CHGCAR和CHG)
#NBANDS= 130 ---指定计算的总能带数。计算需要大量的能带(空带),至少要1个空带
(否则VASP会给出警告)。一般NBANDS=NELECT/2+NIONS/2以上可得到较
精确的结果,如果内存不够就只好减少NBANDS,在牺牲精度和体系大小之
间平衡了。
# ! electric field
EFIELD=0.3
LDIPOL=.TRUE.
IDIPOL=1 -------- calculate monopole/dipole and quadrupole corrections
1-3 只计算第一/二/三晶矢方向适于slab的计算
4 全部计算尤其适于就算孤立分子
DIPOL=0.5 0.5 0.5 ----原胞中心坐标。