vasp计算aggas2能带及态密度及光学性质

《vasp计算aggas2能带及态密度及光学性质》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、VASP计算AgGaS2能带及态密度及光学性质第一步：构型优化1．准备四个输入文件POSCARINCARPOTCARKPOINTSPOSCAR:从ms中导入AgGaS2结构，选择CASTEP，file，save，并保存成原包。这样，得到一隐藏文件.cell,将它用编辑器打开，从中的到vasp所需的POSCAR信息，修改得到POSCAR。AgGaS2bulk1.000000000000000-2.79349994659424102.79349994659423705.20450019836425802.793499

2、9465942370-2.79349994659423905.20450019836425902.79349994659423802.7934999465942380-5.2045001983642580422Direct0.37500000000000000.40720000863075260.53220000863075260.87500000000000000.84279999136924740.4677999913692474............2.POTCAR:用PBE-GGA的赝势，提取，AgGa

3、O的赝势合并成一个赝势。（一般VASP有自带）3.下面是INCARSYSTEM=optimizationofAgGaS2LPLANE=.TRUE.NPAR=8ElecronicminimisationISTART=0LREAL=.FALSE.PREC=LowEDIFF=1e-4EDIFFG=-0.03IALGO=48NELMIN=4ISYM=0GGA=PBEISPIN=1NBANDS=120OUTPUTCONTROLLCHARG=.TRUE.LWAVE=.TRUE.LVTOT=.FALSE.IONICRELAXA

4、TIONNBLOCK=1NSW=1IBRION=-1DOSRELATED(disabled)ISMEAR=0(tetrahedron/gaussian/m-p)SIGMA=0.054.再下面KPOINTSA0M444000至此，四文件已准备好，进行计算，mpiexec-np8~/bin/vasp.4.5-mk-mp-pgivasp.out&（运行命令）（问题：计算前对某些参数的测试，比如截断能，晶格参数等，标准时什么？vasp运行后，不知道结果是否满足要求呢？）第二部：将构型优化后产生的CO

5、NTCAR文件拷贝为POSCAR文件mpiexec-np8~/bin/vasp.4.5-mk-mp-pgivasp.out&（运行命令）第三部：能带的计算.建立新的计算目录mkdirband创建目录cdband进入目录cp../INCAR.复制cp../PO*.cp../CHG*.对INCAR文件作如下修改：ICHARG=11NBANDS=120与上面构型优化时一致建立新的KPOINTS手动定义K点AgGaS251Rec（关键字要对齐第一行）0.50000.50000.50001.00000

6、.45000.45000.45001.0000......K点坐标获得途径布里渊区各特殊点的坐标可利用CASTEP获得然后编辑viinplin1特殊点lin2特殊点之间的取点数lin3第一个K点的坐标lin4第二个K点的坐标....（同上..）....保存运行vasp_kpointsvasp.out&（运行命令）在运行VaspbandDrawthebandstruc

7、tureofVASPTwofiles,EIGENVALandOUTCARareneeded----------------------------------------------------------------Spin-polarizedcalculation(T/F)?(是否自旋极化计算)fNumberoflinesalongtheBZ:(布里渊区特殊线数目)5Kpointsequenceusedtoseparatelines:(各线段始终点序号)11121314151SettheFermilevelt

8、ozero(T/F)?(是否费米设置为零点)t可以得到相应的band.dat用origin绘图然后标注各个特殊点Fermienergy=1.2182vexp=2.76ev(可能是参数设计太粗糙…)不过放大还算清楚..(用MS作图K网格相同)第四部:DOS的计算及其绘图其POSCARPOTCARKPOINTS可以相同于第一步优化设置INCARSYSTEM=optimizat