VASP处理split_dos-推荐下载
- 格式:pdf
- 大小:739.18 KB
- 文档页数:34
VASP计算DOS和能带个人总结一:VASP计算DOS和能带1.计算DOS①POSCAR②POTCAR③KPOINTS(建议以Gamma为中心取点,通常K×a≈45即可)④INCAR(越简洁越好)第一步:结构优化SYSTEM=**ISTART=0ENCUT=500(最好对其进行测试)EDIFF=1E-5EDIFFG=-0.01NSW=100ISIF=2IBRION=2【优化后计算DOS可以一步完成,也可以分为两步来完成,主要是计算量涉及到计算时间的差别】第二步:静态自洽(此时可稍微降低K点数,用第一步优化得到的CONTCAR作为POSCAR进行计算)SYSTEM=**ISTART=0PREC=AccurateEDIFF=1E-5EDIFFG=-0.01ENCUT=500ISMEAR=-5LCHARG=.TRUE.注意:此时得到的E-feimi是准确的,需要记下(grep ‘E-fermi’OUTCAR)第三步:非自洽计算(采用高密度K点)SYSTEM=**ISTART=1ICHARG=11LMAXMIX=2/4/6(VASP官网原话:If ICHARG is set to 11 or 12, it is strongly recommened to set LMAXMIX to twice the maximum l-quantum number in the pseudpotentials. Thus for s and p elements LMAXMIX should be set to 2, for d elements LMAXMIX should be set to 4, and for f elements LMAXMIX should be set to 6)PREC=AccurateEDIFF=1E-5EDIFFG=-0.01ENCUT=500(截断能最好与上一步保持一致)ISMEAR=-5LORBIT=10/11(推荐11,可以得到能级分裂的数据)优化后计算DOS一步完成:(采用高密度K点)SYSTEM=**ISTART=1PREC=AccurateEDIFF=1E-5EDIFFG=-0.01ENCUT=500ISMEAR=-5LORBIT=10/112.计算能带①POSCAR②POTCAR③KPOINTS:使用Line-mode格式,给出高对称性K点之间的分割点数,分割越密,路径积分就越准确。
VASP经典学习教程有用VASP(Vienna Ab initio Simulation Package)是一种用于固体材料计算的第一性原理计算软件包。
它使用密度泛函理论和平面波基组进行计算,可以预测材料的结构、能带、力学性质等基本属性。
本文将介绍VASP的经典学习教程,帮助初学者快速入门。
1.VASP的安装与基本操作-输入文件和输出文件:介绍VASP的常用输入文件和输出文件,以及它们的格式和含义。
-运行VASP计算:教授如何编写VASP运行脚本,以及如何使用命令行界面运行VASP计算。
2.VASP的输入参数和设置-INCAR文件:介绍VASP的主要输入文件INCAR的各种参数和选项,如体系的外部压力、电子迭代的收敛准则等。
-POTCAR文件:讲解VASP的赝势文件POTCAR的作用和用法,以及如何选择合适的赝势。
-KPOINTS文件:讲解KPOINTS文件对计算结果的影响,以及如何选择合适的K点网格。
3.VASP的基本计算-结构优化计算:教授如何进行结构优化计算,寻找稳定的材料晶格参数和原子位置。
-能带计算:讲解如何计算材料的能带结构,以及如何分析能带图和带隙。
-DOS计算:介绍如何计算材料的态密度,以及如何分析态密度图和能带图。
4.VASP的高级计算-弛豫计算:讲解如何进行离子和电子的同时弛豫计算,得到材料的稳定结构和力学性质。
-嵌入原子计算:介绍如何在材料中嵌入原子,并计算嵌入原子的相互作用能。
-软件接口和后处理:讲解VASP与其他软件(如VASPKIT、VESTA等)的接口,以及如何进行后处理分析。
5.VASP的实际应用-表面计算:介绍如何计算材料的表面能和表面形貌。
-催化剂计算:讲解如何通过VASP计算催化剂的吸附能和反应能垒,以预测其催化活性。
-界面计算:讲解如何计算材料的界面能和界面结构。
通过以上内容,初学者可以掌握VASP的基本原理和使用方法,并能在实际应用中进行一些基本的材料计算。
V ASP计算DOSV ASP中电子态密度计算的流程主要分成三步:一、结构优化;二、静态自洽计算;三、非自洽计算以Al-FCC为例子第一步结构优化阿输入文件(INCAR, POTCAR, POSCAR, KPOINT)INCAR文件System=AlISTART=0ISMEAR=1SIGMA=0.2ISPIN=2GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01IBRION=2NSW=50ISIF=2 (OR 3)NPAR=10POTCAR 文件直接在势库中拷贝POSCAR文件Al4.051.0 0.0 0.00.0 1.0 0.00.0 0.0 1.04Direct0.0 0.0 0.0 0.5 0.5 0.0 0.5 0.0 0.5 0.0 0.5 0.5KPOINT 文件Automatic generation 0Mohkorst Pack15 15 150.0 0.0 0.0第二步静态自洽计算INCAR:PREC = Medium,ISTART = 0,ICHARG = 2,ISMEAR = -5 输入文件(INCAR, POTCAR, POSCAR, KPOINT)INCAR文件System=AlISTART=0ISMEAR=1SIGMA=0.2ISPIN=2GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01#IBRION=2#NSW=50#ISIF=2 (OR 3)NPAR=10POTCAR 文件直接在势库中拷贝POSCAR文件Al4.051.0 0.0 0.00.0 1.0 0.00.0 0.0 1.04Selective DynamicDirect0.0 0.0 0.0 T T T0.5 0.5 0.0 T T T0.5 0.0 0.5 T T T0.0 0.5 0.5 T T TKPOINT 文件Automatic generationMohkorst Pack15 15 150.0 0.0 0.0第二步计算是在结构优化的结果上进行的,所以开始第二步的时候,将第一步中的输入文件INCAR, POTCAR, POSCAR, KPOINT 以及C* 文件放入静态自洽计算中去,并且将CONTCAR 拷贝到POSCAR中,然后运行V ASP。
打包压缩命令:tar zcvf 文件名.tar.gz 源文件名采用link方式避免重复的文件浪费内存:ln –s 源文件命名Eg:ln –s ../optic/MME ./EuO.mme交互式绘图工具gnuplot: 命令行打:gnuplot进入格式:plot “文件名”退出:quit求磁矩:getmag如果由于节点掉线在提交任务后秒退,不输出outcar可以指定节点提交任务LJRS -l nodes=c0104:ppn=4chmod +x 名称——使脚本可以执行构造potcar,以A和B元素为例:如果是以Z为拓展名的文件:zcat A/POTCAR.Z B/POTCAR.Z > POTCAR如果是解压后的potcar文件:cat A/POTCAR B/POTCAR > POTCARINCAR中的RWIGS通过POTCAR文件获得单位晶胞体积:grep “vol”OUTCAR自动计算加应力情况下的最优化情况:C/a :vaspcaopt softmode-e -=*(画曲线)如果要看某种材料是FM还是AFM,需要以相同结构计算一次,看能量哪个低。
将POSCAR/CONTCAR/CHGCAR装换成xcrysden(进入xcrysden文件执行./xcrysden)可读取的形式Eg:(BFO)v2xsf CHGCAR -1 83 -2 26 -3 8 -dv2xsf POSCAR -1 83 -2 26 -3 8 -dv2xsf CONTCAR -1 83 -2 26 -3 8 -d晶胞放大时,k点需要等比例的缩小画出曲线图:ISMEAR如果是半导体/绝缘体取-5,如果是金属取1,SIGMA=0.2(一般不改变)如果求DOS,则ISMEAR=-5.如果求band,则ISMEAR=1.1、Relax ISTART=0,ICHARG=2, ISIF=3,NSW=200,EDIFFG=-1*10-3,IBRION=22、Scf NSW=0(关闭结构优化) 用Relax后的CONTCAR替换POSCAR3、DOS 添加scf后的CHGCAR,ISTART=1,ICHARG=11,k放大一倍DOSCAR 第六行:Emax Emin Emin与Emax之间点的数目Ef第七行:能量总的态密度(spin up)总的态密度(spin down)态密度积分(up)态密度积分(down)后面按原子分:能量s轨道态密度(spin up)s轨道态密度(spin down)p轨道态密度(spin up)p轨道态密度(spin down)d轨道态密度(spin up)d轨道态密度(spindown)f轨道态密度(spin up)f轨道态密度(spin down)Split_dos 对dos按原子区分4、band 与DOS所需文件一致,KPOINTS需要使用line模式手动输入。
VASP使用总结VASP(Vienna Ab initio Simulation Package)是一款基于密度泛函理论(DFT)的第一性原理计算软件,主要用于材料科学和凝聚态物理领域的计算。
它提供了丰富的功能和工具,可以用于模拟和研究各种材料的物理和化学性质。
以下是对VASP使用的总结:1.输入文件的准备在进行VASP计算之前,首先需要准备好输入文件。
VASP使用的输入文件包括POSCAR、INCAR、POTCAR等。
POSCAR文件用于定义晶体结构和原子坐标,INCAR文件用于定义计算参数和设置计算方法,POTCAR文件用于定义原子的赝势。
2.材料结构的优化VASP可以通过结构优化计算来确定材料的最稳定结构。
结构优化计算通过改变原子位置和晶胞大小,寻找最低能量的结构。
可以使用ISIF 参数来设置优化类型,如禁止移动原子、禁止改变晶胞大小等。
3.能带结构的计算VASP可以计算材料的能带结构,从而提供关于能带轨道和能带间隙的信息。
能带结构计算需要先进行结构优化计算,然后再进行自洽计算和能带计算。
可以通过设置KPOINTS和NBANDS参数来控制计算的精度和效率。
4.密度状态的计算VASP可以计算材料的密度状态,包括电荷密度、电荷分布和电子态密度等。
通过密度状态计算,可以了解材料的电子结构和性质。
可以通过设置LSORBIT、IALGO和NPAR等参数来控制计算的模式和效率。
5.势能面的计算VASP可以计算材料的势能面,并通过构建势能面图像来显示材料的稳定性和反应性。
势能面计算需要进行结构优化计算,然后通过改变原子位置和晶胞大小来势能面上的最低能量和结构。
6.热力学性质的计算VASP可以通过计算自由能、热容和热膨胀系数等热力学性质来了解材料的热稳定性和热响应。
热力学性质的计算需要进行结构优化计算和自洽计算,然后使用VASP提供的工具和脚本进行热力学性质的分析和计算。
7.计算结果的解析和可视化VASP提供了丰富的工具和脚本,可以用于解析和可视化计算结果。
实例1如何用gcc编译dos-procar.f来处理VASP计算得到的procar如题,VASP计算pdos,含d电子分裂,想在生成的PROCAR中提取数据,然后画图,据说dos-procar.f(/bbs/attachment.php?aid=407340)可以做到。
请问如何编译dos-procar.f?我的linux没有ifort,而有gcc。
以下是我的PROCAR文件的一部分:PROCAR lm decomposed + phase factor# of k-points: 63 # of bands: 184 # of ions: 48k-point 1 : 0.00000000 0.00000000 0.00000000 weight = 0.00165289band 1 # energy -0.44451898 # occ. 2.00000000ion s py pz px dxy dyz dz2 dxz dx2 tot1 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0182 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0183 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0184 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0185 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0186 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0187 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0188 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.0189 0.017 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.002 0.018 ...............................................................................................................edited by zhjh888888 z> which f77which: no f77 in (/shared_scratch/pgi.714/linux86-64/7.1-4/bin/:/home_n22/xqdai2/bin:/usr/local/bin:/usr /bin:/usr/X11R6/bin:/bin:/usr/games:/opt/gnome/bin:/opt/kde3/bin:/usr/lib/mit/bin:/usr/ lib/mit/sbin:/usr/local/maui/sbin:/usr/local/maui/bin)on > gcc -o dos-procar.x dos-procar.fIn file dos-procar.f:1implicit real*8(a-h,o-z)1Error: Non-numeric character in statement label at (1)In file dos-procar.f:1implicit real*8(a-h,o-z)1Error: Unclassifiable statement at (1)In file dos-procar.f:1101000 continue1Error: Bad continuation line at (1)In file dos-procar.f:1101000 continue1Error: Unclassifiable statement at (1)In file dos-procar.f:113do k = 1, nk1In file dos-procar.f:902Error: Variable 'k' at (1) cannot be redefined inside loop beginning at (2) In file dos-procar.f:117do k = 1,nk1In file dos-procar.f:902Error: Variable 'k' at (1) cannot be redefined inside loop beginning at (2) In file dos-procar.f:121101 format(10x,f9.5)1Error: Bad continuation line at (1)In file dos-procar.f:121101 format(10x,f9.5)1Error: Unclassifiable statement at (1)In file dos-procar.f:122102 format(f10.5)1Error: Bad continuation line at (1)In file dos-procar.f:122102 format(f10.5)1Error: Unclassifiable statement at (1)In file dos-procar.f:123103 format(20a4)1Error: Bad continuation line at (1)In file dos-procar.f:123103 format(20a4)1Error: Unclassifiable statement at (1)In file dos-procar.f:124104 format(16x,i3,20x,i5,19x,i4)1Error: Bad continuation line at (1)In file dos-procar.f:124104 format(16x,i3,20x,i5,19x,i4)1Error: Unclassifiable statement at (1)In file dos-procar.f:125105 format(10x,i3,5x,3f11.8,13x,f11.8) 1Error: Bad continuation line at (1)In file dos-procar.f:125105 format(10x,i3,5x,3f11.8,13x,f11.8) 1Error: Unclassifiable statement at (1)In file dos-procar.f:126106 format(4x,i4,9x,f14.8,7x,f12.8) 1Error: Bad continuation line at (1)In file dos-procar.f:126106 format(4x,i4,9x,f14.8,7x,f12.8)1Error: Unclassifiable statement at (1)In file dos-procar.f:127107 format(3x,4f7.3)1Error: Bad continuation line at (1)In file dos-procar.f:127107 format(3x,4f7.3)1Error: Unclassifiable statement at (1)In file dos-procar.f:134do k = 1,nk-11In file dos-procar.f:902Error: Variable 'k' at (1) cannot be redefined inside loop beginning at (2)In file dos-procar.f:148do k=1,nk1In file dos-procar.f:902Error: Variable 'k' at (1) cannot be redefined inside loop beginning at (2)In file dos-procar.f:1899000 continue1 Error: End of nonblock DO statement at (1) is interwoven with another DO loop In file dos-procar.f:192end1Error: END DO statement expected at (1)Error: Unexpected end of file in 'dos-procar.f'请问运行gcc -o dos-procar.x dos-procar.f 后,出现了以上错误是什么原因?谢谢实例2关于V ASP计算PDOS的疑问在计算完TDOS和PDOS之后,不知道怎么将DOS图画出来。
vasp软件包使用入门指南VASP软件包使用入门指南本文档旨在为用户提供VASP软件包的使用入门指南。
以下是本文档的内容概要:1.简介1.1 VASP软件包概述1.2 VASP的功能和应用领域1.3 VASP的系统要求和安装方法2.输入文件准备2.1 结构文件的准备2.2 参数文件的设置2.3 能带计算的准备3.输入文件编辑和提交3.1 INCAR文件的编辑3.2 KPOINTS文件的编辑3.3 POTCAR文件的选择和拼接3.4 批量计算的脚本编写和提交4.输出文件解读4.1 OUTCAR文件解读4.2 VASP计算输出的主要结果解析5.常见问题和错误分析5.1 VASP计算中的常见错误和警告5.2 错误日志文件的解读和分析5.3 VASP计算结果的正确性判断6.进阶使用技巧6.1 计算参数优化方法6.2 并行计算和性能优化6.3 嵌入VASP到自己的代码中7.示例和案例分析7.1 基础输入文件示例7.2 常见VASP计算任务案例分析8.扩展资料和学习资源8.1 VASP官方文档和论坛8.2 相关书籍和资料推荐8.3 其他学习资源和机构推荐附件:1.VASP软件包安装包2.示例输入文件和脚本3.相关代码和案例法律名词及注释:1.VASP - Vienna Ab initio Simulation Package,维也纳从头模拟软件包。
2.INCAR - Input CARd(输入卡)的缩写,VASP输入文件之一,用于设置计算参数。
3.KPOINTS - VASP输入文件之一,用于设置k点和布点方案。
4.POTCAR - VASP输入文件之一,包含了计算所需的原子赝势和波函数。
本文档涵盖了VASP软件包的基本用法和高级技巧,并提供了示例和案例分析帮助用户更好地理解和应用该软件包。
如需进一步了解VASP的详细信息,请参阅附件中的相关资料和学习资源。