pdos计算与画图

pdos计算与画图


1.计算opt
2.计算static         # 以上两步请参考之前计算过程
3.计算pdos
 
#静态计算完成后,拷贝 INCAR KPOINT CONTCAR(拷贝之后需改名为POSCAR) POTCAR vasp.lsf WAVECAR CHGCAR至pdos文件夹
#其中INCAR 参数必须设为下面值,其他参数和static中务必一致!否则后期脚本运行出错非常麻烦
ISTART=1
ICHARG=11
ISMEAR=-5
NEDOS=5000
LORBIT=11
#提交计算
4. pdos数据处理(以下都在pdos目录下处理)

  • 将POSCAR 导出至windows下,用VESTA 打开,找到你想画哪个原子(如第一个Ti原子)的dos图,记住序号为1。pdos计算与画图
  • 将附件中的脚本plotdos1.0脚本拷贝至服务器中pdos目录下,然后运行命令pdos计算与画图
 
chmod 700 plotdos1.0
sed '13,13c Num=1' plotdos1.0 > plotdos
#其中Num=1,中的‘1’是前一步中记录的原子序数,各自更改,请直接复制,避免添加空格等
chmod 700 plotdos
./plotdos 
# 运行完成后会出现 dos1.txt文件,切忌重复运行./plotdos,否则会多次写入数据,无法画图
运行完成如下
pdos计算与画图


  • 拷贝dos1.txt到Origin画图
    pdos计算与画图
  • 成功导入后结果如下
    pdos计算与画图