vasp经验总结.ppt
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_1.gif)
![资源得分’ title=](/images/score_05.gif)
《vasp经验总结.ppt》由会员分享,可在线阅读,更多相关《vasp经验总结.ppt(24页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、vasp计算总结计算总结加加U(LMCO为例)为例)LDAU=.TRUE.LDAUTYPE=2LDAUL=-1 2 2-1 LDAUU=0 5.5 4 0 LDAUJ=0 1.5 1 0#The simplified(rotationally invariant)approach to the LSDA+U#specifies the l-quantum number for which the on-site interaction is added (-1=no on-site terms added,1=p,2=d,3=f,Default:LDAUL=2)#specifies the e
2、ffective on-site Coulomb interaction parameters#specifies the effective on-site Exchange interaction parameters.ISTART=1ICHARG=11ISMEAR=-5LORBIT=11用用VASP计算计算DOS准备好KPOINTS文件,增加k点网格将上一步自洽计算得到的CHG、CHGCAR拷贝至同一目录下用用VASP计算能带计算能带ISTART=1ICHARG=11NBANDS一般可以用默认,有时候可以适当增加ISMARE=-5 不可用。半导体或绝缘体用0,金属用1NAME20Line
3、-modeRec0.0 0.0 0.00.0 0.5 0.00.0 0.5 0.00.5 0.5 0.00.5 0.5 0.00.5 0.5 0.50.5 0.5 0.50.0 0.0 0.0简立方的高对称点能带处理:程序 band.cpp,g+编译,得到可执行文件a.outvasp计算得到EIGENVAL,删除前七行grep-v E EIGENVAL band.in./a.out 运行,得到bnd000.datbnd000.dat文件的第一列数据是k点距离的绝对值,第二列数据是以Ferim level为参考的本征值。简立方的高对称点 60 60 1 1 0.1158924E+02 0.552
4、2500E-09 0.2280000E-08 0.5522500E-09 0.5000000E-15 1.000000000000000E-004 CAR YFTO 492 32 300 0.3742002E-15 0.1262379E-14 0.3742002E-15 0.1388889E-01 1 -46.511682 2 -46.511280 3 -46.451176 4 -46.449273 5 -46.199455 6 -46.197355第8行的前三个数是k点的坐标,第四个数是相应k点的权重。第一行,前三个整数无意义,第四个整数,如果是2,表示是自旋极化的计算,如果是1,表示非自旋
5、极化的计算。第2至5行的数据含义不大明确,可以不管它。第6行的数据表示:第一个数表示体系总的价电子数目,第二个数表示的计算能带时总的k点数目,第三个数表示的是计算能带时计算了多少条能带。第9行给出的是该k点对应的本征值的序号(即第几条能带),及相应的本征值。能量本征值文件能量本征值文件EIGENVALPOSCAR.shchmod+x POSCAR.sh 改为可执行文件./POSCAR.sh 运行#!/bin/shrm WAVECARfor i in 7.6 7.7 7.8 7.9 8.0 8.1 8.2 8.3 8.4 8.5 8.6(改为要计算的值)docat POSCAR display.
6、logE=grep TOTEN OUTCAR|tail-1|awk printf%12.6f n,$5 mv display.log display.log-$imv OUTCAR OUTCAR-$imv DOSCAR DOSCAR-$iecho$i$E commentdone找最低能量的晶格参数找最低能量的晶格参数CVASP中,用中,用Berry Phase的方法计算极化值的方法计算极化值INCAR中要添加的参数:LCALCPOL=.TRUE.打开计算Berry的开关EDIFF=1E-6 精度DIPOL=0.4 0.4 0.4 选取参考点,任意选的,但是不要和离子重合首先,判断一个体系是否有
7、可能存在极化首先,判断一个体系是否有可能存在极化1.看晶胞所属的群是否是极化群,磁性有影响时要考虑进去。MS查看对称性,找到所对应的群。http:/homepage.univie.ac.at/nikos.pinotsis/spacegroup.html#14 找到对应的点群铁电体物理,钟维烈著,P651。附录一,30个晶体点群。10个极性点群是非常有可能产生极化的,11个非极性中兴对称群是可能有极化的,11个中心对称点群是没有极化的。点电荷估计:晶格某一方向所有原子坐标*离子价态相加,与优化后相加得到的值对比。Spin orbital coupling(SOC)的计算的计算(LMCO为例)为例
8、)ISTART=1 要读取WAVCARICHARG=1 最好选取1而不用11,对结果影响是比较大的LSORBIT=.TRUE.LMAXMIX=6 计算soc的时候一定要有此参数,d电子4,f电子6MAGMOM=6*0 0 0 4 0 0 4 0 0 4 0 0-4 0 0-4 0 0-4 0 0 3 0 0 3 0 0 3 0 0 3 0 0-3 0 0-3 0 0-3 0 0-3 72*0 SAXIS=1 1 1 磁矩的方向(见vasp手册,有两种设置方法)NBANDS=408 能带数是线性计算的二倍ISYM=0 计算soc最好去掉对称性GGA_COMPAT=.FALSE.把线性计算的 IB
9、ZKPT文件产生的K点copy成KPONITS文件,作为soc计算的输入文件VASP难收敛的问题难收敛的问题于一些磁性体系、镧系和锕系元素及相关化合物的静态计算(电子迭代),经常会遇到“难收敛”的问题。1、LMAXMIX Default:LMAXMIX=2 An additional flag controls up to which l quantum number the onsite PAW charge densities are passed through the charge density mixer.Higher l-quantum numbers are usually n
10、ot handled by the mixer.In order to obtain fast convergence to the groundstate,you can try the following setting:LMAXMIX=4 for d elements LMAXMIX=6 for f elements这个FLAG对于含d电子和f电子的体系是非常重要的,很大一部分体系的收敛问题可以通过设置合适的LMAXMIX值来解决。2、ALGO,IALGO,LDIAG If the self-consistency loop does not converge within 40 ste
11、ps,it will probably not converge at all.In this case you should reconsider the tags IALGO,LDIAG,and the mixing-parameters.一般情况下,或使用IALGO=48时遇到收敛问题的话,可以考虑设IALGO为38,或设置ALGO=Normal or Fast(in VAS P.4.5 and later versions)。DefaultALGO=Normal3、NELMDLNELMDL gives the number of non-selfconsistent steps at
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- vasp 经验总结
![提示](https://www.deliwenku.com/images/bang_tan.gif)
限制150内