VASP單機(jī)操作流程以Ni111表面為例電子版本

1、Good is good, but better carries it.精益求精，善益求善。VASP單機(jī)操作流程以Ni111表面為例-VASP單機(jī)操作流程（以Ni111表面為例）將生成的POSCAR文件和VASP自帶的相關(guān)元素的贗勢(shì)文件POTCAR導(dǎo)入或者是復(fù)制到工作目錄。打開(kāi)終端，編輯K點(diǎn)文件KPOINTS（在此使用自動(dòng)生成的方法）。CatKPOINTSA0771000!保存退出。二，優(yōu)化截?cái)嗄埽簝?yōu)化過(guò)程的腳本為：#!/bash/shforiin200250300350400450500550600docatINCAREtot.datdone保存退出。輸入命令：shEN.sh點(diǎn)擊Enter鍵

2、運(yùn)行VASP。進(jìn)行截?cái)嗄艿膬?yōu)化。（進(jìn)行并行VASP運(yùn)算前要先輸入命令：mpd&;VASP的并行命令為：mpirunnp2vaspmpi）使用命令：catEtot.dat來(lái)顯示計(jì)算結(jié)果。例如：200-35.941564250-36.653335300-36.667632350-36.627027400-36.617354450-36.611490500-36.614625550-36.619412600-36.616586根據(jù)要求，我選擇400eV的截?cái)嗄堋?yōu)化K點(diǎn)：復(fù)制EN.sh，只是為了方便使用它的格式。CpEN.shKY.sh（自定義文件名）優(yōu)化過(guò)程的腳本為：#!/bash/shforii

3、n5791113151719docatKPOINTSEtot.datdone保存退出。輸入命令：shKY.sh點(diǎn)擊Enter鍵運(yùn)行vasp進(jìn)行K點(diǎn)的優(yōu)化。計(jì)算完成后，輸入命令：catEtot.dat查看結(jié)果。如同截?cái)嗄艿膬?yōu)化一樣，根據(jù)要求選擇合適的K點(diǎn)。（在此我選擇17）優(yōu)化晶格常數(shù)：在優(yōu)化晶格常數(shù)時(shí)，是對(duì)比例系數(shù)的優(yōu)化，不能使用絕對(duì)值。一般在試驗(yàn)值的左右各取一些值，根據(jù)要求選擇最合適的晶格常數(shù)。計(jì)算過(guò)程和上兩步類似。CpEn.shPOS.sh（為方便使用格式）計(jì)算用腳本程序?yàn)椋?!/bash/shforiin0.920.940.960.981.001.021.041.061.081.10do

4、c=echoscale=6;12.821/$i|bccatPOSCAREtot.datdone保存退出。輸入命令：shPOS.sh點(diǎn)擊Enter鍵運(yùn)行VASP,進(jìn)行晶格常數(shù)的優(yōu)化。根據(jù)計(jì)算結(jié)果，選擇對(duì)應(yīng)能量值最低的晶格常數(shù)即可。用計(jì)算所得的晶格常數(shù)代替原有POSCAR文件中的晶格常數(shù)。五，遲豫計(jì)算（relax）：編輯INCAR文件（在原有基礎(chǔ)上加上離子步）輸入命令：viINCAR進(jìn)入編輯模式SYSTEM=Ni(111)ENCUT=400ISTART=0ICHARG=2ISMEAR=0SIGMA=0.05PREC=ANELM=100EDIFF=1E-4NGX=18NGY=18NGZ=162NSW

5、=100IBRION=1ISIF=2EDIFFG=-0.02（保存退出）編輯KPOINTS文件（采用優(yōu)化后的K點(diǎn)）輸入命令：viKPOINTS進(jìn)入編輯模式A0G171710.00.00.0（保存退出）編輯POSCAR文件（采用優(yōu)化后的晶格常數(shù)）將CONTACAR文件的內(nèi)容復(fù)制到POSCAR文件中，原因是CONTCAR文件中包含有結(jié)構(gòu)優(yōu)化的信息。輸入命令：cpCONTCARPOSCAR另外根據(jù)需要固定幾個(gè)原子層Ni(111)1.000000000000001.7619999647140503-1.76199996471405030.00000000000000000.00000000000000

6、001.7619999647140503-1.761999964714050312.821502685546875012.821502685546875011.6554540000000006Ni7SelectivedynamicsDirect0.0000000000000007-0.0000000000000003-0.0000000000000001TTT0.0000000000000002-0.00000000000000020.2748507787801908TTT-0.00000000000000020.00000000000000010.5497015575603815TTT0.6

7、6666666666666660.33333333333333280.1832338525201271FFF0.66666666666666640.33333333333333310.4580846313003180TTT0.33333333333333390.66666666666666660.0916169262600635FFF0.33333333333333330.66666666666666660.3664677050402544TTT（其它內(nèi)容省略）(保存退出)建立relax文件夾：mkdirrelax將CHGCHGCARWAVECARINCARPOSCARKPOINTSPOTCA

8、R復(fù)制到該文件夾：CpCHGCHGCARWAVECARINCARPOSCARKPOINTSPOTCARrelax則在relax文件夾下運(yùn)行命令;Mpirun-np2vaspmpi進(jìn)行遲豫計(jì)算靜態(tài)自洽計(jì)算（scf）：建立文件夾scf：輸入命令：Mkdirscf將POTCARKPOINTSINCAR拷貝到scf文件夾中。將CONTCAR改為POSCAR在scf文件夾下對(duì)相關(guān)文件進(jìn)行編輯。編輯INCAR文件，去掉遲豫計(jì)算的參數(shù)設(shè)置（也就是關(guān)閉離子步）。另外還要控制輸出電荷密度和波函數(shù)。輸入命令：viINCARSYSTEM=Ni(111)ENCUT=400ISTART=0ICHARG=2ISMEAR=

9、0SIGMA=0.05PREC=ANELM=100EDIFF=1E-6（提高精度）NGX=18NGY=18NGZ=162LCHARG=.T.(輸出電荷密度)LWAVE=.T.（輸出波函數(shù)）LVTOT=.T.（保存退出）在scf文件夾下輸入命令：mpirun-np2vaspmpi進(jìn)行遲豫計(jì)算。態(tài)密度計(jì)算（dos）：建立dos文件夾，輸入命令：mkdirdos將KPOINTSINCARPOSCARPOTCARCHGCHGCARWAVECAR復(fù)制到dos文件夾中，輸入命令：cpKPOINTSINCARPOSCARPOTCARCHGCHGCARWAVECARdos(INCAR是遲豫計(jì)算后的文件，POS

10、CAR是晶格常數(shù)優(yōu)化計(jì)算后的文件，KPOINTS是K點(diǎn)優(yōu)化計(jì)算后的文件)編輯INCAR文件，輸入命令：viINCAR修改以下兩個(gè)參數(shù)：ISTART=1IVHARG=11增加以下兩個(gè)參數(shù)：LROBIT=11NPARA=1(控制并行計(jì)算)在進(jìn)行態(tài)密度計(jì)算時(shí)，K點(diǎn)可以不用改變。在dos文件夾下，輸入命令：mpirun-np2vaspmpi運(yùn)行vasp進(jìn)行態(tài)密度的計(jì)算。能帶計(jì)算（bands）：建立bands文件夾，輸入命令：mkdirbands將KPOINTSINCARPOSCARPOTCARCHGCHGCARWAVECAR復(fù)制到bands文件夾中，輸入命令：cpKPOINTSINCARPOSCARPOTCARCHGCHGCARWAVECARbands(INCAR是遲豫計(jì)算后的文件，POSCAR是晶格常數(shù)優(yōu)化計(jì)算后的文件)編輯INCAR文件，輸入命令：viINCAR將dos計(jì)算中的LROBIT=11NPARA=1兩項(xiàng)關(guān)閉或者是刪除編輯KPOINTS文件（在計(jì)