注意：这2个脚本是针对ISPIN=2的计算！！！！！！！！！！！！！

小编写了两个脚本，一个名字叫data-dos

一个名字叫gap.sh

第一个脚本是用来提取态密度数据的，脚本命令如下

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.dat

awk {print $1-$f,$2} DOS.dat DOS-final.dat

注释：

第一行：提取DOS数（就是取了多少个点）

第二行：DOS的终止行，就是我们要提取a到b行的数据

第三行：提取费米能级

第四行：提取a到b行的DOS数据并写入DOS.dat（这时会生成一个DOS.dat文件）

第五行：把提取的DOS.dat做费米修正，输出一个DOS-final.dat

这个DOS-final.dat就是你要的DOS数据，放到origin就可以看到了

第一个脚本的使用方法：

运行：

bash data-dos

你会发现多了两个文件:DOS.dat 和DOS-final.dat

把DOS-final.dat放到origin就可以了

第二个脚本代码如下：

m=`awk {if($10 ($20.05 || $30.05))a=$1}END{print a} DOS-final.dat`

n=`awk {if($10 ($20.05 || $30.05))print $1} DOS-final.dat|head -1`

echo “HOMO = $m”

echo “LUMO=$n”

echo $n - $m|bc

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

注释：第一行提取homo

第二行是LUMO

第三四五行就输出HOMO ，LUMO和能隙了

使用方法：

bash gap.sh

你会看到输出三行内容：

HOMO=

LUMO=

能隙

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