多语言展示
当前在线:1069今日阅读:2今日分享:31

VASP入门到精通[32]态密度和HOMO-LUMO

这一节小编写了个小脚本,提取DOS和HOMO-LUMO能级
工具/原料

vasp linux

方法/步骤
1

注意:这2个脚本是针对ISPIN=2的计算!!!!!!!!!!!!!小编写了两个脚本,一个名字叫data-dos一个名字叫gap.sh

2

第一个脚本是用来提取态密度数据的,脚本命令如下a=`head -6 DOSCAR|tail -1|awk '{print $3}'`b=$((a + 6))f=`awk '{if(NR==6)print $4}' DOSCAR`sed -n '7,'$b' p' DOSCAR > DOS.datawk '{print $1-'$f',$2}' DOS.dat > DOS-final.dat

3

注释:第一行:提取DOS数(就是取了多少个点)第二行:DOS的终止行,就是我们要提取a到b行的数据第三行:提取费米能级第四行:提取a到b行的DOS数据并写入DOS.dat(这时会生成一个DOS.dat文件)第五行:把提取的DOS.dat做费米修正,输出一个DOS-final.dat这个DOS-final.dat就是你要的DOS数据,放到origin就可以看到了

4

第一个脚本的使用方法:运行:bash data-dos你会发现多了两个文件:DOS.dat 和 DOS-final.dat把DOS-final.dat放到origin就可以了

5

第二个脚本代码如下:m=`awk '{if($1<0 &&($2>0.05 || $3>0.05))a=$1}END{print a}' DOS-final.dat`n=`awk '{if($1>0 &&($2>0.05 || $3>0.05))print $1}' DOS-final.dat|head -1`echo “HOMO = $m”echo “LUMO=$n”echo '$n' - '$m'|bc

6

注意这个脚本就五行哦,小编有截图,怕在这里编辑就自动换行,所以截图给大家

7

注释:第一行提取homo第二行是LUMO第三四五行就输出HOMO ,LUMO和能隙了

8

使用方法:bash gap.sh你会看到输出三行内容:HOMO=LUMO=能隙

9

有问题可以留言问小编哦!!

推荐信息