使用vasp计算dos(态密度)的一些思考 |
您所在的位置:网站首页 › k点密度跟平面波函数数量有关吗 › 使用vasp计算dos(态密度)的一些思考 |
一:计算前的准备工作
1,得到一个稳定的弛豫过后的POSCR结构 2,对结构进行一次静态计算 二:计算态密度的参数准备对于态密度计算INCAR参数的一些思考 ISTAT=1 这里我们选择1,但可以选择的还有0,在INCAR中不填,系统应该会默认为1 。他的目的是继承上一步的静态计算当中得到的WAVECAR,以减少计算的耗时。ICHARG=11 这一项还可以选择 0,1,2和12 选择用11进行计算的话就会读取我们在上一步静态计算过程中得到的CHGCAR并进行态密度计算这样,可以节约时间,以便平衡加大K点等其他参数所所带来的耗时增加。IBRION=-1 这里选择 -1 我们在结构优化的时候得到的体系已经是最稳定的这里选择-1就是把所有原子固定住。(思考: 为什么要固定?如果不固定是不是也能算呢?固定是为了节约计算量以节约时间吗?如果不固定也能算,那不固定对态密度峰值精度的影响大吗?)NSW=0 因为不需要调整结构所以不需要弛豫步数LCHCAR=.F. 数据有点大,可以输出,但一般不选择输出节约空间LWAVE=.F. 数据有点大,可以输出,但一般不选择输出节约空间ISMEAR=0 (疑问:我都是用了0,但我也不知道为什么)SIGMA=0.1 这个参数对能带或态密度中间的间隙影响很大,越接近0能带间隙就会越大,但相应计算的耗时也会增加(所以有时也不经怀疑自己,这弄计算的纯粹就是在凑)LORBIT=11 这个参数可以将每一个原子对应的轨道的贡献输出到OUTCAR里ISPIN=2 这个参数用于计算磁性,慎用,开启后会态密度图会成上下分布,对于磁性体系会出现上下分裂的情况NEDOS=3000 加大态密度计算的撒点数,这样可以让态密度的峰值更平顺EMIN=-20 决定态密度中x轴(能量)的最小取值范围EMAX=20 决定态密度中x轴(能量)的最大取值范围 其余的值没有写在INCAR里了,系统会自动使用默认值对于态密度计算KPOINTS取值的一些思考 取值尽量大一点 三:计算态密度需要的文件INCAR KPOINTS POSCAR (上一步静态得到的CONTCAR) POTCAR WAVCAR (上一步静态得到的) CHGCAR (上一步静态得到的) 四:计算态密度将任务投入计算列队 五:处理数据将计算得到的vasprun.xml使用p4vasp.exe处理,导出数据 后使用matlab或origin绘图 最后:留疑静态计算的数据也可以绘制态密度图,但经过态密度加强过的峰值会更平缓, |
今日新闻 |
推荐新闻 |
CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3 |