使用vasp计算dos（态密度）的一些思考 |
您所在的位置：网站首页 › chgcar文件 › 使用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 |