2015源资培训班-VASP上机练习讲解(几何结构优化-表面、吸附).pdf

《2015源资培训班-VASP上机练习讲解(几何结构优化-表面、吸附).pdf》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、VASP上机练习讲解几何结构优化（表面、吸附）Pd(111)表面吸附CO分子2计算流程：（1）对Pd晶胞进行结构优化（2）采用优化好的Pd晶胞切Pd(111)面（建模）（3）优化Pd(111)表面（4）优化CO分子（5）在优化好的Pd(111)表面上吸附CO分子（建模）（6）优化Pd(111)表面吸附CO分子3（1）对Pd晶胞进行结构优化MedeA平台上建模：晶格常数：a=b=c=3.86700α=β=γ=90（1）对Pd晶胞进行结构优化SYSTEM=PdISTART=0ICHARG=2优化晶胞常数P

2、REC=NormalLREAL=.F.IBRION=2ISIF=3！需要改变Pd的晶胞体积和晶格参数NSW=100POTIM=0.5EDIFFG=-0.05ENCUT=250eV！截断能根据相关文献或贋势确定，也可自己测试NELM=60EDIFF=0.1E-04LCHARG=.F.！优化晶胞时，不打开电荷密度，节省计算时间LWAVE=.F.！优化晶胞时，不打开波函数，节省计算时间ISMEAR=1!金属体系，选择1SIGMA=0.2（1）对Pd晶胞进行结构优化Pdbulk1.003.8670000000

3、00000000.00000000000000000.00000000000000000.00000000000000003.867000000000000000.00000000000000000.00000000000000000.00000000000000003.867000000000000004SelectivedynamicsDirect1.000001.000001.00000FFF输入文件POSCAR1.000000.500000.50000FFF0.500001.000000.50

4、000FFF0.500000.500001.00000FFF提示：优化晶胞时，原胞中原子的相对坐标不变，只优化晶胞的晶格参数，因此需要将每个原子坐标固定设为F，并且ISIF=3（1）对Pd晶胞进行结构优化输入文件POTCAR：（略）输入文件KPOINTSmeshauto0G101010000提示：在做晶胞优化时，需要将晶胞多优化几遍，以获得能量最优的结构。比如：逐级精度加高的方法，因此在选择K点的时候，可以先从6×6×6到8×8×8再到10×10×10，优化三遍晶胞，并注意与实验值比较，获得最准确的晶

5、胞结构。7（1）对Pd晶胞进行结构优化输出文件OUTCAR最终晶胞的体积和晶格参数都会变化。晶胞常数优化为3.98391按住shift+G命令，到OUTCAR文件最末端，看到体系优化后的晶胞能量（1）对Pd晶胞进行结构优化SYSTEM=PdISTART=0ICHARG=2弛豫原子PREC=NormalLREAL=.F.IBRION=2ISIF=2！不改变Pd的晶胞体积和晶格参数NSW=100POTIM=0.5EDIFFG=-0.05ENCUT=250eV！截断能根据相关文献或贋势确定，也可自己测试NE

6、LM=60EDIFF=0.1E-04LCHARG=.F.！优化晶胞时，不打开电荷密度，节省计算时间LWAVE=.F.！优化晶胞时，不打开波函数，节省计算时间ISMEAR=1!金属体系，选择1SIGMA=0.2得到的结构为最终优化得到的结构。（2）创建Pd(111)表面4层，2×2，Pd（111）表面,一共16个Pd原子10（2）创建Pd(111)表面将最底下的2层Pd原子fix住，模拟材料的体相结构11（3）优化Pd(111)表面SYSTEM=Pd-111ISTART=0ICHARG=2输入文件INC

7、ARPREC=NormalLREAL=.F.IBRION=2ISIF=2！优化表面NSW=100POTIM=0.5EDIFFG=-0.05ENCUT=250eV！同Pd晶胞优化采用相同设置NELM=60EDIFF=0.1E-04LCHARG=.F.提示：LWAVE=.F.ISIF参数：晶胞优化取3，表面优化取2ISMEAR=1ENCUT参数：由体系决定，可参考文献及SIGMA=0.2ALGO=FastPOTCAR中的值确定。12（3）优化Pd(111)表面Pd1111.05.634100000.000

8、000000.00000000输入文件POSCAR-2.817049884.879273800.000000000.000000430.0000007423.20050000提示：16•Pd（111）表面是简单单斜结构，在定义lattice的时候，Selectivedynamics第一行为晶矢a分别在X，Y，Z方向上的投影，a与X轴平行，所Direct以为0.000000000.000000000.09914000FFF（5.634100000.0000