将vasp结构优化的每一步(XDATCAR文件)转化成POSCAR文件,通过VESTA实现可视化

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)。

  • 2
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值