vasp经验总结.ppt
《vasp经验总结.ppt》由会员分享,可在线阅读,更多相关《vasp经验总结.ppt(24页珍藏版)》请在咨信网上搜索。
1、vasp计算算总结1.加加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 eff
2、ective on-site Coulomb interaction parameters#specifies the effective on-site Exchange interaction parameters.2.ISTART=1ICHARG=11ISMEAR=-5LORBIT=11用用VASP计算算DOS准备好KPOINTS文件,增加k点网格将上一步自洽计算得到的CHG、CHGCAR拷贝至同一目录下3.用用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为参考的本征值。4.简立方的高对称点5.60 60 1 1 0.1158924E+02 0.5
4、522500E-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点对应的本征值的序号(即第几条能带),及相应的本征值。能量本征能量本征值文件文件EIGENVAL6.POSCAR.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 displ
6、ay.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找最低能量的晶格参数找最低能量的晶格参数C7.VASP中,用中,用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个中心对称点群是没有极化的。点电荷估计:晶格某一方向所有原子坐标*离子价态相加,与优化后相加得到的值对比。8.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.把线性计算的
9、 IBZKPT文件产生的K点copy成KPONITS文件,作为soc计算的输入文件9.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
10、not 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值来解决。10.2、ALGO,IALGO,LDIAG If the self-consistency loop does not converge within 40
11、 steps,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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- vasp 经验总结
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,个别因单元格分列造成显示页码不一将协商解决,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【快乐****生活】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【快乐****生活】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。