| 查看: 2499 | 回复: 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,} 教程中第二列显示的是能量值，请问这是什么问题啊？谢谢啊！ |
回复此楼» 猜你喜欢
论文终于录用啦！满足毕业条件了
已经有17人回复
-
不自信的我
已经有5人回复
-
磺酰氟产物，毕不了业了！
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp运行，第一离子步结束后自动停了
已经有7人回复
- vasp中，做截断能测试时，测试收敛的标准
已经有10人回复
- 深圳超算中心vasp编译及任务提交问题
已经有25人回复
- 关于VASP一些琐事的求助（K点，吸附优化等）
已经有12人回复
- vasp并行测试时出错。
已经有7人回复
- vasp并行出错，急于求助！！！
已经有13人回复
- vasp计算：优化能收敛，但同样的设置静态计算不收敛
已经有3人回复
- VASP计算时的几个细节设置问题，希望有经验的人过来指点一下！
已经有23人回复
- VASP 截断能测试 急求
已经有8人回复
- 【求助成功】vasp测试K点收敛性的脚本问题
已经有20人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】求将vasp的输出文件转化成xyz坐标的脚本
已经有11人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【求助】vasp做表面计算时真空层的厚度
已经有13人回复
- 【求助】VASP能对体系加电场计算吗？
已经有18人回复
- 【讨论】vasp结构优化之前的测试
已经有11人回复
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cuicui2397(fzx2008代发): 金币+3 2013-11-14 21:25:06
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:10
cuicui2397: 金币+5, ★★★★★最佳答案 2013-11-23 10:10:40
感谢参与，应助指数 +1
cuicui2397(fzx2008代发): 金币+3 2013-11-14 21:25:06
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:10
cuicui2397: 金币+5, ★★★★★最佳答案 2013-11-23 10:10:40
|
那个，有可能是符号有问题，那个E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}'' 是`（ESC下面那个键，数字1左面那个键），而不是'(单引号)，你试试看 |
» 本帖已获得的红花（最新10朵）
cuicui2397
2楼2013-11-14 20:30:12
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:17
cuicui2397: 金币+12, ★★★很有帮助 2013-11-23 10:11:24
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:17
cuicui2397: 金币+12, ★★★很有帮助 2013-11-23 10:11:24
| 赝势是不是不一样，你有没有选正确的赝势？侯老师的那个教程，有没有说选什么赝势？ |
» 本帖已获得的红花（最新10朵）
4楼2013-11-14 23:06:04
7楼2013-11-15 12:48:07
12楼2013-12-02 12:26:06