源资培训课程-vasp流程输入文件

ViennaAb-initioSimulation输入和输出学习目理解VASP的计算掌握VASP各输入输出文件的掌握VASP各输入文件的创熟悉INCAR中各个熟记INCAR中需要更改的参数意2VASP结构优化输入文

自洽迭

原子3VASP软件包Fortran语言数学库文件主要代码贋势(LDA和一些用来处理计算数据的script和图形界面MedeA（含建模、计算及性质预测、数据处理4VASP基本输入文件POTCAR:文件包含了体系中各类元素的POSCAR:文件描述了所计算的体系的晶胞（包括基矢，晶格常数，原子类型，坐标位置INCAR:文件控制了VASP进行何种性质的5VASP基本输出文件介OUTCAR:计算结果主要输出文 (或MD)的能量信CONTCAR:原 或MD后的体系结构文件，分数坐标，用于续CHG和CHGCAR:电荷密度文WAVECAR:波函数DOSCAR电子态密度EIGENVAL:各k点本征值文件，能带数目，能量XDATCAR:每次原子弛豫后的每一步轨迹文IBZKPT:布里渊区中的kPCDAT:对关联函PROCAR和PROOUT:波函数投影或分解的LOCPOT:总的局域ELFCAR电子局域6在 版本中，VASP贋势库仅提供投影缀加平面波LDA和GGA-A，A_sv，A_pv和A_d（A为元素名称；s、p、d为原子轨道在具体做材料计算时，注意赝势种类的选取VASP贋势文件夹中包含两个压缩文件：potpaw_LDA和potpaw_PBE操作根据泛函及精确度选择需要的贋势文件catfile1file2file3…>POTCAR，将所有贋势文件合并成一个文POTCAR注意赝势的种类要一赝势使用的泛函要与INCAR中选择的泛函一致（PBE的泛函要选PBE的赝势8POSCAR文件示例一Pd（Fm- 1Selective FF

！任意文字注释（一般为体系名称！晶胞缩放系！晶格矢！原子种！原子个！是否要固定部分原子的！输入坐标的格式（direct模式！原子的位置，T为不固定，F为固9POSCAR文件示例二POSCAR1144！分数坐C！绝对标POSCAR文件示例四TTTTFFFFTTTTTTTTFFFFTTTTTTTTFFFFTTTT241TTTTTTTTFFFFTTTTTT提示输入坐标的格式（第八行）共有两种‘Direct’or‘D’or‘d’fordirect‘Cartesian’or‘C’or‘c’or‘K’or‘k’forcartesian各元素原子坐标的排列顺序必须与POSCAR中第6行的原序及POTCAR中的贋势顺序若需要固定部分原子，则可以在原子位置后边添加对x、y、z向是否固定的判断参数（T为不固定，F为固定常见格式（自动生成k点Auto0G

标题或注释行，无特别意0表示自动产生k以字母G开头表示按M-P方法以Γ点为中心自动产生k沿各个基失方向上分割各基失（K点网格）的点对所按网格分割产生的k点进行平移的量（这里不平移提示对六角晶系的结构，强烈推荐采用以Γ为中心按M-P网格产生kK点网格大小k1k2k3按晶格参数a,bc的长度进行合适选一般选取N1,N2和N3为奇数，以便产生的k点包含Γ点，（自动产K点的模式一般来说，k点越密越对于k点的需求，金属>>半导体，绝缘体，不过呢，很多时候主要VASP的INCAR文注释所计算的体系，以示说如果计 中有WAVECAR文件，则默认值为1，否则为0。可赋予值0|1|2|3。决定是否读入0:开始新的计算，按INIWAV初始化波函1:接着计算，通常用在测试ENCUT的收敛性以及计算结合能曲线(体积总能的关系:3如果ISTART=0，则默认值为2，否则为0。决定如何构造初始的电荷密01：从CHGCAR读入，并同原子密度进行 构造原子电子密度线性组合（LCAO），初始的电子密度由定11：读入自洽的CHGCAR，并进行能带计算或态密度的非自洽12ISTART和 推荐的做设置ISTART1ICHARG11一般都设置ISTART0ICHARG22.1自洽迭代步数和收敛标准：NELM,NELMINNELMDL2.3电荷密度混合方:IMIXAMIXAMIN,BMIX,AMIX-电子优化需要设置的参数有EDIFF和ENCUT其他参数使用默认。由NM:60IGODIGEDIFF:EDIFF是电子自洽收敛标准参数。默认值为1E-4。当连续两ENCUT:确定平面波截断能，默认值从POTCAR中读ALGO:确定电子优化的算法Normal则IALGO=38(也就是blockedDavidson方法:FastSMASS,TEBEG,TEEND,POMASS,原子如何移动以及步长和IBRION:决定了原子如何移动 。如果NSW=0或1，则默认值－1，否则为0。可赋予值为－1|0|1|2|3|5。IBRION=0表示进行子动力学模拟。为-1表示原子不移动。为其他值设置结构优化-1 1235用来计算HessianIBRION=0时，默认值为0，否则为2。可赋予值为0|1|2|3|4|5|6|7|0是否是否否1是trace是否否2是是是否否3是是是是是4是是是是否5是是否是否6是是否是是7是是否否是原子如何移动以及步长和POTIM:当IBION=1,2或3时，是力的一个缩放常数()0.5IBION=0MDNSW:子 大步数分子力学的数。认值0在行洽计算并精计算子所受H-力和应力。DIFF:原子 收敛的标准。默认值为EDIF。果的为则表前后次的自能之于EDIF原停止 如果为负，则原子所受的大的力小于EDIFG对，子止迟。分子动力学相关SMASS:确定分子动力学中原子的速度。默认值为-3。可赋予值为-3|2|1|0-3微正则系综(总的自由能守恒:-1在每NBLOCKTEBEG：分子动力学模拟时的初始温度。默认值为0TEEND：分子动力学模拟时的末态温度。默认值为POMASS：每类原子的质PSTRESS：设置加到体系的应力张量上的应力大1.smearing方法和参数：ISMEAR2计算态密度时能量范围和点数：EMIN,EMAX,smearing方法和参确定smearing的方法和展宽的参数。默认值为1和0.2。ISMEAR可赋予为5|4|3|2|1|0|N，它决定了如何确定每个波函数的占有数0Gaussiansmearing方法（氧化物，半导体1—N:Methfessel-PaxtonN阶smearing方法（金属-2:在INCAR中通过FERWE和 手动设-3:采用一个循环的计算来确定smearing的参数，在每一步中的searming设置通过SMEARINGS=smearing1smearing2-4:没有Blöch修正的四面体方-5:带Blöch修正的四面体方smearing方法和参数建ISMEAR=-SIGMA4.2&4.3计算态密度、分波态密度相关参EMIN和EMAX:计算态密度时的能量窗口(或能量范围LORBIT:控制非自洽计算的形S：UPREC:是控制计算精度最重要的一个参数，它决定了为Low,Medium,High,NormalAccurate。它的默认值为Medium(后两个只能在4.5以上版本中才起作用)。Medium。当要提高力和Stresstensor的计算精度时，可以采用High或Accurate，并手动设置ENCUT的值。磁性计MAGMOM每类原子的初始磁距，默认值为Nion*VOSKOWN:确定交换关联函数的关联部分是否采用逻辑参数和体系参LREAL:确定投影算法是否在实空间或倒空间计算。默认值为.FALSE.，可设置的值为.FALSE.|.TRUE|On|Auto|LCHARG确定是否在CHGCAR和CHG5强关联体系LDAU,LDAUTYPE,LDAUL,LDAUU,LDAUTYPE：具体加U的方法1|2|4。通常取2为默认值LDAUL：控制具体在哪个角量子数l（轨道）上加U；-1—不加1—对p轨道加U，2—对d轨道加U，3—对f轨道加LDAUU：U1,U2,……指定有效在位库伦相互作用的LDAUJ：J1,J2,……指定有效在位交换相互作用的参注意：LDAUL，LDAUU，LDAUJ这三个参数必须对体系中所有原都要定义。（包括加U和不加U的原子5杂化计算 ：确定是否进行杂化计算。取HFSCREEN：选择HSE06泛函PRECFOCK=精确，可选择Normal和AccurateINCAR文件——关键参2.电子迭代：NELMEDIFFENCUT；电子部分优化方法：3.原子弛豫：IBRION,NFREEISIFPOTIM,NSWEDIFFG；分子动力学：SMASS,TEBEG,TEEND,PSTRESS4.smearing方法和参数：ISMEAR,SIGMA；计算态密度时能量范围和点数:EMIN,EMAX,NEDOS；计算分波态密度的参数:LORBIT。5.计算精度：PREC；磁性计算：ISPIN,MAGMOM,INCAR文件——SYSTEMCopperISTART=0ICHARG=LREAL=离子弛豫部分

！自己定义 系的名！开始新的！从原子 的！计算精度，一般！实空间投影，一般Auto或者IBRION=2NSW=100POTIM=EDIFFG=-

！离子弛豫的方法，CG方！控制晶格变！离子弛豫的步！离子移动的时间步！离子弛豫的force的收敛标准，＜0.05时，结构收INCAR文件——ENCUT=200eVNELM=EDIFF=1E-04LWAVE=.TRUE.