VASP(计算前的各种测试)

《VASP(计算前的各种测试)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、（计算前的）验证一、检验赝势的好坏：（一）方法：对单个原子进行计算；（二）要求：1、对称性和自旋极化均采用默认值；2、ENCUT要足够大；3、原胞的大小要足够大，一般设置为15Å足矣，对某些元素还可以取得更小一些。（三）以计算单个Fe原子为例：1、INCAR文件：SYSTEM=FeatomENCUT=450.00eVNELMDL=5!makefivedelaystillchargemixing，详细意义见注释一ISMEAR=0SIGMA=0.12、POSCAR文件：atom15.001.000.000.000.001.000.000.000.001.001Direct0003

2、、KPOINTS文件：（详细解释见注释二。）Automatic0Gamma1110004、POTCAR文件：（略）注释一：关键词“NELMDL”:A）此关键词的用途：指定计算开始时电子非自洽迭代的步数（即NELMDLgivesthenumberofnon-selfconsistentstepsatthebeginning），目的是makecalculationsfaster。“非自洽”指的是保持chargedensity不变，由于ChargedensityisusedtosetuptheHamiltonian,所以“非自洽”也指保持初始的哈密顿量不变。B）默认值（defaul

3、tvalue）:NELMDL=-5(当ISTART=0,INIWAV=1,andIALGO=8时)  NELMDL=-12(当ISTART=0,INIWAV=1,andIALGO=48时)  NELMDL=0(其他情况下)NELMDLmightbepositiveornegative.Apositivenumbermeansthatadelayisappliedaftereachionicmovement--ingeneralnotaconvenientoption.（在每次核运动之后）Anegativevalueresultsinadelayonlyforthestart-

4、configuration.（只在第一步核运动之前）C）关键词“NELMDL”为什么可以减少计算所需的时间？ChargedensityisusedtosetuptheHamiltonian,thenthewavefunctionsareoptimizediterativelysothattheygetclosertotheexactwavefunctionsofthisHamiltonian.Fromtheoptimizedwavefunctionsanewchargedensityiscalculated,whichisthenmixedwiththeoldinput-ch

5、argedensity.Abriefflowchartisgivenbelow.（参自ManualP105页）一般情况下，theinitialguessedwavefunctions是比较离谱的，在前NELMDL次非自洽迭代过程中保持chargedensity不变、保持初始的哈密顿量不变，只对wavefunctions进行优化，在得到一个与theexactwavefunctionsofinitialHamiltonian较为接近的wavefunctions后，再开始同时优化chargedensity。这样一来，计算时间要比一开始就同时优化chargedensity和wavef

6、unctions短得多。注释二：为什么这里只需要一个k点？Foratomsandmolecules,theBlochtheoremdoesnotapply,hencethereisnoneedtousemorethanonesinglek-point.Whenmorek-pointsareused,onlytheinteractionbetweentheatoms(whichshouldbezero)isdescribedmoreaccurately.（三）计算任务执行方法：输入：vasp（四）赝势好的判断标准：计算得到的OUTCAR文件中的“energywithoutent

7、ropy”能量值在-0.001～-0.01eV之间。二、筛选合适的ENCUT大小：（一）输入文件：1、用脚本程序optencut.sh代替INCAR文件：rmWAVECARforiin100150200250300350400450500550600docat>INCAR<

8、tail-1

9、a