第一性原理算例.docxVIP

首页| SCI期刊点评|中文期刊点评| 搜索 |论文范文|小木虫论坛 学末科研 学末科研?-站 i/i/wyuXhit￡C h. net 首页 第一原理 查看话题 vasp5.2中GGA=PE, PS, B3计算晶格常数，能带（HSE06杂化） 作者：liliangfang （站内联系 TA）发布：2011-10-20 在vasp5.2中GGA=PE，PS，B3计算晶格常数，能带（HSE06杂化）的对比 下面是其他两个，基于PAW-PBE的赝势，默认的是GGA=PE 在INCAR里要设置GGA=PS，相应的在POTCAR里LEXCH = PS； 在INCAR里要设置GGA=B3，相应的在POTCAR里LEXCH =B3 晶格常数 计算值 误差（实验值5.43） PE 5.472 0.768% PS 5.438 0.147% B3 5.445 0.275% 能带计算（HSE06）Eg （实验值1.12） TOC \o 1-5 \h \z PE 1.26 PS 1.28 B3 1.60 PE 0.61（4.6 的方法） GGA=PE o X— o-5 Aj①u山 GGA=PS o X— o-5 AQOU 山 GGA=B3 o X— o-5 A6JCDU 山 以前4.6的算法得到的 o X— o-5 AJOU 山 下面将计算晶格常数及能带各部的主要文件上传（个人计算总结） 晶格常数计算 INCAR System = Si ENCUT = 400 ISTART = 0 ICHARG = 2 PREC = Accurate EDIFF = 1.0E-05 EDIFFG = -0.01 NSW = 200 IBRION = 2 ISIF = 4 ISPIN = 2 GGA = PS KOPINTS Automatic mesh 0 M POSCAR Si: 5.43 0.0000000000000000 0.5000000000000000 0.5000000000000000 2 0.5000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 S Direct 0.000000000000000 0.000000000000000 0.000000000000000 T T T 0.250000000000000 0.250000000000000 0.250000000000000 T T T POTCAR 优化脚步jvasp #! /bin/bash rm -rf WAVECAR for i in 5.437 5.438 5.439 5.440 5.441 do cat POSCAR Si-Diamond: $i0.00000000000000000.5000000000000000 $i 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 S Direct 0.000000000000000 0.000000000000000 0.000000000000000 T T T 0.250000000000000 0.250000000000000 0.250000000000000 T T T ! echo a= $i ;mpdboot ; mpirun -np 4 vasp E=tail -1 OSZICAR; echo $i $E SUMMARY cp CONTCAR CONTCAR$i done 二次拟合出能量最低点，一般两次足够 能带计算 vasp5.2其中一个计算能带方法是，优化晶格常数后，DFT+HSE06+band计算，见附件 这是其中一个PS的步骤，注意选择POTCAR之后将其中的LEXCH改为LEXCH = PS 1.DFT INCAR SYSTEM = Si ENCUT = 400 ISTART = 0 ISYM = 0 GGA = PS ICHARG = 2 ISMEAR = 0 SIGMA = 0.05 NSW = 0 IBRION = -1 ISIF = 2 EDIFF = 0.00001 EDIFFG = -0.001 KPOINTS Automatic mesh POSCAR Si-Diamond: 5.438000000000000.00000