VASP计算AgGaS2能带及态密度及光学性质

AgGaS2能带、态密度和光学性质的VASP计算步骤1:配置优化1.准备四个输入文件波萨卡角波萨卡关键点POSCAR:从ms导入AgGaS2结构，选择CASTEP，归档，保存，并将其保存为原始包。这样，一个隐藏的文件。单元格是用编辑器获取和打开的，从编辑器中获取vasp所需的POSCAR信息，并对其进行修改以获取POSCAR。AgGaS2批量1.0000-2.42410 2.42370 5.425802.42370 -2.42390 5.425902.42380 2.42380 -5.425804 2 2直接的0.00000 0.07526 0.075260.00000 0.92474 0.92474。2.POTCAR:提取了PBE-GGA的赝势，并将银镓氧的赝势合并为一个赝势。(VASP有自己的)3.以下是肉色系统=AgGaS 2的优化LPLANE=.TRUE。NPAR=8电子最小化ISTART=0LREAL=.FALSE。PREC=低EDIF=1e-4EDIFFG=-0.03IALGO=48NELMIN=4ISYM=0PBEISPIN=1NBANDS=120输出控制LCHARG=.TRUE。LWAVE=.TRUE。LVTOT=.FALSE。离子松弛NBLOCK=1新南威尔士=1IBRION=-1相关操作系统(禁用)ISMEAR=0(四面体/高斯/m-p)西格玛=0.054.以下关键点A0M4 4 4此时，这四个文件已准备好进行计算。mpiexec-np8 /bin/vasp . 4.5-MK-MP-PGI/dev/null vasp . out(运行命令)(问题：在计算前测试一些参数的标准是什么，例如截断能量、晶格参数等。)？vasp运行后，您不知道结果是否符合要求吗？)第二部分：将配置优化后生成的CONTCAR文件复制到POSCAR文件中mpiexec-np8 /bin/vasp . 4.5-MK-MP-PGI/dev/null vasp . out(运行命令)第三部分：能带的计算。新计算目录的建立Mkdir波段创建目录Cd带进入目录Cp.incar.replicationcp./PO*。cp./CHG*。按如下方式修改CORTI文件：ICHARG=11NBANDS=120与上述配置优化一致创建新的关键点来手动定义k个点AgGaS251记录(对齐第一行的键)0.5000 0.5000 0.5000 1.00000.4500 0.4500 0.4500 1.0000。获取k点坐标的方法利用CASTEP可以获得布里渊区特殊点的坐标。然后编辑vi inpLin1特殊点lin2特殊点之间的点数lin3的第一个k点的坐标。lin4第二个k点的坐标。(同上.)。救援运行vasp _ kpoints vasp.out(运行命令)生效Vaspband画出VASP的乐队结构需要两个文件，EIGENVAL和OUTCAR-自旋极化计算？(自旋极化计算)fBZ:沿线的线数5用于分隔：行的点序列1 11 21 31 41 51将费米能级设为零？(费米是否设置为零)tDat可以用原点绘制。然后标记每个特殊点费米能量=1.2182v exp=2.76ev(可能参数设计太粗糙)，但放大仍然清晰。(k网格与毫秒相同)第四部分：DOS计算和绘图它的poscapot carkpoints可以与第一个优化步骤相同。设置角色系统=AgGaS 2的优化NPAR=8LPLANE=.TRUE。电子最小化ISTART=3(从WAVECAR读取波函数)LREAL=.FALSE。PREC=中等EDIF=1e-4EDIFFG=-0.02IALGO=48NELMIN=3ISYM=1ISIF=3ISPIN=1ISMEAR=0西格玛=0.1NBANDS=120LORBIT=11(预测为标准偏差)NEDOS=501(拒绝服务点)ICHARG=11(从CHG读取电荷密度)EMIN=-20EMAX=20输出控制LCHARG=.TRUE。LWAVE=.TRUE。LVTOT=.FALSE。离子松弛NBLOCK=1新南威尔士=1IBRION=-1POTIM=0.05如上运行vasp命令然后写DOSf(自旋极化计算与否)9(轨道数量)6(监督事务司将计算)-1(负数表示被研究原子的DOS)1(原子序数)-1(正态时原子轨道的DOS)2(同上)-13111518保存并运行vaspdos(如果没有设置路径。指定您自己的路径)有25号堡的文件。运行cdos98需要一个附加的dos号你可以得到DOS1 DOS2 DOS3.用原点单独绘制.(由于错误的操作系统编号，未绘制常规(p)和特殊(p).)第5部分：光学性质的计算(a)线性计算1.还执行配置优化(在相同能带计算中的结构优化)。只有计算中的CORTI文件需要包含关键词：ISIF=3在所有计算中必须使用PAW赝势2.配置优化后，将CONTCAR复制到POSCAR。在印度NPAR=1(光学性质必须按能带的顺序计算)LOPTICS=。真(计算光学性质)新南威尔士=1IBRION=-1(不需要配置优化)计算完成后(正常)，光学(用于计算线性光学性质)和动量矩阵被获得。(用于计算NLO，即非线性光学特性)3.编辑OPCTR文件中常用的关键词及其含义。ISYMM=2OMMIN=0OMMAX=20NEDOS=4000NBCON=200LJDOS=.TRUE。LDOS=.TRUE。LKRAMERS=.TURE。=0.002l搜索=.TRUE。EMINSEARCH=0EMAXSEARCH=20剪刀=0.54.创建一个新的目录，并复制光学光学光学控制，关键点，在它下面。命令cp./optic.optctr，kpoints，poscar。修改关键点创建临时目录来复制vasp计算所需的四个输入文件。运行独立vasp。(由于56台计算机未配备独立版本，因此此步骤通过复制到其他计算机来完成)单机操作命令/bin/vasp . 4.5-MK-sp-PGI/nullvasp . out将IBZKPT复制到关键点运行vasp命令：MPI execnp8 /bin/vasp . 4.5-MK-MP-PGI/null vasp . out正常操作将产生2个文件电动助力转向系统-存储介电功能数据JDOS -存储JDOS数据Vasp_lo正在运行用origin软件绘制基于EPS数据的图形在此基础上，可以获得其他线性光学特性：(1)折射率(2)吸附常数(3)能量损失系数：(4)消光系数(5)双折射曲线(6)光电导变化曲线(2)非线性光学性质1.非线的操作需要以下三个文件关键点、本征值和动量矩阵是三个文件。关键点文件的内容与计算的线性光学特性相同。本征值文件计算vasp的输出文件要运行vasp_nlo，请输入以下数据：动量矩阵存储过渡矩阵元素的文件名177.32细胞体积2.76实验带隙(如果为0，则不需要带隙校正)1初始能带(一般为1)38端接频带如果f为假，则只计算静态倍频系数2.运行后将得到静态倍频系数的计算结果。(这些数据似乎不是自己保存的)d(1，1) SHG(实部)=1.60120和(imag部)0.00000(pm/V)d(1，2) SHG(实部)=-2.37106和(imag部)0.00000(pm/V)d(1，3) SHG(实部)=-2.72264和(imag部)0.00000(pm/V)d(1，4) SHG(实部)=10.98424