初学VASP中电子态密度计算设置参考
- 格式:doc
- 大小:34.50 KB
- 文档页数:6
VASP参数设置详解计算材料2010-11-30 20:11:32 阅读197 评论0 字号:大中小订阅转自小木虫,略有增减软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:对所计算的体系进行注释:SYSTEM●定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAV●定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子部分优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX–自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF●定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW –分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFG●定义态密度积分的方法和参数–smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量范围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBIT●其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
VASP参数设置详解计算材料2010-11-30 20:11:32 阅读197 评论0 字号:大中小订阅转自小木虫,略有增减软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:对所计算的体系进行注释:SYSTEM●定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWA V●定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子部分优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX–自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF●定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW–分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFG●定义态密度积分的方法和参数–smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量范围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBIT●其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
VASP参数设置详解计算材料2010-11-30 20:11:32 阅读197 评论0 字号:大中小订阅转自小木虫,略有增减软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:对所计算的体系进行注释:SYSTEM●定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWA V●定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子部分优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX–自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF●定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW–分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFG●定义态密度积分的方法和参数–smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量范围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBIT●其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
VASP参数设置详解计算材料2010-11-30 20:11:32 阅读197 评论0 字号:大中小订阅转自小木虫,略有增减软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:✍对所计算的体系进行注释:SYSTEM✍定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAV✍定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子部分优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX–自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF✍定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW–分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFG✍定义态密度积分的方法和参数–smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量范围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBIT✍其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
VASP参数设置详解要点VASP(Vienna Ab initio Simulation Package)是一种第一原理计算程序,用于计算材料性质和从头计算材料结构。
在进行VASP模拟时,合理设置参数非常重要,它们决定了模拟的准确性和效率。
下面将详细讨论几个关键的VASP参数设置要点。
1.设置能量截断(ENCUT):ENCUT是控制计算中的平面波能量截断的参数。
它应该尽量接近真实波函数的动能截断,以保证计算结果的准确度。
选择合适的ENCUT值非常关键,过低的值可能导致计算不收敛,过高的值则会造成计算时间过长。
一般建议从400eV开始进行尝试,然后根据计算的收敛性和计算结果调整。
2.设置k点密度(KPOINTS):k点密度是控制倒空间采样的参数。
k点密度越高,计算结果越准确,但计算时间也会增加。
为了在准确性和效率之间取得平衡,可以根据材料的对称性和大小进行合理的选择。
一般情况下,对于晶体,k点密度可以使用Reciprocal Space的自动生成程序,对于分子系统,可以使用Gamma Point + Monkhorst Pack方案。
3.设置电子步的最大迭代次数(NELM):NELM是控制电子步迭代收敛性的参数。
它决定了算法进行多少次最大迭代。
在计算过程中,电子步的总数是非常关键的。
如果电子步的迭代次数不足,可能会导致计算不收敛。
通常可以从60次开始进行尝试,如果计算结果不收敛,可以增加NELM的值。
4.设置计算精度(PREC):PREC参数是控制计算精度的参数。
该参数取值从粗到细分别为Low,Medium,High和Accuracy。
选择适当的计算精度可以在减少计算时间和提高计算结果准确性之间取得平衡。
一般情况下,可以从Medium开始尝试。
5.设置自洽迭代的收敛判据(EDIFF):EDIFF是控制自洽迭代收敛性的参数。
当自洽迭代前后两次总能量的变化低于EDIFF时,认为自洽迭代收敛。
合理设置EDIFF可以保证计算结果的准确性。
VASP参数设置详解计算材料 2010-11-30 20:11:32 阅读197 评论0 字号:大中小订阅转自小木虫,略有增减软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:l对所计算的体系进行注释:SYSTEMl定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAVl定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子部分优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX–自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFFl定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW –分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFGl定义态密度积分的方法和参数–smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量范围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBITl其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
VASP中电子态密度计算的流程VASP(Vienna Ab initio Simulation Package)是一种基于密度泛团理论(DFT)的第一性原理计算软件包,适用于从头计算材料的电子结构和相关性质。
电子态密度(Electronic Density of States, DOS)是VASP中一个重要的计算任务,它描述了材料中电子的能量分布情况,可以用来分析材料的能带结构、电导性、磁性等性质。
下面是VASP中计算电子态密度的一般流程:1.构建体系:首先需要确定要研究的体系的晶体结构。
可以通过实验数据、结构数据库或者其他理论方法得到体系的晶体结构,然后使用VASP提供的一些工具生成输入文件。
2.检查和准备输入文件:在进行计算之前,需要检查输入文件的正确性。
输入文件主要包括POSCAR(晶体结构)、POTCAR(势能文件)和KPOINTS(k点网格),还可以包括INCAR(控制参数)和CHGCAR(电荷密度)。
可以使用VASP提供的一些工具来生成这些文件。
3.设置计算参数:在INCAR文件中设置计算参数。
这些参数包括计算方式(GS、NSW等)、电子相关参数(ENCUT、EDIFF、ISMEAR等)和计算资源(NPAR、NCORE等)等。
4.进行自洽计算:运行VASP程序开始自洽计算。
自洽计算是指通过迭代寻找材料中所有电子的基态波函数和电子密度。
5.DOS计算:自洽计算完成后,可以进行DOS计算。
首先需要通过选择一个能量范围,确定所需的DOS信息。
然后在INCAR文件中设置相关参数,如要求计算PDOS(投影态密度)、LORBIT参数(需要计算轨道投影DOS)等。
6.执行DOS计算:运行VASP程序开始DOS计算。
程序会在给定的能量范围内计算电子态密度,并输出相应的结果。
7. 分析结果:根据VASP计算结果,可以通过一些可视化软件(如VESTA、XCrysDen等)绘制电子态密度的能带图、分析能带结构,进而分析材料的电子特性和相关性质。
vasp 态密度计算流程下载温馨提示:该文档是我店铺精心编制而成,希望大家下载以后,能够帮助大家解决实际的问题。
文档下载后可定制随意修改,请根据实际需要进行相应的调整和使用,谢谢!并且,本店铺为大家提供各种各样类型的实用资料,如教育随笔、日记赏析、句子摘抄、古诗大全、经典美文、话题作文、工作总结、词语解析、文案摘录、其他资料等等,如想了解不同资料格式和写法,敬请关注!Download tips: This document is carefully compiled by theeditor. I hope that after you download them,they can help yousolve practical problems. The document can be customized andmodified after downloading,please adjust and use it according toactual needs, thank you!In addition, our shop provides you with various types ofpractical materials,such as educational essays, diaryappreciation,sentence excerpts,ancient poems,classic articles,topic composition,work summary,word parsing,copy excerpts,other materials and so on,want to know different data formats andwriting methods,please pay attention!1. 准备输入文件POSCAR:包含晶体结构信息。
POTCAR:赝势文件。
VASP参数设置详解VASP(Vienna Ab initio Simulation Package)是一种用于计算材料的密度泛净近似(DFT)的第一性原理计算软件包。
它通过解Schrödinger方程来模拟材料的电子结构和相关性质。
在使用VASP进行计算时,合理地设置一些参数对于获得准确和可靠的计算结果至关重要。
下面是一些VASP参数的详细解释和设置建议。
1.ENCUT:电子截断能电子截断能决定了在计算中使用的总能量截断。
它代表了在周期性晶体中采样自由电子波函数所需的最高能量。
对于不同的系统,需要选择一个适当的ENCUT值来平衡计算精度和计算成本。
一般来说,较大的ENCUT可以提供更准确的结果,但也会增加计算的时间和计算资源的需求。
2.KPOINTS:k空间采样KPOINTS参数决定了在倒空间中如何采样能带结构或密度(扩展计算)等。
k点的选择通常取决于晶体结构的对称性和计算的目标。
大多数情况下,使用自动设置的KPOINTS能够提供合理的结果,但对于具有特殊性质(如表面性质或嵌入物性质)的系统,可能需要进行手动调整以获得更准确的结果。
3.ISMEAR:布洛赫函数展宽类型ISMEAR参数定义了用于计算导体体系的占有数矩阵的布洛赫函数的展宽类型。
对于绝缘体,可以使用较小的展宽(如ISMEAR=0),对于导体,使用较大的展宽(如ISMEAR=-5)以便准确地描述费米面附近的行为。
4. SIGMA:Methfessel-Paxton方法的展宽参数SIGMA参数是Methfessel-Paxton方法用来对扩展波函数(导体体系)进行展宽的参数。
较小的SIGMA可以提供更准确的结果,但可能会使计算更加耗时。
合理的SIGMA选择应遵循VASP文档中的建议范围。
5.EDIFF:能量差(收敛标准)EDIFF参数定义了能量收敛的标准。
计算迭代过程中,当两个连续的迭代之间的能量差小于EDIFF时,认为计算已经收敛。
VASP参数设置详解软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:l 对所计算的体系进行注释:SYSTEMl 定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAVl 定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子部分优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX –自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFFl 定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW–分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFGl 定义态密度积分的方法和参数– smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量范围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBITl 其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
VASP参数设置详解软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体l 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型l 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS,下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:对所计算的体系进行注释:SYSTEM定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAV定义电子的优化–平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG–电子部分优化的方法:ALGO,IALGO,LDIAG–电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG,BMIX_MAG,WC,INIMIX,MIXPRE,MAXMIX–自洽迭代步数和收敛标准:NELM,NELMIN,NELMDL,EDIFF定义离子或原子的优化–原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM,NSW –分子动力学相关参数:SMASS,TEBEG,TEEND,POMASS,NBLOCK,KBLOCK,PSTRESS–离子弛豫收敛标准:EDIFFG定义态密度积分的方法和参数–smearing方法和参数:ISMEAR,SIGMA–计算态密度时能量范围和点数:EMIN,EMAX,NEDOS–计算分波态密度的参数:RWIGS,LORBIT其它–计算精度控制:PREC–磁性计算:ISPIN,MAGMOM,NUPDOWN–交换关联函数:GGA,VOSKOWN–计算ELF和总的局域势:LELF,LVTOT–结构优化参数:ISIF–等等。
VASP参数设置详解软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体i 计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型i 计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)丨计算材料的光学性质丨计算材料的磁学性质丨计算材料的晶格动力学性质(声子谱等)丨表面体系的模拟(重构、表面态和STM模拟)丨从头分子动力学模拟丨计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS, 下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:对所计算的体系进行注释:SYSTEM定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAV定义电子的优化-平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG-电子部分优化的方法:ALGO,IALGO,LDIAG-电荷密度混合的方法:IMIX,AMIX,AMIN,BMIX,AMIX_MAG ,BMIX_MAG , WC , INIMIX , MIXPRE , MAXMIX-自洽迭代步数和收敛标准:NELM , NELMIN , NELMDL , EDIFF定义离子或原子的优化-原子位置优化的方法、移动的步长和步数:IBRION,NFREE,POTIM , NSW-分子动力学相关参数:SMASS , TEBEG , TEEND , POMASS , NBLOCK , KBLOCK , PSTRESS-离子弛豫收敛标准:EDIFFG定义态密度积分的方法和参数-smearing 方法和参数:ISMEAR , SIGMA-计算态密度时能量范围和点数:EMIN , EMAX , NEDOS-计算分波态密度的参数:RWIGS , LORBIT其它-计算精度控制:PREC-磁性计算:ISPIN , MAGMOM , NUPDOWN-交换关联函数:GGA , VOSKOWN-计算ELF和总的局域势:LELF , LVTOT-结构优化参数:ISIF-等等。
VASP参数设置详解软件主要功能:采用周期性边界条件(或超原胞模型)处理原子、分子、团簇、纳米线(或管)、薄膜、晶体、准晶和无定性材料,以及表面体系和固体ι计算材料的结构参数(键长、键角、晶格常数、原子位置等)和构型ι计算材料的状态方程和力学性质(体弹性模量和弹性常数)l 计算材料的电子结构(能级、电荷密度分布、能带、电子态密度和ELF)l 计算材料的光学性质l 计算材料的磁学性质l 计算材料的晶格动力学性质(声子谱等)l 表面体系的模拟(重构、表面态和STM模拟)l 从头分子动力学模拟l 计算材料的激发态(GW准粒子修正)计算主要的四个参数文件:INCAR ,POSCAR,POTCAR ,KPOINTS, 下面简要介绍,详细权威的请参照手册INCAR文件:该文件控制VASP进行何种性质的计算,并设置了计算方法中一些重要的参数,这些参数主要包括以下几类:对所计算的体系进行注释:SYSTEM定义如何输入或构造初始的电荷密度和波函数:ISTART,ICHARG,INIWAV定义电子的优化-平面波切断动能和缀加电荷时的切断值:ENCUT,ENAUG-电子部分优化的方法:ALGO,IALGO,LDIAG-电荷密度混合的方法:IMlX,AMIX,AMIN,BMlX,AMIX_MAG ,BMIX_MAG , WC , INlMIX , MlXPRE , MAXMIX-自洽迭代步数和收敛标准:NELM , NELMIN , NELMDL , EDIFF定义离子或原子的优化-原子位置优化的方法、移动的步长和步数:IBRION , NFREE , POTIM , NSW-分子动力学相关参数:SMASS , TEBEG , TEEND , POMASS , NBLOCK , KBLoCK , PSTRESS-离子弛豫收敛标准:EDIFFG定义态密度积分的方法和参数-Smearing 方法和参数:ISMEAR , SIGMA-计算态密度时能量范围和点数:EMIN , EMAX , NEDOS-计算分波态密度的参数:RWIGS , LORBIT其它-计算精度控制:PREC-磁性计算:ISPIN , MAGMOM , NUPDOWN-交换关联函数:GGA , VOSKOWN-计算ELF和总的局域势:LELF , LVTOT-结构优化参数:ISIF-等等。
初学VASP中电子态密度计算基本设置参考主要分成三步:一、结构优化;二、静态自洽计算;三、非自洽计算以Al-FCC为例子第一步结构优化输入文件(INCAR, POTCAR, POSCAR, KPOINT)INCAR文件System=AlISTART=0ISMEAR=1SIGMA=0.2ISPIN=2GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=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.00.5 0.5 0.00.5 0.0 0.50.0 0.5 0.5KPOINT 文件Automatic generationMohkorst 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中,然后运行VASP。
vasp态密度计算流程Title: VASP Density of States Calculation ProcedureTitle: 态密度计算流程Firstly, to calculate the density of states (DOS) using VASP, it is necessary to have performed a structural optimization or a static calculation.The essential inputs required are the POSCAR file which contains information about the structure, and the POTCAR file which includes information about the pseudopotentials.Additionally, the INCAR file should be properly set up with necessary parameters such as the electronic configuration and the energy range for the calculation.首先,使用VASP计算态密度(DOS),必须先进行结构优化或静态计算。
所需的基本输入包括包含结构信息的POSCAR文件和包含伪势信息的POTCAR 文件。
此外,INCAR文件应正确设置所需参数,如电子配置和计算的能量范围。
ext, prepare the KPOINTS file which defines the grid points in the Brillouin zone for the calculation of the DOS.The grid points should be dense enough to obtain a smooth DOS curve, but not too dense to make the calculation computationally expensive.The grid points can be optimized by trial and error or using automatic grid generation methods.接下来,准备KPOINTS文件,它定义了在计算DOS时布里渊区内的网格点。
初学VASP中电子态密度计算基本设置参考主要分成三步:一、结构优化;二、静态自洽计算;三、非自洽计算以Al-FCC为例子
第一步结构优化
输入文件(INCAR, POTCAR, POSCAR, KPOINT)
INCAR文件
System=Al
ISTART=0
ISMEAR=1
SIGMA=0.2
ISPIN=2
GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 IBRION=2
NSW=50
ISIF=2 (OR 3)
NPAR=10
POTCAR 文件直接在势库中拷贝
POSCAR文件
Al
4.05
1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0
4
Direct
0.0 0.0 0.0
0.5 0.5 0.0
0.5 0.0 0.5
0.0 0.5 0.5
KPOINT 文件
Automatic generation
Mohkorst Pack
15 15 15
0.0 0.0 0.0
第二步静态自洽计算
INCAR:PREC = Medium,ISTART = 0,ICHARG = 2,ISMEAR = -5输入文件(INCAR, POTCAR, POSCAR, KPOINT)
INCAR文件
System=Al
ISTART=0
ISMEAR=1
SIGMA=0.2
ISPIN=2
GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01 #IBRION=2
#NSW=50
#ISIF=2 (OR 3)
NPAR=10
POTCAR 文件直接在势库中拷贝
POSCAR文件
Al
4.05
1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0
4
Selective Dynamic
Direct
0.0 0.0 0.0 T T T
0.5 0.5 0.0 T T T
0.5 0.0 0.5 T T T
0.0 0.5 0.5 T T T
KPOINT 文件
Automatic generation
Mohkorst Pack
15 15 15
0.0 0.0 0.0
第二步计算是在结构优化的结果上进行的,所以开始第二步的时候,将第一步中的输入文件INCAR, POTCAR, POSCAR, KPOINT 以及C* 文件放入静态自洽计算中去,并且将CONTCAR 拷贝到POSCAR中,然后运行VASP。
计算结果中的Fermi能是准确的,需要记下。
第三步非自洽计算
INCAR:PREC = Medium,ICHARG = 11,ISMEAR = -5,LORBIT = 10 或者11(这时可不设RWIGS),ISTART = 1
在第二步自洽计算的基础上进行,修改输入文件INCAR, POTCAR, POSCAR, KPOINT。
INCAR文件
System=Al
ISTART=1
ISMEAR=-5
SIGMA=0.2
ICHARG=11
RWIGS=1.402
ISPIN=2
GGA=91; VOSKOWN=1; EDIFF=0.1E-05; EDIFFG=-0.01
#IBRION=2
#NSW=50
#ISIF=2 (OR 3)
NPAR=10
POTCAR 文件直接在势库中拷贝
POSCAR文件
Al
4.05
1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0
4
Selective Dynamic
Direct
1.0 0.0 0.0 T T T
0.5 0.5 0.0 T T T
0.5 0.0 0.5 T T T
0.0 0.5 0.5 T T T
KPOINT 文件
Automatic generation
Mohkorst Pack
21 21 21
0.0 0.0 0.0
VASP中文手册:在进行能带和DOS计算时,ISMEAR 不能使用N阶MP方法。
因为MP方法在空轨道上有负的占据,所以求得的能带和DOS是不正确的。
但是从其它地方看到---“提示:在计算能带结构时,采ISMEAR = 0或1对结果的影响非常小,可以认为是一样的。
但是不能采用ISMEAR = -5 或-4。
”ISMEAR到底多少?
计算能带:ICHARG = 11
导体的话,用ISMEAR=1;
半导体或绝缘体,用ISMEAR=0 。
计算DOS:ICHARG = 11
ISMEAR = -5
计算的时候,金属可用0、1,非金属不要大过0,体材料可用-4、-5(面的话就用-1、0吧?)
参考http://210.45.66.120/doc/vasp5.0/node224.html
设置完成后进行计算,计算完后,得到包含了态密度值的DOSCAR文件,采用split_dos对态密度文件DOSCAR进行分割,得到总态密度DOS0,各个原子的分波态密度DOS1,DOS2……。
另外在运行split_dos程序对DOSCAR文件分割时,要保证当前目录下有对应的OUTCAR和POSCAR文件。
分割后的DOS0,DOS1…等文件的能量值是以费米能级作为能量参考零点。
DOS0的第一列数据是能量值,单位为eV;第二列数据是总态密度的值,单位State/eV.unit cell;第三列数据是总态密度的积分值,也就是电子数,单位为electrons。
DOS1是第一个原子的分波态密度值,其中的第一列数据是能量值,单位为eV;第二、三、四列数据分别对应于s、p、d态的分波态密度值,单位为State/eV.atom。
其他的DOS文件与DOS1类似。