24小时热门版块排行榜    

查看: 1089  |  回复: 1

lulutiantian

新虫 (小有名气)

[求助] 截断能测试脚本出错!急求

大家帮我看看吧,我计算Au的块体,做截断能测试。以下是我的ecut 脚本,为了查看能不能运行,所以截断能暂时设了两个
#!/bin/sh
for i in 360 410
do

cat > INCAR <<!
SYSTEM = Au
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = 1
PREC = Accurate
!
echo "-------------------------------"
echo "ENCUT = $i eV"
mpiexec -n 2 vaspnew >> out &
wait
E=`grep "energy without entropy=" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $4 }'`
echo $i $E >>E-encut.dat

done

运行的时候将INCAR KPOINTS POSCAR POTCAR 还有脚本ENCUT放在一起运行了,但是出来的结果里面没有能量。光有截断能
请做过截断能测试的虫友帮帮忙吧!问题出在哪了?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lulutiantian

新虫 (小有名气)

后来去掉了这一行
echo "-------------------------------"
并将E=`grep "energy without entropy=" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $4 }'`改成
E=`grep "entropy=" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $4 }'`
运行时去掉了INCAR。可以正常运行了
2楼2013-11-21 16:26:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lulutiantian 的主题更新
信息提示
请填处理意见