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。
• 将附件中的脚本plotdos1.0脚本拷贝至服务器中pdos目录下，然后运行命令

 

chmod 700 plotdos1.0

sed '13,13c Num=1' plotdos1.0 > plotdos

#其中Num=1，中的‘1’是前一步中记录的原子序数，各自更改，请直接复制，避免添加空格等

chmod 700 plotdos

./plotdos 

# 运行完成后会出现 dos1.txt文件，切忌重复运行./plotdos,否则会多次写入数据，无法画图

运行完成如下

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