VASP计算实例

1、VASP计算实例目录一、氢气分子H2键长的计算31. 基本文件32. 鹰势细的选择33. 截断能ENCUT 锤的选择44. KPOINTS参数选择55. 对晶格行优化6二、Si晶体晶格常数计算81. 膜势类型选择82. 截断能（ENCUT）参数的选定93. KPOINTS参数的选定114. SIGMA参数的选定125. 晶格常数计算结果13三、Si元素单原子能量计算141. 由内聚能倒推单原子能量142. 基本文件143. 单原子能量计算15四、Si的VASP力学常数计算161. 计算所需文件162. 计算与数据处理163. 计算所用到的公式：18五、SI晶体的电子结构191. 采用VASP计

2、算能带的步骤192. 电荷分布计算结果20能带计算和结果213. 态密度计算和结果21六、Si晶体介电函数和光学性质的计算221. 计算步骤222. 用到的文件233. 计算结果26七、VASP的声子谱计算291. 计算步骤292. 基本文件303. 声子谱、声子态密度计算和结果334. 热学性质计算和结果34八、化合物键长计算351. 计算步骤352. 基本文件35氢气分子H2键长的计算1. 基本文件准备基本文件 INCAR、POTCAR. POSCAR、KPOINT 以及脚本文件 encut、k、optimize2. 鹰势类型的选择输入文件如卜其中参数要靠经验初选INCAR：System

3、= F2ISTART = 0ICHARG = 2NELMDL = 5ISEAR = 0SIGMA = 0. 1PREC = AccurateKPOINTS：Automatic mesh 0M1 1100 0POSCAR：0115.00.000. 000.0014.00.000.000.0013.01D0.000. 000. 00分别选用五个腐势文件进行计算。计算方法：把INCAR.POSCAR.POTCAR.KP0INTS四个文件分别放入同一个文件夹中，在终端中输入“VASP”得到五个能量。GGA-91E0=36620018E+01GGA-CAE0=35466311E+01PAWE0=3441

4、5631E+01PAW-GGAE0=34948621E+01PBEE0=35753116E+01选择单原子能最最低的，品体最稳定，故选择GGA-91鹰势文件POTCARo3.截断能ENCUT参数的选择准备文件KPOINTS、POSCAR、P0TCAR以及脚本文件encut并放入同一个文件夹中encut:#!/bin/shrm WAVECARfor i in 550 600 650 700 750 800 850 900docat > INCAR «1System = H2ISTART=0ISNEAR=0, SIGMAR. 1ICHARG = 2ENCUT=$iPREC=Accu

5、rateecho "ENCUT = $i eV” ; vaspE=' grep "TOTEN" OUTCAR | tail 1 | awk ' (printf "%12.6fn”, $5'、KP= grep 'irreducible" CXJTCAR | tail -1 | awk ' (printf '%5i n”, $ 2广echo $i $KP $E » coirmentdoneKPOINTS：Au I onia I ic mesh 0M1110 00POSCAR：F21.2094

6、. 000. 000.000. 004. 000.000. 000. 004.002Direct0. 000. 000. 000. 000. 000. 25选用上次计算是最低能量的POTCAR运行脚本文件encut»得到comment文件内容如卜:550 1 -3. 654087600 1 -3. 654701650 1 -3.651275700 1 -3. 652794750 1 -3. 652555800 1 -3. 652600850 1 -3. 653120 900 1 -3. 652554 950 1 -3. 652669总能变化在0. OOleV左右就足够了.因此我们选择

7、encut=700eV4. KPOINTS参数选择准备文件POSCA、INCAR. POTCAR以及脚本文件k, POSCAR、POTCAR与上面相同 脚本文件k： #!/bin/sh rm WAVECAR for i in 3 5 7 9 11 13 15 do cat > KPOINTS «! Automatic mesh 0 M $i $i $i 000I echo k mesh 二 $i x $i x $i ; vaspE= grep “TOTEN" OUTCAR | tail 1 | awk ' (printf "%12.6f n”, $5

8、 9 KP= grep "irreducible" OUTCAR | tail T | awk ' (printf "%5i n”,$2 K'echo $i $KP $E »conment done INCAR： System = 02 ENCUT = 700 ISTART = 0 ICHARG = 2 ISNEAR = 0 SIGMA = 0. 1 PREC = Accurate运行脚本文件k,得到comment文件：3 6 -3. 6523825 18 -3. 6524117 40 -3. 6523889 75 -3. 652386总

9、能变化在O.OOleV左右就足够了，因此我们选择5作为k点值。5.对晶格常数进行优化准备文件INCAR、KPOINTS. POTCAR以及脚本文件optimize,其余文件与上样Optimize：#! /bin/shB IN=/ usr /b in/vasprm -rf WAVECARfor i in 1. 16 1.21 1. 26 1.31 1. 36 1.41 1. 46 1.51 1. 56 1.61 1. 66; docat >POSCAR «102$i4. 000.000.000. 004.000.000. 000.004.002Direct0. 000. 000.

10、 000. 000. 000. 25I .echo /zd= $i,$BINE= tail -1 OSZICAR'echo $i $E »RESULT. feedonecat RESULT, fcccat RESULT, fee运行脚本文件optimize,得到RESULT文件如F：1. 16 1F 二15006150E+01 E0=1500615(+01 d E179658E-311.21 1F=-.24084278E+01 E0=-.24084278E+01 d E-.145599ET61. 26 1F=30094870E+01 E0=3009-1870E+01 d E2

11、47133E-211.31 1F=33821346E+O1 E0=33821346E+01 d E-403048E-161. 36 1F=35847561E+01 E0=35847564E+01 d E800387E-691.41 1F=36607958E+01 E0=36607958E+01 d E270236E-211. 46 1F=36424044E+01 E0=36424O44E+O1 d E-.866906E-241.51 1F=35553122E+O1 E0=-35553122E+01 d E289397E-271. 56 1F=34180691E+01 E0=34180694E+

12、01 d E998981E-161. 61 1F=-32473 U2E+01 E0=32473112E+01 d E696853E-2361. 66 1 F= - 30547759E+01 E0= - 30547759E+01 d E 540502ET5利用计算数据拟合二次曲线最低点就是计算出的晶格常数。拟 合 后 的 曲 线 如 下 图计算出晶格常数是1. 41实骑室是同样是1.41故参数选择合话计算准确。二、Si晶体晶格常数计算1. 鹰势类型选择分别带入五个鹰势麻文件带入计算去能量最小值的POTCAR输入文件:INCAR:System = SiI START = 0ENCUT = 240N

13、ELM =100EDIFF = 1E-O4EDIFFG = 1E-O3NSW = 5IBRION =2KPOINTS:Automatic mesh0Monkhorst Pack999000POSCAR:Si5.430.000.500.500.500.000.500.500.500.00Direct0.000.000.000.250.250.25计算结果为：质势类型能量GGA-91E0= -.10765021E+02GGA-CAE0=-.11947589E+02PAWE0= -.11895622E+02PAW-GGAE0=-.10872161E+02PBEE0-108921325E+02GGA-

14、CA鹰势文件POTCAR计算能昂:最低所以o E0=-.11947589E+02eV即优化后Si晶体的总能量。2.截断能（ENCUT）参数的选定利用encut脚本文件计算ENCUT参数，准备POSCAR,POTCAR,KPOINT文件。双击encut进行计算。结果文件在comment文件中。选取变化在0.001的ENCUT虬 所用文件：POSCAR：Si5.430.000.500.500.500.000.500.500.500.002Direct0.000.000.000.250.250.25KPOINTS：Automatic mesh 0Monkhorst Pack999000POTCAR从

15、腰势库GGA-CA里选择所对应的文件脚本文件encut：#!/bin/shrm WAVECARfor iin 150 200 250 300 350docat> INCAR «!System = SiENCUT = $iI START = 0ICHARG = 2ISMEAR = 0SIGMA = 0.1PREC = AccurateIecho "ENCUT = $i eV" vaspE= grep -TOTEN- OUTCAR | tail -1 | awk 'printf ”12.6fn”,$5'、KP='grep "ir

16、reducible" OUTCAR | tail -1 | awk 'printf "%5i n", $2 )'' echo $i $KP $E » comment11done运行脚本文件encut,计算结果位于comment文件，内容如E：150 35 41.900756200 35 -11.937845250 35 -11.945420300 35 41.945545350 35 41.945502选择总能变化在O.OOleV左右就足够了。因此，选择ENCUT=250eV3. KPOINTS参数的选定利用kpoints脚本文件

17、优化KPOINTS文件，准备POSCAR,POTCAR,INCAR文件。双击kpoints 进行计算。结果文件在comment文件中°选取变化在0.001的kpoints值。所用文件： #!/bin/shrm WAVE CARfor i in 5 7 9 1113 15 docat> KPOINTS «!Automatic mesh 0M$i $i $i 000Iecho Hk mesh = $ix $ix $i" ; vaspE= grep ,TOTENt, OUTCAR | tail -1 | awk 'printf ”12.6f n”，$5 &

18、#39;、 KP=、grep “irreducible” OUTCAR | tail -1 | awk '(printf "%5i n", $2 '、echo $i $KP $E »commentdone运行脚本文件k,计算后得到comment文件，内容如下：5 10 -11.8856017 20 -11.9352109 35 -11.94320111 56 -11.95798013 84 41.95865415 120 -11.95895总能变化在O.OOleV左右就足够了。因此，选择k点数目为11。4. SIGMA参数的选定利用sigma脚本文

19、件优化参数sigma,准备POSCAR,POTCAR,KPOINTS文件。双击sigma进行计算。结果文件在comment文件中。选取变化在0.001的sigma值。所用文件：#!/bin/shrm WAVE CARfor iin 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 030docat> INCAR «!System = SiENCUT = 250I START = 0ICHARG = 2ISMEAR = 0SIGMA = $iPREC = Accurate echo "SIGMA = $i eV'

20、1; vaspTS= grep "EENTRO" OUTCAR | tail -1 | awk '(printf “12.6f n”，$5 '、KP='grep 'Irreducible" OUTCAR | tail -1 | awk '(printf "%5i n”，$2 '、echo $i $TS »commentdone计算结果在comment文件中，内容如卜,：0.10 0.0000010.12 -0.0000060.14 -0.0000260.16 -0.0000660.18 -0.00

21、01380.20 -0.0002520.22 -0.0004030.24 -0.0006500.26 -0.0009200.28 -0.001420030 -0.001892计算结果在comment文件，选择SIGMA=0.L5.晶格常数计算结果利用optimize脚本文件优化晶格常数，准备INCAR,POTCAR,KPOINTS文件。双击optimize 进行计算。结果在RESULT.fcc文件中。利用计算数据绘图，图上最低点对应晶格常数：4.9 1 F= -.10663223E+02 E0= -.10663203E+02 d E =-.405173E-045.0 1 F= -.111610

22、17E+02 E0= -.11161004E+02 d E =-.254330E-045.11 F= -.11516661E+02 E0= -.11516661E+02 d E =-.519693E-125.2 1 F= -.11750202E+02 E0= -.11750202E+02 d E =.764632E3053 1 F= -.11879566E+02 E0= -.11879566E+02 d E =-347477E-185.4 1 F= -.11920834E+02 E0= -.11920834E+02 d E =-391193E-615.5 1 F= .11888204E+02

23、E0= .11888204E+02 d E 780930E 465.6 1 F= -.11794236E+02 E0= -.11794236E+02 d E =-.297252E-305.7 1 F= -.11649964E+02 E0= -.11649964E+02 d E =-.366560E-245.8 1 F= -.11465105E+02 E0= -.11465105E+02 d E =-.277191E-225.9 1 F= -.11248137E+02 E0= -.11248137E+02 d E =-.115412E46三、Si元素单原子能量计算1. 由内聚能倒推单原子能量Si

24、的内聚能实验值Ecoh=4.63,计算得到的= 11.92,根据公式£妙=一一 2%）可 计算出Si的单原子能量鸟产-1.3423425。2. 基本文件INCAR、POSCAR、KPOINTS. POTCAR153. 单原子能量计算输入文件如下:INCAR：SYSTEN=SiISTART=OICHARG=2ISMEAR=OSIGMA=0.1ENCUT=250PREC=AccurateISPIN=2POSCAR：15.00.000.000.0014.00.000.000.0013.01Cart0.000.000.00KPOINTS：Automatic mesh 0Monkhorst P

25、ack111把优化的POTCAR INCAR POSCAR KPOINTS放入一个文件夹内，运行vasp,得到计算结果Si单原子能鼠E0=.84449503E+00°Si的VASP力学常数计算1. 计算所需文件defvect.f:这个文件使用FORTRAN语言编写经过编译后得到defvectx文件，它是 被。ptimize文件调用的子程序。这个文件主要进行的功能，对应变的Voigt标记的形式进行 定义，变换基矢，以及生成VASP计算所需的POSCAR文件的数据。对应三个不同参数 C11,C12,C44有三个不同的defvectf <» ifort是Intel公司制作的

26、在linux平台使用的编译软件 其中包含C语言与FORTRAN语言的编译器和MKL函数库编译defvectf的具体过程是：在linux系统卜打开"终端”，在其中输入命令ifort -o defvector.x defvector.foOLDPOS：与POSCAR文件类似包含计算晶体初始的基矢与原子坐标。注意：在第一 行要加上一个数字，数字的值等于计算元素的个数。POTCAR:计算晶体的鹰势文件。KPOINTS:对倒空间K点的选择。optimize:主要功能是每次循环中生成vasp计算所需incar,poscar文件并调用 VASP计算。2. 计算与数据处理新建三个文件夹分别取名C11

27、X12X44,把计算所需的五个文件分别放入其中注意不能 把对应的defvectf搞错。然后分别运行其中的optimize就可执行计算。计算结束后把得到的三个SUMMARY文件记录了系统能最E和相应的应变。，将能最与 当。=0时的能虽想减，得到AE,然后拟合出一条的二次曲线。文件导入Excel绘制二次曲线并拟合二次方程。得到二次项系数AO, Al, A2.A0=72.9651.09E+01 y = 72.965x2 0.0713x 10.857 R2 = 0.9983Al=108.1y = 108.1X2 1.5341X 10.857 R2 = 0.9997A2=9.58-1.09E+01y =

28、 9.5825x2 6E-12x -10.857 R2 = 1在0UTCAR文件中查找volume of cell得到参数VO。代入卜式即可得到力学常数。6%9V0219BG3B-2GE = LL =3B+G 6B + 2G3. 计算所用到的公式：a、剪切模量G：E1G =n G = _(Ci - cp)2(1+/)2 11 Lb、体积模量B；E1B = Ml 、=>B=Z(cil + 2c12) 3(1-2/)3c、弹性模量E：G E2("A = -12G- E EGB =3(1-2/Z) 9G-3E 9BGE =3B + Gd、泊松比|J :6B + 2G五、SI晶体的电子

29、结构1.采用VASP计算能带的步骤第一步计算电荷密度准备INCAR KPOINTS POSCAR POTCAR四个文件。运行vasp计算结束后得到计算结果CAGCARo注意INCAR 文件中要包含参数 ICHARG = 2 ; LCHARG = .TURE.第二步能带结构计算：准备INCAR KPOINTS POSCAR POTCARCAGCAR（其中CAGCAR为上次计算的结果）和三个 脚本文件 sym gk.f pbnd.fo1. 首先从上次计算的结果文件0UTCAR中搜索关键词E-Fermi （费米能级） 和direct lattice vectors （倒格矢）把对应的数字分别填入sy

30、ml文件与gk. f文件的 对应位置。2. 在终端中输入ifort-ogk.xgk.f生成gkx文件：双击运行gkx,得到新 的KPOINTS文件。在终端输入vasp如果计算无错误得到结果文件EIGENVAU3. 在终端输入ifort - o pbnd. x pbnd. f 生成pbnd.x文件，双击运行 pbnd x得到所需计算结果文件bnd. dat; highk. date注意INCAR文件中要包含参数ISTART = 1 ICHARG = 11 LORBIT = 11。POSCAR是上次计算 的结果文件CONTCAR重命名后得到的。第三步态密度计算：准备 INCAR KPOINTS P

31、OSCAR POTCARCAGCAR split_dos 在终端中运行 vasp 得到计算 结果EIGENVALo双击运行split_dos分割EIGENVAL文件得到三个结果文件:DOSO DOS1DOS2.注意INCAR 文件中应该包含参数 ICHARG = 11 ISMEAR = - 5 LORBIT = 11 . KPOINTSPOSCAR是第一步计算所用的。2. 电荷分布计算结果能带计算和结果xr LSi的电子能带结构图3. 态密度计算和结果总态密度图OSOQE-Ef(eV)0.20 -0.15-0.10 -0.05 -0.00-0. 50000000000000000. 00000

32、000000000000. 500000000000000026I15I L II L-15-10-50510E-Ef(eV)第一个原子的分波态密度图六、Si晶体介电函数和光学性质的计 算i.计算步骤第一步DFT自洽能带计算：准备文件 POSCAR、INCAR、POTCAR、KP0INTS。运行 VASP。第二步HSE06计算：此步在DFT基础上计算，用到上步的波函数(WAVECAR)把上步计算得到的WAVECAR 拷贝到此计算文件夹下，其他文件从DPT计算文件夹拷贝过来。运行VASPo 第三步optic计算：这步用到第一步计算的波函数，其他全部文件从HSE06计算文件央中拷贝(包括WAVEC

33、AR),次步的KP0INTS文件中k点的个数应为HSE06计算时的2-3倍，关于这一点加密K点是为保证介电谱收敛。运行VASPo2. 用到的文件第一步用到文件：INCARSystem = diamond SiENCUT = 400ISTART = 0ICHARG = 2ISMEAR = 0SIGMA = 0. 1PREC = AccurateEDIFF = 1. 0E-05EDIFFG = -0. 001NSW = 0IBRION = -1ISIF = 2GGA = PS （此处如果用GGA = PS请在POTCAR中修改LEXCH = PS）KPOINTSAutomatic mesh0Mr*

34、7rill000POSCARSi-Diamond:5. 431000000000000. 00000000000000000. 50000000000000000. 50000000000000000. 50000000000000000. 50000000000000900. 0000000000000000Si2Selective dynamicsDirect0.0000000000000000 T0.2500000000000000 T0. 00000000000000000. 0000000000000000TT0. 25000000000000000. 250000000000000

35、0TT第二步用到文件INCARSYSTEM = Si-DiamondISTART = 1ICHARG = 2EDIFF = 0. 00001EDIFFG = -0. 001GGA = PSISMEAR = 0SIGMA = 0. 1ENCUT = 400ENAUG = 800LREAL = . FALSE.LWAVE =, TRUE.LCHARG = . TRUE.NELM = 200NSW = 0IBRION = -1LMA：MIX = 4NSIM = 4ISYM = 0LHFCALC = . TRUE.HFSCREEN = 0.2ALGO = DampedTIME = 0.4ENCUTFO

36、CK = 0AEXX = 0. 25第三步。Ptic计算用到文件：INCARSYSTEM = Si-DiamondISTART = 1ICHARG = 2EDIFF = 0. 00001EDIFFG = -0. 01GGA = PSENCUT = 400ENAUG = 800NELM = 200NSW = 0IBRI ON = 1LMA：MIX = 4ISMEAR = -5NSIM = 4ISYM = 1 （计算光学时：cubic晶型取1, hep晶型取2,正交晶型取3,其他取6）LHFCALC =, TRUE.HFSCREEN = 0.2ALGO = NTIME = 0. 4ENCUTFOC

37、K - 0AEXX = 0. 2530NELMIN = 5LOPTICS = TRUE （计算光学介电矩阵）CSHIFT = 0. 1NBANDS = 24 （计算能带数量，设置为默认值的2飞倍）NPAR = 1 （按能带依次计算）NEDOS = 2000 （取 2000 个计算点）LRPA =, FALSE.LREAL = . FALSE.LWAVE = . FALSE.LCHARG = . FALSE.KPOINTS：Automatic mesh0M9990003. 计算结果计算结果输出在文件OUTCAR中，搜索关键词frequency dependent IMAGINARY DIELEC

38、TRIC FUNCTION 与 frequency dependent REAL DIELECTRIC FUNCTION, 关键字下有列数据，第一列是能量之后三列对应xx, yy, zzo用excel处理数据， 得出介电函数虚部实部，折射系数和消光系数。以下是计算用到的公式：_ (XX + YY + ZZ)2 _ (XX + YY + ZZ32折射率:n =零2但12+622-1消光系数：k =吸收系数:2 +或一勺）反射系数:r = (”T)2+*2(71 + 1)2+A 2作图如下:1.介电常数虚部£2与能量的关系hw(eV)2 介电常数实部£1与能量的关系he(eV)

39、3.材料光学常数photon energy hw(eV)0.7 nhw(eV)hw(eV)七、VASP的声子谱计算i.计算步骤第一步：准备 INPHON. POSCAR. INCAR （可以不用）、POTCAR. KP0INTS 文件第二步：使用超胞进行VASP计算1）分别复制 POSCAR-OOX 到 POSCARo （本例只 ft POSCAR-OOl）382）进行VASP运算,将输出文件vasprun.xml保存下来。（本例只生成vasprun.xml-001）第三步：执行命令：phonopy -f vasprun.xml-001就会得到文件FORCE_SETS、（这步POSCAR必须为

40、最初的POSCAR）第四步：band.conf计算声子谱。第五步：计算热学性质。执行命令:phonopy -p -t -tmax=1500 （计算的最高温度）“tstep=20（温度步长）得到热学性质，保存在Thermal_properties.yaml中：第六步：计算声子态密度执行命令：phonopy -p mesh.conf得到声子态密度，保存在total_dos.dat中。第七步：计算声子分波态密度执行命令：phonopy -p pdos.conf得到声子分波态密度，保存在partial_dos.dat中2.基本文件INPHONATOM NAME = SiDIM = 2 2 2LSUPE

41、R = .TRUE.POSCARSystem Si5.460.500.000.000.500.500.000.00INCARI H UVM君 SIORvo/VHDq,2nslI，N §2 aOH MSN s.o-nofe&QM WSTW TVSm s.«.onfe&QWI H NIdSI z ho-hvhui 0 nhVLSI 一S n WHJLSASccso 二 UIXUEdsKA luTxusdsccA & 见 dsccAEoAAcc>.30AAdsccA avusod50&<usod &ops.sKOJqsaw 一#

42、 必 do2#s0 0 06 6 6W 0 qsom MmuosvsHNIOd?!7Eri=soad8 8 8 = dWzzz = vaa!S = HKVN- HOJV:juorrsopd9 0 9 0 rO = JjnHSdH2 8 2 = dPV zzz = nia!S = HKVN- HOJV juoj*qsuii? = s±Niod_aNvaso so so oooooo我9097 0 SZ9 0 9 0 0 0 9 0 SE O S7 0 9 0 9 0 0 0 9 0 0 0 0 0 0 0 = QNVa zzz =祝。!S = HWVN-攻 OJV申邮 juoopiiE

43、q3UOp3.声子谱、声子态密度计算和结果3468101214Frequency0 8 6 41 o o oS8SS -0zrGUopraErud2 o o.QFrequencySEES JO 3_suoo1636_2««1Wave vector4.热学性质计算和结果八、化合物功2键长计算1. 计算步骤第一步选定合适的potcar：准备文件INCAR、POTCAR、KPOINTS和五个potcar文件。要把碳元素与氧元 素的potcar文件复制在一起，两个元素应选用相同的potcar文件。分别用五个 potcar文件计算，选取能量最小的potcar。结果文件outcar中的

44、E0代表能昂:。第二步选择合适的截断能：POSCAR,KPOINTS文件的内容同上，POTCAR选择PAWCA,将INCAR文件删 除，并将其内容粘贴到脚本程序encut中，改其中的ENCUT=$i即可。双击运行 文件。结果在comment文件中能量变化在0.001的值就做够了。第三步选择合适的k点：准备基本文件INCAR、POSCAR以及脚本文件k和上步选定的potcar。双击 运行脚本文件k。结果在comment文件中，选取能量变化在0.001范围内的k点 值。第四步计算合适的SIGMA数值：准备基本文件POTCAR. POSCAR、KPOINT以及脚本文件SIGMA。双击运行 即可“第五

45、步计算键长：准备基本文件INCAR、POTCAR. KPOINT以及脚本文件optimize,双击运行 即可。能量最低点对应键长。2. 基本文件8的POTCAR的选择INCARSystem = CO2IS1ART = O;ICHARG = 2ENCUT= 400ISMEAR = 0；SIGMA = 0.1EDIFF = 1E-O4EDIFFG=1E-O3IBRION = 2PREC = Accurate39KPOINTSAutomatic generation0990.0 0.090.0CO24.641.00000.00000.00000.00001.00000.00000.00000.000

46、01.000012Direct0.0 0.0 0.00.0 0.25 0.00.0 0.0 0.5Monhkorst-PackPOSCARco2的截断能优化ENCUT削/bin/shrm WAVECARfor iin 250 300 350 400 450 500 550 600 650docat > INCAR «!System = C02IS7ART = 0;ICHARG = 2ENCUT= $iISMEAR = 0;SIGMA = 0.1EDIFF = 1E-O4EDIFFG = 1E-O3IBRI0N = 2PREC = AccurateI .echo "EN

47、CUT = $i eV" ; vaspE='grep "TOTEN" OUTCAR | tail -1 | awk 'printf ”12.6fn”,$5、KP= grep "irreducible" OUTCAR | tail -1 | awk 'printf ”5i n”，$2 ' echo $i $KP $E » comment doneKPOINTSAutomatic generationMonhkorst-Pack0.0 0.0POSCAR1.00000.00000.00000.00001.00000.00000.00000.00001.0000Direct 0.0 0.00.0 0.0 0.5 0.667 0.333 0.5CQ的SIGMA的优化POSCARcoz4.641.00000.00000.00000.00001.00000.00000.00000.00001.000012Direct0.0 0.0 0.00.0 0.25 0.00.0 0.0 0.5KPOINTSAutomatic generation0Monlikorsl-Pdck9990.0 0.0 0.0其 POTCAR 内容为 PAWGG