第一性原理算例.doc

首页 | SCI期刊点评 | 中文期刊点评 | 搜索 | 论文范文 | 小木虫论坛首页第一原理 查看话题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）PE 1.26PS 1.28B3 1.60PE 0.61（4.6的方法）GGA=PEGGA=PSGGA=B3以前4.6的算法得到的下面将计算晶格常数及能带各部的主要文件上传(个人计算总结）晶格常数计算INCARSystem =SiENCUT = 400ISTART = 0ICHARG = 2PREC = AccurateEDIFF = 1.0E-05EDIFFG = -0.01NSW = 200IBRION = 2ISIF = 4ISPIN = 2GGA = PSKOPINTSAutomatic mesh0M5 5 50 0 0POSCARSi: 5.43 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 2SDirect0.0000000000000000.0000000000000000.000000000000000TTT0.2500000000000000.2500000000000000.250000000000000TTTPOTCAR优化脚步jvasp#! /bin/bashrm -rf WAVECARfor i in5.437 5.438 5.439 5.440 5.441docat POSCAR SUMMARYcp CONTCAR CONTCAR$idone二次拟合出能量最低点，一般两次足够能带计算vasp5.2其中一个计算能带方法是，优化晶格常数后，DFT+HSE06+band计算，见附件这是其中一个PS的步骤,注意选择POTCAR之后将其中的LEXCH 改为LEXCH= PS1.DFTINCARSYSTEM = SiENCUT = 400ISTART = 0ISYM = 0GGA = PSICHARG = 2ISMEAR = 0SIGMA = 0.05NSW = 0IBRION = -1ISIF = 2EDIFF = 0.00001EDIFFG = -0.001KPOINTSAutomatic mesh0G4 4 40 0 0POSCARSi-Diamond: 5.43800000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 0.5000000000000000 0.5000000000000000 0.5000000000000000 0.0000000000000000 Si 2Selective dynamicsDirect0.00000000000000000.00000000000000000.0000000000000000 T T T0.25000000000000000.25000000000000000.2500000000000000 T T T0.00000000E+000.00000000E+000.00000000E+000.00000000E+000.00000000E+000.00000000E+00运行vasp2.HSE06INCAR设置ICHARG =2 ，IBRION =1 ,NELMIN =5 or 其他值如下SYSTEM = SiISTART = 1ICHARG = 2EDIFF = 0.00001EDIFFG = -0.001ENCUT = 400ENAUG = 800GGA = PSLREAL = .FALSE.LWAVE = .TRUE.LCHARG = .TRUE.NELM = 200NSW = 0IBRION = -1LMAXMIX = 4ISMEAR = 0SIGMA = 0.1NSIM = 4ISYM = 0LHFCALC = .TRUE.HFSCREEN = 0.2ALGO = DTIME = 0.4ENCUTFOCK = 0AEXX = 0.25其他文件从DFT文件夹拷贝过来，运行vasp3.band计算INCARSYSTEM = SiISTART = 1ICHARG = 2GGA = PSEDIFF = 0.00001EDIFFG = -0.001ENCUT = 400ENAUG = 800LREAL = .FALSE.LWAVE = .TRUE.LCHARG = .TRUE.NELM = 200NSW = 0IBRION = 1LMAXMIX = 4ISMEAR = 0SIGMA = 0.1NSIM = 4ISYM = 0LHFCALC = .TRUE.HFSCREEN = 0.2ALGO = NTIME = 0.4ENCUTFOCK = 0AEXX = 0.25NELMIN = 5KPOINTS从HSE06文件夹的IBZKPT拷贝到KPOINTS,然后把高对称点拷贝在后面，权重为0更新k点总数Automatically generated mesh 127（36+91后面91为高对称）Reciprocal lattice 0.00000000000000 0.00000000000000 0.00000000000000 1 0.25000000000000 0.00000000000000 0.00000000000000 2 0.50000000000000 0.00000000000000 0.00000000000000 1 0.00000000000000 0.25000000000000 0.00000000000000 2 0.25000000000000 0.25000000000000 0.00000000000000 2 0.50000000000000 0.25000000000000 0.00000000000000 2 -0.25000000000000 0.25000000000000 0.00000000000000 2 0.00000000000000 0.50000000000000 0.00000000000000 1 0.25000000000000 0.50000000000000 0.00000000000000 2 0.50000000000000 0.50000000000000 0.00000000000000 1 0.00000000000000 0.00000000000000 0.25000000000000 2 0.25000000000000 0.00000000000000 0.25000000000000 2 0.50000000000000 0.00000000000000 0.25000000000000 2 -0.25000000000000 0.00000000000000 0.25000000000000 2 0.00000000000000 0.25000000000000 0.25000000000000 2 0.25000000000000 0.25000000000000 0.25000000000000 2 0.50000000000000 0.25000000000000 0.25000000000000 2 -0.25000000000000 0.25000000000000 0.25000000000000 2 0.00000000000000 0.50000000000000 0.25000000000000 2 0.25000000000000 0.50000000000000 0.25000000000000 2 0.50000000000000 0.50000000000000 0.25000000000000 2 -0.25000000000000 0.50000000000000 0.25000000000000 2 0.00000000000000 -0.25000000000000 0.25000000000000 2 0.25000000000000 -0.25000000000000 0.25000000000000 2 0.50000000000000 -0.25000000000000 0.25000000000000 2 -0.25000000000000 -0.25000000000000 0.25000000000000 2 0.00000000000000 0.00000000000000 0.50000000000000 1 0.25000000000000 0.00000000000000 0.50000000000000 2 0.50000000000000 0.00000000000000 0.50000000000000 1 0.00000000000000 0.25000000000000 0.50000000000000 2 0.25000000000000 0.25000000000000 0.50000000000000 2 0.50000000000000 0.25000000000000 0.50000000000000 2 -0.25000000000000 0.25000000000000 0.50000000000000 2 0.00000000000000 0.50000000000000 0.50000000000000 1 0.25000000000000 0.50000000000000 0.50000000000000 2 0.50000000000000 0.50000000000000 0.50000000000000 10.5000000.0000000.5000000.000.4750000.0000000.4750000.000.4500000.0000000.4500000.00。0.0937500.0937500.1875000.000.0750000.0750000.1500000.000.0562500.0562500.1125000.000.0375000.0375000.0750000.000.0187500.0187500.0375000.000.0000000.0000000.0000000.00其他文件从HSE06中拷贝运行vasp4.从EIGENVAL中取权重为0的高对称点做能带图EIGENVAL： 2 2 1 10.2010146E+020.3845247E-090.3845247E-090.3845247E-090.5000000E-151.000000000000000E-004CARSi-Diamond 8127 80.0000000E+000.0000000E+000.0000000E+000.1562500E-01 1 -7.830170 2 5.377682 3 5.377684 4 5.377686 5 8.798907 6 8.798909 7 8.798911 8 9.6933140.2500000E+000.0000000E+000.0000000E+000.3125000E-01 1 -6.949546 2 1.136336 3 4.578605 4 4.578607 5 8.218893 6 9.797289 7 9.797292 8 13.389561。0.0000000E+000.0000000E+000.0000000E+000.0000000E+00 1 -7.830170 2 5.377682 3 5.377684 4 5.377686 5 8.798907 6 8.798909 7 8.798911 8 9.693314去掉前面36个0.0000000E+000.0000000E+000.0000000E+000.1562500E-01 1 -7.830170 2 5.377682 3 5.377684 4 5.377686 5 8.798907 6 8.798909 7 8.798911 8 9.693314留下后面91个高对称点，并改变高对称点数 2 2 1 10.2010146E+020.3845247E-090.3845247E-090.3845247E-090.5000000E-151.000000000000000E-004CARSi-Diamond 891 8。0.0000000E+000.0000000E+000.0000000E+000.0000000E+00 1 -7.830170 2 5.377682 3 5.377684 4 5.377686 5 8.798907 6 8.798909 7 8.798911 8 9.693314然后作图欢迎大家提意见 下载 GGA=B3-VASP5.2.rar 下载 GGA=PE-VASP5.2.rar 下载 GGA=PS-VASP5.2.rar 下载 PE-VASP4.6.rar 下载 Eu.rar举报删除此信息youzhizhe(站内联系TA)谢谢分享。liliangfang(站内联系TA)2楼: Originally posted by youzhizhe at 2011-10-21 15:50:47:谢谢分享。谢谢 大家一起学习kgdu(站内联系TA)谢谢分享！贺仪(站内联系TA)很好的资源 但是有一点没太明白，为很么在算能带的时候，前面取那些权重不为零的k点呢？一般的计算通过line模式直接得出的EIGENVAL,然后直接画图，你的权重为零的k点其实就是line模式产生的k点。能不能解释一下为什么要加那些权重不为零的点呢？liliangfang(站内联系TA)5楼: Originally posted by 贺仪 at 2011-10-21 18:22:54:很好的资源 但是有一点没太明白，为很么在算能带的时候，前面取那些权重不为零的k点呢？一般的计算通过line模式直接得出的EIGENVAL,然后直接画图，你的权重为零的k点其实就是line模式产生的k点。能不能解释一下 .这是我总结了很多人的计算方法，line模式HSE杂化这种方法算不了，通过各种尝试，访问vasp官网论坛和参考说明得到的。如果这些高对称点权重不为0就得不到你想要的能带图。简单的说就是不这样设置，通过杂化计算是算不动的或者是得不到结果。schutsilang(站内联系TA):tiger05:tiger05:tiger05:tiger05:tuhaidan(站内联系TA)顶一下，感谢分享！xx2008(站内联系TA):victory:victory:victory:victory:victory:victory:cavediger(站内联系TA)很好，支持！viebn(站内联系TA)LZ您好，您用的是vasp5.2嘛？请问您能不能把vasp5.2中Eu的PBE势发给我呢？我用的是vasp4.6，计算f电子的时候我的势库好像算不出来！谢谢楼主liliangfang(站内联系TA)11楼: Originally posted by viebn at 2011-10-24 11:36:19:LZ您好，您用的是vasp5.2嘛？请问您能不能把vasp5.2中Eu的PBE势发给我呢？我用的是vasp4.6，计算f电子的时候我的势库好像算不出来！谢谢楼主据我所知没有改变的吧viebn(站内联系TA)12楼: Originally posted by liliangfang at 2011-10-24 14:58:56:据我所知没有改变的吧谢谢了哈:hand:NKDMS(站内联系TA):victory:victory:ppwzq(站内联系TA)顶一下，感谢分享！ppwzq(站内联系TA)顶一下，感谢分享！ppwzq(站内联系TA)顶一下，感谢分享！ppwzq(站内联系TA)顶一下，感谢分享！xiaowenzhi16(站内联系TA)EVERY WEELwhw19850730(站内联系TA)6楼: Originally posted by liliangfang at 2011-10-21 18:48:44:这是我总结了很多人的计算方法，line模式HSE杂化这种方法算不了，通过各种尝试，访问vasp官网论坛和参考说明得到的。如果这些高对称点权重不为0就得不到你想要的能带图。简单的说就是不这样设置，通过杂化计算是 .你好很高兴遇到一个也会使用杂化泛函的，有几个问题要问你ENCUTFOCK和NKRED这两个参数什么物理意义，说明书说高精度计算可以省略，是不是省略了可定不出错，还有看你附件里文件使用杂化泛函时候ISYM = 0为什么去消对称性，还有kpoints文件为什么用G点，我看过一个资料好像不用G发散，想确定一下Automatic mesh0G4 4 40 0 0还有ENAUG = 800这个参数必须设定吗？期待回复谢谢。 查看全部评论相关推荐材料计算模拟实用技巧淘贴专辑VASP and MS淘贴专辑第一性原理淘贴专辑VASP淘贴专辑第一性原理相关文档淘贴专辑第一性原理计算经验淘贴专辑VASP算例淘贴专辑个人收集淘贴专辑VASP淘贴专辑vasp淘贴专辑计算-vasp淘贴专辑纳米及第一性原理方向淘贴专辑HSE & 杂化函数淘贴专辑ms and vasp学习淘贴专辑第一性&纳米&code淘贴专辑第一性淘贴专辑vasp淘贴专辑杂化泛函方面的研究淘贴专辑第一性淘贴专辑VASP 新计算方法淘贴专辑第一性原理和电化学淘贴专辑VASP-study淘贴专辑vasp淘贴专辑个人收藏淘贴专辑vasp淘贴专辑请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例天棚元帅发表于2013-02-17参杂后，晶格常数发生变化，那么计算出来的形成能准确吗xiaoboy185发表于2013-01-17VASP计算单胞GdN的能带结构gougou531发表于2013-01-10各位大侠，vasp下计算完能带如何作图lzlgcdc发表于2012-12-13VASP +U 计算能带redsnowolf发表于2012-