VASP的便捷命令

for i in {1..8};do cp 400 $((400+$i*50)) -r ;sed -i "10s/400/$((400+$i*50))/g" $((400+$i*50))/INCAR;done

上述为创造400到800的截断能的命令，前提需要创造一个400的截断能的文件

下面为配套的提交命令

for i in {1..9};do cd $((350+$i*50));sh run_vasp.sh;cd $OLDPWD;done

for i in *; do cd $i ; qsub sub4; cd $OLDPWD; done

输出截断能，能量，时间的提交命令

for i in *0; do echo -e $i "\t" $(grep User $i/OUTCAR |awk '{print $4}')"\t" $(grep ' without' $i/OUTCAR | tail -n 1 | awk '{print $4}'); done

画图的简要脚本

#!/usr/bin/python3
import matplotlib.pyplot as plt
import numpy as np

x,y,z=np.loadtxt('date1.dat',delimiter=',',usecols=(0,1,2),unpack=True)
plt.xlabel('ENCUT/eV')
plt.ylabel('TIME/s')
plt.plot(x,y,'rs-',linewidth=2.0)
plt.plot(x,z,'bs-',linewidth=2.5)
plt.show()

处理POSCAR文件晶胞大小快捷命令

for i in $(seq 10 2 20);do cp 8 $i -r;sed -i "3,5s/8.0/$i.0/g" $i/POSCAR;done