AgGaS2的計算全过程.docVIP

VASP计算AgGaS2能带及态密度及光学性质 第一步：构型优化 1．准备四个输入文件 POSCAR INCAR POTCAR KPOINTS POSCAR: 从ms中导入AgGaS2结构，选择CASTEP，file，save，并保存成原包。这样，得到一隐藏文件.cell, 将它用编辑器打开，从中的到vasp所需的POSCAR信息，修改得到POSCAR。 AgGaS2 bulk 1.000000000000000 -2.7934999465942410 2.7934999465942370 5.2045001983642580 2.7934999465942370 -2.7934999465942390 5.2045001983642590 2.7934999465942380 2.7934999465942380 -5.2045001983642580 4 2 2 Direct 0.3750000000000000 0.4072000086307526 0.5322000086307526 0.8750000000000000 0.8427999913692474 0.4677999913692474 . . . . . . . . . . . . 2.POTCAR:用PBE-GGA的赝势，提取，Ag Ga O的赝势合并成一个赝势。（一般VASP有自带） 3.下面是INCAR SYSTEM = optimization of AgGaS2 LPLANE=.TRUE. NPAR= 8 Elecronic minimisation ISTART = 0 LREAL = .FALSE. PREC = Low EDIFF = 1e-4 EDIFFG = -0.03 IALGO = 48 NELMIN = 4 ISYM = 0 GGA = PBE ISPIN = 1 NBANDS = 120 OUTPUT CONTROL LCHARG = .TRUE. LWAVE = .TRUE. LVTOT = .FALSE. IONIC RELAXATION NBLOCK = 1 NSW = 1 IBRION = -1 DOS RELATED (disabled) ISMEAR = 0 (tetrahedron/gaussian/m-p) SIGMA = 0.05 4.再下面KPOINTS A 0 M 4 4 4 0 0 0 至此，四文件已准备好，进行计算， mpiexec -np 8 ~/bin/vasp.4.5-mk-mp-pgi /dev/null vasp.out （运行命令） （问题：计算前对某些参数的测试，比如截断能，晶格参数等，标准时什么？vasp运行后，不知道结果是否满足要求呢？） 第二部：将构型优化后产生的CONTCAR文件拷贝为POSCAR文件 mpiexec -np 8 ~/bin/vasp.4.5-mk-mp-pgi /dev/null vasp.out （运行命令） 第三部：能带的计算.建立新的计算目录 mkdir band 创建目录 cd band 进入目录 cp ../INCAR . 复制 cp ../PO* . cp ../CHG* . 对INCAR文件作如下修改： ICHARG =11 NBANDS=120 与上面构型优化时一致 建立新的KPOINTS 手动定义K点 AgGaS2 51 Rec（关键字要对齐第一行） 0.5000 0.5000 0.5000 1.0000 0.4500 0.4500 0.4500 1.0000 . . . . . . K点坐标获得途径 布里渊区各特殊点的坐标可利用CASTEP获得 然后编辑 vi inp lin1 特殊点 lin2 特殊点之间的取点数 lin3 第一个K点的坐标 lin4 第二个K点的坐标 . . . . （同上..） . . . . 保存 运行 vasp_kpoint