24小时热门版块排行榜    

查看: 2501  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

cuicui2397

新虫 (小有名气)

[求助] vasp阶段能测试求助,谢谢啊!急求。

看的侯博士的教程,里面有测试Si的截断能,但我自己去测试老是出现问题,求解啊!谢谢。
POSCAR文件:
Si-Diamond
5.430
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
2
Direct
0.0 0.0 0.0
0.25 0.25 0.25

KPOINTS文件:
Automatic generation
0
Monhkorst-Pack
9 9 9
0.0 0.0 0.0

用来确定ENCUT的脚本程序为run_ecut,其内容是:
#!/bin/sh
rm -rf WAVECAR
for i in 150 200 250 300 350 400
do
cat > INCAR <<!
SYSTEM = Si-Diamond
ENCUT = $i
ISTART = 0; ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo "ENCUT = $i eV"; mpirun -np 2 vasp
E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}''
echo $i $E >>comment
done

同时选择了赝势文件POTCAR,把以上四个文件POSCAR、KPOINTS、run_ecut和POTCAR放入我的计算目录下,在终端输入sh.run_ecut命令,在生成的comment里是这样的:
150 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
200 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
250 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
300 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
350 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
400 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,}
教程中第二列显示的是能量值,请问这是什么问题啊?谢谢啊!
回复此楼

» 猜你喜欢

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

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

微6凉

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by cuicui2397 at 2013-11-14 21:18:24
你好,首先谢谢你,我把符号改过之后,第二列出现了能量值。不过和教程上差别很大,你觉得是什么原因造成的啊?我的结果是这样的:
150 -10.744031
200 -10.795431
250 -10.822115
300 -10.836639
350 -10.84 ...

请问,我也是按照您那样算的,但是不知道为什么出不来结果,我用您那个命令,出现Command not found,我用我们那个提交任务的方法也不行,能告诉我为什么吗
13楼2015-08-03 11:20:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

微6凉

铁虫 (小有名气)

请问,我按照您这样的做法,总会出现这样一些结果Try `rm --help' for more information.
'un_ecut: line 4: syntax error near unexpected token `do
'un_ecut: line 4: `do  这是怎么回事
14楼2015-08-03 11:50:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cuicui2397 的主题更新
信息提示
请填处理意见