材料计算vasp 程序

《材料计算vasp 程序》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、用VASP计算Ｈ原子的能量氢原子的能量为。在这一节中，我们用VASP计算H原子的能量。对于原子计算，我们可以采用如下的INCAR文件PREC=ACCURATE:NELMDL=5makefivedelayscatillchargemixingISMEAR=0;SIGMA=0.05usesmearingmethod采用如下的KPOINTS文件。由于增加K点的数目只能改进描述原子间的相互作用，而在单原子计算中并不需要。所以我们只需要一个K点。MonkhorstPack0MonkhorstPack111000采用如下的POSCAR文件atom115.00000.0

2、0000.00000.0000015.00000.00000.00000.0000015.000001cart000采用标准的H的POTCAR得到结果如下：k-point1:0.00000.00000.0000bandNo.bandenergiesoccupation1-6.31451.000002-0.05270.0000030.48290.0000040.48290.00000我们可以看到，电子的能级不为。Freeenergyoftheion-electronsystem(eV)-------------------------------------

3、--------------alphaZPSCENC=0.00060791EwaldenergyTEWEN=-1.36188267-1/2HartreeDENC=-6.27429270-V(xc)+E(xc)XCENC=1.90099128PAWdoublecounting=0.000000000.00000000entropyT*SEENTRO=-0.02820948eigenvaluesEBANDS=-6.31447362atomicenergyEATOM=12.04670449-------------------------------------

4、--------------freeenergyTOTEN=-0.03055478eVenergywithoutentropy=-0.00234530energy(sigma->0)=-0.01645004我们可以看到也不等于。在上面的计算中有个问题，就是H原子有spin,而在上面的计算中我们并没有考虑到spin。所以如果我们改用LSDA近似，在INCAR中用ISPIN=2的tag，则得到如下结果：k-point1:0.00000.00000.0000bandNo.bandenergiesoccupation1-7.27361.000002-0.12290

5、.0000030.45620.0000040.45620.0000050.45620.00000spincomponent2k-point1:0.00000.00000.0000bandNo.bandenergiesoccupation1-2.41400.000002-0.07010.0000030.51790.0000040.51790.0000050.51790.00000Freeenergyoftheion-electronsystem(eV)---------------------------------------------------alp

6、haZPSCENC=0.00060791EwaldenergyTEWEN=-1.36188267-1/2HartreeDENC=-6.68322940-V(xc)+E(xc)XCENC=2.38615430PAWdoublecounting=0.000000000.00000000entropyT*SEENTRO=0.00000000eigenvaluesEBANDS=-7.27361676atomicenergyEATOM=12.04670449---------------------------------------------------free

7、energyTOTEN=-0.88526212eVenergywithoutentropy=-0.88526212energy(sigma->0)=-0.88526212氢原子的能量约等于。可以看到在LDA中如果限制自旋，使能级大概提高了。但是如何理解所得到的能级，由于用到了赝势，本人并不很清楚如何解释能级意义。用VASP计算Pd金属的晶格常数Pd金属的实验上的晶格常数为。在这里，我们用VASP计算它的晶格常数。首先将Pd所对应的POTCAR文件拷贝到目录下。然后准备好INCAR和KPOINTS文件。POSCAR文件我们将通过一个tcsh的script来产

8、生。KPOINTS文件可以如下：MonkhorstPack0Mon