vasp结构优化后,有时需要提取或者观察优化过程中的特定离子步的结构,可以通过该shell脚本实现,只需要在优化结构所在文件夹下运行下列脚本,即可将每一个离子步输出为POSCAR-x文件!
#!/bin/bash
#zhaoc_chem@126.com
sum_atoms=$(cat XDATCAR | sed -n '7p' | awk '{for (i=1;i<=NF;i++) sum+=$i; print sum}')
sum_steps=$(grep F= OSZICAR | awk '{print $1}')
for j in $sum_steps
do
awk '{if (NR<8) print $0}' XDATCAR > POSCAR-$j
grep -xA $sum_atoms "Direct configuration= $j" XDATCAR >> POSCAR-$j
grep -xA $sum_atoms "Direct configuration= $j" XDATCAR >> POSCAR-$j
grep -xA $sum_atoms "Direct configuration= $j" XDATCAR >> POSCAR-$j
done
使用方法:
1.在优化结构的工作目录新建后缀为 .sh 文件(如 test.sh ),复制上面代码进文件保存退出
2.终端进入当前文件夹:cd ~/....../工作目录/
3.终端输入该命令授权脚本可执行:chmod 777 test.sh
4.终端输入:./test.sh
当前文件夹刷新即出现对应POSCAR文件,继而可以通过VESTA可视化(Download - VESTAhttp://www.jp-minerals.org/vesta/en/download.html)。